Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large system 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A complementarity-based vapor-liquid equilibrium formulation for equation-oriented simulation and optimization

Vapor-Liquid Equilibrium (VLE) is a cornerstone of computer aided process engineering (CAPE). Embedded within process system models, VLE calculations are inherently procedural with non-smooth behavior that frequently require discrete decisions. Traditionally, these features resist the incorporation of VLE within efficient, large-scale equation-oriented (EO) process simulation and optimization strategies. On the other hand, recent reformulation of VLE models through the incorporation of complementarity constraints has broadened its scope to deal seamlessly with phase transitions and even supercritical excursions in process simulation and optimization. In this study, we extend these VLE complementarity models to EO frameworks where procedural thermodynamic property libraries are still required. Here, we develop an efficient, non-intrusive, and intuitive “square-flash” equation system that has been implemented within the IDAES Integrated Platform (IDAES-IP). Further, the effectiveness of this modular approach is demonstrated on case studies for non-ideal flash calculations and distillation optimization with disappearing phases and supercritical transitions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. In this study, two-level generalized Dryja--Smith--Widlund (GDSW)--type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the message passing interface (MPI)--parallel implementation of multilevel Schwarz preconditioners provided by the package FROSch (fast and robust Schwarz) from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To the best of our knowledge, this is the first time two-level Schwarz preconditioners have been applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The preconditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the subdomains. Additionally, several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as nonuniform meshes for the Greenland ice sheet. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32,K processor cores (8,K MPI ranks and 4 OpenMP threads) and 566,M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4,K processor cores (and MPI ranks) and 68,M degrees of freedom for the coupled problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. Here, two-level GDSW (Generalized Dryja–Smith–Widlund) type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the MPI-parallel implementation of multi-level Schwarz preconditioners provided by the package FROSch (Fast and Robust Schwarz)from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To our knowledge, this is the first time two-level Schwarz preconditioners are applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The pre-conditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the sub-domains. Several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored as well. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as non uniform meshes for the Greenland ice sheet are considered. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32 K processor cores (8 K MPI-ranks and 4 OpenMP threads) and 566 M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4 K processor cores (and MPI-ranks) and 68 M degrees of freedom for the coupled problem.

58 GEOSCIENCES↗
DOE OSTI2025

Hardware-Efficient Quantum Phase Estimation via Local Control

Quantum phase estimation plays a central role in quantum simulation as it enables the study of spectral properties of many-body quantum systems. Most variants of the phase estimation algorithm require the application of the global unitary evolution conditioned on the state of one or more auxiliary qubits, posing a significant challenge for current quantum devices. In this work, we present an approach to quantum phase estimation that uses only locally controlled operations, resulting in a significantly reduced circuit depth. At the heart of our approach are efficient routines to measure the complex phase of the expectation value of the time-evolution operator, the so-called Loschmidt echo, for both circuit dynamics and Hamiltonian dynamics. By tracking changes in the phase during the dynamics, the routines trade circuit depth for increased sampling cost and classical postprocessing. Our approach does not rely on reference states and is applicable to any efficiently preparable state, regardless of its correlations. We provide a comprehensive analysis of the sample complexity and illustrate the results with numerical simulations. Our methods offer a practical pathway for measuring spectral properties in large many-body quantum systems using current quantum devices.

