Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled 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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Refined Understanding of the Ice Cloud Longwave Scattering Effects in Climate Model

Abstract Because longwave (LW) absorption by greenhouse gases and clouds is more significant than the LW scattering effect by clouds, most climate models neglect cloud LW scattering to save computational costs. Ignoring cloud LW scattering directly overestimates outgoing longwave radiation (OLR). This study included ice‐cloud LW scattering treatment in the Exascale Energy Earth System Model (E3SM) version 2 and ran fully‐coupled simulations, prescribed sea surface temperature simulations, and offline radiative transfer calculations to comprehensively assess the impact of ice‐cloud LW scattering on global climate simulation. The instantaneous effect due to ice‐cloud LW scattering reduces the OLR by ∼1 W/m 2 on the global average and 2 W/m 2 on the tropical average. Tropospheric warming and high cloud amount reduction act to partially compensate for such instantaneous OLR reduction caused by the inclusion of LW scattering. When the simulation reaches the equilibrium, the surface warms by 0.66 K on average with respect to the simulation without LW scattering, with the Arctic surface temperature differences more than twice as large as that of the global mean. The impact of including LW scattering on the simulated climate change in response to 4 × CO 2 is also assessed. While including the cloud LW scattering does not significantly modify radiative forcing and total radiative feedback under such a scenario, it results in a 10% more positive cloud feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Model‐Based Interpretation of Solute Exports and Carbon Partitioning During Shale Weathering in a Mountainous Hillslope

The weathering of sedimentary rocks in high-elevation catchments influences freshwater quality and the global carbon cycle. While individual biogeochemical mechanisms involved in this process are relatively well understood, quantifying their contributions to solute export and carbon fluxes under natural, transient conditions remains challenging. Here, we implement a numerical multidimensional and multiphase model to simulate coupled hydrological and biogeochemical processes in a shale-underlain, snow-dominated hillslope in the Rocky Mountains, Colorado. The model captures the dynamic interplay between soil respiration, mineral weathering, and climate-driven hydrological forcing, reproducing observed soil CO 2 dynamics, groundwater chemistry, and subsurface flow. Our results reveal that seasonal snowmelt enhances carbonate weathering by promoting the infiltration of CO 2 -rich water to depth, while pyrite oxidation is primarily sensitive to low water saturation that facilitates O 2 diffusion through the regolith. Topography modulates the spatial distribution of shale weathering, as steeper slopes enhance lateral drainage, favoring the delivery of reactants to greater depths. While shale weathering at our site acts as a transient carbon sink, with silicates and carbonates buffering acidity and promoting atmospheric CO 2 consumption (1% of soil-derived CO 2 ), the exported dissolved inorganic carbon is predominantly geogenic (∼73%). Consequently, when accounting for long-term marine carbonate precipitation. The current weathering regime represents a net source of carbon to the atmosphere. The oxidation of pyrite and petrogenic organic carbon together release approximately 0.9 mol·m −2 ·yr −1 of CO 2 . Our findings highlight the role of topography, hydroclimate, and the coupling between acid-base reactions in shaping the carbon balance and the solute exports in mountainous critical zones.

carbon cycling↗
DOE OSTI2024

Active species in chloroaluminate ionic liquids catalyzing low-temperature polyolefin deconstruction

Abstract Chloroaluminate ionic liquids selectively transform (waste) polyolefins into gasoline-range alkanes through tandem cracking-alkylation at temperatures below 100 °C. Further improvement of this process necessitates a deep understanding of the nature of the catalytically active species and the correlated performance in the catalyzing critical reactions for the tandem polyolefin deconstruction with isoalkanes at low temperatures. Here, we address this requirement by determining the nuclearity of the chloroaluminate ions and their interactions with reaction intermediates, combining in situ 27 Al magic-angle spinning nuclear magnetic resonance spectroscopy, in situ Raman spectroscopy, Al K-edge X-ray absorption near edge structure spectroscopy, and catalytic activity measurement. Cracking and alkylation are facilitated by carbenium ions initiated by AlCl 3 - tert -butyl chloride (TBC) adducts, which are formed by the dissociation of Al 2 Cl 7 − in the presence of TBC. The carbenium ions activate the alkane polymer strands and advance the alkylation cycle through multiple hydride transfer reactions. In situ 1 H NMR and operando infrared spectroscopy demonstrate that the cracking and alkylation processes occur synchronously; alkenes formed during cracking are rapidly incorporated into the carbenium ion-mediated alkylation cycle. The conclusions are further supported by ab initio molecular dynamics simulations coupled with an enhanced sampling method, and model experiments using n-hexadecane as a feed.

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

Mesoscopic analyses of the impact of morphology and operating conditions on the transport resistances in a proton-exchange-membrane fuel-cell catalyst layer

