Energy Saving and Charging Discharging Characteristics of Multiple PCMs Subjected to Internal Air Flow

One essential utilization of phase change materials as energy storage materials is energy saving and temperature control in air conditioning and indirect solar air drying systems. This study presents an experimental investigation evaluating the characteristics and energy savings of multiple phase change materials subjected to internal flow in an air heating system during charging and discharging cycles. The experimental tests were conducted using a test rig consisting of two main parts, an air supply duct and a room model equipped with phase change materials (PCMs) placed in rectangular aluminum panels. Analysis of the results was based on three test cases: PCM1 (Paraffin wax) placed in the air duct was used alone in the first case; PCM2 (RT–42) placed in the room model was used alone in the second case; and in the third case, the two PCMs (PCM1 and PCM2) were used at the same time. The results revealed a significant improvement in the energy savings and room model temperature control for the air heating system incorporated with multiple PCMs compared with that of a single PCM. Complete melting during the charging cycle occurred at temperatures in the range of 57–60 °C for PCM1 and 38–43 °C for PCM2, respectively, thereby validating the reported PCMs’ melting–solidification results. Multiple PCMs maintained the room air temperature at the desired range of 35–45.2 °C in the air heating applications by minimizing the air temperature fluctuations. The augmentation in discharging time and improvement in the room model temperature using multiple PCMs were about 28.4% higher than those without the use of PCMs. The total energy saving using two PCMs was higher by about 29.5% and 46.7% compared with the use of PCM1 and PCM2, respectively. It can be concluded that multiple PCMs have revealed higher energy savings and thermal stability for the air heating system considered in the current study.


