Numerical modelling of the turbulent flow developing within and over a 3-d building array, part ii: a mathematical foundation for a distributed drag force approach

In this paper, we lay the foundations of a systematic mathematical formulation for the governing equations for flow through an urban canopy (e.g., coarse-scaled building array) where the effects of the unresolved obstacles on the flow are represented through a distributed mean-momentum sink. This, in turn, implies additional corresponding terms in the transport equations for the turbulence quantities. More specifically, a modified k--ε model is derived for the simulation of the mean wind speed and turbulence for a neutrally stratified flow through and over a building array, where groups of buildings in the array are aggregated and treated as a porous medium. This model is based on time averaging the spatially averaged Navier--Stokes equations, in which the effects of the obstacle--atmosphere interaction are included through the introduction of a volumetric momentum sink (representing drag on the unresolved buildings in the array).The k--ε turbulence closure model requires two additional prognostic equations, namely one for the time-averaged resolved-scale kinetic energy of turbulence,κ, and another for the dissipation rate, ε, of κ . The transport equation for κ is derived directly from the transport equation for the spatially averaged velocity, and explicitly includes additional sources and sinks that arise from time averaging the product of the spatially averaged velocity fluctuations and the distributed drag force fluctuations. We show how these additional source/sink terms in the transport equation for κ can be obtained in a self-consistent manner from a parameterization of the sink term in the spatially averaged momentum equation. Towards this objective, the time-averaged product of the spatially averaged velocity fluctuations and the distributed drag force fluctuations can be approximated systematically using a Taylor series expansion. A high-order approximation is derived to represent this source/sink term in the transport equation for κ . The dissipation rate (ε) equation is simply obtained as a dimensionally consistent analogue of the κ equation. The relationship between the proposed mathematical formulation of the equations for turbulent flow within an urban canopy (where the latter is treated as a porous medium) and an earlier heuristic two-band spectral decomposition for parameterizing turbulence in a plant canopy is explored in detail.

[1]  Fue-Sang Lien,et al.  Simulation of mean flow and turbulence over a 2D building array using high-resolution CFD and a distributed drag force approach , 2004 .

[2]  Francis H. Harlow,et al.  Transport Equations in Turbulence , 1970 .

[3]  Adrian E. Scheidegger,et al.  The physics of flow through porous media , 1957 .

[4]  D. Decroix.  Large-eddy and Gaussian simulations of downwind dispersion from large building HVAC exhaust , 2000 .

[5]  W. L. Ferrar,et al.  Theory and Problems of Real Variables , 1971 .

[6]  A. Thom.  The exchange of momentum, mass, and heat between an artificial leaf and the airflow in a wind‐tunnel , 1968 .

[7]  Michael J. Brown,et al.  COMPARISON OF CENTERLINE VELOCITY MEASUREMENTS OBTAINED AROUND 2D AND 3D BUILDING ARRAYS IN A WIND TUNNEL , 2001 .

[8]  K. B. Oldham,et al.  An Atlas of Functions. , 1988 .

[9]  Fue-Sang Lien,et al.  Numerical Modelling of the Turbulent Flow Developing Within and Over a 3-D Building Array, Part I: A High-Resolution Reynolds-Averaged Navier—Stokes Approach , 2004 .

[10]  J. Finnigan,et al.  A first‐order closure for diturbed plant‐canopy flows, and its application to winds in a canopy on a ridge , 1998 .

[11]  Gerhard Ecker,et al.  Quantum Chromodynamics , 2006 .

[12]  John D. Wilson,et al.  A second-order closure model for flow through vegetation , 1988 .

[13]  A. Chamberlain.  Transport of gases to and from surfaces with bluff and wave‐like roughness elements , 1968 .

[14]  E. Inoue,et al.  On the Turbulent Structure of Airflow within , 1963 .

[15]  John D. Wilson,et al.  Numerical studies of flow through a windbreak , 1985 .

