 Open Access
 Authors : G. Reza VakiliNezhaad, Hamidreza Ziaiefar, Shivani B. Mishra, Hasan Mousa
 Paper ID : IJERTV9IS020305
 Volume & Issue : Volume 09, Issue 02 (February 2020)
 Published (First Online): 16032020
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
An Optimization based Solution for Extended SpieglerKedem Model using Grey Wolf Optimizer

Reza VakiliNezhaad
Department of Petroleum and Chemical Engineering, College of Engineering
Sultan Qaboos University
Shivani B. Mishra
Nanotechnology and Water Sustainability Unit, College of Science, Engineering and Technology
University of South Africa
Abstract Extended SpieglerKedem (ESK) is an irreversible thermodynamics model that has been originally developed to predict the nanofiltration performance of multiple solutes systems by considering the osmotic coefficient in a membrane surface for nonideal solutions. In the ESK model, the membrane
Hamidreza Ziaiefar
Department of Petroleum and Chemical Engineering, College of Engineering
Sultan Qaboos University
Hasan Mousa
Department of Chemical Engineering Jordan University of Science and Technology
the model, Ahmad et al. developed SK for multiple solute systems by considering the solutesolute interactions. This model is named Extended SpieglerKedem (SKE) and can be summarized as the following equations:
is considered as a black box and all of the unknown parameters
=
( ) (1)
must be found by experimental data. Due to difficulties of doing enough precise experiments, traditional mathematical solutions
=1
such as LevenbergMarquardt are not applicable. In this paper,
= [(1) (1 +
)] (2)
at the first step, a modified ESK model for prediction of a high concentration solution is presented. Later, to compensate for the
=1
weakness of traditional mathematical methods, a solution based on optimization algorithms is proposed. Two optimization
=
( ) exp ( )
(3)
algorithms, Genetic Algorithm (as one of the most common algorithms) and Grey Wolf Optimizer (one of the most recent
[ (1 ) ]
(1 )
algorithms), are chosen and the results of ESK solution by them are provided and compared.
=
1
1
exp ( ) (4)
Keywords Nanofiltration; SpieglerKedem; Negative Rejection; Grey Wolf Optimizer; Genetic Algorithm.
= 1
( ) (5)

INTRODUCTION
The process of filtration through a membrane (reverse and forward osmosis) could be explained by the transport mechanism or irreversible thermodynamics. Solution diffusion, KimuraSourirajan and extended NerstPlank are the models based on transport mechanism and the Spiegler Kedem (SK) and KedemKatchalsky are based on irreversible thermodynamics. In irreversible thermodynamics, the membrane is considered as a black box, without considering the structure of membrane or mechanisms of transport. To simplify the modeling, the whole filtration is also considered as a slow process near the equilibrium. In this method of modeling, the flux of solvent and solute are only proportional to the chemical potential which is derived by pressure and concentration gradient.
Among these different methods, SK model has been widely used to predict the performance of single solutes in a Nanofiltration, reverse osmosis[13] and when one of the ions is impermeable to the membrane[4]. SK does not consider the interactions between the solutes and therefore, in many real practical filtrations, like water desalination, which the solution consists of more than two solutes, SK is not useful. To improve
where F is driving force (kWs/m mol), A is osmotic constant (m3Pa/g), Jv is total volumetric flux (m3/m2 s), C is concentration, is reflection coefficient and Ro is observed rejection. In all equations, the subscripts of b, p, s indicates bulk, permeate and i solute respectively.
In the ESK model, the solutes are considered ideal. In a non ideal solution with a concentration of more than 1 molar on the membrane surface, the model must be modified using:
= ( ) (6)
where is the osmotic coefficient and k is the total number of constituent ions in the salt[5]. In this equation, by considering of the concentration on the solution on the membrane, , can be calculated for nonideal solutions.
Ge et al. proposed a threeparameter model to calculate the osmotic coefficient[6]. In their model, the ionion distance, ionsolvent parameter, and solvent parameter are considered. In Ges model, the two parameters of ionion and ionsolvent molecule interaction are not dependent to of the temperature of the solution and solvent[7]. The osmotic coefficient can be obtained from:
= 1
+
1
2
[ 1
1
1 1
2 (1 +
1
1 1
22)
]
problem[8]. The better algorithm can find these optimum values in less time to save time and energy. Nowadays, there
are different optimization algorithms[911]. From a well
2 2 2 2
1 2 2
()
2
+
2
2
known algorithm like Genetic Algorithm (GA) [12] to the most recent ones[1315]. In this paper, GA as a complete
+ 2 + 1 + (7)
famous algorithm in literature and a modern optimization
+
algorithm named Grey Wolf Optimizer (GWO) are implemented. Due to the popularity of GA between
where is the osmotic factor, z is the ion charge, A and B are 6064.613 1/21/23/2 and 56.827 1/21/23/2 respectively. T, a, n and S show temperature, the parameter of ionion distance, the parameter of ionsolvent and salvation parameter respectively.
In equation 7, is the stoichiometric number and can be calculated as:
+ +
+ +
+ =   (8)
I, ionic strength, can be calculated as
1
researchers, there are many papers in the literature that used and discussed GA[16]. Therefore, in this paper GA principles and equations are not discussed and only the results and performance comparison between GA and GWO are provided.
III. GREY WOLF OPTIMIZER
GWO is an intelligent swarmbased algorithm developed by Mirjalili in 2014 [17]. The algorithm is inspired by the natural hunting behavior of the Grey wolf in nature as an apex predator on the top of the food pyramid. The algorithm simulates hunting (optimization of the cost function) with male and female wolves which are followed by second and third top hunters (best solutions) named and . The model of
=
. 2
2
(9)
the optimization consists of searching for prey, encircling, hunting, and attacking.

