skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: High-resolution simulations of downslope flows over complex terrain using WRF-IBM

; ; ; ;
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
OSTI Identifier:
Report Number(s):
DOE Contract Number:
Resource Type:
Resource Relation:
Conference: Presented at: 17th Conference on Mountain Meteorology, Burlington, VT, United States, Jun 27 - Jul 01, 2016
Country of Publication:
United States

Citation Formats

Arthur, R S, Lundquist, K A, Mirocha, J D, Hoch, S W, and Chow, F K. High-resolution simulations of downslope flows over complex terrain using WRF-IBM. United States: N. p., 2016. Web.
Arthur, R S, Lundquist, K A, Mirocha, J D, Hoch, S W, & Chow, F K. High-resolution simulations of downslope flows over complex terrain using WRF-IBM. United States.
Arthur, R S, Lundquist, K A, Mirocha, J D, Hoch, S W, and Chow, F K. Thu . "High-resolution simulations of downslope flows over complex terrain using WRF-IBM". United States. doi:.
title = {High-resolution simulations of downslope flows over complex terrain using WRF-IBM},
author = {Arthur, R S and Lundquist, K A and Mirocha, J D and Hoch, S W and Chow, F K},
abstractNote = {},
doi = {},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Thu Sep 29 00:00:00 EDT 2016},
month = {Thu Sep 29 00:00:00 EDT 2016}

Other availability
Please see Document Availability for additional information on obtaining the full-text document. Library patrons may search WorldCat to identify libraries that hold this conference proceeding.

Save / Share:
  • Presentation given at the 98th AMS Annual meeting on January 8, 2018. Information presented includes the motivation, model configuration, observation info, statistics, summary and conclusion.
  • Boundary layer flows are greatly complicated by the presence of complex terrain which redirects mean flow and alters the structure of turbulence. Surface fluxes of heat and moisture provide additional forcing which induce secondary flows, or can dominate flow dynamics in cases with weak mean flows. Mesoscale models are increasingly being used for numerical simulations of boundary layer flows over complex terrain. These models typically use a terrain-following coordinate transformation, but these introduce numerical errors over steep terrain. An alternative is to use an immersed boundary method which alleviates errors associated with the coordinate transformation by allowing the terrain tomore » be represented as a surface which arbitrarily passes through a Cartesian grid. This paper describes coupling atmospheric physics models to an immersed boundary method implemented in the Weather Research and Forecasting (WRF) model in previous work [Lundquist et al., 2007]. When the immersed boundary method is used, boundary conditions must be imposed on the immersed surface for velocity and scalar surface fluxes. Previous algorithms, such as those used by Tseng and Ferziger [2003] and Balaras [2004], impose no-slip boundary conditions on the velocity field at the immersed surface by adding a body force to the Navier-Stokes equations. Flux boundary conditions for the advection-diffusion equation have not been adequately addressed. A new algorithm is developed here which allows scalar surface fluxes to be imposed on the flow solution at an immersed boundary. With this extension of the immersed boundary method, land-surface models can be coupled to the immersed boundary to provide realistic surface forcing. Validation is provided in the context of idealized valley simulations with both specified and parameterized surface fluxes using the WRF code. Applicability to real terrain is illustrated with a fully coupled two-dimensional simulation of the Owens Valley in California.« less
  • The WRF model version 3.3 is used to simulate near hub-height winds and power ramps utilizing three commonly used planetary boundary-layer (PBL) schemes: Mellor-Yamada-Janjic (MYJ), University of Washington (UW), and Yonsei University (YSU). The predicted winds have small mean biases compared with observations. Power ramps and step changes (changes within an hour) consistently show that the UW scheme performed better in predicting up ramps under stable conditions with higher prediction accuracy and capture rates. Both YSU and UW scheme show good performance predicting up- and down- ramps under unstable conditions with YSU being slightly better for ramp durations longer thanmore » an hour. MYJ is the most successful simulating down-ramps under stable conditions. The high wind speed and large shear associated with low-level jets are frequently associated with power ramps, and the biases in predicted low-level jet explain some of the shown differences in ramp predictions among different PBL schemes. Low-level jets were observed as low as ~200 m in altitude over the Columbia Basin Wind Energy Study (CBWES) site, located in an area of complex terrain. The shear, low-level peak wind speeds, as well as the height of maximum wind speed are not well predicted. Model simulations with 3 PBL schemes show the largest variability among them under stable conditions.« less
  • A finite difference computer code (MATHEW) for ensuring mass-consistency in wind fields, developed by Sherman (1978), has proven successful in providing input for ADPIC (Lange, 1978), a particle-in-cell advection-diffusion code used in pollution transport studies. An analog of MATHEW is described to be used in finite element advection-diffusion codes, finite element planetary boundary models, and for wind energy assessment studies in complex terrain.