Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation models”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Shallow-Water Model for Convective Self-Aggregation

Randomly distributed convective storms can self-aggregate in the absence of large-scale forcings. Here we present a 1D shallow-water model to study the convective self-aggregation. This model simulates the dynamics of the planetary boundary layer and represents convection as a triggered process. Once triggered, convection lasts for finite time and occupies finite length. Here we show that the model can successfully simulate self-aggregation, and that the results are robust to a wide range of parameter values. In the simulations, convection excites gravity waves. The gravity waves then form a standing wave pattern, separating the domain into convectively active and inactive regions. We analyze the available potential energy (APE) budget and show that convection generates APE, providing energy for self-aggregation. By performing dimensional analysis, we develop a scaling theory for the size of convective aggregation, which is set by the gravity wave speed, damping time scale, and number density of convective storms. This paper provides a simple modeling framework to further study convective self-aggregation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

V1G Frequency Regulation: Algorithm Development, Validation & Analysis at Scale

Researchers at Argonne National Laboratory developed and validated a high-fidelity digital twin of a smart charging (V1G) ecosystem to model the participation of up to 1,000 unique electric vehicles (EVs) in the PJM frequency regulation market. Utilizing a discrete-event framework, the simulation models complex interactions, from dynamic grid signals (updated every 2 seconds) to individual EV charging dynamics. The simulation incorporates multiple EV models created from real-world lab test data. Researchers tested multiple control algorithms to balance the dual objectives of maximizing aggregator’s revenue and driver charging needs. Results demonstrate that aggregated EVs function as a controllable, highly effective grid resource, achieving high PJM Performance Scores (80–90%). Additionally, an optimized, market-aware bidding strategy was identified as key to profitability. The platform was shown to provide drivers with an average charging discount of nearly 50%. The algorithm was further validated in the lab using production EVs and charging stations to compare simulation results with real-world performance.

Manne, Nithin↗
DOE OSTI2024

Bimodality in simulated precipitation frequency distributions and its relationship with convective parameterizations

Abstract Bimodality in precipitation frequency distributions is often evident in atmospheric models, but rarely in observations. This study i) proposes a metric to objectively quantify the bimodality in precipitation distributions, ii) evaluates model simulations contributed to the Coupled Model Intercomparison Project (CMIP) phase 5 (CMIP5), phase 6 (CMIP6), and the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) project by comparing them to satellite-based and reanalysis precipitation products, and iii) investigates possible origins of bimodal precipitation distributions. Our results reveal that about 83% (20 out of 24) of CMIP5 and 70% (21 out of 30) of CMIP6 models used in this study exhibit bimodal distributions. The few DYAMOND models that use a deep convective parameterization also show bimodal distributions, while most DYAMOND models do not. Predictably, the bimodality originates from the separation of precipitation process between resolved grid-scale and parameterized subgrid-scale. However, in a larger number of models bimodality arises from the parameterized subgrid-scale convective precipitation alone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Non‐Intrusive Machine Learning Framework for Debiasing Long‐Time Coarse Resolution Climate Simulations and Quantifying Rare Events Statistics

Abstract Due to the rapidly changing climate, the frequency and severity of extreme weather is expected to increase over the coming decades. As fully‐resolved climate simulations remain computationally intractable, policy makers must rely on coarse‐models to quantify risk for extremes. However, coarse models suffer from inherent bias due to the ignored “sub‐grid” scales. We propose a framework to non‐intrusively debias coarse‐resolution climate predictions using neural‐network (NN) correction operators. Previous efforts have attempted to train such operators using loss functions that match statistics. However, this approach falls short with events that have longer return period than that of the training data, since the reference statistics have not converged. Here, the scope is to formulate a learning method that allows for correction of dynamics and quantification of extreme events with longer return period than the training data. The key obstacle is the chaotic nature of the underlying dynamics. To overcome this challenge, we introduce a dynamical systems approach where the correction operator is trained using reference data and a coarse model simulation nudged toward that reference. The method is demonstrated on debiasing an under‐resolved quasi‐geostrophic model and the Energy Exascale Earth System Model (E3SM). For the former, our method enables the quantification of events that have return period two orders longer than the training data. For the latter, when trained on 8 years of ERA5 data, our approach is able to correct the coarse E3SM output to closely reflect the 36‐year ERA5 statistics for all prognostic variables and significantly reduce their spatial biases.

