Couple simulations with CFD and ladybug + honeybee tools for green façade optimizing the thermal comfort in a transitional space in hot-humid climate

ABSTRACT Vertical gardens have become important strategies for urban heat island migration in recent years. However, the simulation tools for green façade (GF) with climbing plants coupling the shading and wind-induced effects are still lacked. This study focused on GF effects on human thermal comfort (HTC) optimization in a transitional space. The main purpose was to develop a couple method on HTC simulations. Field measurements in a hot-humid climate area were conducted to obtain the parameters of GF foliage layers. A couple simulation method with Ansys Fluent, Ladybug, and Honeybee tools was developed. Seven GFs layouts on a standard balcony model were simulated. In field measurements, the Ave. global irradiation was reduced by 73.29%, leaf area index was measured at 1.56–3.61, and foliage cover ratio was measured at 49.74–92.88%. In the simulations, the results of the Ave. wind velocity were reduced by 0.20–0.77 m s−1, the Ave. air temperature was reduced by 0.2°C, the Ave. mean radiant temperature was reduced by 0.5–7°C, and the Ave. physiological equivalent temperature was reduced by 0.3–3.2°C. The couple method could support the architectural and landscape design integrating with the GFs in future. GRAPHICAL ABSTRACT


Introduction
Global climate change and the urban heat island (UHI) effect have become much more serious in recent years. Urban sprawl has also led to a rapid decrease in green urban land. In response to these challenges, urban greenery infrastructure (UGI) strategies, which will improve urban climate, urban health, and urban resilience, have been introduced to policymakers and urban planners. The UGI, which includes a multiscale of green lands and vertical greenery spaces such as urban parks, pocket gardens, vertical gardens, green roofs, and green walls, has been rapidly promoted for UHI migration and healthy city policies (Zölch et al. 2016).
In high-density urban areas, vertical greening (VG) strategies have become important methods of addressing the lack of green land. Different strategies of VGs are promoted depending on the urban environment (Pérez et al. 2014). Green roofs (GR) and green walls (GW) are the two main strategies for VG. GRs have been demonstrated to improve energy saving, rainwater preservation, air quality improvement, and UHI mitigation (Liu et al. 2021;Coma et al. 2017). Although GWs have a long history, they have developed fast in recent years due to UHI mitigation, prefabrication technologies, and urban greening support policies (Aflaki et al., 2017). GWs mainly consist of two systems: the living wall system (LWS) and the green façade (GF). Cross- Transmittance: 11.6-29.8%  Guangzhou, China Sum Cfa -Ave. 4.7 °C (East) (Lin, Xiao, and Musso 012006) Guangzhou, China Sum Cfa Transmittance: Approx. 50% (Zhang et al. 2019) Bari, Italy Al Csa -Max.9.0 °C (Vox, Blanco, and Schettini 2018) Air temperature (T a  climate zones and season studies on GWs have demonstrated the potential for energy saving, improving the thermal environment , noise insulation, dust reduction, and biological diversity promotion (Köhler 2008;Sheweka and Mohamed 2012).
Regarding the thermal environment, human thermal comfort (HTC) is an important issue in high-density urban areas. Recent studies of GWs have demonstrated the positive effect of GWs on HTC via the shading effects of solar radiation, the reduction of surface heat reflection, and the evaporation effects of plant and substrate (Table 1). However, despite the increasing application and various design patterns available for GWs, we still lack meaningful design-supported tools and evaluation methods for GWs, particularly for the determination of building performance in terms of architectural and landscape design.
This study is mainly focused on the GFs' effects on HTC. GFs mainly allow the climbing plants to grow vertically via specific lightweight structures and help them form shading screens with their foliage layers (2010Façade Planting [n] 2010]. With the porous character of the foliage layers, the GFs provide shade and wind-induced effects on HTC optimization (Jim 2015b). GFs are widely grown on the external surfaces of building walls, sky gardens, and balconies. GFs can be implemented in various scenarios because of their low cost of planting, maintenance, and simple structure (Perini and Rosasco 2013). Furthermore, different species of climbing plants support high adaptability in different climates (Isnard and Silk 2009;Paul and Yavitt 2011).
In south China's hot and humid climate, the hot months have been increasing in recent years, highlighting the critical need for HTC optimization. GFs can aid in this endeavor, especially for residential buildings, where the local habitants can add plants and simple shading devices on building facades or balconies, which is particularly relevant in the fast-developing, dense urban environment. The number of residential buildings in dense urban environments increased extremely fast due to China's urbanization trend in the last 30 years. In a residential building, the balcony space is the specific space between outdoor and indoor environments and supports various activities, such as clothes drying, growing gardens, resting, and dining. As climbing plants are widely applied in balcony gardens, this study takes balcony GFs as a climate adaptation strategy for hot and humid climates; to achieve this, we set up a research target to evaluate the effects on HTC. Determining the thermal comfort and temperature control of specific GF design patterns Figure 1 would support the future architectural and landscape design of GFs in a wider background.
In this research area, a balcony could be defined as one of the typologies of transitional space as well, which is a space between the outdoors and indoors with different enclosed elements, such as roofs and walls, reducing the effects of solar radiation on human beings (Chun and Tamura 2005;Chun, Kwok, and Tamura 2004). In hot and humid climate areas, the transitional spaces, such as corridors, arcades, and atriums, support various activities of commerce, leisure, community, and children playing (Zhang, Zhang, and Ding 2017); thus, this combination study on balconies with different GFs patterns would deepen the understanding of the optimization strategies for the transitional spaces and extend the possibility for people staying in these spaces despite the high temperature and strong solar radiation in summer.
This research is based on a previous study of GF effects on HTC in the hot and humid climate in southern China (China Meteorological Administration (in Chinese) 2021). The previous study of GFs had conducted field measurements on thermal comfort indices. A computational fluid dynamics (CFD) method was also developed, validated, and then simulated on different patterns of foliage distribution on a standard building with a transitional space (Lin et al. 2019), testing the GFs' regulations on air temperature (T a ) and wind velocity (V a ); however, the evalua- tion of HTC in a transitional space needs to consider the T a and V a and take into account the global irradiation (G a ), human activity, clothing, and the relative humidity (RH). A better evaluation method for HTC in a transitional space in a hot and humid climate is to introduce the index of physiological equivalent temperature (PET) (Höppe 1999), which is the temperature of a reference environment based on a heat balance model that combines various climatic and physiological variables including T a , RH, solar and environmental thermal radiation, V a , clothing, and metabolic rate to provide a synergetic indicator of HTC (Mayer and Höppe 1987). This index is widely used to understand the thermal comfort environment of outdoor spaces and has been tested and validated in local field measurements and questionnaire surveys in the hot and humid climate in southern China (Zhang, Zhang, and Ding 2017;Zhang, Zhang, and Jin 2018). Furthermore, considering a balcony space with GF would be affected by solar radiation, the index of PET is much more suitable than some HTC indices such as PMV-PPD utilizing in an indoor environment. In this study, we aimed to develop a coupled method to calculate the PET based on CFD simulations and evaluate the HTC level on balconies in summer.

Methodology
In the previous study (Lin et al. 2019), the GFs' effects on HTC had been discussed via field measurements and CFD simulations. The field measurements were taken in a corridor of an apartment in summer in Guangzhou, China. In the CFD simulations, the indices of T a and V a were validated and simulated with different façade layouts, preliminarily discussing the shading and wind-induced effect of the GFs in the transitional space. As introduced above, a much more comprehensive evaluation of HTC should include various indices; therefore, the method of the previous study requires improvement.
As the previous study revealed, T g , T r , and PET indices showed high correlations in the field measurements, which implied that the solar radiation provided a strong effect on HTC and the foliage layer of the GFs provided an effective shading effect improving the HTC (Lin et al. 2019); thus, we conducted a study on the effects of GFs on solar radiation. First, the shading effect was estimated with the coverage indices of the foliage layer of the GFs, which are the key parameters for the solar radiation simulation. Second, the coverage indices were set in the previous standard CFD model to test the shading and evaporation effects separately. At the same time, the models of the various GF layout typologies on a standard balcony model were constructed based on cases studies. Fourth, based on the standard parameters of the GFs' foliage layer and the standard CFD models, a couple simulation method was introduced, tested, and provided a deeper study into HTC.
The methods implemented in this study consist of 4 main steps: 1. The field measurements of GFs' foliage parameters; 2. standard model definitions for GFs simulation; 3. the setting of CFD simulations; 4. the couple simulation with solar radiation and calculation for PET.  This study was performed in Guangzhou, China. Guangzhou is located in southern China (23°08ʹN, 113°16ʹE); in summer, the average high temperature is 32.3 °C, and the average relative humidity is 83% (China Meteorological Administration (in Chinese) 2021). Guangzhou has a humid subtropical climate according to the Köppen climate classification (KCC) (Köppen 2011). The experiment was conducted in June. 2019 and the simulations was set on 22 nd June, which is a typical summer day in hot-humid climate.

Field measurements of GFs' foliage parameters
The shading effect of the GFs was estimated in two aspects. First, the solar radiation was measured behind and outside the foliage layer to compare the reduction ratio. Second, the coverage ratio of the foliage indices were measured to obtain the foliage preference.

Measurements of solar radiation
Field measurements of solar radiation were taken over ten days in June and an index of G a was recorded by the HD 31 with a LP471 SILICON-PYRA probe (Figure 2 (a) and Table 2). The measurement period had covered the different weather characters of cloudy and sunny days, then, we stopped the measurement before arainy day to protect the experiment probes. As the foliage layer of a GF is distributed vertically, the measurements of the GFs were conducted horizontally behind the GFs with a distance of 0.3 m at the height of 0.6 m, 0.9 m, 1.2 m, 1.5 m, 1.8 m, and 2.1 m, recording the foliage character in the height range of a human being. The capacity of the G a adjustment of the foliage layer was analyzed and integrated into the simulation model.

Measurements of the coverage ratio of the GFs' foliage layer
In the literature, the measurements of the foliage distribution could be described as foliage coverage ratio (FCR) or leaf area index (LAI) (Jonckheere et al. 2004;Hunt et al. 2013;Azhar Mat Sulaiman et al. 012025;Casadesús and Villegas 2014;Pérez et al. 2017). The former method is based on applications in agriculture and remote-sensing areas, which compares the color pixels of the plants to the surrounding environment. The color range of the foliage was built up under the empirical measurements (Hunt et al. 2013;Azhar Mat Sulaiman et al. 012025;Casadesús and Villegas 2014); however, the foliage color may be different in different species and seasons, which may also increase the probability of errors in the measurements. Both a photo background with high color contrast and a color range recording in situ for the specific plants' foliage are needed in field measurements to decrease the likelihood of error. The latter method mainly compares the solar radiation values of the space in an open space to values collected behind the foliage layer to calculate the shading level of a crown of a plant (Jonckheere et al. 2004;Pérez et al. 2017). This method is also widely applied in agriculture and forest management, as well as research; however, most LAI measurements were measured on a horizontal surface, as the crowns of the trees are mainly distributed vertically. The climbing plants of the GFs normally grow on a vertical surface, and as such, measurements should be taken on a vertical surface, which is often not performed in field measurements.
To compare both methods, we investigated both indices simultaneously with the same plants. To reduce the inference of the surrounding environment, the GFs were chosen on a building with a bare wall as a contrasting background. The measurements were taken at both positions of the shading area of the plants and the open space. The photos of the FCR of the GFs were processed with a photo pixel analysis program (Adobe Photoshop). The LAI values of the GFs were measured with the equipment of LI-COR LAI-2200C (Figure 2(b) and Table 2). Finally, the data taken from both methods were compared and analyzed.

Standard model for GFs' foliage parameters test
Two standard models based on the validation model in the previous study (Lin et al. 2019) were introduced to test the GFs' effects of shading ratio (SR) and foliage evaporation (FE) (Figure 3(a)). In the SR model, the foliage of the GFs is divided into a difference of 10%. The simulation tempted to test the different effect of shading ratio on thermal comfort; however, it would take too much time to simulate the changes of ratio. Thus, the simulations took an interval of 10% to save the computation resource. Considering that the width of a climbing plant on a single rope net is normally less than 0.4 m, the thickness of the foliage layer was defined as a plant with 0.3 m and the height of a climbing plant was set as the floor height (3 m). As the distance of the ropes supporting the climbing plants is normally arranged in 0.2-0.3 m, the minimum distance of two single plant foliage units was set as 0.3 m as well. Thus, the SR of the GFs was formed with different patterns of the single foliage layer.
In the FE model, the leaf area density (LAD) of the foliage was estimated from 0.5-4.0 m 2 m −3 , following the results from field measurements and literature surveys (Susorova, Azimi, and Stephens 2014;Lin et al. 2019;Gromke et al. 2015). According to Gromke et al. and Shashua et al. (Gromke et al. 2015;Shashua, Pearlmutter, and Erell 2009), the evaporation power (P c ) was estimated as P c = 250 W m −3 for a plant with LAD = 1.0 m 2 m −3 . As calculated in previous studies (Lin et al. 2019;Gromke et al. 2015), the P c was set from 125-1000 W m −3 in the FE models. The study of foliage evaporation tents was performed to test the evaporation effects on temperature adjustments.

GFs layout typologies models
The standard balcony model (2.5 m × 3 m, L × W) was built based on the previous standard single CFD model ( Figure 3(b)) (Lin et al. 2019). The normal width of a balcony is 1.5 m; however, this study expanded the width to evaluate the difference of the T a and V a behind the plant foliage. The previous study had conducted the standard simulations on different orientations (Lin et al. 2019). For this study, the main purpose was to develop a couple simulation method, thus, all models were simplified and only set to face south.
Seven cases of GF typologies on the standard balcony model and one unshaded model (UN) were built to test different layout strategies of GFs. GF typologies included the front, the side, the back, and the railing around the balcony. The front GF models were divided into three typologies: GFs in front of the balcony (Front-1, F1), GFs in front of the balcony with openings (Front-2, F2), and GFs with full shade (Front-all, FA). The side GF models were divided into two typologies: double side GFs (DS) and single side GFs (SS). Furthermore, the GF railing model (GR) and the GW model at the back side of the balconies were also defined for comparison. The models are listed in Table 3. The GF thicknesses were set at 0.3 m, using the field measurements as a guide.

The standard CFD model setting and the grid sensitivity test
The standard process of the CFD simulations had been validated on T a and V a in the previous study (Lin et al. 2019). The software of ANSYS Fluent (Version 16.0, 2014 ANSYS, Inc.) was employed to solve the aerodynamics problems (Gromke et al. 2015;Shashua, Pearlmutter, and Erell 2009;Shih et al. 1995;Fidaros et al. 2010;Richards and Hoxey 1993;López et al. 2012;Mercer 2009;Boulard and Wang 2002;Liu, Chen, and Black 1996;Sanz 2003). The solving model was the threedimensional steady-state Reynolds-averaged Navier-Stokes (RANS) equation with the Boussinesq approximation for the thermal environment and the realizable k-ϵ turbulence model (Shih et al. 1995). The discrete ordinates (DO) model was employed in the radiation model to solve the solar ray-tracing calculations (Fidaros et al. 2010). The second-order scheme was applied in the mean flow, turbulence, and energy equations discretization. The SIMPLEC scheme was applied in the coupling calculation with pressure and velocity. The convergence criteria of the continuity, velocity, energy, and DO model were set to 10 −6 . The computational domain of the standard model was 152 m × 108 m × 96 m. The grid sensitivity of the standard model had been tested with the minimum scale of 0.05 m, 0.1 m, 0.2 m, and 0.4 m, and the results confirmed that the differences of T a and V a of different grid scales were lower than 3%; therefore, in this study, the minimum hex grid size of 0.1 m on the central study building and the maximum size of 1.0 m was applied to the computational domain boundary. The total cell number of the standard model was 3,404,835 (Figures 4 and 5).

The boundary settings of the standard CFD model
In the boundary settings of the model, a vertical profile was defined at the wind-inlet boundary. The V a , k, and ϵ of the neutrally stratified atmospheric boundary layer was set according to (Richards and Hoxey 1993). .
where C u = 0.09; z 0 = 0.5 m; V a was set as 1.0 m s −1 at a height of 10 m in the standard model. Using the previous studies as a guide, the properties of the model boundaries of symmetry, outlet airflow, building walls, and the ground were set (Table 4) (Lin et al. 2019;Gromke et al. 2015;Richards and Hoxey 1993;Boulard and Wang 2002;Liu, Chen, and Black 1996;Sanz 2003). The setting of the building walls include thermal insulation and surface material layer is about 240-270 mm in hot-humid climate and about 300-350 mm in temperate climate, thus, it was simplified to defined as 300 mm in the simulation. Considering the simulation comparison was mainly focused on the different green façade models, thus, the multi-layer of the ground material was simplified to a same setting (Gromke et al. 2015).
The GFs were defined as porous zones in the CFD simulation (López et al. 2012;Mercer 2009), whose transpiration cooling power (P c ) was estimated to be 500 W m −3 in previous studies (Lin et al. 2019;Gromke et al. 2015;Boulard and Wang 2002). The wind-induced effect of the GFs was defined with the momentum, k, and ϵ equations terms (Lin et al. 2019;Gromke et al. 2015;Boulard and Wang 2002). .

Couple simulations
To evaluate the human thermal comfort indices, such as PET, more indices were input into the simulation, including the RH, the mean radiant temperature (T r ), and body characteristics of a human being, which were simulated in a series of algorithms using Ladybug tools   The couple simulations with CFD and LB were based on the same building models and the same coordinate system, which supported the simulation data interaction in GH. In the Ansys software package, a post-processor software named CFD-post (in Ansys 16.0) supported the data of T a and V a , and exported the coordinate system data. The central parts of the simulation zone in CFD, including the building surface and the volumes of the balconies and GFs, were meshed with isometric grids with the same grid length of 0.1 m, which defined as the same grid size in LB and HB launching EP. Consequently, the data of T a and V a were read in LB, and the coordinate system was rebuilt in RN ( Figure 6).
For the simulation models, the balcony areas were defined as HBZones in LB and HB and exported into an OpenStudio® (OS, Version 3.1.0) (SD 2021) file for simulation using the software EnergyPlus TM (EP, Version 8.8, Alliance for Sustainable Energy, LLC) (SE 2021). The GFs were defined as shading components of the HBZones with the transmittance being obtained from the field measurements. Consequently, an open-source weather date file (*.epw file) of Guangzhou, China, was collected from EP (SF 2021) and imported to create a specific site climate for the simulation.
We obtained the standard building model geometric parameters from RN, the indices of the building and environment surface temperatures (T s ), and sky view factor (SVF), and used these values for the simulation, using LB and HB, which launched EP on the environmental parameters. The computation of T s was based on the baffle heat balance assumption that the baffle is sufficiently thin and has high conductivity so that it can be modeled using a single temperature (EnergyPlus™ Version 9.4.0 Documentation 2021). The T s is determined by formulating a heat balance on a control volume that enclosed the baffle surface. According to (EnergyPlus™ Version 9.4.0 Documentation 2021), the heat balance on the baffle surface's control volume is: The T s is: In the above equation, DifShadgRatioIsoSky obs;i = (obstructed sky irradiance on obstruction)/(unobstructed sky irradiance on obstruction). SVF obs;i DifShadgRatioIsoSky obs;i is the sky irradiance at the hit point divided by the horizontal sky irradiance when considering the shadow of the sky diffuse radiation on the obstruction caused by other reflections. The irradiance at one hit point by the reflection of sky radiation from the ground or from other obstructions is not considered.
Basing on the calculation above, both indices of T s and SVF were imported to calculate the T r with the HB algorithm. According to the ISO 7726 (International Organization for Standardization. Available online 2021), the T r could be calculated from the surrounding surface weighted by the surface emissivity: Coupling with the data files of T r , V a , T a , and a group of preset values including RH (70%), obtained from the average value in the field measurements and the parameters of a human being (Table 5), the PET was calculated via the algorithm "ThermalComfortIndices" in LB basing on the PET calculation method according to the Munich Energy-balance Model for Individuals (MEMI) (Höppe 1999;Ng and Ren 2015).

Visualization of the simulations results
Visualizations of the CFD and coupling simulation results on the building model yield a clearer picture of the distribution of different thermal environmental indices. In this study, the CFD results (V a and T a ) were visualized in CFD-post, and the coupling simulations results (T r and PET) were visualized in HB and RN ( Figure 6). The simulation results could be applied to different building floorplans and sections to obtain a much more comprehensive understanding of the effects of the GFs.

Solar radiation results
Results of the measurements between the shaded and unshaded area for ten days in June 2019 revealed that the average reduction of the G a was 73.29% (Table 6). The field measurements conducted over 10 days sufficiently captured different weather characteristics. The G a of the unshaded area varied from 53.27 W m −2 to 403.00 W m −2 on different days due to weather changes; the shaded area stayed below 150 W m −2 in the measurement period Figure 7. These results show that the foliage layer presented a significant and stable shading effect and showed significant potential to reduce building energy consumption. Compared to the previous studies, the result of the G a reduction in this study was higher than Jim's measurement of 67.9% (Jim 2015b)

Foliage parameters results
The measurements of the five GF units are shown in Figure 8. The LAI was 1.56-3.61 and the FCR was 49.74-92.88%. Correlation analysis was performed on two group data and showed that two parameters had a high correlation (R 2 > 0.8). These results imply that both methods could determine the coverage capacity of a GF's foliage layer.
Compared to the results from previous studies, Lee & Jim estimated that the LAI of the scatter foliage of a climbing plant is 0.24 (Köppen 2011); Susorova et al. estimated that the LAI of much denser foliage is 2.0 (Susorova, Azimi, and Stephens 2014). In our study, the foliage thickness of the climbing plants ranged between 0.4 and 0.5 m, and had a significantly increased LAI value. Sulaiman et al. had investigated three climbing plants to determine the FCR value via measurements and found that the FCR values were in the range of 69.4% -86.5% (Sulaiman, Jamil, and Zain 2013;Azhar Mat Sulaiman et al. 012025); however, few studies have compared both parameters for the same plants. In our study, we provide a new perspective to evaluate the coverage ratio of a GF.
Regarding the GF simulations, the FCR value was defined as the transmittance of a shading component in LB and HB, and the LAI value was used to calculate the GF's foliage density with the foliage layer thickness computing the C d and the P c in CFD simulations. Both parameters were shifted to the definition of the GF model in the following simulations. Table 5. Preset values of a human being used in our study (Niu et al. 2015;Jin, Zhang, and Zhang 2017

Simulations of SR Model
In the SR model, SR was defined with ten levels to shade the corridor. The results revealed that T a and V a decreased when the SR increased. The Ave. T a and V a of the corridor were reduced by 0.05-0.2°C and 0.15-0.6 m s −1 at different SR levels (Figure 9(a)). A visualization of the T a and V a distribution at the section of the standard model was added. The data were taken as a series of points in a line along the depth of the balconies at the height of 1.5 m. The results in the corridor section provide a much clearer understanding of the adjusting effects on T a and V a, showing that both indices were reduced more when there was a smaller distance to the GFs (Figure 9 (b1 -b2)). Though the changes of the T a were simulated at a low level, the simulations revealed that the shading and windinduced effects of the GFs strengthened when the SR increased.

Simulations of FE Models
The FE models simulated the P c value from 125 to 1000 W m −3 for the LAD value from 0.5 to 4.0 m 2 m −3 (Gromke et al. 2015) a model without FE was added to compare. The results reveal that the Ave. T a of the  corridor decreased with the increase of the P c . The Max. T a was reduced by 0.75°C when the P c was set to 1000 W m −3 in the experiments (Figure 9(c)). In the correlation analysis, the P c and T a presented a negative correlation (R 2 > 0.9) in the experiment, implying that the increase in the P c , as well as the density of the foliage layer, could improve thermal comfort in the corridor.

Results of CFD simulations on the balcony models
The balcony models were first simulated in Ansys Fluent 16.0 to obtain the T a and V a results. The V a results showed that the model DS, SS, GR, and FA provided much more significant effects on the wind environment and the Ave. V a reduction was 0.77 m s −1 , 0.50 m s −1 , 0.29 m s −1 , and 0.20 m s −1 (Table 7).
Compared to the UN model, the DS model showed that Ave. T a was reduced to 0.2°C, whereas the T a reduction of the other models was lower than 0.1°C (Table 7). A section was added to analyze the V a and T a distribution inside and outside the balconies' area. The T a data showed little difference between the cases with a temperature of 0.1°C. Wind static pressure (P r ) showed a much more significant pressure drop in the balconies' areas, especially in FA, F2, F1, and GR. Thus, the V a of FA, F2, F1, and GR also showed much more reduction behind the GFs. The section data allows for a clearer description of the wind-induced effect of the GFs (Figure 10).To provide a deeper analysis of T a and V a differences in the corridor space, the contour maps were created using the CFD-Post software, and the resulting case data were input into the same coordinate system to compare to the UN model. Three sections were defined in CFD-Post to demonstrate the horizontal and vertical distributions of the data of T a and V a .
In the front-vertical section (Figures 11 and 12), the models of DS, SS, GR, and FA showed clearer effects on V a , especially on the bottom and the edge (the right side of the model) areas of the model because of the air turbulence variation at the building corner. T a was reduced slightly in the models of FA, DS, and GR in the  middle of the balconies that implied the GFs may improve the thermal comfort better in these cases. T a results also showed that the top floor exposing to the solar radiation may be overheated and need much more shading or vegetation devices. Comparably, the GW model presented the lowest effect on both indices.  In the cross-vertical section (Figures 13 and 14), the contour maps showed the GFs' effects on the thermal environment in the middle of the balconies. The changes in T a mainly occurred in the first to the third floor in models of DS, FA, and GR, respectively, in which V a mainly appeared in the edge areas in most models except GW. We noticed that the GF reduced the T a at 0.5°C at a distance of 1 m on the ground in front of the balconies, revealing GFs' optimization effect on the thermal environment of the surrounding area. On the other hand, the reduction of V a around the GFs implied that the building design may need to concern the negative effects of GFs on natural cross ventilation.
In the horizontal section on the third floor (h = 7.5 m, Figures 15 and 16), the contour maps also showed the GF shading and wind-induced effects on the T a and V a , especially in DS, GR, and FA models. The GW model presented the lowest effect on T a and V a at a distance of 0.5 m. We also found that the T a around the building was reduced due to the air turbulence  variation. It also implied that the potential of GFs' effect on thermal comfort may not only existed on the balconies but also in the surrounding environment.

Results of T r on the balcony models
Once we performed the CFD simulations for the balcony models, we performed simulations in the LB and HB with the T a and V a data and a coordinate system.
In the T r comparison, the results revealed that the reduction of T r in FA, F1, F2, DS, GR, SS, and GW models was 7°C, 5.5°C, 3.7°C, 3.5°C, 2.1°C, 2.0°C, and 0.5°C, respectively (Table 8). These results imply that the GFs, as the shading component on the front side, provided a suitable reduction of solar radiation. Furthermore, when comparing the T r results of DS and SS models, the results show that increasing the foliage layer area would positively affect the reduction  of solar radiation. The GW model results also show a reduced Ave. T r by 0.5°C and GW also provided regulation of indirect radiation from vertical surfaces.

Results of the coupling simulations on PET
The results of PET were similar to T r and the reduction in PET in the FA, F1, F2, DS, GR, SS, and GW models was 3.2°C, 2.7°C, 1.9°C, 1.4°C, 0.8°C, 0.6°C, and 0.3°C, respectively (Table 8). In the PET category, thermal comfort was considered to be between the "moderate heat pressure" and "strong heat pressure" level for all models. The PET results at the height of 1.5 m revealed that the GFs of the FA, F1, and F2 models showed much more significant effects on thermal comfort (Figure 17). The adjustment of PET in other models appeared to limit the areas with a 0.5 m distance, which implied that the layout of the GFs would influence the optimization of the thermal comfort as well.   The PET and T r results revealed that the HCT in the balcony models was mainly affected by solar radiation, implying that the shading effect would be important in GF applications.

Measurement and prediction indices for SR of GFs
SR was a key factor for the shading effect of GFs on thermal environmental adjustment; however, it is not easy to measure or predict the SR of a plant that grows vertically. Several factors limit the reliability of field measurements. Firstly, GFs include different species of climbing plants with various growing mechanisms, such as tendrils twining, stems and leaves twining, thorns hooking, and roots adsorbing (Jim 2015b;Herron et al. 2007;Burnham 2021). For the climbing plants adsorbing on the wall surface, there is no space between a foliage layer and a wall to install the measurement equipment. Secondly, some field investigations on GWs have shown that GWs installed on the exterior space of high-rise buildings are difficult to maintain due to problems of accessibility, structural stability, and material durability (Chew, Conejos, and Azril 2019). The long-term monitoring equipment of SR integrating with the façade system is often lacking. Thirdly, there is a lack of synergistic controls on thermal environmental adjustment methods with the variation of GF foliage properties and GF monitoring is normally absent.
In this study, we compared two methods that determine the foliage parameter measurements of GFs to provide a simpler method for evaluating the foliage layer using digital photos, which could be collected with simple cameras. We hope that this method supports a better understanding and maintenance of GFs. The high correlation of the two methods validated the method of photo-pixel post-processing using solar radiation measurements. This simpler image collection method could be used for daily maintenance in many different GFs, even using a user's mobile phone. To measure a GF on a high-rise building, a façadecleaning machine or an uncrewed aerial vehicle could support the collection of images. Future research is still required to test foliage images and color pixel recognition with different plant species, in different seasons, and to distinguish between leaves and stems.

Regression analysis of T a and V a in CFD simulations
The CFD simulation in this study was based on the coupling computation with the air and thermal transmission. The GFs were simulated with a definition of a porous zone, which accounted for the wind-induced and evaporation cooling effects. T a and V a are two main indices useful for discussing the results; a regression analysis was also performed. The analysis revealed that the two indices presented a logarithmic relationship, and the R 2 was higher than 0.7 in most models ( Figure 18). With an increase in V a , the T a reduced due to the cooling effect of the GFs; however, in the FA and DS models, the correlation with the two indices was not obvious. Considering the wind-induced effect and the higher SR of these two models, the V a in the above two models had been reduced to a low level such that the distribution of T a in different ranges became dispersive. Figure 19. Regression analysis with T r and PET in the coupled simulations. In this study, we mainly focused on the balcony areas, which were defined as one typology of transitional spaces; further studies should include other urban environments to investigate the previously mentioned indices.

Correlation and regression analysis of thermal indices with PET
In this study, we conducted a coupled simulation with CFD and LB with HB launching EP. The relationship between the indices is required for the analysis of the different effects on HTC.
The correlation analysis of T r and PET revealed that the correlation index (CI) of F1, F2, FA, and DS models was higher than 0.6, and for T a and PET, it was higher than 0.7 in the Un and GW models. The CI of V a and PET was calculated in the minus zone and revealed that V a harmed thermal comfort regulation; however, only the absolute value of CI of the Un, GR, and GW models was higher than 0.7 (Table 9). In the regression analysis, the R 2 of T r and PET was higher than 0.6 in F1 and DS models, for T a and PET, it was higher than 0.5 in the Un and GW models, and for V a and PET, it was higher than 0.7 in Un and GW models (Figures 19, 20, and 21).
Comparing the CI and R 2 of different thermal indices with PET, the T r yielded a much more obvious effect on thermal comfort, which implies that the shading effect contributes a stronger influence on thermal comfort regulation, which was discussed in previous field measurements in different climates Köppen 2011;Lin, Xiao, and Musso 012006;Lin et al. 2019;Jim 2015b). One could notice that different layouts of GFs also showed different impacts on thermal comfort regulation. The front-side layout (F1, F2, and FA) yielded better effects on the shading. Both regression analyses revealed that the GFs affected T a and V a in CFD simulations but contributed little effects on PET in the coupled simulation.

