Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state fidelity”

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 1,279 records · Page 71

CaloTrilogy: Toward a Breakthrough in One-Step, End-to-End, Physics-Guided Shower Generation for Modern Calorimeters

High-precision calorimeter simulation at current and future colliders imposes rapidly growing computational demands, motivating the development of machine-learning surrogates for traditional Monte Carlo tools such as Geant4. Flow matching and diffusion-based generative models have become leading approaches for high-dimensional fast simulation because of their sample quality, but typically require ${\cal O}(100)$ function evaluations at inference and often rely on auxiliary networks to constrain global observables, compromising streamlined end-to-end generation. We introduce a unified framework that improves the balance between speed, shower quality, and physics fidelity. The method combines: (i) an average velocity field integrator that enables sampling in one or a few evaluations; (ii) a learned generative prior in shower space, constructed from data rather than random noise; and (iii) physics-guided loss terms that impose inductive biases on key observables during training. These elements are training time regularizers, preserving end-to-end inference with no additional cost. With only one or a few evaluation steps, the model achieves shower quality competitive with state-of-the-art flow and diffusion approaches, tested on several public high granularity calorimeter datasets. The results demonstrate inter-layer shower structure consistent with the underlying physics, providing a strong candidate for future fast simulation workflows.

Jiang, Cheng [Edinburgh U.]↗

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. Here, this study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in OpenFOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

42 ENGINEERING↗

mphys-surrogate-model

This repository contains python scripts for building and studying reduced-order-modeling representations of droplet coalescence for eventual use in atmospheric models. The included data are generated from high-fidelity superdroplet methods and are utilized by machine learning pipelines to build data-driven models of droplet size distributions that evolve under coalescence. This repository further includes scripts to determine prediction (uncertainty) intervals on the data-driven model products based on conformal prediction.

Katona, JonasE [Lawrence Livermore National Labora↗

Battery data integrity and usability: Navigating datasets and equipment limitations for efficient and accurate research into battery aging

A tremendous commitment of resources is needed to acquire, understand and apply battery data in terms of performance and aging behavior. There are many state of performance (SOP) and state of health (SOH) metrics that are useful to guide alignment of batteries to end-use, yet how these metrics are measured or extracted can make the difference between usable, valuable datasets versus data that lacks the necessary integrity to meet baseline confidence levels for SOP/SOH quantification. This work will speak to 1) types of data that support SOP and SOH evaluations on mechanistic terms, 2) measurement conditions needed to assure high data integrity, 3) equipment limitations that can compromise data high fidelity, and 4) the impact of cell polarization on data quality. A common goal in battery research and field use is to work from a data platform that supports economical paths of data capture while minimizing down-time for battery diagnostics. An ideal situation would be to utilize data obtained during normal daily use (“pulses or cycles of convenience”) without stopping the daily duty cycles to perform dedicated SOP/SOH diagnostic routines. However, difficulties arise in trying to make use of daily duty cycle data (denoted as cycle-by-cycle, CBC) that underscores the need for standardization of conditions: temperature and duty cycles can vary over the course of a day and throughout a week, month and year; polarization can develop within an immediate cycle and throughout successive cycles as a hysteresis. If CBC data is envisioned as a data source to determine performance and aging trends, it should be recognized that polarization is a frequent consequence of CBC and thus makes it difficult to separate reversible and irreversible components to metrics such as capacity loss and resistance increase over aging. Since CBC conditions can have a major impact on data usability, we will devote part of this paper to CBC data conditioning and management. Differential analyses will also be discussed as a means to detect changing trends in data quality. Our target cell chemistries will be lithium-ion types NMC/graphite and LMO/LTO.

25 ENERGY STORAGE↗

Influence of World and Gravity Model Selection on Surface Interacting Vehicle Simulations

