Engineering Papers⌕ Search

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.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Impacts of Rotational Mixing on the Precipitation Simulated by a Convection Permitting Model

With increased availability of computational resources, regional and global scale convection-permitting model (CPM, Δx ~ 1–10 km) simulations are becoming more common. CPMs have improved accuracy in their representation of deep convection and mesoscale convective systems (MCSs) compared to coarser resolution models. However, CPMs still exhibit convective cloud and precipitation biases relative to observations, notably a lesser frequency of light precipitation rates and greater frequency of heavy precipitation rates. In this work we hypothesize that these CPM biases are related to under-resolved mixing between convective updrafts and their surrounding environment. To test this hypothesis, we introduce a parameterization to the Weather Research and Forecasting model (WRF) that adds a small angular rotation of the grid-scale flow about the axis perpendicular to the plane of convective drafts. This rotated flow is then allowed to alter advection of moisture and hydrometeors. The effects of such mixing on precipitation characteristics are evaluated in month-long 4-km grid spacing simulations over the Amazon. The enhanced mixing transports moisture and condensate from convective cores to other areas including downdrafts. This increases the frequency of low-precipitable water and light precipitation. It also decreases the frequency of intense precipitation from isolated deep convection and MCSs, increases cloud top temperatures, reduces radar echo-top heights, and increases overall precipitation by altering the relationship of precipitation with precipitable water, in better agreement with observations. The results suggest when optimized using multiple observations, such an approach may provide a path toward more accurate representation of convection and precipitation statistics in convection-permitting simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Large-eddy simulation of an atmospheric bore and associated gravity wave effects on wind farm performance in the southern Great Plains

Gravity waves are a common occurrence in the atmosphere, with a variety of generation mechanisms. Their impact on wind farms has only recently gained attention, with most studies focused on wind farm-induced gravity waves. In this study, the interaction between a wind farm and gravity waves generated by an atmospheric bore event is assessed using multiscale large-eddy simulations. The atmospheric bore is created by a thunderstorm downdraft from a nocturnal mesoscale convective system (MCS). The associated gravity waves impact the wind resource and power production at a nearby wind farm during the American Wake Experiment (AWAKEN) in the US southern Great Plains. A two-domain nested setup (Δx=300 and 20 m) is used in the Weather Research and Forecasting (WRF) model, forced with data from the High-Resolution Rapid Refresh model, to capture both the formation of the bore and its interaction with individual wind turbines. The MCS is resolved on the large outer domain, where the structure of the bore and the associated gravity waves are found to be especially sensitive to parameterized microphysics processes. On the finer inner domain, gravity wave interactions with individual wind turbines are resolved; wake dynamics are captured using a generalized actuator disk parameterization in WRF. The gravity waves are found to have a strong effect on the atmosphere above the wind farm; however, the effect of the waves is more nuanced closer to the surface where there is additional turbulence, both ambient and wake-generated. Notably, the gravity waves modulate the mesoscale environment by weakening and dissipating the preexisting low-level jet, which reduces hub-height wind speed and hence the simulated power output, which is confirmed by the observed supervisory control and data acquisition (SCADA) power data. Additionally, the gravity waves induce local wind direction variations correlated with fluctuations in pressure, which lead to fluctuations in the simulated power output as various turbines within the farm are subjected to waking from nearby turbines.

17 WIND ENERGY↗
DOE OSTI2023

New insights on wind turbine wakes from large‐eddy simulation: Wake contraction, dual nature, and temperature effects

