Next Article in Journal
Combined Separator Based on a Porous Ion-Exchange Membrane for Zinc–Halide Batteries
Next Article in Special Issue
Comparison of the Electrodialysis Performance in Tartrate Stabilization of a Red Wine Using Aliphatic and Aromatic Commercial and Modified Ion-Exchange Membranes
Previous Article in Journal
A Straightforward Method to Prepare MOF-Based Membranes via Direct Seeding of MOF-Polymer Hybrid Nanoparticles
Previous Article in Special Issue
Smart Superhydrophobic Filter Paper for Water/Oil Separation and Unidirectional Transportation of Liquid Droplet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Permeate Flux Enhancement in Air Gap Membrane Distillation Modules with Inserting Λ-Ribs Carbon-Fiber Open Slots

1
Department of Chemical and Materials Engineering, Tamkang University, Tamsui, New Taipei 251, Taiwan
2
Department of Water Resources and Environmental Engineering, Tamkang University, Tamsui, New Taipei 251, Taiwan
3
Department of Fundamental and Applied Sciences, HICoE-Centre for Biofuel and Biochemical Research, Institute of Self-Sustainable Building, Universiti Teknologi PETRONAS, Seri Iskandar 32610, Perak Darul Ridzuan, Malaysia
*
Author to whom correspondence should be addressed.
Membranes 2023, 13(1), 66; https://doi.org/10.3390/membranes13010066
Submission received: 16 December 2022 / Revised: 29 December 2022 / Accepted: 31 December 2022 / Published: 4 January 2023

Abstract

:
A novel design of an air gap membrane distillation (AGMD) module was proposed to enhance the permeate flux improvement for the desalination of pure water productivity. The modeling equations for predicting permeate flux in the AGMD module by inserting Λ -ribs carbon-fiber open slots under various hydrodynamic angles were developed theoretically and experimentally. The temperature distributions of both hot and cold feed streams were represented graphically with the hot saline flow rate, inlet saline temperature, and carbon-fiber hydrodynamic angles as parameters. The results showed a good agreement between the experimental results and theoretical predictions. Designed by inserting Λ -ribs carbon-fiber open slots into the flow channel, the membrane distillation module was implemented to act as an eddy promoter and yield an augmented turbulence flow. The effect of Λ -ribs carbon-fiber open slots not only assured the membrane stability by preventing vibration but also increased the permeate flux by diminishing the temperature polarization of the thermal boundary layer. The permeate flux improvement by inserting Λ -ribs carbon-fiber open slots in the AGMD module provided the maximum relative increment of up to 15.6% due to the diminution of the concentration polarization effect. The experimental data was incorporated with the hydrodynamic angle of Λ -ribs carbon-fiber open slots to correlate the enhancement factor with the Nusselt numbers to confirm the theoretical predictions. The accuracy derivation between the experimental results and theoretical predictions was pretty good, within 9.95 E 1.85 . The effects of operating and designing parameters of hot saline flow rate, inlet saline temperature, and hydrodynamic angle on the permeate flux were also delineated by considering both the power consumption increment and permeate flux enhancement.

1. Introduction

The membrane distillation (MD) process becomes one of the crucial issues in arid and semi-arid areas of fresh water scarcity [1,2] because of its simplicity and low-energy requirements [3], which makes the MD system an economic feasibility option. The desalination process is where the volatile species in the hot feed stream are vaporized and transported across the porous hydrophobic membrane in order to produce a high-purity liquid [4]. The vapor pressure gradient across the hydrophobic microporous membrane resulting from a trans-membrane temperature variance to drive the vapor flux permeating from the hot side to the cold side of condensates collected [5]. The existence of temperature differences plays a vital review in their features and device performance [6] and attracts global attention with considerable heat losses [7]. The temperature and concentration polarization effects were attributed to significant limits for the MD permeate flux performance and major assessment factors at industrial application scales [8,9]. AGMD is one of the most energy-efficient MD device among different types of MD modules [10], which can be operated in a relatively low temperature and low-pressure stream driven by low-grade heat from renewable energy sources [11,12]. However, the high heat loss builds up the temperature gradient inside the boundary layers of both membrane-surface sides [8] due to the temperature polarization effect with the significant thermal resistance occurrence. Identification of the magnitude of these effects in improving permeate flux enhancement on the MD system was proposed to reduce the polarization effects by using eddy promoters in the flow channel, including various approaches, e.g., spacers [13], filament [14,15], roughened surface [16], turbulence promoter [17], and porous and non-porous hollow fiber module [18]. Mitigation of the thermal boundary layer by adding eddy promoters could generate vortices and secondary flows in the flow stream, and thus, the permeate flux enhancement is achieved accordingly.
A new design, by inserting Λ-ribs carbon-fiber open slots in a flat-plate AGMD module, was developed theoretically and experimentally to improve distillate flux efficiency. The convective heat transfer coefficient and latent heat across the membrane were incorporated into the energy balance equations derived in the present study to solve the temperature fields, which depend on the turbulence intensity augmentation from inserting eddy promoters. However, the pressure drops with power consumption increments have been taken into account in the economic analysis [19] due to implementing eddy promoters in the flow channel for ultrafiltration and direct contact membrane distillation (DCMD) modules [20]. Many studies offer more insight based on mathematical models to develop strategies for improving permeate flux enhancement [21], and the heat and mass transport modeling approaches were reviewed methodically by Dong [22].
The present work focuses on the AGMD module to improve the device performance by inserting Λ-ribs carbon-fiber open slots into the flow channel of AGMD modules, which were found to disrupt the laminar boundary layer and enhance the local shear stress owing to enhancing the local surface shear stress, and thus, yielding augmentation of driving-force temperature gradients comes out with the increment of trans-membrane permeate flux [8,9]. The new construct of the module, by inserting Λ -ribs carbon-fiber open slots in the flow channel, was studied experimentally and numerically with its performance for saline water desalination, and was obtained in terms of a correlated expression [14,23] of Nusselt numbers. The main purpose of this study is to investigate the effects of designing and operating variables on the device performance, and the comparison was made between the modules with empty channels and the improved module with inserting Λ -ribs carbon-fiber open slots under various operating conditions. An alternative strategy for the membrane distillation modules improves the permeate flux by inserting Λ -ribs carbon-fiber open slots, and adjusting the open slots’ hydrodynamic angles in the flow channel of flat-plate AGMD modules, which disturbs the thermal boundary layer and yields a higher permeate flux improvement for enhanced heat and mass transfer. The suitable selection on economic feasibility was identified and explored theoretically by considering both the permeate flux improvement and power consumption increment.

