Next Article in Journal
Modeling and Optimal Operating Conditions of Hollow Fiber Membrane for CO2/CH4 Separation
Previous Article in Journal
Minimizing Area-Specific Resistance of Electrochemical Hydrogen Compressor under Various Operating Conditions Using Unsteady 3D Single-Channel Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Device Performance of a Tubular Membrane Dialyzer Incorporating Ultrafiltration Effects on the Dialysis Efficiency

1
Department of Chemical and Materials Engineering, Tamkang University, Tamsui, New Taipei 251301, Taiwan
2
Department of Fundamental and Applied Sciences, HICoE-Centre for Biofuel and Biochemical Research, Institute of Self-Sustainable Building, Universiti Teknologi PETRONAS, Seri Iskandar 32610, Perak Darul Ridzuan, Malaysia
*
Author to whom correspondence should be addressed.
Membranes 2023, 13(6), 556; https://doi.org/10.3390/membranes13060556
Submission received: 29 April 2023 / Revised: 19 May 2023 / Accepted: 26 May 2023 / Published: 28 May 2023
(This article belongs to the Section Membrane Applications)

Abstract

:
Membrane dialysis is one of the membrane contactors applied to wastewater treatment. The dialysis rate of a traditional dialyzer module is restricted because the solutes transport through the membrane only by diffusion, in which the mass-transfer driving force across the membrane is the concentration gradient between the retentate and dialysate phases. A two-dimensional mathematical model of the concentric tubular dialysis-and-ultrafiltration module was developed theoretically in this study. The simulated results show that the dialysis rate improvement was significantly improved through implementing the ultrafiltration effect by introducing a trans-membrane pressure during the membrane dialysis process. The velocity profiles of the retentate and dialysate phases in the dialysis-and-ultrafiltration system were derived and expressed in terms of the stream function, which was solved numerically by the Crank–Nicolson method. A maximum dialysis rate improvement of up to twice that of the pure dialysis system ( V w = 0 ) was obtained by employing a dialysis system with an ultrafiltration rate of V w = 2   m L / m i n and a constant membrane sieving coefficient of θ = 1 . The influences of the concentric tubular radius, ultrafiltration fluxes and membrane sieve factor on the outlet retentate concentration and mass transfer rate are also illustrated.

1. Introduction

Membrane extraction [1,2], membrane absorption [3,4], ion-exchange membranes [5,6] and membrane distillation [7,8] are membrane-based separation processes [9] which have been widely applied to many industries over decades in concentration or purification processes, in which the important features of low energy consumption, continuous operation and easy up-scaling [10,11] were the main advantages of these membrane technologies. However, membrane dialysis is unique among the membrane contactors in that the solvent may transport through the membrane accompanied by the solutes from the feed stream (defined as the retentate phase) to the receiving stream (defined as the dialysate phase) by diffusion due to the driving force of the concentration gradient. An early use of membrane dialysis was Donnan dialysis [12] and alcohol reduction in beverages [13], and the most famous applications of hemodialysis (artificial kidney) are conducted to remove waste metabolic end products from the human body [14,15]. Hemodialysis is a process of purifying blood through hemodialyzers containing a hemodialysis membrane, which is the most effective treatment to prolong the life of people suffering from end-stage renal diseases (ESRD) [16].
The dialysate phase concentration in a dialysis process does not remain unchanged and the solvent is inevitably to pass through the membrane, across which exists a pressure gradient between two phases besides the other membrane contactors. Early studies focused on a simplified analysis by ignoring the concentration variance of the dialysate phase and assuming no solvent passed through the membrane [17,18]. A two-dimensional theoretical model was developed [19] to simulate the mass transfer in hemodialyzers considering concentration variations and the overall mass-transfer coefficients were determined in an unsteady-state recycle dialyzer by neglecting the concentration change [20]. The effect of the ultrafiltration operation on the mass-transfer efficiency improvement was studied by simultaneously diffusing and convecting the solute in a dialyzer while assuming there was no solute concentration in the dialysate phase [21] and nonzero dialysate concentration [22]. Henderson [23] and Leber et al. [24] have shown that the ultrafiltration operation not only can improve the mass transfer efficiency in a hemodialysis process but also can reduce the treatment time; in addition, it is tolerated by patients. Furthermore, Yeh et al. [25] investigated the mass transfer in a cross-flow parallel-plate dialyzer coupled with ultrafiltration by using the perturbation method. A one-dimensional theoretical model of a concurrent-flow [26] and countercurrent-flow [27] parallel-plate dialyzer coupled with ultrafiltration was developed and solved by using the Frobenius method. Successful intensifications of membrane dialyzers have been developed and employed in anion-exchange membranes for acid recovery from acidic wastewater [28] and heparin-free hemodialysis [29,30]. Previous studies proved some functionalized hemodialysis membranes were used to overwhelm oxidative stress [31] hemodialysis therapy, in which utilizing a tannic acid (TA) coating as a reactive oxygen species was employed in membrane dialysis indicating an increase in the serum total antioxidant capacity. Moreover, the anion-exchange membranes for acid recovery depend on the dialysis efficiency and selectivity of membrane gas absorption, distribution coefficient and a composition gradient of gas solute in the gas/liquid system [32].
The phenomena of suction or injection through the membrane due to the ultrafiltration operation will lead to a complex flowing pattern in the channels of a membrane dialyzer, in which the velocity distribution can be solved using the stream function coupled with the perturbation method [33,34]. A two-dimensional mathematical model of the concentric tubular dialysis-and-ultrafiltration module is developed theoretically. Though the mass-transfer characteristics in the present study could be analogized from those in our previous work on the flat-plate module [35], the manners of velocity and transport across the membrane are somehow different and the mass-transfer mathematical formulation is relatively complicated. This study actually extends the one-dimensional modeling of the flat-plate module to a two-dimensional analysis of the concentric tubular dialysis-and-ultrafiltration system in the present study to obtain improvements in dialysis efficiency and rate; however, no theoretical study in its simplified expression has been reported in the literature for a concentric tubular dialyzer with an ultrafiltration operation. The general applicability of such a mathematical formulation was investigated in a manner analogous to the parallel-plate dialysis-and-ultrafiltration system [26] to analyze the concentric tubular dialyzer module in this study. The purpose of the present study is to discuss the effects of the retentate phase flow rate, dialysate phase flow rate, ultrafiltration rate, membrane sieving coefficient and channel thickness ratio on the dialysis efficiency in a concentric tubular dialysis-and-ultrafiltration system. This is the value of the present study to delineate a two-dimensional concentric tubular dialyzer under ultrafiltration operation with such a simplified mathematical formulation.

2. Mathematical Formulations

2.1. Dialysis-and-Ultrafiltration Operation in a Concentric Tubular Module

A hydrophilic membrane with thickness δ was inserted in parallel as the inner tube into a circular tube with length L to divide the open conduit into parts: an inner subchannel a (retentate phase) and annular subchannel b (dialysate phase) with inner radius r t and r s , respectively, as shown in Figure 1 which indicates the solutes transporting through the membrane by diffusion for mass transfer. An ultrafiltration flow with flow rate V w and with inlet flow rates V a I and V b I of the retentate and dialysate phases, respectively, was introduced in the concentric tubular dialyzer instead of omitting the ultrafiltration operation, leading to a dialysis rate improvement with channel thickness ratio and membrane sieving coefficient as parameters, as shown in Figure 1. The theoretical analysis is based on the following assumptions: (a) constant physical properties of fluid; (b) a fully developed laminar flow on the entire length in the inner and annular subchannels; (c) negligible entrance length and end effects; (d) negligible longitudinal diffusion; (e) constant membrane sieving coefficient θ .

2.1.1. Velocity Distributions

