Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Modeling Single-Crystal Battery Materials: From Fundamental Understanding to Performance Evaluation

The performance of rechargeable batteries is fundamentally influenced by the physicochemical properties and microstructural features of their key material components. Recent experimental advancements have highlighted the potential of single-crystal (SC) morphologies to address inherent limitations of polycrystalline (PC) electrodes and solid-state electrolytes, offering tunable charge transport kinetics and improved cell cycling performance. Here, this review examines how state-of-the-art computational modeling, from atomistic and mesoscale to continuum-level approaches, including machine learning methodologies, has been utilized to investigate the critical factors governing the electrochemical behavior of SC battery materials. We explore how predictive modeling can elucidate the processing–structure–property–performance relationships of SC cathodes, anodes, and solid-state electrolytes, with a focus on unique SC characteristics such as crystallographic anisotropy, size effects, and facet-dependent properties. Additionally, we identify limitations in commonly used modeling techniques and discuss strategies to address these challenges. By integrating high-fidelity simulations with experimental insights, this review aims to outline a clear path for the rational design and optimization of SC battery components, paving the way for accelerated advancements in energy storage technologies.

Materials science↗
DOE OSTI2020

A thermodynamics-based damage model for the non-linear mechanical behavior of SiC/SiC ceramic matrix composites in irradiation and thermal environments

A damage model is developed and validated with experimental data for the non-linear mechanical behavior of SiC/SiC composite materials in nuclear applications. Cyclic thermal and mechanical loading associated with neutron irradiation effects of these composites leads to wide-spread and progressive micro-cracking that leads to loss of thermal conductivity and further enhancement of thermo-mechanical damage. A physics-based model of wide-spread micro-cracking is developed within the thermodynamic framework of continuum damage mechanics. Evolution equations for damage parameters that describe the growth of continuum damage are developed, where the material variables are obtained from experiments. The model novelty is in coupling mechanical, thermal, and irradiation damage through a consistent thermodynamic framework, including loss of thermal conductivity due to the evolution of mechanically induced micro-cracks. A number of thermo-mechanical experiments were conducted to confirm model assumptions. The model is shown to be validated with out-of-pile experiments, and then implemented using commercial finite element code COMSOL to the fuel cladding problem with normal and high radiation dose cases.

Materials Science↗
DOE OSTI2026

Modeling heat pipe startup and noncondensable gases in Sockeye

For this work, a one-dimensional gas mixture flow model was developed and implemented in the heat pipe code Sockeye to model the effects of noncondensable gases. Additionally, a startup model based on the dusty gas model was implemented to model the transition from rarefied gas dynamics to continuum flow, which occurs during the frozen startup of high-temperature heat pipes. Multiple startup and noncondensable gas models were tested against experimental data for sodium heat pipes, showing excellent agreement. Additionally, the newly developed gas mixture model for modeling noncondensable gas is further tested with a theoretical case study with arbitrary heating configurations. Finally, several recommendations and conclusions are made from the studies in this work to guide future heat pipe modeling efforts.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

LSM and CPT

We study a number of 1+1d lattice models with anti-unitary symmetries that simultaneously reflect space and reverse time. Some of these symmetries are anomalous, leading to Lieb-Schultz-Mattis-type constraints, thus excluding a trivially gapped phase. Examples include a mod 8 anomaly in the Majorana chain and various mod 2 anomalies in the spin chain. In some cases, there is an exact, non-anomalous lattice symmetry that flows in the continuum to CPT. In some other cases, the CPT symmetry of the continuum theory is emergent or absent. Depending on the model, the anomaly of the lattice model is matched in the continuum in different ways. In particular, it can be mapped to an emergent anomaly of an emanant symmetry.

Anomalies in Field and String Theories↗
DOE OSTI2020

Probing UV-sensitive Pathways for CN and HCN Formation in Protoplanetary Disks with the Hubble Space Telescope

The UV radiation field is a critical regulator of gas-phase chemistry in surface layers of disks around young stars. In an effort to understand the relationship between photocatalyzing UV radiation fields and gas emission observed at infrared and submillimeter wavelengths, we present an analysis of new and archival Hubble Space Telescope (HST), Spitzer, ALMA, IRAM, and SMA data for five targets in the Lupus cloud complex and 14 systems in Taurus-Auriga. The HST spectra were used to measure Lyα and far-UV (FUV) continuum fluxes reaching the disk surface, which are responsible for dissociating relevant molecular species (e.g., HCN, N{sub 2}). Semi-forbidden C ii] λ2325 and UV-fluorescent H{sub 2} emission were also measured to constrain inner disk populations of C{sup +} and vibrationally excited H{sub 2}. We find a significant positive correlation between 14 μm HCN emission and fluxes from the FUV continuum and C ii] λ2325, consistent with model predictions requiring N{sub 2} photodissociation and carbon ionization to trigger the main CN/HCN formation pathways. We also report significant negative correlations between submillimeter CN emission and both C ii] and FUV continuum fluxes, implying that CN is also more readily dissociated in disks with stronger FUV irradiation. No clear relationships are detected between either CN or HCN and Lyα or UV-H{sub 2} emission. This is attributed to the spatial stratification of the various molecular species, which span several vertical layers and radii across the inner and outer disk. We expect that future observations with the James Webb Space Telescope will build on this work by enabling more sensitive IR surveys than were possible with Spitzer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗
DOE OSTI2024

