Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy systems model”

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 127 records · Page 7

PINNSTRIPES (Physics-Informed Neural Network SurrogaTe for Rapidly Identifying Parameters in Energy Systems) [SWR-22-12]

Energy systems models typically take the form of complex partial differential equations which make multiple forward calculations prohibitively expensive. Fast and data-efficient construction of surrogate models is of utmost importance for applications that require parameter exploration such as design optimization and Bayesian calibration. In presence of a large number of parameters, surrogate models that capture correct dependencies may be difficult to construct with traditional techniques. The issue is addressed here with the formulation of the surrogate model constructed via Physics-Informed Neural Networks (PINN) which capture the dependence with respect to the parameters to estimate, while using a limited amount of data. Since forward evaluations of the surrogate model are cheap, parameter exploration is made inexpensive, even when considering a large number of parameters.

Hassanaly, Malik↗

Summer Aerosol and Trace Gas Observations in Houston, Texas Using an Adaptable Mobile Facility

An aerosol container featuring a shared inlet system was deployed to Houston, Texas in July 2022, enabling direct, high-time-resolution in situ measurements of aerosols and trace gases. The internal rack system and floorplan was designed for adaptable modularity to elucidate aerosol physicochemical processes at fine scales. The design allowed for the deployment of a core instrument suite and additional customized research grade instruments. A heterogeneous mixture of aerosols was observed during three regimes: (1) intermittent black carbon (BC) and diurnal variations in aerosol chemical composition, (2) observed particle growth associated with SO 2 , (3) transported supermicron dust. The high variability of observed particles and gases in high time resolution indicated a complex urban area with multiple local and regional sources and processes. Particle growth rates of 7–16 nm/hr were observed for submicron particles during periods when SO 2 was >0.5 ppbv. Two periods of multi-day long-range transport events of dust from the African Sahara were observed in the supermicron and submicron particle modes with total mass concentrations up to 30 μg m −3 . Aerosol scattering angstrom exponents and extinction coefficients (B ext ) increased with humidity as a function of particle composition. The measurements demonstrate collaborative capabilities that can be used to increase observations of aerosol processing, microphysical and optical properties, internal mixing state, and supermicron aerosol that are not parameterized or missing in global Earth energy system models.

54 ENVIRONMENTAL SCIENCES↗

NucMesh: nuclear reactor geometry creation and mesh generation module in NEMoSys

NucMesh is a parameterized geometry and mesh generator for nuclear reactors developed within the Nuclear Energy Modeling System NEMoSys at Illinois Rocstar. NEMoSys is a platform developed for mesh generation, adaptive refinement, and solution verification. NucMesh is implemented to be generalized and extensible with a robust computer aided design engine and multiple mesh generation algorithms for unstructured triangular, quad-dominant, and structured quadrilateral meshing. In this paper, we present the geometric and meshing features of NucMesh. Geometrically objects are constructed bottom-to-top and overlaps are addressed automatically. A sophisticated object tracking algorithm prevents data from being lost for segmented objects. We discuss the primitive objects of circle and polygons that constitute the module and show how they are used with example inputs. Arrays of primitives and arrays of arrays are utilized to build large assemblies of objects. The concept of saved objects is discussed to demonstrate how repetitive objects can be reused easily and augmented in place. Three dimensional meshes can be obtained through mesh extrusion where all materials and side sets are extended to three dimensions. We show that side sets can be defined nearly anywhere within the geometry and can then be applied to the mesh. Finally, example reactor meshes are demonstrated for the Idaho National Laboratory Advanced Test Reactor and Los Alamos National Lab Empire reactor, both of which use control drums that NucMesh handles easily. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Quantum-Inspired Bayesian Sampling for Uncertainty Quantification and Machine Learning (Final Technical Report)

