Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗
DOE OSTI2020

Stabilization of Low Valent Zirconium Nitrides in Titanium Nitride via Plasma-Enhanced Atomic Layer Deposition and Assessment of Electrochemical Properties

The increasing awareness of metal nitrides as feasible electrocatalysts in numerous energy-relevant transformations motivates greater synthetic control of mixed metal nitride compositions and their valence states. We report plasma-enhanced atomic layer deposition of titanium zirconium nitride, with an emphasis on the experimental parameters that influence the valence states of the two metal ions. We further examined the consequence of various ratios of Ti and Zr on the electrochemically active surface areas (ECSAs). Films with a Zr-to-Ti ratio of 1:2 or 1:1 yielded the highest ECSA; deviation from these values led to a large fraction of electrically insulating or semiconducting titanium or zirconium oxynitrides/oxides. The systematic examination establishes a fundamental structure-property relationship critical to assessing the electrochemical/electrocatalytic properties of metal nitrides.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Discretization limits of lattice‐Boltzmann methods for studying immiscible two‐phase flow in porous media

Summary Digital images of porous media often include features approaching the image resolution length scale. The behavior of numerical methods at low resolution is therefore important even for well‐resolved systems. We study the behavior of the Shan‐Chen (SC) and Rothman‐Keller (RK) multicomponent lattice‐Boltzmann models in situations where the fluid‐fluid interfacial radius of curvature and/or the feature size of the medium approaches the discrete unit size of the computational grid. Various simple, small‐scale test geometries are considered, and a drainage test is also performed in a Bentheimer sandstone sample. We find that both RK and SC models show very high ultimate limits: in ideal conditions the models can simulate static fluid configuration with acceptable accuracy in tubes as small as three lattice units across for RK model (six lattice units for SC model) and with an interfacial radius of curvature of two lattice units for RK and SC models. However, the stability of the models is affected when operating in these extreme discrete limits: in certain circumstances the models exhibit behaviors ranging from loss of accuracy to numerical instability. We discuss the circumstances where these behaviors occur and the ramifications for larger‐scale fluid displacement simulations in porous media, along with strategies to mitigate the most severe effects. Overall we find that the RK model, with modern enhancements, exhibits fewer instabilities and is more suitable for systems of low fluid‐fluid miscibility. The shortcomings of the SC model seem to arise predominantly from the high, strongly pressure‐dependent miscibility of the two fluid components.

Li, Zhe↗
DOE OSTI2023

Bound-preserving finite element approximations of the Keller–Segel equations

We report this paper aims to develop numerical approximations of the Keller–Segel equations that mimic at the discrete level the lower bounds and the energy law of the continuous problem. We solve these equations for two unknowns: the organism (or cell) density, which is a positive variable, and the chemoattractant density, which is a non-negative variable. We propose two algorithms, which combine a stabilized finite element method and a semi-implicit time integration. The stabilization consists of a nonlinear artificial diffusion that employs a graph-Laplacian operator and a shock detector that localizes local extrema. As a result, both algorithms turn out to be nonlinear and can generate cell and chemoattractant numerical densities fulfilling lower bounds. However, the first algorithm requires a suitable constraint between the space and time discrete parameters, whereas the second one does not. We design the latter to attain a discrete energy law on acute meshes. We report some numerical experiments to validate the theoretical results on blowup and nonblowup phenomena. In the blowup setting, we identify a locking phenomenon that relates the L ∞ (Ω)-norm to the L 1 (Ω)-norm limiting the growth of the singularity when supported on a macroelement.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Dynamics of Intra-Cell Thermal Front Propagation in Lithium-Ion Battery Safety Issues

Thermal runaway (TR), a critical failure mode in lithium-ion batteries (LIBs), poses significant safety risks and hinders wider application of LIBs. TR typically begins at a localized heat source and spreads across the cell. Understanding thermal front propagation (TFP) characteristics, such as front and velocity, is crucial for assessing energy release and temperature distribution for battery hazardous estimation. Recent studies assume that TR within cells propagates at a near-constant velocity, based on the reaction kinetics and thermal properties. Here, in this study, an intra-battery TR model is further proposed and it indicates that TFP velocity stabilizes when the front is distanced from the heat source. Theoretical estimates for propagation velocity and front are developed and validated through numerical simulations and experimental tests from the NREL Battery Failure Databank. The energy release rate during TFP and the impact of preheating based on a point heat source are explored. This work clarifies the long-standing clouds of the thermal font propagation behaviors within the single cell, highlights the power and beauty of mathematics modeling to describe the complicated thermal behaviors, and provides important guidelines for thermal hazardous understanding for next-generation batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Reversible structural transformations in supercooled liquid water from 135 to 245 K