Encircling prey
where m is the molality (mol/kg).


FINDING UNKNOWN PARAMETERS
As it has been mentioned before, in the ESK model, the membrane is considered as a blackbox. This means that there is a need for experimental data to find the unknown parameters. However, it should be considered that the number of unknown parameters in this model is relatively high. For
These steps are mathematically described as:
=  . () () (10)
( + 1) = () . (11)
where t indicates the current iteration, V and W are coefficient vectors, is the position vector of the prey, and X indicates the position vector of a grey wolf. The vectors V and W are
calculated as follows:
instance, to predict the erformance of the membrane in the foursolute solution, the unknown parameters of the model
= 2 . 1
(12)
such as , ks, Pss, Psi has to be found. The term Pss is the solute permeability coefficient of solute s with the consideration of the interaction of solute s and the Psi is the solute permeability coefficient of solute s with the consideration of the interaction of solute i, ss is reflection coefficient and ks is the mass transfer coefficient. In total, by considering n number of solutes in a system, the number of unknown variables will be 2 + 2. By considering p number of experimental data, we will have
Ã— equation that can be used to predict unknown variables. The mathematic solutions for such a problem like Levenberg Marquardt, only work when + 2. However, optimization algorithms can solve this equation. Optimization algorithms test different random values for all of the unknown variables and check the output of the model outputs with the experimental data. A set of random number for variables with
= 2 2 (13)
where components of v are linearly decreased from 2 to 0 over the course of iterations and r1, r2 are random vectors in [0,1].

Hunting mechanism
In the hunting process, alpha, beta, and delta wolves lead the others to the prey position. In the mathematical simulation, we consider that alpha (the current best solution), beta and delta (the second and third top solutions) know the optimum solution (prey location). Therefore, the top best three solutions are saved and considered as the best potential of the optimum solution and the other search agents are obliged to update their position based on them. These equations simulate the search process.
the minimum acceptable difference between real condition and
=  .
 (14)
modeling can be considered as the correct answer. It should be considered that the optimization algorithms do not provide a certain answer for a specific question. They just try to be as close as possible to an optimum answer. This optimum finding is limited by processing time and resources. As much as process power increases or the period of the process be higher, the found answer will be closer to the optimum values. In a practical problem, the accepted closeness of the answer must be determined based on the nature and importance of the
1
=  2.  (15)
=  3.  (16)
1 =  () 1. ( ) (17)
2 =  () 2. ( ) (18)
3 =  () 3. ( ) (19)
and rest of the wolves update their position as the following equation:
TABLE II. THE ESTIMATED FACTOR OF ESK MODEL BY GWO
( + 1) =
1
+ 2
3
+ 3
Factor
Na+
NO3
`Cl
Ca2+
0.9889
0.7489
0.8159
0.8485
4.13e05
7.27e05
4.45e3
6.27e3
1
1.53e05
9.31e08
8.85e08
4.18e08
2
1.26e07
1.20e05
1.01e05
4.10e08
3
2.45e06
1.21e06
1.76e06
8.65e08
4
6.45e07
6.2708
6.5708
1.27e06
R2
0.9623
0.9010
0.9236
0.8277
Factor
Na+
NO3
`Cl
Ca2+
0.9889
0.7489
0.8159
0.8485
4.13e05
7.27e05
4.45e3
6.27e3
1
1.53e05
9.31e08
8.85e08
4.18e08
2
1.26e07
1.20e05
1.01e05
4.10e08
3
2.45e06
1.21e06
1.76e06
8.65e08
4
6.45e07
6.2708
6.5708
1.27e06
R2
0.9623
0.9010
0.9236
0.8277
(20)

