Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation System”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

From Cell to System: Accelerated hpc Simulations of BESS Aging under Frequency Regulation and Arbitrage use cases

Lithium-ion battery energy storage systems (BESS) packs have emerged as a leading solution for grid-scale energy storage, enhancing resiliency and balancing load fluctuations. Yet, experimental characterization of large-format LIB packs-particularly to assess performance and degradation over hundreds of cycles - demands substantial hardware investment and multi-year testing campaigns. In this work, we couple a hierarchical, physics-based modeling framework agnostic to electrode chemistries with high-performance computing to accelerate systems level evaluation by upto two orders of magnitude. Building on the open-source liionpack platform, we implement cell, module, and pack-scale electrochemical models enriched with mechanistic aging mechanisms and deploy them on an HPC cluster to simulate 150−200kWh systems over 500 - 1,000 cycles with in days. We subject these virtual B ESS to both constant-current cycling and realistic grid service profiles spanning frequency regulation, ramp-rate support, and energy arbitrage-and quantify the resulting degradation patterns. Our results reveal that localized cell aging can induce substantial nonuniformity at module and pack levels, with service-specific cycling protocols driving distinct aging modes. This rapid, multiscale modeling approach provides a powerful design-space exploration tool for optimizing electrical architecture, control strategies, and operational schedules to prolong pack lifetime and lower total cost of ownership.