A vehicle simulation is surface-interacting if the state of the vehicle (position, velocity, and acceleration) relative to the surface is important. Surface-interacting simulations perform ascent, entry, descent, landing, surface travel, or atmospheric flight. Modeling of gravity is an influential environmental factor for surface-interacting simulations. Gravity is the free-fall acceleration observed from a world-fixed frame that rotates with the world. Thus, gravity is the sum of gravitation and the centrifugal acceleration due to the world s rotation. In surface-interacting simulations, the fidelity of gravity at heights above the surface is more significant than gravity fidelity at locations in inertial space. A surface-interacting simulation cannot treat the gravity model separately from the world model, which simulates the motion and shape of the world. The world model's simulation of the world's rotation, or lack thereof, produces the centrifugal acceleration component of gravity. The world model s reproduction of the world's shape will produce different positions relative to the world center for a given height above the surface. These differences produce variations in the gravitation component of gravity. This paper examines the actual performance of world and gravity/gravitation pairs in a simulation using the Earth.

Madden, Michael M.↗

Manipulating polymer composition to create low-cost, high-fidelity sensors for indoor CO2 monitoring

Abstract Carbon dioxide (CO 2 ) has been linked to many deleterious health effects, and it has also been used as a proxy for building occupancy measurements. These applications have created a need for low-cost and low-power CO 2 sensors that can be seamlessly incorporated into existing buildings. We report a resonant mass sensor coated with a solution-processable polymer blend of poly(ethylene oxide) (PEO) and poly(ethyleneimine) (PEI) for the detection of CO 2 across multiple use conditions. Controlling the polymer blend composition and nanostructure enabled better transport of the analyte gas into the sensing layer, which allowed for significantly enhanced CO 2 sensing relative to the state of the art. Moreover, the hydrophilic nature of PEO resulted in water uptake, which provided for higher sensing sensitivity at elevated humidity conditions. Therefore, this key integration of materials and resonant sensor platform could be a potential solution in the future for CO 2 monitoring in smart infrastructure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Dynamic Polarization Interactions in the Electrical Double Layer at Calcite (104) Interfaces with Aqueous Solutions

Reactions at mineral interfaces with aqueous solutions control many geochemical and biogeochemical processes in the Earth’s critical zone. At the molecular level, insights into important properties such as the structure of the electrical double layer (EDL) at specific mineral interfaces continue to improve due to the increasing fidelity of laboratory instrumentation and computational approaches. However, molecular simulation approaches suffer from limited reach into relevant scales of time, length, and system complexity. To span this gap, a novel hybrid approach that couples first principle plane-wave density functional theory (DFT) with classical DFT (cDFT) is demonstrated and applied to calcite (104) interfaces with various electrolytes. In this approach, a region of interest described using DFT interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. Benchmarking against experimental observations and entirely first principle DFT simulations demonstrates that this hybrid model efficiently encompasses the key short-range and collective interactions in the EDL. Simulations of calcite (104)/solution interfaces reveal the key static and dynamic polarization interactions that give rise to structuring of ions and water. Ion hydration interactions have the strongest effect on the depth of the first minimum in the density distribution of counterions at the surface, and the position and width of the first density peak is largely determined by the strength of ion-correlation forces. Finer details of ion distributions are controlled by mutual polarization of the calcite surface and interfacial electrolyte. Finally, this new ability to efficiently and rigorously predict EDL structure at mineral surfaces in contact with complex solutions paves the way to accurately modeling sorption, nucleation, dissolution, and growth in realistic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discovering nuclear models from symbolic machine learning

Numerous phenomenological nuclear models have been proposed to describe specific observables within different regions of the nuclear chart. However, developing a unified model that describes the complex behavior of all nuclei remains an open challenge. Here, we explore whether symbolic Machine Learning (ML) can rediscover traditional nuclear physics models or identify alternatives with improved simplicity, fidelity, and predictive power. To address this challenge, we developed a Multi-objective Iterated Symbolic Regression approach that handles symbolic regressions over multiple target observables, accounts for experimental uncertainties and is robust against high-dimensional problems. As a proof of principle, we applied this method to describe the nuclear binding energies and charge radii of light and medium mass nuclei. Our approach identified simple analytical relationships based on the number of protons and neutrons, providing interpretable models with precision comparable to state-of-the-art nuclear models. Additionally, we integrated this ML-discovered model with an existing complementary model to estimate the limits of nuclear stability. These results highlight the potential of symbolic ML to develop accurate nuclear models and guide our description of complex many-body problems.

