3D Modeling of the Thermal Transfer through Precast Buildings Envelopes

: In this paper, a ﬁnite-element-based model is being introduced and developed, using the Cast3m (CEA, Paris, France) simulation tool, to evaluate the thermo-mechanical behavior of a small-scale test bed. In fact, many studies on thermal behavior of cavities have been carried out in literature. However, none of them took into account the co-existence of all thermal phenomena (conduction, convection, internal/external radiation). The work presented in this paper presents a thermo-mechanical model, which aims to combine, in a holistic way, these phenomena. An experimental validation of the thermal model has been ﬁrst carried out using an infrared camera and DS18B20 (Maxim Integrated Products, Dallas, TX, USA) numerical sensors. Results are reported and show the accuracy of the proposed model since both numerical and experimental values of heat transmittance ﬁt together. The main objective is to evaluate heat losses through the walls, by means of heat transmittance calculation, and proposing new functional materials that will help in energy harvesting, as a perspective of this work. As for the mechanical study, it was meant to investigate the distribution of the mechanical stress towards the building envelope submitted to its own weight. Results showed that the stress is uniformly distributed on the lateral walls of the structure as well as on the ﬂoor.


Introduction
Worldwide, the construction and building sector is responsible for more than 40% of total energy consumption (see Figure 1a) and contributes one third of CO 2 emissions [1]. In Morocco, consumption represents about 25%, from which 18% concerns residential buildings while the rest is devoted to the tertiary sector [2]. This expenditure will escalate even more in the coming years given an increasing rate of population increase, as well as the still-improving standards of living [3]. This will not only lead to an increase in the construction rate and in the use of household equipment (e.g., HVAC), which increase electricity consumption, but also to the increase of harmful gas emissions. In order to face the climate change facing the world, reducing energy consumption is becoming an obligation more than a luxurious choice of living. Different countries have adopted laws and construction regulations in order to rationalize construction behavior and register the building sector through energetic sobriety while emphasizing on the comfort of occupants. Indeed, seventy-three countries have their own codes while eight are in the development process. Among these seventy-three, we can find four with voluntary residential codes, twelve with voluntary tertiary codes, fifty-one with imperative tertiary codes and forty-one with obligatory residential codes [4].
Thermal regulation represents a basic document for energy optimization in a building's envelope at the design and conception phases. The appropriate approaches have been then developed to fit each zone according to its corresponding meteorological data [5]. In Other research work focused on the influence of the aspect ratio on the Nusselt number and hence on the heat transfer with respect to the height of the cavity. For example, Yousaf and Usman [18] presented, in their work, an algorithmic study on a natural convection in two-dimensional square cavity with the presence of sinusoidal roughness elements, which are located simultaneously on hot and cold walls. To solve momentum and energy equations, various Rayleigh numbers, varying between 10 3 and 10 6 , were used for a Newtonian fluid with a Prandtl number of 1. Results showed that the presence of these elements has an impact on the hydrodynamic and thermal behavior of the fluid in the cavity. Khatamifar et al. [19] presented numerical simulation results related to combined heat transfer and natural convection flow in a gradually heated square cavity for a large scope of the Rayleigh number (10 5 -10 9 ). In fact, simulations were carried out using the finite volume method following three dimensionless partition thicknesses and positions. Results showed that the Nusselt number increases with the Rayleigh number but decreases with the thickness. Pandey et al., in their paper [20], presented a review of numerical and experimental research works linked to natural convection in enclosures with/out the presence of internal bodies. These latter were taken under different shapes in order to figure out their impact on buoyancy driven-flow among enclosures. The used methods mainly cover spectral element method, finite element approaches, Latice Boltzmann method (LBM), finite volume method and so on. Their paper also discussed the effect of multiple parameters, such as Rayleigh, Grashof, and Prandtl numbers, in order to make the optimal choice of design parameters according to the desired system requirements. Ouakarrouch et al. [21] presented a simulation englobing the conduction heat transfer and natural convection as well as radiative heat through two kinds of alveolar walls used in a recent building's construction (3-and 6-hole numbers). Results showed that heat transfer is increased once radiative transfer is considered into calculations. Further, it was proved that conductive transfer accounts for almost the half of the heat flux in the block of 20 incorporating a thermally conducting vertical partition.
In most of the research work, and to the best of our knowledge, researchers mainly The main idea behind this work is to focus on the building's envelope by investigating the impact of thermo-mechanical aspects on its behavior. In literature, many researchers have carried out lots of work on thermal analyses of building as a means of ensuring thermal comfort. Candau and Piar [11] tracked indoor temperature of a building in response to externally imposed excitations. A comparative study between simulation and experimentation has been held, basing their analysis on, respectively, model identification and spectral analysis. Gargari et al. [12] presented in their work a 'public housing' building type simulation to evaluate energy savings generated by the change of insulation materials for external walls while adopting a green roof solution. Results showed that for a building with insulation that responds to restrictive energy regulation, moderate benefits are received from a green roof's integration. Houda et al. [13] worked on improving the evaluation method of thermal comfort in office buildings through the arid zones with dry and hot climates, by the analysis of the external intervening parameters. It was shown that predictive comfort and perceived comfort are different, since the first can be calculated based on bioclimatic parameters; as for the second, it varies according to occupants.
Not only that, but other researchers have also investigated cavities with various geometries and a variety of construction materials. For instance, Markatos and Perideous [14] refined a calculation method essentially intended to procure solutions of heat transfer, turbulent flow and buoyancy-driven laminar flow in a square cavity with differentially heated side walls. They also took into consideration a refined mesh for fairly high Rayleigh numbers. The outcomes of their research mainly concern rectangular cavities with Rayleigh numbers going up to 10 6 , a Prandtl number of 0.71 and an aspect ratio of 1. Nonetheless, the method is considered to be more general and can be applied to any kind of geometry. Another interesting work has been performed by Cordoba et al. [15]. In this, the authors conducted a numerical study on a laminar, steady, incompressible, free convection with surface radiation in a two-dimensional open cavity. A heating source was imposed on the wall in parallel to the opening, while the rest of the walls were considered adiabatic. The governing equations of, respectively, energy, mass and momentum were resolved based on the usage of a finite volume approach, which was implemented in FORTRAN. Their study revolved around heat transfer characteristics, the orientation of the structure, and the influence of surface radiation on the fluid flow. Reported results showed that the heat transfer decreases when decreasing the cavity tilt angle for all used values of Rayleigh number and emissivity. They also mentioned that, at a certain tilt angle value, thermal radiation exchange between walls is much more significant. De Vahl [16] investigated the bi-dimensional circulation of air in an enclosed square cavity. Desirable movement was generated using a gradient of temperature imposed at both left and right walls as boundary conditions. Accordingly, De Vahl was able to prove that, in terms of vorticity transport and energy equations for Rayleigh numbers, up to 2 × 10 −5 , can converge. Additionally, he also reported the impact of the Prandtl number, which can play the role of a stabilizing factor, especially if it ranges from 10 −1 to 10 3 . Karatas and Derbently [17], in their study, have considered conjugated radiation and natural convection in a rectangular-shaped closed cavity with only one active vertical wall. The cooled surface is made of aluminum and plastic for the rest of the structure. However, since all walls are painted in white, the emissivity is considered to be constant.
Other research work focused on the influence of the aspect ratio on the Nusselt number and hence on the heat transfer with respect to the height of the cavity. For example, Yousaf and Usman [18] presented, in their work, an algorithmic study on a natural convection in two-dimensional square cavity with the presence of sinusoidal roughness elements, which are located simultaneously on hot and cold walls. To solve momentum and energy equations, various Rayleigh numbers, varying between 10 3 and 10 6 , were used for a Newtonian fluid with a Prandtl number of 1. Results showed that the presence of these elements has an impact on the hydrodynamic and thermal behavior of the fluid in the cavity. Khatamifar et al. [19] presented numerical simulation results related to combined heat transfer and natural convection flow in a gradually heated square cavity for a large scope of the Rayleigh number (10 5 -10 9 ). In fact, simulations were carried out using the finite volume method following three dimensionless partition thicknesses and positions. Results showed that the Nusselt number increases with the Rayleigh number but decreases with the thickness. Pandey et al., in their paper [20], presented a review of numerical and experimental research works linked to natural convection in enclosures with/out the presence of internal bodies. These latter were taken under different shapes in order to figure out their impact on buoyancy driven-flow among enclosures. The used methods mainly cover spectral element method, finite element approaches, Latice Boltzmann method (LBM), finite volume method and so on. Their paper also discussed the effect of multiple parameters, such as Rayleigh, Grashof, and Prandtl numbers, in order to make the optimal choice of design parameters according to the desired system requirements. Ouakarrouch et al. [21] presented a simulation englobing the conduction heat transfer and natural convection as well as radiative heat through two kinds of alveolar walls used in a recent building's construction (3-and 6-hole numbers). Results showed that heat transfer is increased once radiative transfer is considered into calculations. Further, it was proved that conductive transfer accounts for almost the half of the heat flux in the block of 20 incorporating a thermally conducting vertical partition.
In most of the research work, and to the best of our knowledge, researchers mainly emphasize characterizing the thermal behaviors of cavities with different shapes. These latter are either heated using sun rays or a simple heating source so that an active wall is always present. They almost never take into account the co-existence of all thermal phenomena (natural convection, conduction, internal/external radiation) to study the thermal behavior of regular cavities. In the work presented in this paper, we shed more light on the thermal characterization of the envelope of a chosen structure using a combination between radiation, conduction among the different materials, and natural convection in a rectangular cavity clogged with air. Not only that, but we also conduct a mechanical study to evaluate stress distribution on the walls of a structure subjected to its own weight. The main objective is to carry out a transient thermo-mechanical characterization of buildings' envelopes. In fact, we are aiming to evaluate heat losses through the walls and especially through potential thermal bridge areas, by means of calculating transmittance coefficients. Once the behavior is known and the simulation models are fitted by means of experimentations in real sitting scenarios, we will be able to introduce (as perspective to this work) new functional materials in order to decrease heat losses and to contribute to energy harvesting.
A small-scale building, called EEBLAB (Touax Maroc, Mohammedia, Morocco) (energy efficient building laboratory), is used as a testbed in order to conduct experimental and simulation evaluations. Its structure is essentially made of galvanized steel for the external walls, roof and flooring base. Expanded polyurethane is injected into the external walls as a means of insulation. Two types of internal insulation are used: chipboard for the floor and an air gap and plaster boards for the roof. Thermal sensors are placed in different areas on the surfaces of the walls in order to obtain real-time temperature values, which are used to evaluate thermal transmittance coefficients. It is also worth to note that some sensors are placed in areas where thermal bridges are expected to appear. These areas were distinguished using a thermal camera [22] (thermography analysis). Based on experimental results, a finite-element-based numerical model was developed. Simulations were directed and the main outputs are presented to show the accuracy and the performance of the developed mechanical and thermal models.
The leading grants of this work compared to our previous published paper [23] are threefold:

•
Extensive experiments have been conducted using thermography analysis to detect thermal bridge placement among the EEBLAB. • A transient numerical model is proposed for buildings' thermal characterization, taking into account its dynamic behavior. • A mechanical numerical model is presented to evaluate the distribution of mechanical stress towards buildings' envelopes in response to the effect of their own weights.
The remainder of this paper is structured as follows. Section 2 presents an overview of the theoretical backgrounds allowing the creation of both thermal and mechanical models. In Section 3, we introduce the materials and methods, including numerical and experimental set-ups, used to develop and validate the numerical models. Results and discussion are presented in Section 4, where a mesh sensibility study is also presented. Section 5 concludes the paper and provides some perspectives on this work.

Thermal Model Equations
The performed thermal model includes conductivity within the materials, natural convection inside and outside the structure as well as both self-surface radiation and external radiation. The thermal transient problem in a certain domain Ω is mainly governed by the heat equation (Equation (1)) and Fourier's law [24] (Equation (2)). This latter describes the transmission caused by the difference in the temperature between two regions of a medium in physical contact; where ρ is the density (kg/m 3 ), C p is the specific heat capacity (J/kg·K), q is the flux density (W/m 2 ), S is the surface of the medium (m 2 ), f is a function which includes all thermal loading sources, k is the conductivity (W/m·K), ∆T is the temperature variation, and n is the normal vector to the surface: To resolve this kind of problem, we usually proceed by using variation methods and the resulting matrix system to solve on the element Ωe is given in the following equations: where: C, K and Q are respectively the capacity (J/kg·K), conductivity (W/m·K) and heat flux (W) matrixes.
The term convection refers to the heat transmission occurring between a certain surface and a fluid in motion when they are at different temperatures. This kind of transfer is usually divided into two categories depending on the nature of the flow. The first one concerns natural convection. It occurs mainly when the flow is induced by a difference in density, which varies with the temperature in the fluid. The second one is forced convection, which occurs when the flow is generated by external means, such as a fan or atmospheric winds.
The calculation of the heat flux, which is induced by natural convection, is estimated in two steps. First, we estimate the non-dimensional numbers of Grashof and Prandtl: where β is the thermal expansion (K −1 ), g is the gravity (m/s 2 ), L c is the characteristic length (m), ρ is the density (kg/m 3 ) and µ is the dynamic viscosity (Pa·s). Following the Grashof number, the correlations in Table 1 allow calculating the Nusselt number as mentioned in Equation (8). Then the heat transfer coefficient is computed by Equation (9): The expression of the heat flux (W) is expressed following Equation (10), where T a is the ambiance (internal or external) temperature (K), T is the temperature of the wall (K), and S represents the area of the wall (m 2 ). In the formulation of the finite element problem adopted in our study, convection is accounted for by adding a convection matrix H, representing the coefficient h (W/m 2 K), to the conductivity matrix, and the addition of a term Q c , representing the term h × T a , to the second member. C· On the other hand, Stefan-Boltzmann's law translates the radiant thermal transfer, where ε is the emissivity of the used material and ϕ is the heat flux density (W/m 2 ).
The condition at the radiation limits is thus considered as an equivalent condition of convection with a nonlinear convective exchange coefficient h = ε·S·T 3 . This factor can be extracted by factorizing by the term (T a − T) in Equation (12).
In addition, thermal transmittance coefficients are estimated on walls, floors, roofs and windows as follows: • Considering a wall separating two different atmospheres, the total thermal resistance is the set of all the layers of materials or air that constitute the wall, in addition to surface-exchange resistances (resulting from the phenomenon of convection) (see Figure 2a). From this resistance, we can obtain the heat transfer coefficient, which represents the amount of heat transmitted through this wall in steady state, per unit time, per unit area and temperature difference of one-degree Kelvin. In other words, this transmittance coefficient represents the opposite of the total thermal resistance of the wall and is given in W/m 2 ·K. • However, if the wall is made of several non-homogeneous materials, a very simplified calculation method can be used. It mainly consists of considering that the total thermal resistance is between upper and lower limits. These limits, for each case, could be computed by dividing the wall into different sections whose layers are homogeneous and then applying the normal rules of calculation.
to the conductivity matrix, and the addition of a term Qc, representing the term h the second member.
.̇+ ( + ). = + On the other hand, Stefan-Boltzmann's law translates the radiant thermal tr where ε is the emissivity of the used material and φ is the heat flux density (W/m 2