2. Experimental Apparatus and Procedure

The experimental setup of the AGMD modules with/without (empty channel) Λ -ribs carbon-fiber open slots is presented schematically in Figure 1, and a photo of the present experimental apparatus is shown in Figure 2 with the acrylic plates as outside walls. The two acrylic plates of the module unit have three holes flowing in and out at both the entrance and exit ends, respectively, and the hot-feed and cold-feed channels are stacked together, as illustrated in Figure 3.
The hot and cold channel’s length, width, and height are 0.21 m, 0.29 m, and 2 mm, respectively. The hot feed saline water of 3.5 wt% NaCl and pure water were transported from two thermostats (G-50, 60 L, 3500 W, DENG YNG, Kaohsiung, Taiwan) through the flat-plate module using conventional pumps (51K40RA-A, ASTK, New Taipei, Taiwan) to regulate steadily two inlet streams at specified temperatures, respectively. The experimental runs were carried out for various hot feed temperatures (40, 45, 50, 55 °C), and both inlet and outlet temperatures were measured using thermometer probes (Type K/J, Tutron, New Taipei, Taiwan) connected to both streams of the flat-plate membrane modules. The permeate flux collected in the cold side was then weighted using an electronic balance (XS 4250C, Precisa Gravimetrics AG, Dietikon, Switerland) and recorded on the PC. The operation conditions of saline feed stream with various flow rates (6.67 × 10−6), 8.33 × 10−6, 11.7 × 10−6 and 15.0 × 10−6 m3/s) were adjusted by using flow meters (FE-091312-D (hot stream), FN-0423112-F (cold stream), Fong-Jei, Hsinchu, Taiwan) and controller (US-2000-40W, ASTK, New Taipei, Taiwan). Hydrophobic polytetrafluoroethylene (PTFE) membrane (J020A330R, Toyo Roshi Kaisha, Ltd., Tokyo, Japan) supported by polypropylene net (PP) with a nominal pore size of 0.2 μm, a porosity of 0.72, and a thickness of 130 μm was conducted in the present experimental runs. Two hydrodynamic angles and two carbon-fiber open slot’s widths were implemented into flow channels to generate vortices with changing direction along the hot saline water feed flows from the grid to grid, and the line carbon-fiber open slots were used in the permeate side to prevent the vibration of membrane acting as a supporter, as shown in Figure 3. The thermal conductivity of carbon-fiber open slots within 374 and 409 W/mK with negligible thermal resistance as compared to that of copper of 400 W/mK.

3. Theoretical Modeling of Mass and Heat Transfer in AGMD Modules

A representation of heat- and mass-transfer models were formulated in the AGMD module and depicted in Figure 4, while the schematic thermal transfer resistances in series are shown in Figure 5.
The permeate flux is dependent on the vapor pressure across the trans-membrane surfaces on the saline feed stream and air gap, as well as the membrane permeation coefficient ( c m ) [24,25] as
N = c T ( P 1 s a t P 3 s a t ) = c T d P d T | T m ( P 1 s a t P 3 s a t ) = c T P a v g λ M w R T m 2 ( P 1 s a t P 3 s a t )
where P 1 s a t   and P 3 s a t   are the saturated pressure on the membrane surface of the hot stream and the condensate film, respectively, and in which c T is as follows:
c T = ( 1 c m + 1 c a ) 1 = [ [ 1.064 ε r τ δ m ( M w R T m ) 1 2 ] 1 + [ 1 | Y m | l n D m ε δ m τ M w R T m ] 1 + [ 1 | Y a | l n D a δ a M w R T a ] 1 ] 1
where the tortuosity τ = 1 ε [26].
Theoretical modeling of both heat- and mass-transfer behaviors for a non-isothermal process inside the AGMD module is schematically illustrated in Figure 4. The heat transfer in the AGMD system caused by the temperature gradient across each component of the module. Using the equations for the conservation of enthalpy flow by considering both conduction heat flow and mass diffusion enthalpy flow through the composite system [27] as follows:
q = h h ( T h T 1 ) = k m δ m ( T 1 T 2 ) + N λ = k a δ a ( T 2 T 3 ) + N λ = h f ( T 3 T 4 ) = k p δ p ( T 4 T 5 ) = h c ( T 5 T c )
Equating the energy conservation of the heat flow and latent heat in each region for two intervals, ( T 1 T 2 ) and ( T 2 T 3 ) , yields the total heat transfer coefficient of the hot feed side and the surface of the condensate film as follows:
q = { ( δ m k m + δ a k a ) - 1 + [ c T ( ( 1 x NaCl ) ( 1 0.5 x NaCl 10 x NaCl 2 ) P w + P 3 ) λ 2 M w 2 R T a v g 2 ] } ( T 1 T 3 ) = H m ( T 1 T 3 )
in which
H m = ( δ m k m + δ a k a ) 1 + [ c T ( ( 1 x NaCl ) ( 1 0.5 x NaCl 10 x NaCl 2 ) P w + P 3 ) λ 2 M w 2 R T a v g 2 ]
where N λ is referred to the latent heat of vaporization and a w = 1 0.5 x NaCl 10 x NaCl 2 [4] is the activity coefficient. In addition, the thermal conductivity of the membrane k m can be determined by the thermal conductivities of vapor in the membrane pore k g and the solid membrane material k s by Warner [28] as:
k m = ε   k g + ( 1 ε ) k s
Similarly, the overall heat transfer coefficient of the cold feed side was obtained as
q = ( 1 h f + δ p k p + 1 h c ) - 1 ( T 3 T c ) = H c ( T 3 T c )
The bulk temperature T h of the hot saline feed stream drops to the membrane surface temperature T 1 , and thus a larger amount of heat was needed to vaporize water at the membrane surface owing to the noticeable temperature polarization effect. The heat transfer rate-limiting step in the thermal boundary layer due to the large difference in temperature between the bulk feed stream and the membrane surface was examined by the temperature polarization coefficient T P C , which is an indicator defined as the ratio of membrane surface temperatures’ gradient to bulk temperatures’ gradient and can be expressed in terms of heat transfer coefficients as follows:
T P C = ( T 1 T 3 ) ( T h T c ) = h h H c h h H c + h h H m + H c H m
Implementing Λ -ribs carbon-fiber open slots into the hot feed stream reduces the thickness of the thermal boundary layer at the membrane surface, as shown in Figure 6, and thus, T P C   increases (i.e., the reduction of temperature polarization effect).
The driving force for permeate flux decreases when the temperature gradient at the membrane surface on the hot feed side is less than the bulk temperature gradient. An attempt was proposed in the last two decades to augment turbulence intensity by inserting eddy promoters into the flow channel to enlarge the turbulence intensity. Restated, the reduction of the temperature polarization effect in decreasing the temperature gradient between the bulk stream and membrane surface indicates that the Λ -ribs carbon-fiber open slots could disturb the thermal boundary layer to strengthen a larger convective heat-transfer coefficient. It comes out with a thinner thermal boundary layer owing to increasing the larger T P C value, and thus, a higher permeate flux is achieved. Various hydrodynamic angles and carbon-fiber widths inserted into the flowing and empty channels (without inserting carbon fiber) were conducted in the experimental work under the same total coverage area, as shown in Figure 7, respectively. A proportion of the hot saline water feed flows across over the carbon-fiber filaments from the open slot to another open slot, and the direction change induced by the hydrodynamic angle due to following a zigzag-like pathway, as seen in Figure 8.
The longitudinal temperature distributions of both hot and cold feed streams were estimated by the one-dimensional modeling equations from the heat balance at steady state, which is presented in a finite fluid element of the energy-flow diagram, as shown in Figure 9.
d T h d z = q W Q h   ρ h C p , h = W Q h   ρ h C p , h H m T P C ( T h T c )
d T c d z = q W Q c   ρ c C p , c = W Q c   ρ c C p , c H m T P C ( T h T c )
The temperature distributions of both hot and cold feed streams were solved in the above two simultaneous first-order ordinary differential equations of Equations (9) and (10) with the use of the estimated convective heat-transfer coefficients and calculated iteratively of Figure 10 by marching the fourth-order Runge–Kutta method numerically along the length of the AGMD module. By following the similar flowchart of calculation procedure adapted from [29], and thus, the theoretical predictions of permeate fluxes and the permeate flux enhancement were obtained accordingly. The temperature distributions were predicted theoretically not only in the hot/cold bulk flows ( T h and T c ) but also on the membrane surfaces ( T 1 and T 3 ) of both hot and cold feed streams, respectively. Meanwhile, the convective heat-transfer coefficients were calculated and validated by the experimental results. Experimental runs were carried out in this study to demonstrate the device performances of permeate fluxes by using the module with inserting Λ -ribs carbon-fiber widths in flow channels and the module using the empty channel.