Searching for and attacking the prey
Searching (exploration) and attacking (exploitation) are two abilities that grey wolves use in hunting. Grey wolves search for prey by diverging and attack prey by converging. Diverge and converge can be simulated mathematically by using factor. is a random number in the interval of [2 , 2].
  1 forces the search agents toward the current optimum
points and   1 forces them to spread for fittest value. Exploration and exploitation make sure that the algorithm moves toward the optimum values and dont be trapped in a local optimum value[18].

RESULTS AND DISCUSSION
The ESK model is tested with experimental data of Wadley et al. [19]. The data consists of the laboratory scale Nano filtration experiments of the solution consists of 111 ternary mixture of NaNO3, NaCl, CaCl2. The total concentration of feeding solution is 0.0015 kg/m3. In total, system contains of Na+ of 0.69 kg/m3, NO3 with 0.93 kg/m3, Cl with concentration of 1.5957 kg/m3 and Ca2+ with the concentration of 0.6 kg/m3. The membrane used as a filter is NF200 and all the tests are done at 35 C with the range of pressure from 2.0 to 3.5 MPa.
The unknown parameters of the model with four solutes,
, , , have to be found by the optimization algorithm. To have a better understanding of the performance of each algorithm in finding the optimum values of unknown parameters, both algorithms are run in Matlab 2017, by a 2670 core i7Â® CPU without pooling process. Both algorithms used the same cost function and the tests are done in the same number generation or process time. However, based on the fundamental differences between the two algorithms, process time and convergence rate are so different for them. While each generation of GA needs a shorter period of time than GWO, the overall trend of convergence in GWO is much faster. To have a better understanding of the performance in solving ESK, both of them are run in an equal period of time [7].
TABLE I. THE ESTIMATED FACTOR OF ESK MODEL BY GA
Factor
Na+
No3
Cl
Ca2+
0.9886
0.7499
0.80877
0.8597
4.12e05
7.26e5
4.44e3
3.62e3
1
1.54e5
9.22e8
8.77e8
4.19e8
2
1.17e7
1.20e5
1.02e5
6.81e8
3
1.96e6
1.21e6
1.79e6
5.60e8
4
7.51e7
6.29e8
6.63e8
9.56e7
R2
0.9423
0.9030
0.8779
0.5374
The results of the optimization show a significant improvement in the coefficient of determination (R2) for Ca2+ and an overall improvement for the others. Fig. 1 shows the ion rejection of the system versus volumetric flux for all four ons. As it can be seen in the figure, the Na+ shows negative rejection which the models and the algorithms could predict completely. The mechanism of this rejection can be explained by the effect of Donnan. Based on this effect, the negative rejection of Na is because of the negative charges of the surface of the membrane[2025].
Fig. 1. Ions rejection of system versus volumetric flux (GWO).

CONCLUSION


SpieglerKedem model is a useful method to predict the behavior of single solutions in the reverse, forward osmosis, and Nanofiltration. An extended version of SpieglerKedem (ESK) can be implemented to estimate the performance of multiple solutes in the filtration. In the ESK model, the membrane is considered as a blackbox and unknown parameters can be found by mathematical calculation over experimental data. Solving the equations using traditional mathematical solutions like LevenbergMarquardt needs an adequate number of experimental data points. In the solutions with a high number of solvents, having an adequate number of experimental points is not easy to access. In this paper, we present an optimization solution to find the unknown parameters of the ESK model. Genetic Algorithm (one of the most common) and Grey Wolf Optimizer (one of the most recent ones) are implemented and the results of the solution with them are compared and presented. The results show that while GA and GWO are both able to move toward the better
answers for the equations, GWO shows better improvement in the accuracy of the answers. Lower process time of this algorithm makes it a more desirable option for finding the solutions in such problems.
ACKNOWLEDGEMENT
The authors appreciate the financial support of Sultan Qaboos University and University of South Africa through the joint research funding (CL/SQUSA/18/04).
REFERENCES

S. Jain and S. K. Gupta, "Analysis of modified surface force pore flow model with concentration polarization and comparison with SpieglerKedem model in reverse osmosis systems," Journal of Membrane Science, vol. 232, no. 12, pp. 4562, 2004.

D. Van Gauwbergen and J. Baeyens, "Modelling reverse osmosis by irreversible thermodynamics," Separation and purification technology, vol. 13, no. 2, pp. 117128, 1998.

