Numerical Assessment of Zebra-Stripes-Based Strategies in Buildings Energy Performance: A Case Study under Tropical Climate

Urban growth has increased the risk of over-heating both in the microclimate and inside buildings, affecting thermal comfort and energy efficiency. That is why this research aims to evaluate the energy performance of buildings in terms of thermal comfort (operative temperature (OP) levels, satisfied hours of natural ventilation SHNV, thermal lag), and energy efficiency (roof heat gains and surface temperatures) in an urban area in Panama City, using superficial-heat-dissipation biomimetic strategies. Two case studies, a base case and a proposed case, were evaluated using the Designbuilder software through dynamic simulation. The proposed case is based on a combined biomimetic strategy; the reflective characteristics of the Saharan ant applied as a coating on the roofs through a segmented pattern such as the Zebra’s stripes (one section with coating, and another without). Results showed that the OP decreased from 8 to 10 °C for the entire urban zone throughout the year. A reduction of 3.13% corresponding to 8790 kWh per year was achieved for cooling energy consumption. A difference of 5 °C in external surface temperature was obtained, having a lower temperature in which the biomimetic strategy was applied. Besides, it was evidenced that a contrasted-reflectivity-stripes pitched roof performed better than a fully reflective roof. Thus, the functionality of Zebra stripes, together with the reflective characteristics of the Saharan ant, provide better performance for buildings’ thermal regulation and energy needs for cooling.


Introduction
In the last century, there has been a high growth of urban areas where deforestation of green areas was not considered, and materials were used in buildings with high thermal retention, soil occupation, and high levels of greenhouse gas emissions. The greenhouse effect is on the other hand; urban areas consume between 60% and 80% of the energy produced, which today has generated a high environmental and energy cost for the planet [1].
Such consequences are reflected in the effect of urban heat islands (UHI) in cities that cause high rates of heat and temperature in urban areas and consequently affect people's comfort both outside and inside buildings and in the high rates of energy consumption of a part of the buildings. Latin America and the Caribbean, according to United Nations studies in 2018, had one of the highest percentages of urban population with around 81%, a trend which will not decrease [2]. Panama is no exception-there is a high uncertainty regarding moregulation mechanisms for hot and cold climates and the models found in nature. Such as in termite mounds, the heat exchange of the prickly pear, and even humans themselves when sweating. Under this concept, an evaporative cooling system for building envelopes (Stoma Brick) was designed, based on four parts that make it up: in which the stoma brick mainly controls the entry and exit of moisture and at the same time retains it for evaporation. The design is adaptable for hot, cold, humid, and dry climates, mainly simulating a plant's stomata, conifers, and human skin.
Moreover, other strategies based on skin characteristics can be found in nature for hot climates, such as the reflective properties of the Saharan silver ant's hair [14], the head position of the Western Reef Heron [15], and the Shell in Snails and slugs [16]. Regarding the reflective strategy for heat dissipation, such as the Saharan ant [17] and the Shell in Snails [18], there are studies that focus on the application of reflective coatings on ceilings and pavements to cool down the surrounding air temperature in urban areas for urban heat island effect mitigation [19,20], a strategy called "cool roofs." A study carried out in [21], based on simulations for cities in California with a Mediterranean climate, where the albedo of all surfaces, both pavement and ceilings, was increased by 0.2, obtaining benefits of up to −0.01$ per year per square meter of energy savings, and CO 2 reductions of less than 1 kg per square meter. However, this type of strategy depends a lot on the climatic zone and also has its disadvantages with respect to the application on floors or ceilings, since it can affect the visibility of people and even increase the temperature of the surfaces of some buildings of greater height to others of less height [22]. Another study applied a reflective coating using light tones combined with phase change materials (PCM) to reduce the monthly cooling energy consumption. Being in a tropical area-Singapore-a reduction of cooling energy between 5% to 12% was achieved throughout the year [23]; supporting results were later obtained in [24].
The thermoregulatory functions associated with the stripes in Zebras' skin also fall within such type of strategy. However, by mimicking the animal body via colored metal barrels, experimental evidence showed no significant difference in the core temperature of barrels covered with horse and Zebra skins with different scratch patterns. Thermographic measurements were employed. It was concluded that the stripes' coloring and patterns do not statistically significantly influence the body thermoregulation in Zebras [25]. In contrast, a study carried out in Kenya also performs thermographic experiments on live Zebras for seven hours a day. It was evidenced that there is a temperature difference of 12 to 15 degrees between contrasted stripes [26], mainly due to sweat that is accelerated by latter, which in turn, with the movement of the air, causes turbulency at the Zebra's hair, increasing cooling by evaporation.
On the other hand, studies have been carried out that involves the biomimetic design in facades specifically applied to insulators, which are simulated in TRNSYS to evaluate the energy reduction potential compared to having no strategy and thus obtaining a significant reduction in the greenhouse gas emissions in the useful life of the building, in which it applies for different types of climates, in different types of buildings [27]. It was recognized that the thermal bioarchitectural framework was valid as a bridge between architecture and biology [28] in order to find designs with efficient thermal performances. In the case of two buildings in New Zealand, in which the buildings were first simulated to recognize the main problems and divide the thermal zones, different scenarios were analyzed in which mainly passive and active techniques were used in combination. The authors concluded that imitating biological forms so far does not seem to be critical to assessing energy efficiency in buildings [29].
It is common to evaluate the potential ideas that biomimetic-based design can offer us. With a solution-based approach [30] these are generated under three main categories (biological domain, transfer, and technological domain) [31], multiple proposals for designs of more efficient buildings. In the review, as a first point, an extensive investigation of all the characteristics and biological behavior is carried out, taking aspects such as: heat control, use of organic material, respiratory control, use of water, among others. This knowledge is transferred and compared with similar processes in construction and engineering. For example, breathing control compared to ductwork, water usage compared to evaporative cooling. Finally, there is the technological domain in which they evaluate which tool or method can be implemented or improved in the design, be it at the level of the building, systems, or components. It was proposed to change load and energy consumption through passive design to retain and release heat, heat sinks emulating the toco-toucan, and develop new heat recovery systems based on aquatic species and products that reduce the range of current detectors by copying the behavior of some insects. The researchers highlighted the importance of bioinspired design and its potential in the future [32].
Finally, while there are multiple proposals focused on the interior conditions of individual buildings, few studies can be found that investigate the biomimetic approach at urban scale via experimental approaches [33][34][35] and even fewer via dynamic simulation [36] at this scale under tropical climates. Based on our particular interest in dynamic simulation and numerical studies, by focusing on the heat island problem mitigation, an extensive investigation of the biological analogies was carried out for a case study of the problem for the Old Town, Casco Antiguo, in Panama City [36]. Here, to evaluate and improve the thermal comfort of pedestrian (outdoor comfort) through the ENVI-met software, considerable changes were made to the roof of the buildings, emulating the physiology of Zebras and the characteristics of the Saharan ant, in addition to changes in the pavement and vegetative growth. The changes made resulted in temperature reductions of up to 4 • C in a specific area of the study and reductions in the Physiological Equivalent Temperature (PET) comfort index indicators for pedestrians. However, this study did not evaluate the impact of such biomimetic strategies in the indoor thermal comfort, envelope performance, and energy consumption.
Therefore, the objective of the present study is based on the research of reference [8], [36]. Such previous studies, performed for the same studied urban zone, concluded the following: 1.
The urban heat island effect, evidenced by the external comfort, is increased by the narrowness of buildings within the urban zone under such a tropical climate [36]; 2.
Natural ventilation (or passive mode) appears not viable, and the extensive use of air conditioning systems (or active mode) may be required to provide acceptable indoor thermal comfort [8]; 3.
Other thermoregulation strategies are needed to improve external comfort [36]; 4.
The use of microclimatic data could significantly influence the estimation of the building energy performance [8], causing about a 10% difference lower for cooling needs than when using standard weather data.
Since [8] evaluated the passive and active modes in the base case study (with no incorporation of biomimicry) indicating that microclimate data should be used for indoor thermal comfort assessment and energy performance, but encountering that the passive mode appeared not viable, the addition of the biomimetic strategies in [36], which improved the microclimate conditions, could increase the viability of the use of natural ventilation, and the buildings' energy performance in this urban zone.
That is why the use of the same biomimetic strategies in [36] are applied to the base case study (namely, proposed case), to assess the thermal and energy performance of the urban zone. The performance is evaluated by means of indoor thermal comfort indicators, such as operating temperature, surface temperature, as well as the evaluation of the satisfied hours of natural ventilation (SHNV), using dynamic simulation through the Designbuilder software v6.1.6.011 [37] and ENVI-met. v4.4.5 [38]. The evaluation was carried out in March (driest) and October (rainiest), comparing the base case (or reference case) and the proposed case. Besides, full reflective coating is applied to the base case and compared with both other cases in terms of performance. Since the combined biomimetic strategy of [36] surpasses in performance the full reflective coating approach, a 2D transient heat transfer model is constructed to clarify this outcome contrasted with those in [25].