4. Permeate Flux Enhancement and Power Consumption Increment

The degree of heat transfer enhancement is commonly expressed by an enhancement factor, α E , which is the ratio of the heat transfer coefficients of the improved module by inserting Λ -ribs carbon-fiber open slots to that of the empty channel (without inserting carbon fiber). The Λ -ribs carbon-fiber open slots inserted in the conduit of the hot feed stream are implemented as compared to the module of the empty channel, and the permeate flux enhancement factor α E depends on the hydrodynamic angles and carbon-fiber widths [30] with respect to the characterization of the turbulence intensity as follows:
α E = N u E N u l a m
where
N u E = h h D h , h k   for   a   module   inserting   Λ - ribs   carbon - fiber   open   slots
N u l a m = 4.36 + 0.036 R e P r ( D h , h L ) 1 + 0.011 [ R e P r ( D h , h L ) ] 0.8   for   the   laminar   flow   in   a   flat   membrane   module .
Applying the regressed correlation approach and quantifying the enhanced permeate flux due to the eddy promoter of the module were solved directly through the permeate flux enhancement factor, which was based on dimensional analysis of Buckingham’s π theorem. Thus, the Nusselt number can be related to dimensionless groups to lump the influence of eddies and vortices created by the turbulent flow due to implementing the Λ -ribs carbon-fiber open slots as:
N u E = f   ( W e D h , h , sin   θ )
where W e and D h , h are the Λ -ribs carbon-fiber width and hydrodynamic angles of the hot stream-side, respectively.
The power consumption increment is required due to inserting Λ -ribs carbon-fiber open slots into the saline feed channel. The friction losses were calculated in determining the power consumption by using the Fanning friction factor f F [31]:
H i = m ˙ h w f , h + m ˙ c w f , c = Q h ρ h w f , h + Q c ρ c w f , c i = c a r b o n   f i b e r ,   e m p t y
w f , j = 2 f F , j v ¯ j 2 L D h , i , j = h , c
in which ( β = d W ) [32]:
f F , j = 24 ( 1 1.3553   β + 1.9467 β   2 1.7012 β   3 + 0.9564 β   4 0.2537 β   5 ) Re   j ,   j = h ,   c .
The power consumption increment I P due to the friction losses in the conduits can be readily derived as follows:
I p = H c a r b o n   f i b e r H e m p t y H e m p t y × 100 % .

5. Results and Discussion

