Next Article in Journal
A Short Review of Advances in MOF Glass Membranes for Gas Adsorption and Separation
Previous Article in Journal
Enrichment of Nutmeg Essential Oil from Oil-in-Water Emulsions with PAN-Based Membranes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enhancing the Permeate Flux Improvement of Direct Contact Membrane Distillation Modules with Inserted S-Ribs Carbon-Fiber Filaments

1
Department of Chemical and Materials Engineering, Tamkang University, Tamsui, New Taipei 251301, Taiwan
2
Department of Chemical Engineering, Faculty of Engineering, Universiti Teknologi Petronas, Seri Iskandar 32610, Perak, Malaysia
3
CO2 Research Center (CO2RES), Institute of Contaminant Management, Universiti Teknologi Petronas, Seri Iskandar 32610, Perak, Malaysia
*
Author to whom correspondence should be addressed.
Membranes 2024, 14(5), 98; https://doi.org/10.3390/membranes14050098
Submission received: 15 March 2024 / Revised: 19 April 2024 / Accepted: 23 April 2024 / Published: 25 April 2024

Abstract

:
Three widths of manufacturing S-ribs carbon-fiber filaments acting as turbulence promoters were implemented into the flow channel of direct contact membrane distillation (DCMD) modules to augment the permeate flux improvement in the present study. Attempts to reduce the disadvantageous temperature polarization effect were made by inserting S-ribs turbulence promoters in improving pure water productivity, in which both heat- and mass-transfer boundary layers were diminished due to creating vortices in the flow pattern and increasing turbulence intensity. The temperature polarization coefficient t t e m p was studied and found to enhance device performance (less thermal resistance) under inserting various S-ribs carbon-fiber thicknesses and operating both cocurrent- and countercurrent-flow patterns. The permeate fluxes in the DCMD modules with inserted S-ribs carbon-fiber turbulence promoters were investigated theoretically by developing the mathematical modeling equations and were conducted experimentally with various design and operating parameters. The theoretical predictions and experimental results exhibited a great potential to considerably achieve permeate flux enhancement in the new design of the DCMD system. The DCMD module with inserted S-ribs carbon-fiber turbulence promoters in the flow channel could provide a relative permeate flux enhancement up to 37.77% under countercurrent-flow operations in comparisons with the module of using the empty channel. An economic consideration on both permeate flux enhancement and power consumption increment for the module with inserted S-ribs carbon-fiber filaments was also delineated.

1. Introduction

Membrane distillation (MD) has gained significant attraction in both industrial and academic settings in recent years due to its technical feasibility. The advancements in separation and water purification techniques have focused on high permeate flux and low energy consumption, particularly in the areas of membrane materials, module fabrication, transport phenomena, and fouling mitigation. MD [1,2] has gained attractive attention due to the existing thermal gradients built up by both microporous hydrophobic membrane surfaces, which are always in contact with both hot and cold bulk solutions at their saturation points to allow the passage of vapors only [3,4]. The MD process involving both heat and mass transfer undergoes temperature and concentration polarization effects [5] near the membrane surface, which diminish the device performance. A higher temperature causes a larger impact on the temperature polarization effect [6], while a higher feed composition induces a more pronounced influence on the concentration polarization effect [7]. The significant thermal resistance across the membrane with a considerable removal of the latent heat associated to water evaporation occurs on the thermal boundary layer due to the vapor pressure gradient across the hydrophobic membrane, which causes the feed temperature decrement, called as the temperature polarization effect [8], and consequently, the net thermal driving force between the bulk temperature and membrane surface temperature to transport the permeate flux thus declines [9].
Augmentation of the temperature driving-force gradient by implementing turbulence promoters to suppress the boundary layer between the bulk and membrane surface results in the increment of transmembrane permeate flux by alleviating the temperature polarization effect [10]. Incorporating suitable flow amendment configurations of turbulence promoters to amplify the hydrodynamic conditions has proven to be beneficial for the simultaneous mass and heat transfer in MD in the boundary layer, such as eddy promoters [11,12] and flow deflectors [13,14,15]. Moreover, the development of advanced composite membranes [16,17] meets the requirements of an ideal membrane with a lower conductive heat loss and a higher mass-transfer rate. Many previous studies focused on the eddies and wakes owing to vortices and secondary flows [18] as the hot stream passes through the spacer strands and proposed the mathematical model to predict heat- and mass-transfer enhancements in the MD desalination system with respect to traditional membrane modules without inserting turbulence promoters [19]. The mathematical heat and mass transport modeling approaches applied in membrane distillation were reviewed [20] and simulated using computational fluid dynamics (CFD) [12,21].
Various approaches, such as the use of spacers [22], filaments [23], and roughened surfaces [24], have been proposed to improve the hydrodynamic conditions by suppressing the thermal boundary layer, thereby enhancing the heat- and mass-transfer rates. The influence affects both heat-transfer rates and heat-transfer coefficients [25,26] which were recognized to identify the relative significance of each heat-transfer mechanism. All of the spacers, filaments, and roughened surfaces play an important role as a turbulence promoter in both heat and mass transport, and thus, heat-transfer improvement is achieved. In MD modules, permeate is transported with simultaneous heat transfer, resulting in complex heat-transfer mechanisms by the temperature difference between the hot and cold feed streams. The resulting turbulence reduces the temperature polarization effect, leading to improved heat and mass-transfer rates [10]. However, the heat- and mass-transfer rate enhancement are at the expense of the pressure drops with power consumption increments due to the augmented turbulent intensity induced by the turbulence promoters. To assess the economic and technical viability, the ratio of transfer rate enhancement to energy consumption increment was evaluated [27]. A recent study has demonstrated the effectiveness of a segment design involving the machining of a carbon-fiber sheet with various geometric shapes to increase turbulent intensity and alleviate the temperature polarization effect in the flow channel of DCMD modules. This disruption of the thermal boundary layer has been shown to significantly increase the convective heat-transfer coefficient by introducing carbon-fiber turbulence promoters into the hot feed stream, thereby changing the flow direction from the curved-shape grid to the grid, as compared to modules with feed channels lacking carbon-fiber filaments. However, the benefits of using turbulence promoters in DCMD are weighed down by their disadvantages associated with the increment of the pressure drop in the feed stream [28].
The present study actually extends the existing study except for inserting S-ribs carved-shape carbon-fiber filaments instead of using straight-line carbon-fiber filaments [29] to accomplish the augmented turbulent intensity due to changing the flow direction by the hydrodynamic angle. This work is based on the principle of the introduction of S-ribs carbon-fiber turbulence promoters (1 mm thickness) with various widths (3 mm, 4 mm, and 5 mm) as feed channel filaments into the DCMD module, which would suppress the boundary layer between the bulk and membrane surface to achieve a smooth transition for the turbulent flow near the membrane surface of fluid flow due to the increased turbulence intensity and decreased thermal boundary layer resistance. Two strata of microscopic and plug-flow descriptions are used to represent the real processes by mathematics, which evolve from transport phenomena principles. These two descriptions are depicted to be related to the complexity of the heat-transfer mechanism in the present DCMD system and conjugated with the dusty-gas model to describe the essential membrane coefficient models, which can estimate the permeate flux across the microporous hydrophobic membrane. In this study, we develop a theoretical model for the DCMD module to estimate the module’s permeate fluxes and analyze its trade-off considerations for energy consumption. The primary objective is to fabricate S-ribs turbulence promoters and insert them into the hot stream channel to mitigate its temperature polarization effect, ultimately leading to an increase in vapor permeate flux. Additionally, we refine and validate a regression equation for the heat-transfer enhancement factor through experimental verification.