Fuzzy gauge theory for quantum computers

Continuous gauge theories, because of their bosonic degrees of freedom, have an infinite-dimensional local Hilbert space. Encoding these degrees of freedom on qubit-based hardware demands some sort of “qubitization” scheme, where one approximates the behavior of a theory while using only finitely many degrees of freedom. We propose a novel qubitization strategy for gauge theories, called “fuzzy gauge theory,” building on the success of the fuzzy σ -model in earlier work. We provide arguments that the fuzzy gauge theory lies in the same universality class as regular gauge theory, in which case its use would obviate the need of any further limit besides the usual spatial continuum limit. Furthermore, we demonstrate that these models are relatively resource-efficient for quantum simulations. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2023

The contribution of Compton ionization to ultrafast x-ray scattering

We investigate the role of Compton ionization in ultrafast non-resonant x-ray scattering using a molecular model system, which includes the ionization continuum via an orthonormalized plane wave ansatz. Elastic and inelastic components of the scattering signal, as well as coherent-mixed scattering that arises from electron dynamics, are calculated. By virtue of a near-quantitative distinction between scattering related to electronic transitions into bound and continuum states, we demonstrate how Compton ionization contributes to the coherent-mixed component. Analogous to inelastic scattering, the contribution to the coherent-mixed signal is significant and particularly manifests at intermediate and high-momentum transfers. Strikingly, for molecules with inversion symmetry, the exclusion of bound or continuum transitions may lead to the prediction of spurious coherent-mixed signals. We conclude that qualitative and quantitative accuracies of predicted scattering signals on detectors without energy resolution require that elements of the two-electron density operator are used. This approach inherently accounts for all accessible electronic transitions, including ionization.

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

Ramp-release experiments for strength measurements: Strain-rate dependence

This paper presents an enhanced analysis method for investigating material properties at high strain rates, extending the capability of established experimental techniques to gain more information. The ramp-release method has been applied to many experiments reported at high (≈10 5 − 10 6 s −1 ) strain-rates. More recent data gathered at the National Ignition Facility (NIF) has enabled higher (≈10 8 s −1 ) strain-rates to be studied. Here, we present an initial application of ramp-release analysis to NIF ramp-compression data, illustrating both the opportunities and the practical challenges of extending these methods to laser-driven platforms. The higher strain-rates accessed at the NIF mean that there is more strain-rate enhancement to strength, and the experimental configuration means that this enhancement is more readily seen in the data. This is enabled by the capability of avoiding peak-compression attenuation through the sample thickness with a designed hold period made possible by the pulse-shaping capability of NIF. We propose that this combination of experimental conditions and an enhanced analysis method enables the strain-rate enhancement to strength to be studied, and potentially for this to inform physics models at smaller scales than the continuum.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Low-energy physics of isotropic spin-1 chains in the critical and Haldane phases

Using a matrix product state algorithm with infinite boundary conditions, we compute high-resolution dynamic spin and quadrupolar structure factors in the thermodynamic limit to explore the low-energy excitations of isotropic bilinear-biquadratic spin-1 chains. Haldane mapped the spin-1 Heisenberg antiferromagnet to a continuum field theory, the nonlinear sigma model (NLσM). In this work, we find that the NLσM fails to capture the influence of the biquadratic term and provides only an unsatisfactory description of the Haldane phase physics. But several features in the Haldane phase can be explained by noninteracting multimagnon states. The physics at the Uimin-Lai-Sutherland point is characterized by multisoliton continua. Moving into the extended critical phase, we find that these excitation continua contract, which we explain using a field-theoretic description. New excitations emerge at higher energies and, in the vicinity of the purely biquadratic point, they show simple cosine dispersions. Using block fidelities, we identify them as elementary one-particle excitations and relate them to the integrable Temperley-Lieb chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A Continuum Dislocation Dynamics Crystal Plasticity Approach to Irradiated Body-Centered Cubic α-Iron

