Passive Strategies to Improve the Comfort Conditions in a Geodesic Dome

Non-conventional thermal zones are low-cost and ecology friendly alternatives to the housing needs of populations in various situations, such as surviving natural disasters or addressing homelessness. However, it is necessary to guarantee thermal comfort for occupants, while aiming to minimize energy consumption and wastage in refrigeration systems. To reduce the cooling requirements in non-conventional thermal zones it is necessary to model the structure and analyze the principal factors contributing to internal temperature. In this paper, a geodesic dome is modellingusing the lumped parameter technique. This structure is composed of a wooden skeleton and wooden floor, with a canvas surface as its exterior. The mathematical model was tuned using experimental data, and its parameters were classified using Monte Carlo sensitivity analysis. The mathematical model was used to evaluate the impact on internal temperature and occupants’ comfort when two strategies are considered. The results obtained indicatee internal temperature reductions down to a range of 7% to 11%; this result is reflected directly in the energy used to refrigerate the thermal zone, contributing to the objective of providing houses with lower energy consumption.


Introduction
The 21st century presents a series of challenges for the current society, such as global warming, depletion of natural resources, and satisfying the requirements of the world population without destroying ecosystems. As a response to these problems, the sustainable development goals were adopted by the member countries of the United Nations (UN) in 2015 [1]. The objectives include inclusive, safe, resilient, and sustainable cities. This is a very difficult goal to fulfill with the increasing growth of the world population, as the more people live in cities, the more resources and energy are needed.
The formulation of the sustainable development goals has affected different research fields, especially those related to reducing the energy consumption in residential buildings [2][3][4]. Energy consumption depends on various factors, such as the nature of the electronic equipment, the number of occupants, the level of activity and the degree of comfort. The last aspects are directly related to the use of heat ventilation air-conditioning (HVAC) systems. The HVAC can represent more than 40% of the energy consumption in buildings [5][6][7], which means that any reduction in the use of air conditioning systems can seriously impact the world energy expenditure.
The scientific literature relates different strategies to minimize the use of air-conditioned systems, which is especially high in tropical countries and during the spring season in countries in temperate zones [8][9][10]. Regularly, these strategies have focused on minimizing the heat flux transfer between buildings and the environment. These strategies are classified as passive strategies and are characterized by the introduction of modifications in the exterior surfaces of the buildings [11][12][13].
The use of green roofs is a passive strategy which is one of the most ancient and frequent technique for minimizing the incoming heat from the exterior, probably because it includes additional environmental benefits such as the oxygen production and runoff water reduction [14,15]. The strategy is based on using vegetation layers over rooftops and walls, taking advantage of the natural ability of plants to process solar radiation, allowing them to function as a thermal barrier [16]. This strategy can be combined with other strategies such as a photovoltaic rooftop that enhances the capabilities of solar panels [17,18]. Another important strategy is the use of a coating of thermal paint, which is known as a cold roof. The use of coating paints in rooftops and walls reduces the incoming heat and directly affects the amount of energy required for cooling [19,20]. Experimental and theoretical research carried out in buildings demonstrate that a thermal coating reduces the internal air temperature by 2.3 • C [21]. In laboratory tests reported in [22], a saving of up to 2.9% was obtained. However, even though these results are promising, new types of insulating materials are constantly being developed [23,24]. This is the case with the new thermal coating being tested in reduced scale models at the Universitat Politècnica de València, Spain. The study yielded reductions of interior temperatures by 7.4% [25].
All these passive strategies can be used to reduce the interior temperature, but one drawback is that in many cases it is not possible to carry out experimental evaluations of multiple thermal coatings, such as in [26], where different products were tested in building facades. An important alternative is to run thermal models to test the impact of passive strategies. Any good mathematical model should be adjusted to real conditions to allow accurate results [3].
The mathematical modeling of thermal zones is important prior to experimental tests of HVAC. For example, in [21], the model was tuned for the specific conditions in a given place and reused for two additional locations. This kind of tests allows for the selection of the best solution for different places with similar conditions, or different kinds of thermal zones and buildings. In [27] a thermal network was used to study the impact of different loads on energy consumption in a building. The model was adjusted to predict the electrical behavior of a campus in Naval Station Great Lakes, Illinois. The investigation presented in [28] studied the thermal response of different structures, such as roofs and decks, which can provide an important heat gain for the solar radiation absorbed. In general, there is plenty of published work on the thermal behavior of buildings, but most are dedicated to traditional structures such as rectangular buildings, which do not consider non-conventional structures and their possible use for low energy consumption [29][30][31].
Different strategies for building modeling have been reported in the literature [32][33][34][35]; each one has advantages and disadvantages. In many cases, the selected strategy is determined by the structure to be studied and by the research capacities [36][37][38]. To reduce the complexity during the construction of the mathematical model, many engineers use simple and adaptable techniques, such as the lumped parameters strategy [39][40][41]. This technique uses an analogy with RC circuits (resistor and capacitor) and, using Kirchhoff's laws, a thermal model for the studied building can be established. This technique allows users to configure the circuit to study different structures with multiple wall layers [42][43][44].
The present work aims to evaluate an autochthonous habitat's adequate behavior in response to the daily thermal oscillation, which was monitored in the mountainous region near Santa Marta. The autochthonous habitat is characterized by a low form factor in its construction (surface/volume ratio) and materiality with low transmission values, due to the use of highly insulating plant fibers. Given the difficulty of reproducing an autochthonous habitat at present, we chose to work on a model that essentially reproduced its thermodynamic behavior. Our model involved two essential elements: a geodesic dome as a form due to the closeness of its form factor to those of vernacular architectures, and a heat-insulating coating that simulated the behavior of the vegetable fibers that make up these enclosures. In this respect, the analysis of a non-conventional thermal zone, namely, a geodesic dome, is presented. The study focusses on modeling the structure and on evaluating the impact on the interior temperature using passive strategies. The goal was to show that by introducing modifications to the envelopem the internal temperature can be reduced, thus minimizing the energy requirements for refrigeration for the habitability of the thermal zone. The article is organized as follows: the second section is dedicated to presenting the mathematical model and its physical and geometric parameters; in the third section, the model is tuned with experimental data and used to analyze the principal factors to guarantee the occupant's comfort; in the fourth section, the impact of passive strategies and adaptations to the dome on the internal temperature are evaluated; and finally, in the fifth section, some conclusions are drawn.

