Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Order parameters”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Thermal Data-driven Model Reduction for Enhanced Battery Health Monitoring

Electric aviation faces a major challenge of avoiding potentially catastrophic consequences of the battery’s thermal runaway while keeping the weight of the battery low. Detection of early warning signals of battery failures requires accurate monitoring of the battery’s health throughout its lifespan. However, identifying the parameters of the battery from field data is notoriously difficult. We investigate this problem within the framework of modeling the temperature dynamics of a Li-ion cell during tests simulating loading in electric aircraft flights. It is found that the parameters of a higher-fidelity physics-based thermal model cannot be identified from the simulated flight data. To resolve this issue, we reduce the higher-fidelity thermal model to a model with fewer parameters. The resulting reduced-order model can predict temperature dynamics accurately and is identifiable throughout the cell’s lifespan which allows using the model’s parameters to monitor the state-of-health of the aging cell and detect anomalies in thermal behavior.

Li ion batteries↗
NASA NTRS1974

Second-order small disturbance theory for hypersonic flow over power-law bodies

A mathematical method for determining the flow field about power-law bodies in hypersonic flow conditions is developed. The second-order solutions, which reflect the effects of the second-order terms in the equations, are obtained by applying the method of small perturbations in terms of body slenderness parameter to the zeroth-order solutions. The method is applied by writing each flow variable as the sum of a zeroth-order and a perturbation function, each multiplied by the axial variable raised to a power. The similarity solutions are developed for infinite Mach number. All results obtained are for no flow through the body surface (as a boundary condition), but the derivation indicates that small amounts of blowing or suction through the wall can be accommodated.

Townsend, J. C.↗
DOE OSTI2025

Imaging Bias‐Driven Domain Wall Motion With Scanning Oscillator Piezoresponse Force Microscopy

Understanding ferroelectric domain wall dynamics at the nanoscale across a broad range of timescales requires measuring domain wall position under different applied electric fields. The success of piezoresponse force microscopy (PFM) as a tool to apply local electric fields at different positions and imaging their changing position, together with the information obtained from associated switching spectroscopies has fueled numerous studies of the dynamics of ferroelectric domains to determine the impact of intrinsic parameters such as crystalline order, defects and pinning centers, as well as boundary conditions such as environment. However, the investigation of sub-coercive reversible domain wall vibrational modes requires the development of new tools that enable visualizing domain wall motion under varying applied fields with high temporal and spatial resolution while also accounting for spurious electrostatic effects. Here, scanning oscillator piezoresponse force microscopy extends the investigation of domain wall dynamics to new regimes, providing direct visualization of domain wall position as a function of an external electric field that varies in time and location. In conclusion, this enables studying the energetics of field-driven ferroelectric domain wall motion, which is shown to obey a thermally activated flow regime in the millisecond timescale.

Piezoresponse Force Microscopy↗
DOE OSTI2024

Lifting of two-mode states in the D1-D5 CFT

We consider D1-D5-P states in the untwisted sector of the D1-D5 orbifold CFT where one copy of the seed CFT has been excited by a pair of oscillators, each being either bosonic or fermionic. While such states are BPS at the orbifold point, they will in general ‘lift’ as the theory is deformed towards general values of the couplings. We compute the expectation value of this lift at second order in the deformation parameter for the above mentioned states. We write this lift in terms of a fixed number of nested contour integrals on a given integrand; this integrand depends on the mode numbers of the oscillators in the state. We evaluate these integrals to obtain the explicit value of the lift for various subfamilies of states. At large mode numbers one observes a smooth increase of the lift with the dimension of the state h; this increase appears to follow a ~ √h behavior similar to that found analytically in earlier computations for other classes of states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Lifting of superconformal descendants in the D1-D5 CFT