The experimental data with empty channel (without implementing carbon-fiber open slots instead of using nylon fiber) and 2 mm and 3 mm Λ -ribs carbon-fiber widths were used to determine the correlation for the permeate enhancement factor α E , as expressed in Equation (19), which was regressed by setting up the normal equations for the least square parameters. The resultant regression analysis from curve-fitting with the squared correlation coefficient ( R 2 = 0.951 ), as shown in Figure 11.
N u E = 1.527 ( W e D h ) 0.427 ( S i n θ ) 0.504
Figure 12 shows the SEM micrographs of the fresh and used membranes of experimental runs. The SEM images indicated that the presence of Λ -ribs carbon-fiber open slots would not be a problem of fouling or scaling in conducting experimental runs.
Implementing carbon-fiber open slots play an important role in the permeate flux improvement due to disturbing the flow stream inside the concentration boundary layer, and thus, the higher turbulence intensity was augmented owing to reducing heat-transfer resistances. The effect of the carbon-fiber open slots in the AGMD module on the longitudinal temperature profiles of both hot and cold feed streams is shown in Figure 13, in which the temperature gradient tapered from the higher value at the entrance to the outlet streams. The temperature profiles show that temperature deviations (say temperature gradients between the membrane surface temperature and bulk temperature of both hot and cold feed streams) are considerably decreased while the AGMD modules with inserting Λ -ribs carbon-fiber open slots for both 120 ° and 90 ° hydrodynamic angles as compared to the module using the empty channel. The temperature polarization effect was reduced for a promising result by inserting Λ -ribs carbon-fiber open slots in the flow channel due to eddy turbulence increment and investigated by computational simulation [33]. One can find that the temperature gradient of inserting Λ -ribs carbon-fiber open slots of the 90 ° hydrodynamic angle was higher than that of the 120 ° hydrodynamic angle.
The temperature polarization coefficients T P C is defined in Equation (8) and calculated by both feed stream temperature distributions, the results investigated and compared the effects of hydrodynamic angles on the temperature polarization effect, as depicted in Figure 14. Inserting Λ -ribs carbon-fiber open slots of the 90 ° hydrodynamic angle was found to be a larger T P C value (a higher heat transfer rate) than those of the Λ -ribs carbon-fiber open slots with the 120 ° hydrodynamic angle and the empty channel as well; this is attributed to more carbon-fiber open slots to strengthen turbulence intensity and shrink the thermal boundary layer with a larger T P C value, and thus, comes out with a higher permeate flux through the hydrophobic membrane. The theoretical results show that a higher inlet temperature of the hot saline feed stream results in a lower T P C values owing to the higher permeate flux associated with the more latent heat of vaporization, T P C increased with a decreasing hydrodynamic angle as well. Improved devices with implementing Λ -ribs carbon-fiber open slots into the hot feed stream could use thinner the thermal boundary layer on the membrane surface due to the strengthened turbulence intensity, which comes out of a reduction of the temperature polarization effect, and thus, the heat transfer enhancements were achieved as compared to the module with using the empty channel.
The good agreement of the theoretical predictions with those obtained from experimental results was achieved. The accuracy deviation [34] of the experimental results from the theoretical predictions was calculated using the following definition as:
E ( % ) = 1 N e x p i = 1 N e x p | N t h e o N e x p | N e x p × 100
where N e x p , N t h e o , and N e x p are the number of experimental runs, theoretical predictions, and experimental results of the permeate fluxes, respectively. The accuracy deviations with two Λ -ribs carbon-fiber widths and two hydrodynamic angles were calculated; the agreement of experimental results deviated from theoretical predictions was pretty good, being within 9.95 E 1.85 .
Inserting Λ -ribs carbon-fiber open slots with two carbon-fiber widths and two hydrodynamic angles produced a higher turbulence intensity that resulted in a higher heat transfer rate, as well as a higher permeate flux. Comparisons were made on theoretical predictions and experimental results of permeate fluxes between the empty channel and the channels by inserting Λ -ribs carbon-fiber open slots of two carbon-fiber widths, say 2 mm and 3 mm, as shown in Figure 15a,b. In general, inserting Λ -ribs carbon-fiber open slots in the AGMD module is a more significant permeate transporting flux through the hydrophobic membrane with the smaller carbon-fiber widths and smaller hydrodynamic angle due to yielding a larger temperature driving-force gradient. Meanwhile, the heat transfer resistance dominating the permeate flux decreased with the increasing inlet volumetric flow rate and inlet saline feed temperature.
The theoretical predictions and experimental results of permeate flux were presented graphically in Figure 16 and Figure 17 for the empty channel and the channels with inserting 2 mm and 3 mm Λ -ribs carbon-fiber open slots under various inlet saline feed temperatures and hydrodynamic angles, respectively. The results show that the permeate flux increased with the saline volumetric flow rate and inlet saline feed temperature but with decreasing hydrodynamic angle and carbon-fiber width.
A percentage increment of permeate flux improvement I E employing a module with inserting Λ -ribs carbon-fiber open slots is best illustrated by calculating in comparisons of the module using the empty channel, as
I E = N c a r b o n   f i b e r N e m p t y N e m p t y × 100 %
The permeate flux enhancements I E of the module with inserting Λ -ribs carbon-fiber open slots into the saline feed stream were significantly achieved in Table 1 for 2 mm and 3 mm carbon-fiber widths, respectively. The theoretical predictions were conducted in comparisons with various hydrodynamic angles, carbon-fiber widths, inlet saline feed temperatures, and inlet volumetric flow rate as parameters. The theoretical predictions show that the permeate flux improvement up to 15.6% is obtained by inserting 2 mm Λ -ribs carbon-fiber open slots as compared to that in the empty channel device, as seen in Table 1. The extent of permeate flux increment is more significant in the smaller hydrodynamic angle and the higher inlet saline feed temperature. Overall, the performance of permeate flux is enhanced by inserting Λ -ribs carbon-fiber open slots into saline feed channel, which plays the important role of eddy promoters in the pressure-driven membrane distillation processes and water treatment technologies.
Although inserting Λ-ribs carbon-fiber open slots into the flow channel of AGMD modules could create resistance for the feed circulation resulting in the power consumption increment, the permeate flux enhancement was achieved, as well as due to enhancing the local surface shear stress with the increment of trans-membrane permeate flux. This study further examines the permeate flux performance by evaluating the desirable permeate flux increment to the undesirable friction loss increment, say I E / I P , due to inserting Λ -ribs carbon-fiber open slots in making the suitable selection from an economic viewpoint. The effects of carbon-fiber widths, hydrodynamic angles, inlet saline feed temperatures, and saline volumetric flow rates are shown in Figure 18 and Table 2, as referred to the ratio of I E / I P .
The value of I E / I P increased with increasing saline volumetric flow rate, which indicates that the expenses of energy consumption increment could compensate by the permeate flux enhancement due to utilizing the driving-force temperature gradient more effectively. Furthermore, the ratio of I E / I P of the module with an inserting 2 mm carbon-fiber width was higher than that of the channel with an inserting 3 mm carbon-fiber width consistent with the trend of the permeate flux obtained with the smaller carbon-fiber widths and smaller hydrodynamic angle, and presented in Table 2. Meanwhile, the increase in the saline volumetric flow rate yields a higher ratio of I E / I P , and the order of the ratio of I E / I P was expected, with the same trend of the permeate flux.