With increasing simulation and measurement data, machine learning and artificial intelligence have been widely used in computational decision-making of complex engineering systems. The resulting tools, such as uncertainty quantification solvers, reinforcement learning, and physics-informed machine learning, have achieved great success in critical DOE tasks such as material discovery and design, energy system modeling and control, and numerical weather and climate prediction. A core topic in scientific machine learning and artificial intelligence is Bayesian inference: given an observed data set, people want to estimate the posterior distribution of a (possibly large) number of hidden parameters. Due to the flexibility and weak assumptions, Bayesian sampling has been the mainstream Bayesian inference solvers despite the rapid progress of approximate Bayesian inference. Classical Bayesian sampling methods such as Markov-chain Monte Carlo suffer from a low-acceptance rate due to the random walk nature, therefore state-of-the-art techniques use Hamiltonian Monte Carlo and its variants to efficiently draw posterior samples in a high dimension. The key idea of Hamiltonian Monte Carlo and its variants is to simulate the Hamiltonian dynamics of a classical particle with a fixed mass, and their performance significantly degrades when the posterior distribution is highly spiky or has multiple modes. Leveraging the idea of quantum physics, this project has investigated new theory, algorithms and applications of Bayesian inference (especially Bayesian sampling). The main results include: (1) novel quantum-inspired Bayesian sampling methods that can lead to better accuracy for challenging multi-modal or spiky distributions, (2) more scalable machine learning framework leveraging tensor-compressed Bayesian inference, and (3) Bayesian and sampling approaches for verifying the robustness of continuous and binary neural networks.

97 MATHEMATICS AND COMPUTING↗

Assessing Two Approaches for Enhancing the Range of Simulated Scales in the E3SMv1 and the Impact on the Character of Hourly US Precipitation

Abstract Improving the representation of precipitation in Earth system models is essential for understanding and projecting water cycle changes across scales. Progress has been hampered by persistent deficiencies in representing precipitation frequency, intensity, and timing in current models. Here, we analyze simulated US precipitation in the low‐resolution (LR) configuration of the Energy Exascale Earth System Model (E3SMv1) and assess the effect of two approaches to enhance the range of explicitly resolved scales: high‐resolution (HR) and multiscale modeling framework (MMF), which incur similar computational expense. Both E3SMv1‐MMF and E3SMv1‐HR capture more intense and less frequent precipitation on hourly and daily timescales relative to E3SMv1‐LR. E3SMv1‐HR improves the intensity over the Eastern and Northwestern US during winter, while E3SMv1‐MMF improves the intensity over the Eastern US and summer diurnal timing over the Central US. These results indicate that both methods may be needed to improve simulations of different storm types, seasons, and regions.

58 GEOSCIENCES↗

Global evaluation of terrestrial biogeochemistry in the Energy Exascale Earth System Model (E3SM) and the role of the phosphorus cycle in the historical terrestrial carbon balance

Abstract. The importance of carbon (C)–nutrient interactions to the prediction of future C uptake has long been recognized. The Energy Exascale Earth System Model (E3SM) land model (ELM) version 1 is one of the few land surface models that include both N and P cycling and limitation (ELMv1-CNP). Here we provide a global-scale evaluation of ELMv1-CNP using the International Land Model Benchmarking (ILAMB) system. We show that ELMv1-CNP produces realistic estimates of present-day carbon pools and fluxes. Compared to simulations with optimal P availability, simulations with ELMv1-CNP produce better performance, particularly for simulated biomass, leaf area index (LAI), and global net C balance. We also show ELMv1-CNP-simulated N and P cycling is in good agreement with data-driven estimates. We compared the ELMv1-CNP-simulated response to CO2 enrichment with meta-analysis of observations from similar manipulation experiments. We show that ELMv1-CNP is able to capture the field-observed responses for photosynthesis, growth, and LAI. We investigated the role of P limitation in the historical balance and show that global C sources and sinks are significantly affected by P limitation, as the historical CO2 fertilization effect was reduced by 20 % and C emission due to land use and land cover change was 11 % lower when P limitation was considered. Our simulations suggest that the introduction of P cycle dynamics and C–N–P coupling will likely have substantial consequences for projections of future C uptake.

54 ENVIRONMENTAL SCIENCES↗

Subsurface hydrological controls on the short-term effects of hurricanes on nitrate–nitrogen runoff loading: a case study of Hurricane Ida using the Energy Exascale Earth System Model (E3SM) Land Model (v2.1)

