SEARCH · Engineering Papers

Results for “WRF simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2026

An Optimized Parameterization of Sub‐Grid Scale Advection for Convection Permitting Models

Convection‐permitting models (CPMs) explicitly resolve deep convection yet under‐resolve the organized lateral exchanges among drafts and their environment that control entrainment/detrainment, precipitation efficiency, and mesoscale structure. In this work, we introduce the Optimized Advection Scheme (OAS), which introduces a small rotation of the Cartesian frame of reference for the horizontal winds relative to other variables used in advection that induces cross‐gradient transport to mimic under‐resolved convective mixing. The rotation angle is selected to minimize the Kullback–Leibler divergence between the simulated and satellite observed precipitation intensity distributions, yielding a physically consistent perturbation that is computationally inexpensive and portable. Optimized Advection Scheme is implemented in WRF and evaluated over Amazon (April 2014). It shifts precipitation–precipitable‐water joint distributions toward lighter rain, reduces overly intense rates, and improves mesoscale convective system (MCS) lifetime and propagation. Mechanistically, the added cross‐gradient transport promotes convective detrainment and environmental mixing, which cools and moistens the mid‐troposphere, weakens downward momentum transport, alleviates excessive downwelling shortwave biases, and warms the surface temperature. The optimized rotation angle yields comparable improvements at 4‐km and 1‐km grid spacing, demonstrating resolution‐independent benefits across the CPM gray zone. By targeting the dynamical root of under‐mixed convective circulations, rather than tuning model microphysics or closures, OAS delivers robust, scale‐aware improvements in precipitation statistics, cloud vertical structure, and characteristics of MCS (MCSs), offering a practical pathway to more reliable CPM simulations for weather and climate applications.

CPM↗

DOE OSTI2021

Challenging and Improving the Simulation of Mid-Level Mixed-Phase Clouds Over the High-Latitude Southern Ocean

Climate models exhibit major radiative biases over the Southern Ocean owing to a poor representation of mixed-phase clouds. In this study, we used the remote-sensing dataset from the Measurements of Aerosols, Radiation and Clouds over the Southern Ocean (MARCUS) campaign to assess the ability of the Weather Research and Forecasting (WRF) model to reproduce frontal clouds off Antarctica. It focuses on the modeling of thin mid-level supercooled liquid water layers which precipitate ice. The standard version of WRF produces almost fully glaciated clouds and cannot reproduce cloud top turbulence. Our work demonstrates the importance of adapting the ice nucleation parameterization to the pristine austral atmosphere to reproduce the supercooled liquid layers. Once simulated, droplets significantly impact the cloud radiative effect by increasing downwelling longwave fluxes and decreasing downwelling shortwave fluxes at the surface. The net radiative effect is a warming of snow and ice covered surfaces and a cooling of the ocean. Despite improvements in our simulations, the local turbulent circulation related to cloud-top radiative cooling is not properly reproduced, advocating for the need to develop a parameterization for top-down convection to capture the turbulence-microphysics interplay at cloud top.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Sensitivity of deep convection and cross-tropopause water transport to microphysical parameterizations in WRF

Deep convective storms can overshoot the tropopause, thus altering the composition of the stratosphere by vertically transporting tropospheric air. The transport of water vapor and ice particles into a sub-saturated environment can hydrate the stratosphere, with implications for radiative forcing and ozone chemistry. Cloud-resolved models, if employed at high spatial resolutions, are used to probe process-level questions about cross-tropopause deep convective hydration and its controls. There is considerable diversity in model representations of processes associated with water transport and transformation, and the choice of a microphysics scheme affects model skill in simulating deep convective events. This motivates our evaluation of state-of-the-art, as well as widely used standard schemes, in a high spatial- and temporal-resolution framework. Six bulk microphysics schemes were employed in a WRF-LES setup, initialized with a sounding profile representative of a tropopause-overshooting storm. Here, we used an idealized framework to isolate the effect of microphysics on the dominant processes that control the reach of deep convection and stratospheric hydration. All schemes produced the highest reaching updrafts 8-12 hours into the simulation but the strength and persistence of updrafts varied across the schemes; maximum storm heights ranged 9.1-12.6 km across the schemes. Varying microphysics produced large differences in the vertical extent and horizontal aggregation of convection, and an order of magnitude spread in above-tropopause water vapor concentrations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Effect of data assimilation in the Polar WRF with 3DVAR on the prediction of radiation, heat flux, cloud, and near surface atmospheric variables over Svalbard