Exploring the origins of local transport resistance and characterizing the oxygen transport resistances in the catalyst layer (RCL) are critical for cost reduction. In this paper, a comprehensive mesoscopic model for simulating coupled transport processes of oxygen and water vapor for different structural parameters under different operating conditions in reconstructed microstructures is proposed. The local transport resistance is calculated after achieving the limiting current density and the effective diffusivity of oxygen. The results demonstrate that RCL increases greatly with decreasing platinum loading (LPt) and the transport resistances in other components of the cell dominate for high-loadings. Both the reduced oxygen permeation coefficient in the ionomer thin-film and the adsorption resistance account for the origins of local transport resistance. The local transport resistance increases with the bare carbon ratio for a constant LPt and Pt/C ratio due to the decreased effective ionomer surface, and increases with the I/C ratio due to the increased ionomer thickness and decreased Knudsen diffusivity. Due to the presence of liquid water, a slight decrease followed by an increase of the local transport resistance versus relative humidity is obtained. The contribution of ionomer thin-films to RCL is more sensitive to liquid saturation compared with that of pores. This journal is

Mu, YT↗
DOE OSTI2023

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

Here, the exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron–ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron–nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

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

Construction of generalized quasilinear diffusion coefficient using neural networks with physical restrictions

The quasilinear diffusion coefficient (D QL ) derived from our machine learning framework shows comparable trends with the ground truth D QL obtained from GENRAY-CQL3D simulations. Additionally, for the strong absorption cases, the radial current drive profiles generated using the D QL from our model exhibit consistent behavior with those obtained from the original simulation. These findings indicate the potential of our surrogate modeling approach with physical restrictions to replicate key wave–plasma interaction characteristics while reducing computational costs. Traditionally, calculating D QL for wave–particle interactions relies on computationally intensive wave simulations coupled with Fokker–Planck solvers. To address this challenge, we developed a machine learning-based surrogate model with physical restrictions derived from cold plasma theory and bounce-averaged damping effects. First, we establish the propagation domain of Lower Hybrid Waves in the (N∥, ρ) space by identifying the accessibility limit and determining the upper and lower bounds of N∥ using the Potential Power Deposition (PPD) method. Subsequently, leveraging a database constructed using Latin hypercube sampling alongside the underlying physical restrictions (e.g. PPD), machine learning methods including U-Net and Recurrent Neural Networks are employed to design a physics-restricted machine learning framework capable of reconstructing D QL .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Correlations between asymmetric compression, burn amplification, and hot-spot velocities in inertial confinement fusion implosions

Here, this manuscript examines the correlations between the hot-spot velocity (an observable signature of residual kinetic energy), low-mode implosion asymmetries, and burn amplification in inertial confinement fusion implosions on the National Ignition Facility (NIF). Using a combination of two-dimensional axis-symmetric and three-dimensional radiation-hydrodynamic simulations coupled to neutronics, we find that for typical NIF implosions, the stagnation asymmetry multiplies the observed hot-spot velocity anywhere from 80% to 120%, while burn amplification always increases it. Additionally, we find stagnation asymmetry typically deflects the observed hot-spot flow. The two mechanisms (low-mode implosion asymmetries and burn amplification) can be decoupled, and application of a simple model to a database of cryogenic implosions on the NIF infers the total hot-spot velocity amplification. This finding modifies the interpretation of data collected from inertial confinement fusion experiments and impacts the magnitude and origin of low-mode asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Core-edge integrated predictive studies of ST40 and NSTX plasmas with the scrape-off layer box model

The ability to model the interplay between the core and edge of tokamak plasmas is crucial to designing both the plasma operating scenario of a fusion pilot plant and the design of the tokamak itself. Scrape-off-layer (SOL) models that are tailored to integrated scenario modeling need to have fast turn-around time and minimal computational burden to enable wide parameter-space coverage for design scoping. The SOL 0-D Box model is a reduced SOL model based on global power and particle balance that captures the essential physics of SOL transport with little computational cost. The usage of the 0-D Box model in core-edge coupled simulations has been demonstrated in both interpretive and predictive modes on a variety of devices. This paper presents a sensitivity study of the 0-D Box model to the input SOL heat-flux width for an ST40 plasma. This study demonstrates that accurate prediction of this width is crucial to predicting global performance parameters of a plasma scenario, such as energy confinement time and flux consumption. We also present an extension of the Box model to 1-D to allow for parallel variation of plasma parameters along the magnetic field lines. The 1-D Box model is then compared with SOLPS-ITER simulations of an NSTX plasma. Advantages and limitations of the Box model are discussed, and future directions are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Conformational free energy landscape of a glutamate transporter and microscopic details of its transport mechanism