Regression analysis of foliage layer area with thermal indices
Considering that the foliage areas of the different GF models are different, an additional regression analysis was performed on foliage areas with T a , V a , T r , and PET. The results revealed that the R 2 was between 0.3-0.5 (Figure 22), implying that the effects of the foliage area were weak. The layout and the shading method of the GFs may present a much more important effect on thermal comfort, as noted in the discussion above.

Couple method for thermal comfort simulations
More and more couple simulations are conducted with CFD and various building energy simulation (BES) methods (Pandey, Banerjee, and Sharma 2021) for building environment and energy research. Internal and external approaches are the two main methods of couple simulation. The internal method requires that the simulation be conducted in the same software and model; however, the external method usually uses different CFD and BES software but shares the same model and boundary information, such as the climate, the building envelope, and the building scale and coordination system (Pandey, Banerjee, and Sharma 2021). Most studies with couple methods are focused on a better prediction of building energy considering the outdoor and indoor air and heat transmission (Tian et al. 2018). Most studies investigated the thermal environments of indoor environments with various methods to transfer boundary conditions from EP to CFD tools (Tian et al. 2018;Shan et al. 2020).
In this study, we followed the external couple method and introduced the open-source LB and HB algorithm to transfer the data from CFD simulation, launch the EP for T r simulation, and calculate the PET. We must use the same model, the same computation grid, the same coordination system, and the same boundary conditions of two series of simulations for data transportation and calculation for the couple method to work properly.
Some aspects may be improved in future studies following the couple method development. This method assumed that the thermal environment was steady and each simulation was only calculated at a specific time; in future investigations, a continuous simulation method would be useful (Pandey, Banerjee, and Sharma 2021). The simulation data transfer and the models and boundary condition settings were manually conducted in this study and could be improved via automation. The inclusion of building energy assessment could also be beneficial for the couple simulation.

