Dynamical Modeling of Water Flux in Forward Osmosis with Multistage Operation and Sensitivity Analysis of Model Parameters

: To mathematically predict the behavior of a forward osmosis (FO) process for water recovery, a model was constructed using an asymmetric membrane and glucose as a draw solution, allowing an examination of both phenomenological and process aspects. It was found that the proposed model adequately described the signiﬁcant physicochemical phenomena that occur in the FO system, including forward water ﬂux, internal concentration polarization (ICP), external concentration polarization (ECP), and reverse solute di ﬀ usion (RSD). Model parameters, namely the physiochemical properties of the FO membrane and glucose solutions, were estimated on the basis of experimental and existing data. Through batch FO operations with the estimated parameters, the model was veriﬁed. In addition, the inﬂuences of ECP and ICP on the water ﬂux of the FO system were investigated at di ﬀ erent solute concentrations. Water ﬂux simulation results, which exhibited good agreement with the experimental data, conﬁrmed that ICP, ECP, and RSD had a real impact on water ﬂux and thus must be taken into account in the FO process. With the Latin-hypercube—one-factor-at-a-time (LH–OAT) method, the sensitivity index of di ﬀ usivity was at its highest, with a value of more than 40%, which means that di ﬀ usivity is the most inﬂuential parameter for water ﬂux of the FO system, in particular when dealing with a high-salinity solution. Based on the developed model and sensitivity analysis, the simulation results provide insight into how mass transport a ﬀ ects the performance of an FO system.


Introduction
Water shortages have become a key issue facing humanity. According to the United Nations World Water Development Report in 2019, over 2 billion people suffer from severe water shortage, and the global demand for fresh water has been increasing by about 1% annually since the 1980s [1]. Consequently, much effort has been made to secure water, in particular safe water, with low energy consumption. One such means is membrane-based water treatment and desalination technology, which is relatively energy-efficient and independent of the water cycle [2]. glucose as draw solute. To validate the proposed model, FO filtration experiments were carried out with a commercialized FO membrane in batch mode using glucose as a draw solute and deionized water and glucose solution as a feed solution. For the batch operation of the FO process, the model prediction of the permeate flux profiles was found to be in line with the experimental data. To obtain insight into the dominant factors affecting permeate water flux, global sensitivity analysis using the Latin-hypercube-one-factor-at-a-time (LH-OAT) method was carried out.

Water Flux in Forward Osmosis Process
The transport of water molecules through FO membranes can be represented by irreversible thermodynamics, with the membrane treated as a black box, as has been done for biological membranes [24]. In this approach, the water flux (J w ) for the FO process (Figure 1), in which hydraulic pressure differences are absent, can generally be described as follows, where A is the water permeability, π i represents the osmotic pressure at the interface between selective layer and support layer, and π Fm represents the osmotic pressure at the interface between selective layer and feed solution as shown in Figure 1. This equation indicates that the water flux is linearly proportional to the effective osmotic pressure difference, π i − π Fm , across the FO membrane selective layer. With the presence of a concentration gradient, draw solute is diffused into feed solution across the boundary layers (δ D and δ F ), across the support layer with thickness t s , and across the selective layer. C Dm , C i , and C Fm represent the solute concentrations at the interfaces between (1) the draw solution and the support layer, (2) the support layer and the selective layer, and (3) the selective layer and the feed solution, respectively. π Dm , π i , and π Fm represent the osmotic pressure of solute concentrations at each interface.