Radiation-induced embrittlement of reactor pressure vessel (RPV) steels can potentially limit the operating life of nuclear power plants. Over extended exposure to radiation doses, these body-centered cubic (BCC) irons demonstrate irradiation damage. Here, we present a continuum dislocation density (CDD) crystal plasticity model to capture the interaction among dislocations and self-interstitial atom (SIA) loops in α-iron. We demonstrate the importance of modeling cross slip using a combined stochastic Monte Carlo approach and the role of slip system strength anisotropy in capturing stochastic cross slip interactions. Through these captured interactions, the CDD crystal plasticity model can capture both the stress response and the physical evolution of dislocations on different slip system planes. Single-crystal verification experiments are used to calibrate the CDD crystal plasticity model, and a set of simplified polycrystalline simulations demonstrates the model’s ability to capture the stress response from tensile experiments on α-iron.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Catalytic resonance of ammonia synthesis by simulated dynamic ruthenium crystal strain

Ammonia affords dense storage for renewable energy as a fungible liquid fuel, provided it can be efficiently synthesized from hydrogen and nitrogen. In this work, the catalysis of ammonia synthesis was computationally explored beyond the Sabatier limit by dynamically straining a ruthenium crystal (±4%) at the resonant frequencies (10 2 to 10 5+ Hz) of N 2 surface dissociation and hydrogenation. Density functional theory calculations at different strain conditions indicated that the energies of NH x surface intermediates and transition states scale linearly, allowing the description of ammonia synthesis at a continuum of strain conditions. A microkinetic model including multiple sites and surface diffusion between step and Ru(0001) terrace sites of varying ratios for nanoparticles of differing size revealed that dynamic strain yields catalytic ammonia synthesis conversion and turnover frequency comparable to industrial reactors (400°C, 200 atm) but at lower temperature (320°C) and an order of magnitude lower pressure (20 atm).

10 SYNTHETIC FUELS↗
DOE OSTI2026

Continuum contribution to charged-current absorption of low-energy $ν_e$ on $^{40}$Ar

Accurate modeling of the absorption of tens-of-MeV $ν_e$ on $^{40}$Ar is needed to enable measurements of astrophysical neutrinos using large liquid argon time projection chamber (LArTPC) detectors, such as those planned for the Deep Underground Neutrino Experiment (DUNE). We revisit the MARLEY neutrino interaction model used in present estimates of DUNE sensitivity to supernova and solar neutrino signals. Multiple theoretical refinements are pursued, especially in the unbound continuum region of nuclear excitation energy. Inclusive charged-current neutrino-argon cross sections are calculated using a hybrid strategy. Nuclear transitions to unbound states are treated using a Hartree-Fock Continuum Random Phase Approximation (HF-CRPA) model, including forbidden contributions. Allowed transitions to low-lying discrete levels are also included using indirect measurements and approximate corrections for the momentum transfer dependence. Exclusive predictions are obtained by coupling these calculations with a statistical nuclear de-excitation model. The impact on observables of interest for DUNE and similar experiments is examined in terms of both total and differential cross sections. Our refined calculations predict a lower allowed portion of the cross section relative to the prior MARLEY model. At neutrino energies appreciably below 100 MeV, the inclusion of forbidden transitions does not fully compensate for the loss of allowed strength. For a representative neutrino burst from a galactic core-collapse supernova, our results suggest that MARLEY 1.2.0 overestimates the event yield in a DUNE-like detector by approximately 20%. However, because this overestimation is more severe at backwards angles, use of the charged-current $ν_e$-$^{40}$Ar reaction for supernova pointing may be more feasible than previously expected.

Gardiner, Steven [Fermilab]↗
DOE OSTI2021

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗
DOE OSTI2021

Perspective—Mass Conservation in Models for Electrodeposition/Stripping in Lithium Metal Batteries

Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. In this paper, using a motivating example of a moving boundary model in one dimension, we show how battery models need proper formulation for mass conservation, especially when simulated over multiple charge and discharge cycles. The article concludes with some thoughts on mass conservation and proper formulation for multiscale models.

25 ENERGY STORAGE↗
DOE OSTI2022

BattPhase—A Convergent, Non-Oscillatory, Efficient Algorithm and Code for Predicting Shape Changes in Lithium Metal Batteries Using Phase-Field Models: Part I. Secondary Current Distribution

Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. Thermodynamics-based equations have been used to study phase changes in lithium batteries using phase-field approaches. However, grid convergence studies and the effect of additional parameters needed to simulate these models are not well-documented in the literature. In this paper, using a motivating example of a moving boundary model in one- and two-dimensions, we show how one can formulate phase-field models, implement algorithms for the same and analyze the results. An open-access code with no restrictions is provided as well. This article concludes with some thoughts on the computational efficiency of phase-field models for simulating dendritic growth.

25 ENERGY STORAGE↗