Comparative Study of the Heat and Mass Transfer Characteristics between Counter-Flow and Cross-Flow Heat Source Towers

: The heat source tower (HST), as a cleaner energy production, which can absorb the low-grade energy from ambient air to drive the heat pump unit without emissions has attracted more and more interest. In addition, HST has excellent economic applicability by using cooling tower equipment, which was idle in winter. However, there are few studies on comparative analysis of thermal behavior between counter-flow and cross-flow HST. A mathematical model suitable for both HST types was developed to identify the performance discrepancies between them. Then a parametric study was carried out in order to investigate the impact of solution and air as well as packing material properties on energy transfer of HSTs. Finally, the characteristics of solution dilution and dehumidification were investigated. As the inlet solution temperature increases, increases first, then decreases gradually, but a transition point occurs in the solution at − 5 °C. Moreover, the transition section of moisture transfer direction for counter-flow HST was located in the 0.78 m and 0.26 – 1.56 m of packing material height, under the condition that the air relative humidity was 50%. In summary, this work intuitively indicates the thermal performance difference between counter-flow and cross-flow HST, also could assist the selection of proper operating conditions in HSTs.


Introduction
A heat source tower (HST), also known as a reversibly used cooling tower (RUCT) [1][2][3], absorbs low-grade energy from the air and supplies the heat to a building. It is generally believed that the heating tower is applicable in multiple countries and regions, especially for areas with high outdoor temperatures and humidity [4]. HST has several advantages compared with traditional heating units. One of the advantages is that HST can absorb low-grade energy from the atmosphere. In addition, the heat pump unit can drive low-level energy, which is more environmentally friendly. Furthermore, the conventional heating unit produces a series of toxic and harmful substances through combustion, which seriously affects the natural environment [5]. Hence, recent concerns about environmental pollution have greatly increased the attention paid to cleaner heating units like HST. In fact, the thermal behavior discrepancies between cooling towers and HSTs are on account of several aspects. First, the direction of energy transfer is reversed. A cooling tower removes the excess The aim of our research is to identify the influence of air and solutions as well as packing material properties on the thermal behavior within HST. A mathematical model validated by experimental data was developed to predict the latent and sensible heat transfer of a HST. The influence of the properties on sensible and latent heat transfer was examined. Additionally, solution dilution and dehumidification characteristics were studied. This research could assist the selection of proper operating conditions for existing HSTs.

Description of the Two Types of HST
The schematic diagrams of the adiabatic cross-flow HST module and the counter-flow HST module are shown in Figure 1, respectively. The solutions are sprayed uniformly onto the packing materias from the top surface to the bottom of the HST by the nozzles. According to the air flow direction in the HST, HSTs are divided into two types. When air enters from the side and flows out from the opposite side it is a cross-flow HST, and if air enters from the bottom side and out of the top side it is a counter-flow HST. The energy transfer arises on the packing material between the solution and the air. When the HST is operating in summer, the air flows through the material in the HST, where the air is humidified and heated, so that the enthalpy of the air increases while the solution temperature decreases. When the HST is operating in winter, the circulating heat transfer medium is replaced by an antifreeze solution with a freezing point below 0 °C. Unlike a HST running in the summer, the air will be cooled and dehumidified, and the solution concentration is diluted while the temperature of the solution rises. Figure 2 illustrates the operating mode of the HST system. The red arrow highlights the heating condition and the blue arrow outlines the cooling condition. When the HST system is operating during the summer, the cooling of the circulating heat transfer medium was carried out by direct evaporation. After pumping the cooling water to the condenser, the heat released by the heat transfer circulating medium is absorbed by the refrigerant in the compressor, and then the circulating cooling water flows out of the condenser.
The HST system collects low-grade heat energy in the outdoor air in winter, the circulating heat transfer medium is an aqueous solution of calcium chloride. The antifreeze solution flows out of the heat pump unit and then flows into the HST. After absorbing heat from the air, the solution enters the evaporator through the water pump, the refrigerant gains heat through the phase change process, and the circulating process occurs again to transfer the heat.