The velocity distributions were derived with the use of the continuity equations and momentum balance equations in both inner subchannel a and annular subchannel b under the assumed conditions as follows:
u a r r + u a r r + u a z z = 0
u b r r + u b r r + u b z z = 0
u a z u a z z + u a r u a z r = 1 ρ p a z + υ 2 u a z r 2 + u a z r r + 2 u a z z 2
u a z u a r z + u a r u a r r = 1 ρ p a r + υ 2 u a r r 2 + u a r r u a r r 2 + 2 u a r z 2
u b z u b z z + u b r u b z r = 1 ρ p b z + υ 2 u b z r 2 + u b z r r + 2 u b z z 2
u b z u b r z + u b r u b r r = 1 ρ p b r + υ 2 u b r r 2 + u b r r u b r r 2 + 2 u b r z 2
The dimensionless forms of the velocity distributions, u a ξ and u a η , are the velocity profiles in the longitudinal and perpendicular directions in the inner subchannel incorporating with the dimensionless groups, = r t r s , η = r r s and ξ = z L , u a I ¯ = V a I π r t 2 , v w 0 = V w 2 π r t L , λ a = Δ r s v w 0 υ , which were derived in Appendix A as follows:
u a ξ = [ 2 Δ 2 u a I ¯ 4 Δ v w 0 L r s ξ ] Δ 2 η 2 Δ 4 + λ a 1 18 Δ 2 η 2 4 Δ 4 + η 4 4 Δ 6 η 6 18 Δ 8 + λ a 2 83 5400 Δ 2 19 ( η 2 ) 270 Δ 4 + 11 ( η 4 ) 144 Δ 6 η 6 36 Δ 8 + η 8 144 Δ 10 η 10 1800 Δ 12
u a η = 2 Δ v w 0 1 Δ 4 Δ 2 η η 3 2 + λ a η 18 Δ 2 η 3 8 Δ 4 + η 5 12 Δ 6 η 7 72 Δ 8 + λ a 2 83 η 5400 Δ 2 19 η 3 540 Δ 4 + 11 η 5 432 Δ 6 η 7 144 Δ 8 + + η 9 720 Δ 10 η 11 10,800 Δ 12
Similarly, the dimensionless forms of velocity distributions in the annulus subchannel, u b ξ and u b η , were obtained using the same procedure except for the stream function, and were derived in Appendix B with 1 = r o r s , u b I ¯ = V b I π r s 2 r o 2 , v w 1 = V w 2 π r o L , λ b = Δ 1 r s v w 1 υ , as follows:
u b ξ = 2 1 Δ 1 2 u b I ¯ + 4 Δ 1 v w 1 L r s ξ f b 0 ( η ) + λ b f b 1 ( η )
u b η = 2 Δ 1 v w 1 η f b 0 ( η ) + λ b f b 1 ( η )

2.1.2. Mass Balance Equations

The mass transfer equations in the retentate, dialysate and membrane phases were derived in the concentric tubular dialyzer with the ultrafiltration system, as shown in Figure 2
u a ξ L ϕ a ξ + u a η r s ϕ a η = D r s 2 2 ϕ a η 2 + 1 η ϕ a η
u b ξ L ϕ b ξ + u b η r s ϕ b η = D r s 2 2 ϕ b η 2 + 1 η ϕ b η
v w ( η ) r s ϕ m η = D m r s 2 2 ϕ m η 2 + 1 η ϕ m η
in which ϕ i = C i C a I , and the subscripts a, b and m refer to the retentate, dialysate and membrane phases, respectively; v w η = 1 η V w 2 π r L , Δ η Δ 1 ; δ is the membrane thickness; D is the solute diffusivity in the retentate and dialysate phases; ε is the porosity of the membrane. The corresponding boundary conditions are
ϕ a ( η , 0 ) = 1 ,   0 η Δ
ϕ a ( 0 , ξ ) η = 0 ,   0 ξ 1
θ ϕ a ( Δ , ξ ) = ϕ m ( Δ , ξ ) ,   0 ξ 1
D ϕ a ( Δ , ξ ) η + v w 0 ϕ a w ( Δ , ξ ) = D m ϕ m ( Δ , ξ ) η + v w 0 ϕ m ( Δ , ξ ) ,   0 ξ 1
ϕ b ( Δ 1 , ξ ) = ϕ m ( Δ 1 , ξ ) ,   0 ξ 1
D ϕ b ( Δ 1 , ξ ) η = D m ϕ m ( Δ 1 , ξ ) η ,   0 ξ 1
ϕ b ( 1 , ξ ) η = 0 ,   0 ξ 1
ϕ b ( η , 0 ) = ϕ b I ,   Δ 1 η 1
where θ is the membrane sieving coefficient in which θ = 1 means the solute can freely pass through the membrane, and θ means that the solutes are partially rejected by the membrane.
The solute in the membrane phase is simultaneously transported by convection and diffusion due to the dialysis-and-ultrafiltration operation, as shown in Figure 2 and indicated in Equation (13), which can be rewritten as
2 ϕ m t 2 + 1 η 1 V w 2 π L D m ϕ m η = 2 ϕ m t 2 α ϕ m t = 0
in which α = V w / 2 π L D m and t = ln η .
The general solution of Equation (22) is
ϕ m η , ξ = a 1 ξ exp α t + a 2 ξ = a 1 ( ξ ) η α + a 2 ( ξ )
One obtains Equations (24) and (25) for a 1 and a 2 , respectively, once Equation (23) is substituted into the boundary conditions of Equations (16) and (18):
a 1 ( ξ ) = θ ϕ a w ( ξ ) ϕ b w ( ξ ) Δ α Δ 1 α
and
a 2 ( ξ ) = Δ α ϕ b w ( ξ ) θ Δ 1 α ϕ a w ( ξ ) Δ α Δ 1 α
Thus, the solute concentration distribution ϕ m η , ξ in the membrane of Equation (23) and the derivative of ϕ m η , ξ were obtained as
ϕ m = θ ( η α Δ 1 α ) ϕ a w ( ξ ) + ( Δ α η α ) ϕ b w ( ξ ) Δ α Δ 1 α
ϕ m η = α θ ϕ a w ( ξ ) ϕ b w ( ξ ) Δ α Δ 1 α η α 1
By substituting Equation (27) into the boundary conditions of Equations (17) and (19), the derivatives of ϕ a and ϕ b were obtained in Equations (28) and (29), respectively:
ϕ a η = α ε θ ϕ a w ( ξ ) ϕ b w ( ξ ) Δ α Δ 1 α Δ α 1 r s v w 0 D θ 1 ϕ a w ( ξ ) ,   0 ξ 1
ϕ b η = α ε θ ϕ a w ( ξ ) ϕ b w ( ξ ) Δ α Δ 1 α Δ 1 α 1 ,   0 ξ 1
Hence, the two-dimensional solute concentration distributions of a concentric tubular dialysis-and-ultrafiltration operation can be solved by the governing equations, Equations (11) and (12), with the use of the boundary conditions of Equations (14), (15), (20), (21), (28) and (29).

2.2. Pure Membrane Dialysis Operation in a Concentric Tubular Module