[16]  S. Whitaker,et al.  The spatial averaging theorem revisited , 1985 .

[17]  S. Belcher,et al.  Adjustment of a turbulent boundary layer to a canopy of roughness elements , 2003, Journal of Fluid Mechanics.

[18]  J. L. Lage,et al.  A modified form of the κ–ε model for turbulent flows of an incompressible fluid in porous media , 2000 .

[19]  W. S. Smith.  A CFD Model Intercomparison and Validation using High Resolution Wind Tunnel Data , 2000 .

[20]  G.P.A. Bot,et al.  Wind-induced airflow through permeable materials, Part I: the motion equation , 2001 .

[21]  E. Takle,et al.  A numerical simulation of boundary-layer flows near shelterbelts , 1995 .

[22]  J. Finnigan,et al.  A Second-Order Closure for Neutrally Stratified Vegetative Canopy Flows , 1999 .

[23]  M. Raupach,et al.  Averaging procedures for flow within vegetation canopies , 1982 .

[24]  T. A. Black,et al.  E-ε modelling of turbulent air flow downwind of a model forest edge , 1996 .

[25]  Hagen Kleinert,et al.  Gauge fields in condensed matter , 1989 .

[26]  Fue-Sang Lien,et al.  Numerical modelling of the turbulent flow developing within and over a 3-d building array, part iii: a istributed drag force approach, its implementation and application , 2005 .

[27]  E. Jaynes.  On the rationale of maximum-entropy methods , 1982, Proceedings of the IEEE.

[28]  Z. Uchijima,et al.  The Energy Budget at the Earth's Surface: An Experimental Study of Air Flow in a Corn Plant-Air Layer. , 1963 .

[29]  Eugene S. Takle,et al.  Boundary-layer flow and turbulence near porous obstacles , 1995 .

[30]  W. Davidon,et al.  Mathematical Methods of Physics , 1965 .

[31]  S. Green.  Modelling Turbulent Air Flow in a Stand of Widely-Spaced Trees , 1992 .

[32]  J. Finnigan,et al.  Turbulent Transport in Flexible Plant Canopies , 1985 .

[33]  Eugene Yee,et al.  Calculation of winds disturbed by an array of fences , 2003 .

[34]  P. Moin,et al.  The basic equations for the large eddy simulation of turbulent flows in complex geometry , 1995 .

[35]  John D. Wilson,et al.  Comments on ’a Numerical Simulation of Boundary-Layer Flows Near Shelterbelts‘ by H. Wang and E. Takle , 1997 .

[36]  Matthias Roth,et al.  Review of atmospheric turbulence over cities , 2007 .

[37]  I. R. Cowan.  Mass, heat and momentum exchange between stands of plants and their atmospheric environment , 1968 .

[38]  P. Moin,et al.  A General Class of Commutative Filters for LES in Complex Geometries , 1998 .

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

[40]  Roland Schiestel,et al.  Multiple‐time‐scale modeling of turbulent flows in one‐point closures , 1987 .

[41]  N. Wilson,et al.  A Higher Order Closure Model for Canopy Flow , 1977 .

[42]  Brian Launder,et al.  Multiple-Time-Scale Concepts in Turbulent Transport Modeling , 1980 .

[43]  M. Raupach,et al.  Experiments on scalar dispersion within a model plant canopy part I: The turbulence structure , 1986 .

[44]  Christophe Sanz.  A Note on k - ε Modelling of Vegetation Canopy Air-Flows , 2003 .

[45]  Adrian E. Scheidegger,et al.  The Physics of Flow Through Porous Media (3rd Edition) , 1957 .

[46]  Brian Launder,et al.  On the computation of buoyancy-driven turbulent flows in rectangular enclosures , 1989 .

[47]  Peter A. Coppin,et al.  Experiments on scalar dispersion within a model plant canopy part II: An elevated plane source , 1986 .