Issue 
Renew. Energy Environ. Sustain.
Volume 7, 2022
Achieving Zero Carbon Emission by 2030



Article Number  1  
Number of page(s)  11  
DOI  https://doi.org/10.1051/rees/2021038  
Published online  06 January 2022 
Research Article
Development and Application of a Thermal Comfort Model in Buildings
Department of Chemical Engineering, School of Engineering, Universidade de Santiago de Compostela, Avda. Lope Gómez de Marzoa s/n. Campus Vida, 15782 Santiago de Compostela, Spain
^{*} email: andres.vilaboa@rai.usc.es
Received:
11
July
2021
Received in final form:
30
September
2021
Accepted:
30
September
2021
Buildings are one of the systems that more energy consumed in the European Union. The study of the thermal envelope is interesting in order to reduce the energy losses. For that, a mathematical model able to predict the system response to external temperature variations is developed. With the mathematical model, different thermal envelope elements of a building based on the lag and the cushioning of the resultant wave can be characterized. In addition, it is important to analyse where the insulation is placed, because when the insulation is outside and the thermal mass is inside, the system produces a response with smooth temperature variations than when the insulation is inside. Therefore, placing the outside insulation generates more steady indoor temperatures, increasing the thermal comfort inside the building. To complete the mathematical model that allows predicting the temperature inside a building taking into account the solar inputs and the thermal inertia of the building. This study will help to establish the optimum design parameters in order to build sustainable and comfortable buildings. Furthermore, it will take one step forward in the construction of nearly ZeroEnergy Buildings.
© A. Vilaboa Díaz and P.M. Bello Bugallo, Published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Energy is a valuable and indispensable resource for social development. Fossil energy resources are finite and the world reserves decrease, so the energy is becoming a most precious asset for the sustainable development of human activity. The continuous emissions of greenhouse gases (GHG) will cause a greater warming and it will increase the likelihood of serious widespread and irreversible impacts on people and the ecosystems [1].
Buildings are one of the major consumers of primary energy in EU. As can be seen in Figure 1, they are responsible of more than 40% of European Union final energy consumption (including household, commercial and public services) and 36% of the total emissions of CO_{2} [2]. The buildings are one of sectors to improve. Among the main efforts is found the renovation of the actual building stock to convert them into buildings with high energy efficiency and decarbonized before 2050 [3].
For increasing energy efficiency in buildings, this study analyses thermal inertia and its implications in the indoor temperature. A good design of the building envelope, taking into account the thermal inertia of its materials allows us to obtain a more stable temperature inside the buildings and therefore comfort is increased.
Ferrari (2007) defined the thermal inertia of the building as the heat storage capacity in its structure and its performance to slow down heat transmission [5]. The heat accumulation capacity of the building envelope has an impact on the final energy consumption. If two building elements with the same thermal transmittance but different energy accumulation capacity are compared, those elements having the greatest heat capacity produce energy savings in heating and cooling [5–7]. On the other hand, if the thermal inertia are combined with other types of energy saving systems (such as solar shadow device), the energy saving in cooling improves notably. Energy savings in cooling in a building with solar shadow device can be 19% lower when using walls with high thermal inertia [8].
On the other hand, the ability of energy accumulation not only depends on its heat capacity, but it also affects its density and its thermal conductivity. By increasing density, a heavier structure is obtained and, as a result, the energy storage capacity increases. This structure delays the heat transfer by increasing the gap of the thermal wave. Furthermore, this increased density produces a greater stability of the indoor temperature. In addition, the optimal density of the material must be selected according to the type of climate varying between 1.800 and 2.500 kg/m^{3} [9]. Finally, the thermal conductivity of the material also influences thermal inertia, since highly conductive materials allow greater heat flow, which facilitates the loss of the energy accumulated in materials. Then, materials with high conductivity cause a minor gap and greater temperature oscillation [10].
Traditional houses, built with materials with high thermal inertia, have an indoor environment more stable than modern houses built with lightweight materials. Specially, in summer and winter, when the outside temperature variation has a value of 25.4 °C and 23.6 °C, respectively, the decrement factor in buildings with heavier materials has a value about 0.20 [11]. This behaviour is also observed in historical buildings. The stone used in this type of building has a high thermal inertia, which allows maintaining the interior temperature at stable values. In addition, ventilation systems at night, allow to eliminate the heat accumulated by the walls [12].
Knowing the influence of the thermal inertia in the indoor temperature evolution helps to improve climate systems. Thus, the control of this property when designing envelope elements improves building energy efficiency. However, in spite of having a low heat capacity, the use of thermal insulation also has special importance. In fact, the placement of thermal insulation on the outside offers better results than placing it on the inside [13]. Also, the insulation helps to decrease the decrement factor and the increase the lag [14].
There are different methods to analyse the thermal behaviour of a wall, allowing the analysis of different scenarios without the need for experimental data collection. Resolution by electrical network analogy is also widely used in scientific literature, which establishes the analogy between the electric and thermal field. Currently, numerical simulation models are already widely used; they allow evaluating the performance of buildings through the use of commercial software or creation of their own code. This type of simulation software allows to evaluate a wide range of physical phenomena [15].
In this paper, the mathematical model is established using the relationship between electrical and thermal analogy. For the development of the equations, the unidirectional heat flow will be assumed, which will allow to avoid the thermal behaviour of the wall in the face of different disturbances. The model will take into account those areas of the buildings in contact with constant temperatures, the influence of ventilation and internal loads. In addition, the model will allow evaluating the influence of solar radiation.
The thermal capacity of materials seems to be a relevant parameter to analyse for moderate climates and intermediate seasons, where it can work as a stabilizing factor of the thermal dynamics of the whole building system [16].
Fig. 1 Final energy consumption in EU by sector in 2017 (Data source: European Environment Agency [4]). 
2 Objectives
The main objective of this paper is to develop a model that allows analyzing the indoor temperature in buildings. The model will be developed based on the analogy between the electric and thermal field. The model will allow to obtain the temperature on the inner surface of a wall. Subsequently, an energy balance will be carried out to obtain the temperature inside a building taking into account different variables. To complete the mathematical model, the solar radiation will be added to achieve a mathematical model more complete.
At the same time, the dynamic response of a singlelayer wall is analysed. The system is subjected to a sinusoidal disturbance, and its response on the indoor surface is also analysed. Thus, the response of enclosures having different thermal properties can be compared.
In addition, this paper tries to show the importance of designing various elements of the thermal envelope properly, in order to achieve more stable indoor conditions; as well as the advantages offered by the correct placement of the insulating material in opaque enclosures with the same thermal transmittance.
3 Methodology
The development of the mathematical model that allows the prediction of the indoor building temperature depending on the outdoor temperature variations is based on the method of Equivalent Electrical Network. This model uses the analogy between an electric circuit and the heat transfer through a multilayer wall. It can be an electric circuit consisting of capacitors and resistors to simulate the heat capacity and envelope thermal resistance, respectively.
In this way, the model that allows predicting the behaviour of a monolayer wall will be developed. This model will allow to evaluate the thermal behaviour of different materials and also to evaluate the optimal placement of the thermal insulation. Subsequently, applying an energy balance, a model will be developed to calculate the interior temperature. Finally, the influence of solar radiation through the walls and windows will be included.
4 Development and presentation of results
4.1 Analysis of a singlelayer wall
To develop the model, the analogy between the electric and thermal field in a monolayer wall will be used. The wall is represented (Fig. 2) by a resistor (R) and a capacitor (C) that represent the thermal resistance and the heat capacity of the wall, respectively. The variation of the external temperature will be represented by the potential difference (V0) and the internal surface temperature by the power difference in the condenser (VC).
The heat flow through the wall is represented by the intensity through the circuit (i = i_{c}).(1)
The voltage on the capacitor (V_{C}), which depends on the time, can be calculated using Kirchhoff's law. So, the intensity through the condenser must be equal to the intensity that circulates through the resistance.(2)
The intensity flowing through the capacitor (i_{C}) is equal to that flows through the resistance (i_{r}). So, combining equations (1) and (2), the following expression is obtained:(3)
By using Laplace transforms, this equation can be solved:(4)
From equation (4), the relationship between the input voltage and the voltage across the capacitor is obtained. Using the electrical analogy, this expression represents the relationship between the inner surface (T_{sur,in}) and external surface (T_{ext}).(5)where C is the heat capacity of the node (J/m^{2} ∙ K); R is the thermal resistance of the node (m^{2} ∙ K/W).
By establishing a monolayer wall as a single RC circuit, inaccurate results can be obtained since the energy accumulation temperature is equal to the temperature of the inner surface. If the wall is divided into a greater number of circuits, it will be possible to obtain a greater precision in the results. In the case of using two nodes to calculate the surface temperature, the equation obtained is the following:(6)
This transfer function has a higher precision than equation (4). The wall is divided in two sections, so the spatial interval is smaller. Considering two nodes, a system of two equations and two unknowns is obtained.(7)
From the matrix obtained, n nodes can be considered to solve a homogeneous wall, obtaining a system of n equations with n unknowns. (8)
with the matrix obtained, the number of calculation nodes is increased and therefore more precise results will be obtained. In this way, the accumulation temperature of the wall will be closer to the inner layer. With this model, it is also possible to analyse multilayer walls where each layer (i) is divided into n/i nodes.
Fig. 2 Representation of a monolayer wall as an electrical circuit. 
4.2 Study of the thermal inertia of materials
Once the model was obtained, the response of different materials to a sinusoidal disturbance with a maximum temperature of 40 °C and a minimum of 20 °C was analysed. The results will be obtained taking as a reference a 20 cm thick wall and 10 nodes will be used for the calculation. This simulation will allow knowing the influence of the properties of the materials (heat capacity, density and conductivity) summarized in Table 1.
The response obtained when a sine wave perturbs the system is shown in Figure 3. As expected, concrete and granite have a similar behaviour because they show similar thermal properties. These materials have a higher internal temperature oscillation and a smaller lag, showing the worst performance, in spite of their high capacity to accumulate energy. This happens because of their high conductivity, which produces a higher loss of the accumulated heat.
On the other hand, the material that shows the best performance is the wood, despite its not being a material with high energy storage. However, it presents less thermal diffusivity, allowing to retain energy during more time, and therefore, more stable temperatures and a greater time lag. Specifically, the temperature oscillation generated would be 1.82 °C (compared to 20 °C external temperature), so it has a factor decrement of 0.0932 (see Tab. 2).
Brick and insulation materials are in the middle. Concerning the insulations, they present high oscillation despite their low thermal diffusivity. Their low capacity of heat accumulation is a disadvantage, so they generate a thermal oscillation of 9.89 °C and have a smaller lag than wood. Finally, brick generates a better response than insulation. It produces waves with an oscillation of 5.88 °C. Bricks have a factor decrement of 0.2952 and a lag of 7.35 h.
The thickness of the wall also has an important influence on the response of the system, since as the thickness increases, the thermal mass of the system increases. Therefore, to analyze the influence of thickness on the response of the system, the material with the worst performance is taken as a reference, which has been granite and concrete. As can be seen in Figure 4, the decrement factor tends to zero as the wall thickness increases and the time lag increases linearly with the thickness. Therefore, an increase in thickness generates very stable indoor temperatures.
Thermal properties of materials used in the simulation (IETCC, 2010).
Fig. 3 Representation the response in the inner surface of different materials with a 20 cm thickness when the system is perturbed from a sine wave. 
Summary of the results obtained in the simulation for different materials with a 20 cm thickness.
Fig. 4 Representation of decrement factor and time lag according to concrete thickness wall. 
4.3 Study on the placement of insulation
The location of thermal insulation in buildings generates a different response in the indoor temperature. To evaluate the influence of the placement of the insulation, the thermal response of two walls composed of 20 cm of concrete (2.5 W/(m ∙ K)) and 8 cm thick thermal insulation (0.032 W/(m ∙ K)) will be compared (details shown in Tab. 3). The difference between the two resides in the placement of the insulation, on the one hand, the evolution of the temperature is analysed when the insulation is in the innermost layer and, therefore, the thermal mass is placed in the outer layer. On the other hand, the opposite case is analysed. That is, the insulation is placed on the outer layer and the thermal mass is placed inside the building.
The climate conditions used correspond with a summer day in the city of Ourense (Spain), specifically on July 8th, 2007. This day the maximum temperature was 40.99 °C at 15:50 hours and the minimum 17.08 °C [17], producing an oscillation of 23.91 °C.
The enclosure has the same thermal resistance (2.62 m^{2} ∙ K/W), and its thermal transmittance has a value of 0.382 W/(m^{2} ∙ K). Nevertheless, the results obtained are different. The decrement factor and the time lag for each case can be obtained through a simulation in Matlab (results shown in Tab. 4).
When the system is disturbed by a sine wave with an oscillation of 20 °C, the placement of the insulation on the inside generates an oscillation of 7.64 °C compared to the oscillation of 0.26 °C that it generates when the insulation is placed on the outside. Introducing the element with greater thermal mass in the outdoor layer (Fig. 5a), energy is stored, but energy is lost more quickly due to its high thermal diffusivity. Therefore there is nothing that prevents the leakage of energy to the outside. Moreover, as the insulation has a low thermal diffusivity, it needs more time to accumulate energy, so in the end the energy flowing to the outside is greater than the energy flowing through the insulation.
If insulation is placed outdoor (Fig. 5b), it keeps the inner surface temperature more stable. This stability is produced because the insulation has a lower thermal diffusivity and prevents energy from flowing outside through the concrete. In this way, the innermost wall accumulates a greater amount of energy, returning the heat inside.
The response of the two walls analysed is shown in Figure 6. Depending on the position of the insulation, great differences can be observed. If the insulation is placed in the outside layer and thus the thermal mass is inside, a constant temperature is recorded in the inner surface. Regarding the time lag, the resultant wave is delayed 0.5 hours more than when the insulation is placed inside.
It is interesting to analyse the time response of each wall, that is, the time that the system takes to stabilize itself when it is perturbed by any disturbance. Both systems are subjected to a unitary step disturbance and the results are obtained through a simulation in Matlab. Figure 7 shows when the insulation is inside and when the insulation is placed in the outside layer. In the first case, the response reaches 63% the perturbation in 0.5 × 10^{5} s (13.9 h) and, in the second case in 1.3 × 10^{6} seconds (15 days). The results obtained between both walls are different, obtaining a greater stability when the thermal mass is placed inside. The thermal mass generates little oscillations in the inner surface temperature and a greater time lag of the thermal wave.
Summary of layer properties of the wall.
Summary of the results obtained through the simulation of the two walls.
Fig. 5 Graphical representation of the cushioning effect of the thermal wave which produces each facade. (a) Indoor insulation and outdoor thermal mass. (b) Outdoor insulation and indoor thermal mass. 
Fig. 6 Variation of inner surface temperature in a wall with outdoor and indoor thermal mass. 
Fig. 7 System response when it is subjected to a unitary step disturbance (Extracted from Matlab). (a) Inside insulation and outside thermal mass. (b) Outside insulation and inside thermal mass. 
4.4 Study of the temperature inside a building
To estimate the indoor temperature, it is necessary to complete the mathematical model. In this case, the electric analogy is discarded and, therefore, the global energy balance of a building will be used.
When the energy balance of a building is made, the know variables are the external temperature and the inner surface temperature, which is calculated previously by the method of Equivalent Electrical Network. In the energy balance we include the energy exchange of the elements of the envelope with outdoors and with zones at a constant temperature. Air renovations and heat contributions of people or electrical equipment are also included. All these heat inputs are accumulated inside a building by the air. Therefore, the energy balance is shown in equation (9).(9)
where h_{in}: is the convection coefficient of air inside (7.69 W/(m^{2} ∙ K)); A_{i} is the area of inner surface of the envelope elements against outdoors (m^{2}); T_{sur,in}: is the inner surface temperature of each envelope element (°C); T_{in} is the temperature inside a building (°C); A_{0,i} is the area of inner surface of each envelope element against a place with constant temperature (m^{2}); T_{0,i} is the temperature of each zones at constant temperature (°C); n is the number of air renovations (h^{−1}); V is the volume of the building analysed (m^{3}); T_{ext} is the external temperature (°C); Q_{in} is the internal loads produced by electronic equipment and by the people living in the building (W/m^{2}); C_{p}_{,air}: is the heat capacity of the air (J/K).
The transfer function to estimate the temperature inside a building is obtained based on the system disturbances. For this, it is necessary to solve equation (9) using the Laplace transformation. This gives as a result equation (10). (10)
By equation (10), the indoor temperature can be calculated in Laplace domain based on their input variables. The equation is solved by Matlab, which allows us to obtain the temperature evolution in the time domain.
With the purpose of analysing a real case, a standard building is taken as a reference. To simplify the operations, the building will be considered as a block. The dwelling has a surface of 80 m^{2}. Two of its facades and its roof are in contact with outdoors, while its other two walls and the concrete deck are in contact with a zone having a constant temperature. These last elements are composed of concrete with a 20 cm thickness.
The elements in contact with outdoors, including the flat roof, are composed of 20 cm of concrete and 8 cm of insulation. This type of elements are analysed from two points of view. Firstly, the thermal mass is placed outside and the insulation in the innermost layer. Secondly, the insulation is placed outside and, therefore, the thermal mass is inside. The thermal properties of the materials are shown in Table 3 and the responses of each facade are in Table 4.
As a reference, the climatic conditions that have been used to the simulation correspond to one of the coldest days recorded in the city of Ourense (the 5th of February 2012). The maximum temperature reached was 8.89 °C and the minimum temperature was −5.85 °C [17].
Figure 8 shows the temperature evolution inside a building when the insulation is inside and outside, disregarding the internal loads and the solar radiation. When insulation is inside, the system response is a sine wave with less amplitude and with a different phase. Specifically, a maximum temperature of 6.0 °C and a minimum of 3.8 °C, so the average temperature is 4.90 °C, with a time lag of 6.5 hours. The average temperature is greater than the outdoor, as the building is in contact with zones with a constant temperature, which is higher than the outdoor average, increasing the inner average temperature.
When the thermal mass is in the innermost layer, the temperature variation of the building is practically constant, reaching a maximum of 5.25 °C and a minimum of 4.55 °C. As for the time lag, the resultant wave is in the same phase than the disturbance. This is due to the influence of the air renovations. As in the inner surface, the temperature is nearly constant. This parameter is influenced by the energy flow of the air renovations, which is a wave with the same phase than external temperature. This produces that the time lag of resultant wave is lost.
If both results are compared, we can observe that the temperature is more stable when the thermal mass is inside and the insulation is placed outdoors. The temperature variation obtained is 0.70 °C versus 2.2 °C, obtained when the thermal mass is outdoors.
Fig. 8 Variation of indoor temperature depending on the placement of the insulation. (a) Inside insulation and outside thermal mass. (b) Outside insulation and inside thermal mass. 
4.5 Study of the influence of solar radiation through walls
For studying the solar radiation through a wall, solair temperature will be included in the model. The solair temperature is the outdoor temperature, which in absence of all radiation exchanges, generates the same heat flow through a wall that is equal to the combination of the solar incident radiation, the exchange of radiant energy with the sky and the exchange of heat convection with outdoor air [18].(11)
Therefore, solair temperature allows knowing the amount of heat conduction through the wall from the outside including solar radiation.
Vertical surfaces receive longwave radiation from the ground, the sky and from other vertical surfaces, so determine the value of ΔR with precision is complex. When solar radiation is high, the temperature of objects are higher than air temperature, thus, the long wave radiation compensates the low emittance of the sky. Therefore, we can assume in vertical surfaces that ε ∙ ΔR = 0. On the other hand, for horizontal surfaces, the long wave radiation is only received from the sky, so the term ΔR has a value of 20 BTU/h ∙ ft^{2} (63.09 W/m^{2}) [18].
Solair temperature will be included in the model. The influence of solar radiation on the outside surface of the wall makes that heat flow that occurs will be equivalent to a maximum temperature of 16.14 °C. The solair temperature must be calculated for every enclosure because the value depends on its tilt and orientation. In the Figure 9, the solair temperature for a horizontal surface (β = 0, γ = 0) is showed:
Figure 10 shows the influence of solar radiation through the walls. In this case, the indoor temperature increases in 0.5 °C, therefore, the average temperature in a building is increased and it reaches 5.5 °C. However, the oscillation is higher when the insulation is inside and the total oscillation reaches the 2.6 °C. On the other hand, when the insulation is outside a constant response is produced. So when solar radiation is introduced is not appreciated an increase of temperature. The maximum temperature reached is 5.66 °C and the minimum is 4.95 °C so the average temperature is 5.30 °C. The average temperature is lower than with the inside insulation because this configuration is less influenced by solar radiation. However, temperatures more stables are achieved with the outside insulation.
Fig. 9 Graphical representation of the outdoor temperature and the solair temperature on a horizontal surface. 
Fig. 10 Variation of the indoor temperature for walls taking into account the solar radiation. On the left with inside insulation. On the right with outside insulation. 
4.6 Study of the influence of solar radiation through windows
The main solar contributions in a building are produced through the windows. Therefore, it is interesting to analyse the influence of this energy in the indoor temperature. These contributions of energy can increase the indoor temperature in winter, however, they can produce an excess of temperature in summer so it is necessary look for a balance between this two situations.
The first step to study the influence of solar radiation is to add the solar radiation through windows to equation (10) and solve it using Laplace transformation. (12)
This new term (Q_{solar}) depends on the orientation, the tilt and the properties of windows. So, solar contributions are equal to the direct solar radiation multiplied by the solar factor of the glass. The solar factor of the window is the relationship between total energy which enters the room through the glass and the incident solar energy. This total energy is the sum of solar energy that enters through direct transmission and energy transferred from the glazing to the indoor space.
To analyse the influence of solar radiation, we will use a glazed surface of 12 m^{2} (20% of the surface of the facade). The window has double glazing with a transmittance of 3.1 W/m^{2} K and a solar factor of 0.75.
Figure 11 shows the evolution during the winter of the indoor temperature taken into account the solar radiation. With the inside insulation, an internal temperature of 8.04 °C is reached, while when placing the insulation on the outside, the temperature reached is 8.37 °C. This difference is produced by the time lag of each wave. When insulation is on the outside and therefore the thermal mass is inside, the indoor temperature is more stable, and the solar radiation increases the indoor temperature during the hours with more sunshine. However, when the insulation is in the outside, the wave lag produces a maximum.
Fig. 11 Variation of the indoor temperature for walls with inside and outside insulation and taking into account the solar radiation through the window during a winter day. 
5 Discussion
The model described in this paper using the electrical analogy, allows predicting the thermal behaviour of the enclosures in a simple way, without the need for complex simulations. The selection of the suitable materials, allows optimizing the thermal behaviour of the buildings and therefore they must be selected properly. The materials that generate a softer internal temperature and greater time lag are those with higher density and low thermal conductivity, due to the high capacity to accumulate heat.
Among the materials used, wood is the one with the best thermal behaviour, generating a more stable temperature inside and producing a time lag of 11.47 h. Wood is a material with great potential, since its use as a substitute for conventional construction materials can reduce CO_{2} emissions by 69% in the construction phase [19]. In the second place is the brick with a decrement factor of 0.2952 and a time lag of 7.35 h.
Granite and concrete are materials with similar properties and therefore a similar response is obtained. These types of materials, despite not having the best thermal response, are widely used materials in old buildings. In this type of construction, the materials generate a lower temperature oscillation and a greater time lag.
The placement of thermal insulation is an important point to take into account when designing a building. In this case, a comparison was made with respect to placing the insulation inside or outside. The results show that placing the insulation on the outside generates more stable temperatures, with a temperature oscillation of 0.26 °C compared to 7.64 °C when placing the insulation inside. In addition, by analysing a variation in the outside temperature over a longer period of time, placing the insulation outside allows to obtain a system that is more stable to changes. The use of insulation inside the walls eliminates the thermal capacity of heat accumulation, generating a greater fluctuation in temperature.
Introducing solar radiation into the model allows an accurate estimation of the indoor temperature of buildings. Solar radiation through the walls of the building causes different disturbances in the indoor temperature depending on the location of the insulation. By introducing solar radiation into the model, there is an increase in temperature inside due to the heat that enters the building by conduction. During winter, a maximum temperature of 7 °C is produced when the insulation is inside, while when the insulation is placed outside, the maximum temperature reached is 5.66 °C.
The inputs of solar radiation through the windows, causes a greater change in the indoor temperature. In the case analysed, a lowquality window is used, in order to observe the influence of solar radiation on the indoor temperature. In this way, it is observed that during the winter, the entry of solar radiation causes a greater oscillation when the thermal mass is placed inside (T_{max} = 8.37 °C) compared to when the thermal mass is outside (T_{max} = 8.04 °C).
In summer solar inputs can generate uncomfortable indoor temperature, so the use of solar shadow devices would control the overheating generated by solar radiation.
6 Conclusions
In this paper, a mathematical model has been developed, which allows for obtaining the temperature variation inside a building, using the analogy between electric field and thermal field and, an energy balance. The results obtained show the following:

The materials with better behaviour, when are perturbed by a sine disturbance, have a high density and heat capacity and a low thermal conductivity.

The thickness of the wall is important to obtain more stable indoor conditions. When the thickness increases, the thermal mass also increases. Thus, when the thickness is high, the factor decrement can be nearly zero.

Introducing thermal mass inside and insulation outside generates more stable temperatures than placing insulation inside, in spite of the thermal transmittance having the same value in both cases.

Placing the thermal mass inside produces a lower heat flow than when it is placed outside. This option allows designing climate systems with less power.

When the building's behaviour is simulated, the indoor temperature is more stable when the thermal mass is inside. However, the influence of the loss of energy, due to air renovations, causes that indoor temperature was practically in the same phase than external temperature.

Solar radiation through the enclosures does not generate large fluctuations in temperature. However, placing the insulation on the inside generates a higher indoor temperature during winter.

Solar contributions through the window generate a greater impact on the indoor temperature with respect to the heat flow through the enclosure.
The design of thermal envelope is basic to achieve nearly ZeroEnergy Buildings. However, achieving enclosures with a low thermal transmittance is not enough. The order of the wall components is also important, systems with the thermal mass inside and the insulation outside showing more advantages.
Being able to predict the temperature inside a building is useful in the design of air conditioning control systems. An approximation of the outside temperature evolution can be obtained from the weather forecast. Thus, it can be designed a control system which anticipates to outside temperature variations, in order to keep the indoor temperature in comfort zone.
References
 IPCC, Climate Change 2014: Synthesis Report, 2014 [Google Scholar]
 EC, COM(2008) 772, Communication from the Commission of 13 November 2008. Energy efficiency: delivering the 20% target. Brussels, 2008 [Google Scholar]
 EU Directive 2018/844 of the European Parliament and of the Council of 30 May 2018 amending Directive 2010/31/EU on the energy performance of buildings and Directive 2012/27/EU on energy efficiency [Google Scholar]
 European Environment Agency: https://www.eea.europa.eu/ [Accessed March 2021] [Google Scholar]
 S. Ferrari, in 28th AIVC and 2nd Palenc Conference: Building Low Energy Cooling and Ventilation Technologies in the 21st Century, 2007 [Google Scholar]
 T. Ihara, A. Gustavsen, B. Petter Jelle, Appl. Energy 158 (2015) [Google Scholar]
 E. Kosseckaa, J. Kosnyb, Energy Build. 34, 4 (2002) [Google Scholar]
 N. Aste, A. Angelotti, M. Buzzetti, Michela. Energy Build. 41, 11 (2009) [CrossRef] [Google Scholar]
 C. AvendañoVera, A. MartinezSoto, V. Marincioni, Renew. Sustain. Energy Rev. 131 (2020) [Google Scholar]
 K.J. Kontoleon, Th.G. Theodosiou, K.G. Tsikaloudaki, Appl. Energy 112 (2013) [Google Scholar]
 S. Martín, F.R. Mazarrón, I. Cañas, Construct. Build. Mater. 24, 5 (2010) [Google Scholar]
 G. Evola, L. Marletta, S. Natarajan, E.M. Patanè, Energy Proc. 133 (2017) [Google Scholar]
 S.A. AlSanea, M.F. Zedan, S.N. AlHussain, Appl. Energy 89, 1 (2012) [CrossRef] [Google Scholar]
 E. Stéphan, R. Cantin, A. Caucheteux, S. TascaGuernouti, P. Michel, Build. Environ. 80 (2014) [Google Scholar]
 S. Verbeke, A. Audenaert, Renew. Sustain. Energy Rev. 82, 3 (2018) [Google Scholar]
 N. Aste, F. Leonforte, M. Manfren, M. Mazzon, Appl Energy 145 (2015) [Google Scholar]
 Meteogalicia: https://www.meteogalicia.gal/web/inicio.action [Google Scholar]
 American Society of Heating, Refrigerating and AirConditioning Engineers, Inc. HVAC Applications, edited by ASHRAE, SI Edition (Atlanta 2011) [Google Scholar]
 A. Himes, G. Busby, Dev. Built Environ. 4 (2020) [Google Scholar]