Water has many anomalous properties compared to “simple” liquids, and these anomalies are typically enhanced in supercooled water.1-3 While numerous models have been proposed, including the liquid-liquid critical point,4,5 the singularity-free scenario,6 and the stability limit conjecture,1 a molecular-level understanding remains elusive. The main difficulty in determining which, if any, of these models is correct is the limited amount of data in the relevant temperature and pressure ranges. For water at ambient pressures, which is the focus of this work, data is largely missing from 160 – 232 K (“No Man’s Land”) due to rapid crystallization.2,3 Whether rapid crystallization is just an experimental obstacle, or a fundamental problem signaling the inability of water to thermally equilibrate prior to crystallization is also a major unanswered question.5,7 Here, we investigate the structural transformations of transiently-heated, supercooled water with nanosecond time resolution using infrared vibrational spectroscopy. The experiments demonstrate three key results. First, water’s structure relaxes from its initial configuration to a “steady-state” configuration prior to the onset of crystallization over a wide temperature range. Second, water’s steady-state structure can be reproduced by a linear combination of two, temperature-independent structures that correspond to a “high-temperature liquid” and a “low-temperature liquid.” Third, the observed structural changes are reversible over the full temperature range. Taken together, these results show that supercooled water can equilibrate prior to crystallization for temperatures from the homogeneous nucleation temperature, TH ~232 K,3 down to the glass transition temperature (Tg ~ 136 K). Second, the results provide support for the hypothesis that supercooled water can be described as a mixture of two, structurally-distinct, interconvertible liquids from 135 K to 245 K.5,8-18

Kringle, Loni M.↗
DOE OSTI2022

Cyber-Physical Secure Observer-Based Corrective Control under Compromised Sensor Measurements

In this paper we introduce the objectives and design principles of corrective control under cyber-physical attacks. We propose two types of observer-based corrective control for both the open-loop stable and the open-loop unstable LTI systems. The basic idea of our corrective control design is to use the observer as the ground-truth during the attack, making the plant dynamics follow the observer behavior. This is the opposite to the no-attack-detected period in which the observer is designed to follow the plant dynamics. We show stability of the proposed control under compromised sensor measurements, and quantify the effects of the discrepancy between the observer and the plant. Numerical examples, with illustrations using microgrid energy dynamics, are presented to show benefits of the proposed corrective control.

Wu, Dan↗
DOE OSTI2020

Multi-frequency General Relativistic Radiation-hydrodynamics with M 1 Closure

We report on recent upgrades to our general relativistic radiation-magnetohydrodynamics code, Cosmos++ , which expands the two-moment, M 1 , radiation treatment from gray to multi-frequency transport, including Doppler and gravitational frequency shifts. The solver accommodates either photon (Bose–Einstein) or neutrino (Fermi–Dirac) statistical distribution functions with absorption, emission, and elastic scattering processes. An implicit scheme is implemented to simultaneously solve the primitive inversion problem together with the radiation–matter coupling source terms, providing stability over a broad range of opacities and optical depths where the interaction terms can be stiff. We discuss our formulations and numerical methods, and validate our methods against a wide variety of test problems spanning optically thin to thick regimes in flat, weakly curved, and strongly curved spacetimes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗
DOE OSTI2020

Direct observation of the formation and stabilization of metallic nanoparticles on carbon supports

Direct formation of ultra-small nanoparticles on carbon supports by rapid high temperature synthesis method offers new opportunities for scalable nanomanufacturing and the synthesis of stable multi-elemental nanoparticles. However, the underlying mechanisms affecting the dispersion and stability of nanoparticles on the supports during high temperature processing remain enigmatic. In this work, we report the observation of metallic nanoparticles formation and stabilization on carbon supports through in situ Joule heating method. We find that the formation of metallic nanoparticles is associated with the simultaneous phase transition of amorphous carbon to a highly defective turbostratic graphite (T-graphite). Molecular dynamic (MD) simulations suggest that the defective T-graphite provide numerous nucleation sites for the nanoparticles to form. Furthermore, the nanoparticles partially intercalate and take root on edge planes, leading to high binding energy on support. This interaction between nanoparticles and T-graphite substrate strengthens the anchoring and provides excellent thermal stability to the nanoparticles. These findings provide mechanistic understanding of rapid high temperature synthesis of metal nanoparticles on carbon supports and the origin of their stability.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Orbital stability analysis and photometric characterization of the second Earth Trojan asteroid 2020 XL5

Trojan asteroids are small bodies orbiting around the L 4 or L 5 Lagrangian points of a Sun-planet system. Due to their peculiar orbits, they provide key constraints to the Solar System evolution models. Despite numerous dedicated observational efforts in the last decade, asteroid 2010 TK 7 has been the only known Earth Trojan thus far. Here we confirm that the recently discovered 2020 XL 5 is the second transient Earth Trojan known. To study its orbit, we used archival data from 2012 to 2019 and observed the object in 2021 from three ground-based observatories. Our study of its orbital stability shows that 2020 XL 5 will remain in L 4 for at least 4 000 years. With a photometric analysis we estimate its absolute magnitude to be ${H}_{r}=18.5{8}_{-0.15}^{+0.16}$, and color indices suggestive of a C-complex taxonomy. Assuming an albedo of 0.06 ± 0.03, we obtain a diameter of 1.18 ± 0.08 km, larger than the first known Earth Trojan asteroid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A novel defoamer for processing nuclear waste: Testing and performance