Although the Polar Weather Research and Forecasting (PWRF) model has been developed for polar environments, simulations of atmospheric states in the polar region still have large uncertainties. Therefore, the effects of data assimilation (DA) to improve forecasts in the polar region were investigated for September 2017 using PWRF and the three-dimensional variational (3DVAR) DA method. The experiments without DA and those with DA assimilating only conventional observations and both conventional observations and satellite radiance data were performed. The forecasts from all experiments both without and with DA underestimated (overestimated) the downward longwave (shortwave) radiation flux due to the underestimation of the amount of Arctic clouds. When satellite radiance data (i.e., the Advanced Microwave Sounding Unit-A (AMSU-A) and Microwave Humidity Sounder (MHS)) were assimilated in addition to conventional observations in the PWRF, the distribution and amount of water vapor became closer to observations, which improves cloud liquid water forecasts. Therefore, when both conventional observations and radiance data were assimilated, the 25–30 h forecast errors of the downward longwave and shortwave radiation fluxes and sensible and latent heat fluxes decreased by 12.7%, 8.1%, 3.2%, and 7.8% with the WSM5 scheme and by 17.1%, 4.7%, 2.5%, and 3.1% with the 2-moment Morrison scheme, respectively, compared to those in the experiments without DA. The forecast errors of the 10 m wind and 2 m temperature with DA were smaller than those without DA at most observation stations. Therefore, the uncertainties of the Arctic forecasts in the PWRF decreased when using DA. To further improve Arctic forecasts, the assimilation of various additional satellite data is necessary.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Simulated wind farm wake sensitivity to configuration choices in the Weather Research and Forecasting model version 3.8.1

Wakes from wind farms can extend over 50 km downwind in stably stratified conditions. These wakes can undermine power production at downwind turbines, adversely impacting revenue. As such, wind farm wake impacts must be considered in wind resource assessments, especially in regions of dense wind farm development. The open-source Weather Research and Forecasting (WRF) numerical weather prediction model includes a wind farm parameterization to estimate wind farm wake effects, but model configuration choices can influence the resulting predictions of wind farm wakes. These choices include vertical resolution, horizontal resolution, and whether or not to include the addition of turbulent kinetic energy generated by the rotating wind turbines. Despite the sensitivity to model configuration, no clear guidance currently exists for these options. Here we compare simulated wind farm wakes produced by varying model configurations with meteorological observations near a land-based wind farm in flat terrain over several diurnal cycles. A WRF configuration comprised of horizontal resolutions of 3 km or 1 km paired with a vertical resolution of 10 m provides the most accurate representation of wind farm wake effects, such as the correct surface warming and elevated wind speed deficit. The inclusion of turbine-generated turbulence is also critical to produce accurate surface warming and should not be omitted.

17 WIND ENERGY↗

DOE OSTI2022

Improved Representation of Horizontal Variability and Turbulence in Mesoscale Simulations of an Extended Cold-Air Pool Event

Abstract Cold-air pools (CAPs), or stable atmospheric boundary layers that form within topographic basins, are associated with poor air quality, hazardous weather, and low wind energy output. Accurate prediction of CAP dynamics presents a challenge for mesoscale forecast models in part because CAPs occur in regions of complex terrain, where traditional turbulence parameterizations may not be appropriate. This study examines the effects of the planetary boundary layer (PBL) scheme and horizontal diffusion treatment on CAP prediction in the Weather Research and Forecasting (WRF) Model. Model runs with a one-dimensional (1D) PBL scheme and Smagorinsky-like horizontal diffusion are compared with runs that use a new three-dimensional (3D) PBL scheme to calculate turbulent fluxes. Simulations are completed in a nested configuration with 3-km/750-m horizontal grid spacing over a 10-day case study in the Columbia River basin, and results are compared with observations from the Second Wind Forecast Improvement Project. Using event-averaged error metrics, potential temperature and wind speed errors are shown to decrease both with increased horizontal grid resolution and with improved treatment of horizontal diffusion over steep terrain. The 3D PBL scheme further reduces errors relative to a standard 1D PBL approach. Error reduction is accentuated during CAP erosion, when turbulent mixing plays a more dominant role in the dynamics. Last, the 3D PBL scheme is shown to reduce near-surface overestimates of turbulence kinetic energy during the CAP event. The sensitivity of turbulence predictions to the master length-scale formulation in the 3D PBL parameterization is also explored. Significance Statement In this article, we demonstrate how a new framework for modeling atmospheric turbulence improves cold pool predictions, using a case study from January 2017 in the Columbia River basin (U.S. Pacific Northwest). Cold pools are regions of cold, stagnant air that form within valleys or basins, and improved forecasts could help to mitigate the risks they pose to air quality, transportation, and wind energy production. For the chosen case study, our tests show a reduction in temperature and wind speed errors by up to a factor of 2–3 relative to standard model options. These results strongly motivate continued development of the framework as well as its application to other complex weather events.