Materials and Methods
For the development of this research, the methodology implemented consists of evaluating microclimatic data within an urban zone in the Old Town of Panama City named Casco Antiguo via ENVI-met and the corresponding standard typical meteorological data. The microclimatic data were adapted via a numerical approach to assessing the energy performance of buildings within the urban zone [8]. The Designbuilder software, together with the microclimatic data, were used to assess such performance of buildings operative in passive and active modes for both the base case and the proposed case. The proposed case is based on biomimicry strategies ( Figure 1).
with both other cases in terms of performance. Since the combined biomimetic strategy of [36] surpasses in performance the full reflective coating approach, a 2D transient heat transfer model is constructed to clarify this outcome contrasted with those in [25].

Materials and Methods
For the development of this research, the methodology implemented consists of evaluating microclimatic data within an urban zone in the Old Town of Panama City named Casco Antiguo via ENVI-met and the corresponding standard typical meteorological data. The microclimatic data were adapted via a numerical approach to assessing the energy performance of buildings within the urban zone [8]. The Designbuilder software, together with the microclimatic data, were used to assess such performance of buildings operative in passive and active modes for both the base case and the proposed case. The proposed case is based on biomimicry strategies ( Figure 1).

Case Study and Problem Identification
The study focuses on the Old Town of Panama City, Caso Antiguo, which is declared a World Heritage Site by UNESCO. Its colonial design is based on the laws of the Indies, in which its buildings are mainly composed of materials such as: wood, calicanto, clay blocks, concrete, tile roofs, and concrete slab roofs [1]. The area is classified as a climatic zone type 3 (LCZ3) according to Oke [39], which consists of high-density, low-level buildings (3 to 4 levels), paved ground, and few trees. To perform the numerical study, part of this urban zone was analyzed: the cutout chosen covers a dimension of 290 m (x) by 226 m (y), as can be seen in Figure 2. This urban zone was selected to analyze further previous studies conducted for the same cutout [1,8,36,40]. Table 1 summarizes the materials of the buildings in the studied zone or cutout. Figure 3 shows the 3D model of the base case study area in Designbuilder and Table  2 presents the transmittance values for each of the buildings' envelope elements. Using the Google Earth Pro tool, the dimensions of the streets and buildings were taken from the cutout, and the Designbuilder blocks tool was used to create the buildings. Only the energy consumption and thermal comfort of gray-colored buildings are considered in the simulation results. The uninhabited buildings, churches, and historical monuments are shown in color pink, which indicates that their individual energy consumption and thermal comfort are not considered in the simulation results but their impact on other buildings performance is considered, through shading, wind blocking, and reflectiveness. In