Introduction
Phase change materials have been widely utilized for thermal energy storage in various heating and cooling applications. One of these essential PCM utilizations is air heating used for buildings' air conditioning and indirect solar dryers. The primary function of thermal energy storage materials is to store the excess energy to be used at relatively low heat supply periods, which results in a reduction in energy consumption. In some applications, it is required to keep the temperature of space at specified limits due to the fluctuations in heat supply resulting from weather variations, as in the case of solar air dryers. To solve this problem, multiple PCMs can be used in air heating systems to control the temperature of the conditioned space at desired values. Various research works [1][2][3][4][5] have been conducted to evaluate the thermal systems' performances based on a single PCM used as a thermal storage material. Several parameters, such as inlet temperature, flow rate, and time of charging and discharging cycles, were investigated. The attained improvements were in the range of 24-57% in the charging and discharging processes and a variation in the range of 4-9 • C in space temperature when PCM was utilized for cooling and heating purposes. Wang et al. [6] and Krzysztof et al. [7] evaluated the heat of fusion, specific heat capacity, and various thermo-physical properties of a mixture of paraffin wax with nanomaterials and microencapsulated PCM slurry at various concentrations, respectively. Significant improvements in PCM performance were observed in these works. Nada et al. [8] experimentally investigated the free air cooling performance using PCM filled in plate containers during charging and discharging cycles. The required number of PCM plates, energy saving, air flow rate, and outdoor temperature were investigated. The results revealed that the discharge time decreases by 17-25% with the increase in air temperature and air flow rate. A solar air heater with a double pass was experimentally and numerically investigated using a solar simulator with and without paraffin wax by Salah et al. [9,10]. The system's charging and discharging process was tested under different air flow rates: 0.6-1.8 kg/min and constant solar irradiation. The results showed that the system's thermal efficiency with and without PCM was 63.53% and 68.23%, respectively. Many researchers [11][12][13][14] have focused on optimizing the performance of solar collectors integrated with PCMs placed under an absorber plate. Thermal energy efficiency and air heater temperature differences with and without PCM were evaluated in natural and forced convection. Modelling of the charging and discharging process of PCM integrated with heat exchangers was investigated by Peter and Cornelia [15] and Wisam et al. [16] for shell and tube and double-pipe heat exchangers, respectively. A significant improvement in heat exchanger performance and a reduction in PCM melting time in the range of 27-46% were obtained with the change in heat transfer fluid temperature from 60 to 80 • C. The numerical modelling of charging and discharging cycles of phase change materials integrated with air heaters and solar systems is covered in several studies [17][18][19][20][21][22]. The charging and discharging behavior of different phase change materials placed in capsules with various geometries is simulated using computational fluid dynamic (CFD) analysis. The numerical results, including PCM thermal performance, melting and solidification times, the configuration of the liquid-solid interface, and temperature field, have revealed an improvement in the storage and recovery of heat energy. The evaluation of the thermal performance of various thermal systems integrated with different phase change materials has been reported by many researchers [23][24][25][26][27]. The results are presented in terms of charging and discharging characteristics, thermal storage capacity, and heat transfer rate and show a significant effect of multiple PCMs on the prevention of the fluctuations in temperature of cool and hot storage systems. The vital features of phase change materials are the melting and solidification times as well as the latent heat of fusion. Therefore, the melting and solidification characteristics during the PCM charging and discharging process are investigated in many experimental and numerical research works [28][29][30][31][32]. Abduljalil et al. [33] presented an experimental study of PCM melting and solidification in a finned triplex tube heat exchanger to examine the charging process under steady and unsteady conditions. The results indicated that the melting time of PCM was reduced by 86% due to the change in the inlet temperature of the heat transfer fluid compared with 58% due to a change in the mass flow rate. Ajay and Vinod [34] analyzed the charging and discharging time of paraffin wax and myristic acid as phase change materials for a comparison. It was found that the melting process of myristic acid required less time (2 h and 25 min) compared with that required for paraffin wax (3 h and 7 min) under similar operating conditions. The behavior of PCM (Plus ICE H190) during melting and solidification process was investigated experimentally using a differential scanning calorimeter (DSC) and thermal gravity (TG) lab test apparatus [35]. The results of the melting and solidification temperatures and latent heat of fusion for the tested PCM revealed that using Plus ICE H190 as a phase change material was feasible in many applications at a temperature in the range of 170 to 200 • C. Kanimozhi et al. [36] introduced an experimental study to investigate the charging and discharging characteristics of paraffin wax as a phase change material placed in a storage tank fabricated from concentric steel cylinders. Heat transfer was enhanced using several copper tubes filled with PCM and inserted in the storage tank. Heat absorption, heat rejection, time and temperature during the melting and solidification of paraffin wax were evaluated under different fluid mass flow rates. It was observed from the reported research works that there is still a need to fill the gap in the research area of melting and solidification characteristics of multiple phase change materials utilized in the air heating systems. The present work focuses on evaluating the charging and discharging characteristics of two PCMs subjected to internal flow in an air heating system. The use of two PCMs will improve the energy savings of the air heaters utilized in building air conditioning systems and indirect solar dryers. The period of supplying the heat energy stored by PCMs for the control of the temperature fluctuations of a room model and keeping it at the desired value is also experimentally investigated in the current study.

Charging and Discharging of PCM Analysis
The air duct illustrated in Figure 1 is thermally insulated, thus, it is assumed that there is no heat loss from the duct to the surrounding air. Based on the energy balance between the hot air in the duct and PCM1 during the charging cycle (melting process), the energy released by the air and gained by PCM1 can be expressed as follows [8]: and: where m . a is the air mass flow rate (kg/s); m pc1 and m Al are the mass (kg) of PCM1 and panels aluminum sheets, respectively, t ch1 is the charging time (s) measured experimentally at a time interval (∆t = 5 min), q pc1 is PCM1 latent heat of fusion (J/kg), T pc1L is the temperature of the PCM1 at the liquid phase; and T pc1M is the melting temperature (melting point) of the PCM1 from the thermo-physical properties of the phase change materials as illustrated in Table 1. T s1i and T s1 f are initial and final aluminum sheets' temperatures of PCM1 panels measured at the start and end of the charging cycle, respectively. with PCM and inserted in the storage tank. Heat absorption, heat rejection, time and temperature during the melting and solidification of paraffin wax were evaluated under different fluid mass flow rates. It was observed from the reported research works that there is still a need to fill the gap in the research area of melting and solidification characteristics of multiple phase change materials utilized in the air heating systems. The present work focuses on evaluating the charging and discharging characteristics of two PCMs subjected to internal flow in an air heating system. The use of two PCMs will improve the energy savings of the air heaters utilized in building air conditioning systems and indirect solar dryers. The period of supplying the heat energy stored by PCMs for the control of the temperature fluctuations of a room model and keeping it at the desired value is also experimentally investigated in the current study.