17 WIND ENERGY↗

DOE OSTI2025

Improved Representations of Land‐Atmosphere Interactions Over the Continental U.S. Through Dynamic Root Modeling

Previous studies have identified the oversimplified root system representation as a key factor leading to inaccuracies in vegetation‐atmosphere feedbacks. In this study, a dynamic root water uptake scheme in the Noah‐MP land surface model has been coupled to the Weather Research and Forecasting (WRF) model to investigate its impact on the surface hydroclimate variables and land‐atmosphere interactions. To evaluate the impact of the dynamic root, two coupled simulations were conducted, one with the dynamic root water uptake scheme (DynRt) and one with the static root water uptake scheme (StcRt), which is based on the default root representation in Noah‐MP, with slight modifications, primarily in vegetation‐related parameters. Both DynRt and StcRt simulations were conducted with a small ensemble of three members to account for variations in physical parameterizations, initial and boundary forcing and model setup. When compared with reference data sets, the DynRt simulations show improved results than the StcRt simulations, reducing biases in the simulated leaf area index, surface energy fluxes, soil moisture and precipitation. Two different mechanisms through which roots affect land‐atmosphere coupling have been identified. Over the transitional climate zone between the dry and wet climate, the dynamic root scheme affects surface climate and land‐atmosphere coupling mainly through changes in soil moisture through hydraulic redistribution by plant root system. Over the energy‐limited mesic zone, the dynamic root affects regional land‐atmosphere coupling mainly through changes in carbon allocation. This work highlights the importance of dynamic root representation in improving vegetation‐atmosphere simulations by enhancing predictions of water, energy, and carbon fluxes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Effects of Horizontal Resolution, Domain Size, Boundary Conditions, and Surface Heterogeneity on Coarse LES of a Convective Boundary Layer

Atmospheric properties in a convective boundary layer vary over a wide range of spatial scales and are commonly studied using large-eddy simulations (LES) in various configurations. We examine how the boundary layer depth and distribution of variability across scales are affected by LES grid spacing, domain size, inhomogeneity of surface properties, and external forcing. Two different setups of the Weather Research and Forecasting (WRF) model are analyzed. A semi-idealized configuration uses a periodic domain, flat surface, prescribed homogeneous surface heat fluxes, and horizontally uniform profiles of large-scale advective tendencies. A nested LES setup employs a larger domain and realistic initial and boundary conditions, including an interactive land surface model with representative topography and vegetation and soil types. Subdomains of identical size are analyzed for all simulations. Characteristic structure sizes are quantified using the variability scales L 50 and L 95 , defined such that features smaller than that contain 50% and 95% of the total variance, respectively. Progressive increase in L 50 from vertical velocity to temperature and moisture structures is systematically reproduced in all simulation configurations. This dependence of L 50 on the considered variable complicates the development of scale-aware parameterizations for models with grid spacing in the “terra incognita”. In simulations using a larger domain with heterogeneous surface properties, the development of internal mesoscale patterns significantly affects variance distributions inside analyzed subdomains. Sizes of boundary layer structures also strongly depend on the LES grid spacing and, in case of heterogeneous surface and topography, on location of the subdomain inside a larger computational domain.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Grid Spacing Sensitivities of Simulated Mid-Latitude and Tropical Mesoscale Convective Systems in the Convective Gray Zone

