Optimization of net power density in Reverse Electrodialysis

Reverse Electrodialysis (RED) extracts electrical energy from the salinity difference between two solutions using selective ion exchange membranes. In RED, conditions yielding a large net power density (NPD) are generally desired, due to the still large cost of the membranes. NPD depends on a large number of physical and geometric parameters. Some of these, for example the inlet concentrations of concentrate and diluate, can be regarded as “scenario” variables, imposed by external constraints (e.g., availability) or chosen by different criteria than NPD maximization. Others, namely the thicknesses HCONC, HDIL and the velocities UCONC, UDIL in the concentrate and diluate channels, can be regarded as free design parameters and can be chosen so as to maximize NPD. In the present study, a simplified model of a RED stack was coupled with an optimization algorithm in order to determine the conditions of maximum NPD in the space of the variables HCONC, HDIL,UCONC, UDIL for different sets of “scenario” variables. The study shows that an optimal choice of the free design parameters for any given scenario, as opposed to the adoption of standard fixed values for the same parameters, may provide significant improvements in NPD

Fluid-structure interaction and flow redistribution in membrane-bounded channels

The hydrodynamics of electrodialysis and reverse electrodialysis is commonly studied by neglecting membrane deformation caused by transmembrane pressure (TMP). However, large frictional pressure drops and differences in fluid velocity or physical properties in adjacent channels may lead to significant TMP values. In previous works, we conducted one-way coupled structural-CFD simulations at the scale of one periodic unit of a profiled membrane/channel assembly and computed its deformation and frictional characteristics as functions of TMP. In this work, a novel fluid-structure interaction model is presented, which predicts, at the channel pair scale, the changes in flow distribution associated with membrane deformations. The continuity and Darcy equations are solved in two adjacent channels by treating them as porous media and using the previous CFD results to express their hydraulic permeability as a function of the local TMP. Results are presented for square stacks of 0.6-m sides in cross and counter flow at superficial velocities of 1 to 10 cm/s. At low velocities, the corresponding low TMP does not significantly affect the flow distribution. As the velocity increases, the larger membrane deformation causes significant fluid redistribution. In the cross flow, the departure of the local superficial velocity from a mean value of 10 cm/s ranges between -27% and +39%

CFD Simulation of Mass Transfer Phenomena in Spacer Filled Channels for Reverse Electrodialysis Applications

Salinity Gradient Power via Reverse Electrodialysis is a topic of primary importance nowadays. It allows to get energy from the \u201ccontrolled\u201d mixing of solutions at different salt concentration. The performance of this technology depends on many factors such as: components properties (i.e. membranes, spacers, electrodes), stack geometry, operating conditions and feeds features. Concentration polarization phenomena may significantly affect the actual membrane potential, thus reducing the gross power produced. On the other hand, C-polarization phenomena may significantly be reduced by suitably choosing the hydrodynamic regime within the stack. Such a choice may in turn significantly require higher pumping power, thus reducing the net power output. In this work, carried out within the EU-FP7 funded REAPower project, CFD simulations were carried out in order to study the fluid flow behaviour and mass transport phenomena within spacer-filled channels for SGP-RE technology. The effect of different parameters (channel geometry, feed flow rate, feed solution concentration and current density) on concentration polarization was assessed. The well known unit cell approach was adopted for the simulations in order to reduce their computational requirements as well as to increase the level of detail. Results show that the electrical potential loss due to polarization phenomena should be regarded as little significant in the case of seawater-brine for the operating conditions and geometrical configurations investigated. Conversely, a great attention should be devoted to such phenomena when very diluted solutions are to be employed (e.g. river water)

Flow and mass transfer in spacer-filled channels for reverse electrodialysis: a CFD parametrical study

In reverse electrodialysis (RED) concentration polarization phenomena and pressure drop affect strongly the power output obtainable; therefore the channel geometry has a crucial impact on the system optimization. Both overlapped and woven spacers are commonly commercialised and adopted for RED experiments; the latter exhibit some potential advantages, such as better mixing and lower shadow effect, but they have been poorly investigated in the literature so far. In this work, computational fluid dynamics was used to predict fluid flow and mass transfer in spacer-filled channels for RED applications. A parametric analysis for different spacer geometries was carried out: woven (w) and overlapped (o) spacers with filaments at 90\ub0 were simulated, and Reynolds number, pitch to height ratio (l/h) and orientation with respect to the main flow (\u3b1=0\ub0 and \u3b1=45\ub0) were made to vary. The filament arrangement was found to be a crucial feature; for any given pumping power, higher Sherwood numbers were provided by the w-arrangement. The influence of flow attack angle and filament spacing depends on Reynolds number and filament arrangement. Only the configuration w-\u3b145 avoids the presence of poorly mixed zones near the wires. Among the cases investigated here, the configuration that provided the best mixing conditions was w, l/h=2, \u3b1=45\ub0

A multi-physics modelling tool for Reverse Electrodialysis

In this work, a multi-physics modelling approach has been developed for the RED process

Water desalination by capacitive electrodialysis: Experiments and modelling

Electrodialysis-related technologies keep spreading in multiple fields, among which water desalination still plays a major role. A new technology that has not yet been thoroughly investigated is capacitive electrodialysis (CED), which couples the standard ED with capacitive electrodes. CED has a number of advantages such as removal of toxic products and system simplification. Little mention is made of this technology in the literature and, to the best of our knowledge, no modelling works have ever been presented. In this work, the CED process has been studied through experiments and modelling. A CED model is presented for the first time. With a simple calibration based on macroscopic membrane properties and the characterisation of electrode behaviour, the model is able to simulate the dynamics of simple as well as more complex layouts. An original experimental characterisation of electrodes is presented, showing how the collected data can be implemented into the model. After a successful validation with experimental data, dynamic simulations of a single pass CED unit have been performed with the aim of assessing the effect of different capacitive electrode properties on process performance. Results show how the impact of these properties is different depending on the number of cell pairs