Limitations of GF model settings
Our model only used a standard balcony model with one orientation and one specific time to compare different shading strategies with the method of couple simulation, which could be studied in the border area of models. In the building scale, different aspects such as indoor environment, outdoor environment, and urban canyon may be considered in further studies to evaluate the effects of GFs. In the boundary condition settings, the building orientation and weather conditions, such as outdoor T a , V a , RH, and solar radiation, may be defined and tested following the local climate data in applications. In the definition of GF, this study only tested a standard setting of the foliage layer, however, the foliage of different plants and in different seasons would be different and need much more specific definition in research and applications. Furthermore, more time steps for the simulation may be useful for continuous climate data but require more computation resources, such as more CPU cores and increased time to process. In this study, we constructed a working process that could support future studies in the future.

Conclusions
This study extends the previous study on thermal comfort assessment in a transitional space in a hot and humid climate area. The main purpose of this study was to develop a couple method on HTC simulations that involved GF SR evaluation, related parameter testing, and a working process construction. The main findings and results are as follows: • In the GFs' shading effects evaluation, Ave. G a was reduced by 73.29% in field measure ments. • In the field measurements of the GFs' foliage layer, the LAI range was 1.56 -3.61 and the FCR range was 49.74 -92.88%. Both indices presented a high correlation (R 2 > 0.8). • Two CFD simulations on a standard model were conducted on the SR and FE simulations. In the SR model, the Ave. T a and V a of the  F air View factor of wall surface to air surface temperature q 00 LWR;cav Net longwave (thermal) radiation flux exchange with the outside face of the underlying surface(s) (W m −2 ), "LWR, cav" denotes to the net longwave radiation and the outside face F gnd View factor of wall surface to ground surface temperature q 00 LWR;Env Net longwave (thermal) radiation flux exchange with the air and surroundings (W m −2 ), "LWR, Env" denotes to the net longwave radiation and the air and surroudings F i Angle factor between person and surface q 00 sky Sky surface longwave radiation flux (W m −2 ) FCR Foliage coverage ratio (%) SVF Sky view factor G a Global irradiation (W m −2 ) SVF obs;i Unobstructed sky view factor of the obstruction ¼ 1 þ cos tilt obs ð Þ=2, "obs, i" denotes the obstruction hit by ray i h co Convection coefficient for the outdoor environment (W m −2 K −1 ) RF Reflection factor h c;cav Convection coefficient for the surfaces facing the plenum (W m −2 K −1 ), "c, cav" denotes the convection coefficient RH Relative humidity (%) h r;atm Linearized radiation coefficient for the surrounding atmosphere (W m −2 K −1 ), "r, atm" denotes the radiation for the surrounding atmosphere