Case Study and Problem Identification
The study focuses on the Old Town of Panama City, Caso Antiguo, which is declared a World Heritage Site by UNESCO. Its colonial design is based on the laws of the Indies, in which its buildings are mainly composed of materials such as: wood, calicanto, clay blocks, concrete, tile roofs, and concrete slab roofs [1]. The area is classified as a climatic zone type 3 (LCZ3) according to Oke [39], which consists of high-density, low-level buildings (3 to 4 levels), paved ground, and few trees. To perform the numerical study, part of this urban zone was analyzed: the cutout chosen covers a dimension of 290 m (x) by 226 m (y), as can be seen in Figure 2. This urban zone was selected to analyze further previous studies conducted for the same cutout [1,8,36,40]. Table 1 summarizes the materials of the buildings in the studied zone or cutout. Figure 3 shows the 3D model of the base case study area in Designbuilder and Table 2 presents the transmittance values for each of the buildings' envelope elements. Using the Google Earth Pro tool, the dimensions of the streets and buildings were taken from the cutout, and the Designbuilder blocks tool was used to create the buildings. Only the energy consumption and thermal comfort of gray-colored buildings are considered in the simulation results. The uninhabited buildings, churches, and historical monuments are shown in color pink, which indicates that their individual energy consumption and thermal comfort are not considered in the simulation results but their impact on other buildings performance is considered, through shading, wind blocking, and reflectiveness. In addition to the information related to the heights of the buildings, the construction materials according to the real characteristics of the pavements, surfaces, and vegetation (only trees without the consideration physiological functions, such as evapotranspiration) were placed.
Biomimetics 2022, 7, x FOR PEER REVIEW 6 of 21 addition to the information related to the heights of the buildings, the construction materials according to the real characteristics of the pavements, surfaces, and vegetation (only trees without the consideration physiological functions, such as evapotranspiration) were placed.

Concrete block
Biomimetics 2022, 7, x FOR PEER REVIEW 6 of 21 addition to the information related to the heights of the buildings, the construction materials according to the real characteristics of the pavements, surfaces, and vegetation (only trees without the consideration physiological functions, such as evapotranspiration) were placed.
(a) (b)  addition to the information related to the heights of the buildings, the construction materials according to the real characteristics of the pavements, surfaces, and vegetation (only trees without the consideration physiological functions, such as evapotranspiration) were placed.
(a) (b)  addition to the information related to the heights of the buildings, the construction materials according to the real characteristics of the pavements, surfaces, and vegetation (only trees without the consideration physiological functions, such as evapotranspiration) were placed.
(a) (b)  addition to the information related to the heights of the buildings, the construction materials according to the real characteristics of the pavements, surfaces, and vegetation (only trees without the consideration physiological functions, such as evapotranspiration) were placed.
(a) (b)    addition to the information related to the heights of the buildings, the construction materials according to the real characteristics of the pavements, surfaces, and vegetation (only trees without the consideration physiological functions, such as evapotranspiration) were placed.

Proposed Designs and Simulation
An evaluation of the study area previously proposed, in which applied biomimetic strategies focused on a problem-based approach, i.e., the basis of an exhaustive search for biological analogies, where the most related were extracted thanks to the biomimetic design methodology [36]. Such biological analogies were applied on the buildings' roofs to evaluate their effects on pedestrians' comfort. Now, the present study concentrates on evaluating their effects on indoor thermal comfort and energy consumption.

Abstraction and Emulation of the Identified Biomimicry-Based Strategies
The biomimicry problem-based approach applied to the case study led to abstract and emulated a two-pinnacle strategy based on the Zebra stripes and the Saharan ant. The Saharan ant for its highly reflective characteristics, and the zebra specifically its black and white appearance, to emulate convective currents and achieve evaporative cooling ( Figure 4). Table 3 shows the summary of such pinnacles' analysis [36]. Table 3. Summary of the Pinnacles Analysis. Own elaboration.

Pinnacles Strategy Mechanism Fundamental Principles Main Feature
Zebra Black and white streaking causes a temperature differential [30,41].
Convective currents are caused by increasing evaporation. High convection and evaporation Arrangement of animal stripes for heat regulation Saharan ant High reflection in the NIR range and emissivity in the NIR [17].
Reflection of thermal radiation and high emissivity to release excess heat High reflectivity and emissivity Silver hairs with triangular structure The emulation of such pinnacle strategies was implemented in the urban zone, as in [36] but in the Designbuilder software on the buildings' roofs. For this, an additional layer of a reflective coating was added to the based-case pitched roof construction (Table 4), following the Zebra-stripes pattern interchanging between the reflective coating and the original roof, as shown in Figure 4. A change in the pitched roof geometry was made from a multi-slopes roof to a two-slope roof for simplicity to adjust and include the stripes pattern. The reflective lines added are 1 cm thick and the width depended on the length of each building ranging from 4 to 6 m wide; no uniform stripe width and shape may be key for Zebras [42]. The buildings' roofs with flat roof constructions were kept as in the base case. The 3D model of the proposed case developed in Designbuilder is shown in Figure 5. To evaluate the base case further, the reflective coating was added to all pitched roof constructions (on gray-colored buildings) for comparison purposes to the proposed case. The value of the coating was based on the value of the reflective properties of the Saharan ant 0.97 [43], commercially a high value of reflective coating of 0.92 [44] was found and therefore that value was taken which was the closest. ant 0.97 [43], commercially a high value of reflective coating of 0.92 [44] was found and therefore that value was taken which was the closest.   ant 0.97 [43], commercially a high value of reflective coating of 0.92 [44] was found and therefore that value was taken which was the closest.  The simulations for both cases were performed using microclimatic data extracted from ENVI-met for the whole year. These microclimatic data were created using a numer-