Modeling a Geodesic Dome
The structure chosen in this work was a geodesic dome built in the Universidad del Magdalena (Colombia). The geodesic dome was modeled from a truncated icosahedron, on whose faces pyramids were built, with their bases on the faces of the polyhedron, and with its vertex on the circumscribed sphere.
The thermal zone had a wooden floor as a base, a wooden structure, and a canvas for the external surface. Figure 1 shows the internal structure of the dome and some of the geometrical data. The left and central pictures show different views of the internal structure, and the right picture shows its building process. Each triangle was formed by wooden segments with dimensions 1.2 m × 0.1 m × 0.025 m. The geodesic dome was modeled and studied using the lumped parameter strategy. Figure 2 represents a cross-section of the dome and shows the internal structure of the thermal zone. This picture also introduces the variables considered in the model, which are the canvas temperature (T l ), the wooden structure temperature (T m ), the floor temperature (T f ), and the internal air temperature (T). Additionally, in Figure 2 other variables are indicated, that is, the environmental temperature (T a ), the ground temperature (T s ), the superficial external and internal canvas temperature (T sle and T sli ), the superficial wooden structure temperature (T sm ), and the superficial floor temperature (T s f ). Finally, the internal loads are represented by I L , and the heat flux transferred by the canvas and the wooden structure (Q l and Q m ) correspond to the incoming heat flux (Q i ). This thermal zone was studied in [45], where a previous version of the mathematical model was presented, which is updated here to include two new aspects. The first modification is the ability to consider the heat transfer by radiation and convection individually. The second aspect is related to using the floor temperature as a state variable. These modifications provide a more accurate mathematical model, according to the circuit diagram used in the literature [39]. Figure 3 shows the circuit diagram used to represent the geodesic dome. This circuit shows the different zones, elements, and heat process in the dome. The resistors and capacitors are calculated according to the physical parameters of the dome. The parameters used are the surface area (A), thickness (L), thermal conductivity (kt), density (ρ) and specific heat (Ce). In the rest of the article the subscripts l, m, f and a are used to indicate the element to which each parameter belongs. The subscript l stands for canvas, m for wooden structure, f for floor, and finally, the subscript a stands for air. Table 1 shows the values for each physical and geometric parameter. This thermal zone was studied in [45], where a previous version of the mathematical model was presented, which is updated here to include two new aspects. The first modification is the ability to consider the heat transfer by radiation and convection individually. The second aspect is related to using the floor temperature as a state variable. These modifications provide a more accurate mathematical model, according to the circuit diagram used in the literature [39]. Figure 3 shows the circuit diagram used to represent the geodesic dome. This circuit shows the different zones, elements, and heat process in the dome. The resistors and capacitors are calculated according to the physical parameters of the dome. The parameters used are the surface area ( ), thickness ( ), thermal conductivity ( ), density ( ) and specific heat ( ). In the rest of the article the subscripts l, m, f and a are used to indicate the element to which each parameter belongs. The subscript l stands for canvas, m for wooden structure, f for floor, and finally, the subscript a stands for air. Table 1 shows the values for each physical and geometric parameter.   represent the experimental In the case of the internal air, it is not possible to use superficial area and thickness; instead we used the internal volume V i = 146 m 3 . According to Figure 3, the exterior of the dome includes the environmental temperature and R 1 resistors, which are used to represent the opposition of the canvas' external surface to transfering heat with convection and radiation process. The two resistors represent the heat transfer processes between the air and the wood or canvas surface. When a single resistor is used, it represents the conduction process.
R 1 is calculated in Equation (1) using h and ε as convection and radiation coefficients for the external surface. These coefficients must be adjusted to environmental conditions to obtain an accurate simulation. Similarly, Equations (3), (5) and (6) represent the convection and radiation resistances for the surfaces on the interior dome, specifically, the canvas' internal surface, the wooden structure, and floor. It is necessary to remark that the coefficient division is one of the improvements in this model which allows for a better characterization of the heat flux process. Previously, the transference between surface and air was handled with one single coefficient. Now, by using two independents coefficients it is possible to evaluate the impact of each and implement solutions directly over one heat transfer process.
Equations (2), (4) and (7) represent the resistance to the conduction heat process in the canvas, wood structure, and floor, respectively. The parameters used in these expressions, such as the thermal conductivity and superficial area are measured or taken from the literature [46]. (1) The constants C l , C m , C f , and C r represent the capacity of each element to store thermal energy and were taken from [45].
The mathematical model is presented in Equations (8)-(11), using canvas, wooden structure, floor, and internal air temperatures as state variables. Sources T a and T s are used to represent the environmental and ground temperatures. The current source I L represent all the internal gains in the system, such as occupants, electronic equipment, refrigeration systems, and the solar gain by skylight. Finally, T c is the canvas limit temperature, which is calculated in Equation (12).
It should be pointed out that the internal loads I L can change according to the quantity and kind of elements inside the dome; in other words, the internal load can take the form I L = N·Q 0 + M·Q e − u·Q c + Q v , using N and M to represent the number of human occupants and electronic equipment working inside the dome, respectively. Q o , Q e and Q v are the heat gain by occupants, electronic equipment, and the skylight, respectively; and finally, Q c is the heat extracted by the cooling system and u is a binary variable used to define the equipment state, u = 1 for on and u = 0 for off. In this work, the dome was considered completely empty; thus, I L = Q v .