Osmotic Pressure of Glucose
To establish the relationship between osmotic pressure and solute concentration, the van't Hoff (1887) equation can be applied for the osmotic pressure (π) for an ideal dilute solution [25], that is, where n represents the van't Hoff factor (n = 1 for glucose), C is the solution molar concentration (molarity), R = 0.0821 L·atm·mol −1 ·K −1 the ideal gas constant, and T is absolute solution temperature (K). Although the van't Hoff equation is well fitted for dilute and ideal solutions in which ions do not affect each other, this is not applicable to the FO process, which deals with highly concentrated draw and feed solutions [25]. In order to more precisely model the osmotic pressure of general solutions, a modification is made in the van't Hoff equation, which is then expressed by a virial expansion to a power series [26,27], as follows: In this equation, the value c represents the solution concentration (g·L −1 ), M w is the solute molecular weight, and B 1 , B 2 , and B 3 are virial osmotic coefficients. It is possible to obtain these virial coefficients empirically by using Equation (3) to fit the experimental osmotic pressure data. It is generally accepted that it is sufficient to determine the coefficients up to the second coefficient (B 1 and B 2 ) for the purpose of reproducing the observed osmotic pressure [28].

Concentration Polarization
Many experimental and theoretical studies have demonstrated that observed water flux values are significantly smaller than those calculated on the basis of difference in bulk concentration. The reason behind this discrepancy is the formation of ICP and ECP; both decrease the effective concentration slope across the FO membrane selective layer; as a result, the value of water flux that is observed, J w , is inevitably lower than expected [29].

Internal Concentration Polarization
It is known that ICP in the support layer can cause severe degradation of cross-membrane water permeation in the FO process [17]. For a porous support layer, in which draw solution is in contact with the porous support, there is a concentration gradient that is steeper than that found in bulk solution. This happens because of decreased solute diffusivity due to porosity and tortuosity. This dilutive ICP brings about a decrease in the net concentration difference across the selective layer and results in a corresponding decline in water flux [29]. According to previous studies [14,23], the solute concentration at the interface between support and selective layers (C i ) can be described as follows: Here, C Dm represents the concentration of the solute on the surface of the support layer membrane and t s indicates thickness of the support layer D s is the reduced solute diffusivity in the porous support layer [30]. S denotes the membrane structural parameter of the support layer, described as follows [14], Equation (5) pertains to the effective distance of the support layer through which a solute must pass to move to the selective layer from the boundary area of the support layer and the draw solution. In experiments using both RO and FO, the average distance of S can be determined and has been described in the literature [14].

External Concentration Polarization
ECP is of great importance in any pressure-driven desalination process, e.g., reverse osmosis [31]. This holds true for the FO system, as the presence of ECP lowers the water movement via a reduced effective concentration gradient through the membrane. This phenomenon was reflected even in some early FO models, such as those of Phillip [14] and Suh and Lee [23]. In fact, because the FO membrane is in contact with two different solutions, ECP can arise on both sides of the membrane surface. To elaborate, the feed side membrane surface faces a concentrated feed solution, a phenomenon termed concentrative ECP; the permeate side membrane surface faces a diluted solution, a phenomenon termed dilutive ECP. According to previous studies [14,23], the solute concentrations at the membrane surface on the support layer (C Dm ) and selective layer (C Fm ) can be estimated as follows, Here, C Db and C Fb are the bulk solute concentration of draw and feed solutions, respectively. D δ can be described in relation to mass transfer coefficient (k); this value has a direct relationship to the Sherwood number (Sh), and also corresponds to the hydrodynamics of the system (d h ), which can be estimated in the previous studies [14,19,23,31].

