An Immersed Boundary Method for the Weather Research and Forecasting Model

2010 ◽  
Vol 138 (3) ◽  
pp. 796-817 ◽  
Author(s):  
Katherine A. Lundquist ◽  
Fotini Katopodes Chow ◽  
Julie K. Lundquist
Keyword(s):  
Immersed Boundary Method ◽  
Complex Terrain ◽  
Wrf Model ◽  
Neumann Boundary ◽  
Weather Research And Forecasting ◽  
Immersed Boundary ◽  
Mesoscale Models ◽  
Boundary Method ◽  
Terrain Following ◽  
Weather Research

Abstract This paper describes an immersed boundary method that facilitates the explicit resolution of complex terrain within the Weather Research and Forecasting (WRF) model. Mesoscale models, such as WRF, are increasingly used for high-resolution simulations, particularly in complex terrain, but errors associated with terrain-following coordinates degrade the accuracy of the solution. The use of an alternative-gridding technique, known as an immersed boundary method, alleviates coordinate transformation errors and eliminates restrictions on terrain slope that currently limit mesoscale models to slowly varying terrain. Simulations are presented for canonical cases with shallow terrain slopes, and comparisons between simulations with the native terrain-following coordinates and those using the immersed boundary method show excellent agreement. Validation cases demonstrate the ability of the immersed boundary method to handle both Dirichlet and Neumann boundary conditions. Additionally, realistic surface forcing can be provided at the immersed boundary by atmospheric physics parameterizations, which are modified to include the effects of the immersed terrain. Using the immersed boundary method, the WRF model is capable of simulating highly complex terrain, as demonstrated by a simulation of flow over an urban skyline.

scholarly journals Mesoscale to Microscale Simulations over Complex Terrain with the Immersed Boundary Method in the Weather Research and Forecasting Model

2020 ◽  
Vol 148 (2) ◽  
pp. 577-595 ◽  
Author(s):  
David J. Wiersema ◽  
Katherine A. Lundquist ◽  
Fotini Katopodes Chow
Keyword(s):  
Boundary Conditions ◽  
Immersed Boundary Method ◽  
Complex Terrain ◽  
Wrf Model ◽  
Weather Research And Forecasting ◽  
Immersed Boundary ◽  
Eddy Simulation ◽  
Boundary Method ◽  
Large Eddy ◽  
Terrain Following

Abstract Improvements to the Weather Research and Forecasting (WRF) Model are made to enable multiscale simulations over highly complex terrain with dynamically downscaled boundary conditions from the mesoscale to the microscale. Over steep terrain, the WRF Model develops numerical errors that are due to grid deformation of the terrain-following coordinates. An alternative coordinate system, the immersed boundary method (IBM), has been implemented into WRF, allowing for simulations over highly complex terrain; however, the new coordinate system precluded nesting within mesoscale simulations using WRF’s native terrain-following coordinates. Here, the immersed boundary method and WRF’s grid-nesting framework are modified to seamlessly work together. This improved framework for the first time allows for large-eddy simulation over complex (urban) terrain with IBM to be nested within a typical mesoscale WRF simulation. Simulations of the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma, are performed using a multiscale five-domain nested configuration, spanning horizontal grid resolutions from 6 km to 2 m. These are compared with microscale-only simulations with idealized lateral boundary conditions and with observations of wind speed/direction and SF6 concentrations from a controlled release from intensive observation period 3. The multiscale simulation, which is configured independent of local observations, shows similar model skill predicting wind speed/direction and improved skill predicting SF6 concentrations when compared with the idealized simulations, which require use of observations to set mean flow conditions. Use of this improved multiscale framework shows promise for enabling large-eddy simulation over highly complex terrain with dynamically downscaled boundary conditions from mesoscale models.

Large-Eddy Simulation over Complex Terrain Using an Improved Immersed Boundary Method in the Weather Research and Forecasting Model