Barthel Sorensen, B.↗
DOE OSTI2025

A dynamic 2D Borehole Thermal Energy Storage (BTES) model for enhanced computational efficiency

Progressing toward a future increasingly reliant on renewable energy sources, the development of effective, durable energy storage solutions becomes essential to balance supply and demand fluctuations. Borehole Thermal Energy Storage (BTES) is a long-duration thermal energy storage technology that captures excess heat generated from renewable energy sources and stores it underground for later use, enabling the efficient utilization of sustainable energy. This approach is particularly valuable in district energy networks when integrated with Ground Source Heat Pumps (GSHP) to provide stable heating and cooling. However, traditional three-dimensional (3D) numerical models of BTES systems demand extensive computational resources, limiting their practicality for real-time and large-scale applications. This study introduces a novel two-dimensional (2D) modeling approach that reduces computational costs while maintaining high accuracy. By employing a radial ring-based discretization method, the model simulates heat injection, retention, and retrieval dynamics over seasonal cycles. A new thermal-mass weighted-average temperature parameter is introduced to evaluate the performance of BTES systems. Model validation against FEFLOW simulations demonstrates a 17-fold improvement in computational speed compared to traditional Computational Fluid Dynamics (CFD) models while achieving a mean absolute percentage error (MAPE) of 2 % during charging and 4 % during discharging. Additionally, a trade-off analysis between computational efficiency and accuracy is conducted, ensuring the model's applicability for real-world scenarios. The findings of this research contribute to the development of computationally efficient BTES models, facilitating better optimization, control, and integration into renewable energy systems. This work provides a foundation for further studies in techno-economic analysis, multi-year performance evaluation, and real-time operational strategies for BTES applications, supporting a more sustainable energy future.

2D modeling↗
DOE OSTI2024

Explicit Consideration of Plant Xylem Hydraulic Transport Improves the Simulation of Crop Response to Atmospheric Dryness in the U.S. Corn Belt

Abstract Atmospheric dryness (i.e., high vapor pressure deficit, VPD), together with soil moisture stress, limits plant photosynthesis and threatens ecosystem functioning. Regions where rainfall and soil moisture are relatively sufficient, such as the rainfed part of the U.S. Corn Belt, are especially prone to high VPD stress. With globally projected rising VPD under climate change, it is crucial to understand, simulate, and manage its negative impacts on agricultural ecosystems. However, most existing models simulating crop response to VPD are highly empirical and insufficient in capturing plant response to high VPD, and improved modeling approaches are urgently required. In this study, by leveraging recent advances in plant hydraulic theory, we demonstrate that the VPD constraints in the widely used coupled photosynthesis‐stomatal conductance models alone are inadequate to fully capture VPD stress effects. Incorporating plant xylem hydraulic transport significantly improves the simulation of transpiration under high VPD, even when soil moisture is sufficient. Our results indicate that the limited water transport capability from the plant root to the leaf stoma could be a major mechanism of plant response to high VPD stress. We then introduce a Demand‐side Hydraulic Limitation Factor (DHLF) that simplifies the xylem and the leaf segments of the plant hydraulic model to only one parameter yet captures the effect of plant hydraulic transport on transpiration response to high VPD with similar accuracy. We expect the improved understanding and modeling of crop response to high VPD to help contribute to better management and adaptation of agricultural systems in a changing climate.

Environmental Sciences & Ecology↗
DOE OSTI2022

Integration of Large-Scale Electrical Imaging into Geological Framework Development and Refinement