Reverse Solute Flux (RSF)
J s can also be described as shown below [16]: Here, B is the salt permeation coefficient. Substituting Equations (4) and (7) for C i and C Fm into Equation (8) can represent the effective solute flux existing between the draw solution and the feed solution, with measurable variables, Note that Equation (9) contains the proportion of RSF to water flux (J s /J w ) as a repetitive term, referred to as the specific RSF [32]. Past research has demonstrated that the specific value of RSF (J s /J w ) can be replaced by a constant [14,27,33], meaning that more water flux leads to more draw solute moving through the membrane. Philip et al. [14] and Suh and Lee [23] also dealt with the selectivity of reverse flux (J w /J s ), which is designated as the reverse form of the specific value of RSF; they also validated the dependency of the reverse flux selectivity (J w /J s ) on the characteristics of the membrane selective layer (water (A) and solute (B) permeability values in the Appendix A), as follows: This result provides insight into the selectivity of the reverse flux, showing that it is independent of the concentration of the bulk draw solute, the crossflow velocity, and structural parameter S, and can be handled as a constant coefficient determined in experiments, provided that these experiments are conducted with the same membrane, the same draw solutes, and identical temperature. Thus, the specific RSF (J s /J w ) can be described as a constant model parameter: Finally, the equation of reverse solute flux and forward water flux for the FO system can be rearranged as follows: It should be noted that for Equation (12), C Db and C Fb are easily measurable independent variables, and A, D, S, and J sw are model parameters. Before simulating FO performance, these model parameters should be determined by experiments. Then, only one unknown dependent variable of the equation remains, namely J w . However, Equation (12) is implicit forms with regard to J w , and cannot be solved explicitly. Consequently, recursive numerical procedures must be used to solve these implicit flux equations.

Mass Balance for Dynamic Modeling (Multistage Operation)
When FO is carried out as a batch operation in such a way that there is an increase of the draw volume and a decrease of the feed volume, such that water molecules continually traverse the FO membrane from feed solution to draw solution, as shown in Figure 2, the time profile for the entire system changes analogously depending on the actual time and can, if the observation time interval (or sampling time interval) of each stage is sufficiently short but not too short (because mass transfer should be reached in quasi-steady state), be expressed as a discrete multistage process, as shown in Figure 3. Considering module length and average crossflow velocity of solution, we set the observation time interval or sampling time interval to 0.5 s because the water molecule takes 0.5 s after entering to exit the FO module. It is assumed that the FO system reaches quasi-steady state at each stage during this sampling time interval. The mass balance of solutes for feed and draw solutions at the j-th stage (or at the j-th sampling in Figure 3) of the FO unit for batch operation can be represented by and Here, C represents the concentration of individual components, V is the volume, J s is the reverse solute flux, A m is the effective area of the membrane, t s is the time interval, the subscript Fb is the bulk feed solution, the subscript Db is the bulk draw solution, the subscript D is the draw solution, the subscript F is the feed solution, and the superscript j is the sampling index.
The total volume of the system, including the feed and draw solutions, remains constant because the volume loss of the feed solution via transverse-FO-membrane water flux permeation is included in the draw solution volume. This relationship is presented below. and Here, J w is forward water flux.

Numerical Simulations
Using the proposed model for the FO process, all iterative calculations for the water flux profile were conducted with MATLAB software under conditions identical to those of the batch operation experiments described above. Figure 4 provides a flow chart of the proposed multistage FO modeling procedure. First, the timer was reset, and invariant variables were initialized. Then, the initial physio-chemical parameters of the feed and draw solutions of the j-th stage are set. In this step, the initial estimate of water fluxĵ w (t) was also properly set. With the given physio-chemical variables and j w (t), calculations of the dynamic viscosity, density, and mass transfer functions of the feed solution and the draw solution were performed, yielding feed and draw solution osmotic pressure values. From the obtained osmotic pressures, the water flux J w (t) was obtained. This calculation was repeated with the updatedĴ w (t) using the mean value ofĴ w (t) and J w (t) until the error tolerance condition, J w (t) −Ĵ w (t) ≤ err, was satisfied. The implicit equation in the proposed model was solved using the bisectional method: The embedded MATLAB 'fzero' function [23]. If this j-stage is not reached by the final stage, as we had hoped, i.e., t = T, the volume and feed and draw solution concentrations were revised using the value of J w (t) determined for the next stage.

