An Immersed Boundary Method Enabling Large-Eddy Simulations of Flow over Complex Terrain in the WRF Model

AbstractThis paper describes a three-dimensional immersed boundary method (IBM) that facilitates the explicit resolution of complex terrain within the Weather Research and Forecasting (WRF) model. Two interpolation methods—trilinear and inverse distance weighting (IDW)—are used at the core of the IBM algorithm. This work expands on the previous two-dimensional IBM algorithm of Lundquist et al., which uses bilinear interpolation. Simulations of flow over a three-dimensional hill are performed with WRF’s native terrain-following coordinate and with both IB methods. Comparisons of flow fields from the three simulations show excellent agreement, indicating that both IB methods produce accurate results. IDW proves more adept at handling highly complex urban terrain, where the trilinear interpolation algorithm fails. This capability is demonstrated by using the IDW core to model flow in Oklahoma City, Oklahoma, from intensive observation period 3 (IOP3) of the Joint Urban 2003 field campaign. Flow in Oklahoma C...

[1]  R. Franke.  Scattered data interpolation: tests of some methods , 1982 .

[2]  Y. Mahrer,et al.  An Improved Numerical Approximation of the Horizontal Gradients in a Terrain-Following Coordinate System , 1984 .

[3]  A. Adcroft,et al.  Representation of Topography by Shaved Cells in a Height Coordinate Ocean Model , 1997 .

[4]  P. M. Gresho,et al.  Projection 2 Goes Turbulent - and Fully Implicit , 1998 .

[5]  J. Ferziger,et al.  A ghost-cell immersed boundary method for flow in complex geometry , 2002 .

[6]  D. Lüthi,et al.  A new terrain-following vertical coordinate formulation for atmospheric prediction models , 2002 .

[7]  Günther Zängl,et al.  An Improved Method for Computing Horizontal Diffusion in a Sigma-Coordinate Model and Its Application to Simulations over Mountainous Topography , 2002 .

[8]  G. Iaccarino,et al.  Immersed boundary technique for turbulent flow simulations , 2003 .

[9]  Oliver Fuhrer,et al.  Numerical consistency of metric terms in terrain-following coordinates , 2003 .

[10]  A Generalized Sigma-Coordinate System for the MM5 , 2003 .

[11]  Olav R. Hansen,et al.  FLACS CFD air quality model performance evaluation with Kit Fox, MUST, Prairie Grass, and EMU observations , 2004 .

[12]  S. Hanna,et al.  Air quality model performance evaluation , 2004 .

[13]  E. Balaras.  Modeling complex boundaries using an external force field on fixed Cartesian grids in large-eddy simulations , 2004 .

[14]  G. Zängl,et al.  Numerical errors above steep topography: A model intercomparison , 2004 .

[15]  Robert L. Lee,et al.  Flow around a Complex Building: Comparisons between Experiments and a Reynolds-Averaged Navier¿Stokes Approach , 2004 .

[16]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

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

[18]  Michael Manhart,et al.  High‐order stable interpolations for immersed boundary methods , 2006 .

[19]  M. Leach,et al.  A Validation of FEM3MP with Joint Urban 2003 Data , 2006 .

[20]  K. Allwine,et al.  Joint Urban 2003: Study Overview And Instrument Locations , 2006 .

[21]  J. Lundquist,et al.  Consequences of Urban Stability Conditions for Computational Fluid Dynamics Simulations of Urban Dispersion , 2007 .

[22]  Y. Tseng,et al.  An improved hybrid Cartesian/immersed boundary method for fluid–solid flows , 2007 .

[23]  E. Hendricks,et al.  Modeling Turbulent Flow in an Urban Central Business District , 2007 .

[24]  Rajat Mittal,et al.  A sharp interface immersed boundary method for compressible viscous flows , 2007, J. Comput. Phys..

[25]  J. Heagy,et al.  Comparisons of Transport and Dispersion Model Predictions of the Joint Urban 2003 Field Experiment , 2008 .

[26]  Roni Avissar,et al.  The Ocean-Land-Atmosphere Model (OLAM). Part II: Formulation and Tests of the Nonhydrostatic Dynamic Core , 2008 .

[27]  B. Stevens,et al.  Monte Carlo Spectral Integration: a Consistent Approximation for Radiative Transfer in Large Eddy Simulations , 2008 .

[28]  Rajat Mittal,et al.  A versatile sharp interface immersed boundary method for incompressible flows with complex boundaries , 2008, J. Comput. Phys..

[29]  Jean-Christophe Golaz,et al.  One-Way Nested Large-Eddy Simulation over the Askervein Hill , 2009 .

[30]  K. Lundquist,et al.  Immersed Boundary Methods for High-Resolution Simulation of Atmospheric Boundary-Layer Flow Over Complex Terrain , 2010 .

[31]  Tamir G. Reisin,et al.  Evaluation of the atmospheric RAMS model in an obstacle resolving configuration , 2010 .

[32]  J. Lundquist,et al.  An Immersed Boundary Method for the Weather Research and Forecasting Model , 2014 .

[33]  Petra M. Klein,et al.  Concentration fluctuations in a downtown urban area. Part I: analysis of Joint Urban 2003 full-scale fast-response measurements , 2011 .

[34]  Inanc Senocak,et al.  A CFD-based wind solver for an urban fast response transport and dispersion model , 2011 .

[35]  Rebecca Vernot,et al.  Comparisons of JU2003 observations with four diagnostic urban wind flow and Lagrangian particle dispersion models , 2011 .

[36]  Sarah-Jane Lock,et al.  Demonstration of a Cut-Cell Representation of 3D Orography for Studies of Atmospheric Flows over Very Steep Hills , 2012 .