We consider D1-D5-P states in the untwisted sector of the D1-D5 orbifold CFT where we excite one copy of the seed CFT with a left-moving superconformal descendant. When the theory is deformed away from this region of moduli space these states can ‘lift’, despite being BPS at the orbifold point. For descendants formed from the supersymmetry $^{Gα}_{Å,–s}$ and R-symmetry $^{Ja}_{–n}$ current modes we obtain explicit results for the expectation value of the lifts for various subfamilies of states at second order in the deformation parameter. A smooth ~ √h behaviour is observed in the lifts of these subfamilies for large dimensions. Using covering space Ward identities we then find a compact expression for the lift of the above $^{Ja}_{–n}$ descendant states valid for arbitrary dimensions. In the large-dimension limit this lift scales as ~ √h, strengthening the conjecture that this is a universal property of the lift of D1-D5-P states. We observe that the lift is not simply a function of the total dimension, but depends on how the descendant level is partitioned amongst modes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Interactions of two and three mesons including higher partial waves from lattice QCD

We study two- and three-meson systems composed either of pions or kaons at maximal isospin using Monte Carlo simulations of lattice QCD. Utilizing the stochastic LapH method, we are able to determine hundreds of two- and three-particle energy levels, in nine different momentum frames, with high precision. We fit these levels using the relativistic finite-volume formalism based on a generic effective field theory in order to determine the parameters of the two- and three-particle K-matrices. We find that the statistical precision of our spectra is sufficient to probe not only the dominant s-wave interactions, but also those in d waves. In particular, we determine for the first time a term in the three-particle K-matrix that contains two-particle d waves. We use three N f = 2 + 1 CLS ensembles with pion masses of 200, 280, and 340 MeV. This allows us to study the chiral dependence of the scattering observables, and compare to the expectations of chiral perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Universal lifting in the D1-D5 CFT

We consider D1-D5-P states in the untwisted sector of the D1-D5 orbifold CFT where one copy of the seed CFT has been excited with a left-moving superconformal primary. Despite being BPS at the orbifold point, such states can ‘lift’ as the theory is deformed away from this point in moduli space. We compute this lifting at second order in the deformation parameter for arbitrary left-moving dimension h of this class of states. This result displays an interesting universality since the lifting does not depend on the details of the superconformal primary; it depends only on the dimension. In the large-dimension limit the lift scales as $\sqrt{h}$; it is observed that such scaling appears to be a universal property of the lift of D1-D5-P states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Exceptional magnetic and magnetoelastic behavior of rare-earth non-centrosymmetric Sm 7 Pd 3

Magnetic compounds possessing an intrinsic combination of near-zero magnetization with high magnetic anisotropy are highly desirable for spinronic applications and memory recording. A comprehensive study of Sm 7 Pd 3 binary compound uncovered a unique combination of strong magnetoelastic behavior, very low net magnetization, and exceptionally high magnetic coercivity. The temperature-dependent X-ray synchrotron powder diffraction study indicates the abrupt changes in the compound's lattice parameters at the magnetic ordering temperature of T C =169 K, although the crystal structure remains non-centrosymmetric hexagonal Th 7 Fe 3 -type down to 6 K. Density functional theory calculations confirm high intrinsic magnetocrystalline anisotropy of Sm 7 Pd 3 , which explains the extremely large coercivity of the polycrystalline sample, up to H cr = 130 kOe at 2 K. In conclusion, this discovery brings to life a novel class of highly anisotropic materials that are distinctly different from known spintronic materials, making them interesting future systems for magnetic memory research.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Higher-order interaction effects among operating conditions and feedstocks shape reactor microbiomes and fatty acid production profiles

