Adaptation of Pressure Based CFD Solvers for Mesoscale Atmospheric Problems

General purpose Computational Fluid Dynamics (CFD) solvers are frequently used in small-scale urban pollution dispersion simulations without a large extent of ver- tical flow. Vertical flow, however, plays an important role in the formation of local breezes, such as urban heat island induced breezes that have great significance in the ventilation of large cities. The effects of atmospheric stratification, anelasticity and Coriolis force must be taken into account in such simulations. We introduce a general method for adapting pressure based CFD solvers to atmospheric flow simulations in order to take advantage of their high flexibility in geometrical modelling and meshing. Compressibility and thermal stratification effects are taken into account by utilizing a novel system of transformations of the field variables and by adding consequential source terms to the model equations of incompressible flow. Phenomena involving mesoscale to microscale coupled effects can be analyzed without model nesting, applying only local grid refinement of an arbitrary level. Elements of the method are validated against an analytical solution, results of a reference calculation, and a laboratory scale urban heat island circulation experiment. The new approach can be applied with benefits to several areas of application. Inclusion of the moisture transport phenomena and the surface energy balance are important further steps towards the practical application of the method.

[1]  János Unger,et al.  Intra-urban relationship between surface geometry and urban heat island: review and new approach , 2004 .

[2]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[3]  D. Reinert,et al.  A new large-eddy simulation model for simulating air flow and warm clouds above highly complex terrain. Part I: The dry model , 2007 .

[4]  C. Montavon,et al.  Validation of a non-hydrostatic numerical model to simulate stratified wind fields over complex topography , 1998 .

[5]  T. Shih,et al.  A New K-epsilon Eddy Viscosity Model for High Reynolds Number Turbulent Flows: Model Development and Validation , 1994 .

[6]  William H. Snyder,et al.  A Laboratory Study of the Urban Heat Island in a Calm and Stably Stratified Environment. Part II: Velocity Field , 1997 .

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

[8]  Antonio Cenedese,et al.  Interaction between an Inland Urban Heat Island and a Sea-Breeze Flow: A Laboratory Study , 2003 .

[9]  John Paul Iselin,et al.  A comparison of contaminant plume statistics from a Gaussian puff and urban CFD model for two large cities , 2005 .

[10]  Hiroshi Yoshik Ado.  Numerical study of the daytime urban effect and its interaction with the sea breeze , 1992 .

[11]  Albert F. Kurbatskii,et al.  Computational Modeling of the Turbulent Penetrative Convection above the Urban Heat Island in a Stably Stratified Environment , 2001 .

[12]  Louis J. Wicker,et al.  Numerical solutions of a non‐linear density current: A benchmark solution and comparisons , 1993 .