Tuning and Analysis of Mathematical Model
The model was adjusted using a database recorded for 5 days, from February 9 to 14, 2019, in Santa Marta city (Colombia) [30]. In this database, the internal and external surface temperatures, environmental temperature, and internal air temperature were registered. Other environmental conditions, such as ground temperature (used in Equation (10)), solar radiation, and wind speed were taken from departmental institutions data [31].
The algorithm used for the tuning was PatternSearch. This algorithm is part of the Optimtool toolbox of Matlab, taking as the optimization function the expression using T exp to represent the experimental temperature recorded at the point i, and T for the simulated temperature at this point [30]. The procedure for adjustion consists in separating the studied period in the day and night and calculating the best value to fit the simulation curves to the experimental data. This method has been tested in similar investigations [25,47]. The tuned coefficients are presented in Table 2. Based on the tuning strategy described, it was possible to obtain a very low percentage difference between the simulation and the experimental data for the day studied. However, using the same parameters for other days can increase the simulation error. In this case, we obtained an error for the studied day of 2.5%, and an error for the rest of days of 3.2%.  Figure 4 depicts the internal simulated temperature with a blue line, the experimental internal temperature with a red line, and the environmental temperature with a green line. The variations in environmental temperature explain the increase in the percentage difference between the experiment and the simulation, because the temperature in later days was greater by 10% than the first-day temperature, which was used to adjust the model. was greater by 10% than the first-day temperature, which was used to adjust the model.  Figure 5 shows the behavior of canvas temperature ( ), wooden structure temperature ( ), and floor temperature ( ). The blue line is used for canvas temperature, the red line for wooden structure, and the green line for the floor temperature. There are noticeable differences between , and . The first two are directly related to environmental conditions and the fast variations in Ta, whereas is related to the ground temperature, , which was not measured. To overcome this situation, we used 28.3 °, which corresponds with the mean annual temperature in the city [48]. Using the tuned mathematical model, we can analyze different aspects such as the incoming heat flux from the exterior, and the impact of each parameter over the internal temperature. The incoming flux is represented by and reaches the interior air by two different routes. The first route is the energy transfer to the canvas, followed by transfer to the wooden structure, and finally to the interior air. The heat flux that occurs in this way is designated by . The second route uses the gaps between the branches in the  Figure 5 shows the behavior of canvas temperature (T l ), wooden structure temperature (T m ), and floor temperature (T f ). The blue line is used for canvas temperature, the red line for wooden structure, and the green line for the floor temperature. There are noticeable differences between T l , T m and T f . The first two are directly related to environmental conditions and the fast variations in T a , whereas T f is related to the ground temperature, T s , which was not measured. To overcome this situation, we used T s = 28.3 • C, which corresponds with the mean annual temperature in the city [48].
Mathematics 2021, 9, x FOR PEER REVIEW 7 of 14 was greater by 10% than the first-day temperature, which was used to adjust the model.  Figure 5 shows the behavior of canvas temperature ( ), wooden structure temperature ( ), and floor temperature ( ). The blue line is used for canvas temperature, the red line for wooden structure, and the green line for the floor temperature. There are noticeable differences between , and . The first two are directly related to environmental conditions and the fast variations in Ta, whereas is related to the ground temperature, , which was not measured. To overcome this situation, we used 28.3 °, which corresponds with the mean annual temperature in the city [48]. Using the tuned mathematical model, we can analyze different aspects such as the incoming heat flux from the exterior, and the impact of each parameter over the internal temperature. The incoming flux is represented by and reaches the interior air by two different routes. The first route is the energy transfer to the canvas, followed by transfer to the wooden structure, and finally to the interior air. The heat flux that occurs in this way is designated by . The second route uses the gaps between the branches in the wooden structure. The transfer process is from the external air to the canvas, and then Using the tuned mathematical model, we can analyze different aspects such as the incoming heat flux from the exterior, and the impact of each parameter over the internal temperature. The incoming flux is represented by Q l and reaches the interior air by two different routes. The first route is the energy transfer to the canvas, followed by transfer to the wooden structure, and finally to the interior air. The heat flux that occurs in this way is designated by Q m . The second route uses the gaps between the branches in the wooden structure. The transfer process is from the external air to the canvas, and then from the canvas to the internal air. The heat flux that occurs in this way is designated by Q l . The relationship between the three heat fluxes is Q i = Q l + Q m . Figure 6 shows the heat coming in by both routes. The blue line is used to represent Q l , and a red line to represent Q m . It can be observed in this figure that most of the incoming heat flux was transferred directly through the canvas, which can be explained by the difference in surface and the greater thickness in the wooden structure. atics 2021, 9, x FOR PEER REVIEW Figure 6 shows the heat coming in by both routes. The blue line is used to , and a red line to represent . It can be observed in this figure that most of th ing heat flux was transferred directly through the canvas, which can be explain difference in surface and the greater thickness in the wooden structure. Another important aspect to mention in relation to Figure 6 is that to min incoming flux to the interior air, the objective should be to reduce the heat trans the canvas. To minimize the heat flux coming in through the canvas, the most i parameters in the heat transfer process should be identified. The classification rameters was performed using Monte Carlo (MC) sensitivity analysis. The MC m based on the generation of random numbers around a central value, which are parameters in this case. The random numbers produce a set of variations in th variable. The dispersions in the output values are used to evaluate the sensitiv model to changes in each parameter.
The variable chosen for the sensitivity analysis was the mean internal tem The studied dispersion was in the range of 0.5%-5% with an increment of 0. 3 shows the dispersion in the mean internal temperature, introducing a dispe tween 0.5% and 5% in each coefficient. According to the results obtained, the m ential parameters were the convection coefficient in the canvas and in the floo internal temperature mean value. Equation (13) shows lineal regressions built to compare the importance of each coefficient. Using to represent the dispersio parameter, to signify the dispersion on the mean internal temperature and resent the linear regression coefficient, the following linear regressions are obtai . Another important aspect to mention in relation to Figure 6 is that to minimize the incoming flux to the interior air, the objective should be to reduce the heat transferred by the canvas. To minimize the heat flux coming in through the canvas, the most important parameters in the heat transfer process should be identified. The classification of the parameters was performed using Monte Carlo (MC) sensitivity analysis. The MC method is based on the generation of random numbers around a central value, which are the tuned parameters in this case. The random numbers produce a set of variations in the output variable. The dispersions in the output values are used to evaluate the sensitivity of the model to changes in each parameter.
The variable chosen for the sensitivity analysis was the mean internal temperature. The studied dispersion (d) was in the range of 0.5%-5% with an increment of 0.5%. Table 3 shows the dispersion in the mean internal temperature, introducing a dispersion between 0.5% and 5% in each coefficient. According to the results obtained, the most influential parameters were the convection coefficient in the canvas and in the floor for the internal temperature mean value. Equation (13) shows lineal regressions built to relatively compare the importance of each coefficient. Using x to represent the dispersion in each parameter, y to signify the dispersion on the mean internal temperature and R 2 to represent the linear regression coefficient, the following linear regressions are obtained:    h s → y = 1.2623x + 0.0607; R 2 = 0.86 h le → y = 0.9989x − 0.036; R 2 = 0.98 h li → y = 0.1543x − 0.0393; R 2 = 0.97 (13) Figure 7 shows the comparison between the linear regressions for the convection coefficients. The horizontal and vertical axes are defined by the dispersion introduced in each parameter (d) and by the dispersion obtained in the internal temperature (σ), respectively. The blue line represents the floor convection coefficient, the red line represents the external convection coefficient, and the green line depicts the coefficient convection for the internal face of the canvas.  Figure 7 shows the comparison between the linear regressions for the convection coefficients. The horizontal and vertical axes are defined by the dispersion introduced in each parameter d and by the dispersion obtained in the internal temperature , respectively. The blue line represents the floor convection coefficient, the red line represents the external convection coefficient, and the green line depicts the coefficient convection for the internal face of the canvas. According to the MC analysis, it is necessary to reduce the convection coefficient of the floor surface in order to reduce the cooling requirements. However, this task can be difficult because this surface is constantly exposed to the footsteps of the occupants and other forms of wear. For this reason, a better strategy would be to minimize the convection coefficients in walls, especially the external surface.