Abstract Legacy radioactive waste from nuclear weapons material processing is the biggest environmental challenge in the state of South Carolina, and one of the biggest in the United States. Although substantial progress has been made in processing and vitrifying high‐level radioactive waste at the Savannah River site, approximately 35 million gallons remains to be treated and dispositioned. In this article, we show the development of a new defoamer for use in the processing of high‐level radioactive waste. The application of the new defoamer is not only more effective at controlling foam but will shorten the batch processing time leading to safer and more efficient processing. This has the potential to help shorten the overall site mission lifetime, saving the federal government hundreds of millions of dollars. The novel aspect of the new defoamer is its effectiveness and stability in the harsh conditions needed for processing high‐level radioactive waste: high temperatures, pH 4–13, strong oxidizing and reducing agents, and numerous metal catalysts. The defoamer will replace an existing antifoam agent that decomposes to form multiple flammable gases. Because the existing antifoam agent decomposes quickly, processing changes increased processing time. Statement of novelty We demonstrate that a new defoamer has been developed for use in processing high‐level radioactive waste. This novel defoamer is effective and stable in the harsh processing conditions with no detectable flammable by‐products. Because the existing antifoam agent decomposes quickly, processing changes, including lowering acid addition rates and evaporation rates, increased processing time. The new defoamer is not only more effective at controlling foam but will shorten the batch processing time and lead to safer and more efficient processing by eliminating the production of flammable gases. Statement of industrial relevance Foaming is of great concern in many industrial processes involving three‐phase gas/liquid/fine‐solids systems, like in water evaporation, mineral floatation, air sparging in situ remediation techniques, and in those found in the paper industry.

Lambert, Dan P.↗
DOE OSTI2023

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗
DOE OSTI2023

Machine learning assisted derivation of minimal low-energy models for metallic magnets

Abstract We consider the problem of extracting a low-energy spin Hamiltonian from a triangular Kondo Lattice Model (KLM). The non-analytic dependence of the effective spin-spin interactions on the Kondo exchange excludes the use of perturbation theory beyond the second order. We then introduce a Machine Learning (ML) assisted protocol to extract effective two- and four-spin interactions. The resulting spin model reproduces the phase diagram of the original KLM as a function of magnetic field and single-ion anisotropy and reveals the effective four-spin interactions that stabilize the field-induced skyrmion crystal phase. Moreover, this model enables the computation of static and dynamical properties with a much lower numerical cost relative to the original KLM. A comparison of the dynamical spin structure factor in the fully polarized phase computed with both models reveals a good agreement for the magnon dispersion even though this information was not included in the training data set.

Sharma, Vikram (ORCID:0000000189105519)↗
DOE OSTI2022

Atomically dispersed iron sites with a nitrogen–carbon coating as highly active and durable oxygen reduction catalysts for fuel cells

Nitrogen-coordinated single atom iron sites (FeN 4 ) embedded in carbon (Fe–N–C) are the most active platinum group metal-free oxygen reduction catalysts for proton-exchange membrane fuel cells. Still, current Fe–N–C catalysts lack sufficient long-term durability and are not yet viable for practical applications. Here we report a highly durable and active Fe–N–C catalyst synthesized using heat treatment with ammonia chloride followed by high-temperature deposition of a thin layer of nitrogen-doped carbon on the catalyst surface. We propose that catalyst stability is improved by converting defect-rich pyrrolic N-coordinated FeN 4 sites into highly stable pyridinic N-coordinated FeN 4 sites. The stability enhancement is demonstrated in membrane electrode assemblies using accelerated stress testing and a long-term steady-state test (>300 h at 0.67 V), approaching a typical Pt/C cathode (0.1 mg Pt cm -2 ). The encouraging stability improvement represents a critical step in developing viable Fe–N–C catalysts to overcome the cost barriers of hydrogen fuel cells for numerous applications.

08 HYDROGEN↗
DOE OSTI2024

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Molecular Brush Surfactants: Versatile Emulsifiers for Stabilizing and Structuring Liquids

Abstract Using amphiphilic molecular brushes to stabilize emulsions usually requires the synthesis of specific side chains, which can be a time‐consuming and difficult challenge to meet. By taking advantage of the electrostatic interactions between water‐soluble molecular brushes and oil‐soluble oligomeric ligands, the in situ formation, assembly and jamming of molecular brush surfactants (MBSs) at the oil‐water interface is described. With MBSs, stable emulsions including o/w, w/o and o/w/o can be easily prepared by varying the molar ratios of the molecular brushes to the ligands. Moreover, when jammed, the binding energy of MBSs at the interface is sufficiently strong to allow the stabilization of liquids in nonequilibrium shapes, i.e., structuring liquids, producing an elastic film at the interface with exceptional mechanical properties. These structured liquids have numerous potential applications, including chemical biphasic reactions, liquid electronics, and all‐liquid biomimetic system.

Wang, Beibei↗