Geologic framework models (GFMs) are critical to the construction of reliable simulation models of groundwater flow and contaminant transport. To support GFM development, direct information (e.g., core samples, fluid samples, hydraulic testing) tends to be sparse and separated by large distances relative to the spatial scales of aquifer heterogeneity. There are additional challenges associated with highly contaminated legacy waste sites, where drilling is particularly costly, and invasive sampling requires specialized handling and disposal of hazardous materials. At these sites in particular, non-invasive geophysical imaging can play an important role in filling spatial gaps between boreholes and reducing characterization costs by optimizing and minimizing the number of necessary boreholes. Here this paper presents a case study demonstrating the use of large-scale (> 30 km 2 ) electrical mapping to identify hydrostratigraphy and potential paleochannels at the Hanford Site, located in Washington State, USA. In two field campaigns, over 36 line-kilometers of electrical resistivity tomography (ERT) data were collected along 14 transects. ERT surveys were sited and performed to image critical aspects (e.g., paleochannels, stratigraphic contacts) of the subsurface, demonstrating a general workflow for integrating ERT with GFM development. Inconsistencies between the GFM and ERT were catalogued to provide a basis for future site characterization using complementary geophysical methods and (or) direct sampling.

58 GEOSCIENCES↗
DOE OSTI2022

Impacts of climate change on estuarine stratification and implications for hypoxia within a shallow subtropical system

Vertical density stratification often plays an important role in the formation and expansion of coastal hypoxic zones through its effect on near-bed circulation and vertical oxygen flux. However, the impact of future climate change on estuarine circulation is widely unknown. Here, we developed and calibrated a three-dimensional hydrodynamic model for Pensacola Bay, a shallow subtropical estuary in the northeastern Gulf of Mexico. Model simulations based on years 2013–2017 were applied to examine changes in salinity, temperature, and density under future climate scenarios, including increased radiative forcing (IR) and temperature (T), increased freshwater discharge (D), sea level rise (SLR), and wind intensification (W). Simulations showed that the impacts of climate change on modeled state variables varied over time with external forcing conditions. The model demonstrated the potential for sea level rise and increased freshwater discharge to episodically increase vertical density gradients in the Bay. However, increased wind forcing destabilized vertical gradients, at times reducing the spatial extent and duration of stable stratification. For time periods with low freshwater discharge, moderate increases in wind speed (10%) can destabilize density gradients strengthened by increased discharge (10%) and sea level rise (0.48 m). In contrast, destruction of strong density gradients that form near the mid-Bay channel following flood events requires stronger wind forcing. Furthermore, these results highlight the importance of considering natural variability in freshwater and wind forcing, as well as local phenomena that are generally unresolved by global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impact of Rainfall on Tropical Cyclone-Induced Sea Surface Cooling

Tropical cyclones (TCs) are often accompanied by strong winds and torrential rains. While the winds associated with TCs tend to enhance mixing in the upper-ocean, the freshwater input from rain can stratify the water column and limit mixing. However, the extent to which the stabilizing effect of rainfall can compete with wind-induced mixing, and to what degree it modulates TC-induced sea surface cooling, remains unknown. Here we show, using a suite of observations, that heavy rains under weak TCs can significantly reduce the magnitude of cold wakes induced by them. Additionally, when compared to storms with low rain rates, the ocean surface under TCs with high rain rates freshens significantly and cools less. High-resolution climate model simulations and idealized experiments with an ocean mixed layer model support these results and reveal that oceanic mixing processes are primarily responsible for the reduced cooling under TCs, with a lesser role for surface fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Summertime low clouds mediate the impact of the large-scale circulation on Arctic sea ice

The rapid Arctic sea ice retreat in the early 21 st century is believed to be driven by several dynamic and thermodynamic feedbacks, such as ice-albedo feedback and water vapor feedback. However, the role of clouds in these feedbacks remains unclear since the causality between clouds and these processes is complex. Here, we use NASA CERES satellite products and NCAR CESM model simulations to suggest that summertime low clouds have played an important role in driving sea ice melt by amplifying the adiabatic warming induced by a stronger anticyclonic circulation aloft. The upper-level high pressure regulates low clouds through stronger downward motion and increasing lower troposphere relative humidity. The increased low clouds favor more sea ice melt via emitting stronger longwave radiation. Then decreased surface albedo triggers a positive ice-albedo feedback, which further enhances sea ice melt. Considering the importance of summertime low clouds, accurate simulation of this process is a prerequisite for climate models to produce reliable future projections of Arctic sea ice.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Characterization of laser-accelerated proton beams from a 0.5 kJ sub-picosecond laser for radiography applications

