Next Article in Journal
Free Energy Analyses of Cell-Penetrating Peptides Using the Weighted Ensemble Method
Next Article in Special Issue
Recent Mitigation Strategies on Membrane Fouling for Oily Wastewater Treatment
Previous Article in Journal
Electrospun Magnetic Nanocellulose–Polyethersulfone-Conjugated Aspergillus oryzae Lipase for Synthesis of Ethyl Valerate
Previous Article in Special Issue
Optimization of a High-Performance Poly(diallyl dimethylammonium chloride)-alumina-perfluorooctanoate Intercalated Ultrafiltration Membrane for Treating Emulsified Oily Wastewater via Response Surface Methodology Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distillate Flux Enhancement of Direct Contact Membrane Distillation Modules with Inserting Cross-Diagonal Carbon-Fiber Spacers

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 2021, 11(12), 973; https://doi.org/10.3390/membranes11120973
Submission received: 13 November 2021 / Revised: 3 December 2021 / Accepted: 6 December 2021 / Published: 9 December 2021
(This article belongs to the Special Issue Special Issue in Honor of Professor Ahmad Fauzi Ismail)

Abstract

:
A new design of direct-contact membrane distillation (DCMD) modules with cross-diagonal carbon-fiber spacers of various hydrodynamic angles in flow channels to promote turbulence intensity was proposed to enhance pure water productivity. Attempts to reduce the temperature polarization coefficient were achieved by inserting cross-diagonal carbon-fiber spacers in channels, which create wakes and eddies in both heat and mass transfer behaviors to enhance the permeate flux enhancement. A simplified equation was formulated to obtain the theoretical predictions of heat transfer coefficients in the current DCMD device. The permeate fluxes and temperature distributions of both hot and cold feed streams are represented graphically with the inlet volumetric flow rate and inlet temperature of the hot saline feed stream as parameters. The higher distillate flux of countercurrent-flow operations for saline water desalination was accomplished as compared to the concurrent-flow operations of various hydrodynamic angles. The results show that the agreement between the theoretical predictions and experimental results is reasonably good. The effects of countercurrent-flow operations and inserting carbon fiber spacers have confirmed technical feasibility and device performance enhancement of up to 45%. The influences of operating and design parameters on the pure water productivity with the expense of energy consumption are also discussed.

1. Introduction

The advantages of membrane distillation (MD) systems to produce pure water by using low grade thermal energy [1] in remote villages or rural areas [2,3] are their simplicity and low operating cost. The DCMD module has been recognized as a potential technology for desalination, solution concentration and wastewater treatment [4]. The DCMD module is the separation process to vaporize the volatile species in the hot feed stream, and the permeate flux collected in the cold feed stream, in which the vapor pressure difference creates a driving force across the hydrophobic membrane surfaces [4] and yields high-purity water [5]. The temperature polarization effect [6] caused the thermal driving-force reduction and transmembrane permeate flux decrement in most previous studies of DCMD modules. Modification of hot fluid channels with various strategies diminish the thermal boundary layer and minimize the temperature polarization effect in DCMD modules, as shown in Figure 1. The mass transfer rate enhancement was accomplished by employing spacer-filled channels [7,8], filaments [9], and rough-surface channels [10], as well as inserting shell side baffles, wavy shape fibers and fibers with a gear-shaped cross section in hollow fiber modules [11] to increase permeate flux up to 30–300%. An alternative configuration improves the permeate flux by inserting cross-diagonal carbon-fiber spacers and adjusting hydrodynamic angles in channels of flat-plate DCMD modules in this study. Mathematical modeling equations were developed to analyze the device performance of the module by inserting cross-diagonal carbon-fiber spacers in channels, and thus, a correlated expression [12] of Nusselt numbers was obtained and validated by the experimental results. The optimal selection’s economic feasibility was investigated theoretically under both concurrent and countercurrent-flow operations. Our objective herein is to determine a permeate flux assessment under various operating conditions. The cross-diagonal carbon-fiber spacers act as eddy promoters to disturb the thermal boundary layer on the hot feed stream, and lead to a trade-off between permeate flux improvement and energy consumption increment. The suitable selection of hydrodynamic angles on the economic analysis for device performance enhancements was identified and explored.

2. Theoretical Modeling of DCMD Modules