RS num
The receiving surface number h r;gnd Linearized radiation coefficient for the ground (W m −2 K −1 ), "r, gnd" denotes the radiation for the ground T a Dry bulb air temperature (°C) h r;sky Linearized radiation coefficient for the sky (W m −2 K −1 ), "r, sky" denotes the radiation for the sky T a;cav Dry bulb temperature for air in the cavity (°C), "a, cav" denotes the air and the cavity Hit obs;i Number of the ray(s) an obstruction, "obs, i" denotes the obstruction hit by ray i corridor were reduced by 0.05 -0.2°C and 0.15 -0.6 m s −1 , respectively, in the SR range of 10 -100%. In the FE model, the Ave. T a of the corridor decreased with the increase of the P c and the Max. T a was reduced by 0.75°C when P c was 1000 W m −3 . • In the CFD simulations on the balcony model, the Ave. V a was reduced by 0.77 m s −1 , 0.50 m s −1 , 0.29 m s −1 , 0.20 m s −1 in the DS, SS, GR, and FA models, respectively, and the Ave. T a was reduced by 0.2°C in the GF models. • In the T r simulations, the reduction of T r in FA, F1, F2, DS, GR, SS, and GW models was 7°C, 5.5°C, 3.7°C, 3.5°C, 2.1°C, 2.0°C, and 0.5°C, respectively. • In the PET calculations, the reduction of PET in FA, F1, F2, DS, GR, SS, and GW models was 3.2°C, 2.7°C, 1.9°C, 1.4°C, 0.8°C, 0.6°C, and 0.3°C, respectively.
In conclusion, the shading effect should be recognized as a key effect on thermal comfort in a transitional space, and the distribution of GF to provide shading on the front side will improve thermal comfort. In this study, the method of couple simulation with CFD and thermal comfort calculation via LB and HB was developed and could support the architectural and landscape design integrating with the GFs.

Disclosure statement
No potential conflict of interest was reported by the author(s).