Charging and Discharging of PCM Analysis
The air duct illustrated in Figure 1 is thermally insulated, thus, it is assumed that there is no heat loss from the duct to the surrounding air. Based on the energy balance between the hot air in the duct and PCM1 during the charging cycle (melting process), the energy released by the air and gained by PCM1 can be expressed as follows [8]: where . is the air mass flow rate (kg/s); m pc1 and m Al are the mass (kg) of PCM1 and panels aluminum sheets, respectively, t ch1 is the charging time (s) measured experimentally at a time interval (Δt = 5 min), q pc1 is PCM1 latent heat of fusion (J/kg), T pc1L is the temperature of the PCM1 at the liquid phase; and T pc1M is the melting temperature (melting point) of the PCM1 from the thermo-physical properties of the phase change materials as illustrated in Table 1. T s1i and T s1f are initial and final aluminum sheets' temperatures of PCM1 panels measured at the start and end of the charging cycle, respectively.   The first term on the right-hand side of Equation (2) represents the latent heat gained by PCM1, while the second and third terms represent the sensible heat gained by the PCM1 and the aluminum sheets of the PCM1 panels, respectively.
Heat energy (Q 1c ) gained by the first energy storage material PCM1 in the air duct during the excess heat energy period of the charging cycle (melting process) represents the energy saving of PCM1, and it can be determined using Equation (3): The heat energy (Q 2c ) gained by the PCM2 in the room model during the charging cycle represents the energy saving of PCM2, and can be estimated as follows [17]: where m ar is the mass of the air (kg) occupied by the room volume, T ari and T ar f are the initial and final room temperatures at the start and end of the charging cycle, respectively, and T s2i and T s2 f are initial and final surface temperatures of PCM2 aluminum panels measured at the start and end of the charging cycle, respectively. During the discharging cycle (solidification process), when the supplied air temperature is less than the desired value for air heating purposes (T ars ) the heat energy (Q 1d ) released by the energy storage material PCM1 and transferred to the air flowing in the duct is approximated using Equation (5). As the heat gained by the aluminum sheets is small compared with that gained by the PCM1, it can be neglected.
where the set point air temperature of the room model T ars was selected in the present study to be in the range of (35-40 • C), which represents the desired temperature for the air heating purpose or food product drying in the solar dryers. The discharging time t dc1 (s), represents the required time to release heat energy from PCM1 at a low air temperature period during the discharging cycle. This time must be determined from the start of the discharging period when the PCM is at the liquid phase and until the PCM temperature reduces and approaches thermal equilibrium with the surrounding air in the PCM solid phase. t dc1 can be estimated based on the energy balance between the air and PCM1 as expressed by Equation (6).
The heat energy (Q 2d ) released by the second energy storage material PCM2 and transferred to the air in the room model can be determined as follows: Fluids 2021, 6, 275 The required time (t dc2 ) to release heat energy by PCM2 during the low room air temperature period of the discharging cycle can be approximated using energy balance Equation (8): where T pc2L is the temperature of the PCM2 during the liquid phase, and T pc2M is the melting temperature (melting point) of the PCM2.

Experimental Setup and Methodology
The experimental work was conducted using a test rig fabricated in the present work to provide the specified conditions for the testing of multiple PCM characteristics. The test rig consisted of two main parts, an air duct of 800 mm in length, 80 mm in width and 100 mm in height, and a room model with dimensions 400 × 250 × 250 mm, as shown in Figures 1 and 2. The air duct is equipped with an electric heating coil of power 1000 W and a variable speed air fan to provide a hot air with a flow of 0.005-0.01 kg/s at temperature in the range of 30-65 • C. The air duct and room model were fabricated from aluminum sheets (k = 232 W/m • C) and insulated using glass wool (k = 0.038 W/m • C) to minimize the heat loss to the surrounding air. Two phase change materials, PCM1 and PCM2, of different melting points and latent heat of fusion were used, as illustrated in Table 1. The first PCM was paraffin wax filled in 4 rectangular panels with dimensions of 310 × 70 × 15 mm, which were fabricated from 0.5 mm thick aluminum sheets. The PCM1 panels were placed horizontally in the duct with equal vertical spacing of 7 mm to allow the air to flow around the panels. The second PCM was paraffin RT-42 filled in 4 aluminum panels of similar dimensions and positioned vertically in the room model as shown in Figure 3. Ten K-type thermocouples (−200 to 1250 • C) with a data logger of the BTM-4208SD model and accuracy of ±0.2 • C were used to monitor and record the temperatures at many locations in the test rig as depicted in Figure 1. Two thermocouples were inserted inside two panels located at different locations for both PCM1 and PCM2 to measure the mean temperatures of PCM1 and PCM2. A digital anemometer (Am802 model) with accuracy of ±3% m/s was used to measure the velocity of air flow in the duct. The specified temperature of the heated air in the duct was controlled using a variable transformer (25 amp. Variac) and a thermostat to regulate the power supplied to the electric heating coil.
where T pc2L is the temperature of the PCM2 during the liquid phase, and T pc2M is the melting temperature (melting point) of the PCM2.

Experimental Setup and Methodology
The experimental work was conducted using a test rig fabricated in the present work to provide the specified conditions for the testing of multiple PCM characteristics. The test rig consisted of two main parts, an air duct of 800 mm in length, 80 mm in width and 100 mm in height, and a room model with dimensions 400 × 250 × 250 mm, as shown in Figures  1 and 2. The air duct is equipped with an electric heating coil of power 1000 W and a variable speed air fan to provide a hot air with a flow of 0.005-0.01 kg/s at temperature in the range of 30-65 °C. The air duct and room model were fabricated from aluminum sheets (k = 232 W/m °C) and insulated using glass wool (k = 0.038 W/m °C) to minimize the heat loss to the surrounding air. Two phase change materials, PCM1 and PCM2, of different melting points and latent heat of fusion were used, as illustrated in Table 1. The first PCM was paraffin wax filled in 4 rectangular panels with dimensions of 310 × 70 × 15 mm, which were fabricated from 0.5 mm thick aluminum sheets. The PCM1 panels were placed horizontally in the duct with equal vertical spacing of 7 mm to allow the air to flow around the panels. The second PCM was paraffin RT-42 filled in 4 aluminum panels of similar dimensions and positioned vertically in the room model as shown in Figure 3. Ten K-type thermocouples (−200 to 1250 °C) with a data logger of the BTM-4208SD model and accuracy of ±0.2 °C were used to monitor and record the temperatures at many locations in the test rig as depicted in Figure 1. Two thermocouples were inserted inside two panels located at different locations for both PCM1 and PCM2 to measure the mean temperatures of PCM1 and PCM2. A digital anemometer (Am802 model) with accuracy of ±3% m/s was used to measure the velocity of air flow in the duct. The specified temperature of the heated air in the duct was controlled using a variable transformer (25 amp. Variac) and a thermostat to regulate the power supplied to the electric heating coil.

Test Procedure
The experimental tests were conducted for three cases: PCM1 (Paraffin wax) panels placed in the air duct were used alone at the first case; PCM2 (RT-42) panels placed in the room model were used in the second case; and in the third case, two PCMs (PCM1 and PCM2) were used at the same time. During the PCM1 charging cycle (melting process), the ambient air entering the duct was heated with an electric coil to 65 °C in the first case (PCM1) and passed through the PCM panels until the PCM1(paraffin wax) melted. The charging cycle was stopped when the temperature of PCM1 exceeded the melting point (57 °C). In the PCM1 discharging cycle (solidification process), the electric heating coil was switched off, and the ambient air entering the duct at a temperature of 29-31 °C was heated by the energy released from the PCM1. The ambient air was heated to the desired temperature in the room model (35-40 °C) via heat exchange between PCM1 panels and entering air until the PCM1 solidified. When the test system approached a steady operation at the start of each charging and discharging cycle, the air and PCMs temperatures, duct inlet air velocity, and power supplied to the electric heater were monitored and recorded at a time interval of 5 min during the testing period of 330 min. These test steps were repeated for the second case (PCM2) and the third case with two PCMs.

Results and Discussion
The results of the three testing cases were considered in the current study to investigate the energy saving and characteristics of PCMs used in the air heating systems. The first test was conducted using PCM1 alone in the air duct to store the excess heat energy during the melting process in the charging cycle when air was supplied at 65 °C. The stored energy was released to the air flowing into the duct at ambient temperature (29-31 °C) in the discharging cycle. The second test was conducted using PCM2 alone in the room model to minimize the fluctuations in room air temperature and energy saving during charging and discharging cycles. In the third case, both PCM1 and PCM2 were used for the same purpose and for the comparison of results. The duct inlet air temperature was maintained at a steady value of 65 °C, while the air temperature of the room model was controlled to be in the range of 35-40 °C. This temperature range is desired in the air heaters and food products drying systems to simulate the output air temperature of solar air heaters working under Iraqi conditions (33.2232° N, 43.6793° E). The solidification steps for 50 g of PCM1 (paraffin wax) and PCM2 (RT-42) heated to 65 °C and cooled at ambient air temperature (Tam = 30.6 °C) are shown in Figure 4. This test was conducted to determine the melting points of PCMs during the discharging cycle under the current work's operating conditions. It can be observed for the PCMs at the liquid phase (65 °C) that the temperatures of complete solidification were in the ranges of 43-45 °C and 35-37 °C for PCM1 and PCM2, respectively.

Test Procedure
The experimental tests were conducted for three cases: PCM1 (Paraffin wax) panels placed in the air duct were used alone at the first case; PCM2 (RT-42) panels placed in the room model were used in the second case; and in the third case, two PCMs (PCM1 and PCM2) were used at the same time. During the PCM1 charging cycle (melting process), the ambient air entering the duct was heated with an electric coil to 65 • C in the first case (PCM1) and passed through the PCM panels until the PCM1(paraffin wax) melted. The charging cycle was stopped when the temperature of PCM1 exceeded the melting point (57 • C). In the PCM1 discharging cycle (solidification process), the electric heating coil was switched off, and the ambient air entering the duct at a temperature of 29-31 • C was heated by the energy released from the PCM1. The ambient air was heated to the desired temperature in the room model (35-40 • C) via heat exchange between PCM1 panels and entering air until the PCM1 solidified. When the test system approached a steady operation at the start of each charging and discharging cycle, the air and PCMs temperatures, duct inlet air velocity, and power supplied to the electric heater were monitored and recorded at a time interval of 5 min during the testing period of 330 min. These test steps were repeated for the second case (PCM2) and the third case with two PCMs.

Results and Discussion
The results of the three testing cases were considered in the current study to investigate the energy saving and characteristics of PCMs used in the air heating systems. The first test was conducted using PCM1 alone in the air duct to store the excess heat energy during the melting process in the charging cycle when air was supplied at 65 • C. The stored energy was released to the air flowing into the duct at ambient temperature (29-31 • C) in the discharging cycle. The second test was conducted using PCM2 alone in the room model to minimize the fluctuations in room air temperature and energy saving during charging and discharging cycles. In the third case, both PCM1 and PCM2 were used for the same purpose and for the comparison of results. The duct inlet air temperature was maintained at a steady value of 65 • C, while the air temperature of the room model was controlled to be in the range of 35-40 • C. This temperature range is desired in the air heaters and food products drying systems to simulate the output air temperature of solar air heaters working under Iraqi conditions (33.2232 • N, 43.6793 • E). The solidification steps for 50 g of PCM1 (paraffin wax) and PCM2 (RT-42) heated to 65 • C and cooled at ambient air temperature (Tam = 30.6 • C) are shown in Figure 4. This test was conducted to determine the melting points of PCMs during the discharging cycle under the current work's operating conditions. It can be observed for the PCMs at the liquid phase (65 • C) that the temperatures of complete solidification were in the ranges of 43-45 • C and 35-37 • C for PCM1 and PCM2, respectively.  Variations in the temperatures of PCM1 and PCM2 with time during the charging cycle at a duct inlet air temperature of 65 °C are shown in Figure 5. A constant trend in the temperatures can be observed in Figure 4 at 60 °C and 43 °C during the periods of 140-260 min and 120-220 min for PCM1 and PCM2, respectively. This behavior of the PCMs as energy storage materials can be attributed to the difference in the heat capacity based on various latent heats of fusion for these PCMs. Complete melting occurred during the charging cycle occurred at temperatures in the ranges of (57-60 °C) and (38-43 °C) for PCM1 and PCM2, respectively, which are identical to the reported PCMs melting-solidification results. A similar trend in the solidification temperatures of the PCM1 and PCM2 can be observed in Figure 6, which shows the variations in the temperatures of PCM1 and PCM2 with time during the discharging cycle at ambient air temperature of 30 °C. A steeper trend for the reduction in the temperatures can be seen before the discharging times of 80 min and 160 min for PCM1 and PCM2, respectively, due to the decrease in sensible heat of PCMs occurring at these discharge periods. Variations in the temperatures of PCM1 and PCM2 with time during the charging cycle at a duct inlet air temperature of 65 • C are shown in Figure 5. A constant trend in the temperatures can be observed in Figure 4 at 60 • C and 43 • C during the periods of 140-260 min and 120-220 min for PCM1 and PCM2, respectively. This behavior of the PCMs as energy storage materials can be attributed to the difference in the heat capacity based on various latent heats of fusion for these PCMs. Complete melting occurred during the charging cycle occurred at temperatures in the ranges of (57-60 • C) and (38-43 • C) for PCM1 and PCM2, respectively, which are identical to the reported PCMs meltingsolidification results. A similar trend in the solidification temperatures of the PCM1 and PCM2 can be observed in Figure 6, which shows the variations in the temperatures of PCM1 and PCM2 with time during the discharging cycle at ambient air temperature of 30 • C. A steeper trend for the reduction in the temperatures can be seen before the discharging times of 80 min and 160 min for PCM1 and PCM2, respectively, due to the decrease in sensible heat of PCMs occurring at these discharge periods.  The complete melting and solidification of PCM1 and PCM2 during the charging and discharging cycles are illustrated in Figure 7. Variations in PCMs and room model air   The complete melting and solidification of PCM1 and PCM2 during the charging and discharging cycles are illustrated in Figure 7. Variations in PCMs and room model air  The complete melting and solidification of PCM1 and PCM2 during the charging and discharging cycles are illustrated in Figure 7. Variations in PCMs and room model air temperatures with the time during the discharging cycle are depicted in Figure 8. The results have revealed steady values of PCMs, room inlet and exit air temperatures for the period from 120 to 240 min. A significant increase of about 26% can be seen in room inlet and exit air temperatures compared with the ambient temperature of 30 • C due to the PCMs energy saving.
Fluids 2021, 6, x FOR PEER REVIEW 9 of 15 temperatures with the time during the discharging cycle are depicted in Figure 8. The results have revealed steady values of PCMs, room inlet and exit air temperatures for the period from 120 to 240 min. A significant increase of about 26% can be seen in room inlet and exit air temperatures compared with the ambient temperature of 30 °C due to the PCMs energy saving.        Figure 9 shows the variations in the room model air temperature with time for four cases: without the use of PCMs, with PCM1 alone, with PCM2 alone and with both PCM1 and PCM2 during the discharge cycle. In the current study, the set point air temperature of the room model was set to be in the range of 35-40 • C, which is suitable for air heating applications in the air conditioning systems and food products indirect solar dryers. It can be observed that, in the case of using a single phase change material, PCM1 or PCM2, the room air temperature was higher than that without PCMs but with undesired fluctuations in the temperature range between 23 • C and 47 • C. In the case of using both PCM1 and PCM2 during a relatively low heat energy supply in the discharge cycle, PCM1 provided heat to increase the air temperature in the duct to about 47 • C. The second material PCM2 worked to control the room temperature by minimizing the fluctuations in the air temperature and stabilizing its value in the range of 35-45.2 • C. The room's average temperature was about 39.2 • C, which is desired in the air heating application. This behavior is attributed to the effect of PCM2, which works as excess heat energy storage when the room air temperature is higher than 40 • C and releases energy below the set room temperature (35-40 • C). The augmentation in the discharging time during energy release and improvement in the room air temperature when using multiple PCMs were about 28.4% higher than those without the use of PCMs.
Fluids 2021, 6, x FOR PEER REVIEW 10 of 15 Figure 9 shows the variations in the room model air temperature with time for four cases: without the use of PCMs, with PCM1 alone, with PCM2 alone and with both PCM1 and PCM2 during the discharge cycle. In the current study, the set point air temperature of the room model was set to be in the range of 35-40 °C, which is suitable for air heating applications in the air conditioning systems and food products indirect solar dryers. It can be observed that, in the case of using a single phase change material, PCM1 or PCM2, the room air temperature was higher than that without PCMs but with undesired fluctuations in the temperature range between 23 °C and 47 °C. In the case of using both PCM1 and PCM2 during a relatively low heat energy supply in the discharge cycle, PCM1 provided heat to increase the air temperature in the duct to about 47 °C. The second material PCM2 worked to control the room temperature by minimizing the fluctuations in the air temperature and stabilizing its value in the range of 35-45.2 °C. The room's average temperature was about 39.2 °C, which is desired in the air heating application. This behavior is attributed to the effect of PCM2, which works as excess heat energy storage when the room air temperature is higher than 40 °C and releases energy below the set room temperature (35-40 °C). The augmentation in the discharging time during energy release and improvement in the room air temperature when using multiple PCMs were about 28.4% higher than those without the use of PCMs. The variations in the energy saving of PCM1 and PCM2 with the charging time are illustrated in Figure 10. The energy saving represents the excess heat energy stored by PCMs during the melting process in the charging cycle. It can be observed that PCM1 started to store the excess energy after a charging time of 130 min, reaching about 250 kJ. The energy saving beyond this time was maintained at a steady value due to the influence of PCM1 latent heat of fusion during the phase change period. Similar behavior could be seen for PCM2 after the charging time of 110 min with steady energy savings of 200 kJ. The total energy saving for PCM1 and PCM2 was about 450 kJ during a charging time of 240 min. Figure 11 shows the variations in the energy saving of PCM1 and PCM2 with the The variations in the energy saving of PCM1 and PCM2 with the charging time are illustrated in Figure 10. The energy saving represents the excess heat energy stored by PCMs during the melting process in the charging cycle. It can be observed that PCM1 started to store the excess energy after a charging time of 130 min, reaching about 250 kJ. The energy saving beyond this time was maintained at a steady value due to the influence of PCM1 latent heat of fusion during the phase change period. Similar behavior could be seen for PCM2 after the charging time of 110 min with steady energy savings of 200 kJ. The total energy saving for PCM1 and PCM2 was about 450 kJ during a charging time of 240 min. Figure 11 shows the variations in the energy saving of PCM1 and PCM2 with the discharging time. During the discharging (solidification) cycle, the energy stored in the PCMs was released to increase the air temperature in the duct and room model; thus, there was a continued reduction in the energy saved by PCMs with the discharging time. The total energy saving in the case of using multiple PCMs was higher by about 29.5% and 46.7% than that achieved when using PCM1 and PCM2 alone, respectively. It can be concluded from the obtained results that multiple PCMs revealed higher energy saving and thermal stability for the considered air heating system. The discharging time required to provide the energy and maintain the room model temperature at the desirable value (35-40 • C) in the case of multiple PCMs was extended to 330 min compared with 240 min and 200 min for PCM1 and PCM2, respectively. discharging time. During the discharging (solidification) cycle, the energy stored in the PCMs was released to increase the air temperature in the duct and room model; thus, there was a continued reduction in the energy saved by PCMs with the discharging time. The total energy saving in the case of using multiple PCMs was higher by about 29.5% and 46.7% than that achieved when using PCM1 and PCM2 alone, respectively. It can be concluded from the obtained results that multiple PCMs revealed higher energy saving and thermal stability for the considered air heating system. The discharging time required to provide the energy and maintain the room model temperature at the desirable value (35-40 °C) in the case of multiple PCMs was extended to 330 min compared with 240 min and 200 min for PCM1 and PCM2, respectively.    discharging time. During the discharging (solidification) cycle, the energy stored in the PCMs was released to increase the air temperature in the duct and room model; thus, there was a continued reduction in the energy saved by PCMs with the discharging time. The total energy saving in the case of using multiple PCMs was higher by about 29.5% and 46.7% than that achieved when using PCM1 and PCM2 alone, respectively. It can be concluded from the obtained results that multiple PCMs revealed higher energy saving and thermal stability for the considered air heating system. The discharging time required to provide the energy and maintain the room model temperature at the desirable value (35-40 °C) in the case of multiple PCMs was extended to 330 min compared with 240 min and 200 min for PCM1 and PCM2, respectively.  . Energy saving of PCM1 and PCM2 during the discharging cycle. Figure 11. Energy saving of PCM1 and PCM2 during the discharging cycle.
The present study's results were validated with the reported experimental results of Wisam et al. [14] and Andrei et al. [30] during the PCM melting process, as shown in Figure 12. The comparison was based on results of PCM2 (RT-42) with a melting point of (38-43 • C) in the current study and the work of Wisam et al. [14], while the results of PCM-RT35 with a melting point of (31-36 • C) were considered in the comparison with the results of Andrei et al. [30]. A similar trend in the current study's results of temperature variation with charging time can be observed in Figure 12, with average deviations of 3.4% and 17.6% when compared with the results of Wisam et al. [14] and Andrei et al. [30], respectively. The deviations in the results can be attributed to some differences in the compared works' operating conditions. The present study's results were validated with the reported experimental results of Wisam et al. [14] and Andrei et al. [30] during the PCM melting process, as shown in Figure 12. The comparison was based on results of PCM2 (RT-42) with a melting point of (38-43 °C) in the current study and the work of Wisam et al. [14], while the results of PCM-RT35 with a melting point of (31-36 °C) were considered in the comparison with the results of Andrei et al. [30]. A similar trend in the current study's results of temperature variation with charging time can be observed in Figure 12, with average deviations of 3.4% and 17.6% when compared with the results of Wisam et al. [14] and Andrei et al. [30], respectively. The deviations in the results can be attributed to some differences in the compared works' operating conditions.

Conclusions
The characteristics and energy savings of multiple phase change materials subjected to internal flow in an air heating system during the charging and discharging cycles are investigated in the present study. The experimental work was conducted using a test rig consisting of two main parts, an air supply duct and a room model equipped with PCMs panels. Investigation of the results was based on three cases using single phase change materials PCM1 and PCM2 in the first and second cases, respectively, while in the third case, both PCM1 and PCM2 were used. The obtained results demonstrate a significant improvement in the energy saving and room model temperature control for the air heating system incorporating two PCMs compared with that of a single PCM. The complete melting during the charging cycle occurred at temperatures in the ranges of 57-60 °C and 38-43 °C for PCM1 and PCM2, respectively, which are identical to the reported experimental results. Using multiple PCMs maintained the room's air temperature at a steady value by minimizing the air temperature fluctuations and stabilizing its value in the range of 35-45.2 °C, which is desirable in air heating applications. The total energy saving in the case of using two PCMs was higher by about 29.5% and 46.7% than that of PCM1 and PCM2, respectively. The augmentation in the discharging time and improvement in the room model air temperature when using two PCMs were about 28.4% higher than those without the use of PCMs.

Conclusions
The characteristics and energy savings of multiple phase change materials subjected to internal flow in an air heating system during the charging and discharging cycles are investigated in the present study. The experimental work was conducted using a test rig consisting of two main parts, an air supply duct and a room model equipped with PCMs panels. Investigation of the results was based on three cases using single phase change materials PCM1 and PCM2 in the first and second cases, respectively, while in the third case, both PCM1 and PCM2 were used. The obtained results demonstrate a significant improvement in the energy saving and room model temperature control for the air heating system incorporating two PCMs compared with that of a single PCM. The complete melting during the charging cycle occurred at temperatures in the ranges of 57-60 • C and 38-43 • C for PCM1 and PCM2, respectively, which are identical to the reported experimental results. Using multiple PCMs maintained the room's air temperature at a steady value by minimizing the air temperature fluctuations and stabilizing its value in the range of 35-45.2 • C, which is desirable in air heating applications. The total energy saving in the case of using two PCMs was higher by about 29.5% and 46.7% than that of PCM1 and PCM2, respectively. The augmentation in the discharging time and improvement in the room model air temperature when using two PCMs were about 28.4% higher than those without the use of PCMs.