When the nutrient level in the soil surpasses vegetation demand, nutrient losses due to surface runoff and subsurface leaching are the major reasons for the deterioration of water quality. The lower Mississippi River basin (LMRB) is one of the sub-basins that deliver the highest nitrogen loads to the Gulf of Mexico. Potential changes in episodic events induced by hurricanes may exacerbate water quality issue in the future. However, uncertainties in modeling the hydrologic response to hurricanes may limit the modeling of nutrient losses during such events. Using a machine learning approach, we calibrated the land component of the Energy Exascale Earth System Model (E3SM), or ELM, version 2.1, based on the water table depth (WTD) of a calibrated 3D subsurface hydrology model. While the overall performance of the calibrated ELM is satisfactory, some discrepancies in WTD remain in slope areas with low precipitation due to the missing lateral flow process in ELM. Simulations including biogeochemistry performed using ELM with and without model calibration showed important influences of soil hydrology, precipitation intensity, and runoff parameterization on the magnitude of nitrogen runoff loss and the leaching pathway. Despite such sensitivities, both ELM simulations produced reduced WTD and increased runoff and accelerated nitrate–nitrogen runoff loading during Hurricane Ida in August 2021, consistent with the observations. With observations suggesting more pronounced effects of Hurricane Ida on nitrogen runoff than the simulations, we identified factors for model improvement to provide a useful tool for studying hurricane-induced nutrient losses in the LMRB region.

54 ENVIRONMENTAL SCIENCES↗

An overview of cloud–radiation denial experiments for the Energy Exascale Earth System Model version 1

Abstract. The interaction between clouds and radiation is a key process within the climate system, and assessing the impacts of that interaction provides valuable insights into both the present-day climate and future projections. Many modeling experiments have been designed over the years to probe the impact of the cloud radiative effect (CRE) on the climate, including those that seek to disrupt the mean CRE effect and those that only disrupt the covariance of the CRE with the circulation. Seven such experimental designs have been added to the Energy Exascale Earth System Model version 1 (E3SMv1) of the US Department of Energy. These experiments include both the first and second iterations of the Clouds On/Off Klimate Intercomparison Experiment (COOKIE) experimental design, as well as the cloud-locking method. This paper documents the code changes necessary to implement such experiments and also provides detailed instructions for how to run them. Analyses across experiment types provide valuable insights and confirm the findings of prior studies, including the role of cloud radiative heating toward intensifying the monsoon, intensifying rain rates, and poleward expansion of the general circulation owing to cloud feedbacks.

Harrop, Bryce E. (ORCID:0000000339524525)↗

Regime-based aerosol–cloud interactions from CALIPSO-MODIS and the Energy Exascale Earth System Model version 2 (E3SMv2) over the Eastern North Atlantic

This study investigates aerosol-cloud interactions in marine boundary layer (MBL) clouds using an advanced deep-learning-driven synoptic-regime-based framework, combining satellite data (CALIPSO vertically resolved aerosol extinction and MODIS cloud properties) with 1° nudged Energy Exascale Earth System Model version 2 (E3SMv2) simulation over the Eastern North Atlantic (ENA; ∼10°×10°, 2006–2014). The E3SMv2 captures observed seasonal variations in cloud droplet number concentrations (N d ) and liquid water path (LWP), though it systematically underestimates N d . We then partition ENA meteorology into four synoptic regimes (Pre-Trough, Post-Trough, Ridge, Trough) via a deep-learning clustering of ERA5 reanalysis fields, enabling regime-dependent aerosol-cloud interactions analyses. Both satellite and E3SMv2 exhibit an inverted-V LWP-N d relationship. In Post-Trough and Ridge regimes, the satellite shows stronger negative LWP-N d sensitivities than in Pre-Trough regime. The Trough regime displays a muted satellite LWP response. In comparison, the model predicts more exaggerated LWP responses across regimes, with LWP increasing too quickly at low N d and decreasing more sharply at high N d , especially in Pre-Trough and Trough regimes. These exaggerated model LWP sensitivities may stem from uncertainties in representing drizzle processes, entrainment, and turbulent mixing. As for N d susceptibility to aerosols, N d increases with MBL aerosol extinction in both datasets, but the simulated aerosol-cloud interactions appear oversensitive to meteorological conditions. Overall, E3SMv2 better captures aerosol effects under regimes that favor stratiform clouds (Post-Trough, Ridge), but performance deteriorates for regimes with deeper, dynamically complex clouds (Trough), highlighting the need for improved representations of those cloud processes in climate models.