Arrested anaerobic digestion (AAD) offers a promising route for producing fatty acids (FAs) from organic residues, yet optimal conditions for selectively generating medium-chain fatty acids (MCFAs) remain poorly defined. Here, we systematically evaluated the main and interaction effects of pH (5, 7, 9), feedstock (food waste, manure), temperature (35 and 45 °C), and inoculum source on microbiome composition and FA production. Anaerobic digester sludge and a novel bison rumen inoculum were compared. Significant higher-order interactions among operating parameters governed FA profiles and microbiome structure. Butyric acid production was driven by a three-way interaction among pH, feedstock, and temperature (p < 0.001), with maximum concentrations achieved in food waste reactors at pH 5.0 and 35 °C (1.2 ± 0.1 g L −1 with sludge and 1.1 ± 0.3 g L −1 with rumen). MCFA production exhibited significant four-way interactions (p < 0.1 to p < 0.001). At 45 °C and pH 5.0, inoculum source tuned MCFA selectivity: sludge favored pentanoic acid (0.4 ± 0.1 g L −1 ), whereas rumen favored hexanoic and heptanoic acids (up to 0.4 ± 0.2 g L −1 ). Manure reactors produced < 0.2 g L −1 MCFAs under all conditions. Genera, including Megasphaera, Prevotella, and Lactobacillus, were associated with production of specific MCFAs. PICRUSt2-based pathway predictions were consistent with MCFA production patterns and suggested a potential role for lactic acid–driven chain elongation pathways. This study provides insights into how interacting operating conditions shape AAD microbiomes, their FA profiles, and advances the trajectory of research aimed at engineering robust and controllable microbiomes for waste valorization.

09 BIOMASS FUELS↗
DOE OSTI2025

Advantages of the Samarskii-type schemes on the Shishkin mesh

The schemes of the Samarskii type are simple modifications of the upwind scheme. We use them on the Shishkin mesh and discuss their advantages over the upwind scheme when applied to the linear one-dimensional singularly perturbed convection–diffusion problem. One of the advantages is that the Samarskii-type schemes have exact first-order accuracy uniform in the perturbation parameter, as opposed to the upwind scheme which is almost first-order uniformly accurate because its accuracy is diminished by logarithmic factors. Although this is not a new result, we re-emphasize it in the paper. We also demonstrate another advantage, that the Samarskii-type schemes are almost second-order uniformly accurate on the layer component of the solution. Motivated by this fact, we present a further improvement of the numerical method.

Convection–diffusion↗
DOE OSTI2022

Goal-oriented a-posteriori estimation of model error as an aid to parameter estimation

In this work, a Bayesian model calibration framework is presented that utilizes goal-oriented a-posterior error estimates in quantities of interest (QoIs) for classes of high-fidelity models characterized by PDEs. It is shown that for a large class of computational models, it is possible to develop a computationally inexpensive procedure for calibrating parameters of high-fidelity models of physical events when the parameters of low-fidelity (surrogate) models are known with acceptable accuracy. The main ingredients in the proposed model calibration scheme are goal-oriented a-posteriori estimates of error in QoIs computed using a so-called lower fidelity model compared to those of an uncalibrated higher fidelity model. The estimates of error in QoIs are used to define likelihood functions in Bayesian inversion analysis. A standard Bayesian approach is employed to compute the posterior distribution of model parameters of high-fidelity models. As applications, parameters in a quasi-linear second-order elliptic boundary-value problem (BVP) are calibrated using a second-order linear elliptic BVP. In a second application, parameters of a tumor growth model involving nonlinear time-dependent PDEs are calibrated using a lower fidelity linear tumor growth model with known parameter values.

A-posterior estimates↗
DOE OSTI2024

Assessing the Sensitivity of Pourbaix Diagrams to Computational Protocols: Electrochemical Stability of Ni Oxides as a Case Study

Pourbaix diagrams stand as a useful tool in assessing and visualizing materials’ electrochemical stability and are widely used for electrocatalyst design. However, their reliability hinges on the accuracy of the chemical potentials of involved phases, which may bear uncertainties and can be significantly impacted by decision-making steps in the computational protocol. Here, this study introduces a robust sensitivity analysis framework, exemplified through a detailed examination of the computational Pourbaix diagram of Ni, the oxides of which are used as high-activity and cost-friendly catalysts for many electrochemical reactions. Quantities of interest derived from the Pourbaix diagram include the appearance and stability domain of the catalytically active Ni oxide phases along with the onset electrochemical potentials of phase transitions. These metrics can guide the design of operational conditions for Ni oxide electrocatalysts. We find that the employed DFT exchange-correlation functional has the most significant influence on the computed Pourbaix diagram. Uncertainties on crystal structures, along with their related ab initio energetics, are also found to affect the size of the phase stability domain. Higher-order coupling among input parameters is found to play a crucial role in influencing the appearance and distribution of Ni phases in the diagram. Our findings suggest a need to consider variations and uncertainties associated with the computational procedures on predicted Pourbaix diagrams for materials design.

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