6. Conclusions

There are two innovation points in this module design, which are easy fabricating of Λ -ribs carbon-fiber open slots and a larger ratio of permeate flux increment to power consumption increment. This is the value of the present study. Two hydrodynamic angles of turbulence promoters were used for the AGMD module and compared to an empty channel under two carbon-fiber widths in the present study. Increasing shear stress in the thermal boundary layer on the membrane surface due to inserting turbulence promoters could considerably reduce the temperature polarization effect. The results show that the permeate flux increases with both saline feed flow rate and inlet saline feed temperature, and demonstrate its feasibility of permeate flux enhancement up to 15.6% by inserting Λ -ribs carbon-fiber open slots in the flow channel of the AGMD system. The comparisons of the permeate flux enhancement with inserting turbulence promoters under various carbon-fiber widths and hydrodynamic angles on the permeate flux in AGMD modules were drawn to the following conclusions:
  • The permeate flux increases with decreasing hydrodynamic angle and carbon-fiber width in the AGMD module by inserting Λ -ribs carbon-fiber open slots.
  • The permeate flux enhancement increases with the increase of the volumetric flow rate by inserting Λ -ribs carbon-fiber open slots where the enhancement of the AGMD module using 2 mm open slots is higher than that of the 3 mm one.
  • The higher the inlet saline feed temperature yields, the higher the permeate flux enhancement.
  • A higher permeate flux enhancement was found in the module with inserting Λ -ribs carbon-fiber open slots compared to that in the empty-channel module under both the smaller hydrodynamic angle and carbon-fiber width.
  • The higher ratio of I E / I P was achieved for the module by inserting Λ -ribs carbon-fiber open slots, and was operated under both the higher saline feed flow rate and inlet saline feed temperature, and both the smaller carbon-fiber width and hydrodynamic angle.
A new design in this study includes the advantage of this membrane distillation module in strengthening the turbulence intensity as an alternative strategy for the permeate flux in the AGMD module. The alternative configurations of inserting Λ -ribs carbon-fiber open slots showed a great potential to considerably diminish the temperature polarization effect, and thus the permeate flux was enhanced.

Author Contributions

C.-D.H.: funding acquisition, conceptualization, methodology, data collection and analysis. L.C.: funding acquisition, writing—original draft, analysis and interpretation of data. Y.-L.Y.; analysis of data, writing—review, editing. S.-T.C.: data curation, investigation. J.W.L.: investigation, writing—review, editing. Z.-Z.C.: data curation, formal analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Science and Technology Council (NSTC) grant number (MOST 104-2815-C-032-038-E).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors wish to thank the National Science and Technology Council (NSTC) of the Republic of China (Taiwan) for the financial support.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

A membrane area (m2)
a w water activity in NaCl solution
C p , c heat capacity of cold fluid (J kg−1K−1)
C p , h heat capacity of hot fluid (J kg−1K−1)
c a mass transfer coefficient of air gap region (kg m−2 Pa−1 s−1)
c m mass transfer coefficient of membrane (kg m−2 Pa−1 s−1)
c T overall mass transfer coefficient of membrane (kg m−2 Pa−1 s−1)
d thickness of empty channel (m)
D equivalent hydraulic diameter of empty channel (m)
D h , h equivalent hydraulic diameter of hot side (m)
D h , c equivalent hydraulic diameter of cold side (m)
Edeviation of experimental results from the theoretical predictions
f F Fanning friction factor
h c convection coefficient of cold stream (W m−2 K−1)
h h convection coefficient of saline stream (W m−2 K−1)
H i hydraulic dissipate energy (J kg−1),   i = p r o m o t e r , e m p t y
H m thermal convection coefficient of membrane (W m−2 K−1)
I E increased percentage of permeate flux
I P raised percentage of hydraulic loss
L axial distance (m)
k a thermal conductivity of the air gap (W m−1 K−1)
k m thermal conductivity of the membrane (W m−1 K−1)
k f thermal conductivity coefficient of aqueous solution (W m−1 K−1)
k g thermal conductivity coefficient of the vapor in the membrane pore (W m−1 K−1)
k s thermal conductivity coefficient of the solid membrane material (W m−1 K−1)
k m thermal conductivity of the membrane (W m−1 K−1)
k p thermal conductivity of cooling plate (W m−1 K−1)
w f friction loss of conduits (J kg−1)
M W molecular weight of water (kg mol−1)
m ˙ mass flow rate (kg s−1)
N permeate flux (kg m−2 s−1)
N u Nusselt number
N u E enhanced Nusselt number with inserting turbulence promoter
N u l a m Nusselt number for laminar flow
P pressure (Pa)
P 1 s a t saturation vapor pressure in the saline feed stream (Pa)
P 3 s a t saturation vapor pressure in the hot feed flow side (Pa)
P w s a t saturated vapor pressure of pure water (Pa)
P r Prandtl number
q heat transfer flux (W)
q heat transfer rate (W/m2)
Q volumetric flow rate (m3 s−1)
R gas constant (8.314 J mol−1 K−1)
R e Reynolds number
r radius of membrane pore (m)
T temperature (°C)
T m mean temperature in membrane (°C)
T P C temperature polarization coefficients
v ¯ average velocity (m s−1)
W e Λ -ribs carbon-fiber width
| Y m | n natural log mean mole fraction of air
x N a C l liquid mole fraction of NaCl
x w liquid mole fraction of water
z axial coordinate along the flow direction (m)
Greek letters
α E heat-transfer enhancement factor
β aspect ratio of the channel
Δ P vapor pressure difference of membrane (Pa)
δ m thickness of membrane (µm)
ε membrane porosity
ε e channel voidage
λ latent heat of water (J/kg)
μ fluid viscosity (kg s−1 m−1)
ρ density (kg m−3)
T P C temperature polarization coefficients
Subscripts
1 on membrane surface in cold feed side
2 on membrane surface in hot feed side
ain the air gap
cin the cold feed flow channel
fin the condensate film
hin the hot feed flow channel
min the membrane
pin the cooling plate
carbon fiberinserting Λ -ribs carbon-fiber open slots
emptymodule with inserting nylon fiber as supporters
expexperimental results
inat the inlet
lammodule with empty channel
outat the outlet
theotheoretical predictions
Superscripts
Ethe channel with inserting Λ -ribs carbon-fiber open slots

