Performance optimization of polymeric porous membrane-based liquid desiccant air dehumidifier used in air conditioning system

In this study an experimental design was developed to optimize the performance and structure of a membrane-based parallel-plate liquid desiccant dehumidifier used in air conditioning regeneration system which operates under high humidity weather conditions. We conducted a series of polymeric porous membranes with different compositions fabricated that were prepared with various weight percentages of polysulfone (PSU), mixed with N-methyl-2-pyrrolidone (NMP) and dimethyl form amide (DMF) solvents. Furthermore, the designed experiments were performed under various operating conditions, indicating that the dehumidification efficiency declines with increasing flow rate, temperature, and humidity. Consequently, a membrane with optimized porosity and moisture permeability was selected which resulted in eliminating the carryover of solution droplets in the air, largely due to separating the flow condition of liquid desiccant (Li Cl) and air. This specific design is also greatly benefited by removing the water vapor from the air stream. The results of mathematical model simulations indicate that the DMF solvent had higher dehumidification capability compared with that of NMP under the optimized operating conditions. Additionally, it can clarify the porosity of the membrane which plays a significant role in the overall performance. Therefore, the fabricated membrane produces fresh cool air, and it can be applied as a guiding sample for designing the membrane-based dehumidifier with improved performance.


Introduction
Growth in the world's population and economy, coupled with rapid urbanization, will result in a substantial increase in energy demand over the coming years, a vast part of which is significantly attributed to the energy consumed by buildings for air conditioning [1]. The desired energy for such a demanding requirement, particularly in humid climates, is provided by burning fossil fuels whose consumption leads to serious energy and environmental problems [2]. Additionally, short-or long-term exposure to air pollution has been associated with a wide range of human health effects. Hence, paying less attention to control measures that reduce fossil fuel consumption could result in human health hazards, discomfort, and productivity decline of human beings [3]. One factor that is also important to be considered is the climatic parameter such as relative humidity. It appears that effective control of humidity stands as a critical matter for maintaining a healthy indoor environment and reducing energy consumptions [4]. Recently, there have been a number of methods to reduce the amount of moisture in the air such as moisture condensing method, desiccant-wheel method, and liquid desiccant method which are considered as conventional air dehumidification procedures [5]. In comparison, the liquid desiccant method which utilizes liquid desiccant such as lithium chloride (Li Cl) [6], lithium bromide (LiBr) [7], and magnesium chloride (MgCl 2 ) [8] for the adsorption of water vapor in air not only results in lower energy consumption, but also provides a faster and more efficient dehumidification due to its lower regeneration temperature compared to other methods [9]. The advantageous liquid desiccant method suffers from the issue of carryover of the liquid desiccant which could lead to the corrosion of operating components [10]. Therefore, it is essential to assess a method or promoter with lower energy consumption and higher regeneration efficiency in air conditioning systems [11].
Nowadays, the membrane-based liquid desiccant air dehumidification which had been mostly applied in wastewater treatment processes has become more interesting for researchers worldwide, who are investigating the field of air conditioning systems [12]. This technology, which is a thermally driven separation process, utilizes the vapor pressure difference across the porous membrane as a driving force to achieve successful separation and dehumidification [13]. The most important difference between this method and the conventional liquid desiccant air dehumidification is the utilization of semi-permeable membranes for the separation of the processing air and liquid desiccant fluids which generally results in the successful inhibition of liquid from penetrating [14]. Due to this phenomenon, the desiccant droplets cannot enter the air stream which further prevents the corrosion of the air conditioner components [15]. Moreover, this separation suppresses the carryover of liquids and the by-product gases of the operation [16]. As a result, this technology has various advantages compared to conventional air conditioning technologies such as (1) higher air dehumidification/purification capability [17], (2) pollution-free structure [18], (3) low energy consumption [19], (4) operability at lower temperatures (40-80 °C) [20], and (5) negligible heat transfer through the membrane. It is clear that the production of membranes for improved and optimized structure and performance has become highly desirable.
The superior performance of polymeric membranes have been proved by the results obtained from the research work of Abdel-Salam et al., who reviewed the application, energy saving, and economic and environmental impact of liquid membrane energy exchangers. They demonstrated that in comparison to the conventional liquid desiccant systems, their energy consumption is reduced by 10-50% [21]. Moreover, Baris Kavasogullari et al. designed and manufactured an eco-friendly liquid desiccant system, which is membrane contained polycarbonate boards considered as innovative packing material in an open-cycle calcium chloride-water system [22]. The proposed system had a higher surface tension compared to other conventional polymeric materials. In further investigation of obtaining more reliable results concerning the performance of the systems, Dec hang Wang et al. introduced a mathematical model that was experimentally validated for a system containing a hollow fiber membrane in its energy storage unit [23]. The project performed by Hongyu Bai et al. was also attributed to the evaluation of the performance of a membrane-based dehumidifier applying an experimental design with a datasheet of operating parameters and their corresponding impacts on the operation [24]. The results indicated that the efficacy of the humidifier increased by increasing the flow rate ratio of the solution to air, but decreased by increasing the operating temperature.
Based on our literature review, the production of membranes with improved and optimized structure and performance has the ability to reduce energy consumption and component corrosion remarkably. However, there are a limited number of studies investigating the impact of the operating conditions and the factors governing the fabrication of membranes on the performance of the dehumidifier that is based on designing a reliable experimental method. Thus, we decided to conduct a research project to primarily optimize the system performance that could also be qualified as a role model for future production. This goal can be achieved by evaluating the control factors with experimental design and perform mathematical model simulations.