Theoretical modeling equations of both heat and mass transfer behaviors for a DCMD module, as shown in Figure 1, were investigated to predict the permeate flux at the membrane/liquid interface of the hot saline feed stream, and diffused through porous hydrophobic membranes, and then condensed at the membrane/liquid interface in the cold stream as distillate flux. Three types of net-like cross-diagonal carbon-fiber spacers with various hydrodynamic angles were implemented in the hot feed side to promote the mass transfer rate, and comparisons of device performance were made between both modules with/without inserting cross-diagonal carbon-fiber spacers.
As hot saline water feed flows from grid to grid, a proportion of the feed stream will change direction by the hydrodynamic angle and follow a zigzag-like pathway, as seen in Figure 2.
Heat and mass transfer models were formulated according to the following assumptions:
(a)
Steady-state operations;
(b)
Physical properties of fluid, frame plates, and membrane are constants;
(c)
Stagnant air within the membrane pore;
(d)
Mass transfer by diffusion and heat transfer by conduction associated with latent heat through the hydrophobic membrane;
(e)
No water transporting through the hydrophobic membrane;
(f)
Good insulation on the entire circumference of modules.
The non-isothermal process inside the DCMD module builds up the temperature gradient to enforce the permeate flux transferring across the hydrophobic membrane, which was condensed in the cold fluid stream as the pure water product. Theoretical modeling of both heat and mass transfer behaviors for a DCMD module was schematically illustrated in Figure 3. Mass-transfer modeling is needed to make balances of permeate flux by vapor diffusion, and then, the enthalpy flow conservation including heat conduction was formulated simultaneously in Equations (1)–(3) with the above assumptions as follows:
q h = h h ( T h T 2 ) the   hot   saline   water   feed   region  
q m = N λ + k m ( T 2 T 1 ) / δ m the   membrane   region
q c = h c ( T 1 T c ) the   cooling   water   region
where N λ is referred to as the latent heat of vaporization and k m ( T 2 T 1 ) / δ m is the conductive heat transfer, 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 is defined, following Warner [13], as:
k m = ε k g + ( 1 ε ) k s
The membrane permeation coefficient ( c m ) and the transmembrane saturation vapor pressure difference ( Δ P ) have been used extensively in mass transfer analysis of permeate flux for membrane distillation processes [2,14] as:
N = c m Δ P = c m [ P 2 s a t ( T 2 ) P 1 s a t ( T 1 ) ] = c m d P d T | T m ( T 2 T 1 ) = c m P m λ M w R T m 2 ( T 2 T 1 )
where P 1 s a t and P 2 s a t are the saturated pressure of water vapor on both membrane surfaces, respectively.
The combinations of the heat flow and latent heat of Equation (2) yields the overall heat transfer coefficient of the membrane as follows:
q m = N λ + k m δ m ( T 2 T 1 ) = ( c m [ a w ( 1 x NaCl ) P 2 + P 1 ] λ 2 M w 2 R T m 2 + k m δ m ) ( T 2 T 1 ) = H m ( T 2 T 1 )
where the membrane permeation coefficient c w is the addition of Knudsen diffusion and Poiseuille flow, a w = 1 0.5 x NaCl 10 x NaCl 2 is the activity coefficient [5] and the tortuosity τ = 1 / ε can be estimated using the porosity of the membrane [15].
The temperature polarization coefficient T P C is an indicator to indicate the extent of the thermal boundary-layer resistance which governs the distillate flux through the membrane. It is used to define as the ratio of membrane surface temperatures’ gradient to bulk temperatures’ gradient as follows:
T P C = ( T 2 T 1 ) / ( T h T c )
Manipulating and solving Equations (1), (3), and (6) by equating all heat transfer regions under steady-state operations, say q = q h = q m = q c and neglecting the heat loss on the outside of the DCMD module leads to the following:
T h = T 2 + H m h h ( T 2 T 1 )
T c = T 1 H m h c ( T 2 T 1 )
Equation (9) is subtracted from Equation (8) to give:
T h T c = ( T 2 T 1 ) + H m h h ( T 2 T 1 ) + H m h c ( T 2 T 1 ) = ( 1 + H m h h + H m h c ) ( T 2 T 1 )
Moreover, thus an alternative form of Equation (7) for T P C expressed in terms of heat transfer coefficients leads to the following:
T P C = h h h c h h h c + h h H m + h c H m
The procedure for calculating of theoretical values of both membrane surface temperatures ( T 1 and T 2 ) and the heat transfer coefficient will be described as follows. First, with the given operation conditions, the heat transfer coefficient is determined from Equations (8) and (9). Next, with the known inlet and outlet temperatures of both hot and cold streams, a temporary value of T 1 (or T 2 ) is estimated from Equation (8) once T 2 (or T 1 ) is assumed in Equation (9). Further, the convective heat transfer coefficient is calculated from Equation (5), with this calculated value of the convective heat transfer coefficient, new values of T 1 and T 2 are then recalculated from Equations (8) and (9). If the calculated values of T 1 and T 2 are different from the assumed values, continuous calculation by iteration is needed until the last assumed values of membrane surface temperatures meet the finally calculated values within a given convergence tolerance, as shown on the right-hand side of Figure 4.
The one-dimensional modeling equations of the energy balances were obtained by making the energy-flow diagram presented in a finite fluid element, as shown in Figure 5, to solve the longitudinal temperature distributions of both hot and cold feed streams as:
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 ) concurrent - flow   operations
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 ) countercurrent - flow   operations
The temperature distributions of both hot and cold feed streams were solved in the above two simultaneous ordinary differential equations of Equations (12) and (13) for concurrent-flow operation (or Equation (14) for countercurrent-flow operation) with the use of the estimated convective heat transfer coefficients, and calculated iteratively in the left-hand side of Figure 4 by marching the fourth-order Runge–Kutta method numerically along the length of the DCMD module, and thus, the theoretical permeate flux and permeate flux enhancement were obtained. The temperature distributions were predicted theoretically not only in the hot/cold bulk flows and on the membrane surfaces of both hot and cold feed streams under concurrent- and countercurrent-flow operations, respectively. Comparisons were made between the channel with cross-diagonal carbon-fiber spacers and the device with the empty channel.