A concentric tubular dialysis system without ultrafiltration operations, that is, V w = 0 , is a pure membrane dialysis without ultrafiltration operation, as shown in Figure 3. The solute is transported by diffusion only in the membrane phase.
The velocity distributions are the same equations, that is, Equations (7) and (8) and Equations (9) and (10), for the inner subchannel a (retentate phase) and annular subchannel b (dialysate phase), respectively, except for incorporating V w = 0 . The governing equations of the pure dialysis system in the retentate, dialysate and membrane phases are
u a ξ L ϕ a ξ = D r s 2 2 ϕ a η 2 + 1 η ϕ a η
u b ξ L ϕ b ξ = D r s 2 2 ϕ b η 2 + 1 η ϕ b η
2 ϕ m η 2 + 1 η ϕ m η = 0
and the corresponding boundary conditions are
ϕ a ( η ,   0 ) = 1 ,   0 η
ϕ a ( 0 , ξ ) η a = 0 ,   0 ξ 1
θ ϕ a ( 1 , ξ ) = ϕ m ( 0 , ξ ) ,   0 ξ 1
D ϕ a ( , ξ ) η = D m ϕ m ( , ξ ) η ,   0 ξ 1
ϕ b ( 1 , ξ ) = ϕ m ( 1 , ξ ) ,   0 ξ 1
D ϕ b ( 1 , ξ ) η = D m ϕ m ( 1 , ξ ) η ,   0 ξ 1
ϕ b ( 1 , ξ ) η = 0 ,   0 ξ 1
ϕ b ( η ,   0 ) = ϕ bI ,   1 η 1
The solute in membrane phase is transported by diffusion only, as shown in Figure 3 and indicated in Equation (32), which is called the Euler equation and can be solved easily by incorporating t = ln η to transform as follows:
2 ϕ m t 2 = 0
Integrating Equation (41) twice to obtain the general solution as
ϕ m ( η , ξ ) = a 1 t + a 2 φ m ( η , ξ ) = a 1 ln η + a 2
Thus, the solute concentration distribution ϕ m η , ξ in the membrane of Equation (42) and the derivative of ϕ m η , ξ were obtained with the use of the boundary conditions of Equations (35) and (37), as
ϕ m ( η , ξ ) = ( ln η ln Δ 1 ) θ ϕ a w ( ξ ) + ( ln Δ ln η ) ϕ b w ( ξ ) ln Δ ln Δ 1
ϕ m ( η , ξ ) η = 1 η θ ϕ a w ( ξ ) ϕ b w ( ξ ) ln Δ ln Δ 1
By substituting Equation (44) into the boundary conditions of Equations (36) and (38), the derivatives of ϕ a and ϕ b were obtained in Equations (45) and (46), respectively:
ϕ m ( η , ξ ) η = 1 η θ ϕ a w ( ξ ) ϕ b w ( ξ ) ln Δ ln Δ 1
ϕ m ( η , ξ ) η = 1 η θ ϕ a w ( ξ ) ϕ b w ( ξ ) ln Δ ln Δ 1
Hence, the two-dimensional solute concentration distributions of a concentric tubular dialysis-and-ultrafiltration operation can be solved by the governing equations, Equations (30) and (31), with the use of the boundary conditions of Equations (33), (34), (39), (40), (45) and (46).

3. Numerical Solutions for Solving Concentration Distributions

The mass balances of Equations (11) and (12) and Equations (30) and (31) were made for membrane dialysis along the flowing direction with and without ultrafiltration operations, respectively. Thus, the solute concentrations in both the inner and annular streams were solved numerically using the Crank–Nicolson method. The availability of computing software facilitates the numerical solution for this problem, in which the node numbers in the η -direction were N and J for the retentate and dialysate phases, respectively, and the node numbers in the ξ -direction were P for both the retentate and dialysate phases, when the iterative procedures reached the convergence tolerance. Hence, dialysis efficiency and dialysis rate improvement were obtained. Comparisons were made for the solute of dialysis efficiency with and without ultrafiltration operations. The partial differential equations have been discretized and transformed with routines incorporating the finite difference algorithm, which are expressed in Appendix C.

4. Dialysis Rate, Dialysis Efficiency and Dialysis Rate Improvement

The dialysis rate of the dialysis system with and without the ultrafiltration operation is defined as
M = V a I C a I V a o C a o
where C a o is the average outlet concentration of the retentate phase and V a o is the net outlet flow rate of the retentate phase ( V a o = V a I V w ). The dialysis efficiency of the dialysis system with the ultrafiltration operation is defined as the ratio of the dialysis rate to the maximum solute concentration gradient at the entrance as follows:
χ = M V aI ( C aI C bI ) = 1 ( 1 V w / V aI ) ϕ ao 1 ϕ bI
Furthermore, the dialysis rate improvement I D ( % ) of the dialysis system with the ultrafiltration operation was defined as the percentage increase in the dialysis efficiency χ relative to the dialysis efficiency of the pure dialysis system without the ultrafiltration operation:
I D ( % ) = χ χ 0 χ 0 × 100 %
where the χ 0 is the dialysis efficiency of the pure dialysis system without the ultrafiltration operation.

5. Results and Discussions

5.1. The Numerical Solutions Validated by Convergence Tolerance

The velocity distributions of u a ξ and u b ξ in the concentric tubular dialysis-and-ultrafiltration operation with the first- and second-order perturbation methods were obtained via comparison of the convergence tolerance, as demonstrated in Table 1 and Table 2, respectively. The selected order of perturbation methods indicated in Table 1 and Table 2 are dominant in this system. It was observed that only the first-order perturbation and zero-order perturbation in the inner subchannel and annular subchannel, respectively, need to be considered due to the convergence tolerance.
The two-dimensional theoretical model of the concentric tubular dialyzer with and without ultrafiltration operations was solved numerically using the Crank–Nicolson method, for which the convergence tolerance of the numerical solutions is shown in Table 3. The node numbers in the η -direction are N and J for the retentate and dialysate phases, respectively, and the node numbers in the ξ -direction are P for both the retentate and dialysate phases. The calculated results shown that two step sizes in the η -direction ( N and J ) and ξ -direction ( P ) agree reasonably well in the concentric tubular dialysis-and-ultrafiltration operation.

5.2. Concentration Distributions

The influences of the ultrafiltration rate V w on the axial velocities in the radial direction of both retentate and dialysate phases at the half axial position ξ = 0.5 with = 0.5 , are shown in Figure 4; both radial velocities vary with V w as a parameter, as indicated in Figure 5. The axial velocity in the retentate phase decreases with V w but the axial velocity in the dialysate phase increases with V w . Meanwhile, both radial velocities in the retentate and dialysate phases increase with V w . The findings shown in Figure 4 had previously been confirmed by Berman [33].
A smaller volumetric flow rate of retentate phase results in a lower concentration distribution due to the increase in the residence time for transporting the solute from the retentate phase to the dialysate phase, as shown in Figure 6. Comparisons were made regarding concentration distributions between the dialysis system with ( V w = 1 mL/min) and without ( V w = 0   m L / m i n , as refers to the pure dialysis system without ultrafiltration) ultrafiltration rates in Figure 7, in which the results demonstrate the concentration distribution of the dialysis-and-ultrafiltration system decreases along the axial direction in the retentate phase, but increases in the dialysate phase.
The influence of ultrafiltration rates and the retentate phase flow rates on the average concentrations of both the retentate and dialysate phases are presented in Figure 8 and Figure 9, respectively. The solute has a longer time to transport from the retentate phase to the dialysate phase with decreasing retentate phase flow rate due to increasing residence time under the lower retentate phase flow rate. Therefore, the average concentration of the retentate phase increases with increasing retentate phase flow rate but with decreasing ultrafiltration rate, as shown in Figure 8. The results indicate that the average concentration of the retentate phase may not decrease with increasing ultrafiltration operation owing to both the solute and solvent transporting through the membrane progressively and simultaneously. Average concentration distributions of the dialysis system in the retentate phase decrease along the axial direction and with the ultrafiltration rate, as shown in Figure 8. Figure 9 illustrates that the average concentration of the dialysate phase increases with increasing ultrafiltration rate and retentate phase flow rate. However, average concentration distributions in the dialysate phase increase along the axial direction as well as with the ultrafiltration rate except for the turning point by the inverse order with the ultrafiltration rate within 0.4 < ξ < 0.6 , as shown in Figure 9. Furthermore, average concentration distributions in both the retentate and dialysate phases increase with the volumetric flow rate of the dialysate phase accordingly.
Similarly, the influences of various ultrafiltration rates and volumetric flow rates of the dialysate phase on both the retentate and dialysate phases are presented in Figure 10 and Figure 11, respectively. Average concentration distributions in the dialysate phase increase along the axial direction as well as with the ultrafiltration rate except for the turning point by the inverse order with the ultrafiltration rate shifting down to 0. 2 < ξ < 0.4 , as shown in Figure 12, and both average concentration distributions in the retentate and dialysate phases decrease with the volumetric flow rate of the dialysate phase. The effects of the membrane sieving coefficient θ and channel thickness ratio on the average concentration distribution of the retentate phase were examined with regard to the dialysis efficiency. Figure 12 shows the average concentration distribution of the retentate phase decreases with increases in both the membrane sieving coefficient and channel thickness ratio.
The average outlet concentrations of the retentate phase ϕ a , o ¯ were studied with various ultrafiltration rates, channel thickness ratios and volumetric flow rates of the retentate phase as parameters, as indicated in Figure 13 and Figure 14. Figure 13 shows the average outlet concentration of retentate phase ϕ a , o ¯ increases with increasing retentate phase flow rate because the residence time is increased in the retentate phase. The results reveal that the average outlet concentrations of the retentate phase ϕ a , o ¯ decrease with increasing ultrafiltration rate and the volumetric flow rate of the dialysate phase, as indicated in Figure 13.
The effect of the channel thickness ratio on the average outlet concentration of the retentate phase ϕ a , o ¯ is demonstrated in Figure 14. The average outlet concentration of the retentate phase ϕ a , o ¯ decreases with decreasing channel thickness ratio, such as enlarging the annulus channel of the dialysate phase.