Simulation of the Based and Proposed Cases
The simulations for both cases were performed using microclimatic data extracted from ENVI-met for the whole year. These microclimatic data were created using a numerical approach as in [8] and results from simulations carried out in ENVI-met using standard weather data (typical meteorological data) obtained from CLIMdata Solargis © ( Table 5). The resulting microclimatic data for the air temperature is presented in Figure 6.   The main activity in the Casco Antiguo is commercial, such as restaurants and bars, equivalent to 50% of the buildings, another 40% is dedicated to hotels and residences, and 10% are offices and public institutions. The cut-out of the old town has a different occupation profile from another area of the city [8]. The main activity in the Casco Antiguo is commercial, such as restaurants and bars, equivalent to 50% of the buildings, another 40% is dedicated to hotels and residences, and 10% are offices and public institutions. The cut-out of the old town has a different occupation profile from another area of the city [8].
The town of San Felipe has a population of approximately 3262 people and is where the Casco Antiguo is located. The population density is presented in [8] (Figure 7), which was considered equivalent to the occupation profile carried out in the simulation (Table 6). The town of San Felipe has a population of approximately 3262 people and is where the Casco Antiguo is located. The population density is presented in [8] (Figure 7), which was considered equivalent to the occupation profile carried out in the simulation (Table  6).

Buildings' Performance Evaluation and Comparison
The buildings' performance was evaluated in terms of the following indicators: indoor thermal comfort, roof heat gains, and the electricity consumption for cooling. This evaluation was performed at the building level, i.e., the average values for the entire studied urban zone. The resulting values for each indicator were analyzed and compared in monthly and hourly averages, depending on the indicator. The urban zone was subjected to two operation modes: passive and active modes. The latter mode corresponds to an operation under the use of air conditioning, and the former mode corresponds to an operation under natural ventilation only. This was completed beforehand for both the base case and proposed case.

Passive Mode Operation
The building thermal performance during passive operation mode was assessed using the following indicators: Operative temperature, roof heat gains, and the Satisfied hours of natural ventilation (SHNV). Results are obtained only for occupied periods for monthly indicators. The SHNV is introduced here as it assesses the possibility of ensuring Figure 7. Population density in the urban zone (adapted from [8]). The studied zone inside black dashed square. Table 6. Occupied periods and energy usages in the urban zone [8].

Buildings' Performance Evaluation and Comparison
The buildings' performance was evaluated in terms of the following indicators: indoor thermal comfort, roof heat gains, and the electricity consumption for cooling. This evaluation was performed at the building level, i.e., the average values for the entire studied urban zone. The resulting values for each indicator were analyzed and compared in monthly and hourly averages, depending on the indicator. The urban zone was subjected to two operation modes: passive and active modes. The latter mode corresponds to an operation under the use of air conditioning, and the former mode corresponds to an operation under natural ventilation only. This was completed beforehand for both the base case and proposed case.

Passive Mode Operation
The building thermal performance during passive operation mode was assessed using the following indicators: Operative temperature, roof heat gains, and the Satisfied hours of natural ventilation (SHNV). Results are obtained only for occupied periods for monthly indicators. The SHNV is introduced here as it assesses the possibility of ensuring acceptable indoor air quality and natural thermal comfort. The latter corresponds to the percentage of the maximum number of hours when the outdoor climate is favorable for natural ventilation inside buildings, based on the thermal comfort requirements established in the ASHRAE 55 standard for naturally conditioned areas [45]. This indicator can be significantly affected by the outdoor climate, effects of building characteristics such as design, internal heat gains, and envelope energy performance. Finally, the criteria for the SHNV indicator are defined based on the average outdoor air temperature and the operative temperature, as the intersection of these two temperatures for an acceptable 80% ambient. If the thermal comfort condition does not meet the established criteria, this indicator equals zero.

Active Mode Operation
The building thermal performance during active operation mode was assessed by focusing on the electricity consumption for cooling purposes, only for occupied periods. The cooling temperature was set to 24 • C considering air conditioning devices with a coefficient of performance of 3.00. Natural ventilation was permitted during unoccupied periods.

Heat Transfer Analysis of the Zebra-Stripe Strategy
Other studies concerning the Zebra striped skin have associated this strategy with the survival way to regulate skin and body temperature. Experimental studies evidenced a significant temperature difference between white and black stripes. Thus, in order to further understand the building energy performance simulation results and for complementary purposes, 2D heat transfer transient simulations were carried out via the Energy2D software [46] based on Designbuilder results. The 2D heat transfer model setup is presented in Figure 8.
acceptable indoor air quality and natural thermal comfort. The latter corresponds to the percentage of the maximum number of hours when the outdoor climate is favorable for natural ventilation inside buildings, based on the thermal comfort requirements established in the ASHRAE 55 standard for naturally conditioned areas [45]. This indicator can be significantly affected by the outdoor climate, effects of building characteristics such as design, internal heat gains, and envelope energy performance. Finally, the criteria for the SHNV indicator are defined based on the average outdoor air temperature and the operative temperature, as the intersection of these two temperatures for an acceptable 80% ambient. If the thermal comfort condition does not meet the established criteria, this indicator equals zero.

Active Mode Operation
The building thermal performance during active operation mode was assessed by focusing on the electricity consumption for cooling purposes, only for occupied periods. The cooling temperature was set to 24 °C considering air conditioning devices with a coefficient of performance of 3.00. Natural ventilation was permitted during unoccupied periods.

