Numerical Simulation and Optimization of Wind Effects of Porous Parapets on Low-Rise Buildings with Flat Roofs

0is paper presents a procedure to optimize the porosity of parapets to improve the aerodynamic behavior of low-rise buildings with flat roofs, by coupling an optimization algorithm and computational fluid dynamics (CFD) simulations. 0e performance of solid parapets to decrease the wind suctions on flat roofs induced by conical vortices was firstly studied, based on four turbulence closure models (standard k-ε, RNG k-ε, SST k-ω, and RSM). 0e simulation results were validated by comparing with the wind tunnel data. Additionally, the porous parapet was treated as a momentum sink in the governing momentum equation, and the RSM turbulence model was employed. As a result, six optimization studies focusing on the highest mean suction minimization that consider parapet height were presented. 0e aim of this paper is to search for the best performing porosity through an automatic CFD-based optimization methodology. At low relative heights (hp/H� 0.01∼0.05, hp is the parapet height, and H is the roof height), the porous parapet with optimal porosity in between 38.2% and 52.3% seems to be more effective than solid parapets in attenuating high corner suctions generated by conical vortices; however, the solid parapet gives the best performance in the reduction of wind suctions when hp/H≥ 0.07.


Introduction
Wind loads on flat or low-sloped roofs of low-rise buildings have been evaluated in many wind tunnel investigations.Of particular concern are the worst mean and peak suctions observed on roof corners near the edge for oblique wind directions, where the existence of conical or delta wing vortices is well established at both model scale and full-scale [1]. is phenomenon becomes increasingly important for the design of low industrial buildings, because of their relatively light dead load and flexibility that making them susceptible to wind loads.Various types of parapets, such as solid parapets, partial or porous parapets, and rooftop spoilers, were widely used in the past to reduce high corner suctions caused by the conical vortices [2][3][4][5][6].
A number of porous parapet configurations have been investigated in wind tunnel and full-scale experiments by previous researches.Surry and Lin [7] found that porous parapets are the most effective parapet configuration, leading to a reduction of the corner suctions of up to about 70%.e porous parapets are expected to (a) disrupt the formation of the corner vortices, (b) absorb some of the flow energy over the roof edges, and (c) generate weakened vortices that formed over the top edge of the parapet by displacing them away from the roof.Banks [8] examined three porous parapet configurations: a mesh of 50% porosity, a porous parapet with circular holes having 50% porosity, and a slatted fence.ey suggested that the effectiveness of the three parapet types is similar, with circular holed parapet providing marginally better results.Pindado and Meseguer [9] conducted tests with six porous parapet types and found that low relative height parapets (h p /H < 0.05) with a medium porosity are more effective than solid parapets to reduce the corner suctions.Kopp et al. [10] also tested a porous perimetric parapet with 50% porosity, which can efficiently mitigate roof suctions in the corner, edge, and interior zones as defined by ASCE 7-16 [11].Lin et al. [12] and Blessing et al. [13] carried out full-scale measurements, utilizing a 6-fan Wall of Wind (WoW) testing apparatus, to investigate the effect of modified roof edge shapes in reducing the peak suctions for the wind direction of 45 °. ey found the cantilever AeroEdge ™ screen with about 50% porosity reduced the peak suction up to 75% in the roof corner compared to the conventional edge shape (i.e., vertical solid parapet).
More recently, several studies by means of computational fluid dynamics (CFD) simulations were carried out to investigate the performance of solid parapets in reducing the roof suctions of low-rise buildings.Poitevin et al. [14] numerically studied surface pressures on open canopy structures with solid parapets.ey found that CFD simulations, by using Reynolds-averaged Navier-Stokes (RANS) equations with proper turbulence models, are in good agreement with the wind tunnel results.Liu et al. [15] applied standard k-ε turbulence model to the wind field of flat roof with solid parapets, and the effect of wind direction was analyzed in detail, incorporated with wind tunnel testing.Li and Liu [16] investigated the unsteady features of conical vortices on a flat roof on the basis of large eddy simulation (LES).e study shows that the addition of solid parapet can mitigate the peak corner suctions, and the dimension of the conical vortices becomes larger with increasing parapet height.Aly and Bresowar [17] carried out an aerodynamic mitigation study of wind-induced lift forces on low-rise buildings using the k-ε and RSM turbulence models, and the mitigation features includes solid parapets (barriers), circular edges, inclined edges, and airfoil edges.
e numerical results revealed that solid parapets can bring a significant reduction in the lift force of about 20% but increase the drag force on the whole structure (building and parapet) by 32.8%.It is noted that none or little attention has been paid on the CFD simulations of wind effects of porous parapets on low-rise buildings.Despite that, porous fences have been widely used for wind protection in many technical applications, ranging from dust erosion control to pedestrian comfort [9].
To the authors' knowledge, experimental methods using wind tunnel testing or CFD techniques provide the basis of traditional "cut and try" approach for the design of parapets located at roof edges.In this approach, several parapets with different porosities or heights are investigated, and the one that yields the best performance in reducing high suctions is identified.With the rapid development of computer hardware performance and CFD techniques in recent years, a number of optimization algorithms have been proposed in the literatures for a variety of applications including CFDbased aerodynamic shape or porosity optimization [18][19][20].One relevant study by San et al. [20] conducted the aerodynamic porosity optimization of wind fences for the shelter effect using the gradient algorithm, and the optimal fence porosities of 10.2% and 22.1% were determined, which provide the best shelter effect in the near and far wake regions, respectively.However, none or little effort has been dedicated to the optimization study of roof-edge parapets, which can help search for an optimal porosity or configuration of parapets for mitigating high corner suctions on building roofs more efficiently.
e main objective of the current paper is to propose a CFD-based optimization methodology to optimize the parapet porosity, for improving the aerodynamic behavior of low-rise buildings with flat roofs.In the first part of this paper, different turbulence closure models (standard k-ε, RNG k-ε, SST k-ω, and RSM) were employed to simulate the effects of solid parapets on the wind pressure distributions on a flat roof and were evaluated using the wind tunnel data of Stathopoulos et al. [5].Additionally, the porous parapet was treated as a momentum sink in the governing momentum equation.Analogical CFD simulation analysis of surface pressures on the ground and a triangular prism behind a porous fence was conducted and compared with the experimental results [21], in order to validate the numerical method used in this study.In the second part, the influence of the height of solid parapets on roof surface pressures has been analyzed numerically, and six parapets with relative heights h p /H from 0.01 to 0.13 were selected for the aerodynamic optimization analysis.e implementation of the optimization procedures using gradient algorithm was then carried out to find the optimal porosity, which aims at improving the corner suction mitigation performance of parapets early in the design.

Numerical Methods
2.1.Governing Equations.Numerical simulations were carried out in this study through the commercial CFD software, ANSYS Fluent 15.0.Airflow is assumed to be steady, incompressible, viscous, and isothermal.e continuity equation and Reynolds-averaged Navier-Stokes (RANS) momentum equation presented in equations ( 1) and (2), respectively, were solved in each of the simulation cases.
where x i is the Cartesian coordinate and u i is the mean velocity component in the ith direction, p denotes the deviation of pressure from its reference value, and ρ and μ are the density and dynamic viscosity of the air, respectively.S i is a source term and assumed to be zero everywhere except inside the porous media [22][23][24], which is treated as a momentum sink as discussed in Section 2.3.−u i ′ u j ′ is the Reynolds stress caused by turbulence [23].
Four turbulence closure models were tested to predict the wind effects of parapets on low-rise buildings: standard k-ε model, Re-Normalization Group (RNG) k-ε model, Shear Stress Transport (SST) k-ω model, and Reynolds-Stress model (RSM).Details of these turbulence models are given in the ANSYS Fluent 15.0 User's Guide.

Numerical Validation by Stathopoulos's Experiment.
Stathopoulos et al. [5], denoted STA hereafter, has experimentally measured the surface pressure coefficients on flat roofs with solid parapets for the wind direction of 45 °. us, in order to validate the possible applicability of the turbulence closure models used later (standard k-ε, RNG k-ε, SST k-ω, and RSM) for solving these cases, the mean pressure study results from STA's experiments are compared.Dimensions of the building model are length (L) × width (W) × height (H) 150 mm × 150 mm × 75 mm.Perimetric parapets with height h p 5 mm are continuous around all four walls of the building, which corresponds to 1 m in full-scale (using a model scale of 1 : 200).STA's experiments were performed in an open country boundary layer (z < 50 m) of that wind tunnel, which has a roughness length, z 0 , of approximately 2 cm at a scale of 1 : 200.
e mean velocity pro le has a power law exponent α 0.15 and is given by where u ref 13 m/s is the wind speed at the reference height (z ref 600 mm). Figure 1 shows the inlet velocity pro le normalized with u ref .Note that we do not have the experimental data about inlet pro les of turbulence kinetic energy (k in ) and its dissipation rate (ε in ) and speci c dissipation rate (ω in ).us, the inlet boundary conditions for k in , ε in , and ω in have been obtained from the assumption of equilibrium boundary layer: where κ is the von Karman constant with the value of 0.4 [25,26]; I u is the free-stream turbulence intensity and is about 2.5% at the building height (H 75 mm) in the simulated boundary layer.No slip condition and standard wall functions have been considered for the ground, building surfaces, and solid parapets.Free slip conditions were assumed for the top and side boundaries.Outlet conditions in terms of zero normal gradients for all variables were adopted, and the outlet ow has been considered fully developed (out ow). Figure 2 presents the schematic illustration of the computational domain, and a multiblock structured grid system (hexahedral grids) has been used for the space discretization.Dimensions of the domain are 42H (length) × 26H (width) × 8H (height), which leads to a blockage ratio around 1.5%.e building was located at 1/4 of the domain length from the inlet boundary.Also, simulations with a domain height of 16H have been conducted and negligible changes in the results (<1% in the maximum of mean pressures on the at roof ) have been found, which reveals that 8H was considered high enough without being a ected by the boundary conditions imposed at z 8H. e domain was split into two subdomains (inner cylindrical region and outer region) as it is shown in the insets in Figure 2. e two subdomains were connected via a nonconformal interface with tetrahedral grids to improve the quality of the mesh in a fully resolved boundary region [27].Nonuniform spacing grid arrangement was employed for the building model, which has 72 × 72 grid points on the roof surface and 48 grid points along the building height.Note that, more dense grid is required near the roof corner to take account of the steep change of turbulent properties due to the existence of parapet.
erefore, a ner grid with Δx/h p 0.0625, Δy/h p 0.0625, and Δz/h p 0.0625 at the roof corner was adopted, and the expansion ratio in the nonuniform grid system along the building and parapet edges did not exceed 1.3. is resolution grid was found to be enough to obtain correct grid-independent results and was also used for fence ow simulations by Santiago et al. [28].e entire computational grids keep the maximum expansion ratio of the grids less than 5, and the estimated values of y + on the roof and parapet are between 30 and 60.
e grid comprised totally about two million individual cells after a few grid independency tests.
e governing equations coupled with the boundary conditions have been solved by using a collocated grid system based on the finite volume method (FVM).Discretized equations are solved by the segregated solution algorithm.e SIMPLEST algorithm has been employed for the pressure-velocity coupling, which is a modification of the SIMPLE (semi-implicit method for pressure-linked equations) algorithm based on the linkage between pressure and velocity [29].e SIMPLEST method has been widely used for simulating wind pressure fields, for example, in applications for irregular-plan shape buildings [30] or singlesided ventilated buildings [31].e solutions were considered converged when the normalized residuals were below 1 × 10 −5 .
Figure 3 compares the mean pressure distributions along the central line of the flat roof obtained from SAT's experiments and present computation results.e data were derived from wind direction of 45 °.Surface pressures have been normalized by the mean dynamic pressure (0.5ρu 2 H ) at the building height and are given by C p � Δp/(0.5ρu 2 H ), where Δp is the mean pressure difference between the surface pressure and a reference pressure.Figure 3 suggests that the performance of any of the four turbulence models against the wind tunnel data is not particularly good.However, the agreement of RSM turbulence model with experimental measurements is acceptable, although C p values in the case of parapet height h p /H � 1/15 are slightly overestimated at x/L > 0.8.Standard k-ε model provides, in general, the highest suction at the windward edge and poor performance of the attached flow simulation in the leeward region.
is turbulence model is generally inaccurate for flows with adverse pressure gradient and underpredicts the size of recirculating zones in some cases [32].RNG k-ε model presents an overprediction of wind suctions in the leeward region for the no parapet case (Figure 3(a)) and near the windward edge for the solid parapet case (Figure 3(b)).Numerical results of SST k-ω model are worse than those of other turbulence models, which underestimates the absolute values of the mean pressure coefficients.As a result, the RSM model was employed to model surface pressure features in the later simulation procedure.Furthermore, it is noted that there are not detailed test measurements close to the windward corner of the roof.us, the turbulence models cannot be verified near the windward corner, and the numerical results obtained there should be observed carefully.

Modelling of Interaction of Wind with a Porous Media.
Most previous computational studies modelled wind fences as a porous jump region (assuming a thin fence), e.g., Guo and Maghirang [22], Song et al. [24], Santiago et al. [28], and Lee and Lim [33].In the present study, porous parapets are also aerodynamically represented in the CFD analyses as porous media.It is assumed that the porous domain of the parapet is homogeneous and isotropic and is treated as a pressure discontinuity surface by using a porous jump model.us, a momentum sink (i.e., pressure loss) is added to the momentum equation ( 2) to model the resistance to the flow following the works of Wilson [34].For steady-state, incompressible flow through a porous media, the pressure loss can be calculated from the relation Δp � 1/2k r ρv 2 , where v is the normal velocity component and k r is the pressure loss coefficient, and an empirical formula k r � 1.04(1 − ϕ 2 )/ϕ 2 provided by Reynolds [35] was employed.Here, ϕ is the porosity defined as the ratio of the open area to the total area of the porous parapet.Note that k r tends to be infinite in the case of solid parapet (ϕ � 0), and as a result, the solid parapet in the simulations was treated as the wall boundary condition.
As far as we are aware, most previous studies focused on the wind effects of solid parapets on roof corners, however, sufficient experimental data are not available for porous ones.It may be the first study to simulate the surface pressures on a flat roof with porous parapets embedded in an atmospheric boundary layer.In order to verify the numerical method used in this study to simulate the porous media flow, the experimental results of Park and Lee [21] were chosen as a benchmark.
ey carried out the point pressure measurements on the triangular prism model and ground surface along the midspan of the porous fence in an atmospheric boundary layer (power law exponent, α � 0.14).
Figure 4 presents the computational domain and coordinate system for the porous fence and triangular prism model used in the simulations.e height (h f ) and length of the numerical fence, were the same as the prism crest height (h � 40 mm) and length (18h), respectively.e prism model, with a span of 95 mm, was fixed at a distance of 60 mm from the fence.Computational domain ranges from x/h � −15 to 45 in the streamwise direction and from z/h � 0 to 15 in the vertical direction with fence located at x/h � 0. e resolution grid (Δx/h f � 0.0625 and Δz/h f � 0.0625) close to the fence, and the types of boundary condition (standard wall functions at ground and prism, outflow at outlet boundary, and free slip conditions at the top and side), are the same used in the parapet cases (see Section 2.2). e porous fence is represented as a momentum sink.e structural grid of 22438 cells was used in the domain discretization.is configuration has been simulated by using the RSM turbulence model mentioned above, with the same simulation set up representing STA's experiment.Inlet profiles of the mean velocity and turbulence intensity were obtained from the experimental data and computed in the same way that for parapet cases.
Figure 5 shows the mean pressure distributions around the fence and prism model obtained from the simulation and wind tunnel experiment.It can be seen that fence porosity has a significant influence on the mean pressure distributions on the ground and prism surface.For porous fences of ϕ � 20% and 65%, the present simulation results are in general in a good agreement with the experimental data. is indicated that the porous jump model and turbulence closure model could better predict the surface pressure and flow field around the porous media.us, the verified numerical method (porous jump model, turbulence model, boundary conditions, etc.,) was applied to predict the wind effects of solid and porous parapets on a flat roof.

Optimization Model.
e essential components of an aerodynamic porosity optimization problem are design 4 Advances in Civil Engineering  Advances in Civil Engineering variables, objective functions, a flow solver, and an optimization algorithm.In the current study, the parapet porosity was chosen as the design variable, which is one of the most effective parameters for attenuating high corner suctions on building roof [9].us, the aerodynamic porosity optimization problem can be described as follows: Min F(ϕ) where F(ϕ) is the objective function, ϕ is the design variable, and l ϕ and u ϕ are, respectively, the lower and upper limits of the design variable.Different goals of the porosity optimization allow the designer to explore more alternatives for the design of roofedge parapets.It is well known that high suctions are the result of conical vortices formed along the roof edge for oblique wind directions (Figure 6(a)), which are the main reason for most of wind-induced damages or wind-borne debris.It has been experimentally demonstrated that the highest mean suction occurs, in general, in the zone under the influence of conical vortices, in particular, beneath the average vortex core [36].In addition, Banks et al. [37] have pointed out that the mean pressure profile is a reasonably accurate reflection of the typical shape of the instantaneous pressure profiles.us, the objective of the optimization in this study is defined to minimize the highest mean suction in the zone affected by conical vortices, and F(ϕ) is given by where C pi,CACV (ϕ) is the mean pressure coefficient of the ith node in the case of porosity ϕ, and the subscript CACV represents the coverage area of conical vortices.CACV can be judged by mean pressure profiles along lines normal to the windward edge, which exhibit bell-shaped suction curves as shown in Figure 6(b).It is generally accepted that the inflection point in mean suction profiles, where the suction gradient reaches its maximal, is considered as the reattachment position (coverage area) of conical vortices.e above objective function will be evaluated through the CFD approach as discussed in Section 2.

Optimization Algorithm.
e aerodynamic porosity optimization used in the current study can be treated as a single-objective optimization problem, and there exists a unique optimal solution.e strategy for solving the porosity optimization problem entails the utilization of the gradient algorithm.One of the key characteristics of this algorithm is that it computes faster than the nongradient based algorithms.
e basic idea of the gradient based methods is that the objective function reaches its minimum or maximum in the gradient direction.e iterative algorithm of the gradient algorithm by means of design variable porosity (ϕ) is given as follows: where i is the iteration step, α i is the search step size derived by the golden section method and d (i) is the negative gradient value and calculated by d (i) � −∇ (i) F(ϕ) � −zF(ϕ)/zϕ.More details on the gradient algorithm can be found by Tsai et al. [38].e convergence criterion of iteration is defined as where τ is the convergence tolerance, and the optimization procedure will terminate as τ ≤ 1 × 10 −4 .Generally, the basic process for an aerodynamic porosity optimization consists of the following steps: (a) assume the initial design variable ϕ (i) and objective function F (i) (ϕ) based on CFD, where i � 0; (b) update the design variable ϕ (i+1) by the gradient algorithm, and compute the new objective function F (i+1) (ϕ); (c) check the convergence criterion, if not satisfy, and (d) repeat the second and third steps until the convergence tolerance is satisfied.A general approach for the aerodynamic porosity optimization utilizing CFD is shown in Figure 7, which was implemented by a FORTRAN program for searching the optimal parapet porosity.

Effect of Parapet Height.
It is clear that parapet height, h p , can significantly influence the roof corner suctions for oblique wind directions.Previous studies have shown that the solid parapet tends to raise the corner vortices above the roof surface, for instance Stathopoulos et al. [5].Data indicates that the addition of low parapets (<0.5 m in height) have been found to increase peak suction on the roof considerably, compared to the roof without a parapet.However, higher parapets (h p ≥ 1.0 m) can lead to a significant reduction in corner suctions.In the present study, CFD simulations of wind pressures on a flat-roof structure in STA's experiment were performed, for solid parapets with relative height h p /H ranging from 0 to 0.2, to verify the above previous findings.
Figure 8 shows the variation of the minimum mean pressure coefficients (i.e., highest mean suction), C pmin , on the roof with the relative parapet height for a 45 °oblique wind direction.Experimental results of Pindado and Meseguer [9] are also included in Figure 8 for comparison, which are obtained from the same roof model used in present simulations.Wind tunnel tests have been carried out in type I atmospheric boundary layer, and the turbulence intensity of the free-stream flow at the roof height is I u � 2.5%.It can be seen from Figure 8 that the present computation, in general, agrees well with the experimental results of Pindado and Meseguer [9].As it was expected, the minimum mean pressure coefficients are significantly affected by the parapet height.It can be observed that the highest mean suction behind solid vertical parapets increases (C pmin varies from −3.2 to −4.1) with relative height ranging from 0 to 0.03.However, in the h p /H range from 0.03 to 0.11, the worst suction load on the roof decreases drastically (C p from −4.1 to about -0.9) as the parapet height increases.Additionally, the highest mean suctions have been found to be nearly invariant when relative height h p /H > 0.11.It should also be noted that, in contrast with roof model 6 Advances in Civil Engineering without parapets (C pmin � −3.2), solid parapets only produced an effective reduction of the wind loads (|C pmin | < 3.2) on the roof when h p /H ≥ 0.07.
To better understand the link between surface pressures and vortex behavior, Figure 9 displays mean pressure distributions on the flat roof exposed to conical vortices, together with mean suction profiles along lines normal to the leading edge.In this figure, numerical results of six solid parapets with relative heights ranging between 0 and 0.13 are presented.As is well known, the formation of conical vortices over a flat roof in oblique flow produces a symmetrical bell-shaped suction curve on the roof surface, with the maximum directly beneath the vortex core [36,37].As a result, the mean vortex core angle can be approximately ascertained as the angle formed between the leading edge and the rays of maximum mean suctions.
For parapets with h p /H � 0∼0.07, as shown in Figures 9(a)-9(e), the highest mean suction generally appears at a point closer to the roof corner than the mean vortex core, and it decreases with increasing distance from the corner.Mean suctions profiles are approximately bellshaped, and the generation and development of conical vortices are closely related with the parapet height.It can be observed that the size of the conical vortices gradually increases with increasing h p /H, whilst the vortex core angle increases from 12 °to 19 °.Additionally, smaller parapet height (h p /H ≤ 0.03) could be less favorable to reduce the highest mean suction induced by conical vortices.e parapet with h p /H � 0.03 (Figure 9(c)) tends to create a more pronounced high-suction region at the leading edge of the roof.Larger reductions in |C pmin | can only be obtained as parapet height increases.
is may be explained by the generation of eddy-turbulence wake downwind the parapets, which weakens the intensity of the vortex.In the case of h p / H 0.13, as shown in Figure 9(f ), it is interesting to note that surface pressures are uniformly distributed near the roof leading edge regardless of parapet height, and |C pmin | can be drastically reduced. is is because higher parapets tend to lift the conical vortices away from the roof surface, and consequently intense shear layers no longer interact with the roof surface.

Optimal Design of Porous Parapets.
Currently, the porous parapets have been consistently found to be the most e ective parapet con guration for mitigating the corner suctions on a at roof [10].erefore, six case studies of porous parapets (h p /H 0.01, 0.03, 0.05, 0.07, 0.09, and 0.13) are considered, as the application of porosity optimization framework proposed in this work.e gradient algorithm was employed to nd out the optimal porosity at a given parapet height, which provides the best performance for alleviating the high suctions on the roof surface.However, the main drawback of gradient algorithm is that it may converge to local minima.us, two initial porosities ϕ (0) 0% and 100% were used simultaneously to search the global optimal solution, which correspond to the lower and upper limits of the design variable, respectively.e optimization procedure will terminate when the change in objective function between iterations reduces below 1 × 10 −4 .
Figure 10 shows the variations with relative parapet height, h p /H, of the highest mean suctions resulting from solid parapets and parapets with optimal porosity.It can be seen that, in the h p /H range between 0.01 and 0.05, |C pmin | for parapets with optimal porosity are at least 26% smaller in magnitude than those for the solid parapet and also smaller than the no parapet case.For parapets with h p /H 0.01, 0.03, and 0.05, the optimal porosities are 52.3%, 47.7%, and 38.2%, respectively.Note that, when h p /H ≥ 0.07, porous parapets are almost useless in reducing high suctions on the roof surface, as compared to the solid parapet.
e optimal porosity for the optimization case of h p /H 0.07 is about 15.5%, but only 2.5% reduction in the magnitude of |C pmin | is obtained compared to the solid parapet case.It should also be noted that solid parapets are found to be the most effective for mitigating corner suctions when h p /H ≥ 0.09, i.e., optimal porosity ϕ 0%, due to the weaken vortices are displaced away from the roof surface.
Figure 11 shows contours of mean pressure coe cients and mean suction pro les on the at roof obtained from the porosity optimizations, with relative parapet height ranging between 0.01 and 0.07.As shown in Figures 9(b)-9(e) and Figure 11, the patterns of mean pressure distributions, along

8
Advances in Civil Engineering with the formation and dimension of conical vortices, between solid parapets and optimized porous parapets are similar.However, porous parapets with optimal porosity are more e ective than solid parapets to reduce the high suctions appearing at the leading corner of the roof, especially when the relative parapet height is lower than 0.05.e di erence in corner suction reduction between solid and optimized porous parapets becomes smaller as the relative parapet height increases up to h p /H 0.07.

Conclusions
e current study introduces an aerodynamic porosity optimization procedure that combines gradient algorithm and 3D CFD simulations based on the RANS equations.e applicability of the proposed optimization procedure was investigated using the example of a at roof with porous parapets located at the roof edges.e goal was to determine an optimal porosity to minimize the high-suction loads on the roof surface caused by conical vortices.e contributions of this paper can be summarized as follows: (i) CFD simulations utilizing four turbulence closure models (standard k-ε, RNG k-ε, SST k-ω, and RSM) were conducted for the prediction of the windinduced pressure distributions on a 3D at roof with solid and porous parapets.e simulation results from the RSM model were in relatively good agreement with the wind tunnel data by

Advances in Civil Engineering
Stathopoulos et al. [5].Additionally, the prediction results of the porous media flow using the pressure loss in the momentum equation were well agreed with the experimental results by Park and Lee [21].As a result, the numerical method was proved to be useful and appropriate to predict the effect of solid or porous parapets on the pressure distribution on the roof surface.(ii) Six optimization case studies of roof-edge parapets with h p /H ranging between 0.01 and 0.13 were presented to demonstrate the potential of the proposed optimization approach.It was found that for low-height parapets with h p /H � 0.01, 0.03, and 0.05, the highest mean suction is reduced by 26% at least for the optimal parapets (optimal porosities are 52.3%, 47.7%, and 38.2% respectively), as compared to the solid parapet, whereas for high-height parapets (h p /H ≥ 0.07), solid parapet seems to be most effective in eliminating the formation of conical vortices and reducing the highest mean suction on the roof surface.Generally speaking, the aerodynamic porosity optimization efficiency coupled with the extensive development in computational capacity is expected to encourage designers to seek for more optimal solutions while alleviating the high suctions appearing on the roof.
e main limitations of this work have been the restricted free-stream flow conditions and limited design height of the parapet.Further study needs to be conducted to evaluate the performance of the optimal parapet designs under different turbulent boundary conditions, and to explore a multiparameter design space to consider both the porosity and height of the parapet.

15 Figure 1 :Figure 2 :
Figure 1: Vertical pro le of the mean velocity at the inlet boundary.

Figure 4 :ϕFigure 5 :
Figure 4: Schematic sketch of the computational domain and coordinate system of the fence and prism.

Figure 3 :
Figure 3: Experimental and numerical mean pressure coefficients along the central line of the flat roof.(a) h p /H � 0, no parapet.(b) h p /H � 1/15.

Figure 6 :Figure 7 :Figure 8 :
Figure 6: Conical vortices on a flat roof and corresponding mean suction profiles in oblique flow.(a) Conical vortices.(b) Mean suction profiles.

Figure 9 :
Figure 9: Contours of mean pressure coe cients and mean suction pro les on at roofs with di erent relative parapet heights.(a) h p /H 0, no parapet.(b) h p /H 0.01.(c) h p /H 0.03.(d) h p /H 0.05.(e) h p /H 0.07.(f ) h p /H 0.13.