Environmental sciences↗

Diurnal rainfall response to the physiological and radiative effects of CO 2 in tropical forests in the Energy exascale earth system model v1

Understanding how the connection between rainfall and tropical forests will respond to increasing CO$_2$ concentrations is a key element in understanding how the tropical water cycle will respond to increasing CO$_2$. The plant physiological and radiative impacts of CO$_2$ on rainfall patterns over tropical forest regions are examined in the Energy Exascale Earth System Model version 1.1 (E3SMv1.1-BGC) biogeochemistry experiments. Composite analysis reveals a dampening of the diurnal cycle of rainfall over the Amazon, Congo, and Maritime Continent in response to rising CO$_2$ levels, regardless of the sign of total rainfall change. A full factorial model experiment confirms that the CO$_2$ radiative and CO$_2$ plant physiological effects can individually or jointly reduce the magnitude of the rainfall diurnal cycle, though the physical pathway giving rise to the reduction differs between the two effects. For the physiological response, stomatal closure reduces evapotranspiration, which dries the boundary layer and raises the lifting condensation level. These effects combine to reduce deep convective rainfall during its peak occurrence in the late daytime to early nighttime period. For the radiative response, a relative reduction in daytime CAPE (consistent with a reduction in the diurnal temperature range) leads to less frequent triggering of deep convection and a reduction of rainfall diurnal amplitude. These diurnal rainfall changes are structurally similar across seasons, and show little sensitivity to representation of nutrient coupling for the land biogeochemistry. In agreement with previous findings, the physiological response has only minor impact on extreme rainfall relative to the radiative response.

54 ENVIRONMENTAL SCIENCES↗

Subseasonal Tropical Convection Characteristics in the Energy Exascale Earth System Model Version 2

Accurate simulation of subseasonal tropical moist convection remains a key challenge for Earth system models. The difficulties stem from the reliance of cumulus cloud processes on model parameterizations and the need to represent the multiscale nature of interactions among clouds, radiation, moisture, circulation, and surface energy fluxes. Equatorial convection drives circulation anomalies that can affect weather patterns and extremes globally, motivating efforts to better understand and simulate these tropical disturbances. Here, a detailed review of subseasonal tropical convective behavior as simulated in the Energy Exascale Earth System Model version 2 (E3SMv2) is presented, with comparison to its predecessor version 1 (E3SMv1) and reference data sets. Model structural changes to the deep convective trigger function and surface fluxes, along with parametric tuning of the cloud and microphysics schemes, together result in an improved depiction of organized tropical convection across scales. In particular, E3SMv2 exhibits a more realistic Madden‐Julian oscillation (MJO) and low‐frequency Kelvin waves—owing to a sharper time mean equatorial meridional moisture gradient and improved convection‐circulation coupling —as well as a better depiction of MJO Northern Hemisphere teleconnections. Despite these improvements, subseasonal precipitation variance continues to be strongly underestimated in E3SMv2. Use of a cloud plume model also reveals that the coupling between daily averaged tropical precipitation and lower tropospheric instability in E3SM is inconsistent with observations, a bias that could potentially impact the simulation of intraseasonal disturbances.

54 ENVIRONMENTAL SCIENCES↗

More Realistic Intermediate Depth Dry Firn Densification in the Energy Exascale Earth System Model (E3SM)

