Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy systems modeling”

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

HybridSystemsSimulations.jl - Solving the Merchant Collocated Facilities with JuMP

The development of new clean-generation technologies also leads to new plant-level architectures that combine several generation and storage assets behind the point of connection. These co-located generation resources (Hybrid Systems) primarily operate as merchant assets that employ automated market bidding models and internal Energy Management Systems (EMS) to comply with the operator's signals. Formulating an optimal bidding model requires embedding the EMS control model into the bidding algorithm, resulting in a bi-level optimization problem. In this presentation, we first showcase using JuMP to formulate and solve this problem effectively for multiple merchant systems and the bidding outcomes considering different model formulations. Second, the bidding outcomes are later integrated into a PowerSimulations.jl (also built with JuMP) simulation to study the system-level effects of the various merchant bidding and the interactions between market-clearing models and the embedded EMS model. We will showcase simulations conducted in the RTS system considering different levels of merchant hybrid systems participation. The presentation provides the following specific insights on JuMP usage: 1) the Formulation of specialized bi-level problems with custom cuts to solve the merchant hybrid system bidding problem; 2) the integration of a modular model within a complex simulation workflow supported by JuMP in PowerSimulations.jl; 3) Accelerating the solution of power systems operations simulation that employ agent optimization problems using JuMP.

energy markets↗

Stochastic Optimization and Uncertainty Quantification of Natrium-based Nuclear-Renewable Energy Systems for Flexible Power Applications in Deregulated Markets

Rapid integration of variable renewable energy sources (VRES) has made modeling and stochastic optimization of hybrid energy systems crucial for studying their long-term performance and viability. However, most studies have focused on just historical data, which may be unreliable for capturing short-term fluctuations, rare events, and long-term patterns of energy demand, price, and the variability of renewable energy sources. For this study, optimal synthetic time series models were developed using Wasserstein distance. The models were validated by comparing the key statistical measures against those of the historical data. They were then used to optimize the integrated Natrium-style advanced energy systems and their long-term (30 years) economics. The stochastic model performs bi-level optimization to find the optimal sizes for the balance of plant and thermal energy storage, while also optimizing energy dispatch to achieve the maximum net present value. In studies of two deregulated markets (California ISO and the Electric Reliability Council of Texas), the integrated Natrium-style system performed better in CAISO than in ERCOT, given higher and more consistent electricity prices during peak-demand periods. The potentially enlarged cost associated with the variable operation and maintenance of the TES system also plays a significant role in driving the system sizing, thus its impacts on the system are investigated in detail through comparison against a baseline case. The study also finds that the bi-level optimization results based on stochastic gradient descent closely match the grid search results. The uncertainty quantification of the stochastic signals provides further NPV-related insights and probability distributions for the case studies. The normal standard error of the mean of NPV for the case with and without TES VOM for CAISO were found to be 7.73M (plus-minus sign) 1.09M USD and 104.99M (plus-minus sign) 1.25M USD, respectively based on a 95% confidence. Given the relatively small NPV variance based on 150 samples, the analysis affords the most robust possible prediction of the techno-economic performance of the integrated Natrium-style energy systems.

25 ENERGY STORAGE↗

Integration of Storage in the DC Link of a Full Converter-Based Distributed Wind Turbine

Energy storage is known to support the dispatchability of variable renewable resources. In this paper, we model a battery energy storage system (BESS) integrated with the DC link of a Type IV full converter-based wind turbine and the necessary controls to achieve efficient dispatch. To support the validation of control methodologies, we build a detailed model of a Type IV research wind turbine at the National Renewable Energy Laboratory (NREL), the Controls Advanced Research Turbine (CART 3), and we integrated a lithium-ion BESS model in grid-following mode into the model. The simulation results illustrate the sizing and control of the DC link-integrated BESS for a given variable wind resource and varying dispatch strategies (i.e., under constant, uncertain, and ramping wind scenarios). The integrated storage can smooth variabilities in distributed wind output, hedge against uncertainties, provide the ramping capability, as well as support stability under voltage and frequency transients. All of these have been illustrated in MATLAB/Simulink.

DC-link voltage↗

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗

Evaluating Impacts of Sustainable Aviation Fuel Production with CO2-to-Fuels Technologies on High Renewable Share Power Grid: Preprint