Removing glutamate from the synaptic cleft is vital for proper function of the brain. Excitatory amino acid transporters mediate this process by uptaking the neurotransmitter from the synaptic cleft back to the cell after its release. The archaeal homolog, Glt Ph , an aspartate transporter fromPyrococcus horikoshii, presents the best structurally characterized model for this family of transporters. In order to transport, Glt Ph undergoes elevator-like conformational changes between inward-facing (IF) and outward-facing (OF) states. Here, we characterize, at an atomic level, the OF⇌IF transition of Glt Ph in differentapo/bound states using a combination of ensemble-based enhanced sampling techniques, employing more than two thousand of coupled simulation replicas of membrane-embedded Glt Ph . The resulting free-energy profiles portray the transition ofapo/bound states as a complex four-stage process, while sodium binding alone locks the structure in one of its states. Along the transition, the transport domain (TD) disengages from the scaffold domain (SD), allowing it to move as a piston sliding vertically with respect to the membrane during the elevator-like motion of TD. Lipid interactions with residues comprising the SD–TD interface directly influence the large-scale conformational changes and, consequently, the energetics of transport. Structural intermediates formed during the transition leak water molecules and may correlate to the uncoupled Cl − ion conductance observed experimentally in both prokaryotic and mammalian glutamate transporters. Mechanistic insights obtained from our study provide a structural framework for better development of therapeutic for neurological disorders.

Science & Technology - Other Topics↗
DOE OSTI2026

Balance of Plant Modeling and Real-Time Hardware-in-the-Loop Integration with the Microreactor Automated Control System

The advent of novel microreactor technology has driven a focused effort to explore safety and efficiency improvements that can be achieved through the use of automated system control. Development of control strategies, especially for initial demonstration, requires an adequate surrogate environment to safely research failure modes and control integration with realistic hardware delay. However, efficiency gains from control strategies are improved when the scope of controller action is expanded to include system-level dynamics such as downstream heat extraction and mass flow. For this reason, a balance-of-plant (BOP) model of a representative microreactor system has been developed using the TRANsient Simulation Framework of Reconfigurable Models library in Modelica. This model captures a reactor and primary NaK coolant loop that represent corresponding system components of the Microreactor Applications Research Validation and EvaLuation (MARVEL) design as well as a secondary coolant loop and heat extraction representative of the Microreactor Agile Non-Nuclear Experimental Test Bed (MAGNET). This model configuration allows for hardware-in-the-loop (HIL) integration with microreactor automated control system (MACS) hardware in real time through a Python-based gRPC client. Real-time simulation of model performance with emulated hardware and communication delay suggests that under independent proportional-integral-derivative control of BOP model drum dynamics and downstream heat extraction, stable power load following is achievable. A slight delay in load following, filtering of high-frequency dynamics, and localized temperature fluctation suggest room for improvement through the development of higher-level control strategies. The simulated coupling of the MAGNET facility lays the groundwork for future digital twin analysis with a coupled MACS-MAGNET HIL demonstration.

McConnell, Jono [ORNL] (ORCID:0000000238984741)↗
DOE OSTI2023

Core-Pedestal Plasma Configurations in Advanced Tokamaks

Here, several configurations for the core and pedestal plasma are examined for a predefined tokamak design by implementing multiple heating/current drive (H/CD) sources to achieve an optimum configuration of high fusion power in a noninductive operation while maintaining an ideally magnetohydrodynamic (MHD) stable core plasma using the IPS-FASTRAN framework. IPS-FASTRAN is a component-based lightweight coupled simulation framework that is used to simulate magnetically confined plasma by integrating a set of high-fidelity codes to construct the plasma equilibrium (EFIT, TOQ, and CHEASE), calculate the turbulent heat and particle transport fluxes (TGLF), model various H/CD systems (TORIC, TORAY, GENRAY, and NUBEAM), model the pedestal pressure and width (EPED), and estimate the ideal MHD stability (DCON). The TGLF core transport model and EPED pedestal model are used to self-consistently predict plasma profiles consistent with ideal MHD stability and H/CD (and bootstrap) current sources. In order to evaluate the achievable and sustainable plasma beta, varying configurations are produced ranging from the no-wall stability to with-wall stability regimes, simultaneously subject to the self-consistent TGLF, EPED, and H/CD source profile predictions that optimize configuration performance. The pedestal density, plasma current, and total injected power are scanned to explore their impact on the target plasma configuration, fusion power, and confinement quality. A set of fully noninductive scenarios are achieved by employing ion-cyclotron, neutral beam injection, helicon, and lower-hybrid H/CDs to provide a broad profile for the total current drive in the core region for a predefined tokamak design. These noninductive scenarios are characterized by high fusion gain (Q ~ 4) and power (P fus ~ 600 MW), optimum confinement quality (H 98 ~ 1.1), and high bootstrap current fraction (f BS ~ 0.7) for Greenwald fraction below unity. The broad current profile configurations identified are stable to low-n kink modes either because the normalized pressure β N is below the no-wall limit or a wall is present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Prehydrodynamic evolution and its signatures in final-state heavy-ion observables