2018 ◽  
Vol 146 (9) ◽  
pp. 2781-2797 ◽  
Author(s):  
Jingyi Bao ◽  
Fotini Katopodes Chow ◽  
Katherine A. Lundquist
Keyword(s):  
Boundary Condition ◽  
Immersed Boundary Method ◽  
Complex Terrain ◽  
Weather Research And Forecasting ◽  
Immersed Boundary ◽  
Flat Terrain ◽  
Boundary Method ◽  
Terrain Following ◽  
Bolund Hill ◽  
Askervein Hill

Abstract The Weather Research and Forecasting (WRF) Model is increasingly being used for higher-resolution atmospheric simulations over complex terrain. With increased resolution, resolved terrain slopes become steeper, and the native terrain-following coordinates used in WRF result in numerical errors and instability. The immersed boundary method (IBM) uses a nonconformal grid with the terrain surface represented through interpolated forcing terms. Lundquist et al.’s WRF-IBM implementation eliminates the limitations of WRF’s terrain-following coordinate and was previously validated with a no-slip boundary condition for urban simulations and idealized terrain. This paper describes the implementation of a log-law boundary condition into WRF-IBM to extend its applicability to general atmospheric complex terrain simulations. The implementation of the improved WRF-IBM boundary condition is validated for neutral flow over flat terrain and the complex terrain cases of Askervein Hill, Scotland, and Bolund Hill, Denmark. First, comparisons are made to similarity theory and standard WRF results for the flat terrain case. Then, simulations of flow over the moderately sloped Askervein Hill are used to demonstrate agreement between the IBM and terrain-following WRF results, as well as agreement with observations. Finally, Bolund Hill simulations show that WRF-IBM can handle steep topography (standard WRF fails) and compares well to observations. Overall, the new WRF-IBM boundary condition shows improved performance, though the leeside representation of the flow can be potentially further improved.

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

2012 ◽  
Vol 140 (12) ◽  
pp. 3936-3955 ◽  
Author(s):  
Katherine A. Lundquist ◽  
Fotini Katopodes Chow ◽  
Julie K. Lundquist
Keyword(s):  
Immersed Boundary Method ◽  
Complex Terrain ◽  
Sonic Anemometer ◽  
Three Dimensional ◽  
Wrf Model ◽  
Oklahoma City ◽  
Immersed Boundary ◽  
Field Campaign ◽  
Boundary Method ◽  
Outer Domain

Abstract This 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 City is simulated concurrently with an outer domain with flat terrain using one-way nesting to generate a turbulent flow field. Results from the IBM-WRF simulation of IOP3 compare well with observations from the field campaign, as well as with results from an urban computational fluid dynamics code, Finite Element Model in 3-Dimensions and Massively Parallelized (FEM3MP), which used body-fitted coordinates. Using the FAC2 performance metric from Chang and Hanna, which is the fraction of predictions within a factor of 2 of observations, IBM-WRF achieves 100% and 71% for velocity predictions using cup and sonic anemometer observations, respectively. For the passive scalar, 53% of the model predictions meet the FAC5 (factor of 5) criteria.

Topographic Effects on Radiation in the WRF Model with the Immersed Boundary Method: Implementation, Validation, and Application to Complex Terrain

2018 ◽  
Vol 146 (10) ◽  
pp. 3277-3292 ◽  
Author(s):  
Robert S. Arthur ◽  
Katherine A. Lundquist ◽  
Jeffrey D. Mirocha ◽  
Fotini K. Chow
Keyword(s):  
Immersed Boundary Method ◽  
Complex Terrain ◽  
Sensible Heat Flux ◽  
Wrf Model ◽  
Atmospheric Modeling ◽  
Surface Energy Budget ◽  
Immersed Boundary ◽  
Topographic Effects ◽  
Grid Spacing ◽  
Boundary Method