3. Experimental Apparatus and Procedures

A detailed configuration schematic of an acrylic parallel-plate channel of length 21 cm, width 29 cm, and 2 mm of each cold and hot stream is illustrated in Figure 6. The hydrophobic membrane surfaces were supported by inserting cross-diagonal carbon-fiber spacers and by winding a 0.1 mm nylon fiber in the hot saline and cold feed sides, respectively, to prevent the membrane bending and wrinkling.
A commercial membrane, made of polytetrafluoroethylene (PTFE) supported by a polypropylene net (PP), was used in the experiments. The principal characteristics specified by the manufacturer (J020A330R, Toyo Roshi Kaisha, Ltd., Tokyo, Japan) are nominal pore size of 0.2 µm, porosity of 0.72, and thickness of 130 µm. PTFE is a hydrophobic membrane of many unique properties with four components, i.e., calcium fluoride, hydrofluoric acid, chloroform and water. The average molecular weight of the membranes ranges from 400,000 to 9,000,000. The manufacturing process involves the synthesis and polymerization of tetrafluoroethylene (TFE) with a series of chemical reactions to create the final product PTFE membrane which is exceptionally resistant to corrosion. The pore size and porosity of the hydrophobic composite membrane made of PTFE/PP will affect the permeation flux. However, the conductance of the distillate flux collected should be monitored and measured during the experimental runs, and it was less than 1.5 μs/cm in the present work. The cross-diagonal carbon-fiber spacers were implemented with various hydrodynamic angles in channels to generate vortices around those net-like carbon-fiber open slots, as shown in Figure 1 and Figure 2. The structure of the spacer is shown by the picture in Figure 7. A silicon rubber with thickness of 2 mm was glued on the acrylic plate to build up a spacer channel and also to prevent leakage.
Experimental runs were conducted under various operating conditions to study the device performances of permeate fluxes for two modules with/without inserting cross-diagonal carbon-fiber spacers in channels. The artificial saline water of 3.5 wt% NaCl was prepared by adding inorganic salts NaCl into distilled water. The experiments were operated by controlling various flow rates (0.1, 0.2, 0.3, 0.4 L/min) for various inlet hot saline temperatures (40, 45, 50, 55 °C) and 25 °C of the cold stream, and read by thermocouples (TM-946, Lutron, New Taipei, Taiwan). The temperatures of both streams were regulated by the thermostat (Water Bath G-50, DENG YNG, Taiwan) and the thermostat (Water Bath D650, DENG YNG, Taiwan), respectively. The distillate flux condensed in the cold side was then collected and weighed using an electronic balance (XS 4250C, Precisa Gravimetrics AG, Dietikon, Switzerland) to measure the distillate flux and recorded on the PC.
Figure 8 presents the SEM micrographs of the fresh membrane and the used membrane after experimental runs. The SEM images indicated that some salts were stuck on the membrane surface, but most of porous channels were not jammed by the salts.

4. Flux Enhancement Factor and Power Consumption Increment

