Enhancing the Permeate Flux of Direct Contact Membrane Distillation Modules with Inserting 3D Printing Turbulence Promoters

Two geometric shape turbulence promoters (circular and square of same areas) of different array patterns using three-dimensional (3D) printing technology were designed for direct contact membrane distillation (DCMD) modules in the present study. The DCMD device was performed at middle temperature operation (about 45 °C to 60 °C) of hot inlet saline water associated with a constant temperature of inlet cold stream. Attempts to reduce the disadvantageous temperature polarization effect were made inserting the 3D turbulence promoters to promote both the mass and heat transfer characteristics in improving pure water productivity. The additive manufacturing 3D turbulence promoters acting as eddy promoters could not only strengthen the membrane stability by preventing vibration but also enhance the permeate flux with lessening temperature polarization effect. Therefore, the 3D turbulence promoters were individually inserted into the flow channel of the DCMD device to create vortices in the flow stream and increase turbulent intensity. The modeling equations for predicting the permeate flux in DCMD modules by inserting the manufacturing 3D turbulence promoter were investigated theoretically and experimentally. The effects of the operating conditions under various geometric shapes and array patterns of turbulence promoters on the permeate flux with hot inlet saline temperatures and flow rates as parameters were studied. The distributions of the fluid velocities were examined using computational fluid dynamics (CFD). Experimental study has demonstrated a great potential to significantly accomplish permeate flux enhancement in such new design of the DCMD system. The permeate flux enhancement for the DCMD module by inserting 3D turbulence promoters in the flow channel could provide a maximum relative increment of up to 61.7% as compared to that in the empty channel device. The temperature polarization coefficient (τtemp) was found in this study for various geometric shapes and flow patterns. A larger τtemp value (the less thermal resistance) was achieved in the countercurrent-flow operation than that in the concurrent-flow operation. An optimal design of the module with inserting turbulence promoters was also delineated when considering both permeate flux enhancement and energy utilization effectiveness.


Introduction
Membrane distillation (MD) has technical feasibility advantages in separation and water treatment with the high permeate flux and low energy consumption. In an MD desalination process [1,2], the existence of temperature differences was built up by the hydrophobic membrane surface temperatures that always contact both hot and cold bulk temperatures. The need for a better temperature driving-force gradient is attracting global attention due to the vapor pressure difference for the permeate passing through the hydrophobic membrane with considerable heat loss. The significant thermal resistance occurs on the thermal boundary layer close to the membrane surface of saline feed stream, which creates a temperature gradient between both membrane surfaces and results in heat loss [3]. Moreover, a higher temperature polarization coefficient yielding augmentation of drivingforce temperature gradients of both bulk feed streams comes out with the increment of trans-membrane permeate flux [4,5]. Development of desalination techniques with a higher permeate flux has been conducted for improved seawater desalination system [6,7]. Despite major advances in developing membrane materials on distillation performance [8][9][10] and pressure driven membrane separation processes are still facing the problems of the temperature polarization effect [11,12]. Previous studies in improving the permeate flux enhancement by incorporating proper flow alteration configurations such as filament [13][14][15], roughened surface [16] and eddy promoter [17] into the flow channel of DCMD modules to boost the turbulence intensity were also investigated. Having spacer filaments in the feed channel has been proven to dynamically change the thermal boundary layer and mitigate temperature polarization. It is well recognized that the turbulence promoters were added individually to the flow channel to generate vortices in enhancing the local shear stress on the membrane surface and to create secondary flows or eddies in the feed stream. Those previous studies confirmed that the mathematical model is capable of predicting heat and mass transfer behaviors in the MD desalination system compared to the experimental results. This study showcases a new design of 3D printing turbulence promoters by machining for generating turbulent intensity in the flow channel of DCMD modules. Extending previous studies, the present work focuses on the overall heat-transfer resistance in which the potential investigation of combining different geometric shapes of 3D turbulence promoters and various array patterns enhance hydrodynamic conditions. The 3D printing turbulence promoters present a higher designing flexibility so as to create multi-faceted geometric shapes under various scales with respect to traditional fabricating techniques [18,19], and this 3D printing technology was used to design spacer-filled flow channel in the membrane separation module to augment the device performance [20,21]. However, the permeate flux enhancement is associated with power consumption increments due to the augmented turbulent intensity attributed to inserting 3D turbulence promoters into the flow channel. The ratio of permeate flux enhancement to energy consumption increment was also assessed to indicate an optimal operation of considering economic feasibility.
Two kinds of 3D printing turbulence promoters, say circle and square shapes, were designed and fabricated to produce secondary flow characteristics and thus alleviate the temperature polarization effect in membrane distillation devices. The turbulence promoters inserted into the flow channel provided a larger convective heat-transfer coefficient due to disrupting the thermal boundary layer on the hot feed stream when compared with that of the device with empty channel. The objective of this study is to fabricate two geometric shapes turbulence promoters and stick onto the membrane surface of the flow channel in hot saline water stream to reduce its temperature polarization effect owing to reducing the heat transfer resistance, and thus, a higher permeate flux was achieved. A theoretical model validated by the experimental data was also developed to predict the pure water productivity for comparisons under various operating conditions, and the fluid velocity distribution was also presented using Ansys software (2018).