Heat Transfer Analysis of the Zebra-Stripe Strategy
Other studies concerning the Zebra striped skin have associated this strategy with the survival way to regulate skin and body temperature. Experimental studies evidenced a significant temperature difference between white and black stripes. Thus, in order to further understand the building energy performance simulation results and for complementary purposes, 2D heat transfer transient simulations were carried out via the En-ergy2D software [46] based on Designbuilder results. The 2D heat transfer model setup is presented in Figure 8. To simplify the transient heat transfer simulation and since the Energy2D software solves the heat conduction equation together with the conservation laws, the following conditions were considered: 1. The thermal boundary conditions for the outside were considered as the same as the air temperature; this was set to 0 °C (default temperature value for any surrounding in the software), except for the bottom boundary, which was set to the same temperature as the body. The mass boundaries were set as open. The thermal properties corresponded to the set-up temperature; 2. The body may represent the body of the Zebra or the roof construction material just below the coated surface. Its temperature was set to 50 °C throughout the entire simulation; 3. The outside airflow speed remained constant at 0.01 m/s (using a default ventilator).
The time step was set to 10 s, and the simulation duration was limited to 30 min; To simplify the transient heat transfer simulation and since the Energy2D software solves the heat conduction equation together with the conservation laws, the following conditions were considered: 1.
The thermal boundary conditions for the outside were considered as the same as the air temperature; this was set to 0 • C (default temperature value for any surrounding in the software), except for the bottom boundary, which was set to the same temperature as the body. The mass boundaries were set as open. The thermal properties corresponded to the set-up temperature; 2.
The body may represent the body of the Zebra or the roof construction material just below the coated surface. Its temperature was set to 50 • C throughout the entire simulation; 3.
The outside airflow speed remained constant at 0.01 m/s (using a default ventilator). The time step was set to 10 s, and the simulation duration was limited to 30 min; 4.
The stripes were also kept at a constant temperature value. The original roof layer was set to the same temperature as the body, and the high-reflective coating was set to 10 • C; 5.
The values for the thermal conductivity, specific heat capacity, and density remained as the default values: 1 W/mK, 1300 J/kgK, and 25 kg/m 3 ; 6.
Any radiation heat transfer interaction was neglected since the stripes' temperatures were set at a constant value, to consider it in advance.
This setup does not intend to replicate the exact situation considered for the energy performance simulations, but it may serve as a reference setup to evidence and understand the heat transfer problem. Finally, since a finite-difference method is used to solve the equations, all physical contact among the stripes and body are considered.

Results Analysis
This section presents the analysis of the results and their discussion, starting with evaluating the indoor thermal performance via operative temperature and SHNV levels. A brief analysis of the thermal lag and damping factor results is presented. This is followed by the energy efficiency evaluation in terms of roof surface temperatures and heat gains and the electricity consumption for cooling. Finally, the 2D transient heat transfer simulation results are presented.

Evaluation of the Thermal Performance via Operative Temperature and SHNV
The average operative temperature results for the entire urban zone can be seen in Figure 9. It can be clearly seen that the proposed case presented a significant reduction concerning the base case, reaching 8.51 • C of reduction for March (from 38.28 • C to 29.77 • C) after applying the Zebra-stripes strategies. The same happened in October with a 9.59 • C reduction. This reduction presented an average of approximately 8.82 • C (from 37.66 • C to 28.84 • C) throughout the year. The main reason for this reduction is associated with the application of the reflective layer, which prevents a large part of the solar radiation that falls on the building from penetrating inside it, which makes the operative temperature lower. Results also showed that when implementing the reflective coating for all pitched roofs in the base case, the monthly operative temperature values obtained are higher than in the proposed case.
4. The stripes were also kept at a constant temperature value. The original roof layer was set to the same temperature as the body, and the high-reflective coating was set to 10 °C; 5. The values for the thermal conductivity, specific heat capacity, and density remained as the default values: 1 W/mK, 1300 J/kgK, and 25 kg/m 3 ; 6. Any radiation heat transfer interaction was neglected since the stripes' temperatures were set at a constant value, to consider it in advance.
This setup does not intend to replicate the exact situation considered for the energy performance simulations, but it may serve as a reference setup to evidence and understand the heat transfer problem. Finally, since a finite-difference method is used to solve the equations, all physical contact among the stripes and body are considered.

Results Analysis
This section presents the analysis of the results and their discussion, starting with evaluating the indoor thermal performance via operative temperature and SHNV levels. A brief analysis of the thermal lag and damping factor results is presented. This is followed by the energy efficiency evaluation in terms of roof surface temperatures and heat gains and the electricity consumption for cooling. Finally, the 2D transient heat transfer simulation results are presented.