2. Experimental Setup and Materials

The carbon fiber used in this study was chosen concerning its lower cost, mechanical strength to preventing membrane vibration, and easy fabrication of the various geometric shapes of the filaments. In addition, incorporating the mesh design (say S-ribs carbon fiber) into the DCMD module was proposed so as to act as turbulence promoters and amplify the hydrodynamic conditions. The fabrication details of S-ribs carbon-fiber filaments and the schematic configuration are depicted in Figure 1. Figure 2 showcases a photograph of the flat-plate DCMD experimental setup, featuring acrylic plates as external walls within a parallel-plate channel. S-ribs carbon-fiber sheets, cut in an S-shape with a 1 mm thickness, serve as turbulence promoters inserted into the hot saline feed stream to induce eddy motion. The S-ribs carbon-fiber filaments adhered to the membrane surface in the hot feed side were inserted in parallel vertically into a parallel-plate channel to conduct a two-stream operation. Meanwhile, the cold feed side was constructed using a 0.1 mm nylon fiber wound as a supportive grid to prevent membrane vibration and wrinkling. Both supporting materials provide mechanical strength to prevent membrane vibration and act as turbulence promoters.
The experimental setup involves two parallel-plate parallel channels (L = 0.21 m, W = 0.29 m, d = 2 mm), separated by a hydrophobic composite membrane made of PTFE/PP (polytetrafluoroethylene and polypropylene, J020A330R, ADVANTEC Toyo Roshi Kaisha, Ltd., Tokyo, Japan) with a thermal conductivity of 5.0 × 10 4   c a l / c m   s e c   ° C (All-Fluoro Co., Ltd. Taoyuan, Taiwan) and a water vapor permeability resistance of 4   s e c / m [30]. The membrane possesses a nominal pore size of 0.2 µm, a porosity of 0.72, and a thickness of 130 µm, serving as the permeating porous medium in this study. A 2 mm-thick silicon rubber was affixed to the acrylic plate to prevent leakage and create two spacer conduits of 2 mm for each channel, respectively. Figure 2 illustrates the top views of S-ribs carbon-fiber filaments with various widths (3 mm, 4 mm, and 5 mm) as a design parameter. These S-ribs carbon-fiber filaments provide mechanical strength to prevent membrane vibration and act as turbulence promoters. The effective permeate areas were partially obstructed by these three widths of S-ribs carbon-fiber filaments, covering approximately 13% of the hydrophobic membrane, a factor considered in the calculation procedure.
The artificial saline water, consisting of 3.5 wt% NaCl, was prepared by adding inorganic salts (NaCl) to both distilled water and pure water, respectively. This saline solution was pumped (51K40RA-A, ASTK, New Taipei, Taiwan) from the thermostat (G-50 and D650, DENG YNG, New Taipei, Taiwan) at specified temperatures (45 °C, 50 °C, 55 °C, and 60 °C). Inlet and outlet temperatures were measured using thermometer probes (TM-946, Lutron, New Taipei, Taiwan) connected to both sides of the flat-plate membrane modules. The operational conditions of saline feed streams, with various flow rates (0.3, 0.5, 0.7, and 0.9 L/min), were adjusted using flow meters (FE-091312-D, Fong-Jei, Hsinchu, Taiwan) and a controller (N12031501PC-540, Protec, Brooks Instrument, Hatfield, PA, USA), while maintaining the temperature at 25 °C for the cold stream (FN-0423112-F, Fong-Jei, Hsinchu, Taiwan). The conductance of the permeate flux was collected and measured, being less than 1.5 μs/cm. Comparisons were made of permeate fluxes under various operation conditions to assess the device performance between two modules with and without inserting S-ribs carbon-fiber turbulence promoters. The experimental run of the permeate flux was collected and weighed using an electronic balance (XS 4250C, Precisa Gravimetrics AG, Dietikon, Switzerland) for measurement and recording on the PC.

3. Theoretical Formulations

3.1. Mass and Heat Transfer

