Stellar Interpretation of Meteoritic Data and PLotting for Everyone (SIMPLE): Isotope Mixing Lines for Six Sets of Core-collapse Supernova Models
Explore the source record for details and available documents.
Physics↗
SEARCH · Search NASA
Results for “numerical modeling”
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.
Dissipative Phase Transition in the Two-Photon Dicke Model
We explore the dissipative phase transition of the two-photon Dicke model, a topic that has garnered significant attention recently. Our analysis reveals that while single-photon loss does not stabilize the intrinsic instability in the model, the inclusion of two-photon loss restores stability, leading to the emergence of superradiant states, which coexist with the normal vacuum states. Using a second-order cumulant expansion for the photons, we derive an analytical description of the system in the thermodynamic limit, which agrees well with the exact calculation results. Additionally, we present the Wigner function for the system, shedding light on the breaking of the 𝑍4 symmetry inherent in the model. These findings offer valuable insights into stabilization mechanisms in open quantum systems and pave the way for exploring complex nonlinear dynamics in two-photon Dicke models.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Usage-based Lifing of Lithium-Ion Battery with HybridPhysics-Informed Neural Networks
Lithium-ion batteries are commonly used to power unmanned aircraft vehicles (UAVs).The ability to model and forecast the remaining useful life of these batteries enables UAV reliability assurance. Building accurate models for battery state of charge and state of health based on first principles is challenging due to the complex electrochemistry that governs battery operations and computational complexity required to solve them. Therefore, reduced order models are often used due to their ability to capture the overall battery discharge. Un-fortunately, these simplifications lead to residual discrepancy between model predictions and observed data. In this paper, we present a hybrid modeling approach merging reduced-order models and neural networks. In this approach, while most of the input-output relationship is captured by Nernst and Butler-Volmer equations, data-driven kernels reduce the gap between predictions and observations. We validate our approach using data publicly available through the NASA Prognostics Center of Excellence repository. Results showed that our hybrid battery prognosis model can be successfully calibrated, even with a limited number of observations.
Lithium-ion Battery↗
SimH 2 : an integrated techno-economic modeling framework for hydrogen pipeline infrastructure and network optimization
Large-scale hydrogen (H 2 ) pipeline transport design and network optimization have seldom been reported due to the lack of a cost model accounting for the relationship between transport cost and hydrogen mass flow rate. Here, this work introduced a system-level cost model for hydrogen pipeline transport at supercritical state and integrated it with an existing CO 2 pipeline network tool, SimCCS, for hydrogen-specific pipeline design and optimization. The Intermountain West (I-West) region of the U.S., historically dependent on fossil fuel-based economies, is chosen to demonstrate the capabilities of our H 2 pipeline cost model and transport network optimization platform called SimH 2 . Two scenarios are examined: one where the pipeline is not allowed to pass through disadvantaged communities and the other where it is permitted. The results highlight that incorporating disadvantaged-community constraints lead to longer pipeline routes and increased transport costs, reflecting the trade-offs involved in equitable infrastructure development. It is demonstrated that the newly developed SimH 2 tool not only enables the efficient design of H 2 transportation pipelines but also optimizes the network by accounting for local terrain and the presence of disadvantaged areas.
08 HYDROGEN↗
Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling
This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.
Fan, Joy [ORNL] (ORCID:0000000229751735)↗
Thermo-Mechanical Phase-Field Modeling of Fracture in High-Burnup UO2 Fuels Under Transient Conditions
This study presents a novel multiphysics phase-field fracture model to analyze high-burnup uranium dioxide (UO2) fuel behavior under transient reactor conditions. Fracture is treated as a stochastic phase transition, which inherently accounts for the random microstructural effects that lead to variations in the value of fracture strength. Moreover, the model takes into consideration the effects of temperature and burnup on thermal conductivity. Therefore, the model is able to predict crack initiation, propagation, and complex morphologies in response to thermal gradients and stress distributions. Several simulations were conducted to investigate the effects of operational and transient conditions on fracture behavior and the resulting cracking patterns. High-burnup fuels exhibit reduced thermal conductivity, elevating temperature gradients and resulting in extensive radial and circumferential cracks. Transient heating rates and temperatures significantly affect fracture patterns, with higher heating rates generating steeper gradients and more irregular crack trajectories. This approach provides critical insights into fuel integrity during accident scenarios and supports the safety evaluation of extended burnup limits.
Chemistry↗
Reverse Logistics Ev Battery Recycling Agent Base Model
This represents the initial regional tier of the electric vehicle (EV) battery recycling agent base model. Through this model, we can ascertain the number of EV purchases at both the state and regional levels. We employ census data to develop a diverse household profile to inform decisions regarding the acquisition of new or used EVs. The number of EV purchases at the state level will affect the future demand for recycling, reuse, and repurposing of end-of-life EV batteries. Additionally, tax credits, EV rebate programs, and the financial capacity of households will influence the number of EV purchases, thereby further impacting the demand for EV battery recycling.
Alam, Lamia [Idaho National Laboratory (INL), Idah↗
Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling
The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.
09 BIOMASS FUELS↗
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↗
A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane
The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.
Chemical looping↗
A multi-region approach for the analysis of porous materials and gas-solid reactions using USAXS-SAXS-WAXS: application to CaO carbonation
The reaction rate of gas-solid non-catalytic reactions is typically investigated using reactant conversion data over time and ex-situ measurements of the porous solid reactant textural properties; these data enable the experimental estimation of the initial intrinsic reaction rate, whereas the evolution of the textural properties and the reaction rate over time are evaluated theoretically using reaction models. In this work, a different methodology is presented, based on: a) in-situ time-resolved USAXS-SAXS-WAXS measurements of the solid sample, and b) a multi-region modeling approach; this methodology allows for the estimation of the textural properties and the intrinsic reaction rate at any time during the reaction. Referring to a reaction in which a solid product is obtained from a solid reactant and a gaseous reactant, the porous particle is described as consisting of several distinct regions with different microstructural properties, and both a gas-solid reaction model and a SAXS model are derived and applied to the carbonation of porous CaO. The proposed reaction model highlights the role of the inaccessible reactant and accurately predicts the solid reactant conversion versus time profile. The SAXS model accurately predicts the measured linear trends of the Porod invariant and of the pre-factor of the power law scattering profile versus CaO mass fraction; in the absence of chemical reactions, the proposed equations extend classical SAXS theory to porous materials containing macropores and nonporous solid phases in addition to standard nanoscale inhomogeneous regions.
CO2 solid sorbents↗
Improvements to RANS Modeling for Aeroheating Predictions on Blunt Bodies
Accurate predictions of aeroheating are critical for designing thermal protection systems for planetary entry vehicles. For larger vehicles, turbulence in the boundary layer can substantially increase convective heating. This turbulence must be accurately modeled to ensure the thermal protection system is sufficient. The majority of hypersonic turbulence model development and validation focuses on boundary layers developing over flat-plates or sharp cones; these cases are substantially different than the boundary layer that develops over the heatshield of a blunt body traveling at hypersonic speeds. Planetary missions often use blunt body geometries, such as the 70-degree sphere-cone favored by Mars missions or the 45-degree sphere-cone planned for the upcoming DAVINCI mission. Due to smaller vehicle size and the lower velocities in the stagnation region, planetary entry vehicles have relatively low Reynolds numbers. Surface curvature and high enthalpy gradients create additional challenges. These difficulties must be addressed to obtain high accuracy needed for the ambitious planetary missions in the upcoming decade. This work focuses on both assessing and improving Reynolds-averaged Navier-Stokes (RANS) turbulence models for blunt-body geometries typical of planetary entry vehicles, with a focus on one-equation and two-equation formulations.
Mars2020↗
Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices
Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.
2-dimensional systems↗
Complex orders and chirality in the classical Kitaev-Γ model
It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Reliability Models and Demonstration of a Fault-Tolerant Motor Concept for Vertical Takeoff and Landing Vehicles
This report documents the completion of the Revolutionary Vertical Lift Technology Project Annual Performance Indicator 24-3.2.4.1: “Apply and document reliability prediction for high reliability motor concept.” Two modeling tools were completed for calculation of reliability of fault-tolerant (FT) motors, and key FT operations of a modular FT motor were demonstrated experimentally. The two models are complementary tools for the stakeholder and user community. Both models employ Markov chain theory. The first model is a time-homogeneous Markov chain model, and the second is a time-inhomogeneous Markov-Weibull model. This report’s main sections are as follows: 1.0 Introduction, 2.0 Theory, 3.0 Motor Reliability Models, 4.0 Validation of FT Operation by Hardware Demonstration, and 5.0 Concluding Remarks. Novel contributions to the field include development of a modular FT motor concept for electrified vertical takeoff and landing (eVTOL) application, solution methods to solve the reliability calculations, development of figures of merit, and the introduction of “linked chains” to formulate a building-block approach for time-inhomogeneous Markov-Weibull modeling of motor reliability. Example case studies have been completed, and results are provided and discussed herein. A four-module FT motor concept was developed to a preliminary-design level of detail. This eVTOL FT motor concept was designed for galvanic, magnetic, and thermal isolation of stator winding faults. The reliability of the concept motor was calculated using a time-inhomogeneous Markov chain model. Employing average failure rate as a metric, 570 times greater reliability was achieved as compared to a baseline motor without fault tolerance. A demonstrator motor was built and tested. The testing demonstrated the key features of FT operation and validated the essential premises of the FT motor concepts presented herein. The experiments included successful demonstration of the feasibility of the following four key FT features: (1) terminal open-circuit operation, (2) thermal isolation after fault, (3) terminal short-circuit operation, and (4) internal short-circuit operation. These works indicate that FT modular motor drives offer promise for addressing the daunting reliability gap that electric aircraft propulsor drives are facing relative to the best conventional motor drive technology that is available today.
Electric Motor↗
Thermodynamically informed priors for uncertainty propagation in first-principles statistical mechanics
Here, this work demonstrates how first-principles statistical mechanics approaches within a Bayesian framework can quantify and propagate uncertainties to downstream thermodynamic calculations. To address the issue of Bayesian prior selection, knowledge of 0 K ground states in the material system of interest is incorporated into the prior. The effectiveness of this framework is shown by creating a phase diagram for the fcc zirconium nitride system, including confidence intervals on order-disorder transition temperatures.
Bayesian methods↗
Rethinking 𝛼−RuCl 3 : Parameters, models, and phase diagram
RuCl 3 was likely the first ever deliberately synthesized ruthenium compound, following the discovery of the 44 Ru element in 1844. For a long time it was known as an oxidation catalyst, with its physical properties being discrepant and confusing, until a decade ago when its allotropic form 𝛼−RuCl 3 rose to exceptional prominence. This “rediscovery” of 𝛼−RuCl 3 has not only reshaped the hunt for a material manifestation of the Kitaev spin liquid, but it has opened the floodgates of theoretical and experimental research in the many unusual phases and excitations that the anisotropic-exchange magnets as a class of compounds have to offer. Given its importance for the field of Kitaev materials, it is astonishing that the low-energy spin model that describes this compound and its possible proximity to the much-desired spin-liquid state is still a subject of significant debate ten years later. In the present study, we argue that the existing key phenomenological observations put strong natural constraints on the effective microscopic spin model of 𝛼−RuCl 3 , and specifically on its spin-orbit-induced anisotropic-exchange parameters that are responsible for the nontrivial physical properties of this material. These constraints allow one to focus on the relevant region of the multidimensional phase diagram of the 𝛼−RuCl 3 model, suggest an intuitive description of it via a different parametrization of the exchange matrix, offer a unifying view on the earlier assessments of its parameters, and bring closer together several approaches to the derivation of anisotropic-exchange models. We explore extended phase diagrams relevant to the 𝛼−RuCl 3 parameter space using quasiclassical, Luttinger-Tisza, exact diagonalization, and density-matrix renormalization-group methods, demonstrating a remarkably close quantitative accord between them on the general structure and hierarchy of the phases, with the zigzag, ferromagnetic, and incommensurate phases that are proximate to each other. As a result, one of the highlights is the detailed agreement on the nature of the incommensurate phases that realize two distinct counterrotating helical states.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Enabling Mission Flexibility to Battery Driven Deep Space Endeavors With Generalized Battery-Health-Monitoring Using Physics-Based and Data-Driven Reduced-Order Models
The needs and requirements for an electrochemical energy storage for deep space exploration is well explored. It is often understood that different mission sites and environmental conditions require different battery chemistries or technologies. Additionally, various engineering solutions are deployed to overcome specific chemical challenges. One often overlooked need is the “health” monitoring of an electrochemical storage system. The term generalized health monitoring, as envisioned in this work, refers to the monitoring of various aspects such as electrode health, electrolyte health, reaction pathway health, cooling system health, sensor health, and BMS health [1]. Generalized health monitoring allows mission leads, engineers, and scientists to incorporate flexibility in mission designs, make on-the-fly mission changes, and extend the duration of science missions. Moreover, it enables automation and data-driven decision-making without compromising safety and performance. Recently, our group developed a hierarchy of thermal reduced-order models (TROM) by combining a physics-based modeling approach and data-driven model reduction techniques applied to flight data [2]. The resulting TROMs were found to be not only accurate but also identifiable from the flight data. Consequently, the coefficient of variance of the model parameters is small over the course of hundreds of flights, allowing for monitoring the parameter evolution trajectories as the battery ages and degrades. These parameters constitute the metrics of the generalized health of a battery. Monitoring their evolution allows such models to be used for anomaly detection and prognostics, improving early detection of abnormal behavior and thus enabling timely maintenance, longer battery life, and enhanced battery safety. For this presentation, the practicality of the thermal model will be validated on a pack of 14cells under various topology configurations such as 1S14P, 2P7S, 7S2P, and 1P14S. It is well known that manufacturing and non-uniform aging lead to variability in the performance of a cell, which is exacerbated by cell balancing during active load. Additionally, in extreme scenarios, the paramount objective is to complete the mission, regardless of the stresses on the battery. Topology-induced balancing issues further stress the battery. The goal of this study is to determine if the noise (identifiability) in the reduced-order thermal model parameters is sensitive to topology, cell spacing, cooling strategy, and manufacturing or age variability. The variability in cells is considered by assuming a multimodal distribution for microscopic parameters of a cell (such as porosity, tortuosity, reaction kinetics, volumetric thermal conductivity, and volumetric heat capacity). The compounded effect of manufacturing variability, topological selection, cooling strategies, and cell balancing ages each cell in a battery differently. The study aims to clarify whether the challenge in extracting maximum information depends on the minimum number of sensors or models used for data extraction.
Automation↗