Evaluation of the Thermal Performance via Operative Temperature and SHNV
The average operative temperature results for the entire urban zone can be seen in Figure 9. It can be clearly seen that the proposed case presented a significant reduction concerning the base case, reaching 8.51 °C of reduction for March (from 38.28 °C to 29.77 °C) after applying the Zebra-stripes strategies. The same happened in October with a 9.59 °C reduction. This reduction presented an average of approximately 8.82 °C (from 37.66 °C to 28.84 °C) throughout the year. The main reason for this reduction is associated with the application of the reflective layer, which prevents a large part of the solar radiation that falls on the building from penetrating inside it, which makes the operative temperature lower. Results also showed that when implementing the reflective coating for all pitched roofs in the base case, the monthly operative temperature values obtained are higher than in the proposed case.  Subsequently, an inspection throughout the urban zone allowed noticing that only a two-story building presented acceptable values for the SHNV calculation (selected building in red framed in Figure 5); the same building as the base case. Note here that Figure 9 shows the operative temperature results for the entire urban zone, while Figure 10 shows the operative temperature for the selected building. Temperature reductions can be observed for six months.
Subsequently, an inspection throughout the urban zone allowed noticing that only a two-story building presented acceptable values for the SHNV calculation (selected building in red framed in Figure 5); the same building as the base case. Note here that Figure 9 shows the operative temperature results for the entire urban zone, while Figure 10 shows the operative temperature for the selected building. Temperature reductions can be observed for six months. Accounting only for the selected building, the results of the operative temperature for the critical day of each month were analyzed in detail, assuming this behavior for the rest of the month: 20 February, 20 October, and 11 November. These dates were chosen because they presented acceptable values for the SHNV in the base case. Hourly temperature results indicate that the maximum temperature is no longer reached as specified by the standard weather data (Table 4); maximum values are reached at 16 h and 17 h for outdoor air temperature and operative temperature, respectively. This evidences the impact of the use of microclimatic data. However, for the proposed case with these operative temperature results and the already recorded outdoor air temperature values (Figure 6), only October presented acceptable values for the evaluation of the SHNV. The corresponding SHNV for October resulted in 37.5% (279 h of the 744 h), which conduce to a 3.18% for the entire year. For the base case, the SHNV for February, October, and November resulted in 4.17% (28 h), 16.67% (124 h), and 16.67% (120 h), respectively, conducing to 3.11% for the entire year. These results confirm the thermal unconformity recorded inside the buildings during the occupancy period of the buildings in passive mode, and such values are considered significantly low according to Causone in [47].
Moreover, two parameters that indicate the envelope's thermal performance were evaluated between the proposed case and the base case shown in the table. Such is the thermal lag, which refers to the time it takes to transfer energy in the form of heat from the exterior to the interior in hours, and the damping, which is a fraction between the difference in max and min temperatures of the interior and exterior of the building [48]. In Table 7, it can be observed that according to the thermal lag in November, the operative temperature reached its peak value (or maximum value) faster in the proposed case than in the base case. No significant changes were observed for other months. Therefore, there is a greater heat transfer penetration in the proposed case than in the base case. 32  Accounting only for the selected building, the results of the operative temperature for the critical day of each month were analyzed in detail, assuming this behavior for the rest of the month: 20 February, 20 October, and 11 November. These dates were chosen because they presented acceptable values for the SHNV in the base case. Hourly temperature results indicate that the maximum temperature is no longer reached as specified by the standard weather data (Table 4); maximum values are reached at 16 h and 17 h for outdoor air temperature and operative temperature, respectively. This evidences the impact of the use of microclimatic data. However, for the proposed case with these operative temperature results and the already recorded outdoor air temperature values (Figure 6), only October presented acceptable values for the evaluation of the SHNV. The corresponding SHNV for October resulted in 37.5% (279 h of the 744 h), which conduce to a 3.18% for the entire year. For the base case, the SHNV for February, October, and November resulted in 4.17% (28 h), 16.67% (124 h), and 16.67% (120 h), respectively, conducing to 3.11% for the entire year. These results confirm the thermal unconformity recorded inside the buildings during the occupancy period of the buildings in passive mode, and such values are considered significantly low according to Causone in [47].
Moreover, two parameters that indicate the envelope's thermal performance were evaluated between the proposed case and the base case shown in the table. Such is the thermal lag, which refers to the time it takes to transfer energy in the form of heat from the exterior to the interior in hours, and the damping, which is a fraction between the difference in max and min temperatures of the interior and exterior of the building [48]. In Table 7, it can be observed that according to the thermal lag in November, the operative temperature reached its peak value (or maximum value) faster in the proposed case than in the base case. No significant changes were observed for other months. Therefore, there is a greater heat transfer penetration in the proposed case than in the base case.

External Surface Temperature for the Proposed Case
The external surface temperature results for 17 March and 20 October, at 15 h, are presented in Figures 11 and 12, for the base case and proposed case, respectively. As it is 3:00 p.m., it can be seen that the temperature is higher on the west side of the roofs, in March the vast majority of the roofs oscillated 55-60 • C and the east side 45-50 • C, by October the west side of the roof oscillates between 60-65 • C and the east side between 45-50 • C-therefore for the base case in October a higher temperature is obtained.