Hydrodynamics and mass transfer in straight fiber bundles with non-uniform porosity

The present study investigates the effects of non-uniformity in a bundle's porosity by considering a model channel made up of "dense" (low porosity) and "loose" (high porosity) regions. In a first, simplified, approach these regions are treated as non-interacting porous media and previously obtained computational results are used for the Darcy permeability and the Sherwood number. In a second, and more complete, approach 3-D CFD simulations are conducted for a checkerboard arrangement of alternately "dense" and "loose" regions with square-arrayed fibers, accounting for entry effects and for interactions between regions. Non-uniformity causes a significant increase of the permeability and a strong reduction of the Sherwood number. These effects are larger, approaching those obtained for non-interacting regions, if the regions' length scale is large. The attainment of fully developed conditions is greatly shifted forward in non-uniform bundles and the mass transfer development length may largely exceed the physical length of most hollow-fiber devices

Investigation of Reverse ElectroDialysis Units by Multi-Physical Modelling

Reverse electrodialysis (RED) is an electrochemical membrane process that converts the salinity gradient energy between two solutions into electric current, by using ion exchange membranes. A novel multi-physical approach for RED modelling is proposed. 2-D simulations of one cell pair with tertiary current distribution (Nernst\u2013Plank equation and local electroneutrality) were performed. Moreover, the Donnan exclusion theory was implemented for simulating double layer phenomena. Transport phenomena and electrochemical behavior were well described. The influence of membrane/channel configuration, dilute concentration and feeds velocity on the process performance was assessed. For a dilute concentration 64 0.01M, stacks with profiled membranes reached lower resistances and higher net powers (up to 4.4 W/m2) with respect to stacks with empty channels, thus suggesting that only in some cases the profiles lead to a performance enhancement

Performance Comparison of Alternative Hollow-Fiber Modules for Hemodialysis by Means of a CFD-Based Model

Commercial hemodialyzers are hollow-fiber cylindrical modules with dimensions and inlet– outlet configurations dictated mostly by practice. However, alternative configurations are possible, and one may ask how they would behave in terms of performance. In principle, it would be possible to depart from the standard counter-flow design, while still keeping high clearance values, thanks to the increase in the shell-side Sherwood number (Sh) due to the cross-flow. To elucidate these aspects, a previously developed computational model was used in which blood and dialysate are treated as flowing through two interpenetrating porous media. Measured Darcy permeabilities and mass transfer coefficients derived from theoretical arguments and CFD simulations conducted at unit-cell scale were used. Blood and dialysate were alternately simulated via an iterative strategy, while appropriate source terms accounted for water and solute exchanges. Several module configurations sharing the same membrane area, but differing in overall geometry and inlet–outlet arrangement, were simulated, including a commercial unit. Although the shell-side Sherwood number increased in almost all the alternative configurations (from 14 to 25 in the best case), none of them outperformed in terms of clearance the commercial one, approaching the latter (257 vs. 255 mL/min) only in the best case. These findings confirmed the effectiveness of the established commercial module design for the currently available membrane properties

Integrated modelling of membrane deformation, fluid dynamics and mass transfer in electromembrane processes

In recent years, water and energy supply issues have drawn the attention of the scientific community to electromembrane processes. Electrodialysis (ED) and Reverse Electrodialysis (RED) are two of the most attractive electromembrane technologies for water desalination and electric energy production from salinity gradients, respectively. In order to gain an important place in the industrial market, technological challenges on various aspects are involved in the optimization of these processes. In this context, profiled membranes exhibit interesting performance. However, the mechanical behavior of the membranes and its interaction with fluid dynamics has been poorly investigated so far. In membrane-based processes, a trans-membrane pressure (TMP) between the different solutions flowing through a module may be a design feature or may arise for various reasons (e.g. flow arrangement, differences in physical properties). This may lead to a local deflection of membranes due to their low mechanical stiffness. As a result, the channel geometry may be modified affecting flow and mass transfer characteristics. In this work, we developed an integrated model for the numerical simulation of local mechanical deformation and of fluid dynamics and associated mass transport phenomena inside deformed channels. Profiled membranes with Round Pillars were simulated under typical REDED working conditions. Membranes were assumed to be with perfectly linear elastic behaviour. 3-D simulations of a couple of membranes and of the interposed fluid were conducted by the unit cell approach (periodic domain). The Ansys Mechanical 18 (Workbench) and the Ansys CFX 18 software was used. The profiled membranes were investigated in both expanded and compressed configurations by varying the TMP of 0.1 bar steps from -0.4 to +0.4 bar. Then, CFD simulations of the deformed channels were performed. The influence of TMP was to increase friction under compression conditions (up to 2.2 times) and to reduce it under expansion conditions (i.e. up to 60%). Overall, compression enhanced mass transfer and expansion reduced it, but with smaller and more complex effects than on friction. Channel deformations affect largely fluid dynamics and mass transfer characteristics. Even mild trans-membrane pressures may produce significant variations in the overall process performance. Results will be implemented in a novel higher-scale simulation tool for studying the distribution of flow in whole channels