5.3. Dialysis Rate

The solute is transported through the membrane by two mechanisms in a dialysis module with ultrafiltration operation: (a) diffusion (caused by the concentration difference across the membrane) and (b) convection (caused by the ultrafiltration operation). The dialysis rate M of the concentric tubular dialysis-and-ultrafiltration system under ultrafiltration operation is calculated, as illustrated in Figure 15 and Figure 16.
The results show that the dialysis rate M increases with the retentate phase flow rate as well as the ultrafiltration rate in comparison to the system without ultrafiltration operation ( V w = 0 when a pure dialysis system is employed). Meanwhile, Figure 15 indicates that the dialysis rate M increases with increasing retentate phase flow rate, which contributes to the average fluid velocity in the retentate phase and results in a higher convective mass transfer coefficient. This is attributed to the concentration gradient between the retentate and dialysate phases owing to the increasing dialysate phase flow rate, which thus results in an increased dialysis rate due to diffusion. Moreover, Figure 15 implies that the dialysis rate M increases with both the retentate and dialysate phase flow rates, which results in an increase in the convective mass-transfer coefficient in both the retentate and dialysate phases. The membrane sieving coefficient plays an important role in increasing the mass transfer rate by diffusing through the membrane, and thus, the dialysis rate M increases with the membrane sieving coefficient as well as the ultrafiltration rate and volumetric flow rate of the dialysate phase, as demonstrated in Figure 16.
For the same reason as refers to the increase in the convective mass-transfer coefficient, the dialysis rate M also increases with increasing channel thickness ratio , as shown in Figure 17.
The experimental work is cost prohibitive at the present time and there is only the available experimental data from the flat-plate dialyzer for confirmation. The present study could be analogized from those in our previous work on the flat-plate module [35]; this is a forthcoming work for our laboratory and will be conducted in the near future. In order to validate the theoretical predictions of the two-dimensional mathematical formulations in the concentric circular membrane dialyzer, a comparison was made with the experimental results conducted by the flat-plate dialyzer under ultrafiltration operation in our previous work [35] which neglected the tube curvature effect in the concentric tube [36]. The diminutive size, that is, = 0.7 , of the annular spacing was relatively smaller than that of the tube diameters. Figure 17 denotes that the dialysis rate M also increases with the channel thickness ratio of both the concentric circular dialyzer (the present study) and flat-plate dialyzer [35] for comparisons. Restated, the geometry of this construction with negligible curvature effect can be approximated as that of a parallel-plate device with = 0.5 . The theoretical predictions and experimental results of the dialyzer rate M are shown in Figure 17 within acceptable experimental deviations.

5.4. Dialysis Efficiency and Dialysis Rate Improvement

The dialysis efficiency χ of the concentric tubular dialysis-and-ultrafiltration system under ultrafiltration operation is expressed in terms of the ratio of the dialysis rate to the maximum solute concentration gradient at the entrance, as illustrated in Figure 18 and Figure 19, respectively.
The dialysis efficiency χ increases significantly with the ultrafiltration rate, especially under the lower volumetric flow rate of the dialysate phase, but decreases slightly with the volumetric flow rate of the retentate phase for the dialysis system under ultrafiltration operation, as shown by Figure 18. Results indicate that the ultrafiltration rate has a significant effect on the clearance of the solute, which has already been demonstrated by Jagannathan and Shettigar [22]. Figure 19 shows that the dialysis efficiency χ also increases with the channel thickness ratio. In addition, Figure 18 and Figure 19 indicate that the dialysis efficiency χ increases with the dialysate phase flow rate but decreases with increasing volumetric flow rate of the retentate phase.
The influence of the membrane sieving coefficient on the dialysis efficiency χ is presented in Figure 20. The results show that the dialysis efficiency χ increases with the membrane sieving coefficient and ultrafiltration rate but decreases with increasing volumetric flow rate of the retentate phase, as confirmed in Figure 20.
The dialysis rate improvements I D % for θ = 1 and = 0.5 were calculated as the relative percentage increment in concentric tubular dialysis-and-ultrafiltration module as compared with the module without ultrafiltration operation. The dialysis rate improvements I D ( % ) are presented in Table 4 and Table 5, which show the results of the dialysis rate improvement of the concentric tubular dialysis-and-ultrafiltration module for various V w and , respectively. The dialysis rate improvements I D ( % ) increase with V w due to the increase in both diffusion and convective transfer through the membrane from the dialysate phase to the dialysate phase. The results also show that the dialysis rate improvements decrease with increasing as well as both V a and V b . However, the maximum dialysis rate improved up to twice as compared with that of the pure dialysis system ( V w = 0 ) by employing a dialysis system with ultrafiltration operation, as shown in Table 4.
The dialysis rate improvements I D % increase with ultrafiltration rate V w and membrane sieving coefficient θ but decrease with channel thickness ratio and both volumetric flow rates ( V a and V b ) of the dialysate and retentate phases. The concentric tubular dialyzer with ultrafiltration operations transports more solute flux through the membrane of the tubular dialyzer with a higher dialysis rate. Generally, operating with the ultrafiltration rate shows a substantial influence which enhances the dialysis rate in the concentric tubular dialyzer. Several results can be concluded from Table 4 and Table 5, as follows. (1) The dialysis rate improvement I D % increases with increasing ultrafiltration rate due to increasing convective mass transfer rate of the solute through the membrane; (2) the dialysis rate improvement I D % decreases with increases in both the retentate phase flow rate and dialysate phase flow rate. The solute is transported through the membrane by diffusion and convection simultaneously in a dialysis-and-ultrafiltration operation, which means the contribution of the ultrafiltration rate on the dialysis rate is weakened when both the retentate phase and dialysate phase flow rates are increased.

6. Conclusions

A concentric tubular dialyzer with ultrafiltration operation to augment the dialysis rate was investigated theoretically. Two-dimensional mathematical equations were developed and formulated to predict the dialysis efficiency and dialysis rate improvement in the dialysis-and-ultrafiltration system as compared with the module without the ultrafiltration operation, which was solved numerically using the Crank–Nicolson method. There were two innovation points of this module design: (1) a two-dimensional mathematical model of the concentric tubular dialyzer was developed theoretically; (2) the simulated results show that the dialysis rate improvement was significantly improved with implementation of the ultrafiltration effect. A more direct method provides a straightforward strategy to the solution for determining the effect of the channel thickness ratio, which allows the specification setting by the designer to be met with economic consideration.
Regarding the pure dialysis system without ultrafiltration operation, its dialysis efficiency could be readily enhanced by employing the module with ultrafiltration operation, especially for operations with lower volumetric flow rate of the retentate phase. A maximum dialysis rate improvement of up to twice was found in the module with the ultrafiltration rate V w = 2   m L / m i n and membrane sieving coefficient θ = 1 in comparisons with the module without ultrafiltration operation. The comparisons of the concentric tubular dialyzer with ultrafiltration operation relative to the module without ultrafiltration operations led to the following conclusions:
  • Average concentration distributions of the dialysis system in the retentate phase decrease along the axial direction and with decreasing ultrafiltration rate whereas the average concentration of the dialysate phase increases with increasing ultrafiltration rate and the retentate phase flow rate.
  • Both average concentration distributions in the retentate and dialysate phases decrease with the volumetric flow rate of the dialysate phase.
  • The average concentration distribution of the retentate phase decreases with increases in both the membrane sieving coefficient and channel thickness ratio.
  • The average outlet concentration of retentate phase increases with increasing retentate phase flow rate because the residence time is increased as well as with increasing ultrafiltration rate, channel thickness ratio and volumetric flow rate of the dialysate phase.
  • The results show that the dialysis rate M increases with increases in both retentate and dialysate phase flow rates, membrane sieving coefficient, ultrafiltration rate and channel thickness ratio relative to the system without ultrafiltration operation.
  • The dialysis efficiency χ increases with increasing dialysate phase flow rate, ultrafiltration rate, membrane sieving coefficient and channel thickness ratio but with decreasing volumetric flow rate of the retentate phase.
  • The dialysis rate improvements I D ( % ) increase with increasing ultrafiltration rate and membrane sieving coefficient but with decreasing channel thickness ratio and both retentate and dialysate phase flow rates.