The main objective of this study is to observationally constrain processes in tropical and midlatitude mesoscale convective systems (MCSs), and to use these constraints for model evaluation. To accomplish this, we leverage MCS observations collected at the U.S. DOE Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site in Oklahoma and ARM's mobile GoAmazon2014/15 site in Manaus, Brazil (MAO). We simulate 13 and 11 of these observed MCSs at the SGP and MAO site, respectively, using the Weather Research and Forecasting model at 12-, 4-, 2-, and 1-km horizontal grid spacing. Observations from radiosondes, surface meteorology, and radar wind profilers are used to characterize MCS properties, such as MCS timing and location, cold pools, and convective drafts, and evaluate these simulations. SGP cases are found in better agreement with observations than MAO cases, and when simulated at 2 km, outperform simulations at 1 km regarding the timing of MCS overpass and the accuracy of surface variable trends. MAO simulations suggest a consistent improvement in model accuracy with increasing model resolution in depicting the downdraft structure, the timing of MCSs, and the surface variables changes, except for the latter two metrics at 2 km. Deficiencies are still evident at km-scales, suggesting the need for higher resolution to simulate tropical MCSs. Overall, location-dependent improvements in MCS representation are obtained with the increasing model resolution, prompting the evaluation of sub-km scale simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

An efficient method to identify uncertainties of WRF-Solar variables in forecasting solar irradiance using a tangent linear sensitivity analysis

Uncertainty in predicting solar energy resources introduces major challenges in power system management and necessitates the development of reliable probabilistic solar forecasts. As the first part of the development of probabilistic forecasts based on the Weather Research and Forecasting model with solar extensions (WRF-Solar), this study presents a tangent linear approach to identify input variables responsible for the largest uncertainties in predicting surface solar irradiance and clouds. A tangent linear analysis is capable of efficiently investigating sensitivities of output variables with respect to various input variables of WRF-Solar because this approach avoids the computational burden of perturbing the initial conditions of individual input variables. We develop tangent linear models (TLMs) for six WRF-Solar physics packages that control the formation and dissipation of clouds and solar radiation, and we evaluate the validity of TLMs using a linearity test. The tangent linear sensitivity analysis is conducted under various scenarios based on satellite observations and model simulations to consider realistic input conditions. A simple method is used to quantify the impact of the uncertainty of input variables on the output variables from the TLMs. The results demonstrate that uncertainties in the output variables that are the focus of this study—including global horizontal irradiance, direct normal irradiance, cloud mixing ratio, cloud tendency, cloud fraction, and sensible and latent heat fluxes—are highly sensitive to uncertainties in 14 input variables. This study indicates that the tangent linear method can identify key variables of physics modules in WRF-Solar that can be stochastically perturbed to generate ensemble-based probabilistic forecasts.

14 SOLAR ENERGY↗

DOE OSTI2021

Simulated aging processes of black carbon and its impact during a severe winter haze event in the Beijing-Tianjin- Hebei region

Black carbon (BC) can mitigate or worsen air pollution through perturbing meteorological conditions. BC aging processes are important for the evolution of particle size, concentration, and optical properties of BC that determines its influence on the meteorology. Here we use the online coupled Weather Research and Forecasting-Chemistry (WRF-Chem) model to quantify the role of BC aging processes, including physical processes (PP) and absorption enhancements (AE), in exerting BC-induced meteorological changes and the associated feedback to PM2.5 (particulate matter less than 2.5 µm in diameter) and O3 concentrations during a severe haze event in Beijing-Tianjin-Hebei (BTH) region during 21-27 February 2014. Our results show that, compared to the simulation without the PP treatment, the simulated near-surface BC concentration and BC mass loading in BTH region is lowered by 6.6 % and 12.1 %, respectively, during the haze event with the PP included. PP increases the proportion of large-size BC (particle diameter greater than 0.312 µm) from 28 %-33 % to 59 %-64 % below 1000 m in BTH region. Both PP and AE enhance the “dome effect” of BC. With both PP and AE considered, a reduction in PBL height due to BC-PBL interaction is 116.3 m (20.7 %), compared to 75.7 m (13.5 %) without AE and 66.6 m (11.9 %) without both PP and AE. However, during this haze event, anomalous northeasterly winds are produced by the direct radiative effect of BC, which further affects the mixing and transport of aerosols. When PBL height is decreased (from 10:00 to 21:00), combining all the impacts on multiple meteorological factors, BC effect without PP and AE, without AE, and with PP and AE, respectively, increases surface concentrations of PM2.5 by 8.3 µg m-3 (6.1 % relative to the mean value), 6.1 µg m-3 (4.5 %) and 9.6 µg m-3 (7.0 %) but decreases surface O3 concentrations by 2.8 ppbv (7.4 %), 4.0 ppbv (9.0 %) and 5.0 ppbv (10.8 %) in BTH region averaged over 21-27 February 2014. Our results highlight the importance of the aging processes and absorption enhancements of BC in simulating weather and air quality.

