Access the full text.
Sign up today, get DeepDyve free for 14 days.
Loading next page...
/lp/springer-journals/simulation-of-a-ground-mounted-prism-in-abl-flow-using-les-on-overview-9nZQ3TDQf0
- Publisher
- Springer Journals
- Copyright
- Copyright © The Author(s) 2023. corrected publication 2023
- eISSN
- 2524-6992
- DOI
- 10.1186/s42774-023-00139-z
- Publisher site
- See Article on Publisher Site
Abstract
jin.xing2@unibo.it The accuracy of wind loading predictions using Large Eddy Simulation (LES) is usually Department of Civil, Chemical, Environmental, and Materials influenced by numerous model parameters, which can influence the obtained results. Engineering (DICAM), The validation of numerical simulations with traditional Wind Tunnel Test ( WTT ) is still University of Bologna, Viale del an important task, necessary to increase our a priori knowledge of possible inaccuracies Risorgimento 2, Bologna 40136, Italy and set up mitigation strategies. In this study, LES is used to simulate the wind fields around an isolated model high-rise building, under seven wind attack angles and vali- dated with WTT results. The influence of various settings and parameters on the model performance is studied. For the angle of attack showing higher inaccuracy, different mesh refinement strategies and turbulence models are tested. Results indicate that LES can accurately predict the mean and local maxima of the pressure coefficients for both perpendicular and skew wind attack angles, as well as reproduce global forces and their envelopes with very good accuracy. Conversely, pronounced errors are found in the prediction of the pressure coefficient standard deviation and the local minima. The highest deviations between LES and WTT are found close to the leading edge in cor- respondence of flow separations which are observed in WTT and not in LES for skew flows. The addition of boundary layer cells and the use of different subgrid models have very limited effectiveness in modifying the obtained results for the analysed case. Keywords: Wind loading, CFD, High-rise buildings, Large eddy simulation, Aerodynamics 1 Introduction Large Eddy Simulation, LES, is known to be a powerful tool for the investigation of wind effects on structures such as high-rise buildings [1–4]. However, comparison with tradi - tional Wind Tunnel Test (WTT) sometimes shows discrepancies, so that their validation and the individuation of critical aspects which might affect results is still an important research activity. Numerous studies can be found in the literature validating LES results in terms of wake flow characteristics [5, 6], global aerodynamic forces [7, 8] and pressure distribution [9–11]. The usefulness of such validations mainly consists in the individu - ation of critical aspects which can be sometimes moderated by appropriate modelling choices. When moderation strategies are not found or lead to excessive computational © The Author(s) 2023. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the mate- rial. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. Xing Advances in Aerodynamics (2023) 5:9 Page 2 of 19 burden, their individuation anyway allows to build up knowledge, which can help in correctly interpreting results obtained by numerical models when WTT results are not available. From the studies above and analogous works [12, 13], it can be found that the model performance of LES is usually good but, in some cases, sensitive to various model parameters. Among them, the accurate reproduction of flow characteristics in the Atmospheric Boundary Layer (ABL) is known to be crucial for the evaluation of wind effects on structures [14, 15]. Nevertheless, numerous other factors can influence the simulation [16–18], which include, as pointed out in [19], ground roughness, subgrid scale model, domain size, near-wall treatment and sensitivity to mesh and timestep. Numerous contributions have been devoted to investigating the influence of some of the aforementioned model parameters on the model accuracy. Unfortunately, due to the variability of the adopted computational setup and analysed flows, contradictory results are sometimes found. In particular, Wang et al. [20] assessed the effects of the mesh type and boundary mesh on the time-averaged and fluctuating wind characteristics around an isolated high-rise building standing in ABL. They found that the case with hexahedral cells has the best agreement with experiments. They also found that the boundary layer mesh does not improve the numerical accuracy in any circumstance, despite the fact that the non- dimensional wall distance ( y ) is often considered to be a very important quantity to be controlled. Wijesooriya et al. [21] showed that the choice of the Sub-Grid Scale (SGS) model does impact the flow field. The wall adapting local eddy viscosity (WALE) model was found to be the most suitable among the tested models including Smagorinsky-Lilly model, Dynamic model, kinetic energy sub grid scale model and wall modeled LES. Liu et al. [22] examined the effect of mesh resolution and different SGS models on the pre - diction of the pedestrian level wind fields around building arrays. They observed that the dynamic Smagorinsky-Lilly method gave the best model performance and recom- mended appropriate values for time step and sampling period. The aim of the present paper is to assess the numerical accuracy of LES in predicting wind loads on a high-rise building and to study the influence of some important model parameters on the accuracy of the obtained results. The high-rise building model is an isolated flat-topped box with an aspect ratio of 1:3:5, standing in an ABL. It is one of the high-rise building models from the public wind tunnel database constructed by Tokyo Polytechnic University (TPU), Japan, which have been considered also in some of the already available contributions [19, 23–25]. Firstly, the numerical results are compared with WTT data in the form of several ◦ ◦ model performance metrics for 7 angles of attack ranging from 0 to 90 . Then, for the worst case, corresponding to a skew flow impinging at 45 , the spatial distribution of the pressure coefficient statistics (i.e. time-average, standard deviation and peak values) and the corresponding Prediction Error (PE), defined as the difference between numerical and experimental results, are analysed. For such worst case, different numerical setups are evaluated varying mesh configuration and testing three SGS models, namely, WALE, k − ω SSTSAS and Smagorinsky. The paper is organized as follows. The descriptions of WTT and LES are provided in the next section, along with the procedure used to generate inflow conditions. Then, the Xing Advances in Aerodynamics (2023) 5:9 Page 3 of 19 numerical results are validated in Section 3. Subsequently, the sensitivity studies of the PE to the mesh configuration and the SGS model are provided in Section 4. Finally, con- clusions are drawn in Section 5. 2 WTT and LES description 2.1 WTT description The WTT were performed in the Boundary Layer Wind Tunnel (BLWT) of TPU, Japan. It is an open circuit tunnel with 2.2 m width and 1.8 m height. The ABL in the experi - ments was generated through spires and square blocks acting as roughness elements. It corresponds to a terrain category IV in the Architectural Institute of Japan (AIJ) stand- ard [26]. The wind speed at building height is around 11.11 m/s, and the stream-wise turbulence intensity ( I ) is approximately 15%. u_WTT The tested model is a flat-topped prismatic box with 0.5 m height (H) and 0.3 m (L) ∗ 0.1 m (W) rectangular cross-section. The length scale reduction factor of the experi - ments is 1:400, leading to a building height equal to 200 m in real scale. 480 wind pres- sure taps are positioned on the four side faces, sampling synchronously with a frequency of 1000 Hz for a duration of 32.8 s. The origin of the coordinate system is the center of the building base. The geometry of the building is presented in Fig. 1 (a). ◦ ◦ ◦ Incidence angles ranging from 0 to 100 with an increment of 5 are considered in WTT. The 0 is defined as the condition in which the flow is moving along the +x direc- tion and it is perpendicularly hitting the long building side, as indicated in Fig. 1 (b). The complete database and detailed description of the WTT are publicly accessible [27]. 2.2 LES description LES are designed to reproduce the WTT using the open source software OpenFoam v6.0. A 3D view of the computational domain is presented in Fig. 2 (a). The domain has the dimensions of 3 m, 2.2 m and 1.8 m in stream-wise, lateral and vertical directions, respectively. The distance between the inlet patch and the building location is 1 m, which equals the double of the building height H. The blockage ratio of the simulations is approxi - ◦ ◦ mately 3.0% for 0 and is 1.0% for 90 , which are both lower than recommendations [28]. ◦ ◦ Fig. 1 The geometry of the building (a) and the wind incidence angles from 0 to 100 (b) Xing Advances in Aerodynamics (2023) 5:9 Page 4 of 19 Fig. 2 The computational domain with boundary conditions (a) and the overviews of the mesh distribution: ◦ ◦ 0 (b) and 45 (c) The turbulent flow at the inflow patch is generated using a synthetic turbulence gen - eration technique called PRFG . It is an extended version of the Prescribed-wavevector Random Flow Generator (PRFG) [29], aiming at giving control over three-dimensional Xing Advances in Aerodynamics (2023) 5:9 Page 5 of 19 spectral densities, in other words, explicitly targeting all integral length and time scales of turbulence. Interested readers are invited to the references [29, 30] for details regard- ing the method. The top and two sides of the domain are treated as symmetry condi - tions, mimicking the effect of wind tunnel walls. The domain bottom and the building surfaces are treated as walls. In particular, the ground of the domain is modeled as a rough wall with a roughness height of 0.015 m (model scale) and a roughness constant of 0.5. The building surfaces are modeled as smooth walls. An inlet-outlet condition is used for the outlet patch. An overview of the mesh is shown in Fig. 2 (b). The snappyHexMesh is used for mesh - ing. A rotor-stator method is adopted. Specifically, the rotor is a cylinder of radius 0.25 m (reduced scale) with an axis oriented along the z-direction passing through the build- ing center. It is rigidly rotated when studying the wind effects at different attack angles without re-meshing the geometry. The stator is the rest of the domain, as indicated in Fig. 2 (b) and (c). As can be seen, a hexahedral mesh is used for the whole domain and cell expansion gradients are used, leading to the fact that the mesh size at the domain top and outlet is around two times bigger than at the domain center. The mesh size is around 0.016 m at the inflow patch and it is around 0.004 m up to a distance of 0.1 m from the building surfaces. This yields the number of cells for the building H, L and W to be 125, 75 and 25, respectively, which fulfills the requirements in [18]. The total number of cells is around 1.8 million. As for the numerical schemes utilized for the simulations, the pressure velocity cou- pling is imposed using the PISO algorithm. The time discretization is performed using the Crank Nicolson scheme with a blending coefficient of 0.85. The face fluxes for all the quantities are calculated using the bounded Gauss linear scheme except for the velocity which is calculated using the second-order LUST scheme. The k-Equation is used as the LES subgrid scale turbulence model. The non-dimensional time step ( �t ∗ (U /H) ) for the calculation is 0.01, yielding the Courant number of the calculation to be around 3.0 on average and only attain 6.0 in some small size cells. 2.3 Inflow As mentioned before, synthetic turbulence is applied at the inflow. In particular, PRFG is used for generating the ABL flow. It targets a Category IV profile in the AIJ standard code, which is in agreement with the one adopted in WTT. The mean velocity and the along-wind turbulence intensity, I , distribution along height are set following the power law in the code prescription, being 0.27 the exponential and 550 m the reference height of ABL. The turbulence intensities in the y and z direc- tions are set to be 0.75 and 0.5 times of the turbulence intensity in the x direction, respectively. We mention that in order to obtain the target wind field characteris- tics at the building location, the turbulence intensity of each velocity component has been increased of about 20% at the inflow to compensate for the energy dissipation between the inflow and the building location. The reference velocity, U , is chosen as the time-averaged wind speed at the reference height, H , which is the height of the building. In order to check the performance of the turbulence generator, an empty domain with the synthesized turbulent inflow condition is tested first. The wind profiles at Xing Advances in Aerodynamics (2023) 5:9 Page 6 of 19 the location where the building will be placed are shown in Fig. 3. As can be seen, the mean wind speed profiles from the AIJ code, WTT and LES agree well with each other. In the same figure the turbulence intensity profiles are reported, being I , u_LES I and I at H equal to 12%, 10% and 8%, respectively. Furthermore, the sim- v_LES w_LES r ulated turbulence integral length scales at H , evaluated from velocity time-histories and assuming frozen turbulence, L , L and L , nondimensionalized to the WTT ux vx wx length scale, are 0.375, 0.2 and 0.125, respectively. The power spectral density of the velocity components u, v and w at the building height are presented in Fig. 4 (a), (b) and (c), respectively, showing a good agreement with the von Karman (indicated as "Tar. Karman" in the figures) wind spectrum up to the cut-off frequency, equal to approximately 65 Hz. Overall, the simulated turbu- lent inflow shows a satisfactory agreement with WTT and with the targeted profiles. 3 Results and discussion 3.1 Mesh As is well known, it is ideally necessary to ensure that results obtained using numeri- cal simulations are independent from the adopted mesh size [28]. When using LES, it Fig. 3 The turbulence characteristics of the wind fields in WTT and LES: (a) mean wind speed, (b) turbulence intensities and (c) integral length scales Fig. 4 The power spectral density of three velocity components at the building location: (a) u, (b) v and (c) w Xing Advances in Aerodynamics (2023) 5:9 Page 7 of 19 is actually more appropriate to state that results of interest shall be independent or vary within acceptable ranges with the mesh size, as a complete independency is actually not expected. We start by considering the case at 0 . For such case, we use three meshes, namely, coarse mesh (CM), medium mesh (MM) and fine mesh (FM), to check the grid dependency of the numerical solution. The number of cells of the three meshes is 0.1, 1.8 and 6.6 million, respectively. In other words, the ratio between the cell size in each direc- tion between two consecutive meshes is approximately 1.5, following the suggestion in [31]. A summary of the characteristics of the three meshes is provided in Table 1. The computational resources needed to perform each simulation are also reported. In order to characterize the obtained results, the time-average, rms, maximum and ¯ ˆ ˇ minimum values of C , respectively denoted as C , C , C and C , are considered. Peak p p p p values are calculated by fitting to the two-minute extreme using a Gumbel cumulative probability distribution function and then, adopting the well-known shifting property of the Gumbel distribution [32], to extract the 10-minute extremes associated with a non- exceedance probability equal to 80%, in agreement with the Cook and Mayne approach [33]. Figure 5 presents the scatter plots of the C statistics from simulations using the afore- mentioned three meshes. As shown in Fig. 5 (a) - (c), no remarkable differences in terms ¯ ˆ ¯ of C , C and C can be observed, except from a very slight underestimation of C in p p p CM (the difference is approximately 0.1). Regarding C , results obtained from CM only arrive at -2.5. However, they attain approximately -4.0 in the simulations with MM and FM, showing no apparent variation between the two meshes. Some localized points have more severe extreme suctions in MM and FM than in CM, while the results from CM are closer to WTT, which seem nonintuitive. However, as indicated in an analogous work [34], a finer grid model might underestimate the turbulence dissipation rate and give a higher turbulent kinetic energy prediction, yielding overall better performance with coarser meshes. In all, results highlight that CM is probably excessively coarse, as it substantially differs from the other two. Results appear to be quite similar with respect to all analyses quantities for MM and FM, so that MM is chosen for the next analyses. 3.2 Pressure coefficients statistics After choosing the mesh size to be used, we proceed at simulating seven attack angles, ◦ ◦ equally spaced between 0 and 90 . The LES results are compared with WTT results for all the simulated wind attack angles in the form of three error metrics, namely, the mean error (ME), the standard deviation of error (SE) and the coefficient of determi - nation ( R ) [12], as reported below. Table 1 Characteristics of the three meshes used for the present investigation Number of cells Number of cells for each t Computational time y building side (million) H L W (-) (CPU - hour) CM 0.1 56 26 9 0.02 25 110.0 MM 1.8 145 73 23 0.01 215 46.0 FM 6.6 220 113 35 0.005 1212 30.0 Xing Advances in Aerodynamics (2023) 5:9 Page 8 of 19 ¯ ˆ ˇ Fig. 5 Comparisons of LES results using different mesh resolutions with WTT: ( a) C , (b) C , (c) C and (d) C p p p ME = (LES − WTT ), (1) i i i=1 SE = σ(LES − WTT ), (2) (WTT − WTT ) ∗ (LES − LES) R = , (3) σ ∗ σ WTT LES where the subscript i represents the pressure tap number and the total number of moni- tors, N, is equal to 480. The WTT (or LES ) represents the time-average and σ (or WTT σ ) represents the standard variation of the WTT (or LES) dataset. The error met - LES rics regarding different C statistics for all the considered cases are reported in Table 2. Notice that ME for C is not reported, as it might be simply related to the choice of the reference pressure, so being irrelevant. Looking at Table 2, an overall correspondence in terms of C is obtained between WTT and LES, with the average R around 0.99. Similarly, good agreement could also be found regarding the C , whose average R is 0.93. p Xing Advances in Aerodynamics (2023) 5:9 Page 9 of 19 ′ ◦ ◦ ¯ ˆ ˇ Table 2 Error metrics of the C , C , C and C for wind attack angles from 0 to 90 p p p C ˆ ˇ p C C p p p 2 2 2 2 ME SE ME SE ME SE ME SE R R R R 0 - 0.08 0.98 0.00 0.03 0.72 0.00 0.15 0.96 -0.14 0.40 0.75 15 - 0.05 0.99 -0.02 0.03 0.65 -0.01 0.14 0.96 0.15 0.41 0.85 30 - 0.06 0.99 -0.02 0.03 0.86 -0.07 0.14 0.95 0.06 0.18 0.93 45 - 0.02 0.99 -0.03 0.02 0.60 -0.04 0.12 0.96 0.16 0.28 0.68 60 - 0.04 0.99 -0.01 0.03 0.67 0.00 0.11 0.95 0.10 0.29 0.72 75 - 0.05 0.98 -0.03 0.03 0.76 -0.08 0.18 0.86 0.12 0.33 0.75 90 - 0.06 0.98 -0.01 0.03 0.89 -0.06 0.15 0.88 0.11 0.42 0.76 mean - 0.05 0.99 -0.02 0.03 0.74 -0.04 0.15 0.93 0.09 0.33 0.78 Considering all the wind attack angles, we observe that no clear trend emerges for the 2 ′ ◦ error metrics when the angle is varied. For instance, R of C for 0 is 0.72, lower than ◦ ◦ ◦ some skew angles such as 0.86 of 30 and 0.76 of 75 . However, at 90 we obtain the highest value of R among all the cases, equal to 0.89. Substantially, a good matching is obtained when the flow is orthogonal to the short side, leading to a reattached flow. In other conditions, results do not follow simple trends. The prediction of fully detached flow expected at 0 is accurate, but it is not more accurate than other skew angles. Figures 6 and 7 show the local peak values of C considering at the same time wind ◦ ◦ incident angles from 0 to 90 , i.e. the extremes from all attack angles are envel- oped together, as for design purposes. In particular, we show scatter plots reporting in abscissa WTT results and in ordinates LES results. In order to allow to individuate taps belonging to different building faces, pressure taps are grouped as face I , face II , face III and face IV , following the indication in Fig. 1 (a). The results in terms of the peak C appear to be characterized by much higher accuracy with respect to suctions. The extreme negative pressures near the building edges on face II , face III and face IV reach values of -4.0 in WTT but only attain -3.5 in LES. Conversely, LES overestimates the negative extremes near the center bottom of face IV , where the C arrives at -5.0 in LES but reaches only -3.5 in WTT. 3.3 Global forces and moments The global forces for all simulated wind attack angles are reported in dimensionless form C , C , C and C , representing the force and moment components in along-wind Fx Fy Mx My (x) and cross-wind (y) directions, as shown in Figs. 8 and 9. Specifically, the time series of C is integrated over the building surfaces to obtain time histories of the non-dimen- sional force coefficients. Then, for each incidence angle, the peak values (including the maximum and minimum) of forces and moments are calculated following the same pro- cedure used for the pressure peaks explained in Section 3.1. Again, for design purposes, the peaks from all the angles are enveloped and indicated as LES env. and WTT env. in the figures. It can be seen that the LES reproduces the global forces and moments very well with only a slight underestimation in terms of the maximum values of positive C and C . Fx Fy The strongest deviation is on 15 , where the extreme values of both forces and moments Xing Advances in Aerodynamics (2023) 5:9 Page 10 of 19 Fig. 6 The peak values of C from all the degrees: (a) scatter plots, (b) surface distribution of WTT and (c) surface distribution of LES Fig. 7 The peak values of C from all the degrees: (a) scatter plots, (b) surface distribution of WTT and (c) surface distribution of LES in LES are about 25% lower than in WTT. Again, discrepancies can be found for 45 , for which the peak values of forces (and moments) have smaller values in LES than in WTT. Similar results are also found in another work [35], which declared that the wind attack Xing Advances in Aerodynamics (2023) 5:9 Page 11 of 19 Fig. 8 The global forces coefficients of all the simulated wind attack angles: (a) C and (b) C Fx Fy Fig. 9 The moments coefficients of all the simulated wind attack angles: (a) C and (b) C Mx My angle of 15 approaches the critical angle or glancing angle and has the minimum drag force coefficient and maximum mean lift magnitude, and the Strouhal number (St) is maximized. 3.4 Prediction errors for 45 case As the 45 case appears to be the one presenting higher discrepancies between WTT and LES, we further investigate it in this section in more detail. Hereafter the Pre- diction Error (PE) is used as an error metric and is defined as PE = LES − WTT , S S S ¯ ˆ ˇ where the subscript S represents the C statistics, i.e., C , C , C and C . p p p p ′ ◦ Figure 10 presents the scatter plots of C and C for the 45 case. The scatter plots ¯ ˆ of C and C are not reported here, since they are satisfactorily accurate. It can be p p seen that C on face II shows systemic underestimations from the numerical model. In fact, monitors with C approaching 0.5 in WTT are lower than 0.3 in LES. As expected, underestimations also regarding C can be found on these probes, with the negative extremes in LES being 60% less pronounced than in WTT. Besides, some ˇ ˇ localized differences regarding C can be found on face III . The values of C of these p p points are around -3.0 in WTT but are lower than -2.0 in LES. Figures 11 and 12 present the surface distributions of C and C as well as their prediction errors, i.e., PE and PE . More specifically, it can be clearly seen that p C p Xing Advances in Aerodynamics (2023) 5:9 Page 12 of 19 ′ ◦ Fig. 10 Scatter plots of (a) C and (b) C for 45 case ′ ◦ Fig. 11 The surface distributions of the C and the PE ′ for 45 case p C the underestimations mentioned above are located near the leading edge of face II . Looking at Fig. 11 (a) and (b), the higher values of C are also located on this region. As for the independent spots near the top edge of face III where the highest PE appears, one possible explanation can be found in a relative study [23]: the extreme value is recorded downstream of the leading edge, from which strong vortices are expected to be shed. Figure 13 presents the frequency distributions of PE ′ and PE . It can be seen that C ˇ p C the PE ′ concentrates around -0.05. It can also be found that deviations between LES and WTT in terms of extreme suctions ( PE ) are mainly inside the range -1.0 and 0.0, which indicates a tendency of the simulation to underestimate such quantity and provides a measure of the expected underestimations. Xing Advances in Aerodynamics (2023) 5:9 Page 13 of 19 Fig. 12 The surface distributions of the C and the PE for 45 case Fig. 13 The frequency distributions of (a) PE and (b) PE for 45 case C ˇ p C 4 Sensitivity study In this section, the influence of parameters including mesh setting and LES subgrid scale (SGS) model on numerical results for the 45 case is investigated, attempting to ameliorate the previously obtained prediction errors. 4.1 Mesh We firstly consider the possibility to ameliorate results by refining the mesh. In Sec - tion 3, we observed that the selected meshing parameters were able to guarantee a low sensitivity of the obtained results to the cell size. Nevertheless, the study was per- ◦ ◦ formed at 0 , so that it cannot be automatically extended to 45 . We here consider two mesh refinement strategies: in the first one we refine all the fluid volume around the building, while in the second one we add boundary layer cells. Then, the two strategies are combined. Xing Advances in Aerodynamics (2023) 5:9 Page 14 of 19 Fig. 14 The configurations of the mesh near the building: (a) BL, (b) BR and (c) BL_BR Consequently, three adjusted meshes based on MM are obtained. The details of the tested meshes are listed below: (I) boundary layer cells, which have 8 layers and the minimum thickness equaling 0.0004 m and scratching ratio of 1.15, are added to the building surfaces, leading to the mesh named BL; (II) the mesh size is refined to 0.004 m up to a distance of 0.2 m from the building and to 0.002 m to a distance of 0.025 m from the building, resulting in the mesh named BR; (III) the two aforementioned strategies are applied at the same time, yielding the mesh named BL_BR. The mesh zones far away from the building are kept unchanged. An overview of the newly con- sidered meshes in the proximity of the building is provided in Fig. 14 and their char- acteristics are provided in Table 3. It is worth noting that the addition of boundary layer cells dramatically increases (double) the computational costs, though a relatively lower (half ) value of y can be obtained at the building surfaces. The error metrics of the three cases are presented in Table 4. Unexpectedly, the mesh with boundary layer cells, i.e., BL, has worse performance in terms of C predic- 2 ′ tion compared to the original mesh, MM. Specifically, the R of C for mesh MM is 0.60 while for the mesh BL it decreases to 0.54. Differently, the BR mesh shows better ′ 2 performance in terms of C , re aching R equal to 0.65. As it can be seen, despite such variations, the overall error does not change substantially, especially if we consider Table 3 Information of all the tested meshes for wind attack angle equal to 45 Number of cells Number of cells for each t Computational time y building side (million) H L W (-) (CPU - hour) BL 2.1 166 100 33 0.005 563 10.0 BR 6.0 250 150 50 0.003 1670 24.0 BL_BR 7.0 250 150 50 0.001 3600 10.0 ¯ ˆ ˇ Table 4 Error metrics of C , C , C and C for LES with different mesh distributions p p p C ˆ ˇ p C C p p p 2 2 2 2 ME SE ME SE ME SE ME SE R R R R BL - 0.02 0.99 -0.03 0.02 0.54 -0.06 0.11 0.96 0.16 0.27 0.71 BR - 0.02 0.99 -0.03 0.02 0.65 -0.04 0.12 0.96 0.14 0.25 0.76 BL_BR - 0.02 0.99 -0.03 0.02 0.63 -0.04 0.13 0.96 0.08 0.26 0.77 Xing Advances in Aerodynamics (2023) 5:9 Page 15 of 19 the high increase in computational resources needed to make the refinements. We also show that the combination of the two refinement strategies BL_BR does not show a more accurate prediction of C . Now we see the simulated results in terms of C . Again, slightly better performance can be found from BR than from BL, while the results from BR and BL_BR are close to each other. Figures 15 and 16 show the surface distributions of the prediction errors PE ′ and PE for the three meshes. The range of PE for BL is close to the original simula- ˇ C C p tion, MM, still showing high error values on face II , which seems to indicate a detachment of the flow at the edge for WTT not predicted by LES. This can be deduced also by Fig. 11 (a) and (b) in which C is shown. The distribution of PE for Fig. 15 The surface distributions of the PE ′ of different mesh resolutions: (a) BL, (b) BR and (c) BL_BR Fig. 16 The surface distributions of the PE of different mesh distributions: (a) BL, (b) BR and (c) BL_BR p Xing Advances in Aerodynamics (2023) 5:9 Page 16 of 19 all the analysed meshes is shown in Fig. 16, which does not show any major differ- ence between the three. 4.2 Subgrid scale model It has been seen that refining the mesh did not yield strong improvements of the results, despite its high increase of computational costs. We thus here investigate the effect of changing the turbulence model. Table 5 reports error metrics of C statis- tics for simulations on mesh MM with different turbulence models (namely, WALE, k − ω SSTSAS and Smagorinsky). Overall, none of them shows better model per- formance than the original turbulence model, i.e., k-Equation, whose ME for C and C are -0.03 and 0.16, respectively. Only marginal improvements can be seen, for instance, regarding R for C , with the value increasing from 0.68 in the original sim- ulation to 0.71 in the simulation with the Smagorinsky model. From the distributions of PE (Fig. 17) and PE (Fig. 18), the differences between C ˇ the three cases are also insignificant. In a summary, the prediction errors of C and C for 45 are not sensitive to changes in the adopted turbulence model. ¯ ˆ ˇ Table 5 Error metrics of C , C , C and C for different turbulence models p p p C ˆ ˇ p p C C p p 2 2 2 2 ME SE ME SE ME SE ME SE R R R R WALE - 0.02 0.99 -0.03 0.02 0.56 -0.04 0.12 0.96 0.16 0.27 0.70 k − ω SSTSAS - 0.03 0.99 -0.04 0.02 0.38 -0.11 0.12 0.94 0.24 0.26 0.63 Smagorinsky - 0.02 0.99 -0.03 0.02 0.58 -0.05 0.12 0.96 0.17 0.26 0.71 Fig. 17 The surface distributions of the PE ′ of different turbulence models: (a) WALE, (b) k − ω SSTSAS and (c) Smagorinsky Xing Advances in Aerodynamics (2023) 5:9 Page 17 of 19 Fig. 18 The surface distributions of the PE of different turbulence models: (a) WALE, (b) k − ω SSTSAS and Cp (c) Smagorinsky 5 Conclusion Wind fields around an isolated high-rise building model with an aspect ratio of 1:3:5 are simulated using LES to evaluate the model performance and investigate the influ - ences of several model parameters on the prediction errors. Simulations are per- ◦ ◦ ◦ formed for seven wind attack angles ranging from 0 to 90 with an increment of 15 . The main conclusions of the study are hereby drawn as below: 1. Globally, LES is capable of capturing the mean and maximum wind pressure on the building for both perpendicular and skew wind attack angles, with the correlation ¯ ˆ ( R ) between the experimental and predicted results being 0.99 for C and 0.93 for C p p on average. The model performance on the prediction of global forces (moments) is also extremely satisfactory, showing very good agreement with WTT results. How- ever, the numerical accuracy in terms of C and C is less satisfactory, especially in a 2 ′ few cases. The R for C and C are 0.74 and 0.78 on average, respectively; 2. For the considered cases, the worst predictions are obtained with an attack angle of 45 . The high values of PE ′ and PE for this case concentrate on the leading edge in C ˇ p C correspondence of the flow separations. It appears that a flow separation occurs in WTT which is not predicted by the adopted numerical models; ′ ◦ 3. The prediction errors of C and C for 45 are not sensitive to the addition of boundary layer cells for the near-wall treatment. Some marginal improvements are obtained by refining the mesh in the surroundings of the building, but this will dra - matically increase the computational time; 4. No major improvement is obtained by changing the adopted turbulence model. Overall, it is difficult to individuate with certainty the cause of the observed discrep - ancies. Actually, many possible causes can be individuated. Such causes range from small deviations of the geometry with respect to WTT to differences in the incoming Xing Advances in Aerodynamics (2023) 5:9 Page 18 of 19 turbulence, despite the fact that a good matching was obtained in terms of turbulence intensity, and length scales were prescribed according to usual practice. Surely, we observe that the considered case is not particularly sensitive to the adopted turbu- lence model and mesh size, once an appropriate mesh has been initially selected. It must be remarked, that the addition of boundary layers proved, beside not ameliorat- ing results, had only a very limited effect on the obtained results. Author’s contributions Jin Xing: Methodology, investigation, conceptualization, writing, and visualization. The author read and approved the final manuscript. Acknowledgements The author is thankful to CINECA for providing the HPC facilities which allowed the completion of the present study. The author is also thankful to Tokyo Polytechnic University for providing the wind tunnel data. The support provided by the China Scholarship Council (CSC) is acknowledged. Availability of data and materials Data used and/or analysed during the current study is included in this published article and are available from the cor- responding author on reasonable request. Part of the data can be openly accessed at TPU. Funding Not applicable. Declarations Competing interests The author declares that she has no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Received: 16 October 2022 Accepted: 21 January 2023 References 1. Jiang G, Yoshie R (2020) Side ratio effects on flow and pollutant dispersion around an isolated high-rise building in a turbulent boundary layer. Build Environ 180:107078 2. Kataoka H, Ono Y, Enoki K (2020) Applications and prospects of CFD for wind engineering fields. J Wind Eng Ind Aerodyn 205:104310 3. Luo Y, Liu H, Xue H et al (2019) Large-eddy simulation evaluation of wind loads on a high-rise building based on the multiscale synthetic eddy method. Adv Struct Eng 22(4):997–1006 4. Vita G, Salvadori S, Misul DA et al (2020) Eec ff ts of inflow condition on RANS and LES predictions of the flow around a high-rise building. Fluids 5(4):233. Accessed 1 March 2022 5. Tominaga Y, Mochida A, Murakami S et al (2008) Comparison of various revised k–ε models and LES applied to flow around a high-rise building model with 1: 1: 2 shape placed within the surface boundary layer. J Wind Eng Ind Aerodyn 96(4):389–411 6. Zu G, Lam KM (2018) LES and wind tunnel test of flow around two tall buildings in staggered arrangement. Compu- tation 6(2):28 7. Alminhana GW, Braun AL, Loredo-Souza AM (2018) A numerical-experimental investigation on the aerodynamic performance of CAARC building models with geometric modifications. J Wind Eng Ind Aerodyn 180:34–48 8. Liu J, Hui Y, Wang J et al (2021) LES study of windward-face-mounted-ribs’ effects on flow fields and aerodynamic forces on a square cylinder. Build Environ 200:107950 9. Buffa E, Jacob J, Sagaut P (2021) Lattice-Boltzmann-based large-eddy simulation of high-rise building aerodynamics with inlet turbulence reconstruction. J Wind Eng Ind Aerodyn 212:104560 10. Ke J (2019) RANS and hybrid LES/RANS simulations of flow over a square cylinder. Adv Aerodyn 1(1):10 11. Yu Y, Yang Y, Xie Z (2018) A new inflow turbulence generator for large eddy simulation evaluation of wind effects on a standard high-rise building. Build Environ 138:300–313 12. Papp B, Kristóf G, Gromke C (2021) Application and assessment of a GPU-based LES method for predicting dynamic wind loads on buildings. J Wind Eng Ind Aerodyn 217:104739 13. Thordal MS, Bennetsen JC, Capra S et al (2020) Towards a standard CFD setup for wind load assessment of high-rise buildings: Part 1 – Benchmark of the CAARC building. J Wind Eng Ind Aerodyn 205:104283 14. Lamberti G, Gorlé C (2020) Sensitivity of LES predictions of wind loading on a high-rise building to the inflow boundary condition. J Wind Eng Ind Aerodyn 206:104370 15. Thordal MS, Bennetsen JC, Koss HHH (2019) Review for practical application of CFD for the determination of wind load on high-rise buildings. J Wind Eng Ind Aerodyn 186:155–168 Xing Advances in Aerodynamics (2023) 5:9 Page 19 of 19 16. Bruno L, Salvetti MV, Ricciardelli F (2014) Benchmark on the aerodynamics of a rectangular 5:1 cylinder: An overview after the first four years of activity. J Wind Eng Ind Aerodyn 126:87–106 17. Schatzman M (2010) COST 732 model evaluation case studies: approach and results. Meteorological Inst 18. Tominaga Y, Mochida A, Yoshie R et al (2008) AIJ guidelines for practical applications of CFD to pedestrian wind envi- ronment around buildings. J Wind Eng Ind Aerodyn 96(10-11):1749–1761 19. Guichard R (2019) Assessment of an improved random flow generation method to predict unsteady wind pressures on an isolated building using large-eddy simulation. J Wind Eng Ind Aerodyn 189:304–313 20. Wang W, Cao Y, Okaze T (2021) Comparison of hexahedral, tetrahedral and polyhedral cells for reproducing the wind field around an isolated building by LES. Build Environ 195:107717 21. Wijesooriya K, Mohotti D, Chauhan K et al (2019) Numerical investigation of scale resolved turbulence models (LES, ELES and DDES) in the assessment of wind effects on supertall structures. J Build Eng 25:100842 22. Liu J, Niu J, Du Y et al (2019) LES for pedestrian level wind around an idealized building array—assessment of sensi- tivity to influencing parameters. Sustain Cities Soc 44:406–415 23. Cao Y, Liu X, Zhou D et al (2022) Investigation of local severe suction on the side walls of a high-rise building by standard, spectral and conditional POD. Build Environ 217:109047 24. Ricci M, Patruno L, Kalkman I et al (2018) Towards LES as a design tool: wind loads assessment on a high-rise build- ing. J Wind Eng Ind Aerodyn 180:1–18 25. Wang Y, Chen X (2020) Simulation of approaching boundary layer flow and wind loads on high-rise buildings by wall-modeled LES. J Wind Eng Ind Aerodyn 207:104410 26. Tamura Y, Ohkuma T, Kawai H, Uematsu Y, Kondo K (2004) Revision of AIJ recommendations for wind loads on build- ings. Structures 2004: Building on the Past, Securing the Future 2004:1–10 27. Tokyo Polytechnic University (2007) TPU aerodynamic database. http:// wind. arch.t- kougei. ac. jp/ system/ eng/ conte nts/ code/ tpu. Accessed 1 March 2022 28. Franke J, Hellsten A, Schlünzen H et al (2007) Best practice guideline for the CFD simulation of flows in the urban environment-a summary. COST Office, Belgium 29. Patruno L, Ricci M (2017) On the generation of synthetic divergence-free homogeneous anisotropic turbulence. Comput Methods Appl Mech Eng 315:396–417 30. Bervida M, Patruno L, Stanič S et al (2020) Synthetic generation of the atmospheric boundary layer for wind loading assessment using spectral methods. J Wind Eng Ind Aerodyn 196:104040 31. Ferziger JH, Perić M (2002) Computational methods for fluid dynamics, 3rd edn. Springer, Berlin, Heidelberg 32. Dyrbye C, Hansen SO (1997) Wind loads on structures. Wiley, Chichester 33. Cook NJ, Mayne JR (1980) A refined working approach to the assessment of wind loads for equivalent static design. J Wind Eng Ind Aerodyn 6(1–2):125–137 34. Gousseau P, Blocken B, van Heijst GJF (2013) Quality assessment of large-eddy simulation of wind flow around a high-rise building: validation and solution verification. Comput Fluids 79:120–133 35. Cao Y, Tamura T, Kawai H (2019) Investigation of wall pressures and surface flow patterns on a wall-mounted square cylinder using very high-resolution Cartesian mesh. J Wind Eng Ind Aerodyn 188:1–18 Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Journal
Advances in Aerodynamics – Springer Journals
Published: Mar 2, 2023
Keywords: Wind loading; CFD; High-rise buildings; Large eddy simulation; Aerodynamics