The enhancement factor α E depends on various hydrodynamic angles, compared to the empty channel, and was correlated to calculate the augmented convective heat transfer coefficients in DCMD modules with implementing cross-diagonal carbon-fiber spacers [16] as:
α E = N u E / N u l a m
where
N u E = h h D h , h k for   a   module   with   inserted   cross - diagonal   carbon - fiber   spacers
N u l a m = 4.36 + 0.036 Re Pr ( D h , h / L ) 1 + 0.011 [ Re Pr ( D h , h / L ) ] 0.8 for   a   module   with   an   empty   channel
A better interpretation of both heat and mass transfer behaviors in the module with cross-diagonal carbon-fiber spacers could be described by a new method based on dimensional analysis of the Buckingham’s π theorem, which expresses the influence of eddies and vortices created by the turbulent flow due to the eddy promoter, and generates more turbulence intensity. The Nusselt number of flow channels with cross-diagonal carbon-fiber spacers can be related to dimensionless groups:
N u E = f   ( W e D h , h , sin θ )
where W e and D h , h are the carbon-fiber spacer width and hydraulic diameter of the hot stream side, respectively. The average velocity [17] and equivalent hydraulic diameter of cold and hot stream sides are defined as follows:
ν ¯ h = Q h d D ε e , ν ¯ c = Q c d W
D h , h = 4 ε e ( 2 / d ) + ( 1 ε e ) S v s p ,   D h , c = 4 d W 2 ( d + W )
S v s p = s u r f a c e v o l u m e = S s p V s p = 2 ( W e + d p ) W e d p
Re h = ρ h v ¯ h d h , h μ h ,   Re c = ρ c v ¯ c d h , c μ c
The power consumption increment is required due to inserting cross-diagonal carbon-fiber spacers into the saline water feed channel as eddy promoters. The friction losses to walls of both hot and cold streams were assumed to be significant, which were calculated in determining the power consumption by using Fanning friction factor f F [18]:
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 = carbon   fiber ,   empty  
w f , j = 2 f F , j v ¯ j 2 L d h , i ,   j = h , c
in which ( β = d / W ) [19]:
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 C = 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 estimated values of membrane surface temperatures and the convective heat transfer coefficients were obtained using a numerical flowchart in Figure 4, which were plugged into Equations (16) and (17) (Equation (18) for countercurrent-flow operations) in solving numerically by the fourth-order Runge–Kutta method of both bulk temperature distributions in hot/cold feed streams along the flowing direction of DCMD module. The theoretical predictions of permeate flux and permeate flux enhancement for the module with inserting cross-diagonal carbon-fiber spacers were thus obtained, and compared to the device with the empty channel. The experimental results with empty channel and 2 and 3 mm carbon-fiber slots were used to regress the correlation for the enhancement factor α E , as expressed in Equation (18). The resultant expression was determined via a regression analysis from curve-fitting with the squared correlation coefficient ( R 2 = 0.952 ), as shown in Figure 9.
α E = N u E N u l a m = 3.163 exp ( W e d h , h ) 0.766 sin θ 0.112
The effect of the carbon-fiber spacers on the longitudinal temperature profiles of both hot and cold feed streams in the DCMD module was shown in Figure 10 with hydrodynamic angles of 120 ° as an illustration. The temperature profiles show that both membrane surface temperatures of both hot and cold feed streams with/without net-like cross-diagonal carbon-fiber spacers are considerably different for both concurrent- and countercurrent-flow operations.
The permeate flux is proportional to the temperature gradient between both membrane surface temperatures T 1 and T 2 in the DCMD system. The temperature gradient between both membrane surfaces is higher in the carbon-fiber spacers channel than that in the empty channel. Reduction of the temperature polarization effect was achieved using carbon-fiber spacers in the channel for a promising result investigated by computational simulation [20]. The more significant temperature gradient results in a larger heat flux of the device with net-like cross-diagonal carbon-fiber spacers; hence more permeate flux or pure water productivity was predicted, as indicated in Equation (5). One can find that the temperature gradient of the carbon-fiber spacer width of 2 mm is higher than that of the 3 mm spacer. The temperature gradient appears a nonuniform profile and tapers from the higher value at the entrance to the outlet in concurrent-flow operations during the comparatively uniform temperature gradient of the countercurrent-flow operations. The descending heat transfer rate and permeate flux along the channel for concurrent-flow operations is thus confirmed as compared to a higher heat transfer rate and permeate flux in countercurrent-flow operations.
The devices with net-like cross-diagonal carbon-fiber spacers were used for promoting eddies and for temperature polarization reduction. This study investigated and compared the effects of carbon-fiber spacer widths on temperature polarization, as depicted in Figure 11.
The theoretical predictions of T P C show that a higher inlet temperature of the hot saline feed stream leads to a lower T P C for both concurrent- and countercurrent-flow operations, because the higher permeate flux needs more latent heat of vaporization, which results in a lower temperature gradient across the membrane surface. Implementing cross-diagonal carbon-fiber spacers into the hot feed stream intensifies turbulence, reduces the thickness of thermal boundary layer at the membrane surface, and increases T P C (i.e., lower temperature polarization). The above procedure leads to the reduction in the thickness of thermal boundary-layer thickness and increment in heat transfer rate, compared to the module without the spacer. In addition, inserting 2 mm cross-diagonal carbon-fiber spacers does have a larger T P C value (a higher heat transfer rate) than the 3 mm cross-diagonal carbon-fiber spacers and the device with the empty channel, which could result from more net-like opening slots to induce turbulence. Furthermore, a higher T P C value (say less thermal resistance) was achieved in countercurrent-flow operations than that in the DCMD system’s concurrent-flow operations. Restated, the reduction of temperature polarization effect in decreasing the temperature gradient between the bulk stream and membrane surface, say the thinner thermal boundary layer with a larger T P C value, and thus, yields a higher permeate flux through the hydrophobic membrane.
The agreement between the experimental results and theoretical predictions is in good consistency, where the accuracy deviation of the theoretical predictions from the experimental results are within an acceptable range, as indicated in Table 1a,b with the definition as follows:
E ( % ) = 1 N i = 1 N exp | N t h e o N exp | N exp × 100
Comparisons were made on theoretical predictions and experimental results of permeate fluxes between the empty channel and the channels with 2 mm cross-diagonal carbon-fiber spacers, as shown in Figure 12a,b as well as in Table 1a,b for illustrations. The results show that the permeate flux increases with the increase of the inlet volumetric flow rate, inlet saline temperature and hydrodynamic angles, and the extent of permeate flux increment is more significant in countercurrent-flow operations.
The theoretical predictions and experimental results of permeate flux were presented graphically in Figure 13 and Figure 14 for the empty channel and the channels with 2 mm and 3 mm inserting cross-diagonal carbon-fiber spacers, respectively. The permeate flux increases with the volumetric flow rate, hydrodynamic angle and the inlet saline temperature. Notice that the effect of hydrodynamic angle on the permeate flux concludes that more permeate flux with the use of the turbulence promoter by inserting cross-diagonal carbon-fiber spacers.
The improvement of the device performance I E from inserting the cross-diagonal carbon-fiber spacers is best illustrated by the percentage increase in the permeate flux compared to the device with an empty channel.
I E = N p r o m o t e r N e m p t y N e m p t y × 100 %
The permeate flux enhancements I E of the device with cross-diagonal carbon-fiber spacers in the hot saline stream were significant achievements under both concurrent- and countercurrent-flow operations, as demonstrated in Table 2a,b for 2 mm and 3 mm cross-diagonal carbon-fiber spacers, respectively. Both experimental results and the theoretical predictions were conducted in comparisons with various hydrodynamic angles, carbon-fiber spacer widths, inlet saline temperatures and inlet volumetric flow rate as parameters. Effects of hydrodynamic angle and the width of carbon-fiber spacer show that the permeate flux increases with the hydrodynamic angle but decreases with the spacer width. Overall, the performance of permeate flux is enhanced by inserting the cross-diagonal carbon-fiber spacers into the channel, which serves as eddy promoters in both concurrent- and countercurrent-flow operations. The results indicate that the permeate flux enhancement is up to 45%. Hence, the device with insertion of cross-diagonal carbon-fiber spacers may also be applied to the pressure-driven membrane distillation processes and water-treatment technologies.
The present work extends the previous study except for inserting cross-diagonal carbon-fiber spacers instead of W-shaped carbon-fiber spacers [21]. The graphical representation for comparisons with theoretical predictions of the permeate flux obtained in the present study and W-shaped carbon-fiber spacers [21] illustrates why the present configurations of inserting cross-diagonal carbon-fiber spacers is preferred, presented in Figure 15 for concurrent-flow operations.
Figure 16 presents the dependence of the Nusselt number on the Reynolds number. The figure shows that the Nusselt number increases as the Reynold number and the inlet saline temperature increase under both concurrent- and countercurrent-flow operations. Notably, the insertion of the diagonal carbon-fiber spacers significantly increases the Nusselt number and the convective heat transfer rate for both 2 and 3 mm open-slot widths. Despite the effect on the flow pattern, the change in the spacer width from 2 to 3 mm only leads to a moderate effect on the change of the Nusselt number.
This study further examines the device performance by evaluating the desirable permeate flux increment to undesirable power consumption increment I E / I P ratio due to inserting cross-diagonal carbon-fiber spacers in the flow channel associated with more flow resistance. The effects of flow configuration, carbon-fiber spacer widths and volumetric flow rate on I E / I P are shown in Figure 17. The value of I E / I P increases with increasing volumetric flow rate, which reveals that the expenses of energy consumption increment can compensate by the permeate flux enhancement. The value of I E / I P in countercurrent-flow configuration with more significant temperature gradient is higher than that in concurrent-flow configuration due to utilizing the driving-force temperature gradient more effectively. The ratio of I E / I P of the channel with 3 mm carbon-fiber spacers is higher than that of the channel with 2 mm carbon-fiber spacers. In other words, inserting 3 mm carbon-fiber spacer gives a higher value of I E / I P , which reflects that a more effective operation in increasing the permeate flux at the expense of energy consumption is expected. Although the permeate flux and the Nusselt number of the channel with 2 mm carbon-fiber spacers in Figure 13, Figure 14 and Figure 16 are higher than that of the channel with 3 mm carbon-fiber spacers, the energy consumption of the former channel is also higher. Comparisons of the economic feasibility among DCMD modules with inserting different widths of cross-diagonal carbon-fiber spacers were examined under both design and operating conditions.