Experimental Setup and Materials
The schematic detailed configuration and the fabrication structure of the flat-plate DCMD module are illustrated in Figures 1 and 2, respectively, and a photo of the present experimental setup is shown in Figure 3 with acrylic plates as outside walls on a parallelplate channel. The DCMD module contains two flow channels with inserting 3D turbulence promoters into the hot feed flow channel and empty channel of the cold stream. The empty channel is constructed with a 0.1 mm nylon fiber routed as a supporting material. The 3D printing turbulence promoters were constructed with a 1mm-thick and glued with Cyanoacrylate Adhesive (Chang Chun Plastics Co., Ltd., Taiwan) on the acrylic plate of the hot feed side in contact with the membrane surface acting as eddy promoters. The printing material of turbulence promoters was made with polyester elastomer and attached onto the hydrophobic membrane surface to offer the eddy motion around those obstacles. The length, width, and clearance of the permeate channel between the two plates were 21 cm, 29 cm, and 1 mm, respectively. The top view of two geometric shapes of turbulence promoters is specified in Figure 4a,b, in which the array pattern of the turbulence promoters was used as a configuration parameter.

Experimental Setup and Materials
The schematic detailed configuration and the fabrication structure of the flat-plate DCMD module are illustrated in Figures 1 and 2, respectively, and a photo of the present experimental setup is shown in Figure 3 with acrylic plates as outside walls on a parallelplate channel. The DCMD module contains two flow channels with inserting 3D turbulence promoters into the hot feed flow channel and empty channel of the cold stream. The empty channel is constructed with a 0.1 mm nylon fiber routed as a supporting material. The 3D printing turbulence promoters were constructed with a 1mm-thick and glued with Cyanoacrylate Adhesive (Chang Chun Plastics Co., Ltd., Taiwan) on the acrylic plate of the hot feed side in contact with the membrane surface acting as eddy promoters. The printing material of turbulence promoters was made with polyester elastomer and attached onto the hydrophobic membrane surface to offer the eddy motion around those obstacles. The length, width, and clearance of the permeate channel between the two plates were 21 cm, 29 cm, and 1 mm, respectively. The top view of two geometric shapes of turbulence promoters is specified in Figure 4a,b, in which the array pattern of the turbulence promoters was used as a configuration parameter.   The schematic detailed configuration and the fabrication structure of the flat-plate DCMD module are illustrated in Figures 1 and 2, respectively, and a photo of the present experimental setup is shown in Figure 3 with acrylic plates as outside walls on a parallelplate channel. The DCMD module contains two flow channels with inserting 3D turbulence promoters into the hot feed flow channel and empty channel of the cold stream. The empty channel is constructed with a 0.1 mm nylon fiber routed as a supporting material. The 3D printing turbulence promoters were constructed with a 1mm-thick and glued with Cyanoacrylate Adhesive (Chang Chun Plastics Co., Ltd., Taiwan) on the acrylic plate of the hot feed side in contact with the membrane surface acting as eddy promoters. The printing material of turbulence promoters was made with polyester elastomer and attached onto the hydrophobic membrane surface to offer the eddy motion around those obstacles. The length, width, and clearance of the permeate channel between the two plates were 21 cm, 29 cm, and 1 mm, respectively. The top view of two geometric shapes of turbulence promoters is specified in Figure 4a,b, in which the array pattern of the turbulence promoters was used as a configuration parameter.   of both sides with using a rotameter while the temperature kept 25 °C of the cold stream. Comparisons were made of permeate fluxes under various operation conditions to study the device performances between two modules with/without inserting 3D printing turbulence promoters. The experimental run of permeate flux condensed in the cold stream side was then collected and weighted using an electronic balance to measure and record the amount of permeate flux on the PC.   °C) regulated by the thermostat, and controlling various flow rates (0.1, 0.2, 0.3, 0.4 L/ of both sides with using a rotameter while the temperature kept 25 °C of the cold str Comparisons were made of permeate fluxes under various operation conditions to s the device performances between two modules with/without inserting 3D printing tu lence promoters. The experimental run of permeate flux condensed in the cold stream was then collected and weighted using an electronic balance to measure and record amount of permeate flux on the PC.   Between the hydrophobic composite membrane and the acrylic plate is inserting turbulence promoters into a 1 mm-thick silicon rubber which sealed and glued upon the acrylic plate to form a channel and to prevent leakage. Nylon fibers of 0.2 mm diameter was implemented to wind upon the hydrophobic membrane surface as a supporter grid in preventing from membrane bending and wrinkling. The artificial saline water of 3.5 wt% NaCl was prepared by adding inorganic salts NaCl into distilled water. The experiments were conducted for adjusting various inlet hot saline temperatures (40, 45, 50, 55 • C) regulated by the thermostat, and controlling various flow rates (0.1, 0.2, 0.3, 0.4 L/min) of both sides with using a rotameter while the temperature kept 25 • C of the cold stream. Comparisons were made of permeate fluxes under various operation conditions to study the device performances between two modules with/without inserting 3D printing turbulence promoters. The experimental run of permeate flux condensed in the cold stream side was then collected and weighted using an electronic balance to measure and record the amount of permeate flux on the PC.

Mass and Heat Transfer
A theoretical formulation coupling both heat and mass transfer mechanisms for the DCMD system, as shown Figure 5, to study how hot saline water was vaporized at the pore entrances of the membrane surface according to the vapor-water equilibrium, then vapor only diffused through porous hydrophobic membranes and condensed in the cold stream as pure water productivity. A mass transfer model is needed to describe the concentration of the hot saline water feed while the heat transfer model is required to determine the temperature gradients at each membrane surface due to the permeate flux transferring. Balances of permeate flux by mass diffusion and enthalpy flow conservation by heat conduction were formulated simultaneously as follows.
Membranes 2021, 11, x FOR PEER REVIEW 6 of 22 determined by the thermal conductivities of vapor in the membrane pore g k and the solid membrane material s k by Warner [24] as: The combinations of the heat flux term from Equations (1), (2) and (7) leads to the overall heat transfer coefficient of the membrane as follows:

Temperature Polarization
The temperature polarization coefficient temp  is an indicator to reveal the magnitude of the thermal boundary-layer resistance which governs the permeate flux of pure water productivity, and is commonly used to express as the ratio of membrane surface temperatures' gradient to bulk temperatures' gradient as follows: The one-dimensional mathematical treatments were developed among all heat transfer regions under steady-state operations of with assuming well insulation on the outside of the DCMD module, as illustrated by the schematic diagram in Figure 5 Both membrane surface temperatures ( Equations (8) and (10) ), respectively, as follows: The membrane permeation coefficient (c m ) and the trans-membrane saturation vapor pressure difference (∆P) have been used extensively in mass transfer analysis of permeate flux for membrane distillation processes [22,23] where the membrane permeation coefficient is the addition of Knudsen diffusion and Poiseuille flow.
The saturation vapor pressure (P sat 1 ) of the hot fluid side can be estimated using water activity coefficient (a w ), which can be determined using a correlation, Membranes 2021, 11, 266 6 of 21 and the tortuosity (τ) can be estimated using the porosity of the membrane [8] Next, the heat balance of the DCMD module could be described on the principle of the non-isothermal process with making the energy balance of enthalpy flow conservation, and the simultaneous occurrences of heat and mass transfer in the DCMD module in each heat transfer region (a) the hot saline water feed; (b) the hydrophobic composite membrane and (c) the cooling water. The energy balance equations were derived for each heat transfer region under the steady-state operation according to the schematic diagram of the DCMD module in Figure 5 and were shown in Equations (6)-(8) as follows: where N λ is characterized as the latent heat of vaporization and is the conductive heat transfer, 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 [24] as: The combinations of the heat flux term from Equations (1), (2) and (7) leads to the overall heat transfer coefficient of the membrane as follows:

Temperature Polarization
The temperature polarization coefficient τ temp is an indicator to reveal the magnitude of the thermal boundary-layer resistance which governs the permeate flux of pure water productivity, and is commonly used to express as the ratio of membrane surface temperatures' gradient to bulk temperatures' gradient as follows: The one-dimensional mathematical treatments were developed among all heat transfer regions under steady-state operations of q = q h = q m = q c with assuming well insulation on the outside of the DCMD module, as illustrated by the schematic diagram in Figure 5 Both membrane surface temperatures (T 1 and T 2 ) and the convective heattransfer coefficients (h h ) were obtained by equating Equations (6) and (10) (q h = q m ) and Equations (8) and (10) (q m = q c ), respectively, as follows: Then, a more simplified form of τ temp expressed in terms of the heat transfer coefficient Membranes 2021, 11, 266 7 of 21 The procedure for calculation of theoretical values of the heat transfer coefficient will be described as follows. First, with the given operation conditions, the heat transfer coefficient is determined from Equations (12) and (13). 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 (12) once T 2 (or T 1 ) is assumed in Equation (13). Further, the convective heat-transfer coefficient is calculated from Equation (1). With this calculated value of the convective heat-transfer coefficient, new values of T 1 and T 2 are then recalculated from Equations (12) and (13). 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. The modeling equations of the energy balances were obtained by making the energy-flow diagram presented in a finite fluid element, as shown in Figure 6, to solve the temperature distributions of both hot and cold stream sides as The 3D printing turbulence promoters inserted into flow channel play a significa role as the eddy promoter, and a better interpretation of both heat and mass transfer haviors in the turbulent flow could be described by a new method based on dimensio analysis of the Buckingham's  theorem. The Nusselt number of flow channels w inserting 3D printing turbulence promoters can be related to dimensionless groups:

Power Consumption Increment
The power consumption increment is unavoidable due to inserting 3D printing t bulence promoters into the hot saline water flowing channel as eddy promoters. T power consumption only the friction losses to walls of both hot and cold streams we significant and may be determined using Fanning friction factor The average velocity and equivalent hydraulic diameter of cold and hot stream sid are defined as follows: The calculated convective heat-transfer coefficients were conveyed to solve the temperature distributions of both hot and cold stream sides in two simultaneous ordinary differential equations of Equations (15) and (16) for concurrent-flow operation (or Equation (17) for countercurrent-flow operation) by marching the fourth order Runge-Kutta method along the length of the module.

Enhancement Factor
The 3D printing turbulence promoters are inserted in the conduit of hot feed stream side instead of using the device of empty channel. The enhancement factor α p depending on the geometric shapes and array patterns was corrected to calculate the augmented convective heat transfer coefficients in DCMD modules with inserting the 3D printing turbulence promoters and were carried out an iterative procedure as [25][26][27] where The 3D printing turbulence promoters inserted into flow channel play a significant role as the eddy promoter, and a better interpretation of both heat and mass transfer behaviors in the turbulent flow could be described by a new method based on dimensional analysis of the Buckingham's π theorem. The Nusselt number of flow channels with inserting 3D printing turbulence promoters can be related to dimensionless groups: where W p and D h,h are the average equivalent width of 3D printing turbulence promoters and hydraulic diameter of the hot stream side, respectively.

Power Consumption Increment
The power consumption increment is unavoidable due to inserting 3D printing turbulence promoters into the hot saline water flowing channel as eddy promoters. The power consumption only the friction losses to walls of both hot and cold streams were significant and may be determined using Fanning friction factor f F , The average velocity and equivalent hydraulic diameter of cold and hot stream sides are defined as follows: The Fanning friction factor can be estimated using a correlated equation based on channel's aspect ratio (β = d/W), The relative extents I P of power consumption increment was illustrated by calculating the percentage increment in the device with inserting turbulence promoters, based on the device of empty channel as where the subscripts promoter and empty represent the channels with inserting turbulence promoters and the empty channel.

The Design of 3D Printed Turbulence Promoters
A smart attempt to advance hydrodynamic conditions in pressure-driven membrane separation processes with gluing 3D printing turbulence promoters on the membrane surface can accomplish the goal of very economical process in terms of implementing cost and speed. The 3D printing technology presents a higher flexibility in designing various complex geometric shapes compared to traditional manufacturing techniques [28]. This is one of the main advantages of 3D printing processes to produce different shapes of turbulence promoters [19] in precisely tailoring the 3D geometric shapes through a layer-by-layer machining process. Two types of turbulence promoters were designed as circular and square by 3D printing technology and empty channel (wound with nylon fiber), respectively, and inserted in the flowing channel throughout the experiments, as shown in Figure 7 Circular type was made of a circle with a diameter of 29.5 mm. Square type was made of a square with each length of 25.69 mm, all of the turbulence promoters had the same area and volume. The height of flowing channel and turbulence promoters are both 1 mm. The turbulence promoters were printed by a 3D printer (ATOM 2.5EX, Taiwan) in making three dimensional solid objects from a digital file of different shapes and incorporated into the flow path and glued to the membrane surface. The printing filament was fabricated from polyester elastomer to durable and hydrophobic promoter structures. Two kinds of array patterns with four half shape and ten full shape turbulence promoters were installed in the hot stream side. Since the printed promoter icon occupied 13% coverage of the entire membrane permeate area, therefore the effective permeate area of the membrane in blocking vapor flux is taken into account in calculation procedure.
Membranes 2021, 11, x FOR PEER REVIEW 1 moters with skewness of 0.8, respectively. The density and viscosity for the fluid w sumed to be 1019 kg/m 3 and 1.06 cp, respectively. The inlet velocity was given by lating the volumetric flow rate and the inlet cross area, and the distribution of cor velocity was obtained from CFD simulation with conserving the voidage change inserting turbulence promoters, as shown in Figure 8.

Numerical Scheme
The temperatures on both sides of membrane surfaces ( 1 T and 2 T ) were calc by equating Equations (6) and (10) and Equations (8) and (10)

CFD Simulation
The detail of the 3D printing protocol of the fabricated promoters is included in Figure 7 The permeate flux area was 6.09 × 10 −2 mA commercial CFD solver was used to solve the velocity profile of the Navier-Stokes equations for the incompressible flow coupled using the CFD solver algorithm under steady state and the volume of fluid model. The convergence criteria were the residuals of continuity and velocity below 10 − The geometry of the flow channel and 3D mesh were illustrated using Ansys software (2018). Hexahedral elements were appropriate for this structure including a wedge element. The mesh sizes were approximately 78334, 77900, 73361 and 72433 for each type of turbulence promoters with skewness of 0.8, respectively. The density and viscosity for the fluid was assumed to be 1019 kg/m 3 and 1.06 cp, respectively. The inlet velocity was given by calculating the volumetric flow rate and the inlet cross area, and the distribution of corrected velocity was obtained from CFD simulation with conserving the voidage change due to inserting turbulence promoters, as shown in Figure 8.

Flux Enhancement by Inserting Turbulence Promoter in DCMD Module
The permeate flux is dependent on the temperature gradient between the membrane surface 1 T and 2 T in the DCMD system. The effect of the geometric turbulence promoters on the longitudinal temperature profiles of both channels in the DCMD module was shown in Figure 9. The devices with inserting turbulence promoters with various geometric shapes and array patterns for eddy promoting come out with temperature polarization reduction. Both membrane-surface and bulk temperatures trend taper along the flowing direction in concurrent-flow operations with decreasing the driving-force temperature gradients shown in Figure 9 while the driving-force temperature gradients in countercurrent-flow operations keep a relative higher average value than those in concurrent-flow operations, and thus, the permeate flux enhancement in countercurrent-flow operations was accomplished. Reduction of the temperature polarization effect was achieved using eddy promoters in flow channel for a favorable result investigated by computational simulation [29]. The heat transfer coefficients of inserting turbulence promoters in flow channel expressed by

Numerical Scheme
The temperatures on both sides of membrane surfaces (T 1 and T 2 ) were calculated by equating Equations (6) and (10) and Equations (8) and (10) with the initial guess of the convective heat-transfer coefficient (h h ) until the iteration procedure reaches the convergence with 0.1% error of accuracy control and validated by the experimental flux N exp under the known inlet and outlet temperatures of both stream sides. The estimated values of membrane surfaces and the convective heat-transfer coefficient were used to solve Equations (13) and (14) (Equation (15) for countercurrent-flow operations) numerically using the fourth order Runge-Kutta method along the flowing direction of 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 but also on the membrane surfaces of both hot and cold feed sides under concurrent-and countercurrent-flow operations, respectively. Comparisons were made for the permeate flux of the channel with inserting turbulence promoters and the empty channel under both concurrent-and countercurrent-flow operations.

Flux Enhancement by Inserting Turbulence Promoter in DCMD Module
The permeate flux is dependent on the temperature gradient between the membrane surface T 1 and T 2 in the DCMD system. The effect of the geometric turbulence promoters on the longitudinal temperature profiles of both channels in the DCMD module was shown in Figure 9. The devices with inserting turbulence promoters with various geometric shapes and array patterns for eddy promoting come out with temperature polarization reduction. Both membrane-surface and bulk temperatures trend taper along the flowing direction in concurrent-flow operations with decreasing the driving-force temperature gradients shown in Figure 9 while the driving-force temperature gradients in countercurrent-flow operations keep a relative higher average value than those in concurrent-flow operations, and thus, the permeate flux enhancement in countercurrent-flow operations was accomplished. Inserting turbulence promoters interrupts the thermal boundary layer on the membrane surface that diminishes heat transfer resistance; hence the permeate flux was boosted. The permeate fluxes were calculated through the enhancement factor, Equation (27), for predicting the Nusselt number, as referred to the heat transfer efficiency, with inserting turbulence promoters in flow channel. A relative increase of permeate flux E I was used to make comparisons between the permeate flux of the channel with turbulence promoters and the empty channel. Reduction of the temperature polarization effect was achieved using eddy promoters in flow channel for a favorable result investigated by computational simulation [29]. The heat transfer coefficients of inserting turbulence promoters in flow channel expressed by Nusselt number and determined by the Buckingham's π theorem, were correlated with the experimental data. The enhancement factor α p derived from the correlation of Nusselt number for the channel with inserting turbulence promoters was determined via a regression analysis as: The Nusselt numbers for the empty channel are in linear diagonal straight line with the experimental data, as shown in Figure 10, while the correlation expression of Nusselt numbers is also applicable for the channel with inserting two geometric shapes and array patterns of turbulence promoters. The Nusselt number in the channels with inserting turbulence promoters is higher than that of the empty channel.
brane surface that diminishes heat transfer resistance; hence the permeate flux was boosted. The permeate fluxes were calculated through the enhancement factor, Equation (27), for predicting the Nusselt number, as referred to the heat transfer efficiency, with inserting turbulence promoters in flow channel. A relative increase of permeate flux E I was used to make comparisons between the permeate flux of the channel with turbulence promoters and the empty channel.  Inserting turbulence promoters interrupts the thermal boundary layer on the membrane surface that diminishes heat transfer resistance; hence the permeate flux was boosted. The permeate fluxes were calculated through the enhancement factor, Equation (27), for predicting the Nusselt number, as referred to the heat transfer efficiency, with inserting turbulence promoters in flow channel. A relative increase of permeate flux I E was used to make comparisons between the permeate flux of the channel with turbulence promoters and the empty channel.
The performance of permeate flux enhancement I E was illustrated by calculating the percentage increase in the device with inserting turbulence promoters, based on the device of empty channel as The comparisons have shown that the increased I E for the channel with inserting turbulence promoters of two geometric shapes and two array patterns in both concurrent-and countercurrent-flow operations. In general, the permeate flux enhanced by the inserting turbulence promoters is more significant in countercurrent-flow operations than those in concurrent-flow operations. As expected, either the increase of hot feed flow rate or the use of countercurrent-flow operations of turbulence promoters results in more permeate flux. The theoretical permeate flux for various hot feed flow rates and inlet temperatures under concurrent-and countercurrent-flow operations are summarized in Tables 1 and 2, respectively. It is also seen from Table 1; Table 2 that the order of the performance of permeate flux enhancement for the device inserting turbulence promoters is the following: Type B Square > Type B Circle > Type A Circle > Type A Square. The Square turbulence promoter with Type B array patterns enhances the permeate flux enhancement by approximately 20-40% compared to Type A, whereas the Circle turbulence promoters of Type A shows a lower permeate flux enhancement than that in Type B by approximately 5-10%, under the same operation conditions. Overall, inserting turbulence promoters into flow channel shows a great potential to significantly enhance the permeate flux in the DCMD module. Thus, temperature polarization reduction is crucial not only in pressure-driven membrane distillation processes but also in other water treatment technologies.   Velocities and vortices in enhancing the local shear stress on the membrane surface to create secondary flows or eddies in the feed stream are two factors to affect the permeate flux enhancement of the device inserting square and circle turbulence promoters. The average equivalent width of square and circle turbulence promoters is 9.72 × 10 −3 m and 1.22 × 10 −2 m, respectively, and hence the velocities in the device with inserting square turbulence promoters are larger than those in the circle ones. Besides, the device with inserting geometric shapes of square turbulence promoters in flowing channel generates more intensive vortices and eddies than those in the device with inserting circle turbulence promoters, resulting the square turbulence promoter is better than circle.
The permeate flux enhancement with inserting turbulence promoters in flow channel is illustrated in Figures 11 and 12 for both experimental results and the theoretical predictions. Two shapes and two array patterns of the inserting turbulence promoters produces the higher turbulence intensity that results in the higher heat transfer or, the higher permeate flux. Comparisons of theoretical predictions of permeate flux were made for both concurrent-and countercurrent-flow operations as depicted in Figure 11 and The agreement of the simulation results with those obtained from experimental runs is pretty good, as indicated in Figure 11 and As expected, either the increase of hot saline water flow rate or the countercurrent-flow operation results in more permeate flux. A relative increase of permeate flux I E was used to compare the permeate flux of the channel with inserting turbulence promoter to the empty channel. In general, the permeate flux enhanced by inserting turbulence promoters is more significant in countercurrent-flow operations than that in concurrent-flow operations. as indicated in Figure 11 and As expected, either the increase of hot saline water flow rate or the countercurrent-flow operation results in more permeate flux. A relative increase of permeate flux E I was used to compare the permeate flux of the channel with inserting turbulence promoter to the empty channel. In general, the permeate flux enhanced by inserting turbulence promoters is more significant in countercurrent-flow operations than that in concurrent-flow operations.  This study also measured, predicted and compared the effects of geometric shapes of turbulence promoters on temperature polarization for countercurrent-flow operations, as depicted in Figure 13. Inserting turbulence promoters into the hot feed stream results in a higher temp  value (a lesser temperature polarization) was found due to the reduction of the thermal boundary-layer thickness as compared to the temp  value of the module without inserting turbulence promoters, as shown in Figure 13. In addition, a larger temp  value (the less thermal resistance) was achieved in the countercurrent-flow operation than that in the concurrent-flow operation.

Analysis of Cake Properties
The precision index of experimental uncertainty of each individual measurement of permeate flux is calculated as referred to Moffat [30] directly from the experimental runs as follows:  (30) and the uncertainty of the reproducibility of the permeate fluxes will be associated with the mean precision index This study also measured, predicted and compared the effects of geometric shapes of turbulence promoters on temperature polarization for countercurrent-flow operations, as depicted in Figure 13. Inserting turbulence promoters into the hot feed stream results in a higher τ temp value (a lesser temperature polarization) was found due to the reduction of the thermal boundary-layer thickness as compared to the τ temp value of the module without inserting turbulence promoters, as shown in Figure 13. In addition, a larger τ temp value (the less thermal resistance) was achieved in the countercurrent-flow operation than that in the concurrent-flow operation.   E . The present work extends the existing study except for inserting 3D turbulence promoters instead of W-shaped carbon-fiber spacers [31]. Besides, to perform additional DCMD tests with the canals filled with supporting mesh, the experiment runs were conducted the channels of DCMD devices with inserting diagonal carbon-fiber spacers to replace the W-shaped carbon-fiber spacers, as shown in Figure 14. Therefore, the graphical representation for comparisons with theoretical predictions of Nusselt numbers obtained in the present study, diagonal carbon-fiber spacers and W-shaped carbon-fiber spacers [31] illustrates why the present design of inserting 3D turbulence promoters is preferred, presented by Figure 15 for countercurrent-flow operations. This is the value and originality of the present study regarding to the implementing cost and technical feasibility.

Energy Consumption
This study further evaluated the device design's effectiveness by considering the ratio of permeate flux increment to power consumption increment P E I I / (Equations (28) and (30)) due to the flow resistance with more power consumptions caused by the insert-

Analysis of Cake Properties
The precision index of experimental uncertainty of each individual measurement of permeate flux is calculated as referred to Moffat [30] directly from the experimental runs as follows: (31) and the uncertainty of the reproducibility of the permeate fluxes will be associated with the mean precision index.
The mean precision index of the experimental measurements of permeate flux enhancement was evaluated within 8.12 × 10 −3 ≤ S N ≤ 2.20 × 10 −2 for both concurrent-flow and countercurrent-flow operations. The permeate flux's experimental results prove the theoretical predictions' validity by defining the accuracy between the numerical solutions and the experimental results as follows: where N theo indicates the theoretical prediction of permeate flux while N exp and N exp are the number of experimental measurements and the experimental data of permeate flux. The error analysis of the experimental measurements determined by Equation (33) for both the theoretical prediction is 0.11 ≤ E ≤ 9.68. The present work extends the existing study except for inserting 3D turbulence promoters instead of W-shaped carbon-fiber spacers [31]. Besides, to perform additional DCMD tests with the canals filled with supporting mesh, the experiment runs were conducted the channels of DCMD devices with inserting diagonal carbon-fiber spacers to replace the W-shaped carbon-fiber spacers, as shown in Figure 14. Therefore, the graphical representation for comparisons with theoretical predictions of Nusselt numbers obtained in the present study, diagonal carbon-fiber spacers and W-shaped carbon-fiber spacers [31] illustrates why the present design of inserting 3D turbulence promoters is preferred, presented by Figure 15 for countercurrent-flow operations. This is the value and originality of the present study regarding to the implementing cost and technical feasibility.
ing turbulence promoters in the flow channel. The effect of feed flow rate, two array patterns of square turbulence promoters, inlet hot feed temperature, and concurrent-/countercurrent-flow operations on P E I I / are illustrated in Figure 16. The 3D printing turbulence promotor inserted in the flow channel is the aim to achieve the augmented turbulence intensity in lessening temperature polarization layers and amplify the convective heat-transfer coefficient as well, and thus, the permeate flux is boosted. The desirable permeate flux improvement and the undesirable power consumption increment are two conflict effects encountered in making economic consideration by inserting turbulence promoters in the flow channel. The percentage of permeate flux enhancement is less than that of energy consumption increment results in a lower value of P E I I / , which indicates that power consumption increment cannot compensate more permeate flux due to a limited vapor transporting rate through the membrane by ing turbulence promoters in the flow channel. The effect of feed flow rate, two array patterns of square turbulence promoters, inlet hot feed temperature, and concurrent-/countercurrent-flow operations on P E I I / are illustrated in Figure 16. The 3D printing turbulence promotor inserted in the flow channel is the aim to achieve the augmented turbulence intensity in lessening temperature polarization layers and amplify the convective heat-transfer coefficient as well, and thus, the permeate flux is boosted. The desirable permeate flux improvement and the undesirable power consumption increment are two conflict effects encountered in making economic consideration by inserting turbulence promoters in the flow channel. The percentage of permeate flux enhancement is less than that of energy consumption increment results in a lower value of P E I I / , which indicates that power consumption increment cannot compensate more permeate flux due to a limited vapor transporting rate through the membrane by

Energy Consumption
This study further evaluated the device design's effectiveness by considering the ratio of permeate flux increment to power consumption increment I E /I P (Equations (28) and (30)) due to the flow resistance with more power consumptions caused by the inserting turbulence promoters in the flow channel. The effect of feed flow rate, two array patterns of square turbulence promoters, inlet hot feed temperature, and concurrent-/countercurrentflow operations on I E /I P are illustrated in Figure 16.
values. Notice that the P E I I / ratio of the countercurrent-flow operation is higher than that of the concurrent-flow operation. Therefore, comparisons on both operations of flow patterns and array patterns were made on considering both the effective utilization of power consumption relative to the increase of permeate flux to indicate the trend of economic feasibility with inserting turbulence promoters for some specific hot feed rates in this study. The introduction of inserting 3D turbulence promoters has positive effects on the heat and mass transfer rate, and outcome assessments of permeate flux improvements requires determining some key operating parameters, such as flow rates, geometric shapes, array patterns and flow patterns. The decrease in temperature polarization (A larger temp  ) with inserting 3D turbulence promoters is accomplished in tapering thermal boundary layers, and hence, establishing a larger trans-membrane mass transfer across the hydrophobic membrane due to the driving-force temperature gradients increment. Meanwhile, a suitable selection of the operating parameters for the device with inserting 3D turbulence promoters would be come out an evident advantage for practical applications with the economic sense which allows the specification setting by the designer. Actually, the present study is an extension of our previous work of DCMD module [31] to apply the device with inserting 3D turbulence promoters instead of carbon-fiber spacers. Though the phenomenon of heat and mass transfer behaviors in the present study could be analogized from The 3D printing turbulence promotor inserted in the flow channel is the aim to achieve the augmented turbulence intensity in lessening temperature polarization layers and amplify the convective heat-transfer coefficient as well, and thus, the permeate flux is boosted. The desirable permeate flux improvement and the undesirable power consumption increment are two conflict effects encountered in making economic consideration by inserting turbulence promoters in the flow channel. The percentage of permeate flux enhancement is less than that of energy consumption increment results in a lower value of I E /I P , which indicates that power consumption increment cannot compensate more permeate flux due to a limited vapor transporting rate through the membrane by the convective heat-transfer coefficient in the hot stream side for the device with inserting turbulence promoters. In other words, the array pattern of Type B gives higher value of I E /I P , which reflects at the expenses of energy consumption are more effective in increasing the permeate flux. Running countercurrent-flow operations with a larger temperature gradient between hot saline feed and cold side than the concurrent-flow operations gives the higher I E /I P value. The comparison reveals that the countercurrent-flow operation obtains permeate flux more effectively than that in the concurrent-flow operation.
As the I E /I P values decrease with the increasing hot feed flow rate larger than 5.0 × 10 6 m 3 /s for both concurrent-and countercurrent-flow operations and array patterns of square turbulence promoters, one may notice that there are existing an optimal I E /I P values. Notice that the I E /I P ratio of the countercurrent-flow operation is higher than that of the concurrent-flow operation. Therefore, comparisons on both operations of flow patterns and array patterns were made on considering both the effective utilization of power consumption relative to the increase of permeate flux to indicate the trend of economic feasibility with inserting turbulence promoters for some specific hot feed rates in this study.
The introduction of inserting 3D turbulence promoters has positive effects on the heat and mass transfer rate, and outcome assessments of permeate flux improvements requires determining some key operating parameters, such as flow rates, geometric shapes, array patterns and flow patterns. The decrease in temperature polarization (A larger τ temp ) with inserting 3D turbulence promoters is accomplished in tapering thermal boundary layers, and hence, establishing a larger trans-membrane mass transfer across the hydrophobic membrane due to the driving-force temperature gradients increment. Meanwhile, a suitable Membrane surface on cold feed side 2 Membrane surface on hot feed side h In the hot feed flow channel c In the cold feed flow channel promoter Inserting turbulence promoters empty Inserting nylon fiber as supporters exp Experimental results in Inlet lam Empty channel out Outlet theo Theoretical predictions Superscripts p The channel with inserting turbulence promoters