Search NASA⌕ Search

SEARCH · Search NASA

Results for “co-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.

530 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Machine Learning for Predicting Multipactor Susceptibility in Planar RF Structures

Multipactor discharge is a persistent challenge in high-power microwave (HPM) and accelerator systems, where secondary electron avalanches can cause heating, vacuum degradation, and failure. This work presents the first supervised machine learning (ML) framework for multipactor prediction, trained on high-fidelity 3D Particle-in-Cell (PIC) simulation data in planar geometries. The model maps operational, geometric, and material-dependent secondary electron yield (SEY) parameters to the time-averaged electron growth rate, enabling rapid reconstruction of susceptibility charts. Among the models evaluated, tree-based ensemble methods such as Random Forest and Extra Trees demonstrate superior generalization to unseen materials compared to neural networks such as multilayer perceptron (MLP). Performance metrics, including Intersection over Union (IoU), Structural Similarity Index Measure (SSIM), and Pearson correlation, show close agreement with simulation benchmarks. Principal Component Analysis attributes generalization limits to material feature-space disjointedness.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Multiscale and multidimensional modeling of particle acceleration and transport in solar flares

Multi-messenger, multi-viewpoint, and time-resolved observations of solar flares are now providing unprecedented constraints on particle acceleration sites, energy conversion, and energy transport. The interpretation of current observations, including microwave imaging spectroscopy from EOVSA, hard x-ray (HXR) imaging from Solar Orbiter/STIX, gamma-ray diagnostics from Fermi, and in situ measurements from Parker Solar Probe and Solar Orbiter, collectively demands modeling frameworks that go beyond traditional spatially unresolved, one-zone models or single-mechanism descriptions. This review surveys multiscale and multidimensional modeling approaches, including kinetic, magnetohydrodynamic (MHD), and macroscopic particle models, that are being developed to meet the need. Kinetic simulations reveal that three-dimensional (3D) effects, including field-line chaos and self-generated turbulence, are essential for sustained power-law particle acceleration. MHD simulations now capture flux-rope eruptions, plasmoid-unstable current sheets, and turbulent flare regions in realistic magnetic topologies. Macroscopic models coupling MHD with energetic-particle models produce spatially resolved electron distributions and synthetic HXR and microwave emissions for direct comparison with observations. Despite these advances, outstanding challenges remain in bridging kinetic and global scales, improving MHD simulations and macroscopic particle models, and achieving quantitative model-observation closure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Machine learning the electric field response of condensed phase systems using perturbed neural network potentials

Abstract The interaction of condensed phase systems with external electric fields is of major importance in a myriad of processes in nature and technology, ranging from the field-directed motion of cells (galvanotaxis), to geochemistry and the formation of ice phases on planets, to field-directed chemical catalysis and energy storage and conversion systems including supercapacitors, batteries and solar cells. Molecular simulation in the presence of electric fields would give important atomistic insight into these processes but applications of the most accurate methods such as ab-initio molecular dynamics (AIMD) are limited in scope by their computational expense. Here we introduce Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) to push back the accessible time and length scales of such simulations. We demonstrate that important dielectric properties of liquid water including the field-induced relaxation dynamics, the dielectric constant and the field-dependent IR spectrum can be machine learned up to surprisingly high field strengths of about 0.2 V Å −1 without loss in accuracy when compared to ab-initio molecular dynamics. This is remarkable because, in contrast to most previous approaches, the two neural networks on which PNNP MD is based are exclusively trained on molecular configurations sampled from zero-field MD simulations, demonstrating that the networks not only interpolate but also reliably extrapolate the field response. PNNP MD is based on rigorous theory yet it is simple, general, modular, and systematically improvable allowing us to obtain atomistic insight into the interaction of a wide range of condensed phase systems with external electric fields.

Science & Technology - Other Topics↗
DOE OSTI2026

Evaluating Polymer Properties with Different Additives for Carbon Capture and Other Applications

