Numerical Study of the Thermal Performance of Three Roof Models in Hot and Dry Climates

The thermal performance of three roofing models: tile, corrugated and earth terrace is numerically analyzed. The mathematical equations which govern the three roofing models are established by the electrical method of analogies. These equations are discretized by an implicit finite difference method and solved by the Gauss-Seidel algorithm. We analyze the influences of geometric parameters (Xlo, Xlarg, α and Ep) on the evolution of the temperatures of the different environments of our three roof models. In particular, we have shown that the effectiveness of a roof in reducing the temperature inside a room is linked to its physical properties. The results obtained that for the same geometric parameters, the earth roof terrace and the earth tile roof compared to the corrugated metal roof improve thermal comfort by lowering the interior temperature of 5 ̊C and 4.6 ̊C.


Introduction
Burkina Faso, a country with low energy resources, is characterized by a hot and dry climate that favours the transfer of heat inside the habitat, hence the thermal discomfort. The rate of sweating exceeds 80% during the hot periods of the year.
To improve thermal comfort, we proceed by active air treatment methods, such as air conditioners. These methods generate excessive energy consumption. Faced with this, it is imperative to develop other methods of ventilation by natural convection. To do this, several possibilities can be considered, notably by taking into account in the design of the roof envelope. Thus, the search for thermal comfort through the improvement of the housing envelope is a very important area of research and has undergone great development in recent years. Several numerical and experimental studies have been carried out on the search for thermal comfort by improving the roof envelope. Among these studies, [1] studied a passive cooling system using a radio-evaporative roof in a hot and arid climate. It has shown that this type of roof reduces the thermal loads of the habitat by 75% during warm periods. [2] through simulations of the thermal behavior of the diode roof has shown that it reduces the interior temperature of the habitat and generates natural ventilation. A modelling of a bioclimatic roof for passive airconditioning of a typical habitat in Burkina has shown that a north-facing roof south inclined at an angle of 45˚ from the horizontal, gives an average indoor air temperature of 31.3˚C during the warmest hours. At an angle of 55˚, the temperature decreases to 29˚C [3]. The experimental study of the thermal rehabilitation of a room in an Arid zone-case of GHARDIA. [4] showed that by placing 8 cm thick polystyrene glue to a 7 cm thick plasterboard on the roof of the building. Room and west walls exposed to the sun and north exposed to the wind reduces the effect of solar radiation on the studio. The author also shows the decrease in temperature between 8˚C and 10˚C for periods of heat. [5] through a theoretical and experimental study of heat transfers in a round hut with a geobeton roof in the city of Abidjan, show that for a closed hut, geobeton improves comfort by lowering the ambient temperature by 0.5˚C to 2˚C. Through an evaluation study of the thermal comfort of public buildings: case of the architecture department of Tamda (Tizi-Ouzou), [6] shows that the interior temperature of the workshops depends on the orientation and the degree of exposure of facades to solar radiation. [7] studied the impact of thermal inertia on hygrothermal comfort and building energy consumption: the case of the colonial era dwelling in Guelma. He shows through this study that the roof alone encompasses 30% of total heat loss. A study of the influence of local insulating roofing materials on the air conditioning loads of a clay-straw detached house was carried out by [8].
The results obtained by these authors show that the clay-straw mixture reduces the air conditioning requirements compared to the clay wall by around 8%. They also show that a 1.5 cm thick insulation for the roof saves around 8.3% on climatic loads. The indoor temperature levels achieved in an un-air conditioned building depend mainly on the outdoor temperature, internal heat gains, solar gains, ventilation and thermal inertia [9] [10] [11]. In this work, we digitally study the thermal behavior of three roof models: tile, sheet metal and roof terrace on the thermal comfort of the habitat in hot and dry climates. At the end of this study, a roof model suitable for countries with hot and dry climates will be proposed in order to improve thermal comfort in buildings.

Description of Roof Models
The tile roof has been used in Burkina Faso for about ten years. It is used in modern constructions, of size 0.3 m × 0.15 m. The interior consists of a 0.005 -0.010 m thick wood ceiling that reduces heat transfer by conduction into the habitat as shown in Figure 1(a). The corrugated iron roof was introduced into the country by French settlers. It is widely used, especially in urban centers. It consists of a sheet of 3 m long and 0.5 m wide with a thickness varying from 0.003 m to 0.015 m. Inside some habitats there is a 0.005 -0.010 m thick wooden ceiling which also reduces heat transfer by conduction into the habitat as shown in Figure 1(b). The earth roof terrace (Figure 1(c)) is a very old technique used in rural buildings. It consists of tree trunks 3 to 4 m long and 5 to 10 m in diameter. It is then covered with a mixture of grass clay and millet fiber. The mixture is then homogenized with the help of pillions and hands. It obeys endogenous techniques which is intended that the roof terrace.

Mathematical Formulation of Equations Governing Heat Transfers in the Three Roof Models
To simulate the thermal behavior of our three roof models, we divide the length (L) of the roof into several dummy sections according to the direction (X). To write the thermal balance of each slice of the length (L), we proceed by the method of analogies that exists between thermal and electrical transfers. Hypothetically, we assume that there is no mass transfer in the three roof models, and the air movement is one-dimensional (X). Also, we assume that the physical properties of air and roofing materials are constant ( Table 1). The application of ohm The inside face Upper side of ceiling