We investigate the effects of prehydrodynamic evolution on final-state observables in heavy-ion collisions using state-of-the art event simulations coupled to different prehydrodynamic scenarios, which include the recently developed effective kinetic transport theory evolution model KøMPøST. Differential flow observables are found to be mostly insensitive to the details of prehydrodynamic evolution. The main effect we observe is in the p T spectra, particularly the mean transverse momentum. However, at least part of this effect is a consequence of the underlying conformal invariance assumption currently present in such approaches, which is known to be violated in the temperature regime probed in heavy-ion collisions. This assumption of early time conformal invariance leads to an artificially large out-of-equilibrium bulk pressure when switching from (conformal) prehydrodynamic evolution to hydrodynamics (using the nonconformal QCD equation of state), which in turn increases the transverse momentum. Our study indicates that a consistent treatment of prehydrodynamic evolution in heavy-ion collisions requires the use of nonconformal models of early time dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Micropinch formation dynamics in X pinches

High temporal resolution x-ray streak camera studies of micropinch formation in Cu hybrid x pinches reveal key plasma conditions. Analysis of Ne-like Cu lines indicate an average electron temperature of about 200 eV and 4.5×10 28 m -3 electron density. Here, the spectra suggest that the electron temperature jumps to about 1 keV, inferred from the continuum and the postcontinuum line emission that includes Li-like Cu lines. There is no sign of a rapid temperature change or a substantial surge in radiation emission during the 200 ps precontinuum x-ray burst, suggesting that the radiative collapse process does not play a major role in micropinch formation. Two-dimensional extended Magnetohydrodynamic (MHD) simulations, coupled to a collisional-radiative spectral analysis code, suggest the significance of the rapid radial implosion of high-temperature, low-density plasma, the axial outflow, and the dynamic plasma pressure in micropinch formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Picosecond spin-orbit torque–induced coherent magnetization switching in a ferromagnet

Electrically controllable nonvolatile magnetic memories show great potential for the replacement of conventional semiconductor-based memory technologies. Here, we experimentally demonstrate ultrafast spin-orbit torque (SOT)-induced coherent magnetization switching dynamics in a ferromagnet. We use an ultrafast photoconducting switch and a coplanar strip line to generate and guide a ~9-picosecond electrical pulse into a heavy metal/ferromagnet multilayer to induce ultrafast SOT. We then use magneto-optical probing to investigate the magnetization dynamics with sub-picosecond resolution. Ultrafast heating by the approximately 9 picosecond current pulse induces a thermal anisotropy torque which, in combination with the damping-like torque, coherently rotates the magnetization to obtain zero-crossing of magnetization in ~70 picoseconds. A macro-magnetic simulation coupled with an ultrafast heating model agrees well with the experiment and suggests coherent magnetization switching without any incubation delay on an unprecedented time scale. Our work proposes a unique magnetization switching mechanism toward markedly increasing the writing speed of SOT magnetic random-access memory devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

PVcircuit [SWR-22-26]

The software contains objects that are building blocks for PV modeling and interactive data fitting based on: Optoelectronic models for tandem/multijunction solar cells including resistive and luminescent coupling; simulation of modules composed of 2T, 3T, and 4T tandem solar cells; and energy yield analysis of PV systems composed of tandem solar cells.

Geisz, John↗
DOE OSTI2022

PVcircuit v0.0.6 [SWR-22-26]

The software contains objects that are building blocks for PV modeling and interactive data fitting based on: Optoelectronic models for tandem/multijunction solar cells including resistive and luminescent coupling; simulation of modules composed of 2T, 3T, and 4T tandem solar cells; and energy yield analysis of PV systems composed of tandem solar cells.

Geisz, John↗
DOE OSTI2026

TEC (ThermoElastoCreep) [SWR-25-167]

This repository provides a Python script to simulate the thermomechanical behavior of Concentrating Solar Power (CSP) receivers under periodic thermal loading. The simulation couples elastic and creep strains to predict stress relaxation, fatigue, and rupture life of the receiver structure. The work was done under the project Light Trapping Planar Cavity Receiver (LTPCR) funded by U.S. Department of Energy. All other relevant work and publications from LTPCR project can be found at https://github.com/NatLabRockies/LTPCR-Public

He, Xin [National Laboratory of the Rockies (NLR),↗