Nuclear structure↗

Grain2Mesh: Mesh Generation for Grain-Scale Nonlinear Elasticity Modeling

The nonlinear hysteretic behavior of rocks under cyclic loading is a crucial area of study in geomechanics. The macroscopic response of a variety of materials has been found to be contingent upon the behavior of the micro-scale structure. This project aims to develop a functional and maintainable software package for generating a multi-phase numerical mesh and accompanying simulation files for finite element modeling used in computational mechanics solvers. Meshes generated from images often lack key preprocessing that reduces noise and prevents mesh element distortion that can increase computational cost. By incorporating user feedback throughout, grain2mesh ensures a high-fidelity mesh that can be used to model grain-scale interactions such as shearing, crack propagation, and interfacial material contrast. Scientific applications of this software include material fracturing, stress-strain analysis for natural and engineered materials, and nonlinear meso-scale analysis.

54 ENVIRONMENTAL SCIENCES↗

A Multi-Qubit Quantum Gate Using the Zeno Effect

The Zeno effect, in which repeated observation freezes the dynamics of a quantum system, stands as an iconic oddity of quantum mechanics. When a measurement is unable to distinguish between states in a subspace, the dynamics within that subspace can be profoundly altered, leading to non-trivial behavior. Here we show that such a measurement can turn a non-interacting system with only single-qubit control into a two- or multi-qubit entangling gate, which we call a Zeno gate. The gate works by imparting a geometric phase on the system, conditioned on it lying within a particular nonlocal subspace. We derive simple closed-form expressions for the gate fidelity under a number of non-idealities and show that the gate is viable for implementation in circuit and cavity QED systems. More specifically, we illustrate the functioning of the gate via dispersive readout in both the Markovian and non-Markovian readout regimes, and derive conditions for longitudinal readout to ideally realize the gate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Impact of Limited Degree of Freedom Drag Coefficients on a Floating Offshore Wind Turbine Simulation

The worldwide effort to design and commission floating offshore wind turbines (FOWT) is motivating the need for reliable numerical models that adequately represent their physical behavior under realistic sea states. However, properly representing the hydrodynamic quadratic damping for FOWT remains uncertain, because of its dependency on the choice of drag coefficients (dimensionless or not). It is hypothesized that the limited degree of freedom (DoF) drag coefficient formulation that uses only translational drag coefficients causes mischaracterization of the rotational DoF drag, leading to underestimation of FOWT global loads, such as tower base fore-aft shear. To address these hydrodynamic modeling uncertainties, different quadratic drag models implemented in the open-source mid-fidelity simulation tool, OpenFAST, were investigated and compared with the experimental data from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project. The tower base fore-aft shear and up-wave mooring line tension were compared under an irregular wave loading condition to demonstrate the effects of the different damping models. Two types of hydrodynamic quadratic drag formulations were considered: (1) member-based dimensionless drag coefficients applied only at the translational DoF (namely limited-DoF drag model) and (2) quadratic drag matrix model (in dimensional form). Based on the results, the former consistently underestimated the 95th percentile peak loads and spectral responses when compared to the OC5 experimental data. In contrast, the drag matrix models reduced errors in estimates of the tower base shear peak load by 7–10% compared to the limited-DoF drag model. The underestimation in the tower base fore-aft shear was thus inferred be related to mischaracterization of the rotational pitch drag and the heave motion/drag by the limited-DoF model.

17 WIND ENERGY↗

Engine Cycle Comparison for Alternative Propellant Nuclear Thermal Propulsion Engines