Passive Strategy Simulation
This section is dedicated to exploring the efficiency of two different passive strategies in reducing the cooling requirements. In each case, the current comfort conditions are compared with those achieved by carrying out the modifications. The limits on the internal temperature are established by using the 55 ASHRAE standard. The typical thermal isolation by clothes in this region is around 0.57 , and the maximum temperature allowed for this standard is 25.72°C. Using the Euclidean norm, we calculated the difference between the experimentally measured internal temperature and this upper According to the MC analysis, it is necessary to reduce the convection coefficient of the floor surface in order to reduce the cooling requirements. However, this task can be difficult because this surface is constantly exposed to the footsteps of the occupants and other forms of wear. For this reason, a better strategy would be to minimize the convection coefficients in walls, especially the external surface.

Passive Strategy Simulation
This section is dedicated to exploring the efficiency of two different passive strategies in reducing the cooling requirements. In each case, the current comfort conditions are compared with those achieved by carrying out the modifications. The limits on the internal temperature are established by using the 55 ASHRAE standard. The typical thermal isolation by clothes in this region is around 0.57clo, and the maximum temperature allowed for this standard is T max = 25.72 • C. Using the Euclidean norm, we calculated the difference between the experimentally measured internal temperature and this upper limit of the comfort zone. We found that the internal temperature was 25.98% above the maximum acceptable value, making it necessary to take cooling actions.

