Modeling and CFD Simulation of Water Desalination Using Nanoporous Membrane Contactors

Numerical evaluation of sweeping gas membrane distillation for desalination of water towards water sustainability and environmental protection

Sweeping gas membrane distillation (SGMD) is considered a membrane distillation configuration. It uses an air stream to collect the water vapour. A 2D mathematical model is prepared in the current study to predict the effect of various operating parameters on the SGMD performance. Also, the temperature distribution in the SGMD was obtained. The effect of air inlet temperature, salt concentration, feed and air flowrate on air and salted solution outlet temperature and vapour flux through the membrane is investigated. There was good agreement between experimental data and modelling outputs. It was found that increase in air inlet temperature from 40 to 72 °C was increased the outlet temperature of air stream and cold solution from 37 to 63 °C and 38 to 65 °C respectively. Furthermore, increase in air inlet temperature led to the enhancement of vapour flux in the membrane distillation. Also, the salt concentration and feed flow rate did not have meaningful influence on the outlet temperatures, however, the flux was increased by increasing feed flowrate.

Pore Size Dependence of Permeability in Bicontinuous Nanoporous Media

In this work, we employ many-body dissipative particle dynamics (mDPD) simulations to investigate the fluid flow process through bicontinuous nanoporous media, which are representative models for a broad class of nanoporous materials. The mDPD formulation includes attractive and repulsive interactions describing accurately fluid-fluid and fluid-solid interactions. As a mesoscale simulation method, mDPD can bridge the length and time scale gap between continuum and atomistic simulations. The bicontinuous nanoporous models are constructed considering a defined morphology, the porosity level, and varying pore sizes in the range from 3.41 to 13.63 nm. All models have a 0.65 porosity level and the same topology. The models provide a stochastic description of the morphology and pore size distribution and allow for a direct investigation of the dependence of permeability on the average pore size. The stationary nanoporous models are filled with fluid particles, and flow is induced by the action of confining pistons. Simulation results, obtained by imposing different pressure differences on the surfaces of the nanoporous media, indicate a linear pressure drop within the nanoporous model. Regardless of the complexities and different scales of the porous media considered, the steady-state fluid flow through the nanoporous models is proportional to the pressure gradient applied, in agreement with Darcy's law. The calculated pore size dependence of permeability is well described by the Hagen-Poiseuille law, considering a single shape correction factor that accounts for the flow resistance due to the complex nanoporous morphology. This work highlights the effect of the average pore size of a complex stochastic bicontinuous nanoporous medium on fluid properties. The results indicate rather a relatively simple dependence of permeability on the average pore size. The novel method we employ to generate the stochastic bicontinuous nanoporous structure allows the control of different geometric features that can be explored in future studies.

Modeling and Simulation of Fabricated Graphene Nanoplates/Polystyrene Nanofibrous Membrane for DCMD

Membrane distillation is an active technique that provides pure water with very good rejection and could be applied to water of extremely high salinity. The low productivity of membrane distillation needs intensive efforts to be competitive with other desalination techniques. In this current study, a composite (PS/GNP) membrane, which is composed of polystyrene (PS) based and 0.25% weight percent graphene nanoplates (GNP) has been fabricated via electrospinning and compared with the blank PS membrane. SEM, FTIR, contact angle and porosity characterization have been performed, and the results show that the validity of the predefined conditions, and the contact angle of the composite membrane, which is found to be 91.68°, proved the hydrophobic nature of the composite membrane. A numerical simulation using Ansys 2020 software has been introduced to study the performance of the fabricated composite membrane when used in direct contact membrane distillation (DCMD). The numerical model has been validated with experimental work from the literature and showed an excellent match. The blank PS and composite PS/GNP membranes have been investigated and compared at different operating conditions, i.e., hot water supply temperature and system flow rate. The results show that the composite PS/GNP membrane outperforms the blank PS membrane at all studied operating conditions.

Exergy Analysis of a Direct Contact Membrane Distillation (DCMD) System Based on Computational Fluid Dynamics (CFD)

Understanding the energy efficiency of direct contact membrane distillation (DCMD) is important for the widespread application and practical implementation of the process. This study analyzed the available energy, known as exergy, in a DCMD system using computational fluid dynamics (CFD). A CFD model was developed to investigate the hydrodynamic and thermal conditions in a DCMD module. After the CFD model was verified, it was used to calculate the temperature polarization coefficient (TPC) and exergy destruction magnitudes under various operating conditions. The results revealed that slight decreases and increases in the TPC occurred with distance from the inlet in the module. The TPC was found to increase as the feed temperature was reduced and the feed and permeate flow rates were increased. The exergy destruction phenomenon was more significant under higher feed temperatures and higher flux conditions. Although the most significant exergy destruction in the permeate occurred near the feed inlet, the effect became less influential closer to the feed outlet. An analysis of exergy flows revealed that the efficiency loss in the permeate side corresponded to 32.9–45.3% of total exergy destruction.

Computational modeling of drug separation from aqueous solutions using octanol organic solution in membranes