Abstract Earth system models account for seasonal snow cover, but many do not accommodate the deeper snowpack on ice sheets (aka firn) that slowly transforms to ice under accumulating snowfall. To accommodate and resolve firn depths of up to 60 m in the Energy Exascale Earth System Model's land surface model (ELM), we add 11 layers to its snowpack and evaluate three dry snow compaction equations in multi‐century simulations. After comparing results from ELM simulations (forced with atmospheric reanalysis) with empirical data, we find that implementing into ELM a two‐stage firn densification model produces more accurate dry firn densities at intermediate depths of 20–60 m. Compared to modeling firn using the equations in the (12 layer) Community Land Model (version 5), switching to the two‐stage firn densification model (with 16 layers) significantly decreases root‐mean‐square errors in upper 60 m dry firn densities by an average of 41 kg m −3 (31%). Simulations with three different firn density parameterizations show that the two‐stage firn densification model should be used for applications that prioritize accurate upper 60 m firn air content (FAC) in regions where the mean annual surface temperature is greater than roughly −31°C. Because snow metamorphism, firn density, and FAC are major components in modeling ice sheet surface albedo, melt water retention, and climatic mass balance, these developments advance broader efforts to simulate the response of land ice to atmospheric forcing in Earth system models.

58 GEOSCIENCES↗

Future Changes in Active and Inactive Atlantic Hurricane Seasons in the Energy Exascale Earth System Model

Abstract North Atlantic tropical cyclones (TCs) have considerable interannual variability, with La Niña and the positive phase of the Atlantic Meridional Mode (AMM) tending to drive active hurricane seasons, and El Niño and the negative AMM often driving inactive seasons. Here, we analyze how active and inactive Atlantic hurricane seasons may change in the future using the high resolution Energy Exascale Earth System Model (E3SM). We performed atmosphere‐only simulations forced by sea‐surface temperature patterns characteristic of La Niña and the positive AMM jointly, and El Niño and the negative AMM jointly, in historical and future climates. Projected Atlantic TCs become more frequent in the future by approximately 34% during El Niño and negative AMM and by 66% during La Niña and positive AMM, with a significant increase in the portion of intense TCs. Warmer SSTs increase TC potential intensity, with reduced wind shear and increased mid‐tropospheric humidity further supporting TC activity.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Water Cycle Over CONUS at the Watershed Scale for the Energy Exascale Earth System Model Version 1 (E3SMv1) Across Resolutions

Abstract The water cycle is an important component of the earth system and it plays a key role in many facets of society, including energy production, agriculture, and human health and safety. In this study, the Energy Exascale Earth System Model version 1 (E3SMv1) is run with low‐resolution (roughly 110 km) and high‐resolution (roughly 25 km) configurations—as established by the High Resolution Model Intercomparison Project protocol—to evaluate the atmospheric and terrestrial water budgets over the conterminous United States (CONUS) at the large watershed scale. The warm season water cycle slows down in the HR experiment relative to the LR, with decreasing fluxes of precipitation, evapotranspiration, atmospheric moisture convergence, and runoff. The reductions in these terms exacerbate biases for some watersheds, while reducing them in others. For example, precipitation biases are exacerbated at HR over the Eastern and Central CONUS watersheds, while precipitation biases are reduced at HR over the Western CONUS watersheds. The most pronounced changes with resolution to the water cycle come from reductions in precipitation and evapotranspiration. The reduction in evapotranspiration reduces the biases across nearly all of the CONUS. Additional exploratory metrics show improvements to water cycle extremes (both in precipitation and streamflow), fractional contributions of different storm types to total precipitation, and mountain snowpack.

54 ENVIRONMENTAL SCIENCES↗

Using a surrogate-assisted Bayesian framework to calibrate the runoff-generation scheme in the Energy Exascale Earth System Model (E3SM) v1

Abstract. Runoff is a critical component of the terrestrial water cycle, and Earth system models (ESMs) are essential tools to study its spatiotemporal variability. Runoff schemes in ESMs typically include many parameters so that model calibration is necessary to improve the accuracy of simulated runoff. However, runoff calibration at a global scale is challenging because of the high computational cost and the lack of reliable observational datasets. In this study, we calibrated 11 runoff relevant parameters in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) using a surrogate-assisted Bayesian framework. First, the polynomial chaos expansion machinery with Bayesian compressed sensing is used to construct computationally inexpensive surrogate models for ELM-simulated runoff at 0.5∘ × 0.5∘ for 1991–2010. The error metric between the ELM simulations and the benchmark data is selected to construct the surrogates, which facilitates efficient calibration and avoids the more conventional, but challenging, construction of high-dimensional surrogates for the ELM simulated runoff. Second, the Sobol' index sensitivity analysis is performed using the surrogate models to identify the most sensitive parameters, and our results show that, in most regions, ELM-simulated runoff is strongly sensitive to 3 of the 11 uncertain parameters. Third, a Bayesian method is used to infer the optimal values of the most sensitive parameters using an observation-based global runoff dataset as the benchmark. Our results show that model performance is significantly improved with the inferred parameter values. Although the parametric uncertainty of simulated runoff is reduced after the parameter inference, it remains comparable to the multimodel ensemble uncertainty represented by the global hydrological models in ISMIP2a. Additionally, the annual global runoff trend during the simulation period is not well constrained by the inferred parameter values, suggesting the importance of including parametric uncertainty in future runoff projections.

