Issue 
Int. J. Simul. Multisci. Des. Optim.
Volume 6, 2015



Article Number  A4  
Number of page(s)  8  
DOI  https://doi.org/10.1051/smdo/2015004  
Published online  29 September 2015 
Research Article
Nonlinear observer to estimate polarization phenomenon in membrane distillation
Process Control Laboratory, Control Engineering Department, Ecole Nationale Polytechnique, 10 avenue Pasteur, HassanBadi, 16200
Algiers, Algeria
^{*} email: khoukhi.billal@gmail.com
Received:
3
April
2015
Accepted:
21
August
2015
This paper presents a bidimensional dynamic model of Direct Contact Membrane Desalination (DCMD) process. Most of the MD configuration processes have been modeled as steadystate onedimensional systems. Stationary twodimensional MD models have been considered only in very few studies. In this work, a dynamic model of a DCMD process is developed. The model is implemented using Matlab/Simulink environment. Numerical simulations are conducted for different operational parameters at the module inlets such as the feed and permeate temperature or feed and permeate flow rate. The results are compared with experimental data published in the literature. The work presents also a feed forward control that compensates the possible decrease of the temperature gradient by increasing the flow rate. This work also deals with a development of nonlinear observer to estimate temperature polarization inside the membrane. The observer gives a good profile and longitudinal temperature estimations and shows a good prediction of pure water flux production.
Key words: Direct contact membrane distillation / Dynamic modeling / Heat and mass transfer / Unknown input observer / Polarization coefficient / Saline water desalination
© B. Khoukhi et al., Published by EDP Sciences, 2015
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Membrane distillation (MD) process is an emerging technology for water treatment. The driving force of the MD process is given by the pressure difference of vapor formed by a difference in temperature of solutions on both sides of a hydrophobic membrane [1]. The advantages of Direct Contact Membrane Desalination (DCMD) lie in its simplicity, the need for only small temperature differences and nearly 100% rejection of dissolved solids [1]. Furthermore, the low energy demand systems in DCMD processes can be equipped with renewable energy equipment such as solar collectors [2] and solar distillers [3]. Many MD configuration processes have been modeled as steadystate onedimensional systems using empirical heat and mass transfer equations [4]. Only few publications use stationary one or twodimensional heattransfer equations to simulate the process more accurately. Although many semiempirical models have been developed, a detailed model for temperature polarization on flatplate MD processes is still lacking. Some studies were interested in a theoretical modeling and experimental analysis of direct contact membrane distillation in steady state such as in [4]. A dynamic modeling of direct contact membrane distillation processes has also been presented in [5, 6]. In another way, to show the interest of using renewable energies in DCMD processes, a performance investigation of a solarassisted direct contact membrane distillation system was conducted in [7].
Our study focuses on a bi dimensional dynamic model to simulate the membrane temperature and the pure water flux. It proposes an algorithm to resolve the temperature partial differential equations (PDE) that describes the process.
Our study also uses a nonlinear observer in order to estimate all temperatures and fluxes inside the membrane from the measurement of the accessible data of the membrane, which are inlet and outlet temperatures.
Because temperature inside the membrane is not accessible for measurement this observer is very useful and can be considered as a software sensor in order to estimate temperature distribution.
Several approaches have been tried to extend the Unknown Input Observer (UIO) design for linear systems to nonlinear systems. The first works have been dedicated to UIOs for bilinear systems [8]. Other approaches based on the transformation of the nonlinear system into canonical forms are proposed in [9] and [10]. The limitation of these approaches is that the required state transformation does not exist for all nonlinear systems but only for a limited class. The observer developed in this paper is designed in a cascade structure and is specific to the model of DCMD.
This paper is organized as follow: in the next section a model of the DCMD is presented. In Section 3 a nonlinear observer is developed. In Section 4 a presentation of the simulation results is given. Finally, the paper close with some concluding remarks.
2. Modeling the membrane distillation process
In Direct Contact Membrane Distillation (DCMD) both sides of the membrane are in direct contact with a liquid stream. On the left of the membrane shown in Figure 1, the hot liquid (i.e. hot seawater) flows in the evaporator channel, whilst on the right a cold liquid (i.e. cooled permeate or distillate) is circulated. Heat transfer (as well as mass transfer) occurs from the hotter to the colder side. The liquid in the evaporator channel is constantly refilled and reheated, whilst the volume of the liquid in the permeate channel increases and heats up. One of the main features of DCMD is that the gas gap between the membrane surface and the condensate stream is very narrow and only exists due to the hydrophobic nature of the membrane. This causes the temperature of the membrane surface in contact with the condensate to be very close to that of the condensate stream itself, thus allowing high temperature drops across the membrane, i.e. high driving forces for mass transfer. Conversely, the Direct Contact configuration causes a relatively high heat loss, as the membrane is the only barrier for the transfer of sensible heat [2]. In the following, we will give the mathematical equations that describe those phenomena.
Figure 1. Schematic diagram of membrane distillation process. 
2.1. Pure water flux in the distillation process
The mass transfer driving force across the membrane is the difference in saturated pressure components on both membrane surfaces due to the temperature gradient. The general mass flux form can be expressed as follows:(1)where and , are the saturated pressure of water on the hot and cold feed membrane surfaces, respectively. For nonideal binary mixtures, the flux can be determined as [4, 11]:(2)where x_{NaCl} is the mole fraction of NaCl in saline solution. Saturated pressures can be determined by the Antoine equation [12]:where T_{ k } is the temperature in °C. The membrane coefficient c_{ m } in (1) can be estimated by a weighted sum (via parameters α(T) and β(T)) of the Knudsen diffusion and the Poiseuille (viscous) flow models [13], i.e.(3)where, M_{ w } is the molecular weight of water, P_{ m } is the mean saturated pressure in membrane, R is the gas constant, r is the pore radius of the membrane, T_{ m } is the mean temperature in membrane, δ_{ m } is the thickness of membrane, ε is the porosity of membrane, η_{ v } is the gas viscosity and τ is the tortuosity factor. The tortuosity of a porous hydrophobic membrane was estimated by [14].
The objective of the modeling is to estimate the temperature distribution in the flow channels and use the equation (2) to calculate the pure water flux.
2.2. The temperature equations in the flow channel
Figure 1 shows a schematic diagram of the parallelflow membrane distillation system. With a hydrophobic, micro porous membrane, liquid water is prevented from entering the pores, while molecular water in the vapor phase can pass through. In order to obtain accurate temperature distribution on membrane surfaces with position, a study of heat transfer and transport phenomena should be done [15, 16].
The threedimensional governing dynamic equation of energy is accurate for temperature distribution in flow channel but the calculation procedure is complex. The main benefit of using a two dimensional computation is saving time and resources cost, but this advantage is followed by losing the information related to scaling the problem from three dimensions to two dimensions.
In order to obtain the approximation temperature distribution on the both channels of the membrane distillation system, the energy equations were simplified with the following assumptions: (1) laminar flow; (2) symmetrical flow and temperature distribution; (3) no internal generation of energy; (4) both zaxis and xaxis directions of velocity exists but only xaxis velocity is explicitly taken into account, the zaxis velocity will be taken into account implicitly as will be demonstrated in resolution algorithm. The equations of energy thus may be obtained as:(4)
The first step to resolve this PDE (Eq. (4)) is to divide the flow channel into M sub systems in cascade so that the variation of the temperature through the zaxe on each sub system is equal to zero. Then, by resolving each sub system separately the summation of flux obtained will give the total pure water flux. The first membrane boundary limit is the inlet flow temperature. The temperatures relaying between the subs systems are the mean temperatures profiles of the previous sub system. The other boundary limits are those in the literature and will be defined in the following. The equation of each sub system is:(5)
This PDE is solved using the finite elements and we can find:(6)
The boundary limits used are those of the heat transfer mechanism [17]:(8)with:
We obtain then, a matrix (M × N) of ordinary differential equations (ODE). M is the number of sub system, which the membrane is divided in cascade in the direction of zaxis. N is the number of equation, which each sub system is discretized in space in the xaxis using finite elements. The diagram in Figure 2 explains the schematic of resolution.
Figure 2. Schematic diagram resolution algorithm. 
Let put T_{a,i }(t) = ψ_{i}^{hot} and T_{b,i}(t) = ψ_{i}^{cold}, Σ_{i} can then be described by equations in stat space. Let’s take:(9)
The output of each sub system is the mean profile temperature and it is also the input of the next sub system. We obtain the stat space equation of the sub system Σ_{i} in the hot side:(10)and(11)and in the cold side we obtain:(12)and(13)with i = 1, M and u_{1} = T_{in} (in the hot and cold side separately).
When seen as a complete system we can write:and(14)where J_{i} is the pure water flux in each sub system.
Figures 3 and 4 show the bi dimensional result of temperature equation in the flow channels of membrane. As can be seen, the temperature in the hot side decrease along the x and zaxis of membrane, in the same way the cold side temperature increase along the x and the zaxis, this phenomenon is called polarization.
Figure 3. Schematic diagram of hot flow temperature in two dimensions coordinates. 
Figure 4. Schematic diagram of cold flow temperature in two dimensions coordinates. 
A case study from literature has been investigated. The data used in our simulations like geometry, physical properties and operating conditions are the same used in [4]. Simulations have been conducted for different inlets temperature. Compared with the stationary experimental data in [4].
As shown in Figure 5, the computed mass flux densities are in good agreement with the measurements.
Figure 5. Pure water flux function of temperature gradient. 
3. Nonlinear observer
In practice, only the inlet and outlet temperatures are accessible for measurement. The need of nonlinear observer arises. The observer will estimate all temperatures inside the membrane and from those temperatures it will be possible to estimate polarization coefficient, polarization ratio, and pure water flux. Because the structure in cascade of the model developed in the previous section, a nonlinear unknown input observer (NUIO) is used. NUIO are used for different systems, a design of unknown input functional observers for nonlinear systems with application to fault diagnosis has been presented in [18]. A nonlinear unknown input observer design for diesel engines is described in [19]. In [20] an observerbased residual design for nonlinear systems with unknown inputs is developed. The structure of the NUIO in our study is in cascade. The advantage of this kind of structure is that we can increase the accuracy of the model by increasing the number of subsystems describing the membrane. By Applying the NUIO to the last subsystem we can estimate all temperatures inside it. We can also estimate the unknown input of these sub system, which is the output of the previous subsystem. Using this output we can estimate the temperatures inside the previous subsystem. Using this method, we can ascend to the first subsystem, which has his input accessible for measurement (inlet temperature). The following scheme (Figure 6) describes the observer used.
Figure 6. Schematic diagram the nonlinear unknown input observer (NUIO) in DCMD process. 
This structure has a lot of advantages; it is possible to increase the accuracy of the model by increasing the number of subsystems. Also, when estimating all temperatures inside the membrane, we will be able to calculate the estimated flux in each part of the membrane, calculate the total flux, and then, calculate different parameters such as polarization ratio and polarization coefficient.
3.1. Development
The NUIO is an observer based on a model, the structure of the hot side observer of the subsystem i is:(15)and y_{i} is the outlet temperature, which is accessible for measurement. Because of the structure of the model in cascade, the estimated unknown input is the output of the previous subsystem this output can be used to estimate all parameters of the previous subsystem (i − 1) and by continuing with same principal to the first subsystem we can estimate all the parameters of the membrane. The same development is applied to the cold side:(16)
4. Simulations results
4.1. Estimation of the pure water flux
The pure water flux is accessible for measurement; it is a good parameter to evaluate the accuracy of the observer by plotting the estimation error between the measured flow and the estimated flow given by the NUIO. The estimated pure water flux is given by the following equation (17):(17)with
By introducing a step in hot water inlet (Figure 7), the pure water flux changes, this is due to the predominance of the temperature gradient in pure water production. The observer follows this changes and the estimating error is given in Figure 8.
Figure 7. Step of hot inlet temperature introduced at t = 400 s. 
Figure 8. Pure water flux production estimation error. 
We can see that the observer gives a good estimation of the pure water flux, and can be used for fault diagnosis to detect abnormal pure water production.
4.2. Estimating polarization ratio
Because temperature polarization is the dominant effect on the pure water flux through membrane, the analysis of temperature polarization phenomenon is presented in Figure 9.
Figure 9. Polarization ratio in the longitudinal axe of the membrane. 
As seen from Figure 9, results show that the temperature polarization ratio, which is defined as T_{a,N+1}/T_{a,0} or T_{b,N+1}/T_{b,0}, increases as the hot stream proceeds in the flowing direction. This increase is due to decrease of the heattransfer coefficient in the flow channel, resulting in increasing both the velocity boundary layer and thermal boundary layer.
4.3. Estimation of the polarization coefficient
Temperature polarization characterizes the temperature distribution between bulk stream (T_{a,0} – T_{b,0}) and membrane interface temperatures, (T_{a,N+1} – T_{b,N+1}), which mainly depends on the ratio of heat transfer in the flow channels and through the membrane. The coefficient characterizing this relation is the temperature polarization coefficient:(18)
It can be seen that, for a very efficient heat transfer between bulk streams and membrane interfaces (T_{a,N+1} ≈ T_{a,0} and T_{b,N+1} ≈ T_{b,0}), the coefficient τ approaches 1. In the case of low heat transfer coefficients in the channels, compared to heat transfer coefficient through the membrane (T_{a,N+1} ≈ T_{b,N+1}), the polarization coefficient τ approaches 0. In MD systems, τ mainly depends on channel setup, module design and operation. The temperature polarization coefficient increases with an increasing heat transfer coefficient, conversely, an increase of the mean temperature, which leads to an increase of vapor transport can also lead to significant reductions in temperature polarization coefficients [21]. Values below 0.1 for flat membranes in laminar flow operation, and above 0.9 for narrow hollow fiber membranes in turbulent flow configuration have been reported [22]. The Figure 10 shows the temperature polarization coefficient estimating error inside the membrane at precise longitudinal position.
Figure 10. Polarization coefficient estimating error inside DCMD. 
We can see also from Figure 10 that the observer gives a good estimation of the polarization coefficient compared with the estimation given by the model. The temperature polarization coefficient must be taken into account in order to calculate the mass transfer from the bulk temperatures, which are usually well known, compared to the membrane interface temperatures, which are almost immeasurable. It can be seen that temperature polarization has a significant influence on flux and on the overall process efficiency, and therefore should be the focus of the MD process design.
4.4. Estimation of the inlet temperature
The observer developed is an observer with an unknown input observer, and the input is estimated using a model. Because of the structure in cascade of the observer, we have to estimate the inputs of all subsystem. The model used in our study is given in the following:

with α > 1 for the hot side and,

with α < 1 for the cold side.
Because the first input (hot and cold inlet temperature) is accessible for measurement, it is important to validate the model of the input by plotting the estimation error (Figure 11) between the estimated input and the measured input.
Figure 11. Hot inlet estimating error. 
We can see that the first input is estimated correctly, and this input is the most difficult to estimate, because the estimation error is cumulative. For accuracy reasons, the first subsystem of the observer is fed with the real measured input.
4.5. Estimation of the longitudinal temperature distribution
The twodimensional model is useful to study the transport mechanism across the micro porous membrane in detail. The membrane surface temperature distributions in the flowing direction are shown in Figure 12. Hot side decrease along flowing direction while the membrane surface temperature of cold side, increase with increasing longitudinal coordinate. Furthermore, the trans membrane pressure affected by the temperature gradient between both membrane surfaces taper along the flowing direction, resulting in decreasing the pure water productivity per unit area of membrane. Accordingly, the pure water flux decreases along the flowing direction due to the temperature difference decreasing. Though increasing the volumetric flow rate is helpful to enhance the productivity of pure water as confirmed from Figure 12.
Figure 12. Schematic diagram of temperature along z axes. 
A comparison with the results given by the observer and those given by the model is shown in Figure 13.
Figure 13. Longitudinal temperature estimating error along DCMD. 
The results show that the estimation error increases in the opposite of the flow direction. This is due to the configuration of the observer, which starts estimations from the last subsystem. The estimation error continues to increase until the introduction of the real measurement of the inlet temperature. The estimation error then decrease until the first input estimation which is measurable.
5. Conclusion
In this paper, a dynamic bi dimensional model has been proposed, a resolution algorithm has been proposed to solve the PDE and simulation results are given. The model is implemented using Matlab/Simulink. A nonlinear observer is used to estimate the major phenomenon inside the membrane. Results show that the model proposed gives good prediction of pure water production. Also, the observer gives a good prediction of the temperatures and flux inside the membrane. The observer can be used to monitor the membrane and the results can be used for fault diagnostic, which is a major step in this kind of MD process. We will continue our efforts to build an experimental setup so that we can perform experiments on our model and observer.
References
 Lawson KW, Lloyd DR. 1996. Membrane distillation. II. Direct contact MD. Journal of Membrane Science, 120, 123–133. [CrossRef] [Google Scholar]
 Koschikowski J, Wieghaus M, Rommel M. 2003. Solar thermaldriven desalination plants based on membrane distillation. Desalination, 156, 295–304. [CrossRef] [Google Scholar]
 de Andrés MC, Doria J, Khayet M, Pena L, Mengual JI. 1998. Coupling of a membrane distillation module to a multieffect distiller for pure water production. Desalination, 115, 71–81. [CrossRef] [Google Scholar]
 Chen T, Ho C, Yeh H. 2009. Theoretical modeling and experimental analysis of direct contact membrane distillation. Journal of Membrane Science, 3(30), 279–287. [CrossRef] [Google Scholar]
 Ashoor BB, Fath H, Marquardt W, Mhamdi A. 2012. Dynamic modeling of direct contact membrane distillation processes. Proceedings of the 11th International Symposium on Process Systems Engineering, 15–19 July, Singapore. [Google Scholar]
 de la Calle A, Bonilla J, Roca L, Palenzuela P. 2014. Dynamic modeling and performance of the first cell of a multieffect distillation plant. Applied Thermal Engineering, 70(1), 410–420. [CrossRef] [Google Scholar]
 Kim YD, Thu K, Ghaffour N, Ng KC. 2013. Performance investigation of a solarassisted direct contact membrane distillation system. Journal of Membrane Science, 427, 345–364. [CrossRef] [Google Scholar]
 Saif M. 1993. A disturbance accommodating estimator for bilinearsystems. IEEE American Control Conference (ACC 93), June 2–4, San Francisco, CA, USA. [Google Scholar]
 Frank PM, Ding X, Guo L. 1990. Nonlinear observer design via extended observer canonical form. Systems and Control Letters, 15, 313–322. [CrossRef] [MathSciNet] [Google Scholar]
 Seliger R, Frank PM. 1991. Robust component fault detection and isolation in nonlinear dynamic system using nonlinear observers, Proceedings of the IFAC/IMACS Symposium SAFE PROCESS, BadenBaden, Germany. [Google Scholar]
 Lawson KW, Lloyd DR. 1997. Membrane distillation. Journal of Membrane Science, 124, 1–25. [CrossRef] [Google Scholar]
 Felder RM, Rousseau RW. 2000. Elementary principles of chemical processes, 3rd ed., John Wiley & Sons: New York. [Google Scholar]
 Schofield RW, Fane AG, Fell CGD. 1990. Gas and vapor transport through microporous membrane. Journal of Membrane Science, 53–159. [Google Scholar]
 Iversen SB, Bhatia VK, DamJphasen K, Jonsson G. 1997. Characterization of microporous membranes for use in membrane contactors. Journal of Membrane Science, 130, 205–217. [CrossRef] [Google Scholar]
 Incropera FP, Dewitt DP. 1996. Fundamentals of heat and mass transfer, 4th ed., John Wiley & Sons: New York. [Google Scholar]
 Bird RB, Stewart WE, Lightfoot EN. 2007. Transport Phenomena, 2nd ed., John Wiley & Sons: New York. [Google Scholar]
 Alklaibi AM, Lior N. 2004. Membranedistillation desalination: status and potential. Desalination, 171, 111–131. [CrossRef] [Google Scholar]
 Teh PS, Trinh H. 2013. Design of unknown input functional observers for nonlinear systems with application to fault diagnosis. Journal of Process Control, 23, 1169–1184. [CrossRef] [Google Scholar]
 Boulkroune B, Djemili I, Aitouche A, Cocquempot V. 2013. Nonlinear unknown input observer design for diesel engines. American Control Conference (ACC), June 17–19, Washington, DC, USA, pp. 1076–1081. [Google Scholar]
 Teh PS. 2011. Observerbased residual design for nonlinear systems with unknown inputs. Australian Control Conference, 10–11 November, Melbourne, Australia, 198–204. [Google Scholar]
 Sherwood TK, Pigford RL, Wilke CR. 1975. Mass transfer, New York, McGrawHill. [Google Scholar]
 Schofield RW, Fane AG, Fell CJD. 1987. Heat and mass transfer in membrane distillation, Journal of Membrane Science, 33, 299–313. [Google Scholar]
Cite this article as: Khoukhi B, Tadjine M & Boucherit MS: Nonlinear observer to estimate polarization phenomenon in membrane distillation. Int. J. Simul. Multisci. Des. Optim., 2015, 6, A4.
All Figures
Figure 1. Schematic diagram of membrane distillation process. 

In the text 
Figure 2. Schematic diagram resolution algorithm. 

In the text 
Figure 3. Schematic diagram of hot flow temperature in two dimensions coordinates. 

In the text 
Figure 4. Schematic diagram of cold flow temperature in two dimensions coordinates. 

In the text 
Figure 5. Pure water flux function of temperature gradient. 

In the text 
Figure 6. Schematic diagram the nonlinear unknown input observer (NUIO) in DCMD process. 

In the text 
Figure 7. Step of hot inlet temperature introduced at t = 400 s. 

In the text 
Figure 8. Pure water flux production estimation error. 

In the text 
Figure 9. Polarization ratio in the longitudinal axe of the membrane. 

In the text 
Figure 10. Polarization coefficient estimating error inside DCMD. 

In the text 
Figure 11. Hot inlet estimating error. 

In the text 
Figure 12. Schematic diagram of temperature along z axes. 

In the text 
Figure 13. Longitudinal temperature estimating error along DCMD. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.