Ayalasomayajula, Surya [ORNL] (ORCID:0009000860788↗
DOE OSTI2023

Constraints on simulated past Arctic amplification and lapse rate feedback from observations

The Arctic has warmed more rapidly than the global mean during the past few decades. The lapse rate feedback (LRF) has been identified as being a large contributor to the Arctic amplification (AA) of climate change. This particular feedback arises from the vertically non-uniform warming of the troposphere, which in the Arctic emerges as strong near-surface and muted free-tropospheric warming. Stable stratification and meridional energy transport are two characteristic processes that are evoked as causes for this vertical warming structure. Our aim is to constrain these governing processes by making use of detailed observations in combination with the large climate model ensemble of the sixth Coupled Model Intercomparison Project (CMIP6). We build on the result that CMIP6 models show a large spread in AA and Arctic LRF, which are positively correlated for the historical period of 1951–2014. Thus, we present process-oriented constraints by linking characteristics of the current climate to historical climate simulations. In particular, we compare a large consortium of present-day observations to co-located model data from subsets that show a weak and strong simulated AA and Arctic LRF in the past. Our analyses suggest that the vertical temperature structure of the Arctic boundary layer is more realistically depicted in climate models with weak (w) AA and Arctic LRF (CMIP6/w) in the past. In particular, CMIP6/w models show stronger inversions in the present climate for boreal autumn and winter and over sea ice, which is more consistent with the observations. These results are based on observations from the year-long Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic, long-term measurements at the Utqiagvik site in Alaska, USA, and dropsonde temperature profiling from aircraft campaigns in the Fram Strait. In addition, the atmospheric energy transport from lower latitudes that can further mediate the warming structure in the free troposphere is more realistically represented by CMIP6/w models. In particular, CMIP6/w models systemically simulate a weaker Arctic atmospheric energy transport convergence in the present climate for boreal autumn and winter, which is more consistent with fifth generation reanalysis of the European Centre for Medium-Range Weather Forecasts (ERA5). We further show a positive relationship between the magnitude of the present-day transport convergence and the strength of past AA. With respect to the Arctic LRF, we find links between the changes in transport pathways that drive vertical warming structures and local differences in the LRF. This highlights the mediating influence of advection on the Arctic LRF and motivates deeper studies to explicitly link spatial patterns of Arctic feedbacks to changes in the large-scale circulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulation of open quantum systems via low-depth convex unitary evolutions

Simulating physical systems on quantum devices is one of the most promising applications of quantum technology. Current quantum approaches to simulating open quantum systems are still practically challenging on NISQ-era devices, because they typically require ancilla qubits and extensive controlled sequences. In this work, we propose a hybrid quantum-classical approach for simulating a class of open system dynamics called random-unitary channels. These channels naturally decompose into a series of convex unitary evolutions, which can then be efficiently sampled and run as independent circuits. The method does not require deep ancilla frameworks and thus can be implemented with lower noise costs. We implement simulations of open quantum systems up to dozens of qubits and with large channel ranks. Published by the American Physical Society 2024

Peetz, Joseph (ORCID:0000000274458028)↗
DOE OSTI2014

Molecular simulation study of role of polymer–particle interactions in the strain-dependent viscoelasticity of elastomers (Payne effect)

The strain-amplitude dependence of viscoelastic behavior of model crosslinked elastomers containing various concentrations of spherical nanoparticles (NPs) was studied by non-equilibrium molecular dynamics simulation. All the filler NPs were in monodispersed state and the interactions between these particles were purely repulsive. The polymer–particle interactions were attractive and their interaction energies were tuned in a broad range. Through the computational study, many important features of the behavior of particle-reinforced elastomers observed in experiments, including the Payne effect, were successfully reproduced. It was shown that the magnitude of the Payne effect was found to depend on the polymer–particle interaction and the filler loading. By examining the microstructures of the simulation systems and their evolution during oscillatory shear, four different mechanisms for the role of the polymer–particle interactions in the Payne effect were revealed that consist of the debonding of polymer chains from NP surfaces, the breakage of polymer-shell-bridged NP network, the rearrangement of the NPs in the network into different layers and the shear-induced yielding of the rigid polymer shell in-between neighboring NPs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

On the feasibility of monitoring carbon monoxide in the lower troposphere from a constellation of northern hemisphere geostationary satellites: Global scale assimilation experiments (Part II)

In this paper, we describe the second phase of an Observing System Simulation Experiment (OSSE) that utilizes the synthetic measurements from a constellation of satellites measuring atmospheric composition from geostationary (GEO) Earth orbit presented in part I of the study. Our OSSE is focused on carbon monoxide observations over North America, East Asia and Europe where most of the anthropogenic sources are located. Here we assess the impact of a potential GEO constellation on constraining northern hemisphere (NH) carbon monoxide (CO) using data assimilation. We show how cloud cover affects the GEO constellation data density with the largest cloud cover (i.e., lowest data density) occurring during Asian summer. We compare the modeled state of the atmosphere (Control Run), before CO data assimilation, with the known “true” state of the atmosphere (Nature Run) and show that our setup provides realistic atmospheric CO fields and emission budgets. Overall, the Control Run underestimates CO concentrations in the northern hemisphere, especially in areas close to CO sources. Assimilation experiments show that constraining CO close to the main anthropogenic sources significantly reduces errors in NH CO compared to the Control Run. We assess the changes in error reduction when only single satellite instruments are available as compared to the full constellation. We find large differences in how measurements for each continental scale observation system affect the hemispherical improvement in long-range transport patterns, especially due to seasonal cloud cover. Finally, a GEO constellation will provide the most efficient constraint on NH CO during winter when CO lifetime is longer and increments from data assimilation associated with source regions are advected further around the globe.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Constraints on methane emissions in North America from future geostationary remote-sensing measurements

The success of future geostationary (GEO) satellite observation missions depends on our ability to design instruments that address their key scientific objectives. Here, an Observation System Simulation Experiment (OSSE) is performed to quantify the constraints on methane (CH 4 ) emissions in North America obtained from shortwave infrared (SWIR), thermal infrared (TIR), and multi-spectral (SWIR+TIR) measurements in geostationary orbit and from future SWIR low-Earth orbit (LEO) measurements. Furthermore, we used an efficient stochastic algorithm to compute the information content of the inverted emissions at high spatial resolution (0.5° × 0.7°) in a variational framework using the GEOS-Chem chemistry-transport model and its adjoint. Our results show that at sub-weekly timescales, SWIR measurements in GEO orbit can constrain about twice as many independent flux patterns than in LEO orbit, with a degree of freedom for signal (DOF) for the inversion of 266 and 115, respectively. Comparisons between TIR GEO and SWIR LEO configurations reveal that poor boundary layer sensitivities for the TIR measurements cannot be compensated for by the high spatiotemporal sampling of a GEO orbit. The benefit of a multi-spectral instrument compared to current SWIR products in a GEO context is shown for sub-weekly timescale constraints, with an increase in the DOF of about 50 % for a 3-day inversion. Our results further suggest that both the SWIR and multi-spectral measurements on GEO orbits could almost fully resolve CH 4 fluxes at a spatial resolution of at least 100 km × 100 km over source hotspots (emissions > 4 × 10 5 kg day -1 ). The sensitivity of the optimized emission scaling factors to typical errors in boundary and initial conditions can reach 30 and 50 % for the SWIR GEO or SWIR LEO configurations, respectively, while it is smaller than 5 % in the case of a multi-spectral GEO system. Our results demonstrate that multi-spectral measurements from a geostationary satellite platform would address the need for higher spatiotemporal constraints on CH 4 emissions while greatly mitigating the impact of inherent uncertainties in source inversion methods on the inferred fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Study of the Effect of Mold Corner Shape on the Initial Solidification Behavior of Molten Steel Using Mold Simulator

The chamfered mold with a typical corner shape (angle between the chamfered face and hot face is 45 deg) was applied to the mold simulator study in this paper, and the results were compared with the previous results from a well-developed right-angle mold simulator system. The results suggested that the designed chamfered structure would increase the thermal resistance and weaken the two-dimensional heat transfer around the mold corner, causing the homogeneity of the mold surface temperatures and heat fluxes. In addition, the chamfered structure can decrease the fluctuation of the steel level and the liquid slag flow around the meniscus at mold corner. The cooling intensities at different longitudinal sections of shell are close to each other due to the similar time-average solidification factors, which are 2.392 mm/s{sup 1/2} (section A-A: chamfered center), 2.372 mm/s{sup 1/2} (section B-B: 135 deg corner), and 2.380 mm/s{sup 1/2} (section D-D: face), respectively. For the same oscillation mark (OM), the heights of OM roots at different positions (profile L1 (face), profile L2 (135 deg corner), and profile L3 (chamfered center)) are very close to each other. The average value of height difference (H{sub D}) between two OMs roots for L1 and L2 is 0.22 mm, and for L2 and L3 is 0.38 mm. Finally, with the help of metallographic examination, the shapes of different hooks were also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Advanced Instrumentation, Information, and Control Systems Technologies Pathway: FY 2016 External Review

This report describes an External Review conducted by the LWRS Program Advanced Instrumentation, Information, and Control (II&C) Systems Technologies Pathway to solicit feedback on the topics and results of the ongoing II&C research program. This review was held in conjunction with the Nuclear Energy Institute (NEI) Digital I&C Working Group meeting that was held at Idaho National Laboratory (INL) on August 9-10, 2016. Given the opportunity to visit INL and see the pathway research projects, NEI agreed that the Working Group would serve as the External Review panel for the purpose of obtaining expert input on the value and timing of the research projects. This consisted of demonstrations in the Human Systems Simulation Laboratory followed by presentations on the II&C research program in general as well as the five technology development areas. Following the meeting, the presentations were sent to each of the attendees so they could review them in more detail and refer to them in completing the feedback form. Follow-up activities were conducted with the attendees following the meeting to obtain the completed feedback forms. A total of 13 forms were returned. The feedback forms were reviewed by the pathway to compile the data and comments received, which are documented in the report. In all, the feedback provided by the External Review participants is taken to be a strong endorsement of the types of projects being conducted by the pathway, the value they hold for the nuclear plants, and the general timing of need. The feedback aligns well with the priorities, levels of efforts allocated for the research projects, and project schedules. The feedback also represents realistic observations on the practicality of some aspects of implementing these technologies. In some cases, the participants provided thoughtful challenges to certain assumptions in the formulation of the technologies or in deployment plans. These deserve further review and revision of plans if warranted. The pathway will take all of the feedback and address the open issues that have been identified by the participants. This includes 11 actionable items for follow up by the II&C Pathway.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Efficient Simulation of Open Quantum Systems on NISQ Trapped‐Ion Hardware

Abstract Simulating open quantum systems, which interact with external environments, presents significant challenges on noisy intermediate‐scale quantum (NISQ) devices due to limited qubit resources and noise. In this study, an efficient framework is proposed for simulating open quantum systems on NISQ hardware by leveraging a time‐perturbative Kraus operator representation of the system's dynamics. This approach avoids the computationally expensive Trotterization method and exploits the Lindblad master equation to represent time evolution in a compact form, particularly for systems satisfying specific commutation relations. The efficiency of this method is demonstrated by simulating quantum channels, such as the continuous‐time Pauli channel and damped harmonic oscillators, on NISQ trapped‐ion hardware, including IonQ Harmony and Quantinuum H1‐1. Additionally, hardware‐agnostic error mitigation techniques are introduced, including Pauli channel fitting and quantum depolarizing channel inversion, to enhance the fidelity of quantum simulations. These results show strong agreement between the simulations on real quantum hardware and exact solutions, highlighting the potential of Kraus‐based methods for scalable and accurate simulation of open quantum systems on NISQ devices. This framework opens pathways for simulating more complex systems under realistic conditions in the near term.

Burdine, Colin [Department of Electrical and Compu↗
DOE OSTI2025

Phase change material integrated core–shell catalyst for in situ thermal control in methanol synthesis from syngas

A model-guided core–shell catalyst design is presented for methanol synthesis, featuring a phase change material (PCM) core encapsulated by a Cu–Zn–Al 2 O 3 (CZA) catalytic shell. The PCM enables in situ thermal management by absorbing reaction heat at its melting point, mitigates the kinetic decline at high temperatures and therefore avoids low conversion, prevents hot spots, and stabilizes the reaction temperature. A two-dimensional axisymmetric, non-isothermal packed-bed reactor model (COMSOL 6.3) was developed for a 10 g system. Simulations evaluate three PCM candidates, that is, LiNO 3 , 9 wt% LiCl + 91 wt% LiNO 3 , and commercial H250, with melting points near 244–250°C. Results indicate that CO conversion can increase from 34.4% to 52.4%, and methanol production can improve by 69% compared to a conventional packed-bed reactor. Beyond methanol synthesis, the PCM-integrated core–shell concept provides a scalable approach for thermal control in exothermic reactions, improving reactor efficiency and safety.

core–shell catalyst↗
DOE OSTI2019

Multi‐model analysis of climate impacts on plant photosynthesis in China during 2000–2015

Abstract Differences, arising from differences in gross primary production (GPP) model structures and driving forces, have fuelled arguments concerning interannual changes of GPP in China since 2000. To better investigate the interannual variability of GPP and its covariance with climate factors in China, this study adopted a multi‐model analysis based on three GPP models (i.e., Terrestrial Ecosystem Carbon flux model [TEC], Breathing Earth System Simulator model [BESS], and MOD17 GPP model). The results show that annual GPP in China increased by 0.021–0.057 Pg C year −1 from 2000 to 2015 attributable to atmospheric‐CO 2 fertilization effects and favourable climate change, that is, increasing precipitation ( P r ) and temperature ( T a ). However, northern China and southern China had a large difference in the amplitude of these GPP changes; annual GPP increased by 0.017–0.039 Pg C year −1 in northern China but only 0.001–0.018 Pg C year −1 in southern China. Northern China and southern China occupy contrasting climate zones and this contrast produced different interannual variability of GPP through different mechanisms. Northern China has a dry climate with GPP changes sensitive to P r . As a result, more P r along with higher T a in northern China produced the strong uptrend of GPP from 2000 to 2015. In contrast, southern China has a wet climate with its GPP sensitive to solar radiation and T a . For the interval of 2000–2015, decreasing radiation plus drought exerted a negative influence on GPP in southern China. This study highlights the diverse mechanisms in which climate change affects GPP in dry and wet climate zones. A robust multi‐model analysis is preferred to reduce uncertainties arising from a single GPP model and its driving data.

Yan, Hao↗
DOE OSTI2020

Permeation of CO 2 and N 2 through glassy poly(dimethyl phenylene) oxide under steady- and presteady-state conditions

Glassy polymers are often used for gas separations because of their high selectivity. Although the dual-mode permeation model correctly fits their sorption and permeation isotherms, its physical interpretation is disputed, and it does not describe permeation far from steady state, a condition expected when separations involve intermittent renewable energy sources. To develop a more comprehensive permeation model, we combine experiment, molecular dynamics, and multiscale reaction–diffusion modeling to characterize the time-dependent permeation of N 2 and CO 2 through a glassy poly(dimethyl phenylene oxide) membrane, a model system. Simulations of experimental time-dependent permeation data for both gases in the presteady-state and steady-state regimes show that both single- and dual-mode reaction–diffusion models reproduce the experimental observations, and that sorbed gas concentrations lag the external pressure rise. The results point to environment-sensitive diffusion coefficients as a vital characteristic of transport in glassy polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Generation of Reversible C++ Code for Optimistic Parallel Discrete Event Simulation

The reversible execution of C/C++ code has been a target of research and engineering for more than a decade as reversible computation has become a central notion in large-scale parallel discrete event simulation (PDES). The simulation models that are implemented for PDES are of increasing complexity and size and require various language features to support abstraction, encapsulation, and composition when building a simulation model. In this paper, we focus on parallel simulation models that are written with user-defined C++ abstractions and abstractions of the C++ Standard Library. We present an approach based on incremental state saving for establishing reversibility of C++ and an evaluation for a kinetic Monte-Carlo simulation implemented in C++. Although a significant runtime overhead is introduced with our technique, it is an enormous win that it allows using the entire C++ language, and has that code automatically transformed into reversible code to enable parallel execution with the Rensselaer’s optimistic simulation system (ROSS).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Tribo-corrosion and Albumin Attachment of Nitrogen Ion-Implanted CoCrMo Alloy During Friction Onset

In this paper, CoCrMo alloy surface was implanted with 100 keV nitrogen ions to modify it. Bovine serum albumin (BSA) adsorption and the initial behavior of tribo-corrosion in the simulated system were studied. Nitrogen ion implantation can promote BSA dynamic adsorption due to the change of friction and wear mechanisms. From the tribo-corrosion test results, the open circuit potential (OCP) increased to about 0.6 V and the coefficient of friction (COF) decreased to about 0.2 for the nitrogen ion-implanted CoCrMo compared with the untreated sample before the modified layer failure. The point when the open circuit potential and the coefficient of friction changed during long wear time (1, 2 and 4 h) is considered the sign of the worn through modified layer. Then, the OCP of the implanted sample rose by 0.3 V compared with the untreated sample, and the COF remained at around 0.3, which is lower compared with the COF of untreated sample after the modified layer has been worn through. Thus, nitrogen ion implantation not only improved wear and corrosion resistance of the CoCrMo alloy, but also promoted BSA adsorption on the CoCrMo alloy surface, which effectively reduced the wear volume and metal ions release.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Design of an experimental facility with a unit cell test section for studies of the lower plenum in prismatic high temperature gas reactors

A proposed design for a Generation IV reactor is the very-high temperature reactor (VHTR), which employs helium as the primary coolant. Although many advantages exist with this design, a number of challenges remain. One such area of concern is the turbulent mixing of non-isothermal flows in the lower plenum, where an array of coolant jets collectively mix and traverse around a series of structural support posts before exiting the core. Neighboring jets in the lower plenum could experience temperature differences as high as 300–400 K, causing potential hot streaking in the support posts. Additional problems could occur further downstream if the level of non-uniformity in the temperatures is too high (e.g, inducing thermal stresses on gas turbine blades). Therefore, it is important that adequate tools exist to predict the turbulent mixing within the VHTR lower plenum. This task remains difficult to achieve due to the complex geometry and flow physics that dictate the mixing. Various modeling approaches can only be assessed using data acquired from a carefully designed experimental facility. This paper presents details of such a facility, where high-fidelity measurements can be acquired for a scaled portion of the lower plenum defined herein as a ‘‘unit cell”, which consists of a hexagonal array of jets whose flow is directed into a test section with support posts positioned between neighboring jets, and subjected to a crossflow. Through adjusting velocity and temperature conditions at the jet inlets, different locations within the lower plenum can be experimentally simulated. System response quantities (SRQs) include both temperature and velocity measurements, captured via multiple thermocouples embedded in each post, and particle image velocimetry (PIV), respectively. These SRQs provide assessment metrics for both the momentum and energy models present in current computational solvers, and therefore the experimental facility is a valuable source for validation data. Another benefit is a more comprehensive understanding of the thermal-hydraulic loading conditions and flow physics present in the VHTR lower plenum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Characteristics of locational uncertainty marginal price for correlated uncertainties of variable renewable generation and demands

With the rapid increase of variable renewable energy sources in power systems, how to manage and price the uncertainty of renewable resources’ power outputs is becoming an urgent issue. Current market designs considering the uncertainties are mainly based on the probabilistic scenario set of demand and renewable energy resources power outputs. This consideration makes market designs vulnerable to three significant challenges when put into practice. First, the accurate probability distribution of renewable generation is hard to obtain in real-time. Second, it is challenging to clear the market timely with many scenarios to guarantee accuracy. Third, generation cost recovery cannot be guaranteed for some scenarios. To overcome these challenges, this paper proposes a locational uncertainty marginal price model to price the uncertainty explicitly based on a scenario-free stochastic market-clearing model. Instead of using the probabilistic scenario set, the uncertainty of renewable energy sources and loads is modeled with distributionally-robust chance constraints. The correlation of uncertainties can be endogenously modeled in both the market-clearing and the locational uncertainty marginal price formation. Furthermore, this paper proves that generation cost recovery, revenue adequacy, and partial market equilibrium can be achieved using the locational uncertainty marginal price model. Numerical results from both the small and large systems simulations validate that the generation cost recovery is maintained no matter the generation participates in uncertainty mitigation or not. The transmission congestion surplus is also allocated appropriately among loads, renewable energy sources, and financial transmission right owners.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2018

Collaborative operational fault tolerant control for stochastic distribution control system

Based on a class of industry processes, a distributed fault diagnosis and a collaborative operational fault tolerant control algorithms are presented for an irreversible interconnected stochastic distribution control system with boundary conditions. Different from the existed collaborative fault tolerant controllers, the controller does not need to change after fault occurs, the system tracking performance can be maintained by adjusting the setpoint of the whole system. Simulations are also included to show the effectiveness of the proposed algorithm and encouraging results have been obtained.

Stochastic distribution control, Fault diagnosis, ↗