Jeffrey D. Mirocha

Development of a Coupled Groundwater–Atmosphere Model

Abstract Complete models of the hydrologic cycle have gained recent attention as research has shown interdependence between the coupled land and energy balance of the subsurface, land surface, and lower atmosphere. PF.WRF is a new model that is a combination of the Weather Research and Forecasting (WRF) atmospheric model and a parallel hydrology model (ParFlow) that fully integrates three-dimensional, variably saturated subsurface flow with overland flow. These models are coupled in an explicit, operator-splitting manner via the Noah land surface model (LSM). Here, the coupled model formulation and equations are presented and a balance of water between the subsurface, land surface, and atmosphere is verified. The improvement in important physical processes afforded by the coupled model using a number of semi-idealized simulations over the Little Washita watershed in the southern Great Plains is demonstrated. These simulations are initialized with a set of offline spinups to achieve a balanced state of ini...

Implementation of a Nonlinear Subfilter Turbulence Stress Model for Large-Eddy Simulation in the Advanced Research WRF Model

Abstract Two formulations of a nonlinear turbulence subfilter-scale (SFS) stress model were implemented into the Advanced Research Weather Research and Forecasting model (ARW-WRF) version 3.0 for improved large-eddy simulation performance. The new models were evaluated against the WRF model’s standard Smagorinsky and 1.5-order turbulence kinetic energy (TKE) linear eddy-viscosity SFS stress models in simulations of geostrophically forced, neutral boundary layer flow over both flat terrain and a shallow, symmetric transverse ridge. Comparisons of simulation results with similarity profiles indicate that the nonlinear models significantly improve agreement with the expected profiles near the surface, reducing the overprediction of near-surface stress characteristic of linear eddy-viscosity models with no near-wall damping. Comparisons of simulations conducted using different mesh sizes indicate that the nonlinear model simulations at coarser resolutions agree more closely with the higher-resolution results ...

Implementation of a generalized actuator disk wind turbine model into the weather research and forecasting model for large-eddy simulation applications

A generalized actuator disk (GAD) wind turbine parameterization designed for large-eddy simulation (LES) applications was implemented into the Weather Research and Forecasting (WRF) model. WRF-LES with the GAD model enables numerical investigation of the effects of an operating wind turbine on and interactions with a broad range of atmospheric boundary layer phenomena. Numerical simulations using WRF-LES with the GAD model were compared with measurements obtained from the Turbine Wake and Inflow Characterization Study (TWICS-2011), the goal of which was to measure both the inflow to and wake from a 2.3-MW wind turbine. Data from a meteorological tower and two light-detection and ranging (lidar) systems, one vertically profiling and another operated over a variety of scanning modes, were utilized to obtain forcing for the simulations, and to evaluate characteristics of the simulated wakes. Simulations produced wakes with physically consistent rotation and velocity deficits. Two surface heat flux values of ...

Large eddy simulation of wind turbine wake dynamics in the stable boundary layer using the Weather Research and Forecasting Model

Recently, an actuator disk parameterization was implemented in the Weather Research and Forecasting (WRF) Model for large eddy simulation (LES) of wind turbine wakes. To thoroughly verify this model, simulations of various types of turbines and atmospheric conditions must be evaluated against corresponding experimental data. In this work, numerical simulations are compared to nacelle-based scanning lidar measurements taken in stable atmospheric conditions during a field campaign conducted at a wind farm in the western United States. Using several wake characteristics—such as the velocity deficit, centerline location, and wake width—as metrics for model verification, the simulations show good agreement with the observations. Notable results include a high average velocity deficit, decreasing from 73% at a downwind distance x of 1.2 rotor diameters (D) to 25% at x = 6.6D, resulting from a low average wind speed and therefore high average turbine thrust coefficient. Moreover, the wake width expands from 1.4D...

Interaction of Nocturnal Low-Level Jets with Urban Geometries as seen in Joint URBAN 2003 Data

Abstract Because accurate modeling of atmospheric flows in urban environments requires sophisticated representation of complex urban geometries, much work has been devoted to treatment of the urban surface. However, the importance of the larger-scale flow impinging upon the urban complex to the flow, transport, and dispersion within it and downwind has received less attention. Building-resolving computational fluid dynamics (CFD) models are commonly employed to investigate interactions between the flow and three-dimensional structures that make up the urban environment; however, such models are typically forced with simplified boundary conditions that fail to include important regional-scale phenomena that can strongly influence the flow within the urban complex and downwind. This paper investigates the interaction of an important and frequently occurring regional-scale phenomenon, the nocturnal low-level jet (LLJ), with urban-scale turbulence and dispersion in Oklahoma City, Oklahoma, using data from the...