6. Conclusions

A parallel-plate direct contact membrane distillation module with net-like cross-diagonal carbon-fiber spacers to enhance the permeate flux was investigated theoretically and experimentally. The theoretical predictions of the permeate flux enhancement by inserting turbulence promoters of cross-diagonal carbon-fiber spacers for various hot feed flow rates, inlet saline temperature and the carbon-fiber spacer widths under both concurrent- and countercurrent-flow operations were examined, and the correlated expression of Nusselt number was obtained as well. Comparisons of the permeate flux enhancements were made and the following conclusions were drawn:
  • The permeate flux increases with the increase of the volumetric flow rate.
  • Higher inlet saline temperature yields higher permeate flux productivity.
  • The permeate flux enhancement is obtained by inserting net-like cross-diagonal carbon-fiber spacers where the enhancement of the 2 mm slot opening is higher than that of the 3 mm one.
  • The permeate flux increases with increasing hydrodynamic angle in the slot of the cross-diagonal carbon-fiber spacers.
  • A maximum of 45.1% permeate flux enhancement was found in the device with cross-diagonal carbon-fiber spacers compared to that in the empty channel device under countercurrent-flow operations of 120° hydrodynamic angle.
  • A more considerable permeate flux was achieved in countercurrent-flow operations than in concurrent-flow operations due to the larger temperature gradient for countercurrent-flow operations.
  • The economic consideration of I E / I P for permeate flux enhancement to power consumption increment concluded that the power utilization is more effective for the channel with cross-diagonal carbon-fiber spacers in higher hot saline water flow rate, and the ratio of I E / I P for the 3 mm slot opening is higher than that of the 2 mm one.
A new contribution of this study is the desirable effect of raising the turbulence intensity as an alternative strategy [10] on the permeate flux in a direct contact membrane distillation module. The advantage of this membrane distillation device is it is somewhat easier to implement the experimental setup and has lower production costs. The alternative configurations of carbon-fiber spacers require further investigations to derive the optimal ratio between the permeate flux enhancement and energy consumption increment to make good economic sense.

Author Contributions

Conceptualization and methodology, C.-D.H.; data curation, P.-H.L.; writing—original draft, L.C. and P.-T.L.; writing—review and editing, J.-W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors wish to thank the Ministry of Science and Technology (MOST) of the Republic of China (Taiwan) for the financial support.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

a w Water activity in NaCl solution
C p , c Heat capacity of cold fluid (J kg−1 K−1)
C p , h Heat capacity of hot fluid (J kg−1 K−1)
c m Mass transfer coefficient of membrane (kg m−2 Pa−1 s−1)
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)
d Height of flow channel (m)
dpHeight of carbon-fiber spacers
EDeviation of experimental results from the theoretical predictions
f F Fanning friction factor
h c Convection coefficient of cold fluid (W m−2 K−1)
h h Convection coefficient of hot fluid (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 Thermal conductivity coefficient of hot saline feed (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)
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 h−1)
N u Nusselt number
N u E Nusselt number of the turbulence promoter
N u l a m Dimensionless Nusselt number for laminar flow
P Pressure (Pa)
P 1 s a t Saturation vapor pressure in the cold feed flow side (Pa)
P 2 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)
Pr Prandtl number
q Heat transfer rate (W/m2)
q c Heat transfer rate between cooling plate and cold fluid (W/m2)
q h Heat transfer rate between hot fluid and membrane surface (W/m2)
q m Heat transfer rate between membrane surface of hot fluid and air gap (W/m2)
Q Volumetric flow rate (m3 s−1)
R Gas constant (8.314 J mol−1 K−1)
ReReynolds number
T Temperature (°C)
T m Mean temperature in membrane (°C)
v ¯ Average velocity (m s−1)
W Width of flow channel
W e Width of carbon-fiber spacers
| 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
1Membrane surface on cold feed side
2Membrane surface on hot feed side
hIn the hot feed flow channel
cIn the cold feed flow channel
carbon fiberInserting cross-diagonal carbon-fiber spacers
emptyInserting nylon fiber as supporters
expExperimental results
inInlet
lamEmpty channel
outOutlet
theoTheoretical predictions
Superscripts
EThe channel with inserting cross-diagonal carbon-fiber spacers

References

  1. Drioli, E.; Ali, A.; Macedonio, F. Membrane distillation: Recent developments and perspectives. Desalination 2015, 356, 56–84. [Google Scholar] [CrossRef]
  2. Lawson, K.W.; Lloyd, D.R. Membrane distillation. J. Membr. Sci. 1997, 124, 1–25. [Google Scholar] [CrossRef]
  3. Martínez-Díez, L.; Vázquez-González, M.; Florido-Díaz, F. Study of membrane distillation using channel spacers. J. Membr. Sci. 1998, 144, 45–56. [Google Scholar] [CrossRef]
  4. Alkhudhiri, A.; Darwish, N.; Hilal, N. Membrane distillation: A comprehensive review. Desalination 2012, 287, 2–18. [Google Scholar] [CrossRef]
  5. Lawson, K.W.; Lloyd, D.R. Membrane distillation. II. direct contact MD. J. Membr. Sci. 1996, 120, 123–133. [Google Scholar] [CrossRef]
  6. Martínez-Díez, L.; Vázquez-González, M.I. Effects of polarization on mass transport through hydrophobic porous membranes. Ind. Eng. Chem. Res. 1998, 37, 4128–4135. [Google Scholar] [CrossRef]
  7. Phattaranawik, J.; Jiraratananon, R.; Fane, A. Effects of net-type spacers on heat and mass transfer in direct contact membrane distillation and comparison with ultrafiltration studies. J. Membr. Sci. 2003, 217, 193–206. [Google Scholar] [CrossRef]
  8. 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]
  9. Wang, T.; Zhang, Z.; Ren, X.; Qin, L.; Zheng, D.; Li, J. Direct observation of single- and two-phase flows in spacer filled membrane modules. Sep. Purif. Technol. 2014, 125, 275–283. [Google Scholar] [CrossRef]
  10. Ho, C.D.; Chang, H.; Chang, C.L.; Huang, C.H. 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]
  11. Yang, X.; Yu, H.; Wang, R.; Fane, A.G. Optimization of microstructured hollow fiber design for membrane distillation applications using CFD modeling. J. Membr. Sci. 2012, 421–422, 258–270. [Google Scholar] [CrossRef] [Green Version]
  12. 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]
  13. Warner, S.B. Fiber Science; Princeton-Hall: Englewood Cliffs, NJ, USA, 1995. [Google Scholar]
  14. 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]
  15. Iversen, S.B.; Bhatia, V.K.; Dam-Jphasen, K.; Jonsson, G. Characterization of microporous membranes for use in membrane contactors. J. Membr. Sci. 1997, 130, 205–217. [Google Scholar] [CrossRef]
  16. Phattaranawik, J.; Jiraratananon, R.; Fane, A.G. Heat transport and membrane distillation coefficients in direct contact membrane distillation. J. Membr. Sci. 2003, 212, 177–193. [Google Scholar] [CrossRef]
  17. Phattaranawik, J.; Jiraratananon, R.; Fane, A.G.; Halim, C. Mass flux enhancement using spacer filled channels in direct contact membrane distillation. J. Membr. Sci. 2001, 187, 193–201. [Google Scholar] [CrossRef]
  18. 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]
  19. Kakac, S.; Shah, R.K.; Aung, W. Handbook of Single-Phase Convective Heat Transfer; Wiley: New York, NY, USA, 1987. [Google Scholar]
  20. Shakaib, M.; Hasani, S.M.F.; 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]
  21. 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. Desalination Water Treat. 2017, 71, 32–44. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Fabrication structure and components of DCMD modules with cross-diagonal carbon-fiber spacers.