Cite this article as: Andrés Vilaboa Díaz, Pastora M. Bello Bugallo, Development and Application of a Thermal Comfort Model in Buildings, Renew. Energy Environ. Sustain. 7, 1 (2022)
All Tables
Summary of the results obtained in the simulation for different materials with a 20 cm thickness.
All Figures
Fig. 1 Final energy consumption in EU by sector in 2017 (Data source: European Environment Agency [4]). 

In the text 
Fig. 2 Representation of a monolayer wall as an electrical circuit. 

In the text 
Fig. 3 Representation the response in the inner surface of different materials with a 20 cm thickness when the system is perturbed from a sine wave. 

In the text 
Fig. 4 Representation of decrement factor and time lag according to concrete thickness wall. 

In the text 
Fig. 5 Graphical representation of the cushioning effect of the thermal wave which produces each facade. (a) Indoor insulation and outdoor thermal mass. (b) Outdoor insulation and indoor thermal mass. 

In the text 
Fig. 6 Variation of inner surface temperature in a wall with outdoor and indoor thermal mass. 

In the text 
Fig. 7 System response when it is subjected to a unitary step disturbance (Extracted from Matlab). (a) Inside insulation and outside thermal mass. (b) Outside insulation and inside thermal mass. 

In the text 
Fig. 8 Variation of indoor temperature depending on the placement of the insulation. (a) Inside insulation and outside thermal mass. (b) Outside insulation and inside thermal mass. 

In the text 
Fig. 9 Graphical representation of the outdoor temperature and the solair temperature on a horizontal surface. 

In the text 
Fig. 10 Variation of the indoor temperature for walls taking into account the solar radiation. On the left with inside insulation. On the right with outside insulation. 

In the text 
Fig. 11 Variation of the indoor temperature for walls with inside and outside insulation and taking into account the solar radiation through the window during a winter day. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.