Abstract Topographic effects on radiation, including both topographic shading and slope effects, are included in the Weather Research and Forecasting (WRF) Model, and here they are made compatible with the immersed boundary method (IBM). IBM is an alternative method for representing complex terrain that reduces numerical errors over sloped terrain, thus extending the range of slopes that can be represented in WRF simulations. The implementation of topographic effects on radiation is validated by comparing land surface fluxes, as well as temperature and velocity fields, between idealized WRF simulations both with and without IBM. Following validation, the topographic shading implementation is tested in a semirealistic simulation of flow over Granite Mountain, Utah, where topographic shading is known to affect downslope flow development in the evening. The horizontal grid spacing is 50 m and the vertical grid spacing is approximately 8–27 m near the surface. Such a case would fail to run in WRF with its native terrain-following coordinates because of large local slope values reaching up to 55°. Good agreement is found between modeled surface energy budget components and observations from the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) program at a location on the east slope of Granite Mountain. In addition, the model captures large spatiotemporal inhomogeneities in the surface sensible heat flux that are important for the development of thermally driven flows over complex terrain.

An Immersed Boundary Method for Simulation of Wind Flow Over Complex Terrain

2011 ◽  
Vol 134 (1) ◽  
Author(s):  
S. Jafari ◽  
N. Chokani ◽  
R. S. Abhari
Keyword(s):  
Immersed Boundary Method ◽  
Complex Terrain ◽  
Stokes Equations ◽  
Ground Level ◽  
Rectangular Domain ◽  
Immersed Boundary ◽  
Test Case ◽  
Wind Flow ◽  
Wall Functions ◽  
Boundary Method

The accurate modeling of the wind resource over complex terrain is required to optimize the micrositing of wind turbines. In this paper, an immersed boundary method that is used in connection with the Reynolds-averaged Navier–Stokes equations with k-ω turbulence model in order to efficiently simulate the wind flow over complex terrain is presented. With the immersed boundary method, only one Cartesian grid is required to simulate the wind flow for all wind directions, with only the rotation of the digital elevation map. Thus, the lengthy procedure of generating multiple grids for conventional rectangular domain is avoided. Wall functions are employed with the immersed boundary method in order to relax the stringent near-wall grid resolution requirements as well as to allow the effects of surface roughness to be accounted for. The immersed boundary method is applied to the complex terrain test case of Bolund Hill. The simulation results of wind speed and turbulent kinetic energy show good agreement with experiments for heights greater than 5 m above ground level.

An Evaluation of a Hybrid, Terrain-Following Vertical Coordinate in the WRF-Based RAP and HRRR Models

2020 ◽  
Vol 35 (3) ◽  
pp. 1081-1096 ◽  
Author(s):  
Jeffrey Beck ◽  
John Brown ◽  
Jimy Dudhia ◽  
David Gill ◽  
Tracy Hertneky ◽  
...  
Keyword(s):  
Complex Terrain ◽  
Wrf Model ◽  
Vertical Motion ◽  
Weather Research And Forecasting ◽  
Mountainous Terrain ◽  
Vertical Coordinate ◽  
Numerical Noise ◽  
Terrain Following ◽  
Model Equations ◽  
Real Time Simulations

Abstract A new hybrid, sigma-pressure vertical coordinate was recently added to the Weather Research and Forecasting (WRF) Model in an effort to reduce numerical noise in the model equations near complex terrain. Testing of this hybrid, terrain-following coordinate was undertaken in the WRF-based Rapid Refresh (RAP) and High-Resolution Rapid Refresh (HRRR) models to assess impacts on retrospective and real-time simulations. Initial cold-start simulations indicated that the majority of differences between the hybrid and traditional sigma coordinate were confined to regions downstream of mountainous terrain and focused in the upper levels. Week-long retrospective simulations generally resulted in small improvements for the RAP, and a neutral impact in the HRRR when the hybrid coordinate was used. However, one possibility is that the inclusion of data assimilation in the experiments may have minimized differences between the vertical coordinates. Finally, analysis of turbulence forecasts with the new hybrid coordinate indicate a significant reduction in spurious vertical motion over the full length of the Rocky Mountains. Overall, the results indicate a potential to improve forecast metrics through implementation of the hybrid coordinate, particularly at upper levels, and downstream of complex terrain.

Sign in / Sign up

Export Citation Format

Share Document