Laser-accelerated ion beams show promise for many applications, including high-resolution flash imaging of static or dynamic objects in next-generation radiography to probe materials and plasmas in extreme environments and inertial confinement fusion. To scale up ion beam production for radiography applications, we conducted experiments using sub-picosecond lasers up to 0.5 kJ at the OMEGA-EP facility to characterize proton beams from solid targets, primarily CH/CD sub-micron thin films from which ion beams were also used for static and dynamic radiography for the first time. For standalone sub-micron thin CH films, the highest detected proton energy is in the range of 72–97 MeV. Proton beams with highest energy near or above 60 MeV at full laser energy and similar beam profiles are also measured from low-density CD foams and flat CH foil target of micrometer-scale thickness. The ~ 700–800 nm CH/CD foils achieve the highest ion yield among the targets tested. For sub-micron thin films, the laser prepulse can expand the target and lead to complex interactions, which is simulated using coupled hydrodynamic and two-step kinetic models. Simulations suggest the presence of a micrometer-scale preplasma plateau with near-critical density and further indicate that target normal sheath acceleration, electron heating from Relativistic transparency in the preplasma plateau, and background proton reflection from carbon ion front at the rear side contribute to the resulting proton spectrum from these sub-micron thin targets at various stages. These proton beams show strong potential for radiography and for production of secondary sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Comparison of methods to estimate aerosol effective radiative forcings in climate models

Abstract. Uncertainty in the effective radiative forcing (ERF) of climate primarily arises from the unknown contribution of aerosols, which impact radiative fluxes directly and through modifying cloud properties. Climate model simulations with fixed sea surface temperatures but perturbed atmospheric aerosol loadings allow for an estimate of how strongly the planet's radiative energy budget has been perturbed by the increase in aerosols since pre-industrial times. The approximate partial radiative perturbation (APRP) technique further decomposes the contributions to the direct forcing due to aerosol scattering and absorption and to the indirect forcing due to aerosol-induced changes in cloud scattering, amount, and absorption, as well as the effects of aerosols on surface albedo. Here we evaluate previously published APRP-derived estimates of aerosol effective radiative forcings from these simulations conducted in the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and find that they are biased as a result of two large coding errors that – in most cases – fortuitously compensate. The most notable exception is the direct radiative forcing from absorbing aerosols, which is more than 40 % larger averaged across CMIP6 models in the present study. Correcting these biases eliminates the residuals and leads to better agreement with benchmark estimates derived from double calls to the radiation code. The APRP method – when properly implemented – remains a highly accurate and efficient technique for diagnosing aerosol ERF in cases where double radiation calls are not available, and in all cases it provides quantification of the individual contributors to the ERF that are highly useful but not otherwise available.

58 GEOSCIENCES↗
DOE OSTI2025

5-minute Wind Power Data based on WFIP2 WRF Simulation

The second Wind Forecast Improvement Project (WFIP2) was a public-private partnership funded by the U.S. Department of Energy and NOAA, aimed at enhancing the forecast skill of numerical weather prediction models for turbine-height winds in regions with complex terrain. An 18-month Weather Research and Forecasting (WRF) model simulation was conducted over the Pacific Northwest, with model outputs validated against observational data collected during WFIP2. Simulated wind speeds were used to estimate wind power generation using reV (the Renewable Energy Potential Model developed by NREL) at ten wind project sites. Two sets of results were produced: one using wind speeds extracted from the model grid cell at the project centroid, and another using wind speeds from the actual turbine locations. For each dataset, power output was calculated using both actual turbine-specific power curves and nine generic power curves to convert wind speed into power.

17 WIND ENERGY↗
DOE OSTI2022

Hamming distance and the onset of quantum criticality