A stochastic perturbation method to generate inflow turbulence in large-eddy simulation models: Application to neutrally stratified atmospheric boundary layers

Despite the variety of existing methods, efficient generation of turbulent inflow conditions for large-eddy simulation (LES) models remains a challenging and active research area. Herein, we extend our previous research on the cell perturbation method, which uses a novel stochastic approach based upon finite amplitude perturbations of the potential temperature field applied within a region near the inflow boundaries of the LES domain [Munoz-Esparza et al., “Bridging the transition from mesoscale to microscale turbulence in numerical weather prediction models,” Boundary-Layer Meteorol., 153, 409–440 (2014)]. The objective was twofold: (i) to identify the governing parameters of the method and their optimum values and (ii) to generalize the results over a broad range of atmospheric large-scale forcing conditions, Ug = 5 − 25 m s−1, where Ug is the geostrophic wind. We identified the perturbation Eckert number, Ec=Ug2/ρcpθpm, to be the parameter governing the flow transition to turbulence in neutrally strat...

Investigating wind turbine impacts on near-wake flow using profiling lidar data and large-eddy simulations with an actuator disk model

Wind turbine impacts on the atmospheric flow are investigated using data from the Crop Wind Energy Experiment (CWEX-11) and large-eddy simulations (LESs) utilizing a generalized actuator disk (GAD) wind turbine model. CWEX-11 employed velocity-azimuth display (VAD) data from two Doppler lidar systems to sample vertical profiles of flow parameters across the rotor depth both upstream and in the wake of an operating 1.5 MW wind turbine. Lidar and surface observations obtained during four days of July 2011 are analyzed to characterize the turbine impacts on wind speed and flow variability, and to examine the sensitivity of these changes to atmospheric stability. Significant velocity deficits ( VD) are observed at the downstream location during both convective and stable portions of four diurnal cycles, with large, sustained deficits occurring during stable conditions. Variances of the streamwise velocity component, σu, likewise show large increases downstream during both stable and unstable conditions, with ...

The Wind Forecast Improvement Project (WFIP): A Public–Private Partnership Addressing Wind Energy Forecast Needs

AbstractThe Wind Forecast Improvement Project (WFIP) is a public–private research program, the goal of which is to improve the accuracy of short-term (0–6 h) wind power forecasts for the wind energy industry. WFIP was sponsored by the U.S. Department of Energy (DOE), with partners that included the National Oceanic and Atmospheric Administration (NOAA), private forecasting companies (WindLogics and AWS Truepower), DOE national laboratories, grid operators, and universities. WFIP employed two avenues for improving wind power forecasts: first, through the collection of special observations to be assimilated into forecast models and, second, by upgrading NWP forecast models and ensembles. The new observations were collected during concurrent year-long field campaigns in two high wind energy resource areas of the United States (the upper Great Plains and Texas) and included 12 wind profiling radars, 12 sodars, several lidars and surface flux stations, 184 instrumented tall towers, and over 400 nacelle anemome...

Comparison of Measured and Numerically Simulated Turbulence Statistics in a Convective Boundary Layer Over Complex Terrain

The Weather Research and Forecasting (WRF) model can be used to simulate atmospheric processes ranging from quasi-global to tens of m in scale. Here we employ large-eddy simulation (LES) using the WRF model, with the LES-domain nested within a mesoscale WRF model domain with grid spacing decreasing from 12.15 km (mesoscale) to 0.03 km (LES). We simulate real-world conditions in the convective planetary boundary layer over an area of complex terrain. The WRF-LES model results are evaluated against observations collected during the US Department of Energy-supported Columbia Basin Wind Energy Study. Comparison of the first- and second-order moments, turbulence spectrum, and probability density function of wind speed shows good agreement between the simulations and observations. One key result is to demonstrate that a systematic methodology needs to be applied to select the grid spacing and refinement ratio used between domains, to avoid having a grid resolution that falls in the grey zone and to minimize artefacts in the WRF-LES model solutions. Furthermore, the WRF-LES model variables show large variability in space and time caused by the complex topography in the LES domain. Analyses of WRF-LES model results show that the flow structures, such as roll vortices and convective cells, vary depending on both the location and time of day as well as the distance from the inflow boundaries.

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