Katherine A. Lundquist

An Immersed Boundary Method for the Weather Research and Forecasting Model

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.

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...

Intercomparison of Mesoscale Model Simulations of the Daytime Valley Wind System

AbstractThree-dimensional simulations of the daytime thermally induced valley wind system for an idealized valley–plain configuration, obtained from nine nonhydrostatic mesoscale models, are compared with special emphasis on the evolution of the along-valley wind. The models use the same initial and lateral boundary conditions, and standard parameterizations for turbulence, radiation, and land surface processes. The evolution of the mean along-valley wind (averaged over the valley cross section) is similar for all models, except for a time shift between individual models of up to 2 h and slight differences in the speed of the evolution. The analysis suggests that these differences are primarily due to differences in the simulated surface energy balance such as the dependence of the sensible heat flux on surface wind speed. Additional sensitivity experiments indicate that the evolution of the mean along-valley flow is largely independent of the choice of the dynamical core and of the turbulence parameteriz...

A New Vertical Grid Nesting Capability in the Weather Research and Forecasting (WRF) Model

AbstractMesoscale atmospheric models are increasingly used for high-resolution (<3 km) simulations to better resolve smaller-scale flow details. Increased resolution is achieved using mesh refinement via grid nesting, a procedure where multiple computational domains are integrated either concurrently or in series. A constraint in the concurrent nesting framework offered by the Weather Research and Forecasting (WRF) Model is that mesh refinement is restricted to the horizontal dimensions. This limitation prevents control of the grid aspect ratio, leading to numerical errors due to poor grid quality and preventing grid optimization. Herein, a procedure permitting vertical nesting for one-way concurrent simulation is developed and validated through idealized cases. The benefits of vertical nesting are demonstrated using both mesoscale and large-eddy simulations (LES). Mesoscale simulations of the Terrain-Induced Rotor Experiment (T-REX) show that vertical grid nesting can alleviate numerical errors due to la...

Implementation of the Immersed Boundary Method in the Weather Research and Forecasting model

Accurate simulations of atmospheric boundary layer flow are vital for predicting dispersion of contaminant releases, particularly in densely populated urban regions where first responders must react within minutes and the consequences of forecast errors are potentially disastrous. Current mesoscale models do not account for urban effects, and conversely urban scale models do not account for mesoscale weather features or atmospheric physics. The ultimate goal of this research is to develop and implement an immersed boundary method (IBM) along with a surface roughness parameterization into the mesoscale Weather Research and Forecasting (WRF) model. IBM will be used in WRF to represent the complex boundary conditions imposed by urban landscapes, while still including forcing from regional weather patterns and atmospheric physics. This document details preliminary results of this research, including the details of three distinct implementations of the immersed boundary method. Results for the three methods are presented for the case of a rotation influenced neutral atmospheric boundary layer over flat terrain.

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

AbstractTopographic 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 w...

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

AbstractThe Weather Research and Forecasting (WRF) Model is increasingly being used for higher-resolution atmospheric simulations over complex terrain. With increased resolution, resolved terrain s...

Assessment of Vertical Mesh Refinement in Concurrently Nested Large-Eddy Simulations Using the Weather Research and Forecasting Model

AbstractTo facilitate multiscale simulation using the Weather Research and Forecasting Model, vertical mesh refinement for one-way concurrent nested simulation was recently introduced. Grid refinement in the vertical dimension removes the requirement of different grid aspect ratios on the bounding versus the nested domain, such that results from refinement are in the horizontal directions only, and thereby can also reduce numerical errors on the bounding domain arising from large aspect ratios in the presence of complex terrain. Herein, the impacts of vertical grid refinement on the evolving downstream flow in nested large-eddy simulations are evaluated in relation to other model configuration choices, including turbulence subfilter-scale (SFS) stress models, mesh configuration, and alternative methods for calculating several near-surface flow parameters. Although vertical nesting requires coarsening of the vertical grid on the bounding domain, leading to a smaller range of resolved turbulence scales in t...