= ( 4 − 4 )
The condition at the radiation limits is thus considered as an equivalent cond convection with a nonlinear convective exchange coefficient h = ε·S·T 3 . This factor extracted by factorizing by the term (Ta−T) in Equation (12).
In addition, thermal transmittance coefficients are estimated on walls, floor and windows as follows: • Considering a wall separating two different atmospheres, the total thermal res is the set of all the layers of materials or air that constitute the wall, in add surface-exchange resistances (resulting from the phenomenon of convectio Figure 2a). From this resistance, we can obtain the heat transfer coefficient, represents the amount of heat transmitted through this wall in steady state, p time, per unit area and temperature difference of one-degree Kelvin. In other this transmittance coefficient represents the opposite of the total thermal res of the wall and is given in W/m 2 ·K. • However, if the wall is made of several non-homogeneous materials, a very fied calculation method can be used. It mainly consists of considering that th thermal resistance is between upper and lower limits. These limits, for eac could be computed by dividing the wall into different sections whose layers mogeneous and then applying the normal rules of calculation.
According to Figure 2b, the transmittance coefficient (in W/m 2 ·K) can be subd into two coefficients (upper and lower) and can be calculated as shown in Equatio and (15). According to Figure 2b, the transmittance coefficient (in W/m 2 ·K) can be subdivided into two coefficients (upper and lower) and can be calculated as shown in Equations (14) and (15). U up = P as U as + P bs U bs + P cs U cs (14) where, U up is the upper limit of the transmission heat coefficient. U low is the lower limit of the transmission heat coefficient. P as , P bs , P cs , P 1 , P 2 and P 3 are respectively the percentage of material constituting layer a s , b s , c s , 1, 2 and 3. The total heat transfer coefficient (U T ) is then expressed as follows: In order to obtain more accuracy in terms of calculations, thermal bridges should be taken into account. They essentially designate the heat losses through the interfaces of constructive elements and must be added to the surface losses. They are largely divided into two categories. The first category represents structural thermal bridges, which are related to the techniques used for insulator placement as well as those of constructions in the frame of wood. The second category represents thermal bridges due to geometric constraints, which are linked to the shape of the envelope.
Broadly, the main thermal bridges of any kind of buildings are located at the articulations of the facades and floors, facades and roofs, facades and opening frames. If the walls are highly insulated, thermal bridges persist more than usual and, hence, thermal losses increase. Therefore, in the case of low-consumption buildings, such as small-scale structure we are considering in this study, having high thermal resistances could ensure low heat losses per junction. If these leakages are not treated, they will cause discomfort for occupants, increase energy expenses, and might cause a deterioration of materials as well.
Taking into account thermal bridges' effect, the global thermal transmittance is given in Equation (17): where ψ i L i represents punctual thermal bridges, χ j represents linear thermal bridges and S is the surface of the layer. As with heat losses through the walls, thermal transmittance of the window needs to be considered. The EEBLAB window consists mainly of two parts: the aluminum frame and the glazing. For each of these parts, we must calculate the heat transfer coefficient. It should be noted that the outline of the glazing is affected by additional heat loss resulting from the combined effects of the frames, panels and spacers (constituting the linear thermal bridges). All of these elements make it possible to determine, by computing, the heat transfer coefficient of the window U w [26].
It is possible to generalize this formula, as shown in Equation (19), if we take into account the panel of the window or the ventilation grid. P g is the visible perimeter of the glazing; S g , S f , S p and S vg are respectively the surfaces of the glazing, frame, panel and ventilation grid.
The "building energy performance" regulation provides a simplified formula for evaluating the effectiveness of glazing and framing, while taking into account good air- tightness [26]. The proportion varies depending on whether the glazing is more thermally efficient or the reverse.
If the transmittance coefficient of the glazing is less than the transmittance coefficient of the frame, the total transmittance coefficient of the window is given as in Equation (20): Otherwise, In the case of absence of ventilation grid and panel (which is our case), the formula becomes simple, as follows: In other words, we consider 70% of glazing (since in most cases it is considered to be the most efficient [26], 30% of framing and 3 m of interlayer per m 2 of window.

Mechanical Model Equations
Considering a general problem of linear mechanics in small deformations, governed by the following equations.

•
The equilibrium equation (Equation (23)), which illustrates the static or dynamic equilibrium of all external and internal forces of the studied system [27]: • The behavior law following three-dimensional Hook elasticity as shown in Equation (24). • The hypothesis of small deformations by assuming that the transformation between the initial equilibrium state and the current equilibrium state is infinitesimal [28], i.e., where ε T , σ d , U and C are respectively the total strain tensor (no unit), the stress tensor (N/m 2 ), the displacement field (m) and the elasticity modulus tensor (Pa).

Numerical Set-Up
This section presents the geometry of our testbed building, together with parameters and tools used.

Geometry
The specimen considered in this study is one of the rectangular cavities, shown in Figure 3, based at our university. It is mainly deployed to investigate energy efficiency approaches for buildings.  This small scaled testbed has 12 m 2 of occupied surface and 30 m 3 of volume. Table 2 gives more information about its dimensions. It is essentially made from galvanized steel and expanded polyurethane injected into the walls for insulation purposes. The roof is made of a stratification of 3 layers (inside-out): plaster, air, galvanized steel. However, the floor is insulated using chipboard. The same geometry of the EEBLAB was re-constructed in respect to the real dimensions (see Figure 4), so that precision and reliability of the model can be expected. The model is implemented on the cast3m tool, based on the finite element method, and using an Inspiron 157,000 gaming laptop with a processor of Intel ® Core™i7-7700HQ CPU@ 2.8 GHz, 2801 MHz, 4 cores and 8 logical processors (Dell Technologies, Round Rock, TX, USA). The assigned number of elements is 19,768 with about 41,030 nodes. Using the Cast3m tool, a density of 0.08 has been imposed in order to acquire the previously mentioned number of elements. The chosen type of meshing elements was cubic, since the geometry is regular and has a parallelepiped shape.  This small scaled testbed has 12 m 2 of occupied surface and 30 m 3 of volume. Table 2 gives more information about its dimensions. It is essentially made from galvanized steel and expanded polyurethane injected into the walls for insulation purposes. The roof is made of a stratification of 3 layers (inside-out): plaster, air, galvanized steel. However, the floor is insulated using chipboard. The same geometry of the EEBLAB was re-constructed in respect to the real dimensions (see Figure 4), so that precision and reliability of the model can be expected. The model is implemented on the cast3m tool, based on the finite element method, and using an Inspiron 157,000 gaming laptop with a processor of Intel ® Core™i7-7700HQ CPU@ 2.8 GHz, 2801 MHz, 4 cores and 8 logical processors (Dell Technologies, Round Rock, TX, USA). The assigned number of elements is 19,768 with about 41,030 nodes. Using the Cast3m tool, a density of 0.08 has been imposed in order to acquire the previously mentioned number of elements. The chosen type of meshing elements was cubic, since the geometry is regular and has a parallelepiped shape.  This small scaled testbed has 12 m 2 of occupied surface and 30 m 3 of volume. Table 2 gives more information about its dimensions. It is essentially made from galvanized steel and expanded polyurethane injected into the walls for insulation purposes. The roof is made of a stratification of 3 layers (inside-out): plaster, air, galvanized steel. However, the floor is insulated using chipboard.  The same geometry of the EEBLAB was re-constructed in respect to the real dimensions (see Figure 4), so that precision and reliability of the model can be expected. The model is implemented on the cast3m tool, based on the finite element method, and using an Inspiron 157,000 gaming laptop with a processor of Intel ® Core™i7-7700HQ CPU@ 2.8 GHz, 2801 MHz, 4 cores and 8 logical processors (Dell Technologies, Round Rock, TX, USA). The assigned number of elements is 19,768 with about 41,030 nodes. Using the Cast3m tool, a density of 0.08 has been imposed in order to acquire the previously mentioned number of elements. The chosen type of meshing elements was cubic, since the geometry is regular and has a parallelepiped shape.

Parameters Identification
In order to perform simulations under Cast3m, complementary information about the construction and insulation materials parameters are needed, especially those related to thermal and mechanical types, as shown in Tables 3 and 4.  GS (Galvanized steel), CB (Chipboard), PB (Plasterboard), Al (Aluminium), Gl (Glazing).

Numerical Tools
In order to numerically compute the thermomechanical behavior of the EEBLAB's envelope, the following procedures were considered. First, we start by the creation of a thermal solving procedure considering all the classical thermal phenomena (see Figure 5). Then, the second procedure is used to solve the mechanical equilibrium problem taking into account the own weight of the EEBLAB as the only force.
The main objective of the thermal procedure is to evaluate the transmitted flux among the EEBLAB's walls in every single node. Hence, the calculation of the global heat transfer coefficient, including all thermal phenomena, is crucial. Figure 5 summarizes the steps of the computing procedures starting with the regeneration of the same EEBLAB's geometry in respect to the real dimensions for more accurate results. Thermal parameters are also required to be declared as input data. For this, a manual thermometer was used to measure the needed information (ambient temperatures, walls' temperatures) and exploit them to calculate the heat transfer coefficients using the classical correlations of Grashof and Nusselt (see Equation (8), Table 1).
It is worth to mention that the flow was considered as being laminar. Table 5 gives complementary information about air characteristics, which are needed for computing heat transfer coefficients.  It is worth to mention that the flow was considered as being laminar. Table 5 gives complementary information about air characteristics, which are needed for computing heat transfer coefficients.
(a) Thermal Boundary and Interface Conditions: In this subsection, thermal boundary and interface conditions (see Figure 6) are presented. Starting with boundary conditions (see Figure 6a), internal and external temperatures have been imposed on each facade of the walls as boundary conditions for convective heat transfer. As for considering the external radiative heat transfer, an infinite temperature of the value 500 K has been imposed. This latter also presents the initial temperature for numerical calculations. It is worth noting that internal radiative heat transfer has also been taken into account by calculating form factors.
For interface conditions, the external walls are made from different construction materials (multi-layers). Hence, it was supposed that during the calculation of heat transfer coefficients that layers in contact with each other have the same temperature at the contact interface. In other words, the contact resistance is neglected.
In addition, thermal conduction at the surface of the galvanized steel is equal to thermal conduction at the surface of the insulation layer. In other words, the interface between layers does not store energy (see Figure 6b). In this subsection, thermal boundary and interface conditions (see Figure 6) are presented. Starting with boundary conditions (see Figure 6a), internal and external temperatures have been imposed on each facade of the walls as boundary conditions for convective heat transfer. As for considering the external radiative heat transfer, an infinite temperature of the value 500 K has been imposed. This latter also presents the initial temperature for numerical calculations. It is worth noting that internal radiative heat transfer has also been taken into account by calculating form factors.
(b) Mechanical Boundary Conditions: To launch the mechanical study of the EEBLAB using the Cast3m tool, we have sup- For interface conditions, the external walls are made from different construction materials (multi-layers). Hence, it was supposed that during the calculation of heat transfer coefficients that layers in contact with each other have the same temperature at the contact interface. In other words, the contact resistance is neglected.
In addition, thermal conduction at the surface of the galvanized steel is equal to thermal conduction at the surface of the insulation layer. In other words, the interface between layers does not store energy (see Figure 6b).
(b) Mechanical Boundary Conditions: To launch the mechanical study of the EEBLAB using the Cast3m tool, we have supposed that the displacement following axes x, y and z are equal to zero. The building is subjected to its own weight (see Figure 7).
(b) Mechanical Boundary Conditions: To launch the mechanical study of the EEBLAB using the Cast3m tool, we have supposed that the displacement following axes x, y and z are equal to zero. The building is subjected to its own weight (see Figure 7).   The EEBLab in both mechanical and thermal models is subjected to weather conditions of Sala El Jadida during the simulated periodof time, mainly in March 2019 (see Figure 8). The resulted numerical data are extracted in nodes where the temperature sensors are placed on both external and internal walls. Figure 9 depicts the orthonormal reference, which allowed us to extract the exact coordinates of the location of the sensors.   The resulted numerical data are extracted in nodes where the temperature sensors are placed on both external and internal walls. Figure 9 depicts the orthonormal reference, which allowed us to extract the exact coordinates of the location of the sensors.  It is worth to mention also that the characteristic length needed in the calculation of the Grashof number was taken as the width of the construction (LC = 4 m). The surface, which is needed to calculate the heat flux afterwards, is equal to 16 m 2 (LC 2 ). This is because the EEBLAB's walls have different surfaces, which may complicate the calculation process and takes more simulation time.
For the mechanical model, the main objective is to evaluate the structure resistance to external forces. However, it is important to first evaluate the mechanical resistance of the structure in response to its own weight. In this study, we mainly focused our interest on this part (see Figure 10). The code requires the geometry and all its mechanical characteristics. Functions and procedures are then integrated in order to evaluate the displacement and stress fields. It is worth to mention also that the characteristic length needed in the calculation of the Grashof number was taken as the width of the construction (LC = 4 m). The surface, which is needed to calculate the heat flux afterwards, is equal to 16 m 2 (LC 2 ). This is because the EEBLAB's walls have different surfaces, which may complicate the calculation process and takes more simulation time.
For the mechanical model, the main objective is to evaluate the structure resistance to external forces. However, it is important to first evaluate the mechanical resistance of the structure in response to its own weight. In this study, we mainly focused our interest on this part (see Figure 10). The code requires the geometry and all its mechanical characteristics. Functions and procedures are then integrated in order to evaluate the displacement and stress fields.

Experimental Set-Up
Experimental setup consists of using two essential devices (see Figure 11a). First, an infrared camera was deployed in order to observe temperature distribution among the studied surfaces and to identify thermal bridges' placements as well. The camera was mainly focused towards the front and left facades, since they contain different materials and components (e.g., beams, door, window, etc.). The thermal analysis leads to areas where temperature sensors, specifically DS18B20 (Maxim Integrated Products, Dallas, Texas, United States) should be placed (see Figure 11b).

Experimental Set-Up
Experimental setup consists of using two essential devices (see Figure 11a). First, an infrared camera was deployed in order to observe temperature distribution among the studied surfaces and to identify thermal bridges' placements as well. The camera was mainly focused towards the front and left facades, since they contain different materials and components (e.g., beams, door, window, etc.). The thermal analysis leads to areas where temperature sensors, specifically DS18B20 (Maxim Integrated Products, Dallas, TX, USA) should be placed (see Figure 11b). infrared camera was deployed in order to observe temperature distribution among the studied surfaces and to identify thermal bridges' placements as well. The camera was mainly focused towards the front and left facades, since they contain different materials and components (e.g., beams, door, window, etc.). The thermal analysis leads to areas where temperature sensors, specifically DS18B20 (Maxim Integrated Products, Dallas, Texas, United States) should be placed (see Figure 11b). These sensors provide temperature measurement in degrees Celsius on the order of 9 to 12 bits and are equipped with a user-defined alarm function. It communicates via a one-bus wire and contains a numerical resolution of 12 bits. This kind of sensor integrates all required functionalities to measure temperature: (i) analog sensor, (ii) analog/digital converter, and (iii) communication electronics and power supply [31]. For our case study, the sensors were linked to an embedded card. To perform the assembly, extra equipment was needed, such as a 4.7 K-Ohm resistance (code color yellow-violet-red), a circuit board and cables.
Depending on an already existing internet of things (IoT) framework, as mentioned in the work of NaitMalek et al. [32], the walls' temperature data were regularly acquired during different periods of time. The data were collected and treated using an Intel ® Xeon ® These sensors provide temperature measurement in degrees Celsius on the order of 9 to 12 bits and are equipped with a user-defined alarm function. It communicates via a one-bus wire and contains a numerical resolution of 12 bits. This kind of sensor integrates all required functionalities to measure temperature: (i) analog sensor, (ii) analog/digital converter, and (iii) communication electronics and power supply [31]. For our case study, the sensors were linked to an embedded card. To perform the assembly, extra equipment was needed, such as a 4.7 K-Ohm resistance (code color yellow-violet-red), a circuit board and cables.
Depending on an already existing internet of things (IoT) framework, as mentioned in the work of NaitMalek et al. [32], the walls' temperature data were regularly acquired during different periods of time. The data were collected and treated using an Intel ® Xeon ® CPU E5-2650 0 @ 2GHz computer (Dell Technologies, Round Rock, TX, USA); 64 Gb of RAM, a GPU of 43 Go-NVIDIA, an operating system of 64 bits and a 64 processor.
This collection of data allows analyzing and comparing both experimental and simulation results. Sensors were, symmetrically, placed on external as well as internal sides of the walls, so that the transmittance coefficients could be correctly evaluated. The obtained data were then transmitted, in real time, via Wi-Fi to our IoT platform for processing and storage.

Results and Discussions
In this part, a mesh model and grid independence checking study is presented. Besides, experimental and numerical results are given and confronted to each other in order to ascertain the finite element model developed on CAST3M.

Mesh Model and Grid Independence Checking
In this subsection, we present a meshing sensibility study. As was previously mentioned in Section 3.1.2, the chosen type of meshing element is cubic since the geometry of the EEBLab is a regular parallelepiped. According to the Cast3m programming, the size of the meshes, and hence the number of elements, can be easily changed by using the command 'dens'. Focusing our interest on the mechanical study, we can notice (see Table 6) that the shear stress value on the external walls slightly varies (from −6.58 × 10 −6 Pa to 3.26 × 10 −5 Pa) until reaching a 'dens' value of 0.06 (SMXX = 3.25 × 10 −5 Pa). As known, reducing mesh sizing leads to increasing the number of elements, which might slow the simulation and certainly requires a lot of memory. Hence, for optimization purposes, we chose to stick with the 0.08 value (19,768 elements with about 41,030 nodes). Table 6. Shear values according to meshing density.

Density of Meshing SMXX (Shear Stress in Pa)
0.03 command 'dens'. Focusing our interest on the mechanical study, we can notice (see Table  6) that the shear stress value on the external walls slightly varies (from −6.58 × 10 −6 Pa to 3.26 × 10 −5 Pa) until reaching a 'dens' value of 0.06 (SMXX = 3.25 × 10 −5 Pa). As known, reducing mesh sizing leads to increasing the number of elements, which might slow the simulation and certainly requires a lot of memory. Hence, for optimization purposes, we chose to stick with the 0.08 value (19,768 elements with about 41,030 nodes).   Table 6. Cont.

Density of Meshing SMXX (Shear Stress in Pa)
0.09 0.08 0.09 Figure 12 shows the temperature evolution given by numerical simulations. According to Figure 12a, the temperature on both the front and right wall evolves exactly in the Figure 12 shows the temperature evolution given by numerical simulations. According to Figure 12a, the temperature on both the front and right wall evolves exactly in the same way, where it slightly decreases from a maximum of 293 K to reach a minimum of 292.8 K. For the sensor placed on the left wall, the temperature evolution may be considered as negligible since the temperature is almost stable at the value of 293 K. Concerning the roof, the temperature evolution behaves differently since it starts from a minimum value of 293 K to reach 293.15 K. All in all, internal temperature variation might be considered as insignificant since the difference between the minimum and the maximum is no larger than 0.3 K. same way, where it slightly decreases from a maximum of 293 K to reach a minimum of 292.8 K. For the sensor placed on the left wall, the temperature evolution may be considered as negligible since the temperature is almost stable at the value of 293 K. Concerning the roof, the temperature evolution behaves differently since it starts from a minimum value of 293 K to reach 293.15 K. All in all, internal temperature variation might be considered as insignificant since the difference between the minimum and the maximum is no larger than 0.3 K. Regarding Figure 12b, the external temperature evolves in the same way for the roof, front and left walls. It varies between a minimum of 293 K and a maximum of 293.4 K; however, the right wall temperature is observed to be stagnant at a value of 293 K. Figure 13 shows the numerical estimated evolution of the thermal transmittance coefficients towards the different walls. According to this figure, this parameter generally alters between a minimum value of 0.04 W/m 2 K and a maximum value of 0.06 W/m2 K. The evolution curves have almost the same qualitative and quantitative aspects. Regarding Figure 12b, the external temperature evolves in the same way for the roof, front and left walls. It varies between a minimum of 293 K and a maximum of 293.4 K; however, the right wall temperature is observed to be stagnant at a value of 293 K. Figure 13 shows the numerical estimated evolution of the thermal transmittance coefficients towards the different walls. According to this figure, this parameter generally alters between a minimum value of 0.04 W/m 2 K and a maximum value of 0.06 W/m 2 K. The evolution curves have almost the same qualitative and quantitative aspects.

Thermal Numerical Results
Regarding Figure 12b, the external temperature evolves in the same way for the roof, front and left walls. It varies between a minimum of 293 K and a maximum of 293.4 K; however, the right wall temperature is observed to be stagnant at a value of 293 K. Figure 13 shows the numerical estimated evolution of the thermal transmittance coefficients towards the different walls. According to this figure, this parameter generally alters between a minimum value of 0.04 W/m 2 K and a maximum value of 0.06 W/m2 K. The evolution curves have almost the same qualitative and quantitative aspects.

Thermal Experimental Results
As previously mentioned, the camera was only used towards the front and left walls because they are the only sides easily accessed, as well as because they contain different components. Regarding the front side, the camera was fixed at 28 m away to have a global view, but for the left side, it was placed 4 m away due to infrastructure constraints. In order to obtain reliable results, the setting parameters, such as values of emissivity, distance and ambient temperature were initialized accounting for the experiment period as well as the studied structure. Horizontal lines were plotted on both sides (see Figure 14),

Thermal Experimental Results
As previously mentioned, the camera was only used towards the front and left walls because they are the only sides easily accessed, as well as because they contain different components. Regarding the front side, the camera was fixed at 28 m away to have a global view, but for the left side, it was placed 4 m away due to infrastructure constraints. In order to obtain reliable results, the setting parameters, such as values of emissivity, distance and ambient temperature were initialized accounting for the experiment period as well as the studied structure. Horizontal lines were plotted on both sides (see Figure 14), so that the temperature distribution can be observed and then the influence of the co-existence of different component on enhancing thermal bridges can be evidenced. Figure 14 shows the temperature distribution on both front and left walls, using thermography analysis. According to Figure 15a, at the level of galvanized steel areas, the temperature varies between 301.5 K and 307 K. However, once Line 1 crosses the window and the door, temperature falls down to reach 297 K and 301 K, respectively. This is due to the variation of the thermal conductivity due to the materials' heterogeneity, which allows the air elusion through the material itself. In addition, in Figure 15b, it is observed that on the level of steel areas, the temperature varies between 298 K and 301 K. A significant offset is noticed when reaching the beam at the middle position and areas labeled IC. This is due to the imperfect contact between the wall and the beam as well as between the combined parts of steel sheets composing the wall. so that the temperature distribution can be observed and then the influence of the coexistence of different component on enhancing thermal bridges can be evidenced. Figure  14 shows the temperature distribution on both front and left walls, using thermography analysis. According to Figure 14a, at the level of galvanized steel areas, the temperature varies between 301.5 K and 307 K. However, once Line 1 crosses the window and the door, temperature falls down to reach 297 K and 301 K, respectively. This is due to the variation of the thermal conductivity due to the materials' heterogeneity, which allows the air elusion through the material itself. In addition, in Figure 14b, it is observed that on the level of steel areas, the temperature varies between 298 K and 301 K. A significant offset is noticed when reaching the beam at the middle position and areas labeled IC. This is due to the imperfect contact between the wall and the beam as well as between the combined parts of steel sheets composing the wall. All in all, thermal bridges mainly persist on areas where imperfect contacts exist in addition to ones where different materials and components coexist. Figure 15 shows the temperature evolution given by the deployed sensors, which are installed symmetrically on each side. Measurements were conducted on the 21st of March 2019 between 10.53 a.m. and 13.44 p.m., with 10 min time step. All in all, thermal bridges mainly persist on areas where imperfect contacts exist in addition to ones where different materials and components coexist. Figure 16 shows the temperature evolution given by the deployed sensors, which are installed symmetrically on each side. Measurements were conducted on the 21st of March 2019 between 10.53 a.m. and 13.44 p.m., with 10 min time step. According to Figure 16a, the temperature on both roof and right wall evolves in the same way quantitatively as well as qualitatively, where it slowly increases from 290 K to reach 294 K. However, for the sensor placed on the left wall, the temperature increases from 290 K to attend 305 K during the same period of time. For the front wall, its evolution behaves differently since it starts from high values (302 K-310 K) between 11 a.m. and All in all, thermal bridges mainly persist on areas where imperfect contacts exist in addition to ones where different materials and components coexist. Figure 16 shows the temperature evolution given by the deployed sensors, which are installed symmetrically on each side. Measurements were conducted on the 21st of March 2019 between 10.53 a.m. and 13.44 p.m., with 10 min time step. 12.30 a.m., and then decreases to reach 295 K afterwards. This can be explained by the fact that this specific wall was directly exposed to sun rays during this period of time.  Figure 16b presents the temperature evolution given by internal temperature sensors. The temperature on each of the four walls evolves in the same manner and starts from a mean value of 293 K to reach 297 K during, approximately, a 3 h period of time. Figure 17 shows the evolution of the thermal transmittance coefficient towards the different walls. It can be noticed that this parameter generally alters between a minimum of 0.05 W/m 2 ·K. and a maximum of 0.06 W/m 2 ·K. It is worth to mention that the transmittance coefficients of all surfaces evolve quantitatively as well as qualitatively in the same way, except for some negligible off-sets noticed at the right wall and roof. According to Figure 16a, the temperature on both roof and right wall evolves in the same way quantitatively as well as qualitatively, where it slowly increases from 290 K to reach 294 K. However, for the sensor placed on the left wall, the temperature increases from 290 K to attend 305 K during the same period of time. For the front wall, its evolution behaves differently since it starts from high values (302 K-310 K) between 11 a.m. and 12.30 a.m., and then decreases to reach 295 K afterwards. This can be explained by the fact that this specific wall was directly exposed to sun rays during this period of time. Figure 16b presents the temperature evolution given by internal temperature sensors. The temperature on each of the four walls evolves in the same manner and starts from a mean value of 293 K to reach 297 K during, approximately, a 3 h period of time. Figure 17 shows the evolution of the thermal transmittance coefficient towards the different walls. It can be noticed that this parameter generally alters between a minimum of 0.05 W/m 2 ·K. and a maximum of 0.06 W/m 2 ·K. It is worth to mention that the transmittance coefficients of all surfaces evolve quantitatively as well as qualitatively in the same way, except for some negligible off-sets noticed at the right wall and roof. The temperature on each of the four walls evolves in the same manner and starts from a mean value of 293 K to reach 297 K during, approximately, a 3 h period of time. Figure 17 shows the evolution of the thermal transmittance coefficient towards the different walls. It can be noticed that this parameter generally alters between a minimum of 0.05 W/m 2 ·K. and a maximum of 0.06 W/m 2 ·K. It is worth to mention that the transmittance coefficients of all surfaces evolve quantitatively as well as qualitatively in the same way, except for some negligible off-sets noticed at the right wall and roof.  Figure 18 presents both numerical and experimental results of transmittance coefficients throughout the different walls. It is noticed that all values alter between a minimum of 0.045 W/m 2 K and a maximum of 0.06 W/m 2 K. Results are shown to be in agreement, with some discrepancies in the order of 1%, 7%, 6% and 8% on the front wall, right wall, left wall and the roof, respectively. This can be explained by the fact that the problem is non-linear and needs more calculation precision (number of iterations, time step, meshing size, etc.) to obtain more accurate results, not to forget the disturbances that numerical sensors may have experienced.  Figure 18 presents both numerical and experimental results of transmittance coefficients throughout the different walls. It is noticed that all values alter between a minimum of 0.045 W/m 2 K and a maximum of 0.06 W/m 2 K. Results are shown to be in agreement, with some discrepancies in the order of 1%, 7%, 6% and 8% on the front wall, right wall, left wall and the roof, respectively. This can be explained by the fact that the problem is non-linear and needs more calculation precision (number of iterations, time step, meshing size, etc.) to obtain more accurate results, not to forget the disturbances that numerical sensors may have experienced. In summary, the thermography analysis allowed us to estimate the temperature distribution on the front and left surfaces as well as to detect thermal bridges through building areas. It was noticed that thermal bridges mainly exist on areas where there is coexistence of different materials making the contact non-perfect. When different materials coexist on the same wall, they lead to heat exhaust since the conductivity changes. Data acquired from temperature sensors allows computing the evolution of the temperature on both internal and external surfaces. According to Figure 16a, the temperature evolution on the front and left walls is more significant compared to the rest of the walls. This is due to the fact that sun rays directly reach these walls during this period of the day. Results In summary, the thermography analysis allowed us to estimate the temperature distribution on the front and left surfaces as well as to detect thermal bridges through building areas. It was noticed that thermal bridges mainly exist on areas where there is co-existence of different materials making the contact non-perfect. When different materials coexist on the same wall, they lead to heat exhaust since the conductivity changes. Data acquired from temperature sensors allows computing the evolution of the temperature on both internal and external surfaces. According to Figure 16a, the temperature evolution on the front and left walls is more significant compared to the rest of the walls. This is due to the fact that sun rays directly reach these walls during this period of the day. Results depicted in Figure 16b show that the temperature evolves mainly in the same way. This is due to the fact that the structure is totally closed, blocking infiltrations and air circulation.

Thermal Results Confrontation
Regarding the experimental thermal transmittance coefficients, Figure 17 shows that these parameters evolve similarly on all of the walls. This is due to the fact that they are mainly made from the same materials with, more or less, the same proportions of thickness.
Indeed, 17 values of the ambient air temperature have been computed, between 10.53 a.m. and 13.44 p.m., and included into the simulation code as input values. Figure 18 confronts simulations and experimental results. Thermal transmittance coefficients turned out to have the same values, which prove the efficiency of our numerical thermal model. Despite that some discrepancies have been noticed, this can be avoided by improving the calculation precision (increasing the number of iterations, optimizing the meshing size, etc.). The small values of the thermal transmittance coefficients are due to the injection of lateral walls with polyurethane. This latter is considered to be a very performant insulation material (k = 0.024 W/m·K), hence, high resistance of the walls.

Mechanical Numerical Results
The implemented code is also used to solve the mechanical equilibrium problem and displays the variation of the building displacement and mechanical stress distribution. Since the structure is stationary, the displacement is considered to be negligible. This can be seen in Figure 19, where the displacement value is uniform through the envelope and is approximatively equal to −1.13 × 10 −7 m. The finite element mechanical model allows us to determine the distribution of the mechanical stress within the building's envelope when it is submitted to its own weight. As shown in Figure 19, the displacement of the overall structure is uniform and very small, so it can be neglected. Regarding the stress distribution, Figure 20 shows that the values are uniform on the lateral walls and do not exceed −1.39 × 10 4 Pa and −5.74 × 10 2 Pa for normal and shear stress, respectively. However, a stress concentration is evidenced through the window and its frame since the construction materials are different from the rest of the envelope. The finite element mechanical model allows us to determine the distribution of the mechanical stress within the building's envelope when it is submitted to its own weight. As shown in Figure 19, the displacement of the overall structure is uniform and very small, so it can be neglected. Regarding the stress distribution, Figure 20 shows that the values are uniform on the lateral walls and do not exceed −1.39 × 10 4 Pa and −5.74 × 10 2 Pa for normal and shear stress, respectively. However, a stress concentration is evidenced through the window and its frame since the construction materials are different from the rest of the envelope. Figure 19, the displacement of the overall structure is uniform and very small, so it can be neglected. Regarding the stress distribution, Figure 20 shows that the values are uniform on the lateral walls and do not exceed −1.39 × 10 4 Pa and −5.74 × 10 2 Pa for normal and shear stress, respectively. However, a stress concentration is evidenced through the window and its frame since the construction materials are different from the rest of the envelope.

Conclusions and Perspectives
In this paper, a numerical finite-element-based model was proposed to investigate the thermo-mechanical behavior of buildings' envelopes. A small building has been chosen as a test specimen for simulation, experimental analysis and evaluation. An infrared camera has been employed to visualize temperature distribution on the external walls and to detect thermal bridges positions as well. Numerical and experimental thermal results have been confronted to evaluate the built-up model. An agreement between results has been observed for the thermal transmittance coefficients validating the numerical model. Some minor differences between experiments and simulation results have been noticed,

Conclusions and Perspectives
In this paper, a numerical finite-element-based model was proposed to investigate the thermo-mechanical behavior of buildings' envelopes. A small building has been chosen as a test specimen for simulation, experimental analysis and evaluation. An infrared camera has been employed to visualize temperature distribution on the external walls and to detect thermal bridges positions as well. Numerical and experimental thermal results have been confronted to evaluate the built-up model. An agreement between results has been observed for the thermal transmittance coefficients validating the numerical model. Some minor differences between experiments and simulation results have been noticed, since the problem is considered to be incremental and needs more precision in order to mimic the real dynamic behavior of the building.
Results of mechanical simulations showed that the displacement values are homogeneous among the lateral wall but still negligible. Nonetheless, a stress concentration has been observed at the level of the window and its frame, since the construction materials are different from the rest of the structure.
As perspectives of this work, we are aiming on running the built-up thermal code on much more significant periods of time. In addition, the mechanical model will include the weight of the installations fixed on its roof to evaluate its actual resistance to external forces. Once both models are validated, a new model coupling both mechanical as well as thermal physics will be developed, so that new materials can be investigated for limiting thermal bridges and for energy harvesting.  and processing (i.e., Temperature data). A special thank is also addressed to Vincent Felix for his guidance and advices concerning the usage of Infrared camera.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.