Abstract Large‐eddy simulation (LES) has been adopted to study wind turbine wakes because it can capture fine‐scale details of turbulent wind flows and interactions with wind turbines. Here, we use the LES version of the Weather Research and Forecasting (WRF) model with an actuator disk model to gain insights on several wake effects that have been traditionally difficult to measure. The first finding is that the wake has a “dual nature,” meaning that the wind speed deficit behaves differently from the added turbulent kinetic energy (TKE) and the two are not co‐located in space. For example, the wind speed deficit peaks at hub height and reaches the ground within 8D (D is the rotor diameter), but added TKE peaks near the rotor tip and generally remains aloft. Second, temperature changes near the ground are driven by the added TKE in the rotor area and by atmospheric stability. The combination of these two factors determines the sign and intensity of the vertical heat flux divergence below the rotor, with convergence and warming associated with stable conditions and weak divergence and modest cooling with unstable conditions. Third, wakes do not expand indefinitely, as suggested by similarity theory applied to the wind speed deficit, but eventually stop expanding and actually contract, at different rates depending on atmospheric stability. The implication of these findings is that, in order to study wakes, it is not sufficient to focus on wind speed deficit alone, because TKE is also important and yet behaves differently from the wind speed deficit.

17 WIND ENERGY↗
DOE OSTI2024

WTK-LED: The WIND Toolkit Long-Term Ensemble Dataset

To satisfy a wide group of stakeholders across various wind energy disciplines, including but not limited to stakeholders in the distributed and utility scale wind industry, the new emerging airborne wind energy field, grid integration, power systems modeling, environmental modeling, and researchers in academia, and to close some of the gaps that current public datasets have, we aimed at developing an updated version of the meteorological WIND Toolkit, named WIND Toolkit Long-term Ensemble Dataset (WTK-LED), which is a meteorological dataset providing time series every 5 min and 2 km, including model uncertainty of wind speed at every modeling grid point so that users are provided with a range of possible wind speeds every 2 km. The data were produced using the Weather Research and Forecasting Model (WRF). The vertical grid used in WTK-LED includes many vertical layers in the atmospheric boundary layer to provide information of atmospheric quantities across the rotor layer of utility scale and distributed wind turbines. The WTK-LED includes: 1) Numerical simulations covering the continental United States, Alaska, and Hawaii, with high-resolution data being available for 3 years (2018-2020). 2) Climate simulations from Argonne National Laboratories covering the North American continent, including Alaska, Canada, and most of Mexico and the Caribbean Islands. These simulations complement the new WTK-LED to offer a 4-km dataset covering 20 years, from 2001-2020. 3) Specific long-term,high-resolution offshore simulations have been conducted separately for the US coasts, Hawaii, and the Great Lakes, leading to the 2023 National Offshore Wind data set. This report focuses on a description of the land-based WTK-LED for CONUS, Hawaii, and Alaska, for the 3-year 2-km/5-min dataset and the 20-year 4-km/hourly dataset, as well as the uncertainty quantification method. We also provide limited validation results. Based on our results to date, we suggest use cases and applications for each dataset of the WTK-LED.

17 WIND ENERGY↗
DOE OSTI2023

Modelling the coupled mercury-halogen-ozone cycle in the central Arctic during spring

Near-surface mercury and ozone depletion events occur in the lowest part of the atmosphere during Arctic spring. Mercury depletion is the first step in a process that transforms long-lived elemental mercury to more reactive forms within the Arctic that are deposited to the cryosphere, ocean, and other surfaces, which can ultimately get integrated into the Arctic food web. Depletion of both mercury and ozone occur due to the presence of reactive halogen radicals that are released from snow, ice, and aerosols. In this work, we added a detailed description of the Arctic atmospheric mercury cycle to our recently published version of the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem 4.3.3) that includes Arctic bromine and chlorine chemistry and activation/recycling on snow and aerosols. The major advantage of our modelling approach is the online calculation of bromine concentrations and emission/recycling that is required to simulate the hourly and daily variability of Arctic mercury depletion. We used this model to study coupling between reactive cycling of mercury, ozone, and bromine during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) spring season in 2020 and evaluated results compared to land-based, ship-based, and remote sensing observations. The model predicts that elemental mercury oxidation is driven largely by bromine chemistry and that particulate mercury is the major form of oxidized mercury. The model predicts that the majority (74%) of oxidized mercury deposited to land-based snow is re-emitted to the atmosphere as gaseous elemental mercury, while a minor fraction (4%) of oxidized mercury that is deposited to sea ice is re-emitted during spring. Our work demonstrates that hourly differences in bromine/ozone chemistry in the atmosphere must be considered to capture the springtime Arctic mercury cycle, including its integration into the cryosphere and ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Quantifying the Resolution Sensitivity of the Kain–Fritsch Scheme Across the Gray Zone by Isolating Interactions: A TWP‐ICE Case Study