3
The objective of this paper is to establish an experimental design for the mathematical simulation of factors governing the structure and performance of porous membranes that are used in the dehumidification of air conditioning systems with improved performance in energy use and environmental impact. Firstly, the structure and performance of the fabricated membrane were designed and then the reasons governing such a capability for avoiding the carryover and corrosion problems were discussed. The reliability of the model was further evaluated using COMSOL. Overall, the system application along with the obtained results was used to provide support for the capability of the fabricated membrane in solution desiccant air conditioning systems, especially in humid areas.

Methodology
The conventional experiments with one variable at a time were applied as an experimental design to optimize the structure and the performance of the membranes in air conditioners. This procedure was highly beneficial for minimizing the variance associated with the estimates of the models with fewer experimental runs. According to the results reported in the literature, the porosity of the membranes, polymer additives, humidity, temperature, and the flow rate are the most important factors governing the performance of the air conditioner. Therefore, the levels of 15, 17.5, and 20% for PSU concentration, 68 and 79% for humidity, 60 and 120 L/h for flow rate, and 312.15, 313.15, and 321.15 K for temperature were specified for designing our experimental procedures. Furthermore, the reliability of the models was evaluated by comparing the results obtained from the experimental runs and the calculated data from the COMSOL analysis. Finally, the most appropriate model was selected for analyzing the structure and performance of the proposed membrane by determining the most optimized operating points of PSU concentration, solvent types, temperature, humidity, and flow rate. This optimization is considered to be highly beneficial for proposing a fabrication route to achieve a membrane with the maximum performance capabilities possible. In the following section, process modeling, membrane fabrication and characterization, and dehumidification test are discussed, respectively. Figure 1 shows the membrane-formed module of parallel-plate channels along with the air and solution stream flow through the channels in a parallel-flow arrangement to facilitate duct sealing. Accordingly, a membrane and two neighboring flow channels were selected for the calculation domain to achieve the desired symmetry and calculation feasibility of the proposed model. Moreover, the COMSOL 5.3 software was used to perform two-dimensional modeling for flat membrane systems employing the laminar flow regimen of a single-phase fluid for the measurement of the velocity and pressure. The incompressible and compressible flows, as well as non-Newtonian fluids, are supported by physics interface at low Mach numbers (< 0.3). The laminar flow interface was solved by using the Navier equations and continuity equation to maintain the momentum and mass, respectively. The Navier equation in its most general form was employed as follows: The characterization of flow variables and fluid features of each point inside the medium is defined by using the

Process modeling
(2) (6) −PI + ∇U + (∇U) T n = −P 0 nP 0 ≤ P 0 . Fig. 1 Module of parallel-plate channels mean and actual values of variables. The obtained properties of a given volume around the designated point in porous domains were also utilized for achieving a significant definition. It is worth mentioning that the control volume should be large enough to the extent that exceeds the volume of pores and elements of the solid matrix. However, it needs to be small compared to the volumes and the dimension of the macroscopic problems as defined in the model. A fraction of the control volume occupied by pores defined as porosity which could vary in a very small magnitude in regions that are purely solid to integrate with free flow domains. Volume averages corresponding to the pore unit volume is defined as the physical property of the fluid such as values of density and viscosity. The specified values that exist in the neighboring free flow were designated to be corresponding physical parameters, which are experimentally measurable. Overall, the flow velocity is defined as the volume average referring to a unit volume of the medium including pores and matrix. It is also referred to Darcy velocity that is defined as the cross-sectional volume per unit of medium flow rate.
Finally, based on the above definitions, the velocity continues throughout the boundaries from porous regions to free flow ones. The Darcy velocity and pressure are the dependent variables of the Brinkman equations. The combination of continuity equation with the momentum equation governs the flow of the porous media. Brinkman equation is derived from the combination of equations as follows: Furthermore, to measure the concentration of the diluted solute, and the reaction rate of gas, liquid, or solid dissolved species, the transportation rate of the diluted species was measured and utilized. The Fick law of convection was further applied to obtain the transport driving force and diffuse migration by coupling with those of fluid flow and electric field, respectively. The utilized equations of Fick law are as follows: Inflow-air equation: Inflow-Li Cl equation: Apart from the convection and diffusion equations, the node has the functionality to measure the evolution of species through adsorption and dispersion mechanisms.
Membrane equation: Air equation: Li Cl equation:

Membrane fabrication and characterization
In this project, the wet-phase inversion procedure was applied to fabricate the designated membranes. Firstly, casting solutions with different solvents such as NMP and DMF were prepared. The prepared solutions were mixed using a magnetic stirrer for 24 h at room temperature. Secondly, a casting machine was employed to disperse the obtained solution on a glass support. Finally, the obtained polymer films floated on a coagulation bath containing deionized water until their inversion process was complete. It is worth mentioning that by knowing the pore volume of the membranes and their membrane volume, it is possible to measure its porosity. Therefore, Density of dry and wet membranes was measured and the results were recorded which demonstrate a weight difference that could be due to presence of eater in its cavities that is the result of the density of water in the membrane. Figure 2 shows the dehumidification setup, which was used for the measurement of outlet air (RH) and the evaluation  Table 1 shows the Membrane properties, contactor specifications and transport properties. As stated earlier, the membrane module was assembled using a membrane with two different flow fields for humid air and liquid Li Cl. The most important control factors such as flow rate, pressure, and temperature of gas and LiCl, were controlled. A humidity indicator (digital thermometer and hygrometer-DT-3) was also employed to measure the gas RH. Since the membranes have a porous structure, the pressure was kept constant on both sides of the membrane, while the transfer of water vapor responded to the difference in chemical potentials. Additionally, a temperature indicator controller was embedded inside, and a mechanical mixer was utilized to keep the temperature constant at the specified temperature points of 312.15, 313.15, and 321.15 K. The gas flow tube was a serpentine channel with a specific length, 6 mm width, and 5 mm surface contact with the membrane. The effective surface of the membrane in the module was 5 × 1800 mm 2 . The LiCl flow field was a simple container with a reticular plate to hold the membrane. In each run, the samples of the membranes with the size of 15 × 15 cm were examined in different designated conditions for mechanical tests. The effects of the operating conditions and the membrane synthesis parameters on dehumidification were evaluated using the apparatus described earlier. The composition of the casting solution was an important factor affecting the position, shape, and orientation of the formed cavities, which had certain impacts on the dehumidification performance. Accordingly, six membranes were prepared, and their attributed experimental data are shown in Table 2.

Experimental validations
As indicated in Fig. 3a, the dehumidification setup was prepared and operated at different operating conditions to evaluate their impact on the performance and characteristics of the membranes (Table 2). Therefore, a module was designed using the COMSOL software, which facilitated validating the experimental results from the mathematical perspective. All the experimental runs were performed isothermally for at least three times to confirm the significance of the results  Table 2 The weight percentage of the polymer and the composition of solvent in different membranes  Table 3. Hence, we obtained the calculated outlet humidity, and mass fraction of the air, as well as solution flows from the humidifier [25]. Overall, there was less than 6.3% discrepancy between the mathematical results and the experimental methods. From this evaluation, we could confirm the reliability of the module in predicting the performance of the membranes [26]. Furthermore, the proposed model by the COMSOL software was able to successfully predict and solve the heat and mass transfer problems in a flat membrane-based liquid desiccant air dehumidification model. The proposed module can be expressed as follows: where pout,ideal is the absolute humidity which would be ideal for the initial airflow of the converter. By assuming that the air is completely dehumidified, most likely pout,ideal would be close to zero.