Anthropogenic climate change is one of this generation’s most pressing concerns, with the potential to completely alter the delicate balance we’ve struck with nature. Already, global temperatures have risen 1.29°C, leading to disrupted weather systems, extinctions, increased risks of wildfires, and sea level rise, to name a few effects. Carbon dioxide emission from the combustion of fossil fuels and other industrial activity is a large driver of this phenomenon, as it absorbs heat before it can be radiated away from Earth, trapping it. Carbon dioxide has reached unprecedented levels in our atmosphere, showing a 50% increase from preindustrial averages to a whopping 430 ppm. Thus, reducing the amount of carbon dioxide via carbon capture technology is an important endeavor that serves to benefit everyone. The Microencapsulated CO 2 Sorbent (MECS) team at Lawrence Livermore National Laboratory (LLNL) has turned to microencapsulation to approach this endeavor. Microcapsules provide an attractive approach to carbon capture, combining large surface areas for more efficient mass transfer, regenerative abilities, reduced solvent loss, and improved handling. Additionally, while existing carbon capture technology relies on industrial plants, capsules could present a modular approach to carbon capture, reducing the need for extensive physical infrastructure. The MECS team’s design consists of a polymer membrane that contains a liquid carbon sequestering sorbent, aqueous sodium carbonate. The carbon capturing reaction occurs in three distinct steps, the first of which is the dissolution of carbon dioxide into the sorbent solution and its conversion into carbonic acid (H 2 CO 3 ), shown in equations 1 and 2 respectively. Because this step hinges upon the ability of carbon dioxide to reach the solution inside the capsule, it is necessary that the microcapsule shell is permeable to carbon dioxide gas. The MECS team produces these microcapsules using the in-air droplet encapsulation apparatus (IDEA) shown in figure 1, which can produce uniform micron-scale droplets at speeds much faster than traditional single-dispersal microfluidic-based techniques. The IDEA Is 100 times faster than these current techniques and can reach up to 1000 times their speed when incorporating a multi-nozzle design. Additionally, because droplets are produced in-air via vibration, IDEA can decrease post-processing times and material waste by 99% and can fabricate microgels that are 10 to 100 times more viscous than can be produced via traditional microfluidics. While this design represents a breakthrough in the throughput, efficiency, and tunability of microcapsule production, it imposes a major constraint on the microcapsule curing process. Because microcapsule shells are crosslinked with UV light while falling 30 cm through the air, this gives them a reaction window of approximately 0.2 seconds. Thus, the system and shell formulations must be optimized such that the shells can be fully crosslinked within this very narrow window, prompting investigations into curing behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effectively Considering the Distribution System in Integrated Resource Plans

While electricity planning practices vary by state and utility based on utility type and market structure, integrated resource planning (IRP) remains a prominent vehicle — even in states with centrally-organized wholesale electricity markets. IRP focuses on meeting forecasted long-term electricity needs. Typically, utilities have not considered impacts of design and operation of the low-voltage distribution network in IRP. With advanced capabilities of grid-edge technologies to generate and store electricity and provide load flexibility, and large utility investments in distribution systems, it's increasingly important to consider at least some distribution planning elements in IRP. This report considers the value proposition for doing so, such as reducing utility costs through resource co-optimization and strategic siting of grid-edge resources, and idenfities the most important touchpoints between planning for bulk power and distribution systems and provide a range of tactics for integrating these two processes.