Schiffer, Benjamin F. [Max Planck Institute of Qua↗
DOE OSTI2020

A global model–measurement evaluation of particle light scattering coefficients at elevated relative humidity

The uptake of water by atmospheric aerosols has a pronounced effect on particle light scattering properties, which in turn are strongly dependent on the ambient relative humidity (RH). Earth system models need to account for the aerosol water uptake and its influence on light scattering in order to properly capture the overall radiative effects of aerosols. Here we present a comprehensive model–measurement evaluation of the particle light scattering enhancement factor f(RH), defined as the particle light scattering coefficient at elevated RH (here set to 85 %) divided by its dry value. The comparison uses simulations from 10 Earth system models and a global dataset of surface-based in situ measurements. In general, we find a large diversity in the magnitude of predicted f(RH) amongst the different models, which can not be explained by the site types. Based on our evaluation of sea salt scattering enhancement and simulated organic mass fraction, there is a strong indication that differences in the model parameterizations of hygroscopicity and model chemistry are driving at least some of the observed diversity in simulated f(RH). Additionally, a key point is that defining dry conditions is difficult from an observational point of view and, depending on the aerosol, may influence the measured f(RH). The definition of dry also impacts our model evaluation, because several models exhibit significant water uptake between RH = 0 % and 40 %. The multisite average ratio between model outputs and measurements is 1.64 when RH = 0 % is assumed as the model dry RH and 1.16 when RH = 40 % is the model dry RH value. The overestimation by the models is believed to originate from the hygroscopicity parameterizations at the lower RH range which may not implement all phenomena taking place (i.e., not fully dried particles and hysteresis effects). This will be particularly relevant when a location is dominated by a deliquescent aerosol such as sea salt. Our results emphasize the need to consider the measurement conditions in such comparisons and recognize that measurements referred to as dry may not be dry in model terms. Recommendations for future model–measurement evaluation and model improvements are provided.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Improving thermodynamic nudging in the E3SM Atmosphere Model version 2 (EAMv2): strategy and hindcast skills on weather systems

Nudging techniques are commonly employed to constrain atmospheric simulations toward observed states, facilitating model evaluation and sensitivity studies. However, if applied improperly – particularly to thermodynamic variables such as temperature and humidity – nudging can distort physical processes and introduce spurious biases, undermining the credibility of the simulations. This study presents an improved nudging implementation that applies vertically modulated tendencies to reduce adverse impacts on model physics. The framework is tested in version 2 of the Energy Exascale Earth System Model (EAMv2) using a suite of hindcast simulations nudged toward ERA5 reanalysis. We systematically evaluate the individual and combined effects of nudging wind, temperature, and humidity fields on the model's ability to represent large-scale atmospheric states and high-impact weather systems. Results show that the revised strategy – particularly when nudging temperature and humidity at selected levels – enhances hindcast skill by improving agreement with ERA5 without degrading the hydrological cycle or precipitation processes. Additional improvements in surface temperature, outgoing longwave radiation, and precipitation biases are achieved through targeted nudging of land-surface variables. The proposed approach strengthens the representation of large-scale conditions relevant to tropical cyclones, atmospheric rivers, and extratropical cyclones in the low-resolution EAMv2. These findings demonstrate that carefully designed thermodynamic nudging, especially of temperature and humidity, improves the realism of constrained simulations and broadens the utility of nudged EAMv2 for atmospheric modeling, machine learning, and high-impact weather research.

Atmospheric river↗
DOE OSTI2025

Computationally efficient models for aqueous organic redox flow batteries

The rising usage of intermittent energy has garnered the need for large scale energy storage systems. Redox flow batteries (RFB) based energy storage system shows promising potential. Numerical simulations and machine learning approaches have been widely used to study RFB performance. The development of autonomous material discovery framework and digital twin of energy storage system usually needs to query cell performance through fast response models. In this study, two computationally efficient models are introduced: a physics-based analytical flow battery model (EZBattery), and a machine learning operator model (Deep Operator Network, denoted by DeepONet). Both models can provide cell performance near instantly, and prediction accuracy was systematically examined on an application of evaluating the performances of a 780 cm 2 aqueous organic redox flow battery (AORFB), using potential anolyte candidates in dihydroxyphenazine (DHP)-based family of organic materials. A validated computationally expansive 3-dimensional multi-physics finite element model by COMSOL was used as the ground truth and provided the training data set for the DeepONet. 1280 samples were generated with 10 properties to mimic the different possible anolyte candidates, and the cell performances were evaluated under 10 different combined operating conditions. The accuracy comparisons for the two computationally efficient models show that both models can provide comparable accuracy in predicting cell charging/discharging voltage curves. DeepONet can provide slightly higher overall accuracy than EZBattery with faster calculation speed, but highly relies on the training dataset. EZBattery does not need a training dataset and can provide interpretable physics-based explanations of the results, while being more flexible to adjust to adapt any different cell designs, flow battery architectures, and electrolyte materials.

Analytical model↗
DOE OSTI2021

Decentralized Data-Driven Estimation of Generator Rotor Speed and Inertia Constant Based on Adaptive Unscented Kalman Filter

This paper proposes an online decentralized and data-driven method for synchronous generator rotor speed and inertia estimation. First, an analytical relationship between the generator terminal voltage/current phasors and its internal rotor speed is derived using the Thevenin equivalent. The parameters of the latter are obtained via a recursive estimator, which allows us to estimate the generator's internal rotor speed from voltage/current phasors without the need for any generator information. Second, we reformulate the equivalent swing equation into the Kalman filter form, which enables us to use a priori-estimated internal rotor speed along with measured real and reactive power injections to obtain an estimate of the inertia for each online synchronous generator. This is achieved through an adaptive unscented Kalman filter. We demonstrate that the use of terminal bus frequency to approximate the rotor speed for inertia estimation by existing methods leads to large estimation biases/errors. Numerical simulations carried out on the IEEE 39-bus and Texas 2000-bus systems reveal that our proposed method is consistently more accurate than existing methods. Moreover, our proposed method is insensitive to both the type and location of system disturbances, and it is robust to different levels of measurement noise.

bus frequency↗
DOE OSTI2021

Stochastic pre-event preparation for enhancing resilience of distribution systems

Extreme weather events are the common causes for power supply interruptions and power outages in electrical distribution systems. Improving the distribution system and enhancing its resilience is becoming crucial due to the increased frequency of extreme weather events. Preparation and allocation of multiple flexible resources, such as mobile resources, fuel resources, and labor resources before extreme weather events can mitigate the effects of extreme weather events and enhance the resilience of power distribution systems. Here, in this paper, a two-stage stochastic mixed-integer linear programming (SMILP) is proposed to optimize the preparation and resource allocation process for upcoming extreme weather events, which leads to faster and more efficient post-event restoration. The objective of the proposed two-stage SMILP is to maximize the served load and minimize the operating cost of flexible resources. The first stage in the optimization problem selects the amounts and locations of different resources. The second stage considers the operational constraints of the distribution system and repair crew scheduling constraints. The proposed stochastic pre-event preparation model is solved by a scenario decomposition method, Progressive Hedging (PH), to ease the computational complexity introduced by a large number of scenarios. Furthermore, to show the impact of solar photovoltaic (PV) generation on system resilience, three types of PV systems are considered during a power outage and the resilience improvements with different PV penetration levels are compared. Numerical results from simulations on a large-scale (more than 10,000 nodes) distribution feeder have been used to validate the effectiveness and scalability of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

A joint role for forced and internally-driven variability in the decadal modulation of global warming

Despite the observed monotonic increase in greenhouse-gas concentrations, global mean temperature displays important decadal fluctuations typically attributed to both external forcing and internal variability. Here, we provide a robust quantification of the relative contributions of anthropogenic, natural, and internally-driven decadal variability of global mean sea surface temperature (GMSST) by using a unique dataset consisting of 30-member large initial-condition ensembles with five Earth System Models (ESM-LE). We present evidence that a large fraction (~29–53%) of the simulated decadal-scale variance in individual timeseries of GMSST over 1950–2010 is externally forced and largely linked to the representation of volcanic aerosols. Comparison with the future (2010–2070) period suggests that external forcing provides a source of additional decadal-scale variability in the historical period. Given the unpredictable nature of future volcanic aerosol forcing, it is suggested that a large portion of decadal GMSST variability might not be predictable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Estimating Lossy Compressibility of Scientific Data Using Deep Neural Networks

Simulation based scientific applications generate increasingly large amounts of data on high-performance computing (HPC) systems. To allow data to be stored and analyzed efficiently, data compression is often utilized to reduce the volume and velocity of data. However, a question often raised by domain scientists is the level of compression that can be expected so that they can make more informed decisions, balancing between accuracy and performance. In this letter, we propose a deep neural network based approach for estimating the compressibility of scientific data. To train the neural network, we build both general features as well as compressor-specific features so that the characteristics of both data and lossy compressors are captured in training. Our approach is demonstrated to outperform a prior analytical model as well as a sampling based approach in the case of a biased estimation, i.e., for SZ. However, for the unbiased estimation (i.e., ZFP), the sampling based approach yields the best accuracy, despite the high overhead involved in sampling the target dataset.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A machine learning-based fast frequency response control for a VSC-HVDC system

An HVDC system can realize a very fast frequency response to the disturbed system under a contingency because its active power control is decoupled from the frequency deviation. However, most of existing HVDC frequency control strategies are coupled with system primary frequency control and secondary frequency control. Since the traditional system frequency control is dominated by the thermal generators, the advantage of the fast response of the HVDC system is not made fully used. The development of a frequency response estimation based on a machine learning algorithm provides another approach to improve the frequency response capability of the HVDC system. Different from other frequency deviation tracking strategies, a machine learning based HVDC frequency response control can directly increase the power flow of a HVDC system by estimation of the system generator or load lost. In this paper, a fast frequency response control using a HVDC system for a large power system disturbance based on the multivariate random forest regression (MRFR) algorithm is proposed. The simulation is carried out with an integrated power system model based on the North American interconnections. The simulation results indicate that the proposed MRFR based frequency response control can significantly improve the frequency low point during an event, while stabilizing the frequency in advance.

42 ENGINEERING↗
DOE OSTI2020

Design and Analysis of a Halo-Measurement Diagnostics

A large dynamical-range diagnostics (LDRD) design at Jefferson Lab will be used at the FAST-IOTA injector to measure the transverse distribution of halo associated with a high-charge electron beam. One important aspect of this work is to explore the halo distribution when the beam has significant angular momentum (i.e. is magnetized). The beam distribution is measured by recording radiation produced as the beam impinges a YAG:Ce screen. The optical radiation is split with a fraction directed to a charged-couple device (CCD) camera. The other part of the radiation is reflected by a digital micromirror device (DMD) that masks the core of the beam distribution. Combining the images recorded by the two cameras provides a measurement of the transverse distribution with over a large dynamical range. The design and analysis of the optical system will be discussed including optical simulation using SRW and the result of a mockup experiment to test the performances of the system will be presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Advancing Understanding of Deep Convective Anvil Clouds

In this project we use ground-based DOE ARM observations, satellite observations, and Weather Research and Forecasting (WRF) model simulations to arrive at a more complete understanding (or, new theory) of how atmospheric heating profiles in thunderstorm cores drives how quickly air rises in the storms. The rate at which air motion increases or decreases with altitude directly relates to how much cloud mass is ejected laterally into the adjacent atmosphere, which then influences the development and maintenance of “stratiform raining cloud anvils”. Our work has led to new equations connecting the heating in these thunderstorms to these stratiform clouds, with the connection modulated by the stability structure of the atmosphere (i.e., if more stable, then a more rapid change in thunderstorm heating with altitude results in more stratiform cloud, independent of whether over land or ocean, or how organized the thunderstorms are). Improving knowledge of stratiform clouds is important, since they contribute to substantial accumulated rainfall (important to agriculture and local hydrology) as well as radiation, given their large areal extents and lifetime. The radiation impacts likely impact simulated cloud feedback in all Earth System models (ESMs).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Large-scale offshore wind farm effects on weather and climate in Puerto Rico (Final Technical Report)

Puerto Rico’s current electricity generation heavily relies on imported fossil fuels. This results in an average cost of electricity higher than the U.S. mainland average in all sectors (residential, commercial, and industrial), despite abundant local offshore wind resources, which have the potential to provide secure, low-cost energy generation and consequent economic prosperity. However, effects on atmospheric and oceanic circulation resulting from large-scale deployments of offshore wind farms have not been previously studied at tropical latitude. This project addressed this knowledge gap through a computational modeling effort designed to capture the coupled dynamics of the atmosphere and the ocean in presence of offshore wind farms. Results indicate that wind farm wakes can alter wind stress, generate Ekman-driven vertical transport, and potentially affect nutrient distribution. While full model coupling remains challenging, progress in parameterization and large-eddy simulations provides a foundation for future research. The project contributes to DOE’s Earth System modeling efforts and supports STEM workforce development.

17 WIND ENERGY↗
DOE OSTI2026

Evaluation of UKESM aerosol size and composition using ATom measurements indicates missing marine aerosol formation mechanisms

Atmospheric aerosols influence climate through their interactions with radiation and clouds, yet large uncertainties remain in their simulation by global models. This study evaluates the United Kingdom Earth System Model version 1.1 (UKESM1.1) using global-scale aircraft observations from the Atmospheric Tomography (ATom) mission, focusing on aerosol lifecycle processes in the remote marine atmosphere. We assess model performance in simulating aerosol precursor vapours, number size distributions, chemical composition, and environmental conditions. Several process improvements are tested, including sulfuric acid-ammonia nucleation, ammonium nitrate scheme, methanesulfonic acid condensation, and low-temperature isoprene-derived secondary organic aerosol formation. Model biases differ significantly between the upper troposphere (UT) and the marine boundary layer (MBL). In the UT, UKESM1.1 overestimates nucleation and Aitken mode particles while underestimating accumulation mode, indicating insufficient growth. In the MBL, the model overestimates primary aerosols (e.g. seasalt) and precursor gases but underestimates nucleation and Aitken mode particles, even after incorporating updated nucleation and ammonium nitrate scheme. The persistence of low aerosol number concentrations, despite overestimated precursors, suggests missing formation pathways likely involving other species such as iodine, amines, and organic vapours. These limitations result in an unbalanced cloud condensation nuclei budget that over-relies on primary emissions. Sensitivity tests reveal that model outputs are strongly influenced by dimethyl sulfide emissions and vapour condensation schemes. Our results highlight the need for future model development to prioritise mechanistic representation of currently missing aerosol sources, rather than relying on empirical tuning, to improve aerosol-climate interaction estimates.

He, Xu-Cheng [Univ. of Cambridge (United Kingdom);↗
DOE OSTI2024

Hemispherically symmetric strategies for stratospheric aerosol injection

Stratospheric aerosol injection (SAI) comes with a wide range of possible design choices, such as the location and timing of the injection. Different stratospheric aerosol injection strategies can yield different climate responses; therefore, understanding the range of possible climate outcomes is crucial to making informed future decisions on SAI, along with the consideration of other factors. Yet, to date, there has been no systematic exploration of a broad range of SAI strategies. This limits the ability to determine which effects are robust across different strategies and which depend on specific injection choices. This study systematically explores how the choice of SAI strategy affects climate responses in one climate model. Here, we introduce four hemispherically symmetric injection strategies, all of which are designed to maintain the same global mean surface temperature: an annual injection at the Equator (EQ), an annual injection of equal amounts of SO 2 at 15° N and 15° S (15N+15S), an annual injection of equal amounts of SO 2 at 30° N and 30° S (30N+30S), and a polar injection strategy that injects equal amounts of SO 2 at 60° N and 60° S only during spring in each hemisphere (60N+60S). We compare these four hemispherically symmetric SAI strategies with a more complex injection strategy that injects different quantities of SO 2 at 30° N, 15° N, 15° S, and 30° S in order to maintain not only the global mean surface temperature but also its large-scale horizontal gradients. All five strategies are simulated using version 2 of the Community Earth System Model with the middle atmosphere version of the Whole Atmosphere Community Climate model, version 6, as the atmospheric component, CESM2(WACCM6-MA), with the global warming scenario, Shared Socioeconomic Pathway (SSP)2-4.5. We find that the choice of SAI strategy affects the spatial distribution of aerosol optical depths, injection efficiency, and various surface climate responses. In addition, injecting in the subtropics produces more global cooling per unit injection, with the EQ and the 60N+60S cases requiring, respectively, 59 % and 50 % more injection than the 30N+30S case to meet the same global mean temperature target. Injecting at higher latitudes results in larger Equator-to-pole temperature gradients. While all five strategies restore Arctic September sea ice, the high-latitude injection strategy is more effective due to the SAI-induced cooling occurring preferentially at higher latitudes. These results suggest trade-offs wherein different strategies appear better or worse, depending on which metrics are deemed important.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Modeling of Microgrid for Critical Data Center Applications

As part of continuing efforts to develop support, understanding, and infrastructure of data center integration onto the electric grid, our work aims to model and analyze the behavior of data center loads in a microgrid power system. Our model consists of renewable energy sources, batteries, and a nuclear reactor-steam Rankine cycle to power data center loads. The simulation studies investigate the electrical behavior of the microgrid system to assess its ability in supporting large data center electrical demands.

14 - SOLAR ENERGY↗