Figure 1. Fabrication structure and components of DCMD modules with cross-diagonal carbon-fiber spacers.
Membranes 11 00973 g001
Figure 2. Saline flow streamlines in the spacer-filled flat-plate channel with hydrodynamic angle.
Figure 2. Saline flow streamlines in the spacer-filled flat-plate channel with hydrodynamic angle.
Membranes 11 00973 g002
Figure 3. Schematic thermal boundary layers and temperature profiles of a DCMD module.
Figure 3. Schematic thermal boundary layers and temperature profiles of a DCMD module.
Membranes 11 00973 g003
Figure 4. Flow chart for solving membrane surface temperatures and heat transfer coefficients.
Figure 4. Flow chart for solving membrane surface temperatures and heat transfer coefficients.
Membranes 11 00973 g004
Figure 5. The energy balance made within a finite fluid element.
Figure 5. The energy balance made within a finite fluid element.
Membranes 11 00973 g005
Figure 6. Schematic diagram of the experimental apparatus for the DCMD system.
Figure 6. Schematic diagram of the experimental apparatus for the DCMD system.
Membranes 11 00973 g006
Figure 7. A photo of the experimental setup.
Figure 7. A photo of the experimental setup.
Membranes 11 00973 g007
Figure 8. SEM micrographs of the fresh membrane and the used membrane after experimental runs. (a) The fresh membrane; (b) the used membrane.
Figure 8. SEM micrographs of the fresh membrane and the used membrane after experimental runs. (a) The fresh membrane; (b) the used membrane.
Membranes 11 00973 g008
Figure 9. Comparison of estimated and experimental Nusselt numbers.
Figure 9. Comparison of estimated and experimental Nusselt numbers.
Membranes 11 00973 g009
Figure 10. Effect of net-like cross-diagonal carbon-fiber spacers on temperature profiles.
Figure 10. Effect of net-like cross-diagonal carbon-fiber spacers on temperature profiles.
Membranes 11 00973 g010aMembranes 11 00973 g010b
Figure 11. Effect of net-like cross-diagonal carbon-fiber spacers on T P C .
Figure 11. Effect of net-like cross-diagonal carbon-fiber spacers on T P C .
Membranes 11 00973 g011aMembranes 11 00973 g011b
Figure 12. Effects of hydrodynamic angles and operation types on permeate flux (2 mm spacers).
Figure 12. Effects of hydrodynamic angles and operation types on permeate flux (2 mm spacers).
Membranes 11 00973 g012
Figure 13. Effects of inlet saline temperatures on permeate flux (2 mm).
Figure 13. Effects of inlet saline temperatures on permeate flux (2 mm).
Membranes 11 00973 g013
Figure 14. Effects of inlet saline temperatures on permeate flux (3 mm).
Figure 14. Effects of inlet saline temperatures on permeate flux (3 mm).
Membranes 11 00973 g014
Figure 15. Comparisons of theoretical predictions and experimental results of permeate flux two type of carbon-fiber spacers (3 mm; concurrent-flow operations).
Figure 15. Comparisons of theoretical predictions and experimental results of permeate flux two type of carbon-fiber spacers (3 mm; concurrent-flow operations).
Membranes 11 00973 g015
Figure 16. Comparisons of theoretical Nusselt numbers of two operations. Dependence of Nu on Re (solid symbols: 60 °C, open symbols: 45 °C).
Figure 16. Comparisons of theoretical Nusselt numbers of two operations. Dependence of Nu on Re (solid symbols: 60 °C, open symbols: 45 °C).
Membranes 11 00973 g016
Figure 17. Effects of flow patterns on the value of I E / I P .
Figure 17. Effects of flow patterns on the value of I E / I P .
Membranes 11 00973 g017
Table 1. (a) Comparison of theoretical and experimental permeate fluxes for concurrent flow. (b) Comparison of theoretical and experimental permeate fluxes for countercurrent flow.
Table 1. (a) Comparison of theoretical and experimental permeate fluxes for concurrent flow. (b) Comparison of theoretical and experimental permeate fluxes for countercurrent flow.
T h ,   i n
( ° C )
Q h × 10 6
(m3 s−1)
(a) Cross-Diagonal Carbon-Fiber Spacers (2 mm)
60 ° 90 ° 120 °
N exp × 10 3
kg m−2 s−1
N t h e o × 10 3
kg m−2 s−1
E (%) N exp × 10 3
kg m−2 s−1
N t h e o × 10 3
kg m−2 s−1
E (%) N exp × 10 3
kg m−2 s−1
N t h e o × 10 3
kg m−2 s−1
E (%)
456.670.690.746.330.690.758.180.720.775.77
8.330.800.833.900.800.855.090.800.807.54
11.70.890.978.120.930.985.600.991.001.01
15.01.001.098.081.061.103.821.111.121.62
506.670.880.968.370.890.988.560.910.998.78
8.331.051.104.641.101.121.371.151.140.94
11.71.291.300.691.361.323.491.421.355.23
15.01.441.461.251.461.480.941.561.504.51
556.671.091.208.901.101.229.871.151.269.26
8.331.291.397.671.331.458.011.401.475.24
11.71.541.688.181.591.717.271.681.754.06
15.01.801.905.531.821.935.742.061.984.00
606.671.411.559.041.441.599.561.501.627.65
8.331.641.798.561.671.838.541.751.866.21
11.72.132.091.772.232.134.772.252.192.92
15.02.372.380.242.492.432.462.522.462.46
T h ,   i n
( ° C )
Q h × 10 6
(m3 s−1)
(b) Cross-Diagonal Carbon-Fiber Spacers (2 mm)
60 ° 90 ° 120 °
N exp × 10 3
kg m−2 s−1
N t h e o × 10 3
kg m−2 s−1
E (%) N exp × 10 3
kg m−2 s−1
N t h e o × 10 3
kg m−2 s−1
E (%) N exp × 10 3
kg m−2 s−1
N t h e o × 10 3
kg m−2 s−1
E (%)
456.670.740.818.740.760.838.750.790.856.79
8.330.840.929.040.870.947.380.910.965.55
11.71.041.062.551.071.070.331.151.104.49
15.01.141.183.261.231.192.661.261.223.66
506.670.971.068.509.781.078.631.041.105.56
8.331.141.269.261.271.270.121.321.311.31
11.71.441.471.691.521.464.091.581.523.82
15.01.681.623.491.741.645.871.801.668.00
556.671.291.418.861.341.457.771.361.509.37
8.331.481.628.731.551.697.861.681.722.60
11.71.801.905.601.851.944.821.951.991.73
15.02.122.072.422.142.102.082.252.163.81
606.671.701.795.271.711.836.891.781.884.88
8.331.942.076.681.982.116.442.072.154.03
11.72.512.356.732.542.405.692.512.462.29
15.02.712.661.852.752.701.972.792.732.07
Table 2. (a) Effects of hydrodynamic angles on flux enhancement for concurrent flow. (b) Effects of hydrodynamic angles on flux enhancement for countercurrent flow.
Table 2. (a) Effects of hydrodynamic angles on flux enhancement for concurrent flow. (b) Effects of hydrodynamic angles on flux enhancement for countercurrent flow.
(a) Effects of Hydrodynamic Angles on Flux Enhancement for Concurrent Flow
T h ,   i n
( ° C )
Q h × 10 6
(m3 s−1)
Empty Channel2 mm3 mm
60 ° 90 ° 120 ° 120 °
N t h e o × 10 3
kg m−2s−1
N t h e o × 10 3
kg m−2s−1
I E N t h e o × 10 3
kg m−2s−1
I E N t h e o × 10 3
kg m−2s−1
I E N t h e o × 10 3
kg m−2s−1
I E
506.670.810.9618.70.9820.40.9922.60.9719.8
8.330.911.1021.61.1223.71.1425.61.1021.7
11.71.031.3026.31.3228.01.3530.71.2824.6
15.01.111.4631.41.4833.11.5034.71.4329.0
606.671.271.5521.71.5925.11.6227.81.5421.3
8.331.431.7925.11.8327.81.8630.31.7824.2
11.71.602.0930.12.1333.22.1936.92.0729.6
15.01.752.3835.62.4338.62.4640.02.3534.1
(b) Effects of Hydrodynamic Angles on Flux Enhancement for Countercurrent Flow
T h ,   i n
( ° C )
Q h × 10 6
(m3 s−1)
Empty channel2 mm3 mm
60 ° 90 ° 120 ° 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 (%)
506.670.881.0621.41.0722.21.1025.91.0621.0
8.331.011.2624.61.2726.11.3129.21.2523.7
11.71.131.4729.71.4629.41.5234.71.4629.5
15.01.201.6235.01.6436.61.6638.71.6033.2
606.671.421.7825.01.8329.01.8832.11.7825.5
8.331.592.0528.92.1031.92.1535.52.0528.7
11.71.742.3535.32.4038.92.4641.12.3233.1
15.01.882.6641.72.7043.62.7345.12.5837.3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/membranes11120973

AMA Style

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(12):973. https://doi.org/10.3390/membranes11120973

Chicago/Turabian Style

Ho, Chii-Dong, Luke Chen, Jun-Wei Lim, Po-Hung Lin, and Pin-Tsen Lu. 2021. "Distillate Flux Enhancement of Direct Contact Membrane Distillation Modules with Inserting Cross-Diagonal Carbon-Fiber Spacers" Membranes 11, no. 12: 973. https://doi.org/10.3390/membranes11120973

APA Style

Ho, C. -D., Chen, L., Lim, J. -W., Lin, P. -H., & Lu, P. -T. (2021). Distillate Flux Enhancement of Direct Contact Membrane Distillation Modules with Inserting Cross-Diagonal Carbon-Fiber Spacers. Membranes, 11(12), 973. https://doi.org/10.3390/membranes11120973

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