A modular Nuclear Thermal Propulsion (NTP) engine modeling suite was coded in Simulink to analyze various engine cycle configurations and propellants. This model was validated and used for previous studies involving alternative propellants for NTP engines. The current study compared hydrogen-based NTP (H-NTP) with ammonia-based Alternative propellant NTP (A-NTP) engines with expander and bleed cycle configurations. Based on prior work, the bleed cycle configuration was modified to yield a feasible result. The assumption that was made was that the Testing Reference Design (TRD) reactor will be used for all cases and that the difference in fuel loading fractions will be enough to offset any differences in the neutronics caused by switching the propellant. This work provided and compared propellant state points throughout the cycles as well as the key performance parameters. The conclusion was that bleed cycles for H-NTP engines are advantageous over expander cycles in terms of maximum system pressure, cycle simplicity, and potentially engine mass. Conversely, expander cycles are more advantageous for A-NTP engines according to the same scope of parameters. Future work will analyze reactor designs that are different from the TRD and use neutronics software integrated with the X-NTP model to provide higher fidelity results.

Nuclear↗

MFIX-Exa: Performance prediction of multiphase energy conversion devices

MFIX-Exa targets the high-fidelity CFD-DEM model in which particles are unresolved by the fluid grid (typically using a fluid mesh approximately twice the particle diameter) but discrete particle dynamics including collisions are resolved with a simple linear-spring dashpot model. The fluid field is solved with a modern low Mach number formation of a cell-centered, nodal-pressure approximate projection method using a Godunov scheme . Physics capabilities include: complex geometries via embedded boundaries (EBs), open and closed system ideal gas equations of state, species transport and heterogeneous chemi- cal reactions. MFIX-Exa is built on the AMReX software framework, the ECP Block-Structured Adaptive Mesh Refinement (AMR) Co-Design Center, which allows the code to be portable and performant.Further integration into the ECP ecosystem includes HYPRE linear solvers and Ascent for in situ visualization. Current project efforts focus on scaling up realistic simulations to the KPP of the challenge problem: NETs 50kW chemical looping reactor, shakedown and scaling on JSLE and OLCF TDS systems and reaching for any remaining performance improvements.

Musser, Jordan↗

Milestone in predicting core plasma turbulence: successful multi-channel validation of the gyrokinetic code GENE

On the basis of several recent breakthroughs in fusion research, many activities have been launched around the world to develop fusion power plants on the fastest possible time scale. In this context, high-fidelity simulations of the plasma behavior on large supercomputers provide one of the main pathways to accelerating progress by guiding crucial design decisions. When it comes to determining the energy confinement time of a magnetic confinement fusion device, which is a key quantity of interest, gyrokinetic turbulence simulations are considered the approach of choice – but the question, whether they are really able to reliably predict the plasma behavior is still open. The present study addresses this important issue by means of careful comparisons between state-of-the-art gyrokinetic turbulence simulations with the GENE code and experimental observations in the ASDEX Upgrade tokamak for an unprecedented number of simultaneous plasma observables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurements of extended magnetic fields in laser-solid interaction

Magnetic fields generated from a laser-foil interaction are measured with high fidelity using a proton radiography scheme with x-ray fiducials. In contrast to prior findings under similar experimental conditions, this technique reveals the self-generated, Biermann-battery fields extend beyond the edge of the expanding plasma plume to a radius of over 3.5 mm by t = + 1.4 ns. An analysis of two monoenergetic proton populations confirms that proton deflection is dominated by magnetic fields far from the interaction ( > 2 mm) and electric fields are insignificant. The results are not captured in state-of-the-art magnetohydrodynamics simulations and suggest the need to consider additional physics mechanisms for the magnetic field generation and transport in laser-solid interactions. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-Impulse to Time Optimal Finite Burn Trajectory Conversion