Sensitivity Analysis Using Latin-Hypercube-One-Factor-at-a-Time (LH-OAT) Method
Sensitivity analyses of model parameters and operating conditions were executed with the Latin-hypercube-one-factor-at-a-time (LH-OAT) method [34], combing Latin-hypercube (LH) sampling [35] and one-factor-at-a-time methods (OAT) [36]. In the LH method, proposed by McKay for use instead of the Monte Carlo method, the parameter space is divided into N intervals with the same probability, and only one variable is randomly extracted from each interval and analyzed by multivariate linear regression. Though this method is advantageous compared to other global sensitivity analyses in that its computational calculation is efficient [34,37], it has limitations in that a linear regression analysis is assumed and the sensitivity to one specific individual variable cannot be identified [34]. To analyze the sensitivity in the OAT method, on the other hand, only one variable in the parameter space is sequentially selected for a small change of the selected parameter, while other parameters are fixed as constant [36,38]. In general, the OAT method has been widely used as an efficient method of local sensitivity analysis. One downside is that it does not give any information on global sensitivity for the whole parameter space. By combining LH sampling and OAT design, therefore, the merits of both methods can be exploited [39]. The process begins with sampling N LH points randomly for N intervals as initial points; this is followed by OAT analysis, in which each of the P parameters is changed [34]. For sensitivity analysis, parameters affecting FO performance were investigated by changing figures with standard deviations of 10%.

Preparation of Feed and Draw Solution
ACS grade glucose (Daejung Chemicals and Metals Co. Ltd., Busan, Korea), which is neutral and highly soluble in water, was utilized as a draw solute. For all empirical procedures, glucose was dissolved in deionized (DI) water (Merck, Daejeon, Korea), which has a resistance value of 18.2 MΩcm at concentrations in the range of 0 to 2.0 M. Viscosity and osmotic pressures for each concentration were measured using a viscosity meter, SV-10 (A&D, Seoul, Korea) and an osmometer, Osmomat 0300-D (Gonotec, Berlin, Germany), respectively. Binary bulk diffusion coefficients for glucose and DI water were maintained at a constant value (6.7 × 10 −10 m 2 /s) [40,41]. The density of each glucose solution was calculated using the Aspen HYSYS ® (Cambridge, MA, USA) chemical database. The characteristics of glucose solutions are summarized in Table 1.

Membrane Preparation and Crossflow Setup
A commercial thin-film composite (TFC) membrane (Porifera, Hayward, CA, USA), known for use in osmotically driven processes, was employed for the FO experiments. This membrane, which has a flat-sheet form, has an asymmetric structure of porous support layer and dense selective layer. This membrane has seen wide use in past research [42][43][44]. Referring to previous studies [42,43], the intrinsic performance parameters of the membrane are summarized in Table 2. The purchased membrane was preserved in DI water at a temperature of 4 • C and was rinsed with DI water prior to use. A widely used bench-scale FO set up was constructed and used to estimate the permeated water flux [11]. The crossflow membrane unit (2.6 cm width × 0.9 cm length × 0.3 cm depth) consisted of an FO cell with channels present on the two membrane sides to promote the flow of feed and draw solutions. As can be seen in Figure 2, the feed solution moved into and out of the feed chamber via a gear pump (Cole-Parmer, Daejeon, Korea). Simultaneously, the draw solution presents in the reservoir circulated into the draw chamber via the same type of pump. Directions of the flow of feed solution and draw solution in the FO module were co-current. The flow rate was maintained at 1.5 L/min by use of a flow meter OM006 (Corea Flow, Daejeon, Korea). The temperature was maintained at 25 ± 1 • C using a heat exchanger (Thermo Fisher, Seoul, Korea A) with a constant-temperature bath. The resulting permeate water flowing through the FO membrane was quantified using a computer-linked analytical balance (CAS, Daejeon, Korea). Operating conditions for simulations and experiments are summarized in Table 3.