Theoretical statements combining both heat- and mass-transfer mechanisms were formulated for the DCMD system with microscopic and plug-flow descriptions, as shown in Figure 3a,b, respectively. Figure 3a,b are two models that evolve from transport phenomena principles. Two levels of microscopic and plug-flow descriptions used to represent real processes by mathematics. These two levels are depicted to be related to the complexity of the heat-transfer mechanism in the present DCMD system. The microscopic description involves a phenomenological approach, representing the system as a continuum, while the plug-flow description deals with only the largest component of the gradient in the balance equation by the bulk flow, neglecting all diffusion terms. A mass-transfer model, coupled with heat-transfer behavior, was developed to illustrate the concentration gradient of the hot saline feedwater, leading to the vapor diffusing exclusively through the porous hydrophobic membrane and condensing in the cold stream, completing the distillation process. This occurs due to the temperature gradients at both membrane surfaces. The theoretical analysis of the DCMD aims to elucidate how the saline water is initially vaporized and eventually condensed at the pore entrances of both membrane surfaces. This is achieved through the permeate/water equilibrium governed by vapor diffusion and enthalpy flow conservation through heat conduction, as detailed below.
The membrane permeation coefficient ( c m ) and the trans-membrane saturation vapor pressure difference ( Δ P ) were evaluated to determine the permeate flux for membrane distillation processes [31,32].
N = c m Δ P = c m P 1 s a t ( T 1 ) P 2 s a t ( T 2 ) = c m d P d T T m T 1 T 2 = c m P m λ M w R T m 2 T 1 T 2
The dusty-gas model was used to describe three essential membrane coefficient models based on a comparison between the mean free path of vapor molecules and the membrane pore sizes, which can estimate the mass flux across the microporous hydrophobic membrane. The influence of the Poiseuille flow can be neglected [33] when the pore size of the membrane is relatively small. Considering that a membrane of 0.2 µm is used in the present study, we may conclude that the addition of the Knudsen and molecular diffusion models is appropriate to estimate the permeate flux through the membrane. The combination of the Knudsen diffusion model (due to the smaller mean free path of vapor molecules than the membrane pore size) and the molecular diffusion model (due to the concentration gradient across the membrane) was investigated [34] to dominate the mass-transfer mechanism and presented by an equation of the membrane permeation coefficient for deaerated microporous membranes to predict the permeate flux through the membrane. The membrane permeation coefficient is the combination of the Knudsen diffusion model and molecular diffusion model as follows:
c m = 1 c K + 1 c M 1 = 2 8 π 3 ε   r τ δ m M w R T m 1 / 2 1 + Y m l n D m ε δ m τ M w R T m 1 1
The temperature differences between the streams near the membrane surfaces and those of the bulk feed streams are utilized to estimate the temperature polarization coefficient. Predictions are made regarding the temperature gradients of membrane surfaces on both the feed stream and permeate side across the entire module, influencing heat transfer. Consequently, the permeate flux is calculated using the mass-transfer modeling equation presented in Equation (1). The saturation vapor pressure ( P 1 s a t ) of the hot saline feed side is estimated using the water activity coefficient ( a w ), determined through a correlation [1]:
P 1 s a t = x w a w P w s a t
a w = 1 0.5 x NaCl 10 x NaCl 2
and the tortuosity ( τ ) can be estimated using the porosity of the membrane [35]
τ = 1 / ε
According to the macroscopic description of the temperature gradient of the DCMD module in Figure 3a, the heat balances of enthalpy flow conservation in the non-isothermal process were made within each heat-transfer region in the DCMD module as follows: (a) the hot saline water stream; (b) the hydrophobic composite membrane; and (c) the cooling water stream. Vapor flux permeating through the hot feed stream, microporous hydrophobic membrane and cold stream for the modules with/without inserted S-ribs carbon-fiber turbulence promoters was expressed in terms of heat-transfer resistances in series due to the temperature gradient, as shown in Figure 4.
The energy balance equations may be derived by balancing energy under the steady-state operation in both feed channels to and from the membrane surfaces. Equating the amount of heat flux by the conservation law among three regions, one may obtain the following:
q h = h h ( T h T 1 )   the   hot   saline   water   feed   region
q m = N λ + k m δ m     ( T 1 T 2 )   the   membrane   region
q c = h c ( T 2 T c )   the   cooling   water   region
where N λ is categorized as the latent heat of vaporization and k m / δ m = h h ( T 1 T 2 ) is the heat conduction, and 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 [36] as
k m = ε   k g + 1 ε k s
Equation (7) can be rewritten by using Equation (1) in terms of the overall heat-transfer coefficient of the membrane Hm including the latent heat across the membrane, which is implicitly related to the temperature polarization coefficient as follows:
q m = N λ + k m δ m     ( T 1 T 2 ) = c m 1 x NaCl 1 0.5 x NaCl 10 x NaCl 2 P 2 + P 1 λ 2 M w 2 R T m 2 + k m δ m T 1 T 2 = H m ( T 1 T 2 )

3.2. Temperature Polarization Coefficient

The temperature polarization coefficient τ t e m p is an indicator to reveal the extent of the thermal boundary-layer resistance between both the hot saline and cooling feed streams, as indicated in Figure 5a,b, which controls the permeate flux through the membrane, and τ t e m p is commonly defined as follows:
τ t e m p = ( T 1 T 2 ) / ( T h T c )
Actually, inserting the S-ribs carbon-fiber filament in the hot saline water compartment would increase h h (for the module without promoter insertion) to h h p (for the module with promoter insertion) but h c   does not change. According to Equations (6)~(8) and (10), all heat-transfer regions of the microscopic description under steady-state operations with assuming q = q h = q m = q c were developed and illustrated by the schematic diagram in Figure 3a.
q = h h T h T 1   = H m   ( T 1 T 2 ) = h c ( T 2 T c )
or
q = T h T c 1 h h + 1 H m + 1 h c = h c T 2 T c   = q c
Equation (13) can be rewritten as
T 2 T c = T h T c h c h h + h c H m + 1
Similarly, we have the increased value of the temperature gradient (say T 2 p T c ) by using the same performed procedure of Equation (14) in the cold feed region for h h p > h h due to inserting the turbulence promoter, as shown in Figure 5b.
T h T c h c h h p + h c H m + 1 = T 2 p T c > T 2 T c
Meanwhile, the temperature polarization coefficient of Equation (11) increases with the promoter insertion, i.e.,
T h T c H m h h p + 1 + H m h c = T 1 p T 2 p = T p > T = T 1 T 2 = T h T c H m h h + 1 + H m h c
The temperature polarization effect was reduced, as demonstrated in Figure 5b, due to the DCMD module with inserted S-ribs carbon-fiber turbulence promoters, which resulted in the temperature difference enlargement of both membrane surface temperatures T p ( T 1 p T 2 p ) > T ( T 1 T 2 ). The gradient of T increases to a higher driving force T p with inserted S-ribs carbon-fiber turbulence promoters, and thus, the vapor flux is increased, as shown in Figure 5b. The S-ribs carbon-fiber filaments have been validated to disrupt the laminar boundary layer and intensify vortices or secondary flow characteristics on the membrane surface, responding to changing hydrodynamic conditions. This disruption plays a crucial role in mitigating the temperature polarization effect and enhancing a higher permeate flux.
The alternative expression of Equation (11) was obtained by equating the energy conservation of heat fluxes in each region, say Equations (6) and (10) ( q h = q m ) and Equations (8) and (10) ( q m = q c ), respectively, to replace both membrane surface temperatures ( T 1 and T 2 ) in terms of the heat-transfer coefficients ( h h , h c and H m ) as follows:
T h = T 1 + H m h h ( T 1 T 2 )
T c = T 2 H m h c T 1 T 2
Then, a more simplified form of τ t e m p in Equation (11) can be rewritten in terms of the heat-transfer coefficient as
τ t e m p = h h h c h h h c + h h H m + h c H m
The procedure for calculating theoretical predictions of the mass-transfer coefficient was performed by continuously iterating T 1 and T 2 from Equations (17) and (18) within the convergence tolerance. The calculated convective heat-transfer coefficients were calculated from Equation (1) and validated by the experimental results, delivered to predict theoretically not only in the hot/cold bulk flows ( T h and T c ) but also those on the membrane surfaces ( T 1 and T 2 ) of both hot and cold feed streams, respectively.

3.3. Governing Equations by Macroscopic Modeling

Both temperature variations in both hot saline and cooling feed streams along the flow direction were developed, as illustrated in Figure 6, by energy conservation in one-dimensional governing equations according to the macroscopic modeling to solve the temperature distributions of both streams in terms of the temperature polarization coefficient τ t e m p as
d T h d z = q W Q h   ρ h C p , h = W Q h   ρ h C p , h H m τ t e m p ( 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 e m p ( T h T c ) Cocurrent-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 e m p ( T h T c ) Countercurrent-flow   operations
Temperature fields of both hot and cold feed streams were solved with the use of two simultaneous ordinary differential equations of Equations (20) and (21) for cocurrent-flow operations (or Equations (20) and (22) for countercurrent-flow operations) by marching the 4th Runge-Kutta method along the length of the module, which can quantify the enhanced permeate flux by the turbulence intensity augmentation with inserted S-ribs carbon-fiber turbulence promoters.

3.4. Hydraulic Consumption Increment

The hydraulic consumption increment is expected and required due to inserting S-ribs carbon-fiber filaments into the hot saline stream, which may be determined using the Fanning friction factor ( f F ) [37], considering only the friction losses to walls of both hot and cold streams as
H = H h + H c = Q h ρ h l w f , h + Q c ρ c l w f , c
l w f , i = 2 f F , i v ¯ i 2 L D h , i , i = h , c
in which
v ¯ h = Q h d W D 1 W 1 N 1 ,   v ¯ c = Q c A  
D e h , h = 4 d W D 1 W 1 N 1 2 d + W + D 1 N 1   , D e h , c = 4 d W 2 d + W
The hydraulic equivalent diameter D h , h of modules with embedding S-ribs carbon-fiber turbulence promoters was calculated by the wetted area A and wetted perimeter P as the average carbon-fiber width of S-ribs carbon-fiber turbulence promoters, say 4 A / P , as shown in Figure 7, which were evaluated by averaging various sections of the S-rib shape.
The Fanning friction factor can be estimated using a correlation based on the aspect ratio of the channel ( α = d / W ) [38]:
f F , h = C R e h , f F , c = C R e c
C = 24 1 1.3553 α + 1.9467 α 2 1.7012 α 3 + 0.9564 α 4 0.2537 α 5
The relative extent of the hydraulic consumption increment I P was illustrated by calculating the percentage increment in the module with inserted S-ribs carbon-fiber filaments, which was based on the module of using empty channels as
I p = H p r o m o t e r H e m p t y H e m p t y × 100 %

3.5. Heat-Transfer Enhancement Factor

S-ribs carbon-fiber turbulence promoters were incorporated into the open conduit of the hot saline feed stream, replacing the device of an empty channel. The correlation of heat-transfer coefficients [39] for the module using filament-filled channels was proposed through multiple linear regressions. The heat-transfer enhancement factor α p was introduced [11] to calculate augmented convective heat-transfer coefficients in DCMD modules with inserted S-ribs carbon-fiber filaments, employing an iterative procedure:
α p = N u p / N u l a m
where the heat-transfer equivalent diameter D e h , h h   was defined [40] as follows:
D e h , h h = 4 d W D 1 W 1 N 1 W D 1 N 1
N u l a m = 4.36 + 0.036 RePr D e h , h h / L 1 + 0.011 RePr D e h , h h / L 0.8 Empty   channel   module
N u p = h h D e h , h h k S-ribs   carbon-fiber-filled   module
The S-ribs carbon-fiber filaments, acting as eddy promoters within the flow channel, play a vital role and offer a comprehensive interpretation of both heat- and mass-transfer behaviors. This interpretation is based on dimensional analysis using Buckingham’s π theorem, relating to the Nusselt number as:
N u p = f W 1 D e h , h h , R e , P r
where W 1 and D e h , h h are the carbon-fiber width and heat-transfer equivalent diameter of channels with S-ribs carbon-fiber filaments, respectively.

4. Results and Discussions

4.1. Lessening Temperature Polarization Effect by Inserting Carbon-Fiber Filaments

Both bulk temperature distributions of hot saline and cold feed streams as well as membrane surface temperatures in DCMD modules were solved numerically using the one-dimensional theoretical model, as presented in Figure 8, along the axial coordinate with the carbon-fiber width as a parameter under both cocurrent- and countercurrent-flow operations, respectively.
In cocurrent-flow operations, theoretical predictions indicate a tapering of the membrane-surface and bulk temperatures along the flowing direction, leading to a reduction in driving-force temperature gradients. Conversely, countercurrent-flow operations maintain a relatively higher average value of the driving-force temperature gradients, resulting in a greater enhancement of permeate flux. The analysis reveals that temperature gradients across both membrane surfaces are higher in flow channels with inserted S-ribs carbon-fiber turbulence promoters compared to those in the empty channel. These increased temperature gradients facilitate greater vapor transport through the membrane, consequently leading to a higher amount of permeate flux condensed in the cooling stream.
Furthermore, theoretical predictions of temperature polarization coefficients τ t e m p were determined according to Equation (19) and are presented in Figure 9. The analysis considered the inlet saline temperature and carbon-fiber widths as parameters. The τ t e m p increased with decreasing both the carbon-fiber width and the inlet saline temperature for both cocurrent- and countercurrent-flow operations. The descending permeate flux along the flowing channel for cocurrent-flow operations is thus expected in contrast to the almost constant temperature gradient of the countercurrent-flow operations. The larger inlet saline temperature associated with the higher vapor pressure creates the greater permeate flux passing through the membrane and results in the lower temperature difference (say T 1 p T 2 p ) across the membrane surfaces; therefore, a lower value of τ t e m p is achieved. This is because the greater permeate flux transferring through the membrane results in reducing the saturation pressure as well as the water activity coefficient of the hot saline stream. However, comparisons of τ t e m p were made on operating the module with inserted S-ribs carbon-fiber filaments of 3 mm and 5 mm with that of the empty channel, as shown in Figure 9. A larger τ t e m p was accomplished substantially when the module included inserted S-ribs carbon-fiber filaments due to disturbing the thermal boundary layer on the membrane surface with a smaller thermal resistance. Suppressing the temperature polarization effect results from inserting S-ribs carbon-fiber filaments in the hot saline feed stream to diminish the thermal boundary layer thickness on the membrane surface. Meanwhile, the effect reduction in τ t e m p is more significant in countercurrent-flow operations, and thus, a higher τ t e m p value is achieved in the countercurrent-flow operation, as demonstrated in Figure 9.

4.2. Permeate Flux Enhancement by Inserting S-Ribs Carbon-Fiber Turbulence Promoters

Applying the regressed correlation approach and quantifying the enhanced permeate flux due to inserting S-ribs carbon-fiber turbulence promoters were solved directly through the permeate flux enhancement factor, which was established based on dimensional analysis of using Buckingham’s π theorem in Equation (34). The experimental runs with empty channels and 3 mm, 4 mm, and 5 mm S-ribs carbon-fiber widths were used to determine the correlation for the permeate enhancement factor α p , which was expressed in Equation (30) of the Nusselt number as a measure of heat-transfer efficiency as follows:
α p = N u p N u l a m = 1.72 ( W 1 D e h , h h ) 0.165 R e 0.04 P r 0.321
The Nusselt number for the empty channel is presented in Equation (32), while the corrected Nusselt number for S-ribs carbon-fiber-filled channels is presented in Equation (33). Furthermore, the correlated Sherwood number was incorporated into the mass-transfer enhancement factor α p , as the improved heat-transfer coefficient reduces temperature polarization, and thus the driving force across the membrane is increased. The heat-transfer enhancement factor was derived from the correlation via a regression analysis and expressed in Equation (35) for implementing S-ribs carbon-fiber turbulence promoters, which results in the augmented convective heat-transfer coefficients in membrane distillation modules and presented in Equation (35) as well as in Figure 10. A regression analysis was set up the normal equations for the least square parameters to obtain the correlated equation and the squared correlation coefficient ( R 2 ) is 0.94. It is concluded from Figure 10 that these corrected Nusselt numbers are more significant in the modules with inserted S-ribs carbon-fiber filaments.

4.3. Accuracy Deviations between Experimental and Theoretical Results

The precision analysis of experimental uncertainty for each individual measurement from the experimental results was determined [41] as follows:
S N e x p = i = 1 N e x p N exp , i N exp , i ¯ 2 N e x p 1 1 / 2
where the mean value of the resulting uncertainty of the experimental measurements was defined by
S N e x p ¯ = S N e x p N e x p
The mean experimental uncertainty in Figure 11, Figure 12 and Figure 13 ranges within 4.81 × 10 3 S N e x p ¯ 9.73 × 10 3 . Meanwhile, the accuracy deviation between the experimental results and theoretical predictions was calculated as follows:
E = 1 N e x p i = 1 N e x p N t h e o , i N exp , i N t h e o , i
where N e x p , N exp , i , and N t h e o , i are the number of experimental data, theoretical predictions, and experimental results of permeate flux, respectively. The accuracy deviation of the experimental results from the theoretical predictions is well minimized within 1.83 × 1 0 2 E 9.97 × 1 0 2 .
Additionally, the predictive capability for permeate flux can be extended to various geometric promoter designs by following the same regression procedure applied to the Nusselt numbers for both modules: one with an empty channel and the others with inserted S-ribs carbon-fiber filaments. In essence, the insertion of S-ribs carbon-fiber turbulence promoters disrupts the thermal boundary layer on the membrane surface, thereby reducing heat-transfer resistance and enhancing permeate flux. The results indicate that higher inlet saline feed temperatures lead to larger N u p numbers, resulting in a higher heat-transfer rate. Both experimental findings and theoretical predictions of permeate flux were visually presented, utilizing inlet saline feed temperatures and feed flow rates as parameters. This representation is outlined in Figure 11, Figure 12 and Figure 13 for both cocurrent- and countercurrent-flow operations, considering 3 mm, 4 mm, and 5 mm carbon-fiber widths, respectively. The fair consistency and agreement between theoretical predictions and experimental data are evident, providing a robust basis for evidence-based validation, as illustrated in Figure 11, Figure 12 and Figure 13. The order of permeate flux magnitude is observed as 3 mm > 4 mm > 5 mm.
It is noteworthy that a larger permeate flux was achieved at a higher hot saline feed temperature of 60 °C compared to 50 °C. This can be attributed to a greater permeate flux emerging due to a higher saturated vapor pressure gradient between both sides of the membrane. Additionally, a higher hot saline feed flow rate results in an increased permeate flux, facilitated by velocities and vortices that effectively mitigate heat-transfer resistance, especially with smaller carbon-fiber widths impacting the thermal boundary layer. Moreover, the module incorporating 3 mm S-ribs carbon-fiber turbulence promoters in the flowing channel generates more intensive vortices and eddies than those using wider S-ribs carbon-fiber filaments under the same total coverage area. As anticipated, countercurrent-flow operations yield a more significant enhancement in permeate flux compared to cocurrent-flow operations when S-ribs carbon-fiber turbulence promoters are inserted.
A relative increase in permeate flux I N was illustrated by calculating the percentage increment in comparisons between the permeate flux of the module using the empty channel and inserted S-ribs carbon-fiber turbulence promoters as
I N = N p r o m o t e r N e m p t y N e m p t y × 100 %
A relative increment of permeate flux I N in the module with inserted S-ribs carbon-fiber turbulence promoters was calculated in comparisons with the module using the empty channel, which were evaluated in terms of the theoretical predictions of permeate flux N , as summarized in Table 1 with carbon-fiber width, inlet saline feed temperature and feed flow rate as parameters.
The module with inserted S-ribs carbon-fiber filaments of 3 mm in width exhibits a relative increment in permeate flux of up to 37.77% under countercurrent-flow operations compared to the module using an empty channel. Furthermore, the analysis from Table 1 reveals that the order of device performance for permeate flux enhancement with S-ribs carbon-fiber turbulence promoters increases with higher inlet saline feed temperatures and feed flow rates but decreases with carbon-fiber width. Overall, the insertion of S-ribs carbon-fiber turbulence promoters into the flow channel demonstrates significant potential for substantially augmenting permeate flux in the DCMD module by mitigating the temperature polarization effect.
The present work extends the existing study except for inserting S-ribs carbon-fiber turbulence promoters instead of using straight-line carbon-fiber filaments [29]. The impact of the more dominating operational parameter of inlet temperatures was monitored in the DCMD system and the experimental results of the permeate fluxes were collected and weighed to evaluate the improved device performance. Therefore, both inlet and outlet temperatures were measured using thermometer probes for each 5 min interval until the outlet temperatures reached an unchanged steady state, and thus comparisons were made of permeate fluxes under various operation conditions for both modules with and without inserted S-ribs carbon-fiber turbulence promoters. The present study shows a graphical representation of comparisons with the theoretical predictions of permeate fluxes obtained in the present study and straight-line carbon-fiber filaments [29], which illustrates why there is a preference for the present design of inserting S-ribs carbon-fiber turbulence promoters, as shown in Figure 14 for both cocurrent- and countercurrent-flow operations. This is the value and originality of the present study with consideration of the economic viewpoint and technical feasibility.

4.4. Energy Consumption Increment

The enhancement of permeate flux is counterbalanced by an increase in energy consumption, constituting trade-offs attributed to heightened turbulent intensity and additional friction losses resulting from the insertion of S-ribs carbon-fiber filaments into the flow channel. An economic consideration on both permeate flux enhancement and power consumption increment for the module with inserted S-ribs carbon-fiber filaments in the present study was also delineated for economic and technical feasibilities, as shown in Figure 15a,b. The impact of hot saline flow rate, inlet saline feed temperature, and carbon-fiber width on the ratio I N / I P is illustrated in Figure 15a,b. This ratio represents the relationship between the power consumption increment and permeate flux enhancement in both cocurrent- and countercurrent-flow operations. An assessment of the effectiveness of membrane turbulence promoters based on economic viewpoint and technical feasibilities was conducted to determine suitable operation and design parameters that balance desirable permeate flux enhancement with undesirable energy consumption increments.
A larger I N / I P value suggests that the increment in energy consumption can compensate for the decrement in permeate flux, owing to the enhancement in convective heat-transfer coefficients increasing permeate flux with the insertion of S-ribs carbon-fiber turbulence promoters. The I N / I P ratio for the channel with wider carbon-fiber filaments exceeds that of the channel with narrower filaments, indicating that adjusting the carbon-fiber width appropriately can achieve more effective permeate flux at the expense of energy consumption. The I N / I P value increases with the hot saline feed temperature and hot saline flow rate but decreases with the S-ribs carbon-fiber width, as revealed in Figure 15a,b. In essence, embedding S-ribs carbon-fiber turbulence promoters in the hot saline feed channel achieves a superior permeate flux enhancement at the cost of a higher friction loss increment under higher inlet saline temperatures. Notably, the values of I N / I P increase with the hot feed flow rate, and larger temperature driving-force gradients in operating countercurrent-flow systems result in a larger value of I N / I P .

5. Conclusions

The theoretical predictions were calculated and validated by experimental results under various hot feed flow rates, inlet saline temperatures, and various carbon-fiber widths for both cocurrent- and countercurrent-flow operations. Thorough comparisons of permeate flux improvement lead to the following conclusions:
  • Inserting S-ribs carbon-fiber filaments of 3 mm in width into the saline feed flow channel results in relative increases in permeate flux up to a maximum permeate flux improvement of 37.77% under countercurrent-flow operations compared to the module using an empty channel.
  • The results show that permeate flux improvement decreases with the width of carbon-fiber filaments, but the ratio of permeate flux improvement to power consumption increment (say I N / I P ) increases with the width of carbon-fiber filaments.
  • Permeate flux improvement is more pronounced in countercurrent-flow operations compared to cocurrent-flow operations due to the attainment of a larger temperature gradient.
The correlated equation of the Nusselt number, derived from the theoretical model, proves valuable for designing a more efficient DCMD for membrane desalination applications. While this paper specifically focuses on assessing permeate flux improvement and energy consumption increment by inserting S-ribs carbon-fiber filaments as turbulence promoters into the saline feed channel, further investigation is needed to explore alternative geometric shapes and array configurations of S-ribs carbon-fiber filaments for optimal operation, taking economic feasibility into account.

Author Contributions

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

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Acknowledgments

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

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

A Cross-sectional area of flow channel (m2)
a w Water activity in NaCl solution
C Friction losses coefficient
C P Heat capacity ( J   k g 1 K 1 )
c k Membrane coefficient based on the Knudsen diffusion model ( mol   m 2 P a 1 s 1 )
c M Membrane coefficient based on the molecular diffusion model ( mol   m 2 P a 1 s 1 )
c m Membrane permeation coefficient ( mol   m 2 P a 1 s 1 )
d Channel height (m)
D 1 Carbon-fiber thickness (m)
D m Diffusion coefficient of air and vapor in the membrane ( m 2 s 1 )
D e h , i Hydraulic equivalent hdiameter of channel (m),   i = h ,   c
D e h , h h The heat-transfer equivalent diameter (m)
E Accuracy deviation of experimental results from the theoretical predictions
f F , i Fanning friction factor,   i = h ,   c
h c Convective heat-transfer coefficient of cold feed (W m 2 K 1 )
h h Convective heat-transfer coefficient of hot saline feed (W m 2 K 1 )
h h p Convective heat-transfer coefficient of hot saline feed with promoter insertion (W m 2 K 1 )
H i Hydraulic dissipate energy (W),   i = h ,   c
I N Permeate flux relative factor
I p Power consumption relative index
k Thermal conductivity of water ( J s 1 m 1 K 1 )
k g Thermal conductivity of gas ( J s 1 m 1 K 1 )
k m Thermal conductivity of membrane ( J s 1 m 1 K 1 )
k s Thermal conductivity of solid membrane ( J s 1 m 1 K 1 )
L Channel length (m)
l w f , j Friction loss (J kg−1), i = h , c
M W Molecular weight of water (kg mol−1)
N 1 Number of carbon-fiber filaments
N Distillate flux ( kg   m 2 s )
N ¯ Average distillate flux ( kg   m 2 s )
N e x p Number of experimental measurements
N u p Enhanced dimensionless Nusselt number
N u l a m Nusselt number for laminar flow
P m Mean saturated pressure in membrane (Pa)
P s a t Saturation vapor pressure (Pa)
Q Volumetric flow rate (m3 s−1)
q Heat flux ( J   m 2 s 1 )
r p Membrane pore radius (m)
R Gas constant (8.314 J mol−1 K−1)
Re   Reynolds number
S N e x p Precision index of an experimental measurements of permeate flux (kg m−2 s−1)
S N e x p ¯ Mean value of S N e x p (kg m−2 s−1)
T 1 Membrane surface temperature in the hot saline feed region (°C)
T 2 Membrane surface temperature in the cold feed region (°C)
T 1 p Membrane surface temperature with promoter insertion in the hot saline feed region (°C)
T 2 p Membrane surface temperature with promoter insertion in the cold feed region (°C)
T m Mean temperature in membrane (°C)
v ¯ Average velocity ( m   s 1 )
W Width of channel (m)
W 1 Carbon-fiber width (m)
x w Liquid mole fraction of water
x NaCl Mole fraction of NaCl in saline solution
Y w Vapor mole fraction of water
Y m l n Natural log mean Vapor mole fraction of water in the membrane
z Axial coordinate along the flow direction (m)
Greek letters
α p Enhancement factor
δ m Thickness of membrane (µm)
εMembrane porosity
η v Gas viscosity ( Ns   m 2 )
λ Latent heat of water ( J   k g 1 )
μ Viscosity ( Ns   m 2 )
ρ Density ( kg   m 3 )
τ Membrane tortuosity
τ t e m p Temperature polarization coefficients
Subscripts
1Membrane surface on hot fluid side
2 Membrane surface on cold fluid side
cCold feed stream
hHot feed stream
cor.Correlated results
emptyChannel without embedding turbulence promoters
exp.Experimental results
inAt the inlet
l a m Laminar flow
outAt the outlet
promoterChannel with embedding turbulence promoters
theoTheoretical predictions

References

  1. Lawson, K.W.; Lloyd, D.R. Membrane distillation. II. Direct contact MD. J. Membr. Sci. 1996, 120, 123–133. [Google Scholar] [CrossRef]
  2. Khayet, M. Membranes and theoretical modeling of membrane distillation: A review. Adv. Colloid Interface Sci. 2011, 164, 56–88. [Google Scholar] [CrossRef] [PubMed]
  3. Al-Obaidani, S.; Curcio, E.; Macedonio, F.; Profio, G.D.; Al-Hinai, H.; Drioli, E. Potential of membrane distillation in seawater desalination: Thermal efficiency, sensitivity study and cost estimation. J. Membr. Sci. 2008, 323, 85–98. [Google Scholar] [CrossRef]
  4. Huang, Y.-X.; Wang, Z.; Horseman, T.; Livingston, J.L.; Lin, S. Interpreting contact angles of surfactant solutions on microporous hydrophobic membranes. J. Membr. Sci. Lett. 2022, 2, 100015. [Google Scholar] [CrossRef]
  5. Martínez-Díez, L.; Vazquez-Gonzalez, M.I. Temperature and concentration polarizationin membrane distillation of aqueous salt solutions. J. Membr. Sci. 1999, 156, 265–273. [Google Scholar] [CrossRef]
  6. El-Bourawi, M.S.; Ding, Z.; Ma, R.; Khayet, M. A framework for better understanding membrane distillation separation process. J. Membr. Sci. 2006, 285, 4–29. [Google Scholar] [CrossRef]
  7. Yun, Y.; Ma, R.; Zhang, W.; Fane, A.; Li, J. Direct contact membrane distillation mechanism for high concentration NaCl solutions. Desalination 2006, 188, 251–262. [Google Scholar] [CrossRef]
  8. Ali, A.; Macedonio, F.; Drioli, E.; Aljlil, S.; Alharbi, O.A. Experimental and theoretical evaluation of temperature polarization phenomenon in direct contact membrane distillation. Chem. Eng. Res. Des. 2013, 91, 1966–1977. [Google Scholar] [CrossRef]
  9. 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]
  10. Ho, C.-D.; Chang, H.; Chang, C.-L.; Huang, C.-H. Theoretical and experimental studies of flux enhancement with roughened surface in direct contact membrane distillation desalination. J. Membr. Sci. 2013, 433, 160–166. [Google Scholar] [CrossRef]
  11. Phattaranawik, J.; Jiraratananon, R.; Fane, A.G. 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]
  12. 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]
  13. Thomas, N.; Sreedhar, N.; Al-Ketan, O.; Rowshan, R.; Abu Al-Rub, R.K.; Arafat, H. 3D printed triply periodic minimal surfaces as spacers for enhanced heat and mass transfer in membrane distillation. Desalination 2018, 443, 256–271. [Google Scholar] [CrossRef]
  14. Lee, J.-Y.; Tan, W.S.; An, J.; Chua, C.K.; Tang, C.Y.; Fane, A.G.; Chong, T.H. The potential to enhance membrane module design with 3D printing technology. J. Membr. Sci. 2016, 499, 480–490. [Google Scholar] [CrossRef]
  15. Tsai, H.-Y.; Huang, A.; Soesanto, J.F.; Luo, Y.-L.; Hsu, T.-Y.; Chen, C.-H.; Hwang, K.-J.; Ho, C.-D.; Tung, K.-L. 3D printing design of turbulence promoters in a cross-flow microfiltration system for fine particles removal. J. Membr. Sci. 2018, 573, 647–656. [Google Scholar] [CrossRef]
  16. 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]
  17. Chen, L.-H.; Huang, A.; Chen, Y.-R.; Chen, C.-H.; Hsu, C.-C.; Tsai, F.-Y.; Tung, K.-L. Omniphobic membranes for direct contact membrane distillation: Effective deposition of zinc oxide nanoparticles. Desalination 2018, 428, 255–263. [Google Scholar] [CrossRef]
  18. 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]
  19. Santos, J.L.C.; Geraldes, V.; Velizarov, S.; Crespo, J.G. 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]
  20. Schofield, R.W.; Fane, A.G.; Fell, C.J.D. Heat and mass transfer in membrane distillation. J. Membr. Sci. 1987, 33, 299–313. [Google Scholar] [CrossRef]
  21. 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]
  22. 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]
  23. 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 2017, 185, 2045–2057. [Google Scholar] [CrossRef]
  24. 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]
  25. Bird, R.B.; Stewart, W.E.; Lightfoot, E.N. Transport Phenomena; Wiley: New York, NY, USA, 1960. [Google Scholar]
  26. Mills, A.F. Mass Transfer, 2nd ed.; Prentice-Hall: Englewood Cliffs, NJ, USA, 2001. [Google Scholar]
  27. Da Costa, A.; Fane, A.; Wiley, D. Spacer characterization and pressure drop modelling in spacer-filled channels for ultrafiltration. J. Membr. Sci. 1994, 87, 79–98. [Google Scholar] [CrossRef]
  28. Taamneh, Y.; Bataineh, K. Improving the performance of direct contact membrane distillation utilizing spacer-filled channel. Desalination 2017, 408, 25–35. [Google Scholar] [CrossRef]
  29. Ho, C.D.; Chang, H.; Tsai, C.H.; Lin, P.H. Theoretical and experimental studies of a compact multi-unit direct contact membrane distillation module. Ind. Eng. Chem. Res. 2016, 55, 5385–5394. [Google Scholar]
  30. Fukazawa, T.; Kawamura, H.; Tamura, T. Water Vapour Resistance of Hydrophobic Microporous Membranes under Reduced Pressures at a Constant Temperature. J. Text. Inst. 1999, 90 Pt 1, 602–615. [Google Scholar] [CrossRef]
  31. Lawson, K.W.; Lloyd, D.R. Membrane distillation. J. Membr. Sci. 1997, 124, 1–25. [Google Scholar] [CrossRef]
  32. 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]
  33. Ding, Z.; Ma, R.; Fane, A. A new model for mass transfer in direct contact membrane distillation. Desalination 2003, 151, 217–227. [Google Scholar] [CrossRef]
  34. Zhang, J.; Gray, S.; Li, J.-D. Modelling heat and mass transfers in DCMD using compressible membranes. J. Membr. Sci. 2012, 387–388, 7–16. [Google Scholar] [CrossRef]
  35. 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]
  36. Warner, S.B. Fiber Science; Princeton-Hall: Englewood Cliffs, NJ, USA, 1995. [Google Scholar]
  37. Welty, J.R.; Wick, C.E.; Wilson, R.E. Fundamentals of Momentum, Heat, and Mass Transfer, 3rd ed.; John Wiley & Sons: New York, NY, USA, 1984. [Google Scholar]
  38. Kakac, S.; Shah, R.K.; Aung, W. Handbook of Single-Phase Convective Heat Transfer; Wiley: New York, NY, USA, 1987. [Google Scholar]
  39. Phattaranawik, J.; Jiraratananon, R.; Fane, A.G. Effect of pore size distribution and air flux on mass transport in direct contact membrane distillation. J. Membr. Sci. 2003, 215, 75–85. [Google Scholar] [CrossRef]
  40. Kern, D.Q. Process Heat Transfer; McGraw-Hill: New York, NY, USA, 1950. [Google Scholar]
  41. Moffat, R.J. Describing the uncertainties in experimental results. Exp. Therm. Fluid Sci. 1988, 1, 3–17. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of fabrication structure and experimental setup for the DCMD system of the S-ribs carbon-fiber-filled DCMD module.
Figure 1. Schematic diagram of fabrication structure and experimental setup for the DCMD system of the S-ribs carbon-fiber-filled DCMD module.
Membranes 14 00098 g001
Figure 2. A photo of the experimental setup of three widths of S-ribs carbon-fiber filaments.
Figure 2. A photo of the experimental setup of three widths of S-ribs carbon-fiber filaments.
Membranes 14 00098 g002
Figure 3. Microscopic and plug-flow descriptions of mass- and heat-transfer mechanisms. (a) Microscopic description; (b) Plug-flow description. (Membranes 14 00098 i001Hot fluid; Membranes 14 00098 i002Cold fluid; Membranes 14 00098 i003Permeate flux; Membranes 14 00098 i004Heat flux).
Figure 3. Microscopic and plug-flow descriptions of mass- and heat-transfer mechanisms. (a) Microscopic description; (b) Plug-flow description. (Membranes 14 00098 i001Hot fluid; Membranes 14 00098 i002Cold fluid; Membranes 14 00098 i003Permeate flux; Membranes 14 00098 i004Heat flux).
Membranes 14 00098 g003
Figure 4. Schematic mass-transfer resistances for three mass-transfer regions of membrane contactor. (Membranes 14 00098 i001Hot fluid; Membranes 14 00098 i002Cold fluid; Membranes 14 00098 i003Permeate flux).
Figure 4. Schematic mass-transfer resistances for three mass-transfer regions of membrane contactor. (Membranes 14 00098 i001Hot fluid; Membranes 14 00098 i002Cold fluid; Membranes 14 00098 i003Permeate flux).
Membranes 14 00098 g004
Figure 5. Reduction in thermal boundary polarization layers in the DCMD module. (a) Empty channel; (b) Channel with inserted turbulence promoters.
Figure 5. Reduction in thermal boundary polarization layers in the DCMD module. (a) Empty channel; (b) Channel with inserted turbulence promoters.
Membranes 14 00098 g005
Figure 6. The energy balance made within a finite fluid element. (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Figure 6. The energy balance made within a finite fluid element. (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Membranes 14 00098 g006
Figure 7. The average carbon-fiber width of S-ribs carbon-fiber turbulence promoters.
Figure 7. The average carbon-fiber width of S-ribs carbon-fiber turbulence promoters.
Membranes 14 00098 g007
Figure 8. Effects of carbon-fiber widths on temperature distributions along the module. (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Figure 8. Effects of carbon-fiber widths on temperature distributions along the module. (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Membranes 14 00098 g008
Figure 9. Effects of saline feed temperature and carbon-fiber widths on τ t e m p . (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Figure 9. Effects of saline feed temperature and carbon-fiber widths on τ t e m p . (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Membranes 14 00098 g009
Figure 10. Comparisons of theoretical Nusselt numbers for the channels with inserted S-ribs carbon-fiber turbulence promoters.
Figure 10. Comparisons of theoretical Nusselt numbers for the channels with inserted S-ribs carbon-fiber turbulence promoters.
Membranes 14 00098 g010
Figure 11. Effect of saline feed temperature and feed flow rate on permeate fluxes (3 mm). (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Figure 11. Effect of saline feed temperature and feed flow rate on permeate fluxes (3 mm). (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Membranes 14 00098 g011
Figure 12. Effect of saline feed temperature and feed flow rate on permeate fluxes (4 mm). (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Figure 12. Effect of saline feed temperature and feed flow rate on permeate fluxes (4 mm). (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Membranes 14 00098 g012
Figure 13. Effect of saline feed temperature and feed flow rate on permeate fluxes (5 mm). (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Figure 13. Effect of saline feed temperature and feed flow rate on permeate fluxes (5 mm). (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Membranes 14 00098 g013
Figure 14. Comparisons of theoretical permeate fluxes of the channels with inserted S-ribs carbon-fiber turbulence promoters and straight-line carbon-fiber filaments [29].
Figure 14. Comparisons of theoretical permeate fluxes of the channels with inserted S-ribs carbon-fiber turbulence promoters and straight-line carbon-fiber filaments [29].
Membranes 14 00098 g014
Figure 15. Effect of inlet saline feed temperature and channel designs on the economic feasibility. (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Figure 15. Effect of inlet saline feed temperature and channel designs on the economic feasibility. (a) Cocurrent-flow operations; (b) Countercurrent-flow operations.
Membranes 14 00098 g015
Table 1. Effects of operation conditions and carbon-fiber widths on N theo   and I N , T c , i n = 25 °C.
Table 1. Effects of operation conditions and carbon-fiber widths on N theo   and I N , T c , i n = 25 °C.
T h , i n
(°C)
Q h
(L/min)
Empty3 mm4 mm5 mm
N empty × 10 3 N theo × 10 3 I N N theo × 10 3 I N N theo × 10 3 I N
(kg m−2 h−1)(kg m−2 h−1)(%)(kg m−2 h−1)(%)(kg m−2 h−1)(%)
50Cocurrent-flow operations
0.30.750.8514.060.8412.580.8311.37
0.59.001.0617.461.0415.681.0314.35
0.71.051.2721.141.2519.241.2317.37
0.91.101.3926.001.3523.101.3320.55
600.31.211.4217.271.3914.921.3713.31
0.51.411.7121.281.6718.161.6416.45
0.71.602.0226.191.9521.941.9219.75
0.91.732.2831.792.1825.952.1323.12
50Countercurrent-flow operations
0.30.871.0015.440.9913.820.9812.73
0.51.011.2221.091.2018.811.1917.52
0.71.121.4226.341.3923.841.3822.77
0.91.171.5532.481.5229.491.5028.03
600.31.391.6518.561.6115.831.5914.39
0.51.571.9524.331.8920.641.8719.30
0.71.752.2930.822.2126.142.1824.31
0.91.882.5937.772.4832.132.4530.27
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.; Wang, Y.-W.; Chao, Y.; Chew, T.L.; Jiang, M.-S.; Chen, J.-H.; Li, C.-Y. Enhancing the Permeate Flux Improvement of Direct Contact Membrane Distillation Modules with Inserted S-Ribs Carbon-Fiber Filaments. Membranes 2024, 14, 98. https://doi.org/10.3390/membranes14050098

AMA Style

Ho C-D, Wang Y-W, Chao Y, Chew TL, Jiang M-S, Chen J-H, Li C-Y. Enhancing the Permeate Flux Improvement of Direct Contact Membrane Distillation Modules with Inserted S-Ribs Carbon-Fiber Filaments. Membranes. 2024; 14(5):98. https://doi.org/10.3390/membranes14050098

Chicago/Turabian Style

Ho, Chii-Dong, Yi-Wun Wang, Yi Chao, Thiam Leng Chew, Ming-Shen Jiang, Jian-Har Chen, and Ching-Yu Li. 2024. "Enhancing the Permeate Flux Improvement of Direct Contact Membrane Distillation Modules with Inserted S-Ribs Carbon-Fiber Filaments" Membranes 14, no. 5: 98. https://doi.org/10.3390/membranes14050098

APA Style

Ho, C. -D., Wang, Y. -W., Chao, Y., Chew, T. L., Jiang, M. -S., Chen, J. -H., & Li, C. -Y. (2024). Enhancing the Permeate Flux Improvement of Direct Contact Membrane Distillation Modules with Inserted S-Ribs Carbon-Fiber Filaments. Membranes, 14(5), 98. https://doi.org/10.3390/membranes14050098

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