58 GEOSCIENCES↗

Status Report on FY2022 Model Development within the Integrated Energy Systems HYBRID Repository

This publication details newly created energy storage models developed within the HYBRID Modelica repository as part of the Department of Energy Office of Nuclear Energy (DOE-NE) Integrated Energy Systems (IES) program, led by Idaho National Laboratory (INL). Model development to-date includes creation of dynamic systems-level models of a pebble bed high temperature gas reactor (HTGR), liquid air energy storage (LAES), and compressed air energy storage (CAES) in the IES-based HYBRID repository. Models are developed using the latest publicly available data and incorporate the possibility of control strategy inclusion for use with the existing IES modeling, analysis, and optimization toolset. Simulations showcase the abilities of each technology to flexibly operate in ways consistent with IES operation expectations. When these models are available, they can be utilized within different integrated energy park concepts to understand optimal system operation, control, and dispatching. Moreover, given the generic nature of the models, industrial partner technologies can be quickly added to the repository using the existing models as a basis. Additional dynamic models for thermal energy storage concepts can be developed and added to the HYBRID repository as needed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An Assessment of Nonhydrostatic and Hydrostatic Dynamical Cores at Seasonal Time Scales in the Energy Exascale Earth System Model (E3SM)

Abstract In global atmospheric modeling, the differences between nonhydrostatic (NH) and hydrostatic (H) dynamical cores are negligible in dry simulations when grid spacing is larger than 10 km. However, recent studies suggest that those differences can be significant at far coarser resolution when moisture is included. To better understand how NH and H differences manifest in global fields, we perform and analyze an ensemble of 28 and 13 km seasonal simulations with the NH and H dynamical cores in the Energy Exascale Earth System Model global atmosphere model, where the differences between H and NH configurations are minimized. A set of idealized rising bubble experiments is also conducted to further investigate the differences. Although NH and H differences are not significant in global statistics and zonal averages, significant differences in precipitation amount and patterns are observed in parts of the tropics. The most prominent differences emerge near India and the Western Pacific in the boreal summer, and the central‐southern Indian Ocean and Pacific in the boreal winter. Tropical differences influence surrounding regions through modification of the regional circulation and can propagate to the extratropics, leading to significant temperature and geopotential differences over the middle to high latitudes. While the dry bubble experiments show negligible deviation between H and NH dynamics until grid spacing is below 6.25 km, precipitation amount and vertical velocity are different in the moist case even at 25 km resolution.

58 GEOSCIENCES↗

Simulating Extreme Precipitation in the United States in the Energy Exascale Earth System Model: Investigating the Importance of Representing Convective Intensity versus Dynamic Structure (Final Report)

Weather events that produce extreme precipitation are associated with severe flooding and winds that result in thousands of deaths and billions of dollars in damages annually in the United States (U.S.). This project aims to improve understanding of the small- and large-scale processes that govern these events and improve our ability to project changes in these extreme events under the influences of natural variability and human activities. We focus on two promising directions in the development of the U.S. Department of Energy’s Energy Exascale Earth System Model (E3SM) to investigate the tradeoffs between resolving the convective-scale processes that control the intensity versus the intermediate-scale processes that control the dynamic structure of the most prominent extreme precipitation events that impact the U.S. throughout the year (i.e., mesoscale convective systems, tropical cyclones, and extratropical cyclone).

58 GEOSCIENCES↗