This paper investigates the impact of Sustainable Aviation Fuel (SAF) production using CO2-to-Fuels technologies on a future power grid with a high share of renewable energy. We focus on understanding the implications of the 2050 SAF production goal on the U.S. power system's long-term planning, encompassing generation, transmission, and cost analysis. Via the Regional Energy Deployment System (ReEDS) model, we developed a detailed SAF electricity demand model based on a low-temperature electrolysis-syngas fermentation-ethanol pathway. Four SAF target scenarios which aimto meet 10%, 15%, 20%, and 27% of SAF demand by 2050. These scenarios are exhaustively simulated to assess their impact on the power grid. Our results reveal that increasing SAF demand will result in higher electricity requirements, as well as expanded generator and transmission capacities, leading to an overall rise in system costs. However, these impacts are manageable within the broader context of U.S. capacity expansion plans. This study provides valuable insights into incorporating the CO2-to-Fuels electricity demand model and other carbon capture technologies into power system planning, emphasizing their significance in shaping a sustainable energy future.

BIOMASS FUELS,ENERGY PLANNING, POLICY, AND ECONOMY↗

MAM4xx

SAND2023-05377O The Modal Aerosol Model with 4 fixed modes (MAM4xx) is a performance-portable, C++ implementation for modeling the chemical reactions of aerosols in the atmosphere. MAM4xx is part of the Energy Exascale Earth System Model (E3SM). The E3SM project is an ongoing Earth system modeling, simulation, and prediction project that optimizes the use of DOE laboratory resources to meet the science needs of the nation and the mission needs of DOE. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Overfelt, James↗

Collaborative Research: Improved Efficiency and Coupling of the Radiation Code in the ACME Earth System Model. Final Report

This final report details all work performed on the project by both project partners. This project provided support to properly couple RTE+RRTMGP, a high-performance broadband radiation code, within DOE’s Energy Exascale Earth System Model (E3SM). RTE+RRTMGP is a successor to the RRTMG radiation code, which has been widely accepted for its speed and accuracy by the global modeling community, and has been in use in the NCAR CESM for many years and was implemented in the initial version of E3SM. However, the computational cost of RRTMG remains high relative to other components in part due to its complexity and to its inefficient use of modern optimization strategies, issues that were rectified by the development of RTE+RRTMGP. Many of the accomplishment in this project necessitated significant collaboration with the E3SM development team. One focus of the project was to enhance the code’s optimization on the limited number of emerging computing systems on which the model is expected be used, including Many Integrated Core (MIC) architectures and Graphics Processing Unit (GPU) hardware. We also developed additional capabilities for RTE+RRTMGP that E3SM scientists identified as important for the planned applications of the model. The result of our project was optimization of a key physical component (radiative transfer calculations) of E3SM, directly supporting E3SM’s overarching global modeling objectives. More broadly, this project provided overall advancements in the use of radiative transfer calculations in atmospheric modeling and simulation, particularly for climate.

54 ENVIRONMENTAL SCIENCES↗

The DOE E3SM Model Version 2: Overview of the Physical Model and Initial Model Evaluation

Abstract This work documents version two of the Department of Energy's Energy Exascale Earth System Model (E3SM). E3SMv2 is a significant evolution from its predecessor E3SMv1, resulting in a model that is nearly twice as fast and with a simulated climate that is improved in many metrics. We describe the physical climate model in its lower horizontal resolution configuration consisting of 110 km atmosphere, 165 km land, 0.5° river routing model, and an ocean and sea ice with mesh spacing varying between 60 km in the mid‐latitudes and 30 km at the equator and poles. The model performance is evaluated with Coupled Model Intercomparison Project Phase 6 Diagnosis, Evaluation, and Characterization of Klima simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate has many realistic features of the climate system, with notable improvements in clouds and precipitation compared to E3SMv1. E3SMv1 suffered from an excessively high equilibrium climate sensitivity (ECS) of 5.3 K. In E3SMv2, ECS is reduced to 4.0 K which is now within the plausible range based on a recent World Climate Research Program assessment. However, a number of important biases remain including a weak Atlantic Meridional Overturning Circulation, deficiencies in the characteristics and spectral distribution of tropical atmospheric variability, and a significant underestimation of the observed warming in the second half of the historical period. An analysis of single‐forcing simulations indicates that correcting the historical temperature bias would require a substantial reduction in the magnitude of the aerosol‐related forcing.

54 ENVIRONMENTAL SCIENCES↗

Energy-conserving physics for nonhydrostatic dynamics in mass coordinate models

