Issue 
Sust. Build.
Volume 2, 2017



Article Number  8  
Number of page(s)  19  
Section  Modelling and Optimisation of Building Performance  
DOI  https://doi.org/10.1051/sbuild/2017006  
Published online  03 October 2017 
Research Article
CFD optimisation of a stadium roof geometry: a qualitative study to improve the wind microenvironment
^{1}
Department of Mechanical Engineering, University of Sheffield,
Sheffield, UK
^{2}
Department of Mechanical and Industrial Engineering, Qatar University,
Doha, Qatar
^{*} email: ben.hughes@sheffield.ac.uk
Received:
16
June
2017
Accepted:
3
July
2017
The complexity of the built environment requires the adoption of coupled techniques to predict the flow phenomena and provide optimum design solutions. In this study, coupled computational fluid dynamics (CFD) and response surface methodology (RSM) optimisation tools are employed to investigate the parameters that determine the wind comfort in a twodimensional stadium model, by optimising the roof geometry. The roof height, width and length are evaluated against the flow homogeneity at the spectator terraces and the playing field area, the roof flow rate and the average interior pressure. Based on nonparametric regression analysis, both symmetric and asymmetric configurations are considered for optimisation. The optimum design solutions revealed that it is achievable to provide an improved wind environment in both playing field area and spectator terraces, giving a further insight on the interrelations of the parameters involved. Considering the limitations of conducting a twodimensional study, the obtained results may beneficially be used as a basis for the optimisation of a complex threedimensional stadium structure and thus become an important design guide for stadium structures.
Key words: computational fluid dynamics / response surface methodology / optimisation / wind comfort / stadiums
© P. Sofotasiou et al., published by EDP Sciences, 2017
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Stadium structures demand particular care during design process to comply with the manufacturing and environmental standards, due to the large number of attendees that they tend to attract. For this reason, progressive design strategies are ordinarily deployed in conceptual phase, in order to combine all the requisite components that will advocate for the optimality of the constructional result. The optimisation of every design should mainly focus on the improvement of users' comfort within the building environment, while aiming at the harmonious coexistence of an environmentally susceptible and climate rational design [1].
This study is focused on the design optimisation of the roof geometry of a twodimensional (2D) stadium model using coupled computational fluid dynamics and response surface methodology optimisation techniques. Based on rational design interventions on the roof configuration, the work pursues the provision of improved or optimum microclimate conditions in the arena bowl, including both spectator terraces and playing field area. The qualitative study includes the exploration of both symmetric and asymmetric roof alterations and the corresponding analytical results on velocity measurements and flow homogeneity at the areas of interest.
2 Previous related work
Early studies on stadium structures highlighted the multidimensional and multiscalar characteristics of thermal comfort and energy consumption in stadiums and accentuated the need to adopt coupled simulation software techniques to obtain a holistic design approach [2,3]. Since then, several studies have been performed to determine the structural characteristics that contribute the most to the alteration of the microenvironment at the stadium bowl.
CFD studies validated with wind tunnel experiments have shown that the roof configuration and the porosity of the stadium envelop have great impact on the wind distribution at both pitch level and spectator terraces and consequently on the induced wind comfort conditions [4–7]. Additional studies revealed the importance of providing ventilation openings to promote natural airflow and ventilation at the occupied areas of the stadium bowl [6,8,9]. The interactions between the urban environment and the stadium structures as well as the wind comfort conditions on pedestrian level have also attracted great interest [10–12]. Finally, thermal comfort studies, employing dynamic thermal modelling and CFD techniques, were also conducted to assess thermal discomfort and temperature distribution on stadia [2,13–15,41] and energy consumption requirements [16,17,42].
To date, only few of the aforementioned research works have been produced with intention to provide a design solution with improved responses towards wind comfort. The completed parametric studies have been conducted under limitations of design parameters and experimental variations. In this study, the design exploration of various roof configurations, within a predefined twodimensional design space, is pursued to provide improved wind environmental conditions and understand the interrelationships of the design parameters involved. To the authors' knowledge, no previous work has used combined CFD and RSM methods to assess and optimise a stadium or its components. This study will detail the method used, which will be useful for optimising stadium structures, since the method can be replicated for the optimisation of any roof and also other structural components.
3 Methodology
For the purpose of this study, coupled CFD modelling and RSM optimisation techniques were employed. CFD is a widely recognised and welldocumented simulation tool that employs numerical methods and algorithms to produce highly accurate predictions of flow movement and distribution inside and outside building structures [18]. On the other hand, RSM is a metamodelingbased optimisation technique that seeks local and/or global optimality, via the evaluation of the assigned design parameters. With substantially limited investment in computational cost, it provides the opportunity to evaluate numerous design cases and assess the influential role of the designated factors.
The simulationbased optimisation study was conducted in four main steps. A 2D stadium design was initially evaluated against velocity and pressure distributions. A series of sequential simulation runs was then performed generating data, regarding predefined geometrical characteristics of the stadium roof (Step 1). The results were used to formulate a high precision mathematical model that would describe the problem and predict relationships among variables (Step 2). The optimisation process was performed according to the selected algorithm and the defined objective sets and thresholds, pursuing design optimality (Step 3). Finally, the proposed solutions were verified with simulation runs and compared with the predicted ones, in order to evaluate the solution and therefore the accuracy of the RSMbased optimisation (Step 4). The flowchart in Table 1 is a graphical representation of the steps followed for the optimumseeking stadium roof configuration to create improved microenvironmental conditions. The methodology followed is fully detailed and validated [19] and it is a computational inexpensive solution to provide quantitative and qualitative results of improved design solutions.
4 Case study description
The 2D stadium design was based on FIFA's technical requirements and recommendations for international games, with a seating capacity of 60 000 [20]. The dimensional characteristics indicate a structure of 60 × 269 m^{2} (H × L). The oculus roof configuration is represented by a 95 m opening, long enough to provide adequate sunlight for the natural grass growth. The stadium consists of two spectator terraces, which are fully covered by the roof configuration, to reduce the direct solar irradiance (Fig. 1). They are also placed within a distance of 11 m far from the playing field to lessen the impact of the expected central vortex in the stadium bowl [5]. Furthermore, a ventilation opening of 1 m was created between the roof and the upper spectator tiers to allow for natural ventilation [6].
Fig. 1
Dimensional characteristics of the 2D stadium model. 
5 CFD settings and parameterisation
For the simulation of the wind distribution, the standard k–ε turbulence model was used, which is accepted and validated for 2D case studies [21–24].
5.1 Geometry and mesh generation
Regardless of the twodimensional geometry, the recommendations for CFD simulations for urban wind studies were followed [25]. Given the height of the stadium structure (H = 60 m), the computational domain was spanning 300 m (5 H) upstream, 900 m (15 H) downstream and 360 m (6 H) on the vertical direction. Thus, the modelled domain would allow for fully developed and homogeneous ABL flow prior approaching the stadium structure, as depicted in Figure 2.
The complexity of the geometry and the large difference in length scales complicated the generation of a fully structured mesh. The highresolution mesh close to the stadium walls, the ground walls and the ventilation openings was achieved by local size control on bodies of influence. As shown in Figure 3, four bodies of influence were created to facilitate the orderly evolution from coarser mesh in the outer domain to a finer mesh in the stadium area. Starting with a cell size of 8 m in the outer domain, the mesh was reduced to 2.5 m cell size for Body 1, 1 m for Body 2 and 0.5 m for Body 3 that included the area on the exact proximity of the stadium. Body 4 included the ground area with a cell size of 0.4 m. On the stadium walls a grid refinement control of 3 was applied. This method ensured the smooth transition from the wind urban scale to the building scale. The hybrid mesh consisted of 237 703 quadrilateral and triangular cells.
Fig. 2
Dimensional characteristics of the 2D computational domain. 
Fig. 3
Creation of bodies of influence and generation of controlled meshing size. 
5.2 Grid verification
Three different mesh sizes were tested to ensure mesh independency of the solution. The initial mesh consisted of 237 703 cells, the coarser mesh with 122 709 cells and the finer one with 576 620 hybrid cells. The selection of the favourable mesh size relied on the correlation of the dimensionless values of velocity (U/U_{freestream}) and the longitudinal direction of the flow (x/L_{Stadium}), along a line passing through the ventilation openings. The correlation of the resultant values is presented in Figure 4 and show good prediction for all three meshes.
Additional comparative results were generated for the flow rate at the front ventilation opening. As presented in Table 2, the maximum error of 5.4% was calculated for the coarse grid size, as opposed to 1.8% for the fine mesh. Thus, the medium mesh was selected for the current simulation study, allowing computationaltime savings and a high degree of solution accuracy.
Fig. 4
Grid sensitivity study among the three different mesh sizes; coarse mesh (122 709), medium mesh (237 703), fine mesh (576 620). 
Flow rate error for the three grid sizes.
5.3 Boundary conditions
The simulation of the wind distribution includes the generation of an atmospheric boundary layer (ABL). The mean ABL velocity profile follows a logarithmic production that when it is fully developed it can be described by equation (1) suggested by [26], followed by the mean turbulent kinetic energy (2) and dissipation rate profiles (3): (1) (2) (3) where is the ABL friction velocity equal to 0.77 m/s calculated from the reference velocity U_{ref} (7 m/s) at reference height y_{ref} (22 m), obtained from the meteorological station of the AlKhor city in Qatar, y_{0} is the aerodynamic roughness length (0.5 m), κ is the von Korman constant equal to 0.42, and C_{µ} is an empirical derived constant, equal to 0.09 in the standard k–ε model.
For the development of the logarithmic inlet velocity profile a userdefined function (UDF) was produced based on equations (1)–(3). It was also decided to improve the case model by applying constant shear stress on the top row of cells equal to [26], which is also suggested by [27]. The parameterisation of the computational model is summarised in Table 3. The roughness height, k_{s}, was calculated based on the equation: k_{s} = 9.793 × y_{0}/C_{s}. The software does not allow the k_{s} value to be greater than y_{p} that is the distance from the ground to the centre of the first cell. If this is the case, k_{s} is automatically set equal to y_{p}. In the current work, the first cell is set equal to 0.4 m (see Body 4, Fig. 3), resulting in y_{p} = 0.2 m, and thus k_{s} is set equal to 0.2 m. Consequently, the roughness constant C_{s} is greater than 1, and thus a UDF is required. The playing field area and the stadium entrances have an aerodynamic roughness equal to 0.03 m (for grass covered area, according to Wieringa [28]), ensuing roughness height equal to 0.58 m, with a set value of C_{s} equal to 0.5.
Solver settings for the 2D ABL flow.
5.4 CFD visualisation
Building on the initial results of the CFD simulation, velocity and pressure outputs allowed the investigation of: (1) the areas that are mostly influenced by the induced wind, (2) the wind speed distribution, (3) the wind behaviour and (4) the potential interventional strategies that could be adopted to improve the stadium's microenvironment. As illustrated in Figure 5, the roof and side ventilation openings seem to have a dominant role on the stadium's microclimate airflow. The incoming air from the roof opening leads to the generation of a primary vortex, located at the centre of the stadium bowl, and clearly shifted towards the leeward terraces. The primary vortex feeds a secondary counterrotating smaller vortex located at the windward spectator tiers, which forces the air to exit the stadium through the front ventilation opening. Additionally, the recirculation zone formed at the rear side of the stadium contributes to the air entrance through the back ventilation opening. Thus, a third vortex is developed, which is spatially limited to the leeward terraces, due to the existence of the primary, and predominant in size, vortex.
The analysis of the pressures acting on the stadium surfaces verifies the wind behaviour at the stadium interior. As shown in Figure 6, the frontal area of the stadium has positive pressure values, as opposed to the rest of the fluid domain. The pressure values at the interior of the stadium vary between −86.8 Pa and −27.5 Pa. The acceleration of the external wind is promoted due to pressure difference, leading to forced directional airflow into the stadium bowl through the roof opening, and thus producing the primary vortex.
In the front of the stadium, positive pressure patterns are observed below the roof level. However, the roof geometry, the extended ventilation opening and the negative angle of attack decelerate the flow leading to adverse pressure gradient and flow separation (Fig. 6, point A). In order to maintain pressure equilibrium, the air exits the stadium via the front ventilation opening. In contrary, on the rear side of the stadium, with pressure values around −20 Pa, the air is forced to enter the stadium interior via the back ventilation opening.
Fig. 5
Normalised velocity vectors and velocity contours of U/U_{ref}, where U_{ref} = 7 m/s at 22 m height. 
Fig. 6
Pressure contours and velocity vectors for the area outside and inside of the stadium. 
5.5 CFD validation
The validation of the CFD simulation for the 2D stadium structure was performed based on the study of Persoon et al. [5], who investigated the wind flow patterns and winddriven rain, regarding different roof structures. The development of a primary central vortex is observed for every case study and its size, and/or penetration level under the leeward spectator terraces, is highly dependent on the roof configuration. Secondary vortices on the windward and leeward spectator tiers were also observed in every design case and it was concluded that the size and length of each one of them is subjective to the geometrical characteristics of the stadium itself and the boundary conditions applied (Fig. 7, top).
Similar results were obtained for the current design case. Figure 7 (bottom) illustrates the interaction of the two vortices on the leeward side of the stadium. The primary vortex, created by the air stream entering through the roof opening, extends vertically over the entire height of the stadium bowl. It is well observed that the primary vortex prevails over the secondary one, penetrating about 20 m below the roof structure (Fig. 7, bottom left). Moving further downstream the roof, the secondary vortex has a more dominant role, directly affecting the wind behaviour on the spectators' level (Fig. 7, bottom right). The rear ventilation opening significantly contributes to the generation of the secondary vortex, allowing the pressure difference to naturally move air to the interior. Longwise comparisons of the dimensionless horizontal component of velocity indicate substantially lower velocity magnitudes to be developed by the secondary vortex, as opposed to the primary one.
Fig. 7
Wind distribution on three different roof configurations, and corresponding graphs of the dimensionless horizontal component of velocity along the vertical direction [5] (top); interaction between the primary and the secondary vortex, based on the horizontal components of velocity magnitude in two vertical lines (bottom); roof edge level (left) and spectator tiers level (right). 
6 Design of experiments
Design of experiments (DoE) is a statistical tool, used to collect and organise data generated from a limited number of experiments. This method enables to quantify how process inputs affect process outputs and study simultaneously several variables in an effort to determine their effects and their interactions. It is the intermediary for the RSM to construct the exploration region of responses, an empirical graph illustrating the response changes, as the continuous variables change under its control.
The three main principles for the implementation of a DoE were described in [29] as follows:

identify the performance assessment of the dependent variables (outputs) to be optimised, considering the influential role of the most significant independent variables (inputs);

designate the design space (constrains), within which the input parameters will vary;

generate the design points and their response values.
6.1 Selection of input parameters
According to literature, stadium structures have been extensively studied in order to determine optimum geometrical configurations to provide comfort conditions for both players and spectators [4,5,30]. The roof dimensional characteristics seem to be the major contributors on the attainment of microenvironmental satisfaction. Thus, the input design parameters were set as: (1) the roof height, (2) the roof radius, and (3) the roof length extension towards the centre of the pitch area.
6.2 Selection of output parameters
Favourable microenvironmental conditions in semioutdoor stadium structures, initially reported by [30], indicated that aerothermal comfort conditions may be achieved with wind threshold speed values between 0.1 m/s, to ensure minimum air movement, and 3.6 m/s, to prevent wind disturbance and performance impairment. On the basis of these considerations, the output parameters were set as the flow homogeneity values on the spectator tiers and the playing field area, the ventilation flow rate through the roof opening, and the overall pressure in the stadium bowl. The maximum and minimum velocity at the occupied areas and the ventilation flow rates through the leeward and windward ventilation openings located between the roof and the spectator terraces were selected as additional derived parameters. In order to control and enable the generation of the results, the areas of interest were separated from the entire fluid domain. As shown in Table 4, zones 2, 3 and 4 comprise an area equal to the zone length times 2 m on the perpendicular direction, to adequately cover the fluid area on human level. Zone 1 was 0.5 m thick, similar to roof thickness.
The flow homogeneity for ventilated areas was calculated based on the equation of the coefficient of spatial variation, presented in equation (4) [31]: (4) where C_{sv}, coefficient of spatial variation, C_{v}, average velocity coefficient obtained by equation (5): (5)σ_{s}(V_{i}/V_{r}) is the standard deviation of C_{v}, V_{i} is the mean velocity at interior location i (m/s), V_{r} is the mean outdoor reference velocity at height 22 m (7 m/s), and n is the number of nodes measured in the area of interest.
The coefficient of spatial variation, C_{sv}, was used as the indicator of the homogeneity at the zonal areas. Values of C_{sv} larger than 1 imply distributions of highvariance, as opposed to values lower than 1 that indicate flow uniformity and velocity magnitudes of lowvariance.
The ventilation flow rate was calculated based on equation (6), for a nonuniform flow of an incompressible fluid (ρ = constant), in order to obtain a direct indicator of the flow direction, either incoming or outgoing flow stream to the stadium bowl. (6) where A is the area of the opening configuration and is the mean velocity in the perpendicular direction to the opening area (A).
In Figure 8, the values of the output parameters for the initial stadium design are presented. The velocity values at the spectator tiers are lower than the velocity threshold of 3.6 m/s, as opposed to the playing field area. The calculated value of 4.41 m/s is 18 % higher than the recommended one that may result to players' performance disruption and compromise. The roof flow rate, with a value of 0.67 m^{3}/s indicates that the air is leaving the stadium area, but this is due to the generated primary vortex that recirculates the air at the stadium interior.
Separation of the fluid domain into four controllable zonal areas.
Fig. 8
Initial 2D stadium structure; illustration of velocity contours, values of homogeneity and maximum velocity speeds and roof flow rate. 
6.3 Generation of design points
The design space, within which different design points are evaluated, is defined as the design space that indicates the range of inputs' variability [32]. In the case of the stadium design, six input design parameters were selected, in order to control both symmetric and asymmetric transformations of the roof configuration, as presented in Table 5. The main difference between the two case studies was that in symmetric case, the input parameters were reduced by three, since P2, P4 and P6 were derived parameters, equal to P1, P3 and P5, respectively. It is worth highlighting at this stage that the design constraints set for the input parameters were subjective, if not arbitrary, selected. The design exploration excluded the option of convex roof configuration, or a roof length that exposes the spectator tiers or partly covers the pitch area.
The selection of the appropriate scheme for the generation of the design points was performed with the Optimal SpaceFilling design (OSF) with nonparametric regression model, because it can describe more accurately the nonlinear model parameters.
Description of the input parameters and the defined constrains for spatial variance.
6.3.1 OSF scheme
The OSF scheme is usually employed in complex metamodelling techniques, such as the nonparametric regression, because it offers an extensive overview of the design space, by equidistributional selection of the input design parameters. The Max–Min Distance design type was selected, which maximises the minimum distance between any two points and it saves computational time. The number of the design points was based on full quadratic samples, by providing the number of the design points needed to generate a full quadratic model [33].
Based on the OSF scheme, 15 and 28 design points were generated in total, for the symmetric and asymmetric cases, respectively. The calculation of the design point responses (output parameters) was performed with sequential simulation runs, contributing to the most computational timeconsuming part of the study, since convergence was achieved for every simulation run. The values of the output parameters will form the basis for the nonparametric regression model, as part of the RSM metamodel technique and the creation of the design space, where different roof configurations, regardless of the input value, may be examined against their responses.
7 Response surface methodology
The response surface methodology (RSM) technique employs mathematical and statistical tools to create approximation functions that can describe complex simulation data. Metamodel algorithms and mathematical relations are used to study the input and output values, their behaviour and interrelations, assisting model simplification, results' interpretation and performance of optimisation studies [34]. RSM determines the relationships among variables by performing regression analysis. Depending on the number of the variables related to the problem, different algorithms may be applied, in order to determine the most accurate regression function that will describe the model. In the current study, for both symmetric and asymmetric roof alterations, the nonparametric regression metamodel algorithm was employed. The regression model enabled the prediction of the highly nonlinear behaviour of the stadium's microclimate, with regard to the defined inputs.
7.1 Nonparametric regression
The nonparametric regression metamodel is based on the categorisation of noisy data using the Support Vector Machine (SVM) technique. SVM enables nonlinear mapping of input–output data (or training data) on a higherdimensional feature space by the aid of kernel maps. The data are split in nonoverlapping classes by the aid of hyperplanes, aiming at the maximization of their inbetween margin. The hyperplanes are used to separate noisy data in groups, by using only those inputs, also known as support vectors that seem to sufficiently represent the outputs. At the same time, as illustrated in Figure 9, the hyperplanes are suitably positioned, so as to minimise the distance of the support vectors that fall outside the margin (slack variables ξ).
Nonparametric regression systems with nonlinear responses and noisy results can be described mathematically by the following equations (7)–(9) [36]: (7) where f(x) is the response of interest, 〈⋅ , ⋅ 〉 is the dot product in x and w is a weight factor given by equation (8): (8)
Threshold b can be described by equation (9): (9) where , α_{i} ≥ 0 are Lagrange variables and x_{r} and x_{s} are support vectors.
The main advantage of nonparametric regression against the parametric one is that the function lacks of rigidity and can be adapted to the equivalent data without compromising form (parabolic, periodic, etc.), generating an inferred function [37].
Fig. 9
Methodology followed by SVM hyperplanes to separate data in classes and minimise the distance of points falling out of the margins [35]. 
7.1.1 Goodness of fit
The goodness of fit gave a first insight into whether the data are adequately represented by the function. As shown in Figures 10 and 11, the predicted values from the response surface and the observed values form the design points perfectly fit the diagonal line, indicating a very good fit of the model.
A commonly used indicator, to show in what percentage a design output can be described by the response surface regression equation, is the coefficient of determination (R^{2}). For linear regression analysis, the R^{2} can be described by equation (10), and lies in the interval [0,1]. (10)
However, in case of nonparametric regression and nonlinear data, the integrity of the function cannot rely on the coefficients of determination (R^{2}), because the more the parameters involved in the system, the higher the value of the regression sum of the squares will be, which might lead to a value of R^{2} greater than 1.0 in some cases (see Eq. (10)) [38]. Indeed, in current case the values of R^{2} for every design output were equal to 1.0, as the software automatically truncates exceeding values to 1.0.
Since nonlinear regression aims at minimizing the sum of squares of the points' distances from the fitting curve, the root mean square error (RMSE) can be considered a good estimator value, given by equation (11): (11)
The closer the value of RMSE is to zero, the better the fit. The values for the output values were varying from 2.2 × 10^{−9} to 1.7 × 10^{−4}, indicating a very good fit and a regression function adequate to represent the model.
Fig. 10
Goodness of fit based on the nonparametric regression for the symmetric case. 
Fig. 11
Goodness of fit based on the nonparametric regression for the asymmetric case. 
7.2 Goal driven optimisation (GDO)
The investigation of the optimum design solution presupposes the designation of the model objectives and constrains. The optimisation goal is to improve the stadium aerocomfort conditions, by providing homogeneity and low wind velocities at the zonal areas 2–4 (Tab. 4), and prevent the high flow rates via the roof opening. For the purpose of that, the optimisation targets set for both symmetric and asymmetric cases are presented in Table 6.
The homogeneity (P8, P9 and P11) in all three occupied areas was set with a minimization goal. The low values of C_{sv} (see Eq. (4)) will indicate lower flow heterogeneity and velocity variance. The average pressure at the stadium interior (P10) was deemed to be maximised, in order to prevent the large pressure differences between the external and the internal environment. The negative pressure at the stadium interior introduces large air streams flowing through the roof opening and the rear ventilation opening, disturbing the stadium's microclimate. At the same time, a negative flow rate through the roof (P7), calculated based on equation (6) using the vertical component of velocity, indicates that more air is entering the stadium bowl than exiting. Thus, a lower bound constraint of 0 m^{3}/s was added, in order to achieve air flow reduction.
The MultiObjective Genetic Algorithm (MOGA) was used to conduct the multiobjective optimisation study. The MOGA belongs to the family of genetic algorithms (GA). By employing Pareto sets, a number of design solutions are determined, after assessing multiple populations of points. The algorithm performs an elitist nondominated sorting according to the set objectives. The theoretical background can be found in [39,41].
Optimisation objectives and constraints set for the output parameters.
8 Results and discussion
8.1 Response surfaces
The response surface methodology gives a plethora of results, regarding the correlation of the input and output parameters, the sensitivity factors, indicating on what extent the input parameters influence the outputs, and the design spaces (or response surfaces), where any parameter can be evaluated against the responses.
Threedimensional graphs of the responses were generated, enabling the investigation of the effect that two inputs have on a single output. Due to the multidimensional character of the system, few response surfaces are selectively presented that show the complexity of the model. For the symmetric case study, the response surfaces of the flow homogeneity on the back spectator tiers (P11) and the play field area (P9) are presented in Figure 12, regarding the alterations in the roof length and radius (P1, P3) and the roof length and height (P3, P5), respectively. On linear correlations between the parameters, the response surface would fit perfectly the design points. However, this is not the case in the current design problem. The system is highly nonlinear and the algorithm adopted, used weighting factors to create a function of high accuracy that enables the prediction of any response within the input design space.
Based on the local sensitivities, by increasing the roof height, the roof flow rate and the homogeneity in the front spectator tiers and playing field area are decreased, as opposed to the average pressure and the homogeneity on the back spectator tiers that are increased. Moreover, the average pressure (P8) is increased by increasing the roof radius, and the homogeneity on the back spectator tiers is also increased by increasing the roof length.
Similar results were obtained for the asymmetric design case. The even larger number of independent parameters, led to a multivariable system, responses of which can be found in Figure 13. In this case, the front roof height (P6) seems to be determinant for the majority of outputs. The increasing value of P6 leads to reduced roof flow rate and poor homogeneity on the front spectator tiers, as opposed to the average pressure and the homogeneity on the back spectator tiers, which are increased. Additional to these, as the back roof height increases, the roof flow rate is increased, and the average pressure and the playing field homogeneity are decreased.
Fig. 12
Response surfaces for selective input–output parameters of the symmetric case. 
Fig. 13
Response surfaces for two sets of input–output parameters of the asymmetric case. 
8.2 Optimisation results
The MOGA metamodel algorithm generated in total six candidate points that performed better based on the defined set of objectives. Table 7 shows the optimisation results generated for the symmetric roof modifications. The geometrical patterns look quite similar for all three candidate points. More specifically, the roof height was reduced by a maximum of 70% and a minimum of 57% for the CP3 and CP1 respectively. The roof radius of CP3 was increased to more than 565%, resulting in an almost flat roof configuration. Finally, the length was reduced between 68.8 m and 71.1 m, providing full protection to the spectator terraces and exposing adequately the pitch area.
The simulated velocity contours indicate how the roof modification may alter dramatically the environment at the stadium interior. The centre of the primary vortex was relocated vertically upwards, resulting in positive flow rate values. Contrariwise, the creation of a new recirculation zone above the rear roof geometry is now considered responsible for the air entering the stadium bowl with highspeed values that are weakened before reaching pitch level. The secondary vortex located at the back spectator tiers was moved towards the pitch level. The air entering through the back ventilation opening resulted in the generation of another vortex limited to the top spectator tiers, for CP1 and CP3. For the CP1, the same vortex is of lower intensity and spatial extension. In any design solution, these vortices promoted flow heterogeneity, resulting to higher values of C_{sv}. The change in roof height and radius resulted in reverse air direction on the front ventilation opening. The impinging air enters the stadium and partly feeds the present secondary vortex and partly enhances the primary vortex to exit the stadium area.
In Table 8, the candidate points for the asymmetric roof optimisation are presented. In order to avoid the air entrance at the stadium interior through the roof opening, the back roof height (P4) was decreased by a maximum of 91.25% for CP2 and a minimum of 61.87% for CP3. The difference in height allows the air to pass above the roof construction, without entering the stadium interior. The recirculation zone that is created behind the stadium structure drives the air through the rear ventilation opening and at the same time contributes to the generation of a recirculation zone above the back roof part. In this case, two factors are responsible for the generation of the primary vortex. The first one is due to the aforementioned recirculation area at the top of the roof surface and the second one is due to the incoming air through the rear ventilation opening. Similarly, as in symmetric case, the secondary vortex in the back spectator tiers is relocated towards the pitch area, providing lower velocity values at zone 3, but slightly higher heterogeneity. The flow homogeneity in the pitch area has been increased by a minimum of 132%, whereas in the front spectator tiers remained relatively stable over the same values. Finally, the front secondary vortex is spatially limited, by the prevailing primary one and with considerable lower velocity magnitude intensity.
Candidate points generated based on MOGA algorithm for symmetric case.
Candidate points generated based on MOGA algorithm for asymmetric case.
8.3 GDO sensitivities
Using correlation analysis, the global sensitivities over the output parameters were generated, by evaluating the entire design space. Figure 14 (top) shows the global sensitivities for the symmetric case study. According to the results, the roof height and length have great impact (either positive or negative) in all the design responses as compared to roof radius. Positive sensitivity denotes that as the value of the design parameter increases, the value of the design response will increase too. Negative sensitivity implies that increasing the value of the design parameter, the response value will be decreased. The roof radius appears to have positive sensitivity in both the average pressure at the stadium interior and the values of homogeneity at the back spectator tiers. It also has negligible impact on the parameters P7, P8 and P9.
In case of the asymmetric roof configuration (Fig. 14, bottom), the front roof height (P6) seems to have the greatest impact both positive and negative in almost all of the output parameters. Considering wind is unidirectional, it is expected that P6 will determine mostly the flow at the stadium interior. In case of the homogeneity in the playing field area, P6 seems to have no impact. This can be explained graphically by the figures in Table 8, where a recirculation zone is created on top of the back roof part and drives the air at the stadium interior, and towards the pitch area.
Fig. 14
Global sensitivity chart of input–output parameters for the symmetric (top) and asymmetric (bottom) case. 
8.4 GDO verification
The verification of the candidate points for optimum aerocomfort conditions at the stadium microenvironment was performed by CFD simulation runs. Figure 15 compares the results of the candidate points obtained from the GDO and the CFD simulations for the symmetric roof modification. The values of homogeneity in the zonal areas 2, 3 and 4 (P8, P9 and P11) show good prediction performance, with a maximum error of 16.35% associated with the homogeneity values at the back spectator tiers (P11) for the CP2, followed by an error of 11.04% for the same output parameter for the CP1. The development of vortices at the windward spectator tiers is susceptible to disturbances caused by the induced wind through the ventilation openings, and thus the computation of velocity magnitudes is highly sensitive. The flow rate values (P7) were predicted with a maximum absolute deviation of 1.84 m^{3}/s for the CP2. Considering the roof length opening of 116.5 m and the presence of the central primary vortex, the difference in flow rate can be consider substantially low. The largest errors were observed for the average pressure (P10) for all three candidate points. With the maximum absolute error to be as high as 29.20 Pa, the pressure output could not be considered an accurate parameter to be considered for the selection of the optimum solution. The generated recirculation zones and wind vortices make the average pressure imprecise and difficult to predict.
In Figure 16, the comparison of the results for the asymmetric roof indicated a GDO of higher performance. The homogeneity values were predicted with a maximum error of 8.65% for the CP3. The maximum absolute deviations of the average pressure (P10) and the flow rate (P7) were equal to 12.4 Pa for the CP1 and 3.4 m^{3}/s for the CP3, respectively.
The verification study has highlighted the inaccuracies in predicting some output parameters, such as the average pressure (P10) and the flow rate (P7), causing some sort of compromise in the optimisation process. Regarding the number of variables involved and the complexity of the model, the verification enabled the assessment of the best roof configuration among the candidates, revealing the global optimum design solutions. For the symmetric case study, the CP1 had the best performance according to the CFD simulation results. Likewise, the CP2 in the asymmetric roof configuration had the less deviation in all parameters against the CFD results.
Fig. 15
RSM vs. CFD results for verification for the three candidate points for the symmetric case. 
Fig. 16
RSM vs. CFD results for verification of the three candidate points for the asymmetric case. 
9 Interpretation of results
The verified optimum solutions are further discussed in terms of maximum velocity magnitudes and internal flow distribution. Figure 17 presents the normalised velocity vectors for the CP1 of the symmetric and the CP2 of the asymmetric roof configuration.
The roof optimisation contributed to a total reduction of 26.5%, 15.4% and 25.9% for the maximum velocities in the front, back spectator tiers and the pitch area for the symmetric case and 76.5%, 62.7% and 55.6% for the asymmetric case, respectively. Thus, the maximum allowed mechanical threshold of 3.6 m/s for stadia was satisfied for both designs. Figure 18 shows a graphical comparison of the dimensionless horizontal component of velocity for the front and back spectator tiers before and after optimisation. It can be concluded that the microenvironment at the spectator terraces is substantially less determined by the primary vortex and consequently by the induced air through the roof opening.
Fig. 17
Normalised velocity vectors and maximum zonal velocity magnitude for the optimum symmetric (top) and asymmetric (bottom) case. 
Fig. 18
Comparison of the dimensionless horizontal component of velocity for symmetric (top) and asymmetric (bottom) case, before (red line) and after (green line) optimisation. 
9.1 Limitations and suggestions for future studies
The twodimensional study poses certain limitations regarding the lack of representing adequately 3D flow characteristics, including turbulent dispersion, buoyancy effects and capture of the separation zones. However, the qualitative study gives a further insight on the design optimisation of the stadium roof geometry for improved aerocomfort conditions in different areas in the stadium. The optimisation method and results presented in this work aim to form the basis for a future study of a 3D stadium design. The definition of design constraints and objectives limits the exploration of solutions of global optimality. The study was additionally formed under subjective selection of input and output parameters. Thus, different design solutions may be obtained under different sets of parameters and objectives.
10 Conclusion
In the present study, a qualitative optimisation study of a 2D stadium roof configuration was performed, with the intention to improve the wind comfort conditions at the stadium bowl. Based on coupled CFD and RSM techniques, the selected design parameters were assessed with reduced computational time and cost, as opposed to traditional deterministic studies. The initial stadium design was simulated using the FLUENT code and validated against the study of Persoon et al. [5] for airflow distribution in football stadia. The investigation of different roof configurations, based on variable height, length and radius, was conducted using design experimentation. Both symmetrical and asymmetrical modifications were assessed against the microclimate conditions developed at the stadium interior, and 6 candidate points were generated in total. CP1 for the symmetric and CP2 asymmetric case were proven to have the best performance against the CFDbased verification.
The results revealed that the optimal symmetric roof design may contribute to a maximum velocity reduction of up to 26.5% for the front spectator tiers, if the roof height is reduced by 57% and the roof radius increased by 835%, over the initial one. In the asymmetric roof configuration, the front roof part seems to be the most crucial for the airflow distribution. With an asymmetric roof difference of 11.4 m the maximum reduction observed for the maximum velocity was equal to 75%. For both designs, wind comfort was achieved that complies with the wind thresholds set for stadiums. Considering the limitations of conducting a twodimensional study, the optimisation method and qualitative results presented in this work will form the preliminary basis for future optimisation study of a complex 3D stadium structure.
Acknowledgments
This paper was made possible by a NPRP grant (64612188) from the Qatar National Research Fund (A member of the Qatar Foundation). The statements made herein are solely the responsibility of the authors.
References
 A. Parr, M. Zaretsky, New directions in sustainable design (Routledge, London, 2011) [Google Scholar]
 K. Lomas, H. Eppel, M. Cook, J. Mardaljevic, Ventilation and thermal performance of design options for Stadium Australia, in: The 5th International IBPSA Conference [online] (International Building Performance Simulation Association, Prague, 1997), Available at: https://dspace.lboro.ac.uk/2134/12393 [Google Scholar]
 D. Fiala, K.J. Lomas, 1999, Application of a computer model predicting human thermal responses to the design of sports stadia, in: CIBSE'99, Conference Proc., Harrogate, UK, pp. 492–499 [Google Scholar]
 J. Bouyer, J. Vinet, P. Delpech, S. Carre, Thermal comfort assessment in semioutdoor environments: application to comfort study in stadia, J. Wind Eng. Ind. Aerodyn. 95, 963 (2007) [CrossRef] [Google Scholar]
 J. Persoon, T. van Hooff, B. Blocken, J. Carmeliet, M. de Wit, On the impact of roof geometry on rain shelter in football stadia, J. Wind Eng. Ind. Aerodyn. 96, 1274 (2008) [CrossRef] [Google Scholar]
 A. Szucs, S. Moreau, F. Allard, Aspects of stadium design for warm climates, Build. Environ. 44, 1206 (2009) [CrossRef] [Google Scholar]
 T. van Hooff, B. Blocken, M. van Harten, 3D CFD simulations of wind flow and winddriven rain shelter in sports stadia: influence of stadium geometry, Build. Environ. 46, 22 (2011) [CrossRef] [Google Scholar]
 P. Biagini, C. Borri, L. Facchini, Wind response of large roofs of stadions and arena, J. Wind Eng. Ind. Aerodyn. 95, 871 (2007) [CrossRef] [Google Scholar]
 T. van Hooff, B. Blocken, Fullscale measurements of indoor environmental conditions and natural ventilation in a large semienclosed stadium: possibilities and limitations for CFD validation, J. Wind Eng. Ind. Aerodyn. 104, 330 (2012) [CrossRef] [Google Scholar]
 B. Blocken, J. Persoon, Pedestrian wind comfort around a large football stadium in an urban environment: CFD simulation, validation and application of the new Dutch wind nuisance standard, J. Wind Eng. Ind. Aerodyn. 97, 255 (2009) [CrossRef] [Google Scholar]
 T. van Hooff, B. Blocken, Coupled urban wind flow and indoor natural ventilation modelling on a highresolution grid: a case study for the Amsterdam ArenA stadium, Environ. Model. Softw. 25, 51 (2010) [CrossRef] [Google Scholar]
 T. van Hooff, B. Blocken, On the effect of wind direction and urban surroundings on natural ventilation of a large semienclosed stadium, Comput. Fluids 39, 1146 (2010) [CrossRef] [Google Scholar]
 A. Stamou, I. Katsiris, A. Schaelin, Evaluation of thermal comfort in Galatsi Arena of the Olympics “Athens 2004” using a CFD model, Appl. Therm. Eng. 28, 1206 (2008) [CrossRef] [Google Scholar]
 M. Ucuncu, D. Woolf, M. Zikri, Thermal comfort of spectators in stadia built in hot climates, paper presented at Adapting to Change: New Thinking on Comfort, Windsor, UK, 9–11April (Network for Comfort and Energy Use in Buildings, London, 2010), p. 8 [Google Scholar]
 A. Matzarakis, D. Fröhlich, Sport events and climate for visitors—the case of FIFA World Cup in Qatar 2022, Int. J. Biometeorol. 59, 481 (2014) [CrossRef] [Google Scholar]
 T. Nishioka, K. Ohtaka, N. Hashimoto, H. Onojima, Measurement and evaluation of the indoor thermal environment in a large domed stadium, Energy Build. 32, 217 (2000) [CrossRef] [EDP Sciences] [Google Scholar]
 P. Sofotasiou, B. Hughes, J. Calautit, Qatar 2022: Facing the FIFA World Cup climatic and legacy challenges, Sustain. Cities Soc. 14, 16 (2015) [CrossRef] [Google Scholar]
 B. Blocken, 50 years of Computational Wind Engineering: Past, present and future, J. Wind Eng. Ind. Aerodyn. 129, 69 (2014) [CrossRef] [EDP Sciences] [Google Scholar]
 P. Sofotasiou, J. Calautit, B. Hughes, D. O'Connor, Towards an integrated computational method to determine internal spaces for optimum environmental conditions, Comput. Fluids 127, 146 (2016) [CrossRef] [Google Scholar]
 FIFA, Football Stadiums: Technical recommendations and requirements, 5th edition [ebook] (FIFA, 2011), p. 108, 283 [Google Scholar]
 F. Lien, E. Yee, Y. Cheng, Simulation of mean flow and turbulence over a 2D building array using highresolution CFD and a distributed drag force approach, J. Wind Eng. Ind. Aerodyn. 92, 117 (2004) [CrossRef] [Google Scholar]
 X. Xiaomin, H. Zhen, W. Jiasong, The impact of urban street layout on local atmospheric environment, Build. Environ. 41, 1352 (2006) [CrossRef] [Google Scholar]
 Y. Huang, X. Hu, N. Zeng, Impact of wedgeshaped roofs on airflow and pollutant dispersion inside urban street canyons, Build. Environ. 44, 2335 (2009) [CrossRef] [Google Scholar]
 L. Yik, S. Salim, A. Chan, C. Cheong, CFD study of flow over parallel ridges with varying height and spacing, in: World Congress on Engineering, London (2010), Available at: http://www.iaeng.org/publication/WCE2010/WCE2010_pp12061211.pdf [Google Scholar]
 F. Joerg, Recommendations of the COST action C14 on the use of CFD in predicting pedestrian wind environment, Fourth Int. Symp. Comput. Wind Eng. 2006, 529 (2006) [Google Scholar]
 P. Richards, R. Hoxey, Appropriate boundary conditions for computational wind engineering models using the kε turbulence model, J. Wind Eng. Ind. Aerodyn. 46–47, 145 (1993) [CrossRef] [Google Scholar]
 D. Hargreaves, N. Wright, On the use of the kε model in commercial CFD software to model the neutral atmospheric boundary layer, J. Wind Eng. Ind. Aerodyn. 95, 355 (2007) [CrossRef] [Google Scholar]
 J. Wieringa, Updating the Davenport roughness classification, J. Wind Eng. Ind. Aerodyn. 41, 357 (1992) [CrossRef] [Google Scholar]
 X. Shen, G. Zhang, B. Bjerg, Investigation of response surface methodology for modelling ventilation rate of a naturally ventilated building, Build. Environ. 54, 174 (2012) [CrossRef] [Google Scholar]
 A. Szucs, S. Moreau, F. Allard, Spectators' aerothermal comfort assessment method in stadia, Build. Environ. 42, 2227 (2007) [CrossRef] [Google Scholar]
 J. Kindangen, G. Krauss, P. Depecker, Effects of roof shapes on windinduced air motion inside buildings, Build. Environ. 32, 1 (1997) [CrossRef] [Google Scholar]
 M. Cavazzuti, Optimization methods (Springer, Berlin, 2013) [CrossRef] [Google Scholar]
 ANSYS, Design Exploration Release 12.1, no. November (ANSYS, Inc, Southpointe 275 Technology Drive Canonsburg, PA 15317, USA, 2009) [Google Scholar]
 L. Friedman, The simulation metamodel (Springer US, Boston, MA, 1996) [CrossRef] [Google Scholar]
 B. Schölkopf, C. Burges, A. Smola, Advances in kernel methods (MIT Press, Cambridge, MA, 1999), p. 11 [Google Scholar]
 V. Vapnik, The nature of statistical learning theory (Springer, New York, 2000) [CrossRef] [Google Scholar]
 M. Wand, M. Jones, Kernel smoothing (Chapman & Hall, London, 1995) [CrossRef] [Google Scholar]
 T. Kvålseth, Note on the R^{2} measure of goodness of fit for nonlinear models, Bull. Psychon. Soc. 21, 79 (1983) [CrossRef] [Google Scholar]
 K. Deb, A. Pratap, S. Agarwal, T. Meyarivan, A fast and elitist multiobjective genetic algorithm: NSGAII, IEEE Trans. Evol. Comput. 6, 182 (2002) [CrossRef] [Google Scholar]
 C. Foncesa, P. Fleming, Genetic algorithms for multiobjective optimization: formulation, discussion and generalization, in: Fifth International Conference on Genetic Algorithms (Morgan Kauffman, San Mateo, CA), 416 (1993) [EDP Sciences] [Google Scholar]
 S.F. Fatima, H.N. Chaudhry, Steadystate CFD modelling and experimental analysis of the local microclimate in Dubai (UAE), Sust. Build. 2, 5 (2017) [CrossRef] [EDP Sciences] [Google Scholar]
 F. Calcerano, C. Cecchini, L. Martinelli, Numerical analysis of passive strategies for energy retrofit of existing buildings in Mediterranean climate: thermal mass and natural ventilation combination, Sust. Build. 2, 4 (2017) [CrossRef] [EDP Sciences] [Google Scholar]
Cite this article as: P. Sofotasiou, B. Hughes and S.A. Ghani: CFD optimisation of a stadium roof geometry: a qualitative study to improve the wind microenvironment. Sust. Build. 2, 8 (2017).
All Tables
Description of the input parameters and the defined constrains for spatial variance.
All Figures
Fig. 1
Dimensional characteristics of the 2D stadium model. 

In the text 
Fig. 2
Dimensional characteristics of the 2D computational domain. 

In the text 
Fig. 3
Creation of bodies of influence and generation of controlled meshing size. 

In the text 
Fig. 4
Grid sensitivity study among the three different mesh sizes; coarse mesh (122 709), medium mesh (237 703), fine mesh (576 620). 

In the text 
Fig. 5
Normalised velocity vectors and velocity contours of U/U_{ref}, where U_{ref} = 7 m/s at 22 m height. 

In the text 
Fig. 6
Pressure contours and velocity vectors for the area outside and inside of the stadium. 

In the text 
Fig. 7
Wind distribution on three different roof configurations, and corresponding graphs of the dimensionless horizontal component of velocity along the vertical direction [5] (top); interaction between the primary and the secondary vortex, based on the horizontal components of velocity magnitude in two vertical lines (bottom); roof edge level (left) and spectator tiers level (right). 

In the text 
Fig. 8
Initial 2D stadium structure; illustration of velocity contours, values of homogeneity and maximum velocity speeds and roof flow rate. 

In the text 
Fig. 9
Methodology followed by SVM hyperplanes to separate data in classes and minimise the distance of points falling out of the margins [35]. 

In the text 
Fig. 10
Goodness of fit based on the nonparametric regression for the symmetric case. 

In the text 
Fig. 11
Goodness of fit based on the nonparametric regression for the asymmetric case. 

In the text 
Fig. 12
Response surfaces for selective input–output parameters of the symmetric case. 

In the text 
Fig. 13
Response surfaces for two sets of input–output parameters of the asymmetric case. 

In the text 
Fig. 14
Global sensitivity chart of input–output parameters for the symmetric (top) and asymmetric (bottom) case. 

In the text 
Fig. 15
RSM vs. CFD results for verification for the three candidate points for the symmetric case. 

In the text 
Fig. 16
RSM vs. CFD results for verification of the three candidate points for the asymmetric case. 

In the text 
Fig. 17
Normalised velocity vectors and maximum zonal velocity magnitude for the optimum symmetric (top) and asymmetric (bottom) case. 

In the text 
Fig. 18
Comparison of the dimensionless horizontal component of velocity for symmetric (top) and asymmetric (bottom) case, before (red line) and after (green line) optimisation. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.