Forward Osmosis Runs for Model Validation
After the membrane was loaded into the unit module with FO mode orientation, the FO system was first operated using deionized water for the feed and draw solutions for 0.5 h; this was done to achieve temperature equilibrium and operational stabilization. Then, to achieve the necessary concentration, a designated quantity of stock glucose solution (4 M) was put into the draw solution tank; the permeated water flux was then quantified. After measuring the water flux, a certain amount of glucose stock solution was put into the draw solution tank to establish the next desired concentration. The aforementioned procedure was continued in a consecutive manner to draw solution concentrations of 0.5, 1, 1.5, and 2 M glucose. After measuring the flux for the last and highest concentration (i.e., 2 M) with a feed solution of DI water, a certain amount of the stock 4 M glucose solution was put into the feed solution to obtain the required feed solution concentration. Flux was measured at 0.5, 1.0, and 1.5 M concentrations of glucose with the draw solution concentration fixed at 2.0 M.

Parameter Estimation
Among the model parameters, the following physiochemical properties of the glucose solution were explored: Osmotic pressure (π), viscosity (µ), and density (ρ). As shown in Figure 5a, the osmotic pressure of the glucose solution, which was measured experimentally using an osmometer, was different from that obtained by van't Hoff equation calculation in the region of the high glucose concentration. A more precise prediction can be achieved by means of a virial equation using a power series, as shown Equation (3). Using the second-order polynomial regression of π cRT and the glucose concentration, the virial coefficients B 1 and B 2 in Equation (3) can be estimated from the experimental osmotic pressure measured at different glucose concentrations, as shown in Figure 5b and Table 1. The viscosity of the solution affects the mass transfer during filtration. The decreased mass transfer rate due to the high viscosity aggravates ECP and thereby ends in a reduction of the water flux. Using a second-order polynomial regression, the viscosity (µ) of the glucose solution can be empirically shown to be a function of the molar concentration (C), as shown in Figure 5c and Table 1. Their relationship can be obtained as follows, in which the coefficients are obtained by the second-order polynomial regression µ = 8.60 × 10 −1 + 2.58 × 10 −1 C + 4.63 × 10 −1 C 2 . According to data sourced from the Aspen HYSYS ® (Cambridge, MA, USA) chemical database, the density (ρ) of the glucose solution also increased slightly in line with its molar concentration (C). This relationship can be expressed as shown below for which the coefficients are obtained by a second-order polynomial regression Among the model parameters, the physio-chemical properties of the FO membrane (Porifera, USA) have been thoroughly studied, with several reports in the literature [42][43][44] focusing on the following: The water permeability (A), the structural parameter (S), and the specific reverse solute (glucose) flux (J sw ) of the membrane. The water permeability and structural parameter were sourced from the literature [43]. The specific reverse solute (glucose) permeability of the membrane J sw was found in another work [42]. These physiochemical properties of FO are summarized in Table 2. Figure 6 shows the model verification process, which relied on comparisons of predicted data obtained using the proposed model and experimental data obtained under identical conditions for FO operations in the FO mode orientation described above. Water flux was plotted against the net concentration difference of the bulk feed and glucose draw solutions (Figure 6a). Dilutive ICP and ECP are on the side of the draw solution. The solid line and the white circle indicate the flux; the feed solution is DI water. When concentrative ECP, dilutive ECP, and ICP are present, the dashed line and black circle show the flux, with glucose solution as the feed. The model predictions of the flux, depicted as solid and dashed lines, agreed well with the empirical results, depicted as white and black circles. More specifically, the experimental flux data were plotted against the data for the estimated flux shown in Figure 6b. The slope of this relationship lies near the dashed line (slope = 1), which means that the model predictions and the experimental results are in perfect agreement. The mean squared error (MSE) between the predicted water flux and the experimental flux data is 0.09, and the R-squared value (R 2 ) is 0.99. The water flux profiles for the batch operation plotted against time are shown in Figure 6c. The blue dashed line represents the process model flux; this is based on the assumption that the glucose solute used is 100% pure. The model prediction of the water flux, represented as the blue dashed line in Figure 6c, initially agreed to a moderate degree, within ±1.0%, with the empirical data; however, over time, this value diverged from the experimental data. This unexpected discrepancy might be due to impurities that were actually contained (around 2%) in the supposedly pure chemical; worse yet, all compounds were ionic compounds, which was not taken into consideration in the modeling process. There is the possibility that the unknown compounds diffused from the draw solution to the feed solution and as a result, were able to build up, in a gradual manner, in the feed solution, impacting considerably the net driving concentration gradient. It is not uncommon that FO models pay too little attention to impurities, and rightly so, seeing that the primary goal is to obtain values for the instantaneous water flux or initial water flux of the FO system, and not the flux change over time. In this study, the 2% impurities were treated as an ionic pure chemical of NaCl, and their behavior was reflected in the modified process model. With the properties of NaCl as listed in Table 4, the modified process model can be depicted by the red solid line in Figure 6c. The modified model prediction of water flux showed better fit to the experimental data than did the unmodified model prediction. All of this demonstrates that the model precisely represents the practical mass transfer phenomenon occurring in the FO membranes.