The resolution sensitivity of the Kain–Fritsch (KF) convection scheme and the role of interactions between the physics and dynamics within the gray zone (<10 km) were investigated using the Separate Physics and Dynamics Experiment (SPADE) framework. Two groups of experiments were conducted using the Weather Research and Forecasting (WRF) model via traditional (Tradition) runs and SPADE runs with resolutions of 1, 2, 4, and 8 km during the wet period of the Tropical Warm Pool–International Cloud Experiment (TWP‐ICE). Results show that the KF scheme simulates the weakened convective processes well as the resolution increases in both groups, and the changes in the convective variables with resolution in SPADE are smaller than in the Tradition group. This indicates the important effects of interactions between model components on convection parameterizations as the resolution changes. Additionally, the microphysics variables remain nearly unchanged with resolution in SPADE and weaken slightly in Tradition as the resolution decreases, suggesting the relatively weaker influences of model interactions for the resolved‐cloud parameterization. Therefore, the scale‐aware behavior of KF scheme is further strengthened in Tradition runs, primarily through inhibiting the strength of stratiform processes through physics–dynamics interactions and physical components.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Impact of Lateral Flow on Surface Water and Energy Budgets over the Southern Great Plains – A Modeling Study

As we see the horizontal grid spacing decrease, treatment of hydrologic processes in land surface models, such as the lateral flow of surface and subsurface flow, need to be explicitly represented. Unlike previous studies that mainly focused on the mountainous regions, in this study the offline WRF-Hydro model is employed to study the impact of lateral flow on soil moisture and energy fluxes over the relatively flat southern Great Plains (SGP). The vast amount of measurements over the SGP provide an unique opportunity to assess the model behavior. In addition, newly developed land surface properties and input forcing are ingested into the model, in an attempt to reduce uncertainties associated with the initial and boundary forcing and help to identify model deficiencies. Our results show that the more realistic inputs (parameters, soil types, forcing) lead to larger underestimation of latent heat flux and dry bias, indicating the existence of model structural uncertainty (embedded errors) in WRF-Hydro that need to be characterized to inform future model development efforts. Including lateral flow processes partly mitigates the model deficiencies in representing hydrologic processes and alleviates the dry bias. In particular, both surface and subsurface lateral flow increase soil moisture mainly over the lower elevations, except that subsurface flow also affects soil moisture over steeper terrains. Additional simulations are performed to assess the effect of routing resolution on model results. When LSM resolution is high, noticeable differences in soil moisture are produced between different routing resolutions especially over steep terrain. Whereas when LSM resolution is coarse, differences between routing resolutions become negligible, especially over flat terrain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Uncertainty quantification and optimization of precipitating hydrometeor parameters for winter precipitation in a cloud microphysics scheme

The precipitating hydrometeor parameters used in cloud microphysics schemes carry inherent uncertainties. The quantification of these uncertainties, together with parameter optimization, can significantly improve precipitation forecasts. This study investigates the effects of 13 parameters in the Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) microphysics scheme, which define the hydrometeor characteristics such as fall velocity–diameter and mass–diameter relationships, as well as the shape parameter of the drop size distribution for precipitating particles such as rain, snow, and graupel on simulated winter precipitation. A comparison between the model's pre-defined parameters and observations from the International Collaborative Experiments for the PyeongChang 2018 Olympic and Paralympic winter games (ICE-POP 2018) field campaign reveals that the fall velocity–diameter relationship for rain, the mass–diameter relationships for snow and graupel, and the shape parameters for all precipitating particles in the WDM6 scheme deviate from the median values observed by the two-dimensional video disdrometer (2DVD). To quantify parameter sensitivities, a perturbed parameter ensemble (PPE) of 256 simulations was conducted within parameter ranges constrained by 2DVD observations for three winter precipitation cases. Bayesian optimization was then applied to identify parameter sets that minimized the root mean square error (RMSE) for each case, achieving reductions of up to 30.2 %. These results demonstrate that ensemble-based uncertainty quantification and parameter optimization can help identify key parameters and provide a pathway to improving precipitation simulation performance. In addition, measurement sites can be strategically selected based on regions that show high sensitivity to variations in hydrometeor characteristic parameters.

