Development and Evaluation of a Reynolds-Averaged Navier-Stokes Solver in WindNinja for Operational Wildland Fire Applications

An open source computational fluid dynamics (CFD) solver has been incorporated into the WindNinja modeling framework widely used by wildland fire managers as well as researchers and practitioners in other fields, such as wind energy, wind erosion, and search and rescue. Here we describe incorporation of the CFD solver and evaluate its performance compared to the conservation of mass (COM) solver in WindNinja and previously published large-eddy simulations (LES) for three field campaigns conducted over isolated terrain obstacles of varying terrain complexity: Askervein Hill, Bolund Hill, and Big Southern Butte. We also compare the effects of two important model settings in the CFD solver and provide guidance on model sensitivity to these settings. Additionally, we investigate the computational mesh and difficulties regarding terrain representation. Two important findings from this work are: (1) the choice of discretization scheme for advection has a significantly larger effect on the simulated winds than the choice of turbulence model and (2) CFD solver predictions are significantly better than the COM solver predictions at windward and lee side observation locations, but no difference was found in predicted speed-up at ridgetop locations between the two solvers.

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

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.

Validation of atmospheric boundary layer turbulence model by on-site measurements

Modeling atmosperic boundary layer with standard linear models does not sufficiently reproduce wind conditions in complex terrain, especially on leeward sides of terrain slopes. More complex models, based on Reynolds averaged Navier-Stokes equations and two-equation k-? turbulence models for neutral conditions in atmospheric boundary layer, written in general curvilinear non-orthogonal co-ordinate system, have been evaluated. In order to quantify the differences and level of accuracy of different turbulence models, investigation has been performed using standard k-? model without additional production terms and k-? turbulence models with modified set of model coefficients. The sets of full conservation equations are numerically solved by computational fluid dynamics technique. Numerical calculations of turbulence models are compared to the reference experimental data of Askervein hill measurements.

Evaluation of Turbulence Closure Models for Large-Eddy Simulation over Complex Terrain: Flow over Askervein Hill

The evaluation of turbulence closure models for large-eddy simulation (LES) has primarily been performed over flat terrain, where comparisons with theory and observations are simplified. The authors have previously developed improved closure models using explicit filtering and reconstruction, together with a dynamic eddy-viscosity model and a near-wall stress term. This dynamic reconstruction model (DRM) is a mixed model, combining scale-similarity and eddy-viscosity components. The DRM gave improved results over standard eddy-viscosity models for neutral boundary layer flow over flat but rough terrain, yielding the expected logarithmic velocity profiles near the wall. The results from the studies over flat terrain are now extended to flow over full-scale topography. The test case is flow over Askervein Hill, an isolated hill in western Scotland, where a field campaign was conducted in 1983 with the purpose of capturing wind data representing atmospheric episodes under near-neutral stratification and steady wind conditions. This widely studied flow provides a more challenging test case for the new turbulence models because of the sloping terrain and separation in the lee of the hill. Since an LES formulation is used, a number of simulation features are different than those typically used in the Askervein literature. The simulations are inherently unsteady, the inflow conditions are provided by a separate turbulent flow database, and (uniquely herein) ensemble averages of the turbulent flow results are used in comparisons with field data. Results indicate that the DRM can improve the predictions of flow speedup and especially turbulent kinetic energy (TKE) over the hill when compared with the standard TKE-1.5 model. This is the first study, to the authors’ knowledge, in which explicit filtering and reconstruction (scale similarity) and dynamic turbulence models have been applied to full-scale simulations of the atmospheric boundary layer over terrain. Simulations with the lowest level of reconstruction are straightforward. Increased levels of reconstruction, however, present difficulties when used with a dynamic eddy-viscosity model. An alternative mixed model is proposed to avoid the complexities associated with the dynamic procedure and to allow higher levels of reconstruction; this mixed model combines a standard TKE-1.5 eddy-viscosity closure with velocity reconstruction to form a simple and efficient turbulence model that gives good results for both mean flow and turbulence over Askervein Hill. The results indicate that significant improvements in LES over complex terrain can be obtained by the use of mixed models that combine scale-similarity and eddy-viscosity components.