Chen, Donglin↗

DOE OSTI2025

Prioritizing urban heat adaptation infrastructure based on multiple outcomes: Comfort, health, and energy

Globally, cities face increasing extreme heat, impacting comfort, health, and energy consumption. Infrastructure-based heat adaptation strategies can improve these outcomes, but each strategy has a unique mix of benefits and drawbacks. Here, we apply an urbanized meteorological model (WRF) with the newly integrated multilayer BEP-Tree street tree model to dynamically downscale Earth System Model projections and a 3-D microclimate model (TUF-Pedestrian) to simulate the street-scale radiation environment impacting pedestrians. We evaluate the performance of five heat adaptation strategies (street trees, cool roofs, green roofs, rooftop photovoltaics (PV), and reflective pavements) during extreme heat events in three cities with contrasting background climates (Toronto, Phoenix, and Miami), under contemporary and end-of-century projected climates, based on three metrics: outdoor heat stress, air conditioning (AC) energy use, and ventilation of vehicular air pollution. No single adaptation strategy improves all three outcomes. While street trees inhibit ventilation, they reduce outdoor heat stress four times more effectively than the next best strategy via shade provision, fully offsetting heat stress increases under a high-emissions end-of-century climate scenario in all cities studied. Cool roofs and green roofs moderately reduce heat stress and energy use. Alternatively, rooftop PV with energy storage can generate sufficient power for space cooling but have marginal effects on heat stress. Reflective pavements are the least effective across metrics. Where the ventilation of street-level emissions is of less concern, our results clearly support the combination of street trees and rooftop PV as a highly complementary and effective means of adaptive mitigation across different climates and neighborhood densities.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

Synthesis of observed and simulated rain microphysics to inform a new Bayesian statistical framework for microphysical parameterization in climate models (Final Scientific Report)

This project aimed to investigate a new approach to bulk microphysics schemes. We identified the need to move beyond the fixed structural assumptions and approximation of existing schemes. For example, most schemes assume some functional form for the rain drop size distribution (e.g. a gamma or exponential distribution). Most schemes then evolve some number of statistical moments of that distribution via various processes, such as evaporation, sedimentation, collision-coalescence, and collisional breakup. These microphysical processes, in turn, are typically some fixed functional form derived from either some other (more detailed) model, or using some single-particle rates that are then integrated over the size distribution. While these process rate functions may have some free parameters to adjust, in other cases doing so is impossible (e.g. when the process rates are analytical or piecewise solutions to some target function). Our goal was to build and test a scheme that assumed no size distribution form, and used a series of power laws as the basis for the process rates. Moments of the size distribution would be predicted, but no underlying size distribution would be specified. Any number or choice of size distribution moments could be used, and any number of power laws could be employed to model the process rates. Thus, our approach is flexible, and can seamlessly scale across levels of complexity, as demanded by the data. Bayesian inference then would provide the formalism to estimate the model parameters and structure, allowing for robust uncertainty quantification. We proposed testing this framework in idealized simulations using bin microphysical schemes as a data source. We also proposed using real data to inform our microphysics scheme, and also that we would integrate our scheme into WRF.

58 GEOSCIENCES↗

DOE OSTI2023

LES ARM Symbiotic Simulation and Observation Scenario for CACTI

The Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) deep-convection scenario for the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. The LASSO-CACTI scenario consists of a library of mesoscale and LES simulations spanning grid spacings from 7.5 km to 100 m for convection near the Sierras de Córdoba mountain range, roughly centered on the ARM Mobile Facility. Supporting skill scores and observation comparisons complement the simulations in addition to subsets of the model output for various categories of model output variables. More information can be found at https://www.arm.gov/capabilities/modeling/lasso.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Developing Metrics to Evaluate the Skill and Credibility of Downscaling