First Strategy: Coating Solutions
The first effort to minimize the cooling requirements was based on apply coating solutions over exterior surfaces. In [25] an experiment was developed to determine the thermal properties of a commercial coating solution. One of the results of this study was the convection coefficient of a surface with a layer of coating solution with a thickness of 0.5 mm. For natural ventilation and high environmental temperature, the convection in the exterior surfaces was tuned according to h le = 5.3 W m 2 K . Using this convection coefficient, the model was adjusted to evaluate the thermal reduction obtained in the case of applying this coating solution over the exterior surface of the canvas. Figure 8 represents the experimental temperature using the blue line, whereas the temperature simulated with the external convection coefficient of the coating solution is represented with red line. The simulation shows that by applying the coating solution, a general thermal reduction of 5.4% can be achieved during the simulation period. However, in individual days, a reduction of 9.4% was obtained.
Mathematics 2021, 9, x FOR PEER REVIEW 10 of 14 thermal properties of a commercial coating solution. One of the results of this study was the convection coefficient of a surface with a layer of coating solution with a thickness of 0.5 mm. For natural ventilation and high environmental temperature, the convection in the exterior surfaces was tuned according to ℎ 5.3 . Using this convection coefficient, the model was adjusted to evaluate the thermal reduction obtained in the case of applying this coating solution over the exterior surface of the canvas. Figure 8 represents the experimental temperature using the blue line, whereas the temperature simulated with the external convection coefficient of the coating solution is represented with red line. The simulation shows that by applying the coating solution, a general thermal reduction of 5.4% can be achieved during the simulation period. However, in individual days, a reduction of 9.4% was obtained. According to the simulations, the use of coating solutions is not enough to achieve results below the comfort limits; however, it is a good strategy to reduce the cooling requirements. In this simulation, the internal simulated temperature was only a 8.12% above the maximum temperature allowed by the comfort standard, which equates to a reduction of 17% of the overheating experienced in normal conditions.