Continuous membrane separation of pharmaceuticals from an aqueous feed was studied theoretically by development of high-performance mechanistic model. The model was developed based on mass and momentum transfer to predict separation and removal of ibuprofen (IP) and its metabolite compound, i.e. 4-isobutylacetophenone (4-IBAP) from aqueous solution. The modeling study was carried out for a membrane contactor considering mass transport of solute from feed to organic solvent (octanol solution). The solute experiences different mass transfer resistances during the removal in membrane system which were all taken into account in the modeling. The model’s equations were solved using computational fluid dynamic technique, and the simulations were carried out to understand the effect of process parameters, flow pattern, and membrane properties on the removal of both solutes. The simulation results indicated that IP and 4-IBAP can be effectively removed from aqueous feed by adjusting the process parameters and flow pattern. More removal was obtained when the feed flows in the shell side of membrane system due to improving mass transfer. Also, feed flow rate was indicated to be the most affecting process parameter, and the highest solute removal was obtained at the lowest feed flow rate.

Molecular separation of ibuprofen and 4-isobutylacetophenone using octanol organic solution by porous polymeric membranes

Molecular separation of pharmaceutical contaminants from water has been recently of great interest to alleviate their detrimental impacts on environment and human well-being. As the novelty, this investigation aims to develop a mechanistic modeling approach and consequently its related CFD-based simulations to evaluate the molecular separation efficiency of ibuprofen (IP) and its metabolite 4-isobutylacetophenone (4-IBAP) from water inside a porous membrane contactor (PMC). For this purpose, octanol has been applied as an organic phase to extract IP and 4-IBAP from the aqueous solution due to high solubility of solutes in octanol. Finite element (FE) technique is used as a promising tool to simultaneously solve continuity and Navier-Stokes equations and their associated boundary conditions in tube, shell and porous membrane compartments of the PMC. The results demonstrated that the application of PMC and liquid-liquid extraction process can be significantly effective due to separating 51 and 54% of inlet IP and 4-IBAP molecules from aqueous solution, respectively. Moreover, the impact of various operational / functional parameters such as packing density, the number of fibrous membrane, the module length, the membrane porosity / tortuosity, and ultimately the aqueous solution flow rate on the molecular separation efficiency of IP and 4-IBAP is studied in more details.

Developing and validating a dynamic model of water production by direct-contact membrane distillation

We consider the development and fitting of a dynamic model for desalinated water production by a direct-contact membrane distillation (DCMD) unit. Two types of dynamic-model structures, namely, lumped parameter and spatial, were evaluated. Both the models were validated using experimental response data generated by step testing the inlet hot stream temperature of a DCMD pilot plant. Both the model structures failed to follow the dynamic response adequately. However, a modification of the model by adding a heat loss term resulted in enhanced predictions for both model structures. The overall relative error in the model–plant mismatch was approximately 3%. This is reasonable considering the random uncertainties associated with the plant operation. This observation also improves our understanding of the importance of using better correlations for heat-transfer coefficients, to develop a more reliable and accurate predictive model for a wide range of operating conditions.

Molecular insights on NaCl crystal formation approaching PVDF membranes functionalized with graphene

Atomistic simulations of graphene–PVDF membranes speeding up NaCl crystal nucleation and growth in comparison to the pristine PVDF membranes.

Engineering Water and Solute Dynamics and Maximal Use of CNT Surface Area for Efficient Water Desalination

While polymer-based membranes and the consistent plants and elements have long been considered and optimized, there are only few studies on optimization of the new generation of carbon-based porous membranes for water desalination. By modeling the elements and their corresponding parameters in a vertical configuration via COMSOL Multiphysics software, an experimental setup was modified that contained various bare and carbon nanotube (CNT)-covered microprocessed porous membranes in parallel and in series. Several design parameters such as inlet pressure, length of outlet, vertical distance of the parallel membranes, and horizontal distances of the series membranes were optimized. Taking advantage of the uttermost surface area of CNTs and the engineered particle trajectory, almost 90% NaCl rejection and 97% Allura red rejection were obtained with very high permeation values. Considering microsized outlets, the results of particle rejections are outstanding owing to the smart design of the setup. The results of this work can be extended to larger and smaller scales up to the point where the governing equations still hold.

Development of a mechanistic model for prediction of CO2 capture from gas mixtures by amine solutions in porous membranes

A mechanistic model was developed in order to predict capture and removal of CO₂ from air using membrane technology. The considered membrane was a hollow-fiber contactor module in which gas mixture containing CO₂ was assumed as feed while 2-amino-2-metyl-1-propanol (AMP) was used as an absorbent. The mechanistic model was developed according to transport phenomena taking into account mass transfer and chemical reaction between CO₂ and amine in the contactor module. The main aim of modeling was to track the composition and flux of CO₂ and AMP in the membrane module for process optimization. For modeling of the process, the governing equations were computed using finite element approach in which the whole model domain was discretized into small cells. To confirm the simulation findings, model outcomes were compared with experimental data and good consistency was revealed. The results showed that increasing temperature of AMP solution increases CO₂ removal in the hollow-fiber membrane contactor.