References

  1. Alkhudhiri, A.; Darwish, N.; Hilal, N. Treatment of saline solutions using Air Gap Membrane Distillation: Experimental study. Desalination 2013, 323, 2–7. [Google Scholar] [CrossRef]
  2. Elhenawy, Y.; Elminshawy, N.A.S.; Bassyouni, M.; Alhathal Alanezi, A.; Drioli, E. Experimental and theoretical investigation of a new air gap membrane distillation module with a corrugated feed channel. J. Membr. Sci. 2020, 594, 117461. [Google Scholar] [CrossRef]
  3. Dehesa-Carrasco, U.; Perez-Rabago, C.A.; Arancibia-Bulnes, C.A. Experimental evaluation and modeling of internal temperatures in an air gap membrane distillation unit. Desalination 2013, 326, 47–54. [Google Scholar] [CrossRef]
  4. Lawson, K.W.; Lloyd, D.R. Membrane distillation. II. direct contact MD. J. Membr. Sci. 1996, 120, 123–133. [Google Scholar] [CrossRef]
  5. Alklaibi, A.M.; Lior, N. Membrane-distillation desalination: Status and potential. Desalination 2005, 171, 111–131. [Google Scholar] [CrossRef]
  6. Olatunji, S.O.; Camacho, L.M. Heat and Mass Transport in Modeling Membrane Distillation Configurations: A Review. Front. Energy Res. 2018, 6, 130. [Google Scholar] [CrossRef]
  7. Al-Obaidani, S.; Curcio, E.; Macedonio, F.; Di Profio, G.; Al-Hinai, H.; Drioli, E. Potential of membrane distillation in sea-water desalination: Thermal efficiency, sensitivity study and cost estimation. J. Membr. Sci. 2008, 323, 85–98. [Google Scholar] [CrossRef]
  8. Schofield, R.; Fane, A.; Fell, C. Heat and mass transfer in membrane distillation. J. Membr. Sci. 1987, 33, 299–313. [Google Scholar] [CrossRef]
  9. Martinez-Diez, L.; Vázquez-González, M. Effects of polarization on mass transport through hydrophobic porous membranes. Ind. Eng. Chem. Res. 1998, 37, 4128–4135. [Google Scholar] [CrossRef]
  10. Alkhudhiri, A.; Darwish, N.; Hilal, N. Membrane distillation: A compressive review. Desalination 2012, 287, 2–18. [Google Scholar] [CrossRef]
  11. Gryta, M.; Tomaszewska, M.; Grzechulska, J.; Morawski, A. Membrane distillation of NaCl solution containing natural organic matter. J. Membr. Sci. 2001, 181, 279–287. [Google Scholar] [CrossRef]
  12. Drioli, E.; Ali, A.; Macedonio, F. Membrane distillation: Recent developments and perspectives. Desalination 2015, 356, 56–84. [Google Scholar] [CrossRef]
  13. Chang, H.; Hsu, J.-A.; Chang, C.-L.; Ho, C.-D.; Cheng, T.-W. Simulation study of transfer characteristics for spacer-filled membrane distillation desalination modules. Appl. Energy 2015, 185, 2045–2057. [Google Scholar] [CrossRef]
  14. Phattaranawik, J.; Jiraratananon, R.; Fane, A. Effects of net-type spacers on heat and mass transfer in direct contact mem-brane distillation and comparison with ultrafiltration studies. J. Membr. Sci. 2003, 217, 193–206. [Google Scholar] [CrossRef]
  15. Santos, J.; Geraldes, V.; Velizarov, S.; Crespo, J. Investigation of flow patterns and mass transfer in membrane module channels filled with flow-aligned spacers using computational fluid dynamics (CFD). J. Membr. Sci. 2007, 305, 103–117. [Google Scholar] [CrossRef]
  16. Ho, C.; Chang, H.; Chang, C.; Huang, C. Theoretical and experimental studies of performance enhancement with roughened surface in direct contact membrane distillation desalination. J. Membr. Sci. 2013, 433, 160–166. [Google Scholar] [CrossRef]
  17. Yang, X.; Yu, H.; Wang, R.; Fane, A.G. Analysis of the effect of turbulence promoters in hollow fiber membrane distillation modules by computational fluid dynamic (CFD) simulations. J. Membr. Sci. 2012, 415–416, 758–769. [Google Scholar] [CrossRef] [Green Version]
  18. Singh, D.; Sirkar, K.K. Desalination by air gap membrane distillation using a two hollow-fiber-set membrane module. J. Membr. Sci. 2012, 421–422, 172–179. [Google Scholar] [CrossRef]
  19. Da Costa, A.R.; Fane, A.G.; Wiley, D.E. Spacer characterization and pressure drop modelling in spacer-filled channels for ultrafiltration. J. Membr. Sci. 1994, 87, 79–98. [Google Scholar] [CrossRef]
  20. Ho, C.-D.; Chen, L.; Lai, J.-Y.; Ng, C.A. Theoretical and experimental studies of direct contact membrane distillation modules with inserting W-shaped carbon-fiber spacers. Desalin. Water Treat. 2017, 71, 32–44. [Google Scholar] [CrossRef]
  21. Yang, X.; Wang, R.; Fane, A.G. Novel designs for improving the performance of hollow fiber membrane distillation modules. J. Membr. Sci. 2011, 384, 52–62. [Google Scholar] [CrossRef] [Green Version]
  22. Dong, Y.; Dai, X.; Zhao, L.; Gao, L.; Xie, Z.; Zhang, J. Review of Transport Phenomena and Popular Modelling Approaches in Membrane Distillation. Membranes 2021, 11, 122. [Google Scholar] [CrossRef] [PubMed]
  23. Lee, J.-G.; Kim, Y.-D.; Kim, W.-S.; Francis, L.; Amy, G.; Ghaffour, N. Performance modeling of direct contact membrane distillation (DCMD) seawater desalination process using a commercial composite membrane. J. Membr. Sci. 2015, 478, 85–95. [Google Scholar] [CrossRef] [Green Version]
  24. Lawson, K.; Lloyd, D. Membrane distillation. J. Membr. Sci. 1997, 124, 1–25. [Google Scholar] [CrossRef]
  25. Schofield, R.W.; Fane, A.G.; Fell, C.J.D. Gas and vapor transport through microporous membranes. I. Knudsen-Poiseuille transition. J. Membr. Sci. 1990, 53, 159–171. [Google Scholar] [CrossRef]
  26. Iversen, S.; Bhatia, V.; Dam-Johansen, K.; Jonsson, G. Characterization of microporous membranes for use in membrane contactors. J. Membr. Sci. 1997, 130, 205–217. [Google Scholar] [CrossRef]
  27. Izquierdo-Gil, M.; Garcia-Payo, M.; Pineda, C.F. Air gab membrane distillation of sucrose aqueous solutions. J. Membr. Sci. 1999, 155, 291–307. [Google Scholar] [CrossRef]
  28. Warner, S.B. Fiber Science; Princeton-Hall: Englewood Cliffs, NJ, USA, 1995. [Google Scholar]
  29. Ho, C.D.; Chen, L.; Lim, J.W.; Lin, P.-H.; Lu, P.-T. Distillate flux enhancement of direct contact membrane distillation modules with inserting cross-diagonal carbon-fiber spacers. Membranes 2021, 11, 973. [Google Scholar] [CrossRef]
  30. Phattaranawik, J.; Jiraratananon, R.; Fane, A. Heat transport and membrane distillation coefficients in direct contact membrane distillation. J. Membr. Sci. 2003, 212, 177–193. [Google Scholar] [CrossRef]
  31. Welty, J.R.; Wicks, C.E.; Wilson, R.E. Fundamentals of Momentum, Heat, and Mass Transfer, 3rd ed.; John Wiley & Sons: New York, NY, USA, 1984. [Google Scholar]
  32. Kakac, S.; Shah, R.K.; Aung, W. Handbook of Single-Phase Convective Heat Transfer; Wiley: New York, NY, USA, 1987. [Google Scholar]
  33. Shakaib, M.; Hasani, S.; Ahmed, I.; Yunus, R.M. A CFD study on the effect of spacer orientation on temperature polarization in membrane distillation modules. Desalination 2012, 284, 332–340. [Google Scholar] [CrossRef]
  34. Moffat, R.J. Describing the uncertainties in experimental results. Exp. Therm. Fluid Sci. 1988, 1, 3–17. [Google Scholar] [CrossRef]