Effects of Changes in the Concentration Difference on Concentration Polarization
The developed FO model was used to investigate how changes in the bulk concentration differences influence the concentration polarization and water flux. Using deionized water as a feed solution and glucose as a draw solution, with concentrations in the range of 0.5 to 2.0 M, Figure 7 shows the proportion of each concentration polarization and resultant water flux. All model conditions were identical to those in the previous parameter estimation and model validation. The side view of Figure 7a shows a profile of solution concentration change across the membrane when 2 M glucose is used for the draw solution; it can be interpreted that the water flux increment was reduced as the draw solution concentration increased. The percentage of effective concentration difference for water transfer across the selective layer decreased from 27.3% to 11.3% when the concentration of the draw solution was increased; on the other hand, the proportions of dilutive ICP and ECP increased. The decline in the increase of the slope of the water flux can be explained in terms of the comparatively increased ratio of ECP on the support layer surface and support layer ICP.
The most likely reason behind the driving force reduction is probably the presence of internal ICP [22,29,45]. In addition, the ratio of ICP increased from 52.2% to 56.8% with the increase of the draw solution concentration shown in Figure 7a. Surprisingly, there was a significant increment of the dilutive ECP in the draw solution side boundary layer; the increase was in the range of 20.4% to 31.9% as the solute concentration of the draw solution increased. Dilutive ECP may have a crucial role under certain circumstances, e.g., when a mixed solution is applied and/or when low crossflow velocities or high water flux levels are used [46,47]. However, the effect of ECP in a concentrative form on the side of the feed solution, with deionized water as a feed solution, was found to be insignificant for all conditions of the draw solution concentration because the simulation also clearly revealed that the concentrative ECP decreased from 0.07% to 0.03%, which is not shown in Figure 7.  As presented in Figure 7b, when the concentration of the feed solution surpasses a certain level, the ratio of the concentrative ECP increased from 0.1% to 12.9% under the condition of the ideal net concentration difference in bulk phase, which implies that the concentrative ECP should not be assumed to have little influence on the water flux; rather, it is because the absolute concentrations of the two solutions, both the feed and the draw, was increased. Furthermore, the ICP was increased from 52.2% to 68.1% when the draw solution concentration was increased from 0.5 to 2.0 M, suggesting that the concentration of the draw solution has a direct effect on the severity of the ICP. Resultantly, both processes end up slowing down the diffusional movement of water molecules. When the feed solution is highly concentrated, the draw solution concentration has to be correspondingly high in order for water permeation to take place, all because of natural osmotic pressure. However, Figure 7b reveals that the high salinity of both solutions worsens the level of ICP and concentrative ECP. This result is consistent with the results of previous studies involving NaCl as a draw solution [23].

