ISRN ThermodynamicsVolume 2012 (2012), Article ID 878279, 4 pagesdoi:10.5402/2012/878279

Research Article

A Case Study of Nanofiltration Using an Irreversible Thermodynamics-Based Model

G. R. Vakili-Nezhaad1,2 and Z. Akbari²

¹Department of Petroleum & Chemical Engineering, College of Engineering, Sultan Qaboos University, Muscat 123, Oman²Institute for Nanoscience and Nanotechnology, University of Kashan, Kashan 87317-51167, Iran

Received 7 December 2011; Accepted 9 January 2012

Academic Editor: Z. Slanina

Copyright © 2012 G. R. Vakili-Nezhaad and Z. Akbari. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The extended Spiegler-Kedem model that can be solved by using genetic algorithm (GA) has been used for study of the rejection of charged and uncharged solutes in a nanofiltration process. In this work, the accuracy and applicability of this model were tested with the data obtained from the nanofiltration of a real solution. The permeated solution was obtained from the whey ultrafiltration system (ultrafiltered whey). The studied process separates lactose and mineral salts from ultrafiltered whey. Moreover, the ultrafiltered whey nanofiltration system has been simulated with higher precision compared to the other sources of existing data.

1. Introduction

One of the applications of the nanofiltration membranes is the separation of lactose from the whey. Lactic acid is one of the most valuable retentive food material which has a high price. The whey is the main byproduct which is obtained from the cheese production. It contains high concentration of organic matter, mainly proteins and lactose. Study of the mathematical mass transfer modeling of lactose and other ions through the nanofiltration membranes is very important, because applying this modeling in experimental tests can lead to obtaining a product with higher purity. Some researchers have studied the nanofiltration process in situation where ions and uncharged solutes are present in the feed solution. They have predicted amount of solute retention by using the Donnan steric partitioning model (DSPM). Van Der Horst et al. [1] only described salt retention, and not lactose retention, when whey was used as feed. Cuartas-Uribe et al. [2, 3] used the DSPM and the Spiegler-Kedem model (SK) to analyze and predict the retention of only uncharged solutes (lactose) in the nanofiltration process.

The SK model is based on irreversible thermodynamics. This model, at first, was developed for reverse osmosis membranes and for single-solute nonelectrolyte solutions [3]. In this model the membrane is considered as a black box, and it is characterized in terms of solute permeability and the reflection coefficient.

Ahmad et al. [4] proposed extended Spiegler-Kedem (ESK) model by incorporating solute-solute interactions in the nanofiltration of multiple solutes systems. The global equations of ESK model are as follows: is the total volumetric flux (m³/m2 s). is the solute permeability coefficient of solute with the consideration of the interaction of solute . is the solute permeability coefficient of solute with the consideration of the interaction of solute . is reflection coefficient of solute . is the mass transfer coefficient of solute . is the observed rejection of solute , and is true rejection of solute where is the permeate concentration of solute . is the solute concentration in the bulk solution, and is the solute concentration on the membrane surface (kg/m³). is hydraulic permeability constant of membrane (m/Pas).   is transmembrane pressure (Pa), and is ideal gas constant (8.314 m³ Pa/mol K). is the temperature (K). is the molar mass of solute (g/mol), and is the total number of solutes. Vakili-Nezhaad and Akbari [5] have proposed a novel method for solving the ESK model by simulating the multiple solutes nanofiltration systems via genetic algorithm (GA). Their method is independent of the number of solutes and data points (experimental data of observed rejection () for each solute verses flux () are taken at different and constant feed rate), and it is able to predict the performance of multiple solutes systems with high precision. In the present work, the performance of UF-whey (ultrafiltered-whey) nanofiltration system was predicted by using the ESK model and GA method. The UF-whey nanofiltration system is a pilot plant equipped with flat and spiral wound membrane modules. The used membranes were Desal 5 DL and NF 200. Whey was ultrafiltered to avoid fouling problems in the nanofiltration process [2].

2. Results and Discussion

2.1. Feed Characterization

The results of analytical characterization of UF whey are shown in Table 1. It was observed by Cuartas-Uribe et al. that the UF whey did not contain fat and the amount of protein was low [2].

Table 1: Analytical characterization of the UF whey [2].

Solute
Concentration in UF whey (NF-Feed) Kg/m³

Cl⁻
1.29

Na⁺
0.37

K⁺
1.48

0.6

Mg²⁺
0.072

0.093

Lactose
41.9

Ca²⁺
0.289

2.2. Membrane Characterization

Cuartas-Uribe et al. used two aromatic polyamide composite membranes. These membranes were NF 200 and Desal 5 DL, supplied by Dow Chemical (USA) and GE Osmonics (USA), respectively [2].

2.3. Results of Solving ESK Model with GA Method

2.3.1. Genetic Algorithm

In order to find global optimization method, a more powerful method compared to the traditional ones is necessary. In the recent years genetic algorithm (GA) has been found to be an efficient algorithm for the complex engineering calculations. GA is based on the natural selection and the process that drives biological evolution. Those optimization problems that are not well suited for standard optimization algorithms, including problems in which the objective function is discontinuous, nondifferentiable stochastic, or highly nonlinear can be solved by GA. GA differs from a classical, derivative-based, optimization algorithm in two main aspects. First, classical algorithm generates a single point at each iteration and the sequence of points approaches an optimal solution while GA generates a population of points at each generation and the best point in the population approaches to an optimal solution. The second difference is that classical algorithm selects the next point in the sequence by a deterministic computation while GA selects the next population using random number of generators. GA modifies a population of individual solutions between the values of lower and upper bounds of variables successively, in which the initial range only restricts the range of the points in the initial population. At each step, GA selects individuals from the current population stochastically to be the parents of the children for the next generation based on genetic operators including selection, crossover, and mutation. The population size determines the size of population at each generation. Increasing the population size enables GA to search more points, thereby obtaining better results. Fixing the correct population size in a problem is an important factor for proper working of GA. Over a wide range of function optimization problems, a population size of 10 times of the number of variables is recommended. GA moves from generation to generation until the termination criterion is met. The most frequently used stopping criterion is specified maximum number of generations.

2.3.2. Solution of the ESK Model Using GA

The ESK model parameters (solute permeability (, ), solute reflection coefficient (), and mass transfer coefficient ()) were obtained by the analysis of real rejection and flux data of Cuartas-Uribe et al. [2, 3]. For this purpose the procedure proposed by Vakili-Nezhaad and Akbari [5] was considered (multiobjective optimization method in genetic algorithm by using Matlab 7.6.1 (R2008a), the software package). The results obtained for NF 200 membrane are summarized in Table 2 and for Desal 5 DL in Table 3. Figures 1 and 2 compare lactose and ions rejection predicted by the ESK model with the experimental data for two nanofiltration membranes. As can be observed, the results obtained with the ESK model are in a good agreement with the experimental results. In addition, the shape of the curve for lactose and ions rejection versus permeate flux predicted by the ESK model is more similar to the experimental results than the one obtained with the DSPM and SK model by Cuartas-Uribe et al. [3] for two membranes.

Table 2: Parameters estimated from solving ESK model with GA method for NF 200 membrane.

Parameter/s
Cl⁻(1)
Na⁺(2)
K⁺(3)
(4)
Mg⁺(5)
(6)
Lactose(7)
Ca²⁺(8)

0.8099
0.9994
0.85978
0.9908
0.9990
0.9950
0.9900
1.0000

5.81
1.51
9.90
9.00
9.00
9.50
8.51
9.50

6.86
9.04
2.01
2.00
3.00
5.00
1.01
8.01

1.01
1.20
2.01
2.00
8.91
5.01
2.00
9.00

1.26
1.00
9.55
2.01
7.61
5.01
2.00
2.00

1.27
1.01
2.01
2.07
5.51
5.00
2.01
2.01

3.58
1.06
2.00
2.01
9.00
9.00
2.00
2.01

1.40
1.03
2.01
2.01
6.90
7.01
4.01
2.01

1.20
5.06
2.01
2.01
1.01
1.51
9.99
1.00

1.29
5.01
2.01
2.00
3.01
2.51
2.01
5.06

Table 3: Parameters estimated from solving ESK model with GA method for Desal 5 DL membrane.

Parameter/s
Cl⁻(1)
Na⁺(2)
K⁺(3)
(4)
Mg⁺(5)
(6)
Lactose(7)
Ca²⁺(8)

0.999
0.6498
0.4802
0.9043
0.8912
0.9502
0.9428
1.0000

5.80
7.00
2.00
5.04
3.24
2.50
8.99
6.52

5.56
2.00
2.01
2.00
2.99
5.00
2.48
8.01

3.14
3.75
2.01
2.00
8.99
5.01
6.52
9.00

1.06
2.02
4.06
2.01
6.99
5.01
5.08
2.00

1.27
2.00
1.99
1.01
5.50
5.00
1.98
2.01

1.58
1.99
1.99
2.01
1.00
9.00
3.00
2.01

3.97
1.99
2.07
2.01
6.98
7.01
7.01
2.01

2.02
4.99
1.01
2.01
1.41
7.99
8.99
1.00

1.29
2.08
2.00
2.00
3.01
2.51
2.01
1.51

(a)

(b)

(c)

Figure 1: Ions and lactose rejection of UF whey versus volumetric flux for NF 200 membrane. Experimental results and ESK model by using GA method simulations. (a) Cl⁻, Na⁺, , and ; (b) Mg²⁺, K⁺, and Ca²⁺; (c) lactose.

(a)

(b)

(c)

Figure 2: Ions and lactose rejection of UF whey versus volumetric flux for Desal 5 DL membrane. Experimental results and ESK model by using GA method simulations. (a) Cl⁻, Na⁺, , and ; (b) Mg²⁺, K⁺, and Ca²⁺; (c) lactose.

3. Conclusions

The simulations which were made under the extended Spiegler-Kedem model and were solved by genetic algorithm method are reliable and suitable for predicting the performance of lactose and other solutes in UF-whey nanofiltration system. In addition, the model is able to predict lactose and other ions retention in presence or absence of ions without any information about the membrane structure. In comparison with similar works, this simulation is more perfect and outperforms the precision in the prediction of the performance of membranes in rejection and flux.

References

H. C. Van Der Horst, J. M. K. Timmer, T. Robbertsen, and J. Leenders, “Use of nanofiltration for concentration and demineralization in the dairy industry: model for mass transport,” Journal of Membrane Science, vol. 104, no. 3, pp. 205–218, 1995.

B. Cuartas-Uribe, M. I. Alcaina-Miranda, E. Soriano-Costa, and A. Bes-Piá, “Comparison of the behavior of two nanofiltration membranes for sweet whey demineralization,” Journal of Dairy Science, vol. 90, no. 3, pp. 1094–1101, 2007.

B. Cuartas-Uribe, M. C. Vincent-Vela, S. Álvarez-Blanco, M. I. Alcaina-Miranda, and E. Soriano-Costa, “Nanofiltration of sweet whey and prediction of lactose retention as a function of permeate flux using the Kedem-Spiegler and Donnan Steric Partioning models,” Separation and Purification Technology, vol. 56, no. 1, pp. 38–46, 2007.

A. L. Ahmad, M. F. Chong, and S. Bhatia, “Mathematical modeling and simulation of the multiple solutes system for nanofiltration process,” Journal of Membrane Science, vol. 253, no. 1-2, pp. 103–115, 2005.

G. Vakili-Nezhaad and Z. Akbari, “Modification of the extended Spiegler-kedem model for simulation of multiple solute systems in nanofiltration process,” Desalination and Water Treatment, vol. 27, no. 1–3, pp. 189–196, 2011.