Figure 1. Experimental configuration of the air gap membrane distillation (AGMD) system: (A) hot fluid thermostat; (B) cold fluid thermostat; (C) pump; (D) flow meter; (E) beaker; (F) electronic balance; (G) temperature probe; (H) temperature indicator; (I) the AGMD module.
Figure 1. Experimental configuration of the air gap membrane distillation (AGMD) system: (A) hot fluid thermostat; (B) cold fluid thermostat; (C) pump; (D) flow meter; (E) beaker; (F) electronic balance; (G) temperature probe; (H) temperature indicator; (I) the AGMD module.
Membranes 13 00066 g001
Figure 2. A photo of the experimental apparatus of the AGMD system.
Figure 2. A photo of the experimental apparatus of the AGMD system.
Membranes 13 00066 g002
Figure 3. Components of the inserting Λ -ribs carbon-fiber open slots of the AGMD module. (1) Acrylic plate; (2) silicon rubber sealing; (3) Λ -ribs carbon-fiber open slots; (4) PTFE membrane; (5) line carbon-fiber open slots; (6) silicon rubber sealing; (7) aluminum plate; (8) silicon rubber sealing; and (9) acrylic plate.
Figure 3. Components of the inserting Λ -ribs carbon-fiber open slots of the AGMD module. (1) Acrylic plate; (2) silicon rubber sealing; (3) Λ -ribs carbon-fiber open slots; (4) PTFE membrane; (5) line carbon-fiber open slots; (6) silicon rubber sealing; (7) aluminum plate; (8) silicon rubber sealing; and (9) acrylic plate.
Membranes 13 00066 g003
Figure 4. Schematic thermal boundary layers and temperature profiles of the AGMD module.
Figure 4. Schematic thermal boundary layers and temperature profiles of the AGMD module.
Membranes 13 00066 g004
Figure 5. Schematic heat transfer resistances in the AGMD module.
Figure 5. Schematic heat transfer resistances in the AGMD module.
Membranes 13 00066 g005
Figure 6. Schematic mass and heat transfer resistances in the AGMD module.
Figure 6. Schematic mass and heat transfer resistances in the AGMD module.
Membranes 13 00066 g006
Figure 7. Components of inserting various Λ -ribs hydrodynamic angles and carbon-fiber widths.
Figure 7. Components of inserting various Λ -ribs hydrodynamic angles and carbon-fiber widths.
Membranes 13 00066 g007
Figure 8. Saline feed flow streamlines in the Λ -ribs carbon-fiber channel with hydrodynamic angles.
Figure 8. Saline feed flow streamlines in the Λ -ribs carbon-fiber channel with hydrodynamic angles.
Membranes 13 00066 g008
Figure 9. The energy balance made within a finite fluid element.
Figure 9. The energy balance made within a finite fluid element.
Membranes 13 00066 g009
Figure 10. Flow chart for solving membrane surface temperatures and heat transfer coefficients.
Figure 10. Flow chart for solving membrane surface temperatures and heat transfer coefficients.
Membranes 13 00066 g010
Figure 11. Comparisons of both calculated and experimental Nusselt numbers.
Figure 11. Comparisons of both calculated and experimental Nusselt numbers.
Membranes 13 00066 g011
Figure 12. SEM micrographs of the fresh and used membranes of experimental runs.
Figure 12. SEM micrographs of the fresh and used membranes of experimental runs.
Membranes 13 00066 g012
Figure 13. The effect of various hydrodynamic angles with a 2 mm carbon-fiber width on temperature profiles ( Q h = 8.33 × 10 6 m 3 s ,   Q c = 1.5 × 10 5 m 3 s ).
Figure 13. The effect of various hydrodynamic angles with a 2 mm carbon-fiber width on temperature profiles ( Q h = 8.33 × 10 6 m 3 s ,   Q c = 1.5 × 10 5 m 3 s ).
Membranes 13 00066 g013
Figure 14. The effect of various hydrodynamic angles of 2 mm carbon-fiber width on T P C   .
Figure 14. The effect of various hydrodynamic angles of 2 mm carbon-fiber width on T P C   .
Membranes 13 00066 g014
Figure 15. Effects of various hydrodynamic angles and carbon-fiber widths on permeate fluxes.
Figure 15. Effects of various hydrodynamic angles and carbon-fiber widths on permeate fluxes.
Membranes 13 00066 g015
Figure 16. The effects of inlet saline feed temperatures and hydrodynamic angles on permeate fluxes with inserting 2 mm carbon-fiber widths.
Figure 16. The effects of inlet saline feed temperatures and hydrodynamic angles on permeate fluxes with inserting 2 mm carbon-fiber widths.
Membranes 13 00066 g016
Figure 17. The effects of inlet saline feed temperatures and hydrodynamic angles on permeate fluxes with inserting 3 mm carbon-fiber widths.
Figure 17. The effects of inlet saline feed temperatures and hydrodynamic angles on permeate fluxes with inserting 3 mm carbon-fiber widths.
Membranes 13 00066 g017
Figure 18. The effects of hydrodynamic angles and inlet saline feed temperatures on the value of I E / I P .
Figure 18. The effects of hydrodynamic angles and inlet saline feed temperatures on the value of I E / I P .
Membranes 13 00066 g018
Table 1. The effects of hydrodynamic angles and carbon-fiber widths on permeate flux enhancement.
Table 1. The effects of hydrodynamic angles and carbon-fiber widths on permeate flux enhancement.
T h ,   i n ( ° C ) Q h × 10 6 ( m 3   s 1 ) Empty Channel2 mm3 mm
90 ° 120 ° 90 ° 120 °
N t h e o × 10 3 kg   m 2   s 1 N t h e o × 10 3 kg   m 2   s 1 I E N t h e o × 10 3 kg   m 2   s 1 I E N t h e o × 10 3 kg   m 2   s 1 I E N t h e o × 10 3 kg   m 2   s 1 I E
406.670.400.437.800.436.420.425.370.424.35
8.330.480.5310.20.528.400.526.850.515.11
11.70.530.5911.90.5810.30.578.010.565.94
15.00.550.6312.70.6211.40.608.770.596.35
456.670.540.598.280.587.070.575.660.574.77
8.330.660.7310.70.729.190.717.660.705.62
11.70.720.8112.70.8011.00.788.910.776.66
15.00.760.8713.90.8612.40.839.480.827.20
506.670.690.758.740.747.760.736.260.725.00
8.330.840.9411.50.9210.10.918.560.896.09
11.70.921.0514.01.0311.81.0110.00.997.21
15.00.971.1214.81.1013.11.0810.61.057.68
556.670.830.929.930.908.230.896.960.885.30
8.331.021.1512.71.1310.61.129.021.096.68
11.71.121.2914.61.2612.41.2410.51.217.72
15.01.191.3715.61.3513.61.3211.01.298.16
Table 2. The effects of hydrodynamic angles and carbon-fiber widths on I E / I P .
Table 2. The effects of hydrodynamic angles and carbon-fiber widths on I E / I P .
T h ,   i n ( C ) Q h × 10 6 ( m 3   s 1 ) 2 mm3 mm
90 ° 120 ° 90 ° 120 °
I E / I P I E / I P I E / I P I E / I P
406.670.180.150.130.11
8.330.230.200.160.13
11.70.270.240.190.15
15.00.290.260.210.16
456.670.190.170.140.12
8.330.250.220.180.14
11.70.290.260.210.17
15.00.310.280.220.18
506.670.200.190.150.13
8.330.270.240.200.15
11.70.320.280.240.18
15.00.330.300.250.19
556.670.230.200.170.14
8.330.290.250.220.17
11.70.330.290.250.19
15.00.350.320.260.20
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ho, C.-D.; Chen, L.; Yang, Y.-L.; Chen, S.-T.; Lim, J.W.; Chen, Z.-Z. Permeate Flux Enhancement in Air Gap Membrane Distillation Modules with Inserting Λ-Ribs Carbon-Fiber Open Slots. Membranes 2023, 13, 66. https://doi.org/10.3390/membranes13010066

AMA Style

Ho C-D, Chen L, Yang Y-L, Chen S-T, Lim JW, Chen Z-Z. Permeate Flux Enhancement in Air Gap Membrane Distillation Modules with Inserting Λ-Ribs Carbon-Fiber Open Slots. Membranes. 2023; 13(1):66. https://doi.org/10.3390/membranes13010066

Chicago/Turabian Style

Ho, Chii-Dong, Luke Chen, Yan-Ling Yang, Shih-Ting Chen, Jun Wei Lim, and Zheng-Zhong Chen. 2023. "Permeate Flux Enhancement in Air Gap Membrane Distillation Modules with Inserting Λ-Ribs Carbon-Fiber Open Slots" Membranes 13, no. 1: 66. https://doi.org/10.3390/membranes13010066

APA Style

Ho, C. -D., Chen, L., Yang, Y. -L., Chen, S. -T., Lim, J. W., & Chen, Z. -Z. (2023). Permeate Flux Enhancement in Air Gap Membrane Distillation Modules with Inserting Λ-Ribs Carbon-Fiber Open Slots. Membranes, 13(1), 66. https://doi.org/10.3390/membranes13010066

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop