Thermal Dynamic Behavior in Bi-Zone Habitable Cell with and without Phase Change Materials †

: The thermal dynamic behavior of buildings represents an important aspect of the energy e ﬃ ciency and thermal comfort of the indoor environment. For this, phase change material (PCM) wallboards integrated into building envelopes play an important role in stabilizing the temperature of the human comfort condition. This article provides an assessment of the thermal behavior of a “bi-zone” building cell, which was built based on high-energy performance (HEP) standards and heated by a solar water heater system through a hydronic circuit. The current study is based on studying the dynamic thermal behavior, with and without implantation of PCMs on envelope structure, using a simpliﬁed modeling approach. The evolution of the average air temperature was ﬁrst evaluated as a major indicator of thermal comfort. Then, an evaluation of the thermal behavior’s dynamic proﬁle was carried out in this study, which allowed for the determination of the PCM rate anticipation in the thermal comfort of the building cell.


Introduction
The sustained growth in the energy demand's dependence on fossil fuel and fluctuating prices, as well as environmental constraints, have prompted African countries to explore potential sources of energy savings. One of the promising responses is energy efficiency. In Morocco, since the launch of the national energy strategy in 2009, energy efficiency has become one of the main pillars for the sustainability and energy independence policy of the kingdom.
The building sector is among the most energy-demanding sectors in Morocco, with an energy consumption of up to 33% of the total energy in the country, including 26% for the residential sector and 7% for the tertiary sector. The building sector today has the greatest potential for improving energy efficiency by representing an energy-saving potential of 40% [1]. An important rate is dominated by heating, ventilation, and air conditioning (HVAC) systems. The Moroccan thermal standard was launched in 2015 to reduce the energy bill and environmental impact. Thus, high-energy performance (HEP) in building standards is a key element to achieve the Moroccan thermal regulation of buildings (RTCM) [2]. The thermal performances required in the RTCM differ according to the type of building (residential or tertiary) and according to the climatic zone by location [3].
Besides, several research works have been developed to improve building energy efficiency, the thermal performance of buildings, and energy consumption savings. In this context, several solutions have been presented, among them the incorporation of phase change materials (PCMs) as promising systems with the ability to improve the imbalance between energy supply and demand and provide high-energy storage density using latent heat.
PCMs are substances that release or absorb energy during melting and freezing processes at a constant or semi-constant temperature based on the phase change temperature of the material.
Therefore, one of the uses of variable phase materials is the storage of latent thermal energy, whereby the variable phase materials combined with building materials can increase the thermal breakdown of buildings to contribute to reducing fluctuations in internal temperature and achieving internal thermal comfort, reducing the energy consumption needed for heating and air conditioning.
In the Mediterranean context, Mandilaras et al. [4] thermally characterize a typical two-story family house embedded with PCM gypsum pasteboard panels. They analyze the building walls orientation effect for meteorological conditions and with conducted experimental measurements, the authors conclude that the thermal mass enhancement for early summer and autumn results in a decrease in decrement factor by 30 to 40% and an increase in time lag by approximately 100 min. These two cited factors are key output parameters to evaluate thermal performance inside buildings; they were also determined and discussed in the works of Kharbouch et al. [5] who investigated the thermal performance of PCM-enhanced walls and roofs in northern Moroccan buildings.
Most supplementary studies have also shown that PCM integration into building envelopes has major benefits. Mourid et al. [6,7] have been experimentally concluding that the application of PCMs on the envelope of a cavity living space is an effective technique for improving thermal comfort by storing heating losses, which reduces the electrical energy consumption, compared to a similar cavity without PCMs. Other research works [8,9] revealed experimentally that the optimal site to achieve greater efficiency by integrating PCMs was located closer to the interior surface of the wall. In addition, Li et al. [10] construed that building envelopes containing PMCs manifest different levels of enhancement by reducing peak temperature and temperature fluctuations. Li et al. [11] improved heat transfer reduction in multi-layer walls of conventional buildings by using PCMs.
This work aims to numerically investigate the thermal dynamic behavior of all the components of a whole building partially covered by PCMs and that of one without PCMs. The outcomes of the study may contribute to clearing the way for energy analysis and comparison with standard and conventional buildings in matters of thermal conditioning.
In the presented paper, the first part describes the thermal model and the proposed bi-zone building description, while the second part discusses briefly the developed scheme of numerical dynamic simulation as well as the different adopted scenarios that were shown to evaluate the use of PCMs in the building envelope. Finally, the results obtained from the dynamic simulation are developed to discover the optimal scenario in the proposed building.

Building Description
The building studied is a prototype that was built at the Abdelmalek Essaadi University Faculty of Sciences in Tetouan (Figure 1), located in the climate Z2 zone [2], which is characterized by moderate climate in the Tangier-Tetouan region in the north of Morocco (N 35 • 34 W 5 • 22 ). For Moroccan conditions, the location is relatively rain-laden and is located on a plateau, about 90 m above mean sea level and 10 km south of the Mediterranean Sea. The city is surrounded by the Rif Mountains in the northern and southern directions.
The building has an entire area of 53 m 2 and consists of two rooms, which are connected by a simple wooden door (4 cm thick, 83 cm × 210 cm) as illustrated by Figure 2. The entrance door on the north side of the building has a heat transfer coefficient of 1.627 W/m 2 K; its structural characteristics are the same as those of the internal door. Furthermore, the building has six simple windows (76 cm × 117 cm, 4 mm thick) with aluminum frames. Two are located in the south of the building, while four face westward. The thermo-physical characteristics of the building are summarized in Table 1. The building has an entire area of 53 m 2 and consists of two rooms, which are connected by a simple wooden door (4 cm thick, 83 cm × 210 cm) as illustrated by Figure 2. The entrance door on the north side of the building has a heat transfer coefficient of 1.627 W/m 2 K; its structural characteristics are the same as those of the internal door. Furthermore, the building has six simple windows (76 cm × 117 cm, 4 mm thick) with aluminum frames. Two are located in the south of the building, while four face westward. The thermo-physical characteristics of the building are summarized in Table 1.
In this work, the PCM wall board was proposed to be embedded in the inside part of the eastern and southern building faces, because those faces were assumed to receive the highest amount of daily solar radiation. Furthermore, this study intended to limit the integration of the PCM in the whole building for economic reasons and to investigate firstly the partial PCM integration in building walls ( Figure 2).    The building has an entire area of 53 m 2 and consists of two rooms, which are connected by a simple wooden door (4 cm thick, 83 cm × 210 cm) as illustrated by Figure 2. The entrance door on the north side of the building has a heat transfer coefficient of 1.627 W/m 2 K; its structural characteristics are the same as those of the internal door. Furthermore, the building has six simple windows (76 cm × 117 cm, 4 mm thick) with aluminum frames. Two are located in the south of the building, while four face westward. The thermo-physical characteristics of the building are summarized in Table 1.
In this work, the PCM wall board was proposed to be embedded in the inside part of the eastern and southern building faces, because those faces were assumed to receive the highest amount of daily solar radiation. Furthermore, this study intended to limit the integration of the PCM in the whole building for economic reasons and to investigate firstly the partial PCM integration in building walls (Figure 2).    In this work, the PCM wall board was proposed to be embedded in the inside part of the eastern and southern building faces, because those faces were assumed to receive the highest amount of daily solar radiation. Furthermore, this study intended to limit the integration of the PCM in the whole building for economic reasons and to investigate firstly the partial PCM integration in building walls (Figure 2).

Heat and Mass Transfer Modeling
The chosen mathematical approach to model the thermal behavior in the bi-zone cell was the heat capacity and thermal conductivity variation concerning the PCM temperature method [12].
The heat equation in the solid-state is defined by the Laplace equation and reduced to the developed form described by the equation; this model describes formally the thermal balance in a transient regime without heat term sources for different sub-layers of the system in a wall sample.
where ρ i , C i , k i , and e i are respectively the density, the heat capacity, the thermal conductivity, and the thickness of the layer i, while i + 1 leads to the adjacent layer. It can be noted in this work that the reference for the indoor chosen direction is from indoor to outdoor. The integration of the PCM layer can be modeled by the following equation: where φ subscript denotes the phase change material layer, while φ − 1 and φ + 1 corresponds to the adjacent layers in front (indoor air) and behind the PCM layer, respectively. The approach used in this work consists to extend the wall model with and without PCM pasteboard into full-building bi-zone scale. Figure 3 represents the multi-scale process RC (resistance capacity) equivalent network modeling for a 1D single wall to a mono and then for a full bi-zone dynamic model with and without PCMs. The heat capacity and thermal conductivity variation were developed in the work of Kharbouch et al. [5]; the PCM used in this work was composed of paraffin-wax microcapsule material and its phase change temperature was 21 • C, which is near thermal comfort conditions. Proceedings 2020, 63, 41 5 of 9 Figure 3. The multi-scale process RC equivalent network modeling.

Meteorological Data (Hourly TMY)
The meteorological data used were extracted from the hourly typical meteorological years) (TMY) Meteonorm database where different levels of global solar radiation that rose into the four facades (north, south, east, and west) were checked. Additionally, the ambient temperature and the The heat transfer between the building envelopes and its outdoor and indoor environments is realized by convection and radiation. The thermo-physical parameters are estimated by referring to the literature and physical assumptions. Figure 4 represents 38 thermal points' predictions that are determined by the simplified dynamic model in the case of PCM integration.

Meteorological Data (Hourly TMY)
The meteorological data used were extracted from the hourly typical meteorological years) (TMY) Meteonorm database where different levels of global solar radiation that rose into the four facades (north, south, east, and west) were checked. Additionally, the ambient temperature and the wind speed were also useful in this work [13].

Finite Difference Method
The finite difference method (FDM) is one of the performing numerical simulation methods that were employed to resolve differential equations, by substituting the derivatives in the equation using

Meteorological Data (Hourly TMY)
The meteorological data used were extracted from the hourly typical meteorological years) (TMY) Meteonorm database where different levels of global solar radiation that rose into the four facades (north, south, east, and west) were checked. Additionally, the ambient temperature and the wind speed were also useful in this work [13].

Finite Difference Method
The finite difference method (FDM) is one of the performing numerical simulation methods that were employed to resolve differential equations, by substituting the derivatives in the equation using differential quotients, principally approximating the boundary conditions value problems by appropriate finite-difference operators [14]. FDM consists of approximating the derivatives of the heat equation-employing Taylor expansions and is deducted directly from the definition of the derivative [15].
The governing equation was solved along with the corresponding boundary and initial conditions using a finite difference method. The simulations were performed using MATLAB Software. The thermal transfer and storage processes of PCMs were taken into account by employing a one-dimensional conduction finite-difference solution as the heat balance algorithm. To account for phase change energy, this algorithm gives an implicit finite difference scheme combined with heat capacity and thermal conductivity function.

Results and Discussion
Dynamic simulation was conducted to determine the thermal behavior in the unheated bi-zone habitable cell with and without phase change materials (PCMs). It is worthwhile to know that the modeling approach predicted the temperature dynamics state for about 41 components of the system with PCMs and 38 components for the non-PCM case. Figure 5 represents the temperature dynamic behavior with PCM integration phase for the 10th to the 12th day of January 2012. It can be seen that the PCMs amplify the natural heat sources in the indoor area. The southern side of the building that is most exposed to the sun causes high thermal behavior in indoor surfaces. Figure 5 represents the temperature dynamic behavior with PCM integration phase for the 10th to the 12 th day of January 2012. It can be seen that the PCMs amplify the natural heat sources in the indoor area. The southern side of the building that is most exposed to the sun causes high thermal behavior in indoor surfaces.  Figure 6 represents the temperature dynamic behavior without the PCM integration cases, and it can be seen that the temperature profile is lower than the temperature profile of the PCM case integration.

Thermal Profile in Summer Period
In this section, we show temporal evolutions of the temperature for the cells with and without PCMs. During the trees first day of July in the summer period, it can be seen from Figures 7 and 8 that  Figure 6 represents the temperature dynamic behavior without the PCM integration cases, and it can be seen that the temperature profile is lower than the temperature profile of the PCM case integration. Figure 5 represents the temperature dynamic behavior with PCM integration phase for the 10th to the 12 th day of January 2012. It can be seen that the PCMs amplify the natural heat sources in the indoor area. The southern side of the building that is most exposed to the sun causes high thermal behavior in indoor surfaces.  Figure 6 represents the temperature dynamic behavior without the PCM integration cases, and it can be seen that the temperature profile is lower than the temperature profile of the PCM case integration.

Thermal Profile in Summer Period
In this section, we show temporal evolutions of the temperature for the cells with and without PCMs. During the trees first day of July in the summer period, it can be seen from Figures 7 and 8 that

Thermal Profile in Summer Period
In this section, we show temporal evolutions of the temperature for the cells with and without PCMs. During the trees first day of July in the summer period, it can be seen from Figures 7 and 8 that the PCMs reduce significantly the thermal behavior in the building cells for the first room and has negative impact for the second room in the summer period.
the PCMs reduce significantly the thermal behavior in the building cells for the first room and has negative impact for the second room in the summer period.   the PCMs reduce significantly the thermal behavior in the building cells for the first room and has negative impact for the second room in the summer period.

Conclusions
The aim of this study was to investigate the thermal behavior of PCM integration in a bi-zone building cells. The dynamic simulation was conducted for 31 subcomponents. The represented results show the dynamic thermal profile along the daytime period. The study was limited to the dynamic simulation and the comparison with the partial PCM integrations into the indoor walls. The prospect of this study was to investigate deeply the effect of the parameters that control the mathematical model, and to optimize the PCM integration in the buildings.