Assessment of the Total Volume Membrane Charge Density through Mathematical Modeling for Separation of Succinic Acid Aqueous Solutions on Ceramic Nanoﬁltration Membrane

: Nanoﬁltration of aqueous solutions of succinic acid with the addition of sodium hydroxide or magnesium hydroxycarbonate has been investigated experimentally and modeled with the comprehensively described Donnan–Steric partitioning model. The experimental retentions of acid at the same pH varied between 16% and 78%, while the estimated total volume membrane charge densities were in the range of − 35.73 and + 875.69 mol / m 3 . This work presents a novel insight into the modeling of nanoﬁltration and investigates the relations between the estimated total volume membrane charge densities, ionic strength, and component concentration on the performance of ceramic membrane. In addition, this study takes into consideration other parameters such as pH regulation and viscosities of solutions.


Introduction
Diminishing resources of fossil fuels, vulnerability of their prices, global warming, and environmental pollution result in a large interest in renewable and unconventional energy sources. Hence, one of these resources is biodiesel, which is currently under rigorous investigation [1]. Biodiesel is composed from renewable biological sources, such as vegetable oils and animal fats. It is biodegradable and nontoxic, and is characterized by low emission profiles and therefore is environmentally benign [2,3], unfortunately significant amounts of waste are generated during its production, which mainly consists of the glycerol phase. According to data presented by BP [4], the global biofuel production increased on the average by 5.13% (+4.1 Mtoe) in 2017 compared to 2014. In 2017, in the Europe region, biodiesel production reached approximately 14.167 million tons of oil equivalent (toe) whereas in North America nearly 38.190 Mtoe. Additionally, global primary energy consumption was projected to grow by 1.5 % per year between 2012 and 2035, whereas the energy from renewable resources, including biofuels, was expected to grow by 6.4% per year [5]. The biodiesel industry produces approximately 600 million tons of crude glycerol [6]. Therefore, due to the increase of glycerol waste, the technologies which allow for its beneficial reuse are of great interest. Many researchers focused on the biotechnological processing of glycerol by microbial bioconversion, delivering additional benefit to the environment in view of lower pollution. Through appropriate determine the total volume membrane charge density through parameter estimation based on the retention experiments.
Modeling of a nanofiltration membrane performance comprises two aspects: flux and rejection predictions [27]. Generally, the NF transport description models should consider the interaction between the charged membrane and ionic solutes [38]. Many charged membrane transport theories have been proposed which account for electrostatic effects as well as diffusive and convective flow to describe the solute separation [39], such as the steric-hindrance pore model (SHP), electric-steric-hindrance pore model (ESHP) [40], Teorell-Mayer-Sievers model [41], frictional model or space-charge model [42]. However, the most popular and widely adopted mathematical model for NF process is based on the extended Nernst-Planck (eN-P) equation [39]. Bowen and co-authors [25,[43][44][45] proposed the Donnan-steric-partitioning model (DSPM) which arose from the eN-P equation. The DSPM has been also used by many authors [32,[46][47][48][49] with fairy good results. However, prediction of separation or membrane charge is very limited, and modeling is used to describe and correlate rejection results.
In order to predict the separation performance, it is important to evaluate the membrane charge density in well-defined solutions. Therefore, the authors aim to consider each ion and water in the modeling, and therefore to obtain values of the total volume membrane charge densities through mathematical modeling. As a case study, the nanofiltration of aqueous solutions of succinate sodium and magnesium salts in pH equal to 9 and in two different concentrations (3.6 and 36.0 g/L) were studied experimentally and modeled with a comprehensively described DSPM model. Aside from sodium hydroxide, magnesium hydroxycarbonate (4MgCO 3 × Mg(OH) 2 × 5H 2 O) was also selected for pH adjustment of the feeds. It was because the magnesium hydroxycarbonate was used in the fermentation process reported in [50] to keep pH at approximately 8-9 in the bioreactor and at around 8.5 in the post-fermentation broth.

Experimental and Operating Conditions
The pilot plant presented in Figure 1 was used in the study of the separation of succinic acid by the NF technique. The exploited system was equipped with two membrane modules, whereas all research was carried out on the only one of them. The used membrane module was equipped with a tubular ceramic mono-channel membrane (TiO 2 active layer, TAMI Industries, Nyons, France) with a cut-off of 450 Da, support mean porosity 38.5 mm 2 , membrane mean porosity 3 µm, open porosity 30-40%, external diameter 10 mm, channel diameter and filtration area per 0.6 m tube equal to 0.0125 m 2 [51]. The transmembrane pressure (TMP) was set to 1.5 MPa. The process temperature for all experiments was fixed to 300 ± 2.0 K. The system was operated in a continuous mode, which means that the permeate and retentate were recirculated to the feeding tank continuously. The total volume of each investigated solution was equal to 12 L.
After each experiment, the ceramic membrane was cleaned to recover its initial permeance according to the procedure described in details elsewhere [51]. All cleaning steps were carried out at a flow system, the deionized water was filtered for 60 min, and the resulting flux was compared to the initial water flux. The decline in water flux after cleaning was not observed.

Analysis
The concentration of sodium ions was measured based on microwave induced plasma with optical emission spectrometry (MIP-OES) using a PLASMAQUANT 100 (Carl Zeiss, Jena, Germany), which included a microwave energy generator and resonator delivered by Plazmatronika, Wroclaw. The resonant cavity was combined with a microwave energy generator operating with frequency equal to 2.45 GHz. The generator has two cooling systems: aqueous (magnetron cooling) and pneumatic (external cavity wall cooling). More details regarding the MIP-OES used in that analysis was published elsewhere [52]. The analysis was performed using solutions with a concentration equal to 0.5000 g/L NaOH (microgranules, analytically pure, POCH, Gliwice, Poland), 0.500 g/L 4MgCO3 × Mg(OH)2 × 5H2O, 3.6, and 36.0 g/L C4H6O4 (purum p.a., >99.0%, Fluka Analytical, Sigma-Aldrich), which were appropriately diluted in order to prepare the analytical curve in the range of 0.5-30.0 mg/L. Standard solutions and samples were introduced to the excitation source (argon-helium plasma) with a rate of 1.5·10 −3 L/min delivered by peristaltic pump (SPETEC, Erding, Germany) and Figure 1. Pilot installation of the NF process: 1-feed tank; 2-rotameter; 3-heat exchanger; 4-pump; 5-membrane module; 6-temperature and flow rate controller; 7-permeate measuring cylinder; 8-pressure gauge; 9-heater.

Materials
Three one-main-component model solutions were tested in the reported study. All of them were prepared by dissolving succinic acid, either 3.6 or 36.0 g/L, in deionised water with an electrolytic conductivity not exceeding 3 µS/cm. The pH of each of the solutions was adjusted by addition of sodium hydroxide in model solutions 1 (MS1) and 2 (MS2), and addition of magnesium hydroxycarbonate in model solution 3 (MS3), with precise control of added amounts of each substance. The pH of all model solutions was equal to 9 ± 0.7. The compositions and concentrations of all model solutions were presented in Table 1. Analytical grade succinic acid was purchased from Sigma-Aldrich (pure p.a.) while sodium hydroxide and magnesium hydroxycarbonate from Avantor Performance Materials Poland S.A. (pure p.a).

Analysis
The concentration of sodium ions was measured based on microwave induced plasma with optical emission spectrometry (MIP-OES) using a PLASMAQUANT 100 (Carl Zeiss, Jena, Germany), which included a microwave energy generator and resonator delivered by Plazmatronika, Wroclaw. The resonant cavity was combined with a microwave energy generator operating with frequency equal Processes 2019, 7, 559 5 of 19 to 2.45 GHz. The generator has two cooling systems: aqueous (magnetron cooling) and pneumatic (external cavity wall cooling). More details regarding the MIP-OES used in that analysis was published elsewhere [52]. The analysis was performed using solutions with a concentration equal to 0.5000 g/L NaOH (microgranules, analytically pure, POCH, Gliwice, Poland), 0.500 g/L 4MgCO 3 × Mg(OH) 2 × 5H 2 O, 3.6, and 36.0 g/L C 4 H 6 O 4 (purum p.a., >99.0%, Fluka Analytical, Sigma-Aldrich), which were appropriately diluted in order to prepare the analytical curve in the range of 0.5-30.0 mg/L. Standard solutions and samples were introduced to the excitation source (argon-helium plasma) with a rate of 1.5 × 10 −3 L/min delivered by peristaltic pump (SPETEC, Erding, Germany) and concentric atomizer (Meinhard Glass Products, Golden, CO, USA) equipped with cyclone cloud chamber (EPOND, Vevey, Switzerland).
The contents of succinic acid and succinates in feed solutions, retentate and permeate obtained during the NF processes were determined by high-performance liquid chromatography method using HP Agilent 1100 Series system (Germany). The apparatus was equipped with an autosampler, interface (HP 35900), RI Detector (HP 1047A), pump (HP 1050), and column Rezex ROA-OrganicAcid H+ (8%), Phenomenex ® . The carrier phase was a 2.5 mM H 2 SO 4 solution with the flow of 0.5 × 10 −3 L/min. The column temperature at the input to the detector was equal to 50 • C and pressure of 5.6 MPa. Before measurement, all samples were acidified to pH ≤ 2 by addition of 25% H 2 SO 4 in an amount of 100 × 10 −6 L to 1000 × 10 −6 L of investigated solution. The pH of solutions was measured by laboratory pH-meter CP-505 (Elmetron, Poland) equipped with EPS-1 electrode.
The obtained concentrations of components were used to calculate retentions (R) according to Equation (1): All rheological studies were conducted using a rotational stress rheometer Physica MCR 501 (Anton Paar) with a double gap concentric cylinder measuring geometry (DG26.7, cup diameter 27.59/23.83 mm, bob diameter 26.66/24.66 mm). The experiments were conducted in the scanning shear rate measurements with increasing shear rates in order to obtain flow curves of succinic acid solutions in the range from 10 to 500 s −1 .

Theory of the Utilized Model
In order to describe the ion transport through the NF ceramic membrane, the comprehensively described DSPM (ddDSPM) was proposed. The ddDSPM is based on the eN-P equation, extensively utilized by Bowen and co-authors [25,30,33,44]. Description of solute fluxes in the DSPM considers convection, diffusion, and electromigration mechanism. Convection occurs due to the applied pressure difference over the membrane, diffusion due to the concentration gradient across the membrane and finally charge effects related to electrostatic repulsion between the charged membrane and a charged organic compound [53]. The proposed ddDSPM explicitly takes all ions, solutes, and solvent into account. The whole ddDSPM model used in this work consists of Equations (2)- (19) set. Equation (2) describes the solvent velocity (V), which depends on membrane properties such as pore size (r p ), porosity (A k ) and thickness of active layer (∆x), properties of separated solutions-i.e., osmotic pressures (∆π) and feed viscosity (η)-and process parameter such as transmembrane pressure (∆P). Difference of osmotic pressures defined by Equation (3) is calculated according to Equations (6) and (7) which are based on feed (x f,i ) and permeate (x p,i ) molar fractions (Equations (4) and (5)). Equation (8) defines ratio of solute (r s,i ) to pore radius, which is used along with Equation (9) to compute diffusive (K d,i ) and convective (K c,i ) hindrance factors for each component i present in the mixture. The main equations of the ddDSPM model describe the gradient of individual ion concentration (c m,i ) expressed by Equation (12) and electric potential gradient (ψ) across the membrane active layer thickness presented by Equation (13). The equations described above are solved under the condition that the membrane poses effective membrane charge density X d , which is present in Equation (14), and that the separated mixture is electroneutral (Equation (15)). The presented model is constructed under the assumption of the Donnan exclusion mechanism on the feed-membrane interface which is expressed in the form of Equation (16). The components retentions are expressed by Equation (17). In order to solve the model equations set discussed above, there is a need to define a set of boundary conditions representing component concentrations at the membrane feed (c m(0 + ),i ) and permeate (c p,i ) sides; the boundary conditions are presented by Equations (18) and (19). In summary, the ddDSPM consists of 17 + 14NC variables (listed in Table 2) in 7 + 11NC equations which were provided with appropriate descriptions in Table 3.
and boundary conditions

Differential variables Unit Number
Concentration of ion in the membrane c m,i mol/m 3 NC

Algebraic and implicit variables Unit Number
Potential gradient inside the membrane pore ψ V 1 Ratio of solute to pore radius

Total number of variables: 17 + 14NC
Although, the eN-P equations (Equations (12) and (13)) have been commonly used for the calculation of ion rejection by RO and NF membranes, they have rarely been applied to organic solutes [31]. Moreover, the ddDSPM model presented above in this work was not only used for process simulation but also for the estimation of parameter X d which is present in Equation (14), i.e., the total volume membrane charge density across the membrane active layer. Thus, the estimated X d is evaluated under the assumption of the constant surface charge and constant surface potential at the interface of pore entrance. In detail, the constant surface charge means that the pore wall surface charge density is identical to the free surface charge density, e.g., measured on particles. Whereas constant surface potential means that the pore wall surface potential at the pore entrance is identical to the free surface potential, which is related by some researchers [54] to the Donnan potential at the feed-membrane interface and therefore it is equal to the zeta potential ζ. Additionally, since X d is present in electroneutrality condition (Equation (14)), it actually combines all electrochemical interactions in close membrane neighborhood, i.e., those between solutes, solvents, and membrane material. Therefore, naming X d as the total volume charge density across the membrane active layer is justified.

Determination of Total Volume Membrane Charge Density Values in Nanofiltration
The degree of freedom (DOF) of the presented model is equal to 10 + 3NC, where NC stands for number of solutes present in the mixture. In order to solve the derived ddDSPM, the DOF must be equal to zero, therefore values of all parameters and known variables need to be provided. As already mentioned, during the parameter estimation in the ddDSPM, each ion existing in the solution is considered, even those originating from the sodium hydroxide or magnesium hydroxycarbonate, which were used for regulation of the pH of separated solutions. Values of diffusion coefficient D i , ions charge z i and radius of ions r i,s used in all calculations were presented in Table 4. Due to lack of the data, the Stokes-Einstein Equation (20) was used to determine the ionic radius of the succinate anion.
In Figure 2, the values of viscosities of investigated solutions were presented in comparison to pure water viscosities. It is important to notice that the difference between viscosities of model solutions and water varied between 4% and 24%.
The modeling in this study is considering each ion presented in the system, even ions originating from solutions used to set the desired values of pH. Such detailed approach is innovative in modeling of NF processes. Until now, researchers dealing with modeling with the DSPM model, did not consider ions originating from solutions used for regulating pH or at least had not shown it explicitly. The solutes dissociate in aqueous solutions, then they deliver specific ionic forms to the separated feed. Authors are convinced that the presence of additional ions (such as Na + , OH − , Mg 2+ , or CO 3 2− ) may influence the total volume membrane charge density. It was also assumed in the ddDSPM, that the concentrations of the components in the feed are constant (i.e., steady state model), transmembrane pressure for the entire duration of the process is constant, pores are straight cylindrical in shape and of length equal to the effective membrane layer thickness. Due to the crossflow velocity and achieved Reynolds numbers within the membrane module in experiments equal to 2.3 m/s and 19,293, respectively, it was also assumed that concentration polarization effect and fouling phenomena are negligible. Additionally, as it is well known that the NF ceramic TiO 2 membrane has a support layer (Al 2 O 3 ), the influence of that layer was neglected in the view of ratio of ions radii to support layer pore radii and the assumption that support layer is uncharged (neutral). cylindrical in shape and of length equal to the effective membrane layer thickness. Due to the crossflow velocity and achieved Reynolds numbers within the membrane module in experiments equal to 2.3 m/s and 19,293, respectively, it was also assumed that concentration polarization effect and fouling phenomena are negligible. Additionally, as it is well known that the NF ceramic TiO2 membrane has a support layer (Al2O3), the influence of that layer was neglected in the view of ratio of ions radii to support layer pore radii and the assumption that support layer is uncharged (neutral).  The parameter estimations were conducted in the gPROMS software, which employs a rigorous optimization-based approach for model validation by offering parameter estimation capabilities, i.e., fitting model parameters to experimental data. Parameter estimation in gPROMS is based on the maximum likelihood formulation which provides simultaneous estimation of parameters in the physical model of the process [60]. Assuming independent, normally distributed measurement errors, with zero means and standard deviations, that maximum likelihood goal can be achieved through the objective function presented by Equation (21) [60]. In cases discussed in this study, the parameters estimation problems had the following values of parameters following Equation (21)  The parameter estimations were conducted in the gPROMS software, which employs a rigorous optimization-based approach for model validation by offering parameter estimation capabilities, i.e., fitting model parameters to experimental data. Parameter estimation in gPROMS is based on the maximum likelihood formulation which provides simultaneous estimation of parameters in the physical model of the process [60]. Assuming independent, normally distributed measurement errors, with zero means and standard deviations, that maximum likelihood goal can be achieved through the objective function presented by Equation (21) [60]. In cases discussed in this study, the parameters estimation problems had the following values of parameters following Equation (21)

Results and Discussion
Generally, it is well known that the following key parameters primarily affect the rejection of organic solutes during NF separations and are related to: • Solute parameters: molecular size, acid dissociation constant, hydrophobicity/hydrophilicity, and diffusion coefficient, • Membrane properties: cut-off, pore size, surface charge, • Feed composition: pH, ionic strength, hardness, and the presence of organic matter.
It has to be clearly stated that the negatively hydrophilic solutes can be rejected by electrostatic repulsion through negatively charged membrane surfaces. Electrostatic interactions between charged solutes and a porous membrane have been frequently reported to be an important rejection mechanism. Ions such as Na + , K + , Ca 2+ , and Mg 2+ in feed water reduced the negative zeta potential of a membrane.
The experimentally obtained retention rates in the NF process for all model solutions were presented in Figure 3. The obtained retentions in process time slightly changed only initially and reached constant values after 40 min for whole processes which lasted up to 240 min. It has to be noticed, that two model solutions containing 3.6 g/L of succinic acid have displayed totally different values of retention. When sodium hydroxide was used as a pH regulator, the retention rate varied between 67% and 77% (see MS1 in Figure 3). When magnesium hydroxycarbonate was used, the retention decreased from significantly reaching a range of 20-22% (see MS3 in Figure 3). On the other hand, when the concentration of succinic acid increased from 3.6 g/L (MS1) to 36 g/L (MS2), the retention decreased form 67-77% to 3-16%. When retention profiles of MS1 and MS2 were analyzed separately, the initial increase of retention rate could be related to the fouling caused by concentration polarization, internal pore blocking or cake formation as it was described in [61]. However, since the MS2 has 10 times higher concentration than MS1, then it should be expected that retention rate would be kept at the same level, or that it should increase if there fouling would occur; but this was not the case. Therefore, it can be assumed that there is no typical fouling but a type of electrokinetic saturation of surface charges, the mathematical description of which is unknown to authors and requires further detailed studies.
The concentrations of sodium in the feed and permeate have changed significantly in the end of MS2 separation experiments, which is not the case in MS1 (Figure 4). That observation indicates that there is a mechanism of ion binding by the pore wall at higher concentration of sodium. Based on the obtained results, the ion binding for lower concentrations of succinic acid and pH-regulators (MS1 and MS3) can be excluded.

Impact of Dynamic Viscosity on Modeled Permeate Flux
The obtained dynamic viscosity values for all model solutions were presented in relation to the temperature in Figure 2. In modeling, it is generally assumed that viscosities of aqueous solutions with low solute concentrations are equal to water viscosity. For the MS1 at 303 K, the obtained viscosity value differed by approximately 10% in comparison with reference water viscosity. Such difference can influence the modeling results, which is highlighted by comparison of volumetric fluxes calculated with experimentally obtained viscosities and with viscosity being equal to water at certain temperatures which was presented in Figure 5. In general, the experimental viscosities in calculations of volumetric fluxes for investigated solutions resulted in a decrease of calculated volumetric fluxes from 7.5% for MS1 and up to 66.4% for MS3 to those calculated with pure water viscosities. It has to be stated that including the experimentally obtained viscosities resulted in better description of the experimental volume flux for MS2. However, for MS1 and MS3 cases the discrepancies between utilization of water and experimental viscosities were rather far from the experimental, but still of the same order. Generally, for the investigated solutions, the use of experimental viscosities with combination of the Hagen-Poiseuille equation resulted in an underestimation of the calculated volume flux.   Permeate @ t=240 Feed @ t=15 Feed @ t=240   Permeate @ t=240 Feed @ t=15 Feed @ t=240

Comparison of the ddDSPM Model with the Standard Approach
The standard approach considers only concentrations of solutes and ions coming from dissolved components, in the presented cases from succinic acid. The ddDSPM model takes into account all solutes, ions, and solvent, which include succinic acid, pH regulating solutions and water. In the standard approach, the model consists of 29  Based on the obtained results, it is difficult to postulate which modeling approach is the best, since it is impossible to experimentally obtain the overall membrane charge density Xd. However, it is evident that the obtained values of Xd for higher concentration of solutes (MS2) differ significantly depending on the applied modeling approach. At first glance, it seems that there is no clear reason to use the ddDSPM model for diluted solution such as MS1; however, when feed contains a higher number of components or polyanions like in MS3, the ddDSPM would be recommended. The presented results clearly show that use of experimentally obtained volumetric flux in comparison to the ddDSPM with experimentally obtained viscosities does not significantly influence the Xd when computed volumetric flux is in good comparison to the experimental one (MS2). Otherwise, use of experimental volumetric flux in estimation of Xd is recommended.

Variation of the Overall Volume Charge Densities in Relation to Used pH Regulator
Comparison of the estimated values of the total volume membrane charge densities with standard and detailed models was presented in Figure 6 along with ionic strength and obtained

Comparison of the ddDSPM Model with the Standard Approach
The standard approach considers only concentrations of solutes and ions coming from dissolved components, in the presented cases from succinic acid. The ddDSPM model takes into account all solutes, ions, and solvent, which include succinic acid, pH regulating solutions and water. In the standard approach, the model consists of 29 equations and 46 variables (NC = 2), while in the ddDSPM there are 62 equations with 88 variables (NC = 5), which were solved with the aim to estimate X d . Additionally, estimation of X d was also conducted under direct consideration of values of experimental fluxes (V exp ) instead of volume fluxes calculated with the Hagen-Poiseuille equation. All estimation results were compared in Figure 6. In the standard approach (i.e., DSPM model, ions coming from pH regulator not included), X d is changing between −35.59 and +278.09 mol/m 3 , while in the detailed approach (i.e., ddDSPM model, ions coming from pH regulator included) the changes ranged between −35.73 and +875.69 mol/m 3 . The use of the water viscosity in the ddDSPM resulted in X d values ranging between −34.98 and +939.67 mol/m 3 and use of experimental values of volume fluxes (V exp ) in ddDSPM (i.e., ddDSPM model, ions coming from pH regulator included) resulted in X d values ranging between −19.57 and +871.74 mol/m 3 . It is important to highlight that the retention obtained with the models overlaid with the experimental values, regardless of which model was used.
Based on the obtained results, it is difficult to postulate which modeling approach is the best, since it is impossible to experimentally obtain the overall membrane charge density X d . However, it is evident that the obtained values of X d for higher concentration of solutes (MS2) differ significantly depending on the applied modeling approach. At first glance, it seems that there is no clear reason to use the ddDSPM model for diluted solution such as MS1; however, when feed contains a higher number of components or polyanions like in MS3, the ddDSPM would be recommended. The presented results clearly show that use of experimentally obtained volumetric flux in comparison to the ddDSPM with experimentally obtained viscosities does not significantly influence the X d when computed volumetric flux is in good comparison to the experimental one (MS2). Otherwise, use of experimental volumetric flux in estimation of X d is recommended.

Variation of the Overall Volume Charge Densities in Relation to Used pH Regulator
Comparison of the estimated values of the total volume membrane charge densities with standard and detailed models was presented in Figure 6 along with ionic strength and obtained experimental rejections for each investigated solution. The ionic strength I for each model solution was calculated according to Equation (22).  As presented in Figure 6, the estimated values of Xd for each model solutions were clearly different, but these values were adequate to retention rates obtained in the NF processes. In general, the higher ionic strength of the separated solutions resulted in a higher value of the total volumetric charge of the membrane. It should be noted that Xd values for all model solutions were values obtained for the experiments which reached the steady state and these values should not be confused with the fixed membrane charge because the Xd relates all electrokinetic phenomena present in the vicinity of membrane surface and not only those related to the membrane surface as it is in the case of fixed membrane charge [62][63][64]. The achieved Xd value for MS1 (3.6 g/L succinic acid of pH = 9.7 regulated with granulate NaOH) equal to −35.73 mol/m 3 may indicate the presence of strong electrostatic repulsion between succinate anions present in aqueous solution and membrane surface functional groups, hence the Xd reached a negative value and retention achieved 78%. For two other cases, the Xd values were positive. Such behavior could be related to the amphoteric characteristics of TiO2 active membrane layer, although it is generally explained by the charge change at isoelectric point at specific pH [46]. 3 In the studied cases, when pH was equal to approximately 9, it was expected that the active layer would possess negatively charged groups. The succinate anion has the lowest diffusion coefficient among other ions present in the system (Table 4), therefore succinate anions should grind and penetrate the membrane as the last anions. Although such strongly positive Xd value obtained for As presented in Figure 6, the estimated values of X d for each model solutions were clearly different, but these values were adequate to retention rates obtained in the NF processes. In general, the higher ionic strength of the separated solutions resulted in a higher value of the total volumetric charge of the membrane. It should be noted that X d values for all model solutions were values obtained for the experiments which reached the steady state and these values should not be confused with the fixed membrane charge because the X d relates all electrokinetic phenomena present in the vicinity of membrane surface and not only those related to the membrane surface as it is in the case of fixed membrane charge [62][63][64]. The achieved X d value for MS1 (3.6 g/L succinic acid of pH = 9.7 regulated with granulate NaOH) equal to −35.73 mol/m 3 may indicate the presence of strong electrostatic repulsion between succinate anions present in aqueous solution and membrane surface functional groups, hence the X d reached a negative value and retention achieved 78%. For two other cases, the X d values were positive. Such behavior could be related to the amphoteric characteristics of TiO 2 active membrane layer, although it is generally explained by the charge change at isoelectric point at specific pH [46].
In the studied cases, when pH was equal to approximately 9, it was expected that the active layer would possess negatively charged groups. The succinate anion has the lowest diffusion coefficient among other ions present in the system (Table 4), therefore succinate anions should grind and penetrate the membrane as the last anions. Although such strongly positive X d value obtained for MS2 (36.0 g/dm 3 succinic acid at pH = 8.8 regulated with granulate NaOH) may reflect high concentration of sodium cations, their selective adsorption on membrane surface, which results in formation of additional surface layer and appearance of the electrostatic attraction of the succinic anions, resulted in very low retention (16%). In cases of basic salts, it is very difficult to determine which hypothetical mechanisms might influence the separation and could result in such charge values. Although the estimated value for the MS3 (3.6 g/dm 3 of succinic acid with pH = 8.7 regulated with 4MgCO 3 × Mg(OH) 2 × 5H 2 O) was equal to +163.27 mol/m 3 , which suggests the selective adsorption of magnesium cations or other cations on the membrane active layer surface, also within the pores.
In general, the increase of the membrane charge densities with pH increase might be caused by the selective adsorption of ions at the membrane-separated mixture interface and additional adsorption in membranes' pores in active layer [65]. In this work, authors consider the total volume membrane charge density as the sum of the fixed membrane charge density and the number of adsorbed ions in the whole active membrane volume and its close vicinity at the feed side and membrane pores. The possible mechanism for the formation of membrane charge assumes that ions are partitioned from the bulk solution into the membrane pore under the influence of the Donnan potential. Among the partitioned ions in the membrane pores, either cations or anions are adsorbed selectively by the pore walls. Next, the adsorbed ions are bound on the pore wall and provide the electric charge to the membrane, which is specific to the investigated cases. In view of all that was stated above, the values of total volume membrane charge densities X d will always be different depending on the type of solution subjected to the NF process, although the pH values are the same. Despite that pH for all studied solutions was equal to approximately 9, the X d values varied between −35.73 and +875.69 mol/m 3 ( Figure 6). Therefore, it can be assumed that, the mechanism of selective ion adsorption plays a significant role despite similar pH values, the schematic explanation of which was presented in Figure 7. The overall volume membrane charge density in MS1 is negative because not all negative membrane surface group are associated with cations present in the solution, therefore succinic ions are repulsed by the membrane interface. The MS3 contains a similar amount of added magnesium hydroxide as sodium hydroxide in MS1. Since the magnesium ion possesses two positive charges, it can interact with two negative membrane sites and more negative membrane sites can be associated. Therefore, in MS3, the succinic ions can permeate more easily than in case of MS1. The lowest retention observed for MS2 could be related to the highest overall volume membrane charge due to high concentration of dissociated succinic acid and the highest amount of added sodium hydroxide, which in that case could compensate all present negative membrane charges. even reporting the amount of added pH regulator; (7) including all components in mathematical modeling can allow for better understanding of nanofiltration separation.  Table 4, Xdtotal volume membrane charge density; Csuc2 − -molar concentration of succinate anion; ISf-m-feedmembrane interface; ISf-m-permeate-membrane interface).
Based on the obtained results and authors knowledge, it is impossible to state which modeling approach is the best, however it is clear to recommend the comprehensively described Donnan-Steric partitioning model to analyze the separation of multicomponent mixtures as it gives reasonable results and takes into account all components present in the mixture.

Conflicts of Interest:
The authors declare no conflict of interest.   Table 4, X d -total volume membrane charge density; C suc2 − -molar concentration of succinate anion;

Conclusions
The main goal of the performed computer-aided simulations was to estimate the membrane surface charge densities with the use of the comprehensively described Donnan-Steric partitioning model, which is derived from the extended Nernst-Planck equation with Donnan partitioning assumption. The obtained total volume membrane charge densities X d are consistent with experimental values of retention. Values such as charge density of the membrane are very important for explanation of the mechanism of ions transport across the membrane, defining retention, and describing influences on electrostatic repulsion between ions and membrane.
The presented work can be summarized in the following points: (1) separation of organic acids, such as succinic acid, on the ceramic membrane is affected by the acid concentration and the type of used pH regulator; (2) experiments clearly show that increase of the concentration of components can decrease the retention of separated components (MS1 and MS2); (3) use of a different pH regulator can dramatically change the performance of the separation (MS1 and MS3); (4) in case of low succinic acid concentration, the membrane functional group (-TiO 2 ) dissociation has the greatest impact on the charge formation, whereas in case of higher succinic acid concentrations, it is postulated that the charge formation is determined by adsorption of specific ions on the membrane active layer; (5) defining separation performance at certain pH can be totally different in the vicinity of other components or impurities; (6) detailed description of experiments is required in order to compute X d , even reporting the amount of added pH regulator; (7) including all components in mathematical modeling can allow for better understanding of nanofiltration separation.
Based on the obtained results and authors knowledge, it is impossible to state which modeling approach is the best, however it is clear to recommend the comprehensively described Donnan-Steric partitioning model to analyze the separation of multicomponent mixtures as it gives reasonable results and takes into account all components present in the mixture.

Conflicts of Interest:
The authors declare no conflict of interest.  Boltzmann constant (1.38×10 −23 J/K) r p pore radii, m r s,i ion radii, m x f,i molar fraction on the feed side, mol/mol x p,i molar fraction on the permeate side, mol/mol z i charge of individual ion Greek Letters θ set of model parameters to be estimated Ψ potential gradient inside membrane pore, V Ψ D Donnan potential, V ∆π osmotic pressure difference, Pa ∆P transmembrane pressure, Pa ∆x thickness of membrane active layer, m η viscosity, Pa × s ηs solvent viscosity, Pa × s λi ratio of solute to pore radius π permeate osmotic pressure on the permeate side, Pa π feed osmotic pressure on the feed side, Pa σ 2 ijk variance of the k-th measurement of variable j in experiment i Φ i steric term