Relf, Grace [Lawrence Berkeley National Laboratory↗
DOE OSTI2026

Precise Modeling of a Complex Solenoidal Magnetic Field Using a Combination of Analytic Functions and a PINN

We demonstrate an iterative approach to modeling a sparsely measured magnetic field in a large-bore solenoid. This approach uses a hybrid of traditional and machine learning techniques. The traditional technique is a linear least-squares fit using a series solution to Laplace's equation, while the machine learning technique involves the training of a physics-informed neural network (PINN) on the least-squares fit residuals. We use a newly defined activation function "DELTAsnake," a modification to the snake activation function proposed by Ziyin et al. that allows for stronger curvature and non-monotonicity. The combined model approximately obeys Maxwell's equations to a level sufficient for producing high quality physics simulations and analysis. Our approach is applied to a highly realistic calculation of the expected magnetic field in the Mu2e experiment's Detector Solenoid which includes a simple model for the expected statistical measurement uncertainties. Using ten toy measurement simulations, we demonstrate the capabilities of our model in comparison to the least-squares method alone; the least-squares method alone results in a reduced chi-squared statistic of ${2.15 \pm 0.01}$, while our approach improves the reduced chi-square to ${1.034 \pm 0.005}$. Furthermore, for an average toy simulation, we show that the range of the RMS of the three field component residuals reduces from ${0.07-0.37}$ Gauss to ${0.05-0.07}$ Gauss. We find that this novel method is robust against a realistic systematic uncertainty deriving from Hall probe calibration bias and can be used to significantly reduce the number of measurements required to achieve an accurate model.

Kampa, Cole [Caltech] (ORCID:0000000192972920)↗
DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗
DOE OSTI2026

Modeling of a Four-Stage Linear Ionization Cooling Channel for a Muon Collider in G4Beamline

A previous study of an eight-stage rectilinear ionization cooling channel in the ICOOL software demonstrated a five-order-of-magnitude reduction in a muon beam’s 6D emittance. In this study, we look to compare the ways ICOOL and Muons, Inc.’s g4Beamline software model ionization cooling by comparing their modeling of the first four stages to this optimized cooling channel constructed in ICOOL. We begin by identifying the parameters used to construct the optimized ionization cooling channel in ICOOL. We then reconstruct this beam in g4Beamline with identical parameter specifications and simulate the cooling of an identical input beam. Finally, we compare the two simulations based on their beam transmission, longitudinal emittance, and transverse emittance along the channel length. Through this process, we demonstrate that G4Beamline accurately reproduces transverse cooling results but predicts systematically different longitudinal emittance evolution while maintaining similar overall cooling performance, reproducing a 97.9% reduction in 6D emittance over four stages.

Keeler, Dominic [Purdue U., West Lafayette] (ORCID↗
DOE OSTI2026

Selective depolymerization of polyolefin copolymer via high-entropy alloy catalyzed pyrolysis in a double-layered ultrafast heating reactor

Chemical recycling of polyolefin copolymers (PCPs) is essential for mitigating plastic pollution, yet achieving selective depolymerization into C 2 –C 3 monomers at high productivity remains challenging. Here, we introduce a double-layered ultrafast heating reactor (DUHR) that integrates Joule heating with a high-entropy alloy (HEA) catalyst composed of near-equimolar Co, Cu, Fe, Mn, and Ni. We achieve highly selective depolymerization of PCP into ethylene and propylene, yielding a 70% non-condensable fraction. Significant aromatic fractions form via radical-driven pathways—a behavior distinct from conventional externally heated pyrolyzers—highlighting DUHR's ability to generate reactive species and modulate reaction pathways. Furthermore, these findings demonstrate an effective strategy for chemical recycling of real-world plastic waste for polyolefin circularity and sustainable resource recovery.

Chemical recycling↗
DOE OSTI2026

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗
DOE OSTI2026

Assembly bias and local Primordial non-Gaussianity from DESI DR1 quasars

The analysis of the large-scale clustering of quasars (QSO) observed by the Dark Energy Spectroscopic Instrument (DESI) represents a promising avenue for constraining local Primordial non-Gaussianity (PNG), parameterized by f NL . The signal to be constrained is the scale-dependent bias induced in the 2-point clustering of the considered tracer sample. The resulting constraints on f NL , however, are fully degenerate with the local PNG bias parameter b ϕ , dependent on the assembly bias parameter p. Using IllustrisTNG hydrodynamical simulations, we select a QSO sample reflecting the selection criteria and properties of DESI QSOs, and provide a robust prior for p, and thus for b ϕ , building on the findings of Fondi et al. 2024. We find a distribution with mean p̅ ≃ 1.4 with weak redshift dependence, stable to selection noise and consistent with the expected recent merger history typical of quasar-hosting halos. By comparing with the CAMELS simulations we demonstrate that this prior is robust to astrophysical assumptions and cosmic variance. Finally, applying this prior to the DESI DR1 dataset, we derive updated constraints on local PNG, obtaining f NL = -3.3±9.2.

cosmological parameters from LSS↗
DOE OSTI2026

Anomalous lattice expansion of bcc Ta in helium-loaded diamond anvil cells

Anomalous lattice expansion of bcc Ta is observed in helium-loaded diamond anvil cell (DAC) experiments near 6 GPa and above 400 K, while co-loaded W, Mo, and Au exhibit normal thermal expansion. The Ta volume expansion increases with temperature, reaching ~14% near 650 K before plateauing, and is largely retained after recovery to ambient conditions. Time-resolved measurements at 500 K reveal progressive expansion over several hours, consistent with gradual, thermally activated He incorporation. The largely retained expansion during subsequent cooling and holding at 400 K indicates that the incorporated He is predominantly kinetically trapped. The data reveal two distinct regimes: below ~550 K, modest and spatially uniform expansion consistent with interstitial He trapping at isolated, randomly distributed vacancies; above ~550 K, He-vacancy cluster growth producing dramatically larger expansion, peak broadening, and a core–shell microstructure revealed by micron-resolution X-ray diffraction mapping. These results demonstrate that He is not always a passive pressure medium in DAC experiments and reveal a previously unrecognized regime of He-mediated lattice modification in Ta, complementing conventional ion-implantation studies relevant to nuclear and fusion structural materials.

Diamond anvil cell↗
DOE OSTI2024

Crack‐Free Single‐Crystalline LiNiO 2 for High Energy Density All‐Solid‐State Batteries

Abstract Single‐crystalline layered oxide (LiNi 1‐ x ‐ y Mn x Co y O 2 ) cathodes have been found to exhibit exceptional electrochemical properties when coupled with various inorganic solid electrolytes (ISEs) in all‐solid‐state batteries (ASSBs). Their advantages stem from the robust morphological integrity with the absence of grain boundaries and the high electrochemical oxidative stability. Here, ASSBs featuring single‐crystalline LiNiO 2 (LNO) with the highest Ni content are reported, offering a high theoretical specific capacity of 275 mAh g ‐1 alongside a high average discharge voltage (3.7 V vs Li + /Li). Through a careful investigation, it is demonstrated that micron‐sized single‐crystalline LNO (µSC‐LNO) composite cathodes with a halide ISE exhibit a high initial discharge capacity of 205 mAh g ‐1 with an outstanding cycle performance over 200 cycles in room‐temperature ASSBs. The significance of engineering parameters is emphasized, such as particle size and specific density, in promoting a homogeneous and fast Li + transport within the composite cathodes. Furthermore, the formation of undesirable interphase between the halide ISE in the cathode and sulfide ISE separator is elucidated, which may be a critical factor impeding long‐term cyclability of ASSBs. This work provides insights into the design of composite cathodes for high‐energy‐density ASSBs.

Chemistry↗
DOE OSTI2026

A Probabilistic Approach to Load Modeling for Central HVAC Systems in Large Commercial Buildings for Retrofit Decisions Under Uncertainty

Retrofitting central HVAC systems in large commercial buildings with advanced technologies like heat recovery chillers (HRCs) offers a significant opportunity to enhance energy efficiency. However, analyzing these retrofits is challenging with traditional whole-building simulation tools, which require intensive calibration and struggle to model innovative system configurations and controls. To overcome these limitations, this study proposes a load profilebased retrofit analysis framework that provides better decisions under uncertainty. The main focus of this paper is the development of a probabilistic load profile model that can be used in the framework by using exploratory data analysis (EDA) of measured building data to properly quantify its inherent variability. A non-parametric Gaussian Process (GP) model was employed to capture the time- and weather-dependent characteristics of the heating load while explicitly modeling its uncertainty. The model's effectiveness is demonstrated through strong predictive performance on unseen data and physically interpretable insights into load behavior. This data-driven, probabilistic load profile serves as a robust and flexible input for subsequent system simulations, enabling a more confident and statistically sound analysis of retrofit potential.

Ham, S W↗
DOE OSTI2025

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran↗
DOE OSTI2026

Final Technical Report: Transport of Complex Mixtures in Ion-Containing Polymer Membranes

Permselective ion-containing membranes are an integral component for many applications from water treatment, fuel cells, and solar fuels devices where the selective transport of molecules and ions is desired. In solar fuels devices, ion-containing polymer membranes are responsible for permitting selective transport of ions between electrodes to maintain overall charge neutrality yet limit transport of reaction products produced at the electrodes. While the transport of single solutes through such membranes has been fairly well described, binary and multicomponent transport is poorly understood due to the myriad of interactions that occur in these systems (i.e. between co-permeants and between permeants and the membrane). Solar fuels devices are just one example of an application where understanding the transport of multiple simultaneous species is critically important to improving device performance as product crossover leads to reductions in overall device performance. The objectives of this research was to improve our understanding of the complex array of factors that influence transport behavior of multiple solutes within ion-containing polymer membranes. This experimental project addressed the lack of fundamental understanding of multicomponent transport behavior by synthesizing ion exchange membranes with varied incorporation of comonomers (ionic and neutral moieties) to investigate fundamental relationships between membrane structure, membrane physiochemical properties, and transport behavior of solutes and complex solute mixtures through dense, hydrated membranes.

25 ENERGY STORAGE↗
DOE OSTI2026

Superdiffusion resilience in Heisenberg chains with two-dimensional interactions on a quantum processor

Superdiffusive spin transport in the one-dimensional (1D) Heisenberg model is a key theoretical discovery in nonequilibrium quantum many-body physics. Although extensively studied in 1D systems, the breakdown and sustenance of superdiffusion in two-dimensional (2D) lattices with integrability-breaking terms, as found in real materials, remains an open question. To address this, we develop a toy model that extends the 1D Heisenberg model with a representative set of 2D interaction types and tunable strengths. Our model exhibits varying degrees of superdiffusion breakdown depending on the interaction type, spanning ballistic to diffusive regimes. We establish and justify a hierarchy of 2D interactions based on their resilience against superdiffusion breakdown: Heisenberg >𝑋⁢𝑋 > Ising. This precise control over the superdiffusive behavior also enables rigorous benchmarking of quantum hardware, and our simulations on IBM's Heron devices confirm the hardware's ability to accurately capture these many-body nonequilibrium phenomena. Overall, our results are relevant not only to simulating superdiffusion in real materials, such as the 1D Heisenberg compound KCuF3, which contains modest nonintegrable 2D terms, but also to extending superdiffusive behavior to larger 2D qubit lattices and other 2D materials.

Alagarsamy Manikandan, Keerthi Kumaran [ORNL]↗
DOE OSTI2026

Temperature as a control knob on spin-orbit coupling

Spin-orbit coupling (SOC) governs many physical phenomena. Through ab initio molecular dynamics simulations, Lu and Sun demonstrated that structural disorder at elevated temperatures substantially reduces the effective SOC contribution that stabilizes band inversion, driving a topological-to-normal-insulator transition in Bi2Se3. Their work identifies temperature as a meaningful control parameter for SOC-mediated topology and other properties.

Liang, Liangbo [ORNL] (ORCID:0000000311990049)↗