Improvement of k-epsilon turbulence model for CFD simulation of atmospheric boundary layer around a high-rise building using stochastic optimization and Monte Carlo sampling technique

The accuracy of the computational fluid dynamics (CFD) to model the airflow around the buildings in the atmospheric boundary layer (ABL) is directly linked to the utilized turbulence model. Despite the popularity and their low computational cost, the current Reynolds Averaged Navier-Stokes (RANS) models cannot accurately resolve the wake regions behind the buildings. The default values of the RANS models’ closure coefficients in CFD tools such as ANSYS CFX, ANSYS FLUENT, PHOENIX, and STAR CCM+ are mainly adapted from other fields and physical problems, which are not perfectly suitable for ABL flow modeling. This study embarks on proposing a systematic approach to find the optimum values for the closure coefficients of RANS models in order to significantly improve the accuracy of CFD simulations for urban studies. The methodology is based on stochastic optimization and Monte Carlo Sampling technique. To show the capability of the method, a test case of airflow around an isolated building placed in a non-isothermal unstable ABL was considered. The recommended values for this case study in accordance with the optimization method were thus found to be 1.45 ≤

[1]  S. Murakami,et al.  COMPARISON OF VARIOUS TURBULENCE MODELS APPLIED TO A BLUFF BODY , 1993 .

[2]  Abdelkader Frendi,et al.  Uncertainty Quantification of Turbulence Model Coefficients via Latin Hypercube Sampling Method , 2010 .

[3]  Alex Van der Velden,et al.  Isight Design Optimization Methodologies , 2009 .

[4]  W. Rodi,et al.  Calculation of the flow past a surface-mounted cube with two-layer turbulence models , 1997 .

[5]  P. Kastner-Kleina,et al.  A wind tunnel study of organised and turbulent air motions in urban street canyons , 2001 .

[6]  K. Hibi,et al.  Turbulent measurments of the flow field around a high-rise building , 1998 .

[7]  Kenny C. S Kwok,et al.  Wind tunnel study of pedestrian level wind environment around tall buildings: Effects of building dimensions, separation and podium , 2012 .

[8]  T. Shih,et al.  A new k-ϵ eddy viscosity model for high reynolds number turbulent flows , 1995 .

[9]  P. Richards,et al.  Appropriate boundary conditions for computational wind engineering models using the k-ε turbulence model , 1993 .

[10]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[11]  P. Spalart.  A One-Equation Turbulence Model for Aerodynamic Flows , 1992 .

[12]  John C. LaRue,et al.  The decay power law in grid-generated turbulence , 1990, Journal of Fluid Mechanics.

[13]  J. Hunt,et al.  Kinematical studies of the flows around free or surface-mounted obstacles; applying topology to flow visualization , 1978, Journal of Fluid Mechanics.

[14]  Yoshihiko Hayashi,et al.  Examining the κ-ϵ model by means of a wind tunnel test and large-eddy simulation of the turbulence structure around a cube , 1990 .

[15]  I. Tezaur,et al.  Uncertainty Quantification , 2011, Encyclopedia of Parallel Computing.

[16]  Parham A. Mirzaei,et al.  Pollution removal effectiveness of the pedestrian ventilation system , 2011 .

[17]  W. Rodi.  Comparison of LES and RANS calculations of the flow around bluff bodies , 1997 .

[18]  Stavros Tavoularis,et al.  Further experiments on the evolution of turbulent stresses and scales in uniformly sheared turbulence , 1989, Journal of Fluid Mechanics.

[19]  Yoshihide Tominaga,et al.  Cooperative project for CFD prediction of pedestrian wind environment in the Architectural Institute of Japan , 2007 .

[20]  Anita Schöbel,et al.  Algorithm Engineering in Robust Optimization , 2016, Algorithm Engineering.

[21]  Jacques Periaux,et al.  Uncertainty based robust optimization method for drag minimization problems in aerodynamics , 2012 .

[22]  Peter Richards,et al.  Appropriate boundary conditions for computational wind engineering models revisited , 2011 .

[23]  P. Duynkerke.  Application of the E – ε Turbulence Closure Model to the Neutral and Stable Atmospheric Boundary Layer , 1988 .

[24]  Yoshihide Tominaga,et al.  Development of a new k–ε model to reproduce the aerodynamic effects of snow particles on a flow field , 2015 .

[25]  Ryuichiro Yoshie,et al.  CFD simulations of gas dispersion around high-rise building in non-isothermal boundary layer , 2011 .

[26]  Yun Kyu Yi,et al.  Dynamic integration between building energy simulation (BES) and computational fluid dynamics (CFD) simulation for building exterior surface , 2013 .

[27]  S. Murakami,et al.  Development of a new k−ε model for flow and pressure fields around bluff body , 1997 .

[28]  S. Orszag,et al.  Development of turbulence models for shear flows by a double expansion technique , 1992 .

[29]  Peter Richards,et al.  Pedestrian level wind speeds in downtown Auckland , 2002 .

[30]  Shuzo Murakami,et al.  Environmental design of outdoor climate based on CFD , 2006 .

[31]  Michael Schatzmann,et al.  Numerical Model Inter-comparison for Wind Flow and Turbulence Around Single-Block Buildings , 2011 .

[32]  Serge Guillas,et al.  STATISTICAL CALIBRATION OF CFD MODELLING FOR STREET CANYON FLOWS , 2011 .

[33]  Bje Bert Blocken,et al.  Quality assessment of Large-Eddy Simulation of wind flow around a high-rise building : validation and solution verification , 2013 .

[34]  H. Panofsky,et al.  Atmospheric Turbulence: Models and Methods for Engineering Applications , 1984 .

[35]  Parham A. Mirzaei,et al.  A novel approach to enhance outdoor air quality: Pedestrian ventilation system , 2010 .

[36]  Megdouda Tari,et al.  Refined descriptive sampling: A better approach to Monte Carlo simulation , 2006, Simul. Model. Pract. Theory.

[37]  Jonas Allegrini,et al.  Coupled CFD, radiation and building energy model for studying heat fluxes in an urban environment with generic building configurations , 2015 .

[38]  Nyuk Hien Wong,et al.  Integrated urban microclimate assessment method as a sustainable urban development and urban design tool , 2011 .

[39]  Slawomir Koziel,et al.  Multi-Fidelity Robust Aerodynamic Design Optimization under Mixed Uncertainty , 2015 .

[40]  Hanqing Wu.  Pedestrian-level wind environment around buildings , 1994 .

[41]  Erik Dick,et al.  Prediction of the pressure distribution on a cubical building with implicit LES , 2010 .

[42]  C. Mann,et al.  A Practical Treatise on Diseases of the Skin , 1889, Atlanta Medical and Surgical Journal (1884).

[43]  Mathias W. Rotach,et al.  A wind tunnel study of organised and turbulent air motions in urban street canyons , 2001 .

[44]  Sai Hung Cheung,et al.  Bayesian uncertainty analysis with applications to turbulence modeling , 2011, Reliab. Eng. Syst. Saf..

[45]  Parham A. Mirzaei,et al.  A procedure to quantify the impact of mitigation techniques on the urban ventilation , 2012 .

[46]  B. Blocken,et al.  CFD simulation of pollutant dispersion around isolated buildings: on the role of convective and turbulent mass fluxes in the prediction accuracy. , 2011, Journal of hazardous materials.

[47]  Parham A. Mirzaei,et al.  Impact of non-uniform urban surface temperature on pollution dispersion in urban areas , 2011 .

[48]  Takeshi Ohkuma,et al.  Numerical prediction of wind loading on buildings and structures : Activities of AIJ cooperative project on CFD , 1997 .

[49]  Koichi Hishida,et al.  Buoyancy effects on the wake behind a heated obstacle immersed in a turbulent boundary layer , 1995 .

[50]  Parham A. Mirzaei,et al.  Dynamical computational fluid dynamics modeling of the stochastic wind for application of urban studies , 2013 .

[51]  R N Meroney,et al.  On the use of numerical modelling for near-field pollutant dispersion in urban environments--A review. , 2016, Environmental pollution.

[52]  Sergio Teggi,et al.  Analysis of the urban heat island effects on building energy consumption , 2015 .

[53]  P. Moin,et al.  Turbulence statistics in fully developed channel flow at low Reynolds number , 1987, Journal of Fluid Mechanics.

[54]  Serge Guillas,et al.  Bayesian calibration of the constants of the k–ε turbulence model for a CFD model of street canyon flow , 2014 .

[55]  D. Etling,et al.  Application of the E-ε turbulence model to the atmospheric boundary layer , 1985 .

[56]  Y. Tominaga,et al.  Numerical simulation of dispersion around an isolated cubic building: Model evaluation of RANS and LES , 2010 .

[57]  Hester Bijl,et al.  Bayesian estimates of parameter variability in the k-ε turbulence model , 2014, J. Comput. Phys..

[58]  P. A. Mirzaei,et al.  Recent challenges in modeling of urban heat island , 2015 .

[59]  Edna Shaviv,et al.  Climatic aspects in urban design—a case study , 2003 .

[60]  Ren-Jye Yang,et al.  Design for six sigma through robust optimization , 2004 .

[61]  S. Orszag,et al.  Renormalization group analysis of turbulence. I. Basic theory , 1986, Physical review letters.

[62]  Todd A. Oliver,et al.  Bayesian uncertainty quantification applied to RANS turbulence models , 2011 .

[63]  B. Launder,et al.  The numerical computation of turbulent flows , 1990 .

[64]  Abdelkader Frendi,et al.  Uncertainty Quantification of Turbulence Model Coefficients via Latin Hypercube Sampling Method , 2011 .

[65]  Yoshihide Tominaga,et al.  AIJ guidelines for practical applications of CFD to pedestrian wind environment around buildings , 2008 .

[66]  Y. Tominaga.  Flow around a high-rise building using steady and unsteady RANS CFD: Effect of large-scale fluctuations on the velocity statistics , 2015 .

[67]  Parham A. Mirzaei,et al.  Toward design and fabrication of wind-driven vehicles: Procedure to optimize the threshold of driving forces , 2013 .

[68]  M. Kato.  The modeling of turbulent flow around stationary and vibrating square cylinders , 1993 .

[69]  K Schittkowski.  NLPQLP: A Fortran Implementation of a Sequential Quadratic Programming Algorithm with Distributed and Non-Monotone Line Search , 2005 .

[70]  Jonas Allegrini.  Coupled CFD, Radiation and Building Energy Model for Studying the Impact of Building Height Topology and Buoyancy on Local Heat Island Formation in Urban Environments , 2017 .

[71]  Yoshihide Tominaga,et al.  Cross Comparisons of CFD Results of Wind Environment at Pedestrian Level around a High-rise Building and within a Building Complex , 2004 .