The results demonstrate the technical feasibility of dialysis rate improvement in the tubular membrane dialyzer with ultrafiltration operation. It was also found that the concentrations of the components (minerals and contaminants) in the retentate could be removed to reduce solutes in wastewater treatment processes, which depend on the separation technique and the operational parameters in membrane-based separation processes. In this paper, both the dialysis rate and ultrafiltration operation were examined from an economic perspective by implementing various ultrafiltration rates in the tubular dialyzer. Therefore, the alternative membrane sieving coefficient, membrane material and ultrafiltration rate require further investigation regarding the economic considerations of the tubular dialyzer. It is believed that the availability of such a simplified mathematical formulation as developed here is the value in the present work in designing a concentric tubular dialyzer and will be an important contribution to the design and investigation of multi-stream or multi-phase problems with coupling mutual boundary conditions. One may follow the present theory and develop a mathematical formulation to deal with multi-stream or multi-phase heat- or mass-transfer devices for each particular application with various geometries.

Author Contributions

C.-D.H.: funding acquisition, conceptualization, methodology, data collection and analysis. J.-W.T.: formal analysis, data curation, investigation. J.-W.L.: investigation, writing—review, editing. W.-C.L.: writing—original draft, analysis and interpretation of data. 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 are contained within the article.

Acknowledgments

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

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

C Solute concentration in the stream (mol m−3)
D Solute diffusivity in both retentate and dialysate phases (m2 s−1)
D m Solute diffusivity in membrane (m2 s−1)
I D Dialysis rate improvement
J Node numbers in η -direction of dialysate phase
L Length of membrane dialyzer (m)
M Dialysis rate (mol s−1)
N Node numbers in η -direction of retentate phase
P Node numbers in ξ -direction
p ξ -direction step sizes of the Crank–Nicolson method
p a Pressure in retentate phase (N m−2)
p b Pressure in dialysate phase (N m−2)
r Transversal coordinate (m)
r o Outside radius of inner tube (m)
r s Inside radius of the shell tube (m)
r t Inside radius of inner tube (m)
R e Reynolds number
u a Velocity distribution of the retentate phase (m s−1)
u b Velocity distribution of the dialysate phase (m s−1)
u a I ¯ Averaged velocity of the retentate phase (m s−1)
u b I ¯ Averaged velocity of the dialysate phase (m s−1)
V a I Inlet volumetric flow rate of the retentate phase (m3 s−1)
V a o Outlet volumetric flow rate of the retentate phase (m3 s−1)
V b I Inlet volumetric flow rate of the dialysate phase (m3 s−1)
V b o Outlet volumetric flow rate of the dialysate phase (m3 s−1)
V w Ultrafiltration rate (m3 s−1)
u w Velocity distribution of the dialysate phase (m s−1)
v w Ultrafiltration flux distribution (m s−1)
v w 0 Average ultrafiltration flux on the inner surface of the membrane (m s−1)
v w 1 Average ultrafiltration flux on the outer surface of the membrane (m s−1)
z Axial coordinate along the flow direction (m)
Greek letters
Ratio of inside radius of inner tube to shell radius, r t / r s
1 Ratio of outside radius of inner tube to shell radius, r o / r s
δ Membrane thickness (m)
ε Membrane porosity
η Dimensionless transversal coordinate
θ Membrane sieving coefficient
λ a Wall Reynolds number of the retentate phase, λ a = Δ r s v w 0 / υ
λ b Wall Reynolds number of the dialysate phase, λ b = Δ 1 r s v w 1 / υ
μ Wall Reynolds number of the dialysate phase ( kg   m 1 s 1 )
ξ Dimensionless longitudinal coordinate
l ρ Density ( k g m 3 )
l υ Kinetic viscosity ( m 2 s 1 )
ϕ Dimensionless solute concentration
ϕ Dimensionless average solute concentration
χ Dialysis efficiency
φ Stream function
Subscripts
a In the retentate phase
b In the dialysate phase
I At the inlet
o At the outlet
w On the membrane surface

Appendix A

The corresponding boundary equations of Equations (3)–(6) are
u a z ( 0 , z ) r = 0 ,   0 z L
u a z ( r t , z ) = 0 ,   0 z L
u a z ( r , 0 ) = 2 u a I ¯ ( 1 r 2 r t 2 ) ,   0 r r t
p a + μ u a z ( r , 1 ) z = 0 ,   0 r r t
u a r ( r t , z ) = u b r ( r o , z ) ,   0 z L
u a r ( 0 , z ) = 0 ,   0 z L
u a r ( r , 0 ) = 0 ,   0 r r t
u a r ( r , 1 ) z = 0 ,   0 r r t
u b z ( r o , z ) = 0 ,   0 z L
u b z ( r s , z ) = 0 ,   0 z L
u b z ( r , 0 ) = 2 u b I ¯ r s 2 r 2 + r s 2 r o 2 ln r ln r s ln r s ln r o r s 4 r o 4 r s 2 r o 2 r s 2 r o 2 ln r s ln r o ,   r o r r s
p b + μ u b z ( r ,   1 ) z = 0 ,   r o r r s
u b r ( r s , z ) = 0 ,   0 z L
u b r ( r , 0 ) = 0 ,   r o r r s
u b r ( r , 1 ) z = 0 ,   r o r r s
in which u a I ¯ = V a I / π r t 2 and u b I ¯ = V b I / π ( r s 2 r o 2 ) are the average velocities for both inner subchannel a and annular subchannel b, respectively. Meanwhile, the average flow rate through the membrane was calculated as
v w ( r ) = V w 2 π r L ,   r t r r o
Then, the average ultrafiltration flux on the inner and outer surfaces of the membrane were evaluated accordingly as v w 0 = V w / 2 π r t L and v w 1 = V w / 2 π r o L , respectively. Thus, Equation (A5) can be rewritten as follows:
u a r ( r t , z ) = v w 0 = constant ,   0 z L
u b r ( r o , z ) = v w 1 = constant ,   0 z L
The stream function [25] of the inner subchannel was simplified by using the reservation of mass in the inner subchannel in terms of the dimensionless group η * = r / r s 2 as
φ a ( η * , z ) = [ r t 2 u a I ¯ 2 r t v w 0 z ] f a ( η * )
The derivatives of Equation (A19) with respect to r and z , respectively, were obtained to give
φ a r = r u a z   or   u a z = 2 Δ 2 u a I ¯ 4 Δ v w 0 r s z f a ( η * )
φ a z = r u a r   or   u a r = 2 Δ v w 0 η * f a ( η * )
in which = r t / r s . Substituting Equations (A20) and (A21) into Equations (3) and (4), Equations (3) and (4) can be rewritten as
4 r s 2 2 Δ 2 u a I ¯ 4 Δ v w 0 r s z Δ v w 0 r s f a f a ( f a ) 2 υ η * f a + f a = 1 ρ p a z
and
2 Δ v w 0 r s Δ v w 0 r s 2 f a f a η * f a 2 ( η * ) 2 2 υ f a = 1 ρ p a η *
Combinations of the derivative of Equation (A22) with respect to r and the derivative of Equation (A23) with respect to z will obtain
η * Δ v w 0 r s f a f a ( f a ) 2 υ η * f a + f a = 0
or
η * f a + f a λ a f a f a ( f a ) 2 = k a
in which λ a = Δ r s v w 0 / υ is the wall Reynolds number and k a is the integration constant. When the stream function of Equation (A19) was introduced into the boundary conditions, Equations (A1), (A2), (A6) and (A17) can be rewritten as
f a ( 0 ) = 0
f a ( Δ 2 ) = 0
f a ( 0 ) = 0
f a ( Δ 2 ) = 1 / 2
The perturbation method [28] was used to solve the expressions of both k a and f a when λ a 1 as follows:
f a = f a 0 + λ a f a 1 + λ a 2 f a 2 + . . . + λ a n f a n + . . .
k a = S a 0 + λ a S a 1 + λ a 2 S a 2 + . . . + λ a n S a n + . . .
Applying the perturbation method to Equations (A30) and (A31) with the use of boundary conditions yields f a and the derivative f a once all constants ( f a i and k a i , i = 0 , 1 , 2 , ) were obtained as follows:
f a = 1 Δ 4 Δ 2 η * 1 2 ( η * ) 2 + λ a η * 18 Δ 2 ( η * ) 2 8 Δ 4 + ( η * ) 3 12 Δ 6 ( η * ) 4 72 Δ 8 + λ a 2 83 η * 5400 Δ 2 19 ( η * ) 2 540 Δ 4 + 11 ( η * ) 3 432 Δ 6 ( η * ) 4 144 Δ 8 + ( η * ) 5 720 Δ 10 ( η * ) 6 10,800 Δ 12
f a = Δ 2 η * Δ 4 + λ a 1 18 Δ 2 η * 4 Δ 4 + ( η * ) 2 4 Δ 6 ( η * ) 3 18 Δ 8 + λ a 2 83 5400 Δ 2 19 ( η * ) 270 Δ 4 + 11 ( η * ) 2 144 Δ 6 ( η * ) 3 36 Δ 8 + ( η * ) 4 144 Δ 10 ( η * ) 5 1800 Δ 12
The dimensionless forms of velocity distributions, u a ξ and u a η , of Equations (7) and (8) are obtained after substituting Equations (A32) and (A33) into Equations (A20) and (A21) accordingly.