Bayesian optimization↗
DOE OSTI2022

SOWFA Large-Eddy Simulation of a Diurnal Cycle at the SWiFT Site

These are large-eddy simulation (LES) data output from the Simulator for Wind Farm Applications (SOWFA) with internal mesoscale-to-microscale coupling. A diurnal cycle was simulated over the Scaled Wind Farm Technology (SWIFT) facility near Lubbock, Texas, from 2013-11-08 12:00:00 UTC to 2013-11-09 12:00:00 UTC. Indirect profile assimilation was used to enforce nonstationary mesoscale conditions. Mesoscale forcing data were extracted from WRF numerical weather prediction outputs or derived from observations at SWiFT.

17 WIND ENERGY↗
DOE OSTI2021

Improving the prediction of DNI with physics-based representation of all-sky circumsolar radiation

Direct normal irradiance (DNI) is often interpreted differently in ground measurements and forecasts using numerical weather prediction (NWP) models, leading to substantial bias in DNI forecasting especially under cloudy-sky conditions. To mitigate the bias, we use the Fast All-sky Radiation Model for Solar applications with DNI (FARMS-DNI) to conduct physics-based simulations of solar radiation in the circumsolar region. The DNI is quantified using the sum of the radiation along the sun direction and the scattered radiation within the circumsolar region. FARMS-DNI is coupled with the Weather Research and Forecasting model with solar extensions (WRF-Solar) to predict day-ahead DNI in 9-km pixels over an extended domain covering the contiguous United States. The DNI forecasts for the entire year of 2018 are validated using surface-based observations at the Atmospheric Radiation Measurement (ARM) and Surface Radiation Budget Network (SURFRAD) sites and a satellite-based solar radiation product developed by the procedure of the National Solar Radiation Data Base (NSRDB). The results show that the WRF-Solar coupled with the conventional Beer-Bouguer-Lambert law underestimes the overall DNI forecasts by 30%-60%, which is a joint effect of the uncertainties in computing the circumsolar radiation and forecasting cloud properties. This bias is significantly reduced by using FARMS-DNI.

14 SOLAR ENERGY↗
DOE OSTI2021

Evaluating the performance of WRF urban schemes and PBL schemes over Dallas-Fort Worth during a dry summer and a wet summer

This study evaluated the Weather Research and Forecasting (WRF) model sensitivity to different planetary boundary layer (PBL) schemes (the YSU and MYJ schemes) and urban schemes including the bulk scheme (BULK), single-layer urban canopy model (UCM), multi-layer building environment parameterization (BEP) model, and multi-layer building energy model (BEM). Daily reinitialization simulations were conducted over Dallas-Fort Worth during a dry summer month (July 2011) and a wet summer month (July 2015) with weaker (stronger) daytime (nocturnal) UHI in 2011 than 2015. All urban schemes overestimated the urban daytime 2m temperature in both summers, but BEP and BEM still reproduced the daytime urban cool island in dry summer. All urban schemes reproduced the nocturnal urban heat island, with BEP producing the weakest one due to its unrealistic urban cooling. BULK and UCM overestimated the urban canopy wind speed, while BEP and BEM underestimated it. The urban schemes showed prominent impact on daytime PBL profiles. UCM+MYJ showed a superior performance than other configurations. The relatively large (small) aspect ratio between building height and road width in UCM (BEM) was responsible for the overprediction (underprediction) of urban canopy temperature. The relatively low (high) building height in UCM (BEM) was responsible for the overprediction (underprediction) of urban canopy wind speed. Improving urban schemes and providing realistic urban parameters were critical for improving urban canopy simulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Improved Prediction of Cold-Air Pools in the Weather Research and Forecasting Model Using a Truly Horizontal Diffusion Scheme for Potential Temperature