Establishment of Mathematical Equations
A mathematical model was developed to investigate the thermal performances of cross-flow and counter-flow HSTs. The major simplifications and assumptions for the mathematical model are as follows: (1) The specific heat and density are constant with temperature [12,20]; (2) The heat and mass exchange through the tower surface to the environment are negligible; (3) Heat and mass transfer perpendicular to the flow direction is ignored; (4) The thermal resistance between the solution and the air occurs on the airside; (5) The properties of the packing material are uniform and stable; (6) There is no loss of solute during HST operation; (7) Both inlet air and solution flow rates are steady and uniform.
Based on the above assumptions, the cross-flow HST can be simplified into a two-dimensional model, and the counter-flow HST can be simplified into a one-dimensional model. Because both the inlet air and solution flow rates are steady and uniform, any one of the control volumes can be selected for analysis throughout the entire model. The control volumes of a cross-flow and counterflow HST are shown in Figure 3.

Development of Governing Equations
The control volume is shown in Figure 3. Every control volume follows the conservation of energy. Therefore, for any of the control volumes, there are the following rules: The value of variables that fluids enter the control volume + the increment of variables produced by the interaction between control volume and fluids = the value of variables that fluids leave the control volume + the cumulative value of variables in the control volume. The packing material is thermally insulated and has no internal heat source inside, and then the energy is completely released between the air and the solution. Therefore, the energy generated and accumulated by the control volume is none [6].
Extracted to the conservation of energy, the enthalpy of the inlet and outlet should be equal, thus, which can be written as Equation (1): The input enthalpy of the control volume can be written as below: The output enthalpy of the control volume can be written as Equation (3) Combining Equations (2) and (3), the energy conservation equation of the control volume can be rewritten as Equation (4): The mass of water brought into the control volume by air and solution is equal to that of water brought out of the control volume by air and solution. Similarly, the water vapor mass balance equation can be written as Equation (5) Simplifying Equation (5) we get Equation (6): It is assumed that the solute of the antifreeze solution does not change, that is, the solute has no loss [24], then Equation (7) can be written as: where ζ is the solute mass fraction of the antifreeze solution, %.
As can be seen from above, the properties of the packing material are uniform and stable. Assume that A is packing material heat and mass transfer surface. According to the literature, the value can be calculated as Equation (8): where β is the packing material porosity, it is the ratio of the void volume of the total volume. D is the packing material spacing, m. During the energy transfer process, the energy transfer of the air is caused by two effects. On the one hand, convective heat transfer occurs between air and solution. On the other hand, the water vapor undergoes a phase change process, the air condenses and liquefies and rapidly releases the latent heat of condensation. So, the energy transfer equation of the control volume can be written as Equation (9) Simplifying Equation (9) we get Equation (10): Similarly, during the mass transfer process, the difference in air humidity ratio is due to the convective mass transfer process, so the control volume moisture transfer equation can be written as Equation (11): Simplifying Equation (11) we get Equation (12): The enthalpy difference between moist air can be written as Equation (13): where Cp is specific heat of air (kJ·kg −1 ·K −1 ) and r is water latent heat of condensation at 0 °C (kJ·kg −1 ·K −1 ). The mass transfer unit number NTUm is a dimensionless number that measures the heat transfer capacity within different fluids. The Lewis number Le representing the ratio of thermal diffusion coefficient to the mass diffusion coefficient. It describes the extent of heat transfer and mass transfer during convective heat transfer. The NTUm and the Le can be written as Equations (14) and (15): The transfer coefficient varies linearly. Cp represents the specific heat of air (kJ·kg −1 ·K −1 ), hm is the evaporation coefficient. Combined with Equations (12)-(17) this becomes: Equations (4), (6), (7), (16) and (17) establish the major energy transfer model in a cross-flow HST. Therefore, through similar steps, the counter-flow conservation of energy Equation ( The heat transfer process can be written as Equation (20): The mass transfer process can be written as Equation (21): Taking Equation (16) into equations (20) and (21), and then simplifying it with equations (18) and (19), the conservation of energy can be rewritten as Equation (22): The mass balance equation can be rewritten as Equations (23) and (24): The heat transfer process can be rewritten as Equation (25): The mass transfer process can be rewritten as Equation (26):

Numerical Methods
In order to solve the equation using the finite difference method, the governing equation is rewritten into the form of a differential equation [25,26]. Figure 4 explains method intuitively. The section can be seen as a two-dimensional mathematical model, which divides the packing material into N parts of the x direction and into K parts of the y direction. For any control volume (i, j) in this model, the input parameters are the enthalpy of the solution hs i,j , air enthalpy ha i,j , solution mass flow rate ms i,j , air humidity ratio ωa i,j , and solution concentration ζ i,j , respectively. The output parameters are hs i + 1,j , ha i,j + 1 , For the (i + 1, j) or (i, j + 1) control volumes, the same calculation process is performed to obtain the import of the next volume, and then repeating the calculation process [25,26]. When the values of K and N are large enough, the result can be considered reliable [6]. For mesh (i, j), the Equations (22)-(26) can be discretized as follows: This model uses an aqueous solution of calcium chloride as fluid. The boundary layer generated while the solution contact with air, the energy transfer process occurs on the boundary layer, simultaneously. The humidity ratio of moist air in the boundary layer can be written as Equation (32): The air enthalpy can be written as Equation (33) The enthalpy of solution can be written as Equations (34)-(36) [27]: The values of the correlation coefficients are shown in Table 1. Combining equations (34)-(36) and Table 1, Equation (34) can be rewritten as below:

Determination of Other Parameters
For both forms of HSTs, the solution specific heat and the packing NTUm are the same, so the two forms of packing materials can be discussed together.

Solution Specific Heat
For aqueous calcium chloride solutions in the non-crystalline region, the specific heat can be calculated using Equations (38)-(42) [27]: where, ξi (i = 1, 2, ..., 6), a-e and f are all constant value, and their values areas below in table 2.

Packing Material NTUm
The mass transfer coefficient is difficult to determine in the Equation (14). And the relation between Sherwood number (Sh) and hm can be obtained from the following formula: in which ds is the equivalent diameter of the channel as defined by Al-Farayedhi et al. [28]. Therefore, the Sherwood number (Sh) was introduced into Equation (14). The Sherwood number of the packing material can be obtained by fitting the experimental data. So, NTUm can be rewritten as Equation (44)

Heat and Mass Transfer Coefficient
According to literature [26], the heat transfer coefficient hc can be calculated as follows: A program written in MATLAB (R2018b, The MathWorks, Inc., Natick, MA, USA) is used to solve the mathematical models. The computational flow of the program is illustrated below in Figure  5.

Experiment and Model Verification of the HST
To validate the developed mathematical model, field experimental work has been carried out by our group in Chongqing City. Figure 6 shows the major structure of the experimental platform. The experimental device used in the experiment is shown in Table 3.  The packing material used was a structured PVC (polyvinyl chloride) material with dimensions of 5000 mm long × 3300 mm wide × 2900 mm high for counter-flow, and 1460 mm long × 970 mm wide × 1060 mm high for cross-flow. The counter-flow HST experiment was carried out during the 11:50-18:30 period within five days under outdoor inlet air conditions, while the cross-flow experiment was carried out during the 14:30-18:15 period.
Ts,o, Ta,o and ωa,o were selected for verification. Figure 7 demonstrates the difference between the simulation data and the experimental data. The R 2 , r and RMSE (root mean square error) were carried during the model verification. The model verification of the cross-flow model is based on the experimental data from Appendix A, and the model was also used for counter-flow verification using the data from Lu [29]. The parameters used to quantify the packing material are mainly the porosity and packing space. The influence of these parameters on the thermal behaviour of the filler was discussed in Section 4.3. Different flow directions cannot directly affect these parameters, but they will cause different heat and mass transfer processes between the air and the solution inside the packing material, which indirectly affects the thermal performance of HST. Therefore, the thermal behaviour of a counterflow HST and a cross-flow HST are different.

Parametric Study
Through the mathematical model verified in the previous section, the mass and energy transfer of HST can be predictable, to compare the thermodynamic characteristics between the cross-flow and counter-flow HST, the same inlet boundary conditions are used for the calculation.
It can be reflected in Figure 8 that the air temperature decreases as the solution temperature increases. This figure shows the distribution of air and solution temperature inside the packing material. The form of this figure is a box diagram, where the box is used to represent the parameter distribution of cross-flow at the same height, the upper and lower boundaries are the maximum and minimum values of a set of data, and the upper and lower limits of the box are a set of 1/4 quantile of data. This phenomenon mainly due to the temperature difference between the air and the solution. Therefore, SHT (sensible heat transfer) occurs due to the temperature difference, and latent heat transfer (LHT) occurs simultaneously. The air temperature decreases from 6.0 °C to 1.88 °C with the solution temperature increases from −3.0 °C to −1.126 °C. In addition, for the cross-flow HST, the air temperature decreases from 3.82 °C to 3.09 °C as the solution temperature increases from −3.0 °C to −0.924 °C. Table 4 shows the Simulation boundary conditions.   Figure 9a derives the temperature contour of the air in the HST. The air temperature decreases along the airflow direction (X-direction), and the air temperature changes slightly following the Yaxis. The maximum temperature appears in the lower-left area of the X-Y plane, where the temperature difference between air and solution is the minimum. However, in the upper right corner of the X-Y plane, the air temperature taken lowest value here, which is caused by the maximum temperature difference between the solution and air in this area. Further, a solution temperature contour can be obtained as shown in Figure 9b. Compared with Figure 9a, the solution temperature decreases along the Y-direction, and the air temperature changes slightly along the X-axis. A rapid change of the solution temperature appears in the lower-left part of the X-Y plane. The temperature difference between the two fluids on the left side of the positive X-axis region is higher, which can be the major cause of the solution temperature evolution. The variation of the air humidity ratio is revealed in Figure 10. The air humidity ratio in both HSTs decreased along the fluid flow channel. The counter-flow air humidity ratio decreases from 4.50 g/kg to 3.82 g/kg and the cross-flow HST decrease is from 4.024 g/kg to 4.157 g/kg. Part of the reason for this discovery is that a pressure difference exists between the water vapor surface and the solution surface, so consequently the moisture in the air is continuously transmitted to the solution side under the pressure difference. In addition, the LHT is spontaneously generated in the process. Moreover, the cross-flow humidity ratio decreased faster than in the counter-flow configuration. The evolution of the solution mass flow rate is displayed in Figure 11. The solution mass flow rate varies from 37.70 kg/s to 37.73 kg/s for cross-flow, while the counter-flow solution mass flow increased to 37.726 kg/s. A major source for this is that the solution is diluted under the water vapor pressure difference, the solution temperature is far below the freezing point, so the solution will continuously absorb moisture from the air. However, as the freezing point approaches 0 °C, the system will have a risk of freezing, so the characteristics of solution dilution should be investigated.

Impact of Inlet Air Parameters on Energy Transfer
A parametric study is carried out using the above constant boundary conditions as benchmarks. In order to compare the thermodynamic properties of two HST, the SHT, LHT and ξ of two HST under different air, solution and material boundary conditions will be investigated.

Impact of Inlet-Air Temperature
The energy transfer inside the HST is indicated in Figure 12. With the inlet air temperature elevates, the SHT in the HST increases, while the LHT decreases slightly and the difference of SHT between the two HSTs also increases. The counter-flow SHT increases from 73.2 kW to 281.2 kW, while the SHT varies from 79.81 kW to 312.5 kW for the cross-flow setup. However, the counter-flow LHT decreases from 69.16 kW to 64.37 kW when the cross-flow varies from 78.7 kW to 70.45 kW. This suggests that the air inlet temperature has a significant effect on the SHT, but has little influence on the LHT. ξ indicates the proportion of LHT caused by moisture transfer in THT. With the inlet air temperature increases, not solely counter-flow ξ decreases from 0.4858 to 0.1863 significantly, but also cross-flow ξ decreased from 0.4965 to 0.184. As the air inlet temperature increases, resulting in a reduction in the air humidity ratio, the capacity of the air moisture content also decreases, following, a decrease in the LHT further triggered by a decrease in the mass transfer between the air and the solution, finally the ξ decreases gradually. Accordingly, the inlet air temperature has a majority effect on ξ, the cross-flow ξ is more sensitive to the inlet air temperature subsequently.

Impact of Inlet air Humidity Ratio
The effect of inlet air humidity ratio on energy transfer in the HST is presented in Figure 13 when the inlet air humidity ratio varies from 2.6 g/kg to 5.4 g/kg, the LHT in the HST increases, while the SHT slightly attenuate, the difference of LHT between two HSTs also increases. It is contrary to inlet air temperature, the counter-flow SHT decreases from 170.4 kW to 166.6 kW, with the cross-flow SHT varies from 190.5 kW to 183.2 kW. Furthermore, the counter-flow LHT increases quickly from −8.858 kW to 102.5 kW when cross-flow varies from −12 kW to 116.3 kW, in contrast to the inlet air temperature, the inlet air humidity ratio causes a pivotal contribution to the LHT but has little amelioration on the SHT. It is distinct from the ξ of air inlet temperature, the counter-flow ξ rises from −0.055 to 0.381 along with the cross-flow ξ increases from −0.067 to 0.389. The LHT between the fluids enhanced as the inlet air humidity ratio increased, However, the increase in air humidity doesn't drive the air temperature increases. In other words, the temperature difference is still constant, thus the SHT remains almost unchanged and even a decline is obtained. Nevertheless, the water vapor partial pressure highly increased with the air humidity ratio amplified, then the increase of LHT results in an ξ increment, especially for the cross-flow HST. Hence, the inlet air humidity ratio has a significant influence on ξ. In spite of the LHT increase, the solution is also diluted, therefore, the HST system is at risk of freezing.

Impact of Inlet Air Mass Flow Rate
With the inlet air mass flow rate rises, both SHT and LHT increase, the difference between SHT and HST among the two HSTs also increased, which is shown in Figure 14. The counter-flow SHT increases from 68.12 kW to 244.1 kW, as well as the cross-flow SHT, which varies from 73 kW to 271 kW. Besides, the counter-flow LHT increases from 27.75 kW to 87.54 kW when the cross-flow varies from 31.64 kW to 106.2 kW. This conveys that the air inlet mass flow rate leads to a functional increase in both SHT and LHT. By contrast, the counter-flow ξ decreased from 0.287 to 0.265 with the crossflow ξ reduced from 0.302 to 0.282. Different from inlet air temperature and humidity ratio, there is no direct effect on the air enthalpy, on the contrary, the heat and mass transfer coefficients are indirectly changed by the air mass flow rate, which can enhance the SHT and LHT, owing to the SHT, which increases much more than LHT, so ξ decreases, particularly in a counter-flow configuration. Importantly, the energy transfer times were reduced even though the energy transfer coefficients were strengthened, as a consequence, decreases in the increase rates of SHT and HST were obtained.

Impact of Inlet Solution Temperature
The variation of SHT and LHT with solution inlet temperature varies is observed in Figure 15.  As the inlet solution temperature increases, ξ increases first, then decreases gradually. A transition temperature occurs at −5 °C, ξ is 0.284 and 0.289, especially for counter-flow and crossflow. Additionally, with the increase in solution temperature, LHT decreases rapidly. Although it can weaken the solution dilution, conversely, the energy transfer was reduced. As a result, the solution inlet temperature should be selected according to the HST system operating conditions.

Impact of Inlet Solution Mass Flow Rate
It was previously found that both SHT and LHT can strengthen significantly with the air mass flow rate is increased. The influence of the solution inlet mass flow rate on energy transfer characteristics is shown as Figure 16  When the inlet mass flow reaches 25 kg/s, the value of ξ almost remains constant. As a matter of fact, both SHT and LHT stay at a high value. In addition, a higher solution mass flow rate brings about a slower increment of energy transfer. Therefore, the solution mass flow rate should be selected precisely.

Impact of Material Packing Space
The influence of material packing space with SHT and LHT is presented in Figure 17. The counter-flow SHT decreases from 257 kW to 117 kW with the cross-flow SHT varies from 264.4 kW to 133.9 kW. Moreover, the counter-flow LHT decreases from 102.8 kW to 42.39 kW as the cross-flow varies from 112.8 kW to 52.98 kW. Furthermore, it is also seen that the value of ξ keeps decreasing slightly with the increase of material packing space. The cross-flow ξ decreases from 0.2991 to 0.2835 and from 0.2856 to 0.2666 for counter-flow. The above results indicate that the solution material packing space should be as small as possible under demand.

Impact of Material Porosity
The influence of material porosity with SHT and LHT is proposed in Figure 18. The counterflow SHT increases from 109.4 kW to 168 kW with the cross-flow SHT varies from 125.8 kW to 185.6 kW. Moreover, the counter-flow LHT increases from 39.6 kW to 62.52 kW when the cross-flow varies from 49.63 kW to 75.1 kW. Furthermore, it is also seen that the value of ξ keeps increasing slightly with the increase of material porosity. The counter-flow ξ increases from 0.265 to 0.272 and from 0.282 to 0.288 for cross-flow.

Investigation on Solution Dilution Characteristics
The above parametric study about air, solution and packing material revealed the effect on energy transfer, then the previous research has specifically shown a strong link between mass flow as well as humidity ratio and LHT. Nevertheless, LHT derives from moisture transfer between air and solution which may cause the solution freezing point to raise and even the HST system to break down. Thus, there are further reasons to investigate solution dilution characteristics thoroughly.
When the air relative humidity varies from 10% to 100%, the direction of moisture transfer also changes, which as shown in Figure 19. First, when the air relative humidity within 10%-49%, the solution in both cross-flow and counter-flow HST is dehumidified, while, the air is humidified. Second, when the air relative humidity is between 49%-60%, the direction of moisture transfer will be reversed in counter-flow HST, in other words, the solution will be diluted and dehumidified for counter-flow and cross-flow, respectively. Third, when the air relative humidity higher than 60%, solution in both cross-flow and counter-flow HST are diluted, while the air is dehumidified. Therefore, it can be seen from the previous research that although cross-flow HST has a better energy transfer coefficient than the counter-flow, the counter-flow HST can operate under 49-60% air relative humidity, while the cross-flow HST can only operate above 60% air relative humidity in winter. Further analysis needs to be performed on the moisture transfer of internal HST under 50% air relative humidity conditions, in order to analyze the reasons why the mass transfer direction changed. The variation of air humidity ratio in the air flow ωa and air humidity ratio at the air-solution interface ωs along the fluid flow channel length under 50% air relative humidity are shown in Figure  20. As the air flows upward, the difference in humidity ratio between the solution boundary layer and the air decreases below 0.78 m height, for counter-current HST, meanwhile, the solution humidity ratio is always higher than the air humidity ratio, this suggests that the air absorbs moisture from the solution. Conversely, the solution absorbs moisture from the air above the 0.78 m. For crossflow, there is a transition zone, which lies between 0.26-1.56 m. Ultimately, the solution absorbed moisture from the air for whole packing material, and the moisture transfer rate is 0.8933 g/s and −11.01 g/s for counter-flow and cross-flow respectively. It can be concluded that moisture transfer occurs in different directions at different heights, the transition surface is located in the 0.78 m for counter-flow and 0.26-1.56 m for cross-flow of the material packing height under 50% air relative humidity.

Conclusions
This paper numerically investigated the energy transfer characteristics in the HST with packing material. A mathematical model validated by experiments was developed for predicting the HST, considering counter-flow and cross-flow simultaneously. The main conclusions are as follows: • A higher air inlet temperature can improve the energy transfer performance of the HST. The above studies illustrated that as the inlet air temperature increases, the SHT increases, while the LHT and ξ are decreased, especially for cross-flow. Besides, the air inlet temperature leads to a significant effect on the SHT and ξ, but has little influence on the LHT. Conversely, when the air humidity ratio increases, the LHT increases, while the SHT and ξ decrease. Moreover, the inlet air humidity ratio causes a pivotal contribution to the LHT and ξ but has little amelioration effect on the SHT. However, the air inlet mass flow rate leads to a functional increase in both SHT and LHT.

•
The results have shown that inlet solution temperature greatly decreases the ability of SHT and LHT, and both SHT and LHT decrease when the solution inlet temperature increases. Further, as the inlet solution temperature increases, ξ increases first, then decreases gradually, but a transition point occurs in the solution at −5 °C. In addition, the solution mass flow rate has little influence on SHT and LHT.

•
The variation of packing space is opposite to that of packing porosity. The SHT and the LHT are simultaneously reduced as the packing space increases. Meanwhile, the SHT and the LHT increase as the packing porosity increases. The influence of the packing spacing on the energy transfer is higher than that of the packing porosity. Additionally, counter-flow HST is less affected than cross-flow. The results have found that solution material packing space should be as small as possible under demand conditions. • The mass transfer amount and direction vary with the air relative humidity. The transition air relative humidity of the mass transfer direction is 50% and 60%, respectively, for the counterflow and cross-flow HST. It also can be concluded that moisture transfer occurs in different directions at different heights, the counter-flow transition surface is located at 0.78 m, while the cross-flow transition area lies between 0.26-1.56 m of the material packing height under 50% air relative humidity.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A