Appendix B

The stream function of the annulus subchannel was simplified by using the reservation of mass as
φ b ( η * , z ) = [ ( r s 2 r o 2 ) u b I ¯ + 2 r o v w 1 z ] f b ( η * )
The derivatives of Equation (A34) with respect to r and z , respectively, were obtained to give
φ b r = r u b z   or   u b z = 2 ( 1 Δ 1 ) 2 u b I ¯ 4 Δ 1 v w 1 r s z f b ( η * )
φ b z = r u b r   or   u b r = 2 Δ 1 v w 1 η * f b ( η * )
in which f b ( η ) = f b 0 ( η ) + λ b f b 1 ( η ) and
f b 0 η = S b 0 1 2 + W 1 + η 4 2 η 2 + W 1 η 2 ln η 2 η 2
and
f b 1 η = S b 1 2 η 4 + I 0 η 2 l n η 2 η 2 + I 1 η 2 + I 2 + S b 0 2 [ η 8 72 η 6 72 6 + 5 W 1 ] W 1 η 2 2 1 + 2 W 1 + W 1 12 η 2 l n η 2 ( 6 12 η 2 + η 4 + 12 W 1 36 W 1 η 2 ) + W 1 4 η 2 ( l n η 2 ) 2 ( 1 2 W 1 + 2 W 1 η 2 ) + W 1 η 2 4 ( 1 + 10 W 1 + 24 W 1 2 )
in which I 0 , I 1 and I 2 are integration constants, and S b 0 and S b 1 are obtained as follows:
S b 0 = 1 l n 1 2 ( 1 1 4 ) + 2 ( 1 1 2 ) 2
S b 1 = S b 0 2 11 36 + W 1 + 9 W 1 2 + S b 0 2 W 2 9 W 1 2 36 4 1 2 39 + 11 1 6 36 1 8 8 + W 1 ( 1 2 8 1 4 4 + 1 6 8 ) + l n 1 2 S b 0 2 W 2 1 36 + 4 1 4 72 1 6 6 + 1 8 24 W 1 1 144 + 3 1 2 2 + 29 1 4 16 + 11 1 6 36 + W 1 2 ( 1 2 4 1 2 + 7 1 4 2 ) + ( l n 1 2 ) 2 S b 0 2 W 2 W 1 1 4 1 2 4 1 4 + 1 6 6 + W 1 2 ( 1 2 + 1 2 2 + ( l n 1 2 2 3 ) 1 4
where W 1 = 1 1 2 l n 1 2 and W 2 = 1 2 1 2 + 1 4 + l n 1 2 2 ( 1 1 4 )

Appendix C

The partial differential equations have been discretized and transformed with routines incorporating the finite difference algorithm, as shown in Figure A1, as follows:
ϕ ξ = ϕ p + 1 n ϕ p n s
ϕ η = 1 2 ϕ p + 1 n + 1 ϕ p + 1 n 1 2 k + ϕ p n + 1 ϕ p n 1 2 k
2 ϕ η 2 = 1 2 ϕ p + 1 n + 1 2 ϕ p + 1 n + ϕ p + 1 n 1 k 2 + ϕ p n + 1 2 ϕ p n + ϕ p n 1 k 2
Figure A1. The scheme meshes of the Crank–Nicolson method. ( s and k present the step sizes of ξ -direction and η -direction, respectively).
Figure A1. The scheme meshes of the Crank–Nicolson method. ( s and k present the step sizes of ξ -direction and η -direction, respectively).
Membranes 13 00556 g0a1

References

  1. Chen, L.; Wu, Y.L.; Dong, H.J.; Meng, M.J.; Li, C.X.; Yan, Y.S.; Chen, J. An overview on membrane strategies for rare earths extraction and separation. Sep. Purifi. Technol. 2018, 197, 70–85. [Google Scholar] [CrossRef]
  2. Overmans, S.; Ignacz, G.; Beke, A.K.; Xu, J.J.; Saikaly, P.E.; Szekely, G.; Lauersen, K.J. Continuous extraction and concentration of secreted metabolites from engineered microbes using membrane technology. Green Chem. 2022, 24, 5479–5489. [Google Scholar] [CrossRef]
  3. Zhang, Z.; Chen, F.; Rezakazemi, M.; Zhang, W.X.; Lu, C.F.; Chang, H.X.; Quan, X.J. Modeling of a CO2-piperazine-membrane absorption system. Chem. Eng. Res. Des. 2018, 131, 375–384. [Google Scholar] [CrossRef]
  4. Ochedi, F.O.; Yu, J.L.; Yu, H.; Liu, Y.X.; Hussain, A. Carbon dioxide capture using liquid absorption methods: A review. Environ. Chem. Lett. 2021, 19, 77–109. [Google Scholar] [CrossRef]
  5. Ran, J.; Wu, L.; He, Y.B.; Yang, Z.J.; Wang, Y.M.; Jiang, C.X.; Ge, L.; Bakangura, E.; Xu, T.W. Ion exchange membranes: New developments and applications. J. Membr. Sci. 2017, 522, 267–291. [Google Scholar] [CrossRef]
  6. Epsztein, R.; Shaulsky, E.; Qin, M.H.; Elimelech, M. Activation behavior for ion permeation in ion-exchange membranes: Role of ion dehydration in selective transport. J. Membr. Sci. 2019, 580, 316–326. [Google Scholar] [CrossRef]
  7. Wang, P.; Chung, T.S. Recent advances in membrane distillation processes: Membrane development, configuration design and application exploring. J. Membr. Sci. 2015, 474, 39–56. [Google Scholar] [CrossRef]
  8. Thomas, N.; Mavukkandy, M.O.; Loutatidou, S.; Arafat, H.A. Membrane distillation research & implementation: Lessons from the past five decades. Sep. Purif. Technol. 2017, 189, 108–127. [Google Scholar]
  9. Alsarayreh, A.A.; Al-Obaidi, M.A.; Patel, R.; Mujtaba, I.M. Scope and limitations of modelling, simulation, and optimisation of a spiral wound reverse osmosis process-based water desalination. Processes 2020, 8, 573. [Google Scholar] [CrossRef]
  10. Ezugbe, E.O.; Rathilal, S. Membrane Technologies in Wastewater Treatment: A Review. Membranes 2020, 10, 89. [Google Scholar] [CrossRef]
  11. Asif, M.B.; Zhang, Z. Ceramic membrane technology for water and wastewater treatment: A critical review of performance, full-scale applications, membrane fouling and prospects. Chem. Eng. J. 2021, 418, 129481. [Google Scholar] [CrossRef]
  12. Lake, M.A.; Melsheimer, N.J. Mass transfer characterization of Donnan dialysis. AIChE J. 1978, 24, 130–137. [Google Scholar] [CrossRef]
  13. Moonen, H.; Niefind, N.J. Alcohol reduction in beer by means of dialysis. Desalination 1982, 41, 327–335. [Google Scholar] [CrossRef]
  14. Klein, E.; Holland, E.F.; Lebeouf, A.; Donnand, A.; Smith, J.K. Transport and mechanical properties of hemodialysis hollow fibers. J. Membr. Sci. 1976, 1, 371–396. [Google Scholar] [CrossRef]
  15. Ho, W.S.W.; Sirkar, K.K. Handbook of Membranes; Van Nostrand Reinhold: New York, NY, USA, 1992. [Google Scholar]
  16. Meng, F.; Seredych, M.; Chen, C.; Gura, V.; Mikhalovsky, S.; Sandeman, S.; Ingavle, G.; Ozulumba, T.; Miao, L.; Anasori, B.; et al. Mxene sorbents for removal of urea from dialysate: A step toward the Wearable Artificial Kidney. ACS Nano 2018, 12, 10518–10528. [Google Scholar] [CrossRef]
  17. Colton, C.K.; Smith, K.A.; Stroeve, P.; Merrill, E.W. Laminar flow mass transfer in a flat duct with permeable walls. AIChE J. 1971, 17, 773–780. [Google Scholar] [CrossRef]
  18. Walker, G.; Davis, T. Mass transfer in laminar flow between parallel permeable plates. AIChE J. 1974, 20, 881–889. [Google Scholar] [CrossRef]
  19. Cooney, D.O.; Kim, S.S.; Davis, E.J. Analyses of mass transfer in hemodialyzers for laminar blood flow and homogeneous dialysate. Chem. Eng. Sci. 1974, 29, 1731–1738. [Google Scholar] [CrossRef]
  20. Özdural, A.R.; Alkan, A. Derivation of a new explicit equation for the determination of overall mass transfer coefficients in continuous dialyzers. J. Membr. Sci. 2003, 223, 49–57. [Google Scholar] [CrossRef]
  21. Popovich, R.; Christopher, G.; Babb, A.L. The effect of membrane diffusion and ultrafiltration properties on hemodialyzer design and performance. Chem. Eng. Prog. Symp. Ser. 1971, 67, 105–115. [Google Scholar]
  22. Jagannathan, R.; Shettigar, U.R. Analysis of a tubular hemodialyser—Effect of ultrafiltration and dialysate concentration. Med. Biol Eng. Comput. 1977, 15, 500–512. [Google Scholar] [CrossRef]
  23. Henderson, L.W. Current status of hemofiltration. Artif. Organs. 1978, 2, 120–124. [Google Scholar] [CrossRef]
  24. Leber, H.W.; Wizemann, V.; Goubeaud, G.; Rawer, P.; Schutterle, G. Hemodiafiltration: A new alternative to hemofiltration and conventional hemodialysis. Artif. Organs. 1978, 2, 150–153. [Google Scholar] [CrossRef]
  25. Yeh, H.M.; Chen, T.W.; Chen, Y.J. Mass transfer for dialysis with ultrafiltration flux declined in cross-flow membrane modules. J. Chem. Eng. Jpn. 2000, 33, 440–448. [Google Scholar] [CrossRef]
  26. Tu, J.W.; Ho, C.D.; Yeh, H.M. The analytical and experimental studies of the parallel-plate concurrent dialysis system coupled with ultrafiltration. J. Membr. Sci. 2006, 281, 676–684. [Google Scholar] [CrossRef]
  27. Tu, J.W.; Ho, C.D.; Chuang, C.J. Effect of ultrafiltration on the mass-transfer efficiency improvement in a parallel-plate countercurrent dialysis system. Desalination 2009, 242, 70–83. [Google Scholar] [CrossRef]
  28. Lin, J.Y.; Huang, J.M.; Wang, J.; Yu, J.W.; You, X.Q.; Lin, X.C.; der Bruggen, B.V.; Zhao, S.F. High-performance porous anion exchange membranes for efficient acid recovery from acidic wastewater by diffusion dialysis. J. Membr. Sci. 2021, 624, 119116–119126. [Google Scholar] [CrossRef]
  29. Wang, J.X.; Liu, Z.Y.; Qiu, M.; He, C.J. Heparin-mimicking semi-interpenetrating composite membrane with multiple excellent performances for promising hemodialysis. J. Membr. Sci. 2021, 618, 118740–118750. [Google Scholar] [CrossRef]
  30. Liu, Y.; Li, G.L.; Han, Q.; Lin, H.B.; Li, Q.; Hu, J.; Liu, F. Anticoagulant dialyzer with enhanced Ca2+ chelation and hydrophilicity for heparin free hemodialysis. J. Membr. Sci. 2020, 604, 118082–118092. [Google Scholar] [CrossRef]
  31. Chen, Q.; He, Y.; Zhao, Y.P.; Chen, L. Intervening oxidative stress integrated with an excellent biocompatibility of hemodialysis membrane fabricated by nucleobase-recognized co-immobilization strategy of tannic acid, looped PEtOx brush and heparin. J. Membr. Sci. 2021, 625, 119174–119188. [Google Scholar] [CrossRef]
  32. Du, X.M.; Wang, Z.; Zhang, H.Y.; Yuan, Y.J.; Wang, H.; Zhang, Z.G. Prepared poly(aryl piperidinium) anion exchange membranes for acid recovery to improve dialysis coefficients and selectivity. J. Membr. Sci. 2021, 619, 11805–118812. [Google Scholar] [CrossRef]
  33. Berman, A.S. Laminar Flow in Channels with Porous Walls. J. Appl. Phys. 1953, 24, 1232. [Google Scholar] [CrossRef]
  34. Yuan, S.W.; Finkelstein, A.B.; Brooklyn, N.Y. Trans. ASME 1956, 78, 719–724. [Google Scholar]
  35. Tu, J.W.; Ho, C.D. Two-dimensional mass-transfer model of a flat-plate dialyzer with ultrafiltration operation. Chem. Eng. Technol. 2010, 33, 1358–1368. [Google Scholar] [CrossRef]
  36. Yeh, H.M.; Ward, H.C. The improvement in separation of concentric tube thermal diffusion columns. Chem. Eng. Sci. 1971, 26, 937–947. [Google Scholar] [CrossRef]
Figure 1. Concentric tubular dialysis-and-ultrafiltration module.
Figure 1. Concentric tubular dialysis-and-ultrafiltration module.
Membranes 13 00556 g001
Figure 2. Concentric tubular dialyzer module with ultrafiltration operation.
Figure 2. Concentric tubular dialyzer module with ultrafiltration operation.
Membranes 13 00556 g002
Figure 3. Concentric tubular dialyzer module without ultrafiltration operation.
Figure 3. Concentric tubular dialyzer module without ultrafiltration operation.
Membranes 13 00556 g003
Figure 4. Axial velocity distributions in both retentate and dialysate phases.
Figure 4. Axial velocity distributions in both retentate and dialysate phases.
Membranes 13 00556 g004
Figure 5. Radial velocity distributions in both retentate and dialysate phases.
Figure 5. Radial velocity distributions in both retentate and dialysate phases.
Membranes 13 00556 g005
Figure 6. Concentration distributions of both retentate and dialysate phases with V w and V a as parameters.
Figure 6. Concentration distributions of both retentate and dialysate phases with V w and V a as parameters.
Membranes 13 00556 g006
Figure 7. Concentration distributions of both phases with various longitudinal positions.
Figure 7. Concentration distributions of both phases with various longitudinal positions.
Membranes 13 00556 g007
Figure 8. Average concentration distribution of the retentate phase with various ultrafiltration rates and volumetric flow rates of the dialysate phase.
Figure 8. Average concentration distribution of the retentate phase with various ultrafiltration rates and volumetric flow rates of the dialysate phase.
Membranes 13 00556 g008
Figure 9. Average concentration distribution of the dialysate phase with various ultrafiltration rates and volumetric flow rates of the dialysate phase.
Figure 9. Average concentration distribution of the dialysate phase with various ultrafiltration rates and volumetric flow rates of the dialysate phase.
Membranes 13 00556 g009
Figure 10. Average concentration distribution of the retentate phase with various ultrafiltration rates and volumetric flow rates of the retentate phase.
Figure 10. Average concentration distribution of the retentate phase with various ultrafiltration rates and volumetric flow rates of the retentate phase.
Membranes 13 00556 g010
Figure 11. Average concentration distribution of the dialysate phase with various ultrafiltration rates and volumetric flow rates of the retentate phase.
Figure 11. Average concentration distribution of the dialysate phase with various ultrafiltration rates and volumetric flow rates of the retentate phase.
Membranes 13 00556 g011
Figure 12. Average concentration distribution of the retentate phase with the membrane sieving coefficient and channel thickness ratio as parameters.
Figure 12. Average concentration distribution of the retentate phase with the membrane sieving coefficient and channel thickness ratio as parameters.
Membranes 13 00556 g012
Figure 13. Average outlet concentration distribution of the retentate phase with the ultrafiltration rate and the retentate phase flow rate as parameters.
Figure 13. Average outlet concentration distribution of the retentate phase with the ultrafiltration rate and the retentate phase flow rate as parameters.
Membranes 13 00556 g013
Figure 14. Average outlet concentration distribution of the retentate phase with the channel thickness ratio and the retentate phase flow rate as parameters.
Figure 14. Average outlet concentration distribution of the retentate phase with the channel thickness ratio and the retentate phase flow rate as parameters.
Membranes 13 00556 g014
Figure 15. Dialysis rate with the ultrafiltration rate and the dialysate phase flow rate as parameters.
Figure 15. Dialysis rate with the ultrafiltration rate and the dialysate phase flow rate as parameters.
Membranes 13 00556 g015
Figure 16. Dialysis rate with the membrane sieving coefficient and ultrafiltration rate as parameters.
Figure 16. Dialysis rate with the membrane sieving coefficient and ultrafiltration rate as parameters.
Membranes 13 00556 g016
Figure 17. Comparisons of dialysis rate between theoretical results and experimental data [35] with the channel thickness ratio and the dialysate phase flow rate as parameters.
Figure 17. Comparisons of dialysis rate between theoretical results and experimental data [35] with the channel thickness ratio and the dialysate phase flow rate as parameters.
Membranes 13 00556 g017
Figure 18. Dialysis efficiency with the ultrafiltration rate and the dialysate phase flow rate as parameters.
Figure 18. Dialysis efficiency with the ultrafiltration rate and the dialysate phase flow rate as parameters.
Membranes 13 00556 g018
Figure 19. Dialysis efficiency with the channel thickness ratio and volumetric flow rate of the dialysate phase as parameters.
Figure 19. Dialysis efficiency with the channel thickness ratio and volumetric flow rate of the dialysate phase as parameters.
Membranes 13 00556 g019
Figure 20. Dialysis efficiency with the membrane sieving coefficient and ultrafiltration rate as parameters.
Figure 20. Dialysis efficiency with the membrane sieving coefficient and ultrafiltration rate as parameters.
Membranes 13 00556 g020
Table 1. Comparisons of the first- and second-order perturbation methods for u a ξ ( = 0.5 , ξ = 0.5 ,   V a = 20   m L / m i n ).
Table 1. Comparisons of the first- and second-order perturbation methods for u a ξ ( = 0.5 , ξ = 0.5 ,   V a = 20   m L / m i n ).
η V w = 1   m L / m i n V w = 2   m L / m i n
First OrderSecond OrderFirst OrderSecond Order
0.054.054.053.943.94
0.253.073.072.992.99
0.450.780.780.760.76
Table 2. Comparisons of the zero- and first-order perturbation methods for u b ξ ( = 0.5   , ξ = 0.5 , V b = 100   m L / m i n ).
Table 2. Comparisons of the zero- and first-order perturbation methods for u b ξ ( = 0.5   , ξ = 0.5 , V b = 100   m L / m i n ).
η V w = 1   m L / m i n V w = 2   m L / m i n
Zero OrderFirst OrderZero OrderFirst Order
0.551.161.161.171.17
0.755.465.465.485.48
0.951.821.821.831.83
Table 3. Comparisons of the convergence tolerance with the use of the Crank–Nicolson method ( = 0.5 , θ = 1 , V w = 2   m L / m i n , V a = 40   m L / m i n ,   V b = 300   m L / m i n ) .
Table 3. Comparisons of the convergence tolerance with the use of the Crank–Nicolson method ( = 0.5 , θ = 1 , V w = 2   m L / m i n , V a = 40   m L / m i n ,   V b = 300   m L / m i n ) .
Dimensionless
Concentration
P = 100P = 200
N = J = 400N = J = 500N = J = 400N = J = 500
ϕ a w 0.5430.5430.5430.543
ϕ a , o ¯ 0.9450.9450.9450.945
ϕ b w 0.2780.2780.2780.278
ϕ b , o ¯ 0.0170.0170.0170.017
Table 4. Dialysis rate improvement I D ( % ) of the concentric tubular dialysis-and-ultrafiltration operation ( = 0.5 and V b = 300   m L / m i n ) .
Table 4. Dialysis rate improvement I D ( % ) of the concentric tubular dialysis-and-ultrafiltration operation ( = 0.5 and V b = 300   m L / m i n ) .
V a
( m L / m i n )
θ = 1 θ = 0.6
V w = 1   m L / m i n V w = 2   m L / m i n V w = 1   m L / m i n V w = 2   m L / m i n
20103.55208.7987.2676.73
3084.34169.6463.3143.47
4072.69145.8348.9323.81
Table 5. Dialysis rate improvement I D ( % ) of the concentric tubular dialysis-and-ultrafiltration operation ( V a = 40   m L / m i n and V w = 2   m L / m i n ) .
Table 5. Dialysis rate improvement I D ( % ) of the concentric tubular dialysis-and-ultrafiltration operation ( V a = 40   m L / m i n and V w = 2   m L / m i n ) .
θ = 1 θ = 0.6
V b = 100   m L / m i n V b = 200   m L / m i n V b = 100   m L / m i n V b = 200   m L / m i n
0.3212.23176.21317.25274.24
0.5164.40145.83305.13278.61
0.7136.52127.00277.15261.17
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.; Tu, J.-W.; Lim, J.-W.; Lai, W.-C. Device Performance of a Tubular Membrane Dialyzer Incorporating Ultrafiltration Effects on the Dialysis Efficiency. Membranes 2023, 13, 556. https://doi.org/10.3390/membranes13060556

AMA Style

Ho C-D, Tu J-W, Lim J-W, Lai W-C. Device Performance of a Tubular Membrane Dialyzer Incorporating Ultrafiltration Effects on the Dialysis Efficiency. Membranes. 2023; 13(6):556. https://doi.org/10.3390/membranes13060556

Chicago/Turabian Style

Ho, Chii-Dong, Jr-Wei Tu, Jun-Wei Lim, and Wei-Chi Lai. 2023. "Device Performance of a Tubular Membrane Dialyzer Incorporating Ultrafiltration Effects on the Dialysis Efficiency" Membranes 13, no. 6: 556. https://doi.org/10.3390/membranes13060556

APA Style

Ho, C. -D., Tu, J. -W., Lim, J. -W., & Lai, W. -C. (2023). Device Performance of a Tubular Membrane Dialyzer Incorporating Ultrafiltration Effects on the Dialysis Efficiency. Membranes, 13(6), 556. https://doi.org/10.3390/membranes13060556

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