Simulation results
The relative humidity of the output stream that was obtained from the mathematical simulation is listed in Table 3. Additionally, Fig. 4a shows the channel layout in the simulation in which its sequence type was selected as the main physical mesh control, while its element size was considered to be relatively coarse [27]. Figure 4b shows the variations of velocity along the channel. It is clearly revealed that the velocity of the channel in the middle was maximal and the velocity along the walls was zero [28]. Figure 5a also depicts the pressure variations along the channel, which due to the presence of the porous membrane remained constant across the channel with no noticeable variations. As shown in Fig. 5b, the effect of membrane composition induced variations in both air and soluble concentration along the length of the channel [29]. Accordingly, air concentration decreased by increasing the airflow through the membrane surface, while the solution concentration increased, indicating the transfer of moisture from the air to liquid [30].  Figure 6a shows that the impact of solvent type and porosity of the membranes varied due to the presence of DMF or NMP with different weight percentages on the porosity of the membranes. According to the obtained results, the membranes with DMF solvent had high dehumidification, and the relative humidity of the outlet air substantially lower compared to the membranes containing the NMP. This clearly indicates that the dehumidification efficiency of the former samples is higher than that of the previous one. Figure 6b also shows the positive impact of solvent porosity on the efficiency of the dehumidification procedure, as the amount of relative humidity of the outlet air was decreased significantly in the samples containing the DMF [31]. We further attempted to assess the observed  Figs. 7 and 8, the median pore diameter in the DMF containing membranes was smaller, and the number of pores was higher. Therefore, it can be concluded that the DMF had more interaction with moist air and thus resulted in a more enhanced air dehumidification yield [32]. According to Fig. 7, the bottom layer of membrane containing the DMF solvent was larger in size, which allowed less liquid desiccant to diffuse deeper into the membrane. It should be mentioned that unlike air, the membrane was not nearly capable of letting the liquid desiccant to diffuse into it, thus allowing an easier transmittance of the DMF solvent [33].

Effect of air flow rate
The impact of the inlet flow rate on the humidity of air and consequently the efficiency of the membranes in the dehumidification process are summarized in Table 3. It is noted in the literature that with increasing AIFR from 60 to 120 L/h and the stability of the discharge rate for all the synthesized membranes, the humidity of the air will be less affected and the absorption of moisture from the air besides the dehumidification rate will decrease as well [34]. When the AIFR increases, due to the constant cross section, air velocity increases as well; thus, the air stagnation time period decreases and air mass remains constant over an area for an extended period. The simulation results are presented in the following sections for different flows (Fig. 9). According to the findings, by increasing the airflow percentage, the air concentration decreases, and the efficiency of dehumidification also drops [35].

Effects of humidity and temperature
The impact of the increasing moisture while the temperate and input flow rate stay constant is shown in Fig. 10 and 11. Additionally, increasing the relative humidity results in the increase of the air concentration [36]. Equation (14) is defined as follows: The membrane was expected to increase the moisture flux as verified by the implication of the modules. It should be = X W * P∕P sat . noted that in Eqs. 14 and 15, N W is the moisture flux, K m is the overall mass transfer coefficient in the membrane. X W stands for the mole fraction of moisture in the air, W is the activity coefficient of moisture, Y W is the mole fraction of moisture in liquid, P and P sat are the total pressure and water saturated pressure, respectively, and φ is the relative humidity.
Also, if the relative humidity is kept constant and the temperature is raised, according to the diagram of the pyrometer, we noticed that the percentage of humidity in the air also increases. The same phenomenon is also applied to Eq. (14), as the transferred flux also increases following the same path [37].
According to Table 3, the overall results indicate that in all the membranes evaluated in this study, increasing the Fig. 9 Concentration of air outlet across the channel a and concentration variations along the membrane b for different air flow rates temperature and humidity from 39 to 48 °C and 68 to 79% caused the efficiency of dehumidification to rise from 42 to 48% in all the membranes [38]. Thus, it was concluded that in any area with a warm and humid climatic condition, the most appropriate efficiency could be achieved. Additionally, by increasing the amount of intake air, the relative humidity of the outlet air was also diminished in the synthesized membranes. Under other conditions, increasing the air velocity results in the reduction of the amount of air outlet [39]. In conclusion, the results of this study indicate that the synthesized membranes most likely demonstrate remarkable performance in hot and humid climates as compared with other fabricated membranes (Figs. 12, 13). Therefore, they can be considered as reliable membranes for the dehumidification in areas with extreme climactic conditions in terms of humidity and temperature [40].

Conclusion
Recently, increasing attention has been paid to membranebased technology due to the fact that air dehumidification plays an important role in improving air quality and maintaining thermal comfort. Hence, membrane-based air dehumidification has been applied in ventilation and air conditioning. In this study, a conventional experimental design was developed to optimize the performance and structure of six samples of fabricated membranes for utilizing under diverse climatic conditions.
We conducted a series of polymeric porous membranes fabricated with different compositions that were prepared with various weight percentages of polysulfone (PSU), mixed with N-methyl-2-pyrrolidone (NMP) and dimethyl form amide (DMF) solvents. Additionally, the operating control factors such as flow rate, temperature, and humidity were controlled and optimized. Overall, the results indicate that the membranes with higher porosity are capable of dehumidifying the air to a higher extent. The selected membranes further demonstrated an improved performance in dry and humid areas with lower intake flow rates. Therefore, the fabricated membrane produces fresh cool air, and it can be applied as a guiding sample for designing the membrane-based dehumidifier with improved performance.