Sensitivity Analysis
In the proposed model, parameters that can possibly affect the permeated water flux can be categorized into operating conditions, membrane properties, and solution properties, as listed in Tables 1-3. In this study, the bulk diffusion coefficient of glucose (D), the first and second virial coefficients (B 1 and B 2 ), the viscosities of the feed and draw solutions (µ F and µ D ), the densities of the feed and draw solutions (ρ F and ρ D ), the water permeability (A), the structural parameter (S), the specific reverse solute flux for glucose (J sw ), the linear velocities of the feed and draw solutions (ν F and ν D ), the hydraulic diameters of the feed and draw channels (d h F and d h D ), and the system temperature (T°C) were selected for the sensitivity analysis. This sensitivity analysis was meant to provide relative and quantified indices of influential factors and thereby insight into the characteristics of the FO process. Figure 8 shows the sensitivity indices of the 15 parameters under different conditions. Similar to previous studies [23,48], four different feed concentrations, for which the glucose concentration of the draw solution varies from 0.5 to 2.0 M while maintaining the 0.25 M glucose concentration of the feed solution, were applied while the draw concentration was kept constant at 2.0 M (Figure 8a), and four different draw concentrations were used when the feed concentration was fixed at 0.25 M (Figure 7b). Under all different combinations of the feed and draw concentrations, the diffusivity of the draw solute (D) had the highest sensitivity index, followed by S and A. An increase of S negatively affects the permeated water flux, while increases in D and A positively improve it. Similarly, because D over S ( D S ) is closely tied to the solute mass transfer in the support layer, as described in Equation (4), this value can have a dominant effect on the ICP; that is, an increase of D S alleviates the ICP. As discussed above based on the simulation results, the ICP was a dominant hindrance that caused the water flux to drop. Thus, these sensitivity analysis results, showing that D and B are the parameters that have the most influence on the permeated water flux, are indeed in line with previous simulation results and other studies [23,48]. These results prove that not only developing a support layer but also choosing an appropriate draw solute is important when it comes to FO performance.
When either the feed solution concentration was raised while keeping the draw solution fixed, or the draw solution concentration was increased while keeping the feed solution fixed, the sensitivity indices D and S rose while that of A decreased, as shown in Figure 7. As mentioned for the previous simulation, the ICP and ECP have specific effects on the water flux that increase according to the concentrations of the feed and draw solutions, implying that D and B are more related to the ICP and ECP compared to A. Thus, when dealing with high salinity solutions, the augmentation of D S is probably the best method of maximizing the water flux. ν is average crossflow velocity, d h is hydraulic diameter of channel, T°C is temperature (°C ) of system, D is diffusivity, µ is solution viscosity, ρ is solution density, S is structural parameter of membrane, A is water permeability of membrane, and J sw is specific reverse solute flux. The subscripts F and D indicate feed solution and draw solution, respectively.

Conclusions
In this paper, a numerical dynamic model was established, capable of predicting results for a process of forward osmosis (FO). After a parameter estimation process, the model was found to be capable of describing significant physio-chemical phenomena during the FO process, such as the internal (ICP) and external concentration polarization (ECP), as well as diffusion of the reverse draw solute. The proposed model was verified through comparisons with experimental data. The results of our simulation agree well with the empirical data. Furthermore, the influences of the ICP and ECP on the FO system water flux were investigated with different solute concentrations. The simulation results indicate that the influences of the ICP, ECP, and reverse draw solute flux must be taken into account for FO systems. It was also observed that the concentrative ECP on the selective layer surface does not need to be taken into account when deionized water is used as a feed solution. However, the concentrative ECP is not negligible if the salinity of the feed solution exceeds a certain point. Similarly, the high salinity of a draw solution can affect the dilutive ECP on the support layer surface, because the changes in the ECP can have an effect on the support layer dilutive ICP. Furthermore, with the verified model, a global sensitivity analysis was used to consider the effects of certain model parameters on the FO performance. The simulation results confirmed solute diffusivity is the most influential parameter for water flux; the reason is that the solute diffusivity directly affects both ECP and ICP.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.