The terrain-following vertical coordinate system used by many atmospheric models, including the Weather Research and Forecasting (WRF) Model, is prone to errors in regions of complex terrain. These errors stem, in part, from the calculation of horizontal gradients within the diffusion term of the momentum or scalar evolution equations. In WRF, such gradients can be calculated along coordinate surfaces, or using metric terms that help account for grid skewness. However, neither of these options ensures a truly horizontal gradient calculation, especially if a grid cell is skewed enough that the heights of the neighboring grid points used in the calculation fall outside the vertical range of the cell. In this work, an improved scheme that uses Taylor series approximations to vertically interpolate variables to the level necessary for a truly horizontal gradient calculation is implemented in WRF for the diffusion of potential temperature. The scheme is validated using an atmosphere-at-rest configuration, in which spurious flows develop only as a result of numerical errors and can thus be used as a proxy for model performance. Following validation, the method is applied to the simulation of cold-air pools (CAPs), which occur in regions of complex terrain and are characterized by strong near-surface temperature gradients. Using the truly horizontal scheme, idealized simulations demonstrate reduced numerical mixing in a quiescent CAP, and a realistic case study in the Columbia River basin shows a reduction in positive wind speed bias by up to roughly 20% compared to observations from the Second Wind Forecast Improvement Project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Impacts of cloud microphysics parameterizations on simulated aerosol–cloud interactions for deep convective clouds over Houston

Aerosol–cloud interactions remain largely uncertain with respect to predicting their impacts on weather and climate. Cloud microphysics parameterization is one of the factors leading to large uncertainty. Here, we investigate the impacts of anthropogenic aerosols on the convective intensity and precipitation of a thunderstorm occurring on 19 June 2013 over Houston with the Chemistry version of Weather Research and Forecast model (WRF-Chem) using the Morrison two-moment bulk scheme and spectral bin microphysics (SBM) scheme. We find that the SBM predicts a deep convective cloud that shows better agreement with observations in terms of reflectivity and precipitation compared with the Morrison bulk scheme that has been used in many weather and climate models. With the SBM scheme, we see a significant invigoration effect on convective intensity and precipitation by anthropogenic aerosols, mainly through enhanced condensation latent heating. Such an effect is absent with the Morrison two-moment bulk microphysics, mainly because the saturation adjustment approach for droplet condensation and evaporation calculation limits the enhancement by aerosols in (1) condensation latent heat by removing the dependence of condensation on droplets and aerosols and (2) ice-related processes because the approach leads to stronger warm rain and weaker ice processes than the explicit supersaturation approach.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Observations and simulations of a wind farm modifying a thunderstorm outflow boundary

Abstract. On 18 June 2019, National Weather Service (NWS) radar reflectivity data indicated the presence of thunderstorm-generated outflow propagating east-southeastward near Lubbock, Texas. A section of the outflow boundary encountered a wind farm and then experienced a notable reduction in ground-relative velocity, suggesting that interactions with the wind farm impacted the outflow boundary progression. We use the Weather Research and Forecasting model and its wind farm parameterization to address the extent to which wind farms can modify the near-surface environment of thunderstorm outflow boundaries. We conduct two simulations of the June 2019 outflow event: one containing the wind farm and one without. We specifically investigate the outflow speed of the section of the boundary that encounters the wind farm and the associated impacts on near-surface wind speed, moisture, temperature, and changes to precipitation features as the storm and associated outflow pass over the wind farm domain. The NWS radar and nearby West Texas Mesonet surface stations provide observations for validation of the simulations. The presence of the wind farm in the simulation clearly slows the progress of the outflow boundary by over 20 km h -1 , similar to what was observed. Simulated perturbations of surface wind speed, temperature, and moisture associated with outflow passage were delayed by up to 6 min when the wind farm was present in the simulation compared to the simulation without the wind farm. However, impacts on precipitation were localized and transient, with no change to total accumulation across the domain.

