New Surrogate Model for Wind Pressure Coefficients in a Schematic Urban Environment with a Regular Pattern

Natural ventilation and the use of fans are recognized as sustainable design strategies to reduce energy use while reaching thermal comfort. A big challenge for designers is to predict ventilation rates of buildings in dense urban areas. One significant factor for calculating the ventilation rate is the wind pressure coefficient (Cp). Cp values can be obtained at a high cost, via real measurements, wind tunnel experiments, or high computational effort via computational fluid dynamic (CFD) simulation. A fast surrogate model to predict Cp for a schematic urban environment is required for the integration in building performance simulations. There are well-known surrogate models for Cp. The average surface pressure coefficient model integrated in EnergyPlus considers only a box-shaped building, without surrounding buildings. CpCalc, a surrogate model for Cp, considers only one height of neighbouring buildings. The Toegepast Natuurwetenschappelijk Onderzoek (TNO) Cp Generator model was available via web interface, and could include several box-shaped buildings in the surrounding area. These models are complex for fast integration in a natural ventilation simulation. For optimization processes, with thousands of simulation runs, speed is even more essential. Our study proposes a new surrogate model for Cp estimation based on data obtained from the TNO CP Generator model. The new model considers the effect of different neighbouring buildings in a simplified urban configuration, with an orthogonal street pattern, box-shaped buildings, and repetitive dimensions. The developed surrogate model is fast, and can easily be integrated in a dynamic energy simulation tool like EnergyPlus for optimization of natural ventilation in the urban areas.


Introduction
Building energy consumption is responsible for about 40% of the global energy demand [1][2][3].Due to increasing comfort requirements, global warming, and the urban heat island effect, the demand for air conditioning by individual households has risen [4][5][6].To reduce energy consumption, "passive architecture" or "nearly zero energy buildings", together with sustainable occupant behavior, have become essential.
Natural ventilation is crucial for passive design in warm and hot humid climates, and also important in other conditions [7].Studies in the context of Malaysia and Singapore have proven that thermal comfort can be reached for a large period of the year only with ventilation [8].Airflow through building openings is a critical factor that influences heat and moisture exchange between thermal zones and the outdoor environment [9][10][11][12].The wind pressure is the driving force for airflow and should be predicted carefully [13].In the next section, an airflow network model and existing surrogate models for wind pressure coefficient (Cp) values are reviewed, in order to develop an improved model.The aim is to have a surrogate model that can easily be integrated into dynamic energy simulations in schematic urban environments.

Multi-Zone Airflow Network Model
Airflow networks have been applied in almost all building performance simulations.The airflow network model, embedded in the standard version of building performance simulation as EnergyPlus, provides the ability to simulate airflow, which is driven by wind pressure, through multiple zones.The air flow rate highly influences the predicted energy requirement [14].In order to obtain the air flow rates, the building performance simulation takes the wind's velocity and direction from the hourly weather data, and predicts the wind pressure coefficients at different external nodes (Figure 1).
Atmosphere 2018, 9, x FOR PEER REVIEW 2 of 16 thermal zones and the outdoor environment [9][10][11][12].The wind pressure is the driving force for airflow and should be predicted carefully [13].In the next section, an airflow network model and existing surrogate models for wind pressure coefficient (Cp) values are reviewed, in order to develop an improved model.The aim is to have a surrogate model that can easily be integrated into dynamic energy simulations in schematic urban environments.

Multi-Zone Airflow Network Model
Airflow networks have been applied in almost all building performance simulations.The airflow network model, embedded in the standard version of building performance simulation as EnergyPlus, provides the ability to simulate airflow, which is driven by wind pressure, through multiple zones.The air flow rate highly influences the predicted energy requirement [14].In order to obtain the air flow rates, the building performance simulation takes the wind's velocity and direction from the hourly weather data, and predicts the wind pressure coefficients at different external nodes (Figure 1).

Wind Pressure Coefficient Estimation Method
The wind pressure distribution on the envelope of a building is described by Cp, which is defined as the ratio of the dynamic pressure at a point on the surface over the dynamic pressure in the undisturbed flow pattern, measured at a reference height.Cp values, used to simulate multi-zone airflow network models for natural ventilation, can be obtained from many sources.The first source is full-scale measurements when an existing building is being studied.Precise pressures on a particular building in a specific environment can be measured.However, those pressures are only applicable to that specific building-layout and to that unique environment.Thus, they offer less relevant results for new designs.Moreover, these real scale measurements require a long measurement period and generate a high cost.Another source is wind tunnel tests, which can give more relevant results, because changes in building layout and urban layout are made easily.The limitations of wind tunnels include the requirement of special tools and large wind tunnels for investigating urban models.Finally, the computational fluid dynamic (CFD) approach has the same advantages as wind tunnel tests.CFD analysis is done on powerful computers, but is still timeconsuming, requiring validation and user expertise.
The Cp Generator of TNO (Netherlands Organisation for Applied Scientific Research), a metamodel available via the web, was developed in the Netherlands [15].This meta-model is based on finite element calculations, and has been verified through wind tunnel experiments [16] and measured data [17].This approach offers a rather good correlation between measurements and meta-

Wind Pressure Coefficient Estimation Method
The wind pressure distribution on the envelope of a building is described by Cp, which is defined as the ratio of the dynamic pressure at a point on the surface over the dynamic pressure in the undisturbed flow pattern, measured at a reference height.Cp values, used to simulate multi-zone airflow network models for natural ventilation, can be obtained from many sources.The first source is full-scale measurements when an existing building is being studied.Precise pressures on a particular building in a specific environment can be measured.However, those pressures are only applicable to that specific building-layout and to that unique environment.Thus, they offer less relevant results for new designs.Moreover, these real scale measurements require a long measurement period and generate a high cost.Another source is wind tunnel tests, which can give more relevant results, because changes in building layout and urban layout are made easily.The limitations of wind tunnels include the requirement of special tools and large wind tunnels for investigating urban models.Finally, the computational fluid dynamic (CFD) approach has the same advantages as wind tunnel tests.CFD analysis is done on powerful computers, but is still time-consuming, requiring validation and user expertise.
The Cp Generator of TNO (Netherlands Organisation for Applied Scientific Research), a meta-model available via the web, was developed in the Netherlands [15].This meta-model is based on finite element calculations, and has been verified through wind tunnel experiments [16] and measured data [17].This approach offers a rather good correlation between measurements and meta-model predictions.The TNO Cp Generator approach was also used in another study, in order to obtain the Cps for a large urban fragment [18].
The CpCalc model was developed by analysing data from wind tunnel tests using a parametrical approach [19].The CpCalc model considers only one height of neighbouring buildings.The model considers climate, environment, and building parameters.These independent parameters were varied to generate the data that were used to obtain polynomial functions via regression.Beside the climate parameters, the model allows variance in parameters, such as plan area density, relative building height, building layout, and relative position on façades.The model uses many regression functions.For each of them, the appropriate parameters have to be loaded from a data table.Using many data tables makes it complex to integrate the model into a whole-building energy simulation program.
Blocken and his group have compared the Cp values estimated through these different methods: CpCalc, Cp Generator, and Cp integrated in EnergyPlus [19].Cp integrated in EnergyPlus is the average surface wind pressure coefficient based on the method developed by Swami and Chandra in 1988.Results of the three methods are similar [20].
To date, only a few studies consider surrogate models for the wind pressure coefficient of terraced buildings in urban areas, such as CpCalc, TNO Cp Generator, and free-standing average Cp [19,21].A new parametric equation for the wind pressure coefficient for low-rise buildings fits the Tokyo database values with a goodness-of-fit value of R 2 = 0.992 [22].Although existing surrogate models for Cp are good for analysing individual cases without obstructions, or considering average height of neighbouring box buildings, these models are not appropriate to be integrated into a dynamic simulation with parametric models or an optimisation process.Therefore, as a novelty in this research, a new surrogate model for Cp is developed.This model offers time saving and easy incorporation into the optimization processes, whereby a large number of simulations are required and both building geometry and urban layout can be changed.

Data Description and Methods
The method consists of five steps: (1) defining a schematic urban environment; (2) collecting Cp data from Cp Generator of TNO; (3) multiple linear regression; (4) CFD simulation; and (5) sensitivity analysis.Cp data for 36 orientations for each urban layout are collected to obtain regression functions for each orientation (in steps of 10 • ).

Schematic Urban Environment
The simplified urban layout model, as shown in Figure 2, can be considered as a representation for major urban areas with terraced housing around the world.This urban layout is described by the following key parameters: the height of each of the surrounding buildings (H1 to H14), road width in two directions (Rd, Rl), back garden (Bg), the building depth (D), and building length (L).Different urban layouts are generated based on the ranges and step size mentioned in Table 1.These ranges are defined based on case studies [23,24].

Collecting Cp Data from the TNO Cp Generator
The "Latin Hypercube Sampling method" is applied to generate 200 combinations of 19 independent parameters (14 building heights and five urban parameters) of the building and urban layouts.Input data of those combinations for the Cp Generator consists of defining the terrain's roughness for the wind's different flow directions (kept constant as "urban") as well as "obstacles".In the TNO Cp Generator, each obstacle is referred to via a unique name.The following characteristics are associated: locations (x, y) of a corner point, orientation (β), and size of the buildings.Those data have to be transmitted to the programme as a text file with the correct formatting.In a next step, the text files are uploaded.The results of the Cp values are returned, and consist of a data table for each evaluated point on facades and roofs, and of graphical files, as can be seen in Figures 3 and 4. For our meta-model, seven points are considered: red points on the front facade, blue points on the rear façade, and the green point on the roof.Other positions on the facades will be considered in further research.Points on each facade are equally distributed over the height in the middle of the building.The point on the roof is in the middle.1).

Collecting Cp Data from the TNO Cp Generator
The "Latin Hypercube Sampling method" is applied to generate 200 combinations of 19 independent parameters (14 building heights and five urban parameters) of the building and urban layouts.Input data of those combinations for the Cp Generator consists of defining the terrain's roughness for the wind's different flow directions (kept constant as "urban") as well as "obstacles".In the TNO Cp Generator, each obstacle is referred to via a unique name.The following characteristics are associated: locations (x, y) of a corner point, orientation (β), and size of the buildings.Those data have to be transmitted to the programme as a text file with the correct formatting.In a next step, the text files are uploaded.The results of the Cp values are returned, and consist of a data table for each evaluated point on facades and roofs, and of graphical files, as can be seen in Figures 3 and 4.
For our meta-model, seven points are considered: red points on the front facade, blue points on the rear façade, and the green point on the roof.Other positions on the facades will be considered in further research.Points on each facade are equally distributed over the height in the middle of the building.The point on the roof is in the middle.Wind direction (degree)

Surrogate Model: Multiple Linear Regression
The multiple linear regression approach is a well-known technique for modelling the relationship between two or more parameters.It was selected to develop the surrogate model of wind pressure coefficient.The general multiple linear regression model with parameter X is represented by the following equation.
where, Y is the predicted value; a o , a 1 , . . .a n are constant coefficients of the regression function; X 1 , X 2 , . . .X n are values of n parameters.
To develop the regression model for the wind pressure coefficients, 200 generated design cases were randomly separated into two sub-samples: the training samples (90% of the cases) and testing samples (10% of the cases).The training dataset was used for model development.The testing data was used to evaluate the model.Initial trials showed a very low correlation between Cp values and urban parameters, including the orientations.However, Cp values had a very high correlation with the urban parameters when the regression function was developed for each orientation.Therefore, the surrogate model algorithms consisted of 36 functions for 36 orientations (in steps of 10 • ).The subdivision into 36 regression functions for 36 orientations was also based on the concept of the wind shadow model [25,26].The surrogate model for the wind pressure coefficients is described by the following equations.

•
Cps in middle of front facade: • Cps in middle of rear facade: • Cps in middle of the roof: where j is the orientation, j = 1, 2, . . ., 36 (wind orientations); i is the number of neighbouring buildings (i = 1, 2, . . ., 14 buildings in Figure 2); Cp frontj , Cp rearj and Cp froofj are wind pressure coefficients of front, rear and roof surfaces, respectively; H i , L, D, Rw, RL and Bg were defined in the Table 1.

Computational Fluid Dynamic Simulation
In this section, the CFD is used to provide qualitative and quantitative results to understand the effect of the different wind directions and the urban layout on Cp values.The CFD simulations use OpenFoam with Reynolds-Averaged Navier-Stokes (RANS) equations and a simple algorithm.The computational domain is (600 × 600 × 110) m 3 , full scale.All buildings are 9 m high and 24 m for Rl, Rd, and Bg.The depth and length of the terraced building are 24 m and 120 m, respectively.The basic square grids size is 5 m × 5 m × 5 m, but the mesh (grid) is resized near the buildings, in order to increase the quality of the calculation, Figure 5.The dimensions of the computational domain were chosen based on best practice guidelines by Franke et al. [27] and Tominaga et al. [28].The inlet flow with fixed velocity and the outlet with free pressure were applied for this simulation.

Sensitivity Analysis
In this study, the standardized regression coefficient (SRC) was used to identify the impact of each urban parameter on Cp values.SRCi is calculated by dividing the standard deviation of each parameter (σxi) by the standard deviation of the Cp (σCp), and multiplying the regression coefficient (ai) of each parameter in the multiple linear regression model.The building and urban parameters, together with Cp results of 200 urban layouts, were used to obtain SRC.The parameters with a high absolute value of SRC have an important impact on the Cp values.

Sensitivity Analysis
In this study, the standardized regression coefficient (SRC) was used to identify the impact of each urban parameter on Cp values.SRC i is calculated by dividing the standard deviation of each parameter (σ xi ) by the standard deviation of the Cp (σ Cp ), and multiplying the regression coefficient (a i ) of each parameter in the multiple linear regression model.The building and urban parameters, together with Cp results of 200 urban layouts, were used to obtain SRC.The parameters with a high absolute value of SRC have an important impact on the Cp values.

Cp Surrogate Model Based on the TNO Cp Generator
The model was developed by using 90% case data and 10% data for testing.The predicted Cp results of 20 cases from the testing data were good predictions compared with the TNO Cp Generator values, with R 2 = 0.93 as shown in Figure 6.The testing process is repeated with random data sets, and similar R 2 results were always obtained.The few cases with larger errors when predicting Cp are the cases whereby many neighbouring buildings reach the maximum height (36 m).
ameter (σxi) by the standard deviation of the Cp (σCp), and multiplying the regression coefficient of each parameter in the multiple linear regression model.The building and urban parameters, ether with Cp results of 200 urban layouts, were used to obtain SRC.The parameters with a high olute value of SRC have an important impact on the Cp values.

Cp Surrogate Model Based on the TNO Cp Generator
The model was developed by using 90% case data and 10% data for testing.The predicted Cp ults of 20 cases from the testing data were good predictions compared with the TNO Cp Generator ues, with R² = 0.93 as shown in Figure 6.The testing process is repeated with random data sets, similar R² results were always obtained.The few cases with larger errors when predicting Cp are cases whereby many neighbouring buildings reach the maximum height (36 m).The factors of the multiple linear regression are reported in the Appendix 5.The regression functions show a high correlation (Figure 7).For each urban layout, the R 2 value for 36 Cp values is calculated, in order to obtain the correlation between Cp values obtained from the TNO Cp Generator and estimated Cp values.The R 2 values of 200 urban layouts are shown in Figure 8.
The factors of the multiple linear regression are reported in the Appendix.The regression functions show a high correlation (Figure 7).For each urban layout, the R² value for 36 Cp values is calculated, in order to obtain the correlation between Cp values obtained from the TNO Cp Generator and estimated Cp values.The R² values of 200 urban layouts are shown in Figure 8.

Computational Fluid Dynamic Results and Sensitivity Analysis of Cp Values for the Schematic Urban Layout
The neighbouring buildings modify the external wind flow around the evaluated building.The velocity and pressure of the wind is related to the urban parameters and different angles of attack.The CFD result visualized in Figure 9 shows the effect of the urban layout on the wind pressure and wind direction.The factors of the multiple linear regression are reported in the Appendix.The regression functions show a high correlation (Figure 7).For each urban layout, the R² value for 36 Cp values is calculated, in order to obtain the correlation between Cp values obtained from the TNO Cp Generator and estimated Cp values.The R² values of 200 urban layouts are shown in Figure 8.

Computational Fluid Dynamic Results and Sensitivity Analysis of Cp Values for the Schematic Urban Layout
The neighbouring buildings modify the external wind flow around the evaluated building.The velocity and pressure of the wind is related to the urban parameters and different angles of attack.The CFD result visualized in Figure 9 shows the effect of the urban layout on the wind pressure and wind direction.

Computational Fluid Dynamic Results and Sensitivity Analysis of Cp Values for the Schematic Urban Layout
The neighbouring buildings modify the external wind flow around the evaluated building.The velocity and pressure of the wind is related to the urban parameters and different angles of attack.The CFD result visualized in Figure 9 shows the effect of the urban layout on the wind pressure and wind direction.The detailed graphical and numerical results are provided in Figures 10-13.The Cp values of the middle point on the front façade of 200 cases were used to calculate SRC (the right part of the Figures 10-13).The results shows the air pressure in a horizontal plane at half the building height above the ground, and in the Section A-A (Figure 10).On the right, the SRC of Cp results from the TNO model indicate the impact of each parameter on the Cp values for a certain direction.

•
If the wind angle is 0° (Figure 10), the height of building in font (H5) has the strongest effect on the Cp values.The width (D), length of the building (L) and the back garden (Bg) also have a strong impact on the Cp.They impact the turbulence of the wind flow.Buildings H2 and H5 block the wind flow and generate negative Cp on the front façade of the evaluated building.

•
If the wind angle is 30° (Figure 11), building H5 has also a high impact on the Cp value.

•
If the wind angle is 60° (Figure 12), building H10 has a very high effect on the Cp value, because the back garden varies from 8 m to 24 m.

•
When the wind flows parallel to the façade (Figure 13), building H10 and the width of road parallel the terraced row strongly impact on the Cp value.
The four cases illustrate that the dimensions of the neighbouring buildings and the width of the roads will generate various Cp values.This physical effect depends on the interplay of several parameters, such as wind direction, building layout, and obstructions.Figures 10-13 support the understanding for different wind directions, wind pressures, and the wind speed at different locations in 3D.
The Figure 14 shows that the height of the buildings in front (H10) and behind (H5) the evaluated building impact strongly the Cp of the front façade.The width of the road parallels the evaluated building (Rl) also has a strong effect on Cp values.The detailed graphical and numerical results are provided in Figures 10-13.The Cp values of the middle point on the front façade of 200 cases to calculate SRC (the right part of the Figures 10-13).The results shows the air pressure in a horizontal plane at half the building height above the ground, and in the Section A-A (Figure 10).On the right, the SRC of Cp results from the TNO model indicate the impact of each parameter on the Cp values for a certain direction.

•
If the wind angle is 0 • (Figure 10), the height of building in font (H5) has the strongest effect on the Cp values.The width (D), length of the building (L) and the back garden (Bg) also have a strong impact on the Cp.They impact the turbulence of the wind flow.Buildings H2 and H5 block the wind flow and generate negative Cp on the front façade of the evaluated building.

•
If the wind angle is 30 • (Figure 11), building H5 has also a high impact on the Cp value.

•
If the wind angle is 60 • (Figure 12), building H10 has a very high effect on the Cp value, because the back garden varies from 8 m to 24 m.

•
When the wind flows parallel to the façade (Figure 13), building H10 and the width of road parallel the terraced row strongly impact on the Cp value.The four cases illustrate that the dimensions of the neighbouring buildings and the width of the roads will generate various Cp values.This physical effect depends on the interplay of several parameters, such as wind direction, building layout, and obstructions.Figures 10-13 support the understanding for different wind directions, wind pressures, and the wind speed at different locations in 3D.
The Figure 14 shows that the height of the buildings in front (H10) and behind (H5) the evaluated building impact strongly the Cp of the front façade.The width of the road parallels the evaluated building (Rl) also has a strong effect on Cp values.This surrogate model is very significant in a case where a huge number of simulations is required, like in an optimization process.However, the new Cp model considers the middle terraced unit of the terraced row.Estimating Cp of a house located closer to the edges of the terraced row requires further research.A more robust model will have to consider the Cp along the evaluated building, instead of only at the middle points of the building.

Figure 1 .
Figure 1.Plan view of a simple airflow network and computational fluid dynamic (CFD) results show possible airflow pattern.

Figure 1 .
Figure 1.Plan view of a simple airflow network and computational fluid dynamic (CFD) results show possible airflow pattern.

16 Figure 3 .
Figure 3.The geometry of the building and obstacles as submitted to the Toegepast Natuurwetenschappelijk Onderzoek (TNO) Cp Generator TNO.Results are obtained for points 1, 2, and 3 on the front façade; and points 4, 5, and 6 on the rear facade.with heights from the ground 1.5 m, 4.5 m, and 7.5 m, respectively.Point 7 is located in the middle on the roof.

Figure 3 .
Figure 3.The geometry of the building and obstacles as submitted to the Toegepast Natuurwetenschappelijk Onderzoek (TNO) Cp Generator TNO.Results are obtained for points 1, 2, and 3 on the front façade; and points 4, 5, and 6 on the rear facade.with heights from the ground 1.5 m, 4.5 m, and 7.5 m, respectively.Point 7 is located in the middle on the roof.

Figure 3 .
Figure 3.The geometry of the building and obstacles as submitted to the Toegepast Natuurwetenschappelijk Onderzoek (TNO) Cp Generator TNO.Results are obtained for points 1, 2, and 3 on the front façade; and points 4, 5, and 6 on the rear facade.with heights from the ground 1.5 m, 4.5 m, and 7.5 m, respectively.Point 7 is located in the middle on the roof.

Figure 4 .
Figure 4. Wind pressure coefficient values of the middle point in the front façade of the terraced house, two buildings in front of the terraced building largely influence the wind pressure coefficient of direction 150° to 210°.

Figure 4 .
Figure 4. Wind pressure coefficient values of the middle point in the front façade of the terraced house, two buildings in front of the terraced building largely influence the wind pressure coefficient of direction 150 • to 210 • .

Figure 5 .
Figure 5.The computational domain (600 × 600 × 110) m³ full scale with the grid is resized near the buildings.

Figure 5 .
Figure 5.The computational domain (600 × 600 × 110) m 3 full scale with the grid is resized near the buildings.

Figure 6 .
Figure 6.Predicted Cp and Cp obtained by the TNO Cp Generator of 20 urban forms of the testing data.A total of 720 points ( 36 orientations × 20 cases) are represented in the Scatter chart.

Figure 6 .
Figure 6.Predicted Cp and Cp obtained by the TNO Cp Generator of 20 urban forms of the testing data.A total of 720 points ( 36 orientations × 20 cases) are represented in the Scatter chart.

Figure 7 .
Figure 7. Cp value estimation for terraced house patterns by using regression approach, based on Cp results of 200 scenarios from the TNO Cp Generator web base.3D models at the left and right show two examples of randomly-generated urban layouts.The analysed building is always in the centre.

Figure 8 .
Figure 8. Goodness-of-fit value (R²) values of the linear link between the predicted Cp values and Cp values of the TNO Cp Generator of 200 urban patterns; average R² of the front façade = 0.93 and average R² of the rear façade = 0.95.

Figure 7 .
Figure 7. Cp value estimation for terraced house patterns by using regression approach, based on Cp results of 200 scenarios from the TNO Cp Generator web base.3D models at the left and right show two examples of randomly-generated urban layouts.The analysed building is always in the centre.

Figure 7 .
Figure 7. Cp value estimation for terraced house patterns by using regression approach, based on Cp results of 200 scenarios from the TNO Cp Generator web base.3D models at the left and right show two examples of randomly-generated urban layouts.The analysed building is always in the centre.

Figure 8 .
Figure 8. Goodness-of-fit value (R²) values of the linear link between the predicted Cp values and Cp values of the TNO Cp Generator of 200 urban patterns; average R² of the front façade = 0.93 and average R² of the rear façade = 0.95.

Figure 8 .
Figure 8. Goodness-of-fit value (R 2 ) values of the linear link between the predicted Cp values and Cp values of the TNO Cp Generator of 200 urban patterns; average R 2 of the front façade = 0.93 and average R 2 of the rear façade = 0.95.

Figure 9 .
Figure 9.The external air pressure of the case when the wind angle is 0° (perpendicular the front façade).

Figure 9 .
Figure 9.The external air pressure of the case when the wind angle is 0 • (perpendicular the front façade).

Atmosphere 2018, 9 , 16 Figure 10 .
Figure 10.Left: visualization of pressure (colour) and wind direction (arrows); Right: Standardized regression coefficient (SRC) for orientation at 0° for the 19 variables, based on input data and Cp values of 200 cases.

Figure 10 .
Figure 10.(Left) visualization of pressure (colour) and wind direction (arrows); (Right) Standardized regression coefficient (SRC) for orientation at 0 • for the 19 variables, based on input data and Cp values of 200 cases.

Figure 10 .
Figure 10.Left: visualization of pressure (colour) and wind direction (arrows); Right: Standardized regression coefficient (SRC) for orientation at 0° for the 19 variables, based on input data and Cp values of 200 cases.

Figure 11 .
Figure 11.Left: visualization of pressure (colour) and wind direction (arrow); Right: SRC for orientation at 30° for the 19 variables, based on input data and Cp values of 200 cases.

Figure 11 . 16 Figure 12 .
Figure 11.(Left) visualization of pressure (colour) and wind direction (arrow); (Right) SRC for orientation at 30 • for the 19 variables, based on input data and Cp values of 200 cases.Atmosphere 2018, 9, x FOR PEER REVIEW 11 of 16

Figure 13 .
Figure 13.Left: visualization of pressure (colour) and wind direction (arrows); Right: SRC for orientation at 90° for the 19 variables, based on input data and Cp values of 200 cases.

Figure 12 . 16 Figure 12 .
Figure 12. (Left) visualization of pressure (colour) and wind direction (arrows); (Right) SRC for orientation at 60 • for the 19 variables, based on input data and Cp values of 200 cases.

Figure 13 .
Figure 13.Left: visualization of pressure (colour) and wind direction (arrows); Right: SRC for orientation at 90° for the 19 variables, based on input data and Cp values of 200 cases.

Figure 13 .
Figure 13.(Left) visualization of pressure (colour) and wind direction (arrows); (Right) SRC for orientation at 90 • for the 19 variables, based on input data and Cp values of 200 cases.

Figure 13 .
Figure 13.Left: visualization of pressure (colour) and wind direction (arrows); Right: SRC for orientation at 90° for the 19 variables, based on input data and Cp values of 200 cases.

A
computationally intensive step to calculate Cp values of various urban patterns and building geometries was solved by developing a surrogate model.This multiple linear regression model, with different functions for different orientations, is based on other surrogate models that were developed based on measurement in situ, wind tunnel measurement, and CFD simulations.The developed meta-model of wind pressure coefficients allows for the inclusion of the effect of the schematic urban environment.The regression functions of this meta-model show a high correlation with the values of the existing surrogate model.This meta-model is able to predict the Cp values on different surfaces very quickly.

Table 1 .
Numerical variables and their design options for urban layout.