Abstract. Motivated by reducing errors in the energy budget related to enthalpy fluxes within the Energy Exascale Earth System Model (E3SM), we study several physics–dynamics coupling approaches. Using idealized physics, a moist rising bubble test case, and the E3SM's nonhydrostatic dynamical core, we consider unapproximated and approximated thermodynamics applied at constant pressure or constant volume. With the standard dynamics and physics time-split implementation, we describe how the constant-pressure and constant-volume approaches use different mechanisms to transform physics tendencies into dynamical motion and show that only the constant-volume approach is consistent with the underlying equations. Using time step convergence studies, we show that the two approaches both converge but to slightly different solutions. We reproduce the large inconsistencies between the energy flux internal to the model and the energy flux of precipitation when using approximate thermodynamics, which can only be removed by considering variable latent heats, both when computing the latent heating from phase change and when applying this heating to update the temperature. Finally, we show that in the nonhydrostatic case, for physics applied at constant pressure, the general relation that enthalpy is locally conserved no longer holds. In this case, the conserved quantity is enthalpy plus an additional term proportional to the difference between hydrostatic pressure and full pressure.

58 GEOSCIENCES↗

Southern Ocean polynyas and dense water formation in a high-resolution, coupled Earth system model

Antarctic coastal polynyas produce dense shelf water, a primary source of Antarctic Bottom Water that contributes to the global overturning circulation. This paper investigates Antarctic dense water formation in the high-resolution version of the Energy Exascale Earth System Model (E3SM-HR). The model is able to reproduce the main Antarctic coastal polynyas, although the polynyas are smaller in area compared to observations. E3SM-HR also simulates several occurrences of open-ocean polynyas (OOPs) in the Weddell Sea at a higher rate than what the last 50 years of the satellite sea ice observational record suggests, but similarly to other high-resolution Earth system model simulations. Furthermore, the densest water masses in the model are formed within the OOPs rather than on the continental shelf as is typically observed. Biases related to the lack of dense water formation on the continental shelf are associated with overly strong atmospheric polar easterlies, which lead to a strong Antarctic Slope Front and too little exchange between on- and off-continental shelf water masses. Strong polar easterlies also produce excessive southward Ekman transport, causing a build-up of sea ice over the continental shelf and enhanced ice melting in the summer season. This, in turn, produces water masses on the continental shelf that are overly fresh and less dense relative to observations. Our results indicate that high resolution alone is insufficient for models to properly reproduce Antarctic dense water; the large-scale polar atmospheric circulation around Antarctica must also be accurately simulated.

54 ENVIRONMENTAL SCIENCES↗

Early Evaluation of Fugaku A64FX Architecture Using Climate Workloads

The Energy Exascale Earth System Model (E3SM) Project is an ongoing, state-of-the-science Earth system modeling, simulation, and prediction project that targets efficient utilization of U.S. Department of Energy’s (DOE) supercomputers to meet the science needs of the nation and the mission needs of DOE. This work focuses on our early evaluation of the A64FX architecture on Fugaku supercomputer using E3SM benchmarks. We will present results that track hardware trends, facilitate architecture comparison and the specific impact on our workload using an atmospheric model benchmark. We have two variants of the code written in Fortran and C++/Kokkos respectively which were used to collect data on a variety of CPU and GPU platforms. Furthermore, we have conducted a comparative evaluation of the compilers on the A64FX architecture and found GNU to be the best performer for our workload. Our experience so far indicates that Fugaku/A64FX shows promising energy efficiency (performance/Watt) with further performance gains possible through architecture-aware optimization efforts.

Sreepathi, Sarat↗

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model con figurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

58 GEOSCIENCES↗

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Abstract. Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model configurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

54 ENVIRONMENTAL SCIENCES↗

Analysis of Stationary, Photovoltaic-based Surface Power System Designs at the Lunar South Pole

Combinations of solar arrays and either batteries or regenerative fuel cells are analyzed for a surface power system module at the lunar south pole. The systems are required to produce 5 kW of net electrical power in sunlight and 2 kW of net electrical power during lunar night periods for a 10-year period between 2020 and 2030. Systems-level models for energy conservation, performance, degradation, and mass are used to compare to various systems. The sensitivities of important and/or uncertain variables including battery specific energy, fuel cell operating voltage, and DC-DC converter efficiency are compared to better understand the system. Switching unit efficiency, battery specific energy, and fuel cell operating voltage appear to be important system-level variables for this system. With reasonably sized solar arrays, the regenerative fuel cell system has significantly lower mass than the battery system based on the requirements and assumptions made herein. The total operational time is estimated at about 10,000 hours in battery discharge/fuel cell mode and about 4,000 and 8,000 hours for the battery charge and electrolyzer modes, respectively. The estimated number of significant depth-of-discharge cycles for either energy storage system is less than 100 for the 10-year period.