In this project, our overarching goal was to evaluate the skill and credibility of statistical and dynamical downscaling techniques for reconstructing the historical climate record and projecting future climate change. This report summarizes our approach and accomplishments. In the first year of the project, the priority was to generate new data sets. For Year 2, the overall goal was to focus comparisons during the historical period and for a subset of variables, especially temperature and precipitation. In Year 3 and beyond, the focus was mainly on comparisons and metrics for future climate change simulations. To quantify the uncertainty associated with gridded dataset choice, the UCLA group compared eight gridded temperature datasets, examining climatologies, variability, and trends. The Scripps group completed statistical downscaling of NARR with the LOCA method, and UCLA and Scripps collaborated on a comparison of historical and end-of-century temperature and precipitation patterns over California from WRF, LOCA-WRF, and LOCA-Livneh. In addition to downscaling temperature and precipitation, the Scripps group has spent considerable time working on near-surface wind and humidity downscaling. The Santa Clara group assembled the Variable Infiltration Capacity (VIC) for nine snow-dominated river basins draining the western slopes of the Sierra Nevada range and assessed the role of snow albedo feedback in their VIC simulations. The results will help inform future downscaling studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Developing Metrics to Evaluate the Skill and Credibility of Downscaling (Final Report)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Evaluation of precipitation indices in suites of dynamically and statistically downscaled regional climate models over Florida

Abstract The present work evaluates historical precipitation and its indices defined by the Expert Team on Climate Change Detection and Indices (ETCCDI) in suites of dynamically and statistically downscaled regional climate models (RCMs) against NOAA’s Global Historical Climatology Network Daily (GHCN-Daily) dataset over Florida. The models examined here are: (1) nested RCMs involved in the North American CORDEX (NA-CORDEX) program, (2) variable resolution Community Earth System Models (VR-CESM), (3) Coupled Model Intercomparison Project phase 5 (CMIP5) models statistically downscaled using localized constructed analogs (LOCA) technique. To quantify observational uncertainty, three in situ-based (PRISM, Livneh, CPC) and three reanalysis (ERA5, MERRA2, NARR) datasets are also evaluated against the station data. The reanalyses and dynamically downscaled RCMs generally underestimate the magnitude of the monthly precipitation and the frequency of the extreme rainfall in summer. The models forced with CanESM2 miss the phase of the seasonality of extreme precipitation. All models and reanalyses severely underestimate both the mean and interannual variability of mean wet-day precipitation (SDII), consecutive dry days (CDD), and overestimate consecutive wet days (CWD). Metric analysis suggests large uncertainty across NA-CORDEX models. Both the LOCA and VR-CESM models perform better than the majority of models. Overall, RegCM4 and WRF models perform poorer than the median model performance. The performance uncertainty across models is comparable to that in the reanalyses. Specifically, NARR performs poorer than the median model performance in simulating the mean indices and MERRA2 performs worse than the majority of models in capturing the interannual variability of the indices.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Evaluating Implementations of the Immersed Boundary Method in the Weather Research and Forecasting Model

The terrain-following coordinate system used by many atmospheric models can cause numerical instabilities due to discretization errors as resolved terrain slopes increase and the grid becomes highly skewed. The immersed boundary (IB) method, which does not require the grid to conform to the terrain, has been shown to alleviate these errors, and has been used successfully for high-resolution atmospheric simulations over steep terrain, including vertical building surfaces. Since many previous applications of IB methods to atmospheric models have used very fine grid resolution (5 m or less), the present study seeks to evaluate IB method performance over a range of grid resolutions and aspect ratios. Two classes of IB algorithms, velocity reconstruction and shear stress reconstruction, are tested within the common framework of the Weather Research and Forecasting (WRF) Model. Performance is evaluated in two test cases, one with flat terrain and the other with the topography of Askervein Hill, both under neutrally stratified conditions. WRF-IB results are compared to similarity theory, observations, and native WRF results. Despite sensitivity to the location at which the IB intersects the model grid, the velocity reconstruction IB method shows consistent performance when used with a hybrid RANS/LES surface scheme. The shear stress reconstruction IB method is not sensitive to the grid intersection, but is less consistent and near-surface velocity errors can occur at coarse resolutions. This study represents an initial investigation of IB method variability across grid resolutions in WRF. Future work will focus on improving IB method performance at intermediate to coarse resolutions.

54 ENVIRONMENTAL SCIENCES↗