Discretization of the Equations of Our Roof Models
The discretization of Equations (1) to (4) is expressed in the form: External wall of the roof Inside the roof face t  t  ce  ce  ci  ci  ce  ce  ci  i  a  ci  ce   M Cp T  T  T T  hc T T  S  t  Ep   λ   +   ⋅  −  ⋅  =  −  +  −  ∆   1   1   t  t  t  be  ce  i  ci  be  ci  ce  ce  ce  be  ce  ce  ce Air in the roof Upper side of ceiling

Numerical Resolution Method
The Equations (1)-(4) are discretized (5)-(8) using an explicit method of finite differences. This method transforms the transfer equations into a system of algebraic equations which are solved by the Gauss-Seidel method. The coefficients of heat transfer by natural convection and radiation depend on the temperatures of the different environments which are unknown [12]. Also, an iterative calculation is required. An arbitrary value, taken equal to the ambient temperature, is assigned to the temperatures of the different environments for the calculation of the transfer coefficients which allows the resolution of the system of algebraic equations. The obtained temperatures are then compared with the arbitrary values thus formed by the distribution of the temperatures in the different sections of the roof. If the deviation satisfies the desired precision 3 10 err − = , the calculated temperature is displayed otherwise, the arbitrary value is replaced by the calculated temperatures and the calculation is repeated until the satisfaction of the precision.

Numerical Results and Discussions
We present the numerical results of our calculations for the three roof models, in the form of temperature curves in the different roof environments. Thus, Figure  2 shows the temperature curves of the different environments of the roof tile. For our typical day which is from 8 h to 20 h, the curve of the external wall exposed to the sun grows at 10 h from 311.7 K until it reaches its maximum at 13 h or 314.6 K; after this hour, we observe a drop in temperature which reaches 306.4 K at 6 pm. For the same hours, the internal wall curve shows the values 310.9 K, 313.8 K and 305.3 K. As for the air contained in the roof, the temperature also increased from 310.6 K to 313.5 K and decreased to 305.1 K. These differences in temperature between the roof environments may be due to the physical properties of the materials. Figure 3 shows the evolution of the temperatures of the different environments of the corrugated roof during the same typical day. We observe an increase in the temperature of the outer wall from 10 h to 13 h respectively 316.6 K until 319.5 K; after 13 h, we notice a decrease in temperature which reaches 311.6 K at 18 h. We also observe an increase in the temperature curve of the internal wall from 316.2 K at 10 h to reach 319.1 K at 13 h its maximum value, time at which the curve begins to decrease and reaches 311 K at 18 h. We also notice a temperature growth of the confined air in this model of 10 h and 13 h, being 315.7 K and 318.6 K; and a temperature decrease after 13 h until reaching 310.4 K at 18 h. These small temperature differences between the corrugated roof environments are due to the high thermal conductivity of aluminum. Figure 4 illustrates the evolution of the temperatures of the different    Influence of calculation parameters on temperature distributions:  To find the optimal parameters of our three roof models, we vary the following parameters: length(Xlo), width(Xlarg), inclination(α) and thickness of materials(Ep). The curves in Figure 6 and Figure 9 show the temperatures of the various environments of the corrugated iron roof. The temperature of its wall exposed to solar radiation is identical for the two figures of maximum value 319.4 K at 13 h.       At the level of the roof terrace Figure 7 and Figure 10, we note that the temperature of the external wall is identical for the three figures of maximum value 314.6 K at 13 h. This value is similar to that of Figure 5 and Figure 8 of roof tile. The internal wall temperature differs from roof to roof and is 313.4 K for Figure  7, 313.5 K for Figure 10. Examination of the curves of the Figures 5-10 shows that the geometric parameters influence the temperatures of the different environments of our three roofs. The temperatures of the internal walls and the confined air decrease for ever greater parameters.
By comparing the curves of the Figures 5-10, we find that the curves of Figure 6, Figure 9 of the roof corrugated sheet shows the highest temperatures.

Conclusions
The obvious interest of the roof on the reduction of the internal temperature in a room, we have proceeded to the development of mathematical models based on the electric analogy that can predict the temperature of the different roof environments. These models integrate both meteorological data, geometric and physical parameters of the materials of our roof models. Through this study, we have provided an analysis on the thermal behavior of three roof models used in the Burkinabè habitat. A good use of the numerical method allows the thermal characterization of roofing materials.
At the end of this study, we show that the main causes of thermal discomfort in the habitat remain certain roofing materials with high thermal conductivity such as sheet metal, as well as the geometric shape. Thus, the corrugated roof shows a maximum temperature of 318.8 K for a thermal conductivity of Whatever the variations in the geometric parameters of the three roof models, materials with low thermal conductivity are those that generate the lowest temperatures within the habitat. The results show that the compressed earth roof terrace shows a satisfactory thermal behavior. However, the temperatures of the different media of each roof decrease as geometric parameters increase. The effectiveness of a roof to reduce the interior temperature of a room is related to its geometric shape and physical properties. This digital study results in an endogenous solution on heat transfer of three roof models in hot and dry climates.
Our results also show that the choice of roof of a habitat has a significant impact on thermal comfort. It is always considered as the main element of the thermal regulation of the heat exchange between the interior and the exterior. This leads us to say that the roof terrace is the most suitable for the hot and dry climate that is that of Burkina. As a result, this work contributes to the awareness and promotion of local materials.