Freeh, Joshua E.↗

High quality superconducting Nb co-planar resonators on sapphire substrate

We present measurements and simulations of superconducting Nb co-planar waveguide resonators on sapphire substrate down to millikelvin temperature range with different readout powers. In the high temperature regime, we demonstrate that the Nb film residual surface resistance is comparable to that observed in the ultra-high quality, bulk Nb 3D superconducting radio frequency cavities while the resonator quality is dominated by the BCS thermally excited quasiparticles. At low temperature both the resonator quality factor and frequency can be well explained using the two-level system models. Through the energy participation ratio simulations, we find that the two-level system loss tangent is $\sim 10^{-2}$, which agrees quite well with similar studies performed on the Nb 3D cavities.

Zhu, S.↗

Effects of coupling a stochastic convective parameterization with the Zhang–McFarlane scheme on precipitation simulation in the DOE E3SMv1.0 atmosphere model

Abstract. A stochastic deep convection parameterization is implemented into the US Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1.0 (EAMv1). This study evaluates its performance in simulating precipitation. Compared to the default model, the probability distribution function (PDF) of rainfall intensity in the new simulation is greatly improved. The well-known problem of “too much light rain and too little heavy rain” is alleviated, especially over the tropics. As a result, the contribution from different rain rates to the total precipitation amount is shifted toward heavier rain. The less frequent occurrence of convection contributes to suppressed light rain, while more intense large-scale and convective precipitation contributes to enhanced heavy total rain. The synoptic and intraseasonal variabilities of precipitation are enhanced as well to be closer to observations. The sensitivity of the rainfall intensity PDF to the model vertical resolution is examined. The relationship between precipitation and dilute convective available potential energy in the stochastic simulation agrees better with that in the Atmospheric Radiation Measurement (ARM) observations compared with the standard model simulation. The annual mean precipitation is largely unchanged with the use of the stochastic scheme except over the tropical western Pacific, where a moderate increase in precipitation represents a slight improvement. The responses of precipitation and its extremes to climate warming are similar with or without the stochastic deep convection scheme.

54 ENVIRONMENTAL SCIENCES↗

Machine Learning Calibration of Groundwater Table Depth in ELM: Impact on Land Surface Hydrology and Land‐Atmosphere Fluxes

Accurate representation of groundwater table depth (GWTD) is crucial for simulating hydrological cycling in Earth system models (ESM). Nevertheless, there is a notable gap in the literature regarding the validation of GWTD simulations in ESMs and their subsequent impact on downstream hydrological components. This study explores the calibration of parameterization of global GWTD using machine learning within the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Despite achieving significant gains in simulating GWTD through calibration, offline ELM simulations unexpectedly show that these improvements do not translate to substantial enhancements in model performance for other key hydrological variables, including soil moisture (SM), runoff, groundwater contribution to runoff or base flow index (BFI), and evapotranspiration and its partitioning. The performance in SM and runoff was even degraded in some regions, while BFI was mostly overestimated. Although there is significant improvement in GWTD within the critical range of 1–5 m, where groundwater traditionally influences land surface energy fluxes, these improvements occurred mostly in humid areas where the impact of GWTD on surface processes is minimal. Although the impacts of model calibration are generally small in offline ELM simulations, coupled land-atmosphere simulations exhibit much stronger responses to GWTD calibration, highlighting the role of land-atmosphere feedbacks in Earth system modeling. These findings underscore the need for integrated calibration strategies that simultaneously optimize multiple hydrological variables. However, if a single-variable approach is necessary, it is crucial to establish clear priorities for calibration, identifying the most critical variables that have the greatest impact on overall model performance.

Fang, Yilin [Pacific Northwest National Laboratory↗

A Model for Hybrid Systems for Production Cost Modeling Studies Considering Ancillary Services: Preprint

This paper introduces a model for simulating hybrid plants participating in energy and ancillary services for bulk power system studies. The model considers a hybrid plant comprised of a renewable energy source, a thermal power unit, a storage unit, fixed power loads, or any combination of these technologies. The model focuses on Production Cost Modeling (PCM) studies under the assumption of centralized dispatch. We present an example case study to illustrate the use of the model in a single-stage production cost model similar to those conducted by planning agencies. We explore the allocation of behind-themeter ancillary services products and total energy participation to hybrid plant sub-assets and the resulting impacts on the system's ancillary service allocation. The model is implemented and simulated in a unit commitment problem in the RTS test system, which was modified to include a hybrid plant asset.

ancillary services↗