External Surface Temperature for the Proposed Case
The external surface temperature results for 17 March and 20 October, at 15 h, are presented in Figure 11 and Figure 12, for the base case and proposed case, respectively. As it is 3:00 p.m., it can be seen that the temperature is higher on the west side of the roofs, in March the vast majority of the roofs oscillated 55-60 °C and the east side 45-50 °C, by October the west side of the roof oscillates between 60-65 °C and the east side between 45-50 °C-therefore for the base case in October a higher temperature is obtained. As expected, significant differences between the reflective coating and the origi roof construction surface temperature can be observed. Similarly, low surface tempe ture values were obtained when implementing the reflective coating for all pitched ro in the base case. Regarding the monthly average surface temperature values, for Mar the difference between the roof of the same building using the Zebra-stripes biomime strategy (ZSBS) oscillates in a range of 27.5 °C and 30 °C in the reflective segment and the non-reflective segment between 32.5 °C and 35 °C, that is, a difference of 5 °C betwe both segments. For flat buildings in which the ZSBS were not applied, the temperatu difference of these buildings concerning those with strategies is about 5 °C.
For October, the difference between the roof of the same building having the ZS oscillates in a range of 25-27.5 °C in the reflective segment and in the non-reflective s ment between 32.5 °C and 35 °C, corresponding to a difference of 7.5 °C between bo As expected, significant differences between the reflective coating and the original roof construction surface temperature can be observed. Similarly, low surface temperature values were obtained when implementing the reflective coating for all pitched roofs in the base case. Regarding the monthly average surface temperature values, for March, the difference between the roof of the same building using the Zebra-stripes biomimetic strategy (ZSBS) oscillates in a range of 27.5 • C and 30 • C in the reflective segment and in the non-reflective segment between 32.5 • C and 35 • C, that is, a difference of 5 • C between both segments. For flat buildings in which the ZSBS were not applied, the temperature difference of these buildings concerning those with strategies is about 5 • C.
For October, the difference between the roof of the same building having the ZSBS oscillates in a range of 25-27.5 • C in the reflective segment and in the non-reflective segment between 32.5 • C and 35 • C, corresponding to a difference of 7.5 • C between both segments. In buildings not having the ZSBS, the temperature difference concerning those that have the ZSBS is 3.75 • C. For March, considering the complete area of the model, the proposed case decreased 7.5 • C with respect to the base case. The same happened for October. Figure 13 shows the monthly averages for the roof-ceiling heat gains for the urban zone during the occupied hours for the base and proposed cases. Besides, the base case when implementing the reflective coating on all pitched roof constructions is shown. As expected, the lowest heat gains are encountered for the base case with fully reflective coating on all pitched roofs. The proposed case presented slight reductions concerning the base case. Figure 13 shows the monthly averages for the roof-ceiling heat gains for the urban zone during the occupied hours for the base and proposed cases. Besides, the base case when implementing the reflective coating on all pitched roof constructions is shown. As expected, the lowest heat gains are encountered for the base case with fully reflective coating on all pitched roofs. The proposed case presented slight reductions concerning the base case.

Comparison of Roof-Ceiling Heat Gains and Cooling Electricity Consumption
On the other hand, the electricity consumption for cooling needs can be observed in Figure 14. Slight differences are encountered between the base and proposed cases, with an average difference of 3.13%. The latter corresponds to an annual reduction of 8790 kWh, which is significant when considering electricity costs, but this should be addressed further contrasted with the retrofit costs to implement the proposed case at a real scale. The active operation mode was not simulated when implementing the reflective coating on all pitched roofs for the base case. On the other hand, the electricity consumption for cooling needs can be observed in Figure 14. Slight differences are encountered between the base and proposed cases, with an average difference of 3.13%. The latter corresponds to an annual reduction of 8790 kWh, which is significant when considering electricity costs, but this should be addressed further contrasted with the retrofit costs to implement the proposed case at a real scale. The active operation mode was not simulated when implementing the reflective coating on all pitched roofs for the base case.
segments. In buildings not having the ZSBS, the temperature difference concerning those that have the ZSBS is 3.75 °C. For March, considering the complete area of the model, the proposed case decreased 7.5 °C with respect to the base case. The same happened for October. Figure 13 shows the monthly averages for the roof-ceiling heat gains for the urban zone during the occupied hours for the base and proposed cases. Besides, the base case when implementing the reflective coating on all pitched roof constructions is shown. As expected, the lowest heat gains are encountered for the base case with fully reflective coating on all pitched roofs. The proposed case presented slight reductions concerning the base case.

Comparison of Roof-Ceiling Heat Gains and Cooling Electricity Consumption
On the other hand, the electricity consumption for cooling needs can be observed in Figure 14. Slight differences are encountered between the base and proposed cases, with an average difference of 3.13%. The latter corresponds to an annual reduction of 8790 kWh, which is significant when considering electricity costs, but this should be addressed further contrasted with the retrofit costs to implement the proposed case at a real scale. The active operation mode was not simulated when implementing the reflective coating on all pitched roofs for the base case.

Heat Transfer Analysis of the Zebra-Stripes Strategy
Considering that the proposed case has presented significant results with respect to the base case and base case with fully reflective pitched roof construction. The 2D transient heat transfer simulation results when implementing the ZSBS are presented in Figure 15 (case A). Gray color flashes indicate the temperature range where the highest temperature zones (50 • C) correspond to the clearer zones (in white) and the lowest temperature zones (0 • C) to the darkest. The blue lines represent the heat flux lines. These results were compared with the case where the reflective coating covered the entire surface (case B) in terms of a temperature measurement at the center of the "body". At the end of the simulation time, this temperature measurement reported a lower value for case A than case B, with a difference of 3.8 • C. The graph in Figure 15 shows that in the one with total reflectivity (case B, orange line), the temperature decreases more slowly than the one with the biomimetic strategy (case A, blue line). In more detail, case A shows distortions of the heat flow lines around the contrasting segments, which indicates that the temperature differential in the segments increases heat dissipation, explaining its speed compared to case B.

Heat Transfer Analysis of the Zebra-Stripes Strategy
Considering that the proposed case has presented significant results with respect to the base case and base case with fully reflective pitched roof construction. The 2D transient heat transfer simulation results when implementing the ZSBS are presented in Figure 15 (case A). Gray color flashes indicate the temperature range where the highest temperature zones (50 °C) correspond to the clearer zones (in white) and the lowest temperature zones (0 °C) to the darkest. The blue lines represent the heat flux lines. These results were compared with the case where the reflective coating covered the entire surface (case B) in terms of a temperature measurement at the center of the "body". At the end of the simulation time, this temperature measurement reported a lower value for case A than case B, with a difference of 3.8 °C. The graph in Figure 15 shows that in the one with total reflectivity (case B, orange line), the temperature decreases more slowly than the one with the biomimetic strategy (case A, blue line). In more detail, case A shows distortions of the heat flow lines around the contrasting segments, which indicates that the temperature differential in the segments increases heat dissipation, explaining its speed compared to case B.