Simulating models for quantum correlated matter unveils the inherent limitations of deterministic classical computations. In particular, in the case of quantum Monte Carlo methods, this is manifested by the emergence of negative weight configurations in the sampling, that is, the sign problem (SP). There have been several recent calculations which exploit the SP to locate underlying critical behavior. Here, utilizing a metric that quantifies phase-space ergodicity in such sampling, the Hamming distance, we suggest a significant advance on these ideas to extract the location of quantum critical points in various fermionic models, in spite of the presence of a severe SP. Combined with other methods, exact diagonalization in our case, it elucidates both the nature of the different phases as well as their location, as we demonstrate explicitly for the honeycomb and triangular Hubbard models, in both their U(1) and SU(2) forms. Furthermore, our approach exemplifies a possible path allowing the exploration of the phase diagram of a variety of fermionic quantum models hitherto considered to be impractical via quantum Monte Carlo simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Comparison of simulated neutrino emission models with data on Supernova 1987A

Here we compare models of supernova (SN) neutrino emission with the Kamiokande II data on SN 1987A using the Bayesian approach. These models are taken from simulations and are representative of current one-dimensional SN models. We find that models with a brief accretion phase of neutrino emission are the most favored. This result is not affected by varying the overall flux normalization or considering neutrino oscillations. We also check the compatibility of the best-fit models with the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The impact of coupled air–sea interaction on extreme East Asian summer monsoon simulation in CMIP5 models

Abstract In this study, the relationship between the ability to simulate air–sea interactions over the western North Pacific (WNP), and to reproduce the extreme East Asian summer monsoon (EASM), were investigated by comparing the performances of several global climate models (GCMs). High ranked in air–sea interaction simulation (HRA) and low ranked in air–sea interaction simulation (LRA) models were selected, according to their performance in simulating relations between sea surface temperature (SST) and precipitation over the WNP, from the ensemble of models that participated in the third and fifth phases of the Coupled Model Intercomparison Project (CMIP3, CMIP5). Compared with CMIP3 models, CMIP5 models exhibited improved simulations of the distinctive air–sea interaction over the WNP, namely, the strong atmospheric forcing on the ocean. Among CMIP5 models, HRA models, which reproduced intrinsic negative correlations between precipitation and SST over the WNP, could simulate the extreme EASM better than LRA models. In particular, HRA models generated a more realistic spatial distribution of the extreme EASM compared with LRA models. The defects of the LRA models resulted from distorted synoptic fields, including underestimated geopotential height and overestimated low‐level wind over the WNP, inducing unrealistic moisture supply and convection due to the exaggerated SST forcing. In contrast, reasonable air–sea interactions represented in HRA models lead to realistic synoptic fields over the WNP, and proper simulation of the extreme EASM.

Kim, Taehyung↗
DOE OSTI2026

Equivalent Point Sources for Simulating Nuclear/Radiological Search Operations

Modeling and simulation can be used to estimate the detection of different types and strengths of radioactive sources. Simulations of shielded gamma-ray sources in complex environments, such as searches in urban areas using vehicle-mounted detectors, can demonstrate what source strengths and shielding combinations should be able to be detected. For faster calculations, gamma-emitting sources that could have different shielding layers can be simulated ahead of time, resulting in equivalent point source emission distributions that can be used in larger simulation models. Larger-scale models will then run faster using equivalent point sources because the shielding layers around the sources have already been accounted for. Accompanying this paper are 260 equivalent point source emission files for 72 sources with different shielding layers, formatted for three different Monte Carlo codes. Users wanting to request versions of the SCALE and/or MCNP software can find details here: https://rsicc.ornl.gov/

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Calibration of the 2-Phase Bubble Tracking Model for Liquid Mercury Target Simulation with Machine Learning Surrogate Models

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory is one of the most powerful accelerator-driven neutron sources in the world. The intense protons strike on SNS's mercury target to provide bright neutron beams, which also leads to severe fluid-structure interactions inside the target. Prediction of resultant loading on the target is difficult particularly when helium gas is injected into mercury to reduce the loading and mitigate the pitting damage on vessel walls. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this multi-physics problem. However, several uncertain parameters in the R-P model require intensive simulations to determine their optimal values. With the help of machine learning and the measured target strain, we have studied the major uncertain parameters in this R-P model and developed a framework to identify optimal parameters that significantly reduce the discrepancy between simulations and experimental strains. The preliminary results show the possibility of using this mercury/helium mixture and surrogate models to predict a better match of target strain response when the helium gas is injected.

Lin, Lianshan↗