An Atom-Precise Approach to Damp First-Order Phase Transitions and Its Implications for Neuromorphic Signal Processing

Neuromorphic computing inspired by mammalian intelligence aims to emulate the nonlinear dynamics of biological neurons and synapses to achieve fast, low-energy, and highly efficient information processing. Brain-inspired computing relies on the design and discovery of materials exhibiting nonlinear current–voltage profiles, frequently underpinned by electronic state transitions, to achieve spiking neurons and dynamically tunable synapses. A signature challenge in the design of artificial neurons is controlling the steepness of first-order transitions in active elements, as abrupt transitions are at risk of driving unstable voltage and temperature oscillations, which result in catastrophic device failure. A critical knowledge gap is the lack of structure–function correlations mapping the composition and atomistic structure of crystalline solids to nonlinear dynamical response characteristics. Here, we address the key question of how modification of atomistic structure correlates with alteration of neuron-like functionality. Constructing oscillator circuits from millimeter-scale single crystals enables high-resolution atomic structure solutions, which we use to demonstrate that the selective positioning of Pb cations modifies charge ordering along a one-dimensional CuxV2O5 framework even at low insertion stoichiometries, thereby providing an atom-precise design parameter for damping first-order transitions. We use temperature-variant X-ray diffraction and X-ray spectroscopy to elucidate the suppression of Cu-ion shuttling based on the precise positioning of Pb ions in seven-coordinated tunnel interstitial sites as the mechanistic basis for transition broadening, thus bridging a critical gap between statistical mechanics and quantum chemical descriptions of phase transitions. Such mechanistic understanding thus paves the way to site-selective modification strategies for modulating the sharpness of first-order transitions, with an exemplary demonstration here in tuning neuronal signal processing.

Crystal structure↗
DOE OSTI2020

On a variational formulation of the weakly nonlinear magnetic Rayleigh–Taylor instability

The magnetic-Rayleigh–Taylor (MRT) instability is a ubiquitous phenomenon that occurs in magnetically-driven Z-pinch implosions. It is important to understand this instability since it can decrease the performance of such implosions. In this work, I present a theoretical model for the weakly nonlinear MRT instability. I obtain such a model by asymptotically expanding an action principle, whose Lagrangian leads to the fully nonlinear MRT equations. After introducing a suitable choice of coordinates, I show that the theory can be cast as a Hamiltonian system, whose Hamiltonian is calculated up to the sixth order in a perturbation parameter. The resulting theory captures the harmonic generation of MRT modes. It is shown that the amplitude at which the linear magnetic-Rayleigh–Taylor instability exponential growth saturates depends on the stabilization effect of the magnetic-field tension. Overall, the theory provides an intuitive interpretation of the weakly nonlinear MRT instability and provides a systematic approach for studying this instability in more complex settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics

In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.

Chemistry↗
DOE OSTI2020

Cosmological angular trispectra and non-Gaussian covariance

Angular cosmological correlators are infamously difficult to compute due to the highly oscillatory nature of the projection integrals. Motivated by recent development on analytic approaches to cosmological perturbation theory, in this paper we present an efficient method for computing cosmological four-point correlations in angular space, generalizing previous works on lower-point functions. This builds on the FFTLog algorithm that approximates the matter power spectrum as a sum over power-law functions, which makes certain momentum integrals analytically solvable. The computational complexity is drastically reduced for correlators in a "separable" form—we define a suitable notion of separability for cosmological trispectra, and derive formulas for angular correlators of different separability classes. As an application of our formalism, we compute the angular galaxy trispectrum at tree level, with and without primordial non-Gaussianity. This includes effects of redshift space distortion and bias parameters up to cubic order. We also compute the non-Gaussian covariance of the angular matter power spectrum due to the connected four-point function, beyond the Limber approximation. Finally, we demonstrate that, in contrast to the standard lore, the Limber approximation can fail for the non-Gaussian covariance computation even for large multipoles.

79 ASTRONOMY AND ASTROPHYSICS↗