17 WIND ENERGY↗
DOE OSTI2020

Simulation of Continental Shallow Cumulus Populations Using an Observation-Constrained Cloud-System Resolving Model

Continental shallow cumulus (ShCu) clouds observed on 30 August 2016 during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign are simulated by using an observation-constrained cloud-system resolving model. On this day, ShCu forms over Oklahoma and southern Kansas and some of these clouds transition to deeper, precipitating convection during the afternoon. We apply a four-dimensional ensemble-variational (4DEnVar) hybrid technique in the Community Gridpoint Statistical Interpolation (GSI) system to assimilate operational data sets and unique boundary layer measurements including a Raman lidar, radar wind profilers, radiosondes, and surface stations collected by the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) atmospheric observatory into the Weather Research and Forecasting (WRF) model to ascertain how improved environmental conditions can influence forecasts of ShCu populations and the transition to deeper convection. Independent observations from aircraft, satellite, as well as ARM's remote sensors are used to evaluate model performance in different aspects. Several model experiments are conducted to identify the impact of data assimilation (DA) on the prediction of clouds evolution. The analyses indicate that ShCu populations are more accurately reproduced after DA in terms of cloud initiation time and cloud base height, which can be attributed to an improved representation of the ambient meteorological conditions and the convective boundary layer. Extending the assimilation to 18 UTC (local noon) also improved the simulation of shallow-to-deep transitions of convective clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

SSP‐Based Land‐Use Change Scenarios: A Critical Uncertainty in Future Regional Climate Change Projections

Abstract To better understand the role projected land‐use changes (LUCs) may play in future regional climate projections, we assess the combined effects of greenhouse‐gas (GHG)‐forced climate change and LUCs in regional climate model (RCM) simulations. To do so, we produced RCM simulations that are complementary to the North‐American Coordinated Regional Downscaling Experiment (NA‐CORDEX) simulations, but with future LUCs that are consistent with particular Shared Socioeconomic Pathways (SSPs) and related to a specific Representative Concentration Pathway (RCP). We examine the state of the climate at the end of the 21st century with and without two urban and agricultural LUC scenarios that follow SSP3 and SSP5 using the Weather Research and Forecasting (WRF) model forced by one global climate model, the MPI‐ESM, under the RCP8.5 scenario. We find that LUCs following different societal trends under the SSPs can significantly affect climate projections in different ways. In regions of significant cropland expansion over previously forested area, projected annual mean temperature increases are diminished by around 0.5°C–1.0°C. Across all seasons, where urbanization is high, projected temperature increases are magnified. In particular, summer mean temperature projections are up to 4°C–5°C greater and minimum and maximum temperature projections are increased by 2.5°C–6°C, amounts that are on par with the warming due to GHG‐forced climate change. Warming is also enhanced in the urban surroundings. Future urbanization also has a large influence on precipitation projections during summer, increasing storm intensity, event length, and the overall amount over urbanized areas, and decreasing precipitation in surrounding areas.

Bukovsky, Melissa S.↗
DOE OSTI2022

Simulation of the effects of biomass burning in a mesoscale convective system in the central amazon

During the dry season of 2014, the formation of a mesoscale convective system (MCS) caused intense precipitation and strong winds in the central Amazon region. In this period, cases of MCS that occurred during the days when there were higher concentrations of CO were analyzed. Through this criterion, a case of MCS occurred on August 16th, 2014 was selected. We used the chemical-atmospheric model WRF-Chem to assess the influence of biomass burning aerosol on the intensity of precipitation, winds, vertical and horizontal transport associated with this convective system. We show that biomass burning aerosol reduces the strength of the mesoscale convective system, with less vertical and horizontal transport of carbon monoxide and ozone. Here, in the absence of biomass burning aerosol downdrafts and precipitation rate are more intense, and cloud development up to 5km is more developed, making the horizontal flow and vertical transport of ozone more intense, however, much more efficient in reducing the concentrations of other gases.

54 ENVIRONMENTAL SCIENCES↗