A novel conversion algorithm is presented that combines the fidelity of indirect optimization methods with the generality of direct methods to more easily solve time-optimal, finite-burn pseudo-rendezvous problems. An algorithm is described that converts a set of multiple-impulses, representing the entirety or a portion of a high- or low-thrust maneuver, to an exact time optimal finite-burn trajectory for a thrust limited, constant exhaust velocity spacecraft. A pseudo-rendezvous problem is one that yields a solution whose final time, position and velocity state is equal to that of the original post-impulsive trajectory. An iterative adjoint-control transformation is used to initialize the optimal control two-point boundary value problem. Examples are shown for both high and low-thrust non-coplanar Earth orbit transfers, as well as a low-thrust Hohmann-type Earth-Mars transfer.

Fogel, J.↗

Autonomous Navigation of a Lunar Relay Using GNSS and Other Measurements

Many of the highest priority destinations at the Moon lack a continuous view of Earth, such as the lunar poles or lunar far side. Exploration of these sites will require spacecraft in cislunar space to relay communications and provide position, navigation, and timing (PNT) services. Accurate knowledge of relay position, velocity, and time is essential to these services. This paper describes a concept for a PNT Instrument being developed for the Lunar Communications Relay and Navigation Systems (LCRNS) Project. The instrument is intended as a payload that would enable autonomous, on-board, real-time navigation and timing using Global Navigation Satellite System (GNSS), optical navigation, and one-way measurements from Earth-based ground stations. Hardware-in-the-loop simulations using flight software are used to realistically characterize performance on hardware platforms with a path to flight. These results provide preliminary validation of the proposed PNT Instrument, demonstrate the benefits of augmenting GNSS with other measurements, and serve as an insightful reference for the design of future lunar missions, including those that will operate within the LunaNet framework of standards. This instrument concept relies on several technologies developed at NASA Goddard Space Flight Center (GSFC). For GNSS observables, the instrument relies on the high-altitude NavCube 3 mini (NC3m) GNSS receiver specifically designed for cislunar applications. The autoNGC system, which consists of flight software and a hardware platform, is responsible for fusing the observables using its extended Kalman filter, the Goddard Enhanced Onboard Navigation System (GEONS). Optical navigation observables are processed within autoNGC (“autonomous Navigation, Guidance, and Control”) using the Goddard Image Analysis & Navigation Tool (GIANT) which is also responsible for simulating high-fidelity images for test and analysis. In addition to describing the PNT Instrument and its components, the paper will present predicted performance based on simulation results. As a baseline, it will present GNSS-only hardware-in-the loop results using a NC3m test unit to process Spirent-simulated GPS signals in a potential lunar relay trajectory: a 12-hour elliptical frozen lunar orbit (ELFO). GEONS then processes the GPS pseudorange and time differenced carrier phase measurements to estimate and propagate the relay state (position, velocity, and time). These results extend previously published work that showed preliminary ELFO performance. Previous work has shown the importance of other measurement types, so additional simulations are performed which augment GNSS with ground station observables and several methods of optical navigation, including celestial navigation, limb-finding (e.g., observations of the lunar horizon), and terrain relative navigation (TRN). TRN involves correlating simulated predicted images of the lunar surface with actual imagery; misalignments of landmarks identified in each image are translated into relay state updates. TRN is valuable as a measurement of the relay’s state relative to the Moon, especially during GNSS outages or after maneuvers. One-way Pseudorange and Doppler measurements from Earth-based ground stations are also simulated. The full set of observables is processed using autoNGC. These simulations make use of autoNGC and NC3m test units, a lab atomic clock, and a pulse-per-second (PPS) generation and distribution system. This combination of subsystems, and the hardware platforms used in this analysis, represents a PNT Instrument that could be flown on a lunar relay. Results from the hardware-in-the-loop simulations presented in this paper provide a preliminary assessment of the achievable navigation performance of this instrument concept. PNT Instrument performance is compared to the GPS-only performance, and a discussion is provided on the apparent merits and challenges of each measurement type.

Ben Ashman↗