I. Sutzkover, D. Hasson, and R. Semiat, "Simple technique for measuring the concentration polarization level in a reverse osmosis system," Desalination, vol. 131, no. 13, pp. 117127, 2000.

A. Hassan et al., "A new approach to membrane and thermal seawater desalination processes using nanofiltration membranes (Part 1)," Desalination, vol. 118, no. 13, pp. 3551, 1998.

V. K. Gupta, S.T. Hwang, W. B. Krantz, and A. R. Greenberg, "Characterization of nanofiltration and reverse osmosis membrane performance for aqueous salt solutions using irreversible thermodynamics," Desalination, vol. 208, no. 13, pp. 118, 2007.

X. Ge, X. Wang, M. Zhang, and S. Seetharaman, "A new three particleinteraction model to predict the thermodynamic properties of different electrolytes," The Journal of Chemical Thermodynamics, vol. 39, no. 4, pp. 602612, 2007.

G. VakiliNezhaad and Z. Akbari, "Modification of the extended SpieglerKedem model for simulation of multiple solute systems in nanofiltration process," Desalination and Water Treatment, vol. 27, no. 13, pp. 189196, 2011.

K. A. Hatmi, A. a. A. Mashrafi, S. U. A. Balushi, H. AlKalbani,
M. Shaikh, and G. R. VakiliNezhaad, "Ways for Optimizing and Increasing Reliability of Centrifugal Compression Systems in Oil and Gas Industry: Case Study in Sultanate of Oman," in Abu Dhabi International Petroleum Exhibition & Conference, 2017: Society of Petroleum Engineers.

A. A. Esmin, R. A. Coelho, and S. Matwin, "A review on particle swarm optimization algorithm and its variants to clustering high
dimensional data," Artificial Intelligence Review, vol. 44, no. 1, pp. 2345, 2015.

S. SalcedoSanz, "A review on the coral reefs optimization algorithm: new development lines and current applications," Progress in Artificial Intelligence, vol. 6, no. 1, pp. 115, 2017.

B. Shahriari, K. Swersky, Z. Wang, R. P. Adams, and N. De Freitas, "Taking the human out of the loop: A review of Bayesian optimization," Proceedings of the IEEE, vol. 104, no. 1, pp. 148 175, 2015.

J. R. Koza, "Genetic programming," 1997.

S. Mirjalili, "The ant lion optimizer," Advances in Engineering Software, vol. 83, pp. 8098, 2015.

S. Mirjalili and A. Lewis, "The whale optimization algorithm,"
Advances in engineering software, vol. 95, pp. 5167, 2016.

S. Mirjalili, "Mothflame optimization algorithm: A novel nature inspired heuristic paradigm," KnowledgeBased Systems, vol. 89, pp. 228249, 2015.

L. Davis, "Handbook of genetic algorithms," 1991.

S. Mirjalili, S. M. Mirjalili, and A. Lewis, "Grey wolf optimizer,"
Advances in engineering software, vol. 69, pp. 4661, 2014.

S. Mirjalili, S. Saremi, S. M. Mirjalili, and L. d. S. Coelho, "Multi objective grey wolf optimizer: a novel algorithm for multicriterion optimization," Expert Systems with Applications, vol. 47, pp. 106 119, 2016.

S. Wadley, C. Brouckaert, L. Baddock, and C. Buckley, "Modelling of nanofiltration applied to the recovery of salt from waste brine at a sugar decolourisation plant," Journal of Membrane Science, vol. 102, pp. 163175, 1995.

X.L. Wang, T. Tsuru, S.i. Nakao, and S. Kimura, "The electrostatic and sterichindrance model for the transport of charged solutes through nanofiltration membranes," Journal of Membrane Science, vol. 135, no. 1, pp. 1932, 1997.

A. E. Yaroshchuk, "Negative rejection of ions in pressuredriven membrane processes," Advances in colloid and interface science, vol. 139, no. 12, pp. 150173, 2008.

J. Gilron, N. Gara, and O. Kedem, "Experimental analysis of negative salt rejection in nanofiltration membranes," Journal of Membrane science, vol. 185, no. 2, pp. 223236, 2001.

S. Mandale and M. Jones, "Interaction of electrolytes and non electrolytes in nanofiltration," Desalination, vol. 219, no. 13, pp. 262271, 2008.

I. Koyuncu and D. Topacik, "Effect of organic ion on the separation of salts by nanofiltration membranes," Journal of Membrane Science, vol. 195, no. 2, pp. 247263, 2002.

A. Yaroshchuk, "Asymptotic behaviour in the pressuredriven separations of ions of different mobilities in charged porous membranes," Journal of Membrane Science, vol. 167, no. 2, pp. 163185, 2000.