Discussion
A proposed model was made based on the emulation of the characteristics of the Saharan ant employing a high reflectivity coating applied to the buildings' roofs and the contrasted stripes of the Zebra by applying the reflective coating in a segmented distribution. All this is using the Designbuilder software in which the microclimatic data were taken for the dynamic simulation carried out in the ENVI-met software. From the simulation results obtained, the following could be highlighted:

•
In the proposed case, the indoor operative temperature of the entire urban zone under study was reduced by about 8-10 °C on average for the year, compared to the base case, due to the application of the reflective coating on the roof, preventing heat gain in inside;

•
Regarding the external surface temperature of a building with biomimetic strategies for March, the reflective and non-reflective segments presented a temperature difference of about 5 °C, and it was also found that the building with strategies is 5 °C colder than a building without the strategies. In October, the temperature difference

Discussion
A proposed model was made based on the emulation of the characteristics of the Saharan ant employing a high reflectivity coating applied to the buildings' roofs and the contrasted stripes of the Zebra by applying the reflective coating in a segmented distribution. All this is using the Designbuilder software in which the microclimatic data were taken for the dynamic simulation carried out in the ENVI-met software. From the simulation results obtained, the following could be highlighted:

•
In the proposed case, the indoor operative temperature of the entire urban zone under study was reduced by about 8-10 • C on average for the year, compared to the base case, due to the application of the reflective coating on the roof, preventing heat gain in inside; • Regarding the external surface temperature of a building with biomimetic strategies for March, the reflective and non-reflective segments presented a temperature difference of about 5 • C, and it was also found that the building with strategies is 5 • C colder than a building without the strategies. In October, the temperature difference between the reflective and non-reflective segments was 7.5 • C, while buildings with strategies have 3.75 • C less (colder) without applied strategies; • The interior temperature for damping for November in the proposed case reached its peak value faster than the base case due to a reduction in the hours of the thermal lag indicator; • Significant reduction in the electricity consumption for cooling was achieved by the proposed case, with respect to the base case, with an average of 3.13% (about 8.8 MWh); • The simplified setup for 2D transient heat transfer simulations evidenced enhanced heat transfer dissipation for the proposed case compared to a fully reflective coating on the surface.
Although the operative temperatures reached lower values than the base case, results from the SHNV indicator confirm the thermal unconformity recorded inside the buildings during the occupancy periods in passive mode, which is considered significantly low just as in [47].
Moreover, the proposed case performed better than the base case and the base case when implementing the reflective coating over all roof constructions in terms of indoor operative temperature values and electricity consumption for cooling. The average reduction in electricity consumption for cooling of 5.07% (about 505 kWh) when using microclimate data as in [8], compared to the average reduction of 3.13% (8.8 MWh), here in the proposed case, might not represent high enough energy savings to undergo the microclimate data generation process from standard weather data or in situ data collection for simulation purposes only, as suggested in [8,49]. However, the use of microclimate data could strengthen the evaluation-verification process of new urban designs [49] before implementation at a real scale by giving more precise values for energy needs forecast. Furthermore, the difference in the roofs with the applied Zebra-stripes biomimetic strategy is comparable as it occurs with the Zebra, which through its black and white coloration, obtains a temperature differential that causes convective currents that accelerate heat dissipation [26]. The effectiveness of this biomimetic strategy is also observed at the surface temperature at the roof west slopes of both cases for both months (Figures 11 and 12). For instance, above the park on the left (surrounding with trees), the buildings' roofs west slope without the reflective coating (color yellow in Figure 12a) appeared 5 • C lower than the same slopes in Figure 11a (in orange). The same can be spotted on the buildings' roofs above the cathedral and those on the righthand side (below the park).
To corroborate, the 2D transient heat transfer simulations helped realize that the heat flow perturbance at the surface of the body (heat flow lines in Figure 15) caused by the stripes allows a faster temperature drop than the surface with homogenous heat flow. However, the impact of this strategy could be enhanced by combining different patterns since the latter may be as fundamental as the stripes' pattern distribution and shape [25]. It is worth highlighting, based on the study [25] that denies the fact that Zebra stripes do not have a significant effect on the decrease in body temperature, our study is based on biomimetics, which does not take the behavior of the analogy directly. Rather, it is adapted to the applications and needs of the design, for example, the reflective coating that also promotes heat dissipation. Finally, using a commercial coating with a reflectivity value of 0.92 [44] on roofs limits the Saharan ant reflectivity value of 0.96, which brings an opportunity for future coating developments. Although significant energy saving for cooling was attained with the combined biomimetic strategy used here, a cost-effective analysis of retrofitting the roofs' construction contrasted with energy saving costs, needs to be performed, as in as in [21], before implementation.

Conclusions
The main objective of this research was to evaluate the indoor thermal comfort and energy efficiency of buildings in an urban area in the Casco Antiguo of Panama City through biomimetic strategies from a previously carried out study. Among the findings this investigation provided, the following can be retained: The functionality of Zebra stripes, together with the reflective characteristics of the Saharan ant, provide better performance for buildings' thermal regulation and energy needs for cooling. Further investigation through experimental studies might support the effectiveness of this combined biomimetic strategy, but a cost-effective analysis of retrofitting the roofs' construction contrasted with costs energy saving needs to be performed before implementation.
It is important to highlight that biomimetics has been applied throughout history by nature without generating negative impacts on the environment. That is why it must be learned from, and we must seek alternatives to solve the challenges currently faced by society through the application of different approaches so that they may be solved sustainably.
A recommendation for future work would be to conduct more in-depth studies verifying the effectiveness of the strategies applied in this study, such as convective currents that influence heat dissipation through experimental tests. The evaluation of the application of biomimetic strategies in walls added to the strategies of this study.