Second Strategy: Reinforcing the Structure
The second strategy used to reduce cooling requirements includes a wood reinforcement in each gap between the wooden structure. Figure 9 shows the position of the wood reinforcements. The reinforcements were in the form of equilateral triangles with an area of 6235.4 using 120cm as the length of each side. Additionally, the thickness of each reinforcement used was 5mm. According to the simulations, the use of coating solutions is not enough to achieve results below the comfort limits; however, it is a good strategy to reduce the cooling requirements. In this simulation, the internal simulated temperature was only a 8.12% above the maximum temperature allowed by the comfort standard, which equates to a reduction of 17% of the overheating experienced in normal conditions.

Second Strategy: Reinforcing the Structure
The second strategy used to reduce cooling requirements includes a wood reinforcement in each gap between the wooden structure. Figure 9 shows the position of the wood reinforcements. The reinforcements were in the form of equilateral triangles with an area of A tm = 6235.4 m 2 using 120 cm as the length of each side. Additionally, the thickness of each reinforcement used was L tm = 5 mm. According to the reinforcement incorporated to the dome, the resistance is modified to include the convection and conduction process, as shown inEquation (14). In this expression, ℎ and ℎ are the internal and external convection coefficients, which in this case were adjusted to ℎ 2.7 and ℎ 3.1 , respectively.
The mathematical model was simulated using the new , and the reduction in internal temperature was 11.8%. Figure 10 represents the experimental data with a blue line, and the simulation results using a red line. The wood reinforcements worked as a new isolation layer, reducing the heat coming into internal air, and allowing for a reduction of the cooling requirements. The simulation with reinforcements showed that the internal temperature was only 6.26% above the upper limit of the comfort zone according to the 55 standards [49].  According to the reinforcement incorporated to the dome, the resistance R 3 is modified to include the convection and conduction process, as shown in Equation (14). In this expression, h te and h ti are the internal and external convection coefficients, which in this case were adjusted to h te = 2.7 W m 2 K and h ti = 3.1 W m 2 K , respectively.
The mathematical model was simulated using the new R 3 , and the reduction in internal temperature was 11.8%. Figure 10 represents the experimental data with a blue line, and the simulation results using a red line. The wood reinforcements worked as a new isolation layer, reducing the heat coming into internal air, and allowing for a reduction of the cooling requirements. The simulation with reinforcements showed that the internal temperature was only 6.26% above the upper limit of the comfort zone according to the 55 standards [49]. According to the reinforcement incorporated to the dome, the resistance is modified to include the convection and conduction process, as shown inEquation (14). In this expression, ℎ and ℎ are the internal and external convection coefficients, which in this case were adjusted to ℎ 2.7 and ℎ 3.1 , respectively.
The mathematical model was simulated using the new , and the reduction in internal temperature was 11.8%. Figure 10 represents the experimental data with a blue line, and the simulation results using a red line. The wood reinforcements worked as a new isolation layer, reducing the heat coming into internal air, and allowing for a reduction of the cooling requirements. The simulation with reinforcements showed that the internal temperature was only 6.26% above the upper limit of the comfort zone according to the 55 standards [49].  The inclusion of new surfaces in the geodesic dome envelope increases the thermal resistance and reduces the heat flux entering the thermal zone. This produces an effect of thermal inertia, in which the dome's internal temperature is reduced during the day. However, at night, it remains a little higher than the ambient temperature.

Conclusions
The study of thermal zones involves the improvement of previous models. In this work, the mathematical model of a geodesic dome was readjusted to achieve a percentage difference of about 3.2% in comparison with experimental data. The geodesic dome model was studied using a Monte Carlo (MC) sensitivity analysis to determine the factors with a major impact on the internal temperature. The MC analysis showed that the convection coefficients were the most important factors influencing the thermodynamic process and therefore needed be limited using passive strategies. Two different passive strategies were tested in our simulations. The first consisted of applying a coating solution to the exterior faces of the structure. This solution provided a general reduction on the internal temperature of around 5.4%. The second strategy consisted of reinforcing the gaps left by the wooden structure. The model was tuned to simulate the new element, showing a reduction of 11.8% in the internal temperature. The passive strategies simulated in this paper were insufficient to guarantee thermal comfort in the geodesic dome, but were proven to contribute to the cooling requirements on the order of 17% to 19%. Both strategies are viable alternatives to reducing the energy consumption of an HVAC system. In future works, the results obtained in this investigation will be verified experimentally, by applying one or both passive strategies in order to reduce refrigeration requirements.

Conflicts of Interest:
The authors declare no conflict of interest.