SEARCH · Search NASA

Results for “first order”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Revisiting separation: Algorithms and complexity

Abstract Linear temporal logic (LTL) with Since and Until modalities is expressively equivalent, over the class of complete linear orders, to a fragment of first-order logic known as FOMLO (first-order monadic logic of order). It turns out that LTL, under some basic assumptions, is expressively complete if and only if it has the property, called separation, that every formula is equivalent to a Boolean combination of formulas that each refer only to the past, present or future. Herein we present simple algorithms and their implementations to perform separation of the LTL with Since and Until, over discrete and complete linear orders, and translation from FOMLO formulas into equivalent temporal logic formulas. We additionally show that the separation of a certain fragment of LTL results in at most a double exponential size growth.

Oliveira, Daniel↗

DOE OSTI2025

FedOSAA: Improving Federated Learning with One-Step Anderson Acceleration

Federated learning (FL) is a distributed machine learning approach that enables multiple local clients and a central server to collaboratively train a model while keeping the data on their own devices. First-order methods, particularly those incorporating variance reduction techniques, are the most widely used FL algorithms due to their simple implementation and stable performance. However, these methods tend to be slow and require a large number of communication rounds to reach the global minimizer. We propose FedOSAA, a novel approach that preserves the simplicity of first-order methods while achieving the rapid convergence typically associated with second-order methods. Our approach applies one Anderson acceleration (AA) step following classical local updates based on first-order methods with variance reduction, such as FedSVRG and SCAFFOLD, during local training. This AA step is able to leverage curvature information from the history points and gives a new update that approximates the Newton-GMRES direction, thereby significantly improving the convergence. We establish a local linear convergence rate to the global minimizer of FedOSAA for smooth and strongly convex loss functions. Numerical comparisons show that FedOSAA substantially improves the communication and computation efficiency of the original first-order methods, achieving performance comparable to second-order methods like GIANT.

Feng, Xue [University of California, Davis]↗

DOE OSTI2023

Peculiarities of beta functions in sigma models

In this paper we consider perturbation theory in generic two-dimensional sigma models in the so-called first-order formalism, using the coordinate regularization approach. Our goal is to analyze the first-order formalism in application to β functions and compare its results with the standard geometric calculations. Already in the second loop, we observe deviations from the geometric results that cannot be explained by the regularization/renormalization scheme choices. Moreover, in certain cases the first-order calculations produce results that are not symmetric under the classical diffeomorphisms of the target space. Although we could not present the full solution to this remarkable phenomenon, we found some indirect arguments indicating that an anomaly similar to that established in supersymmetric Yang-Mills theory might manifest itself starting from the second loop. We discuss why the difference between two answers might be an infrared effect, similar to that in β functions in supersymmetric Yang-Mills theories.

Sigma Models↗

DOE OSTI2025

Addendum to SAND2023-09604 Xyce lumped-element transmission line model verification to support Empire-Cable cable SGEMP analyses

This report supplements the Verification of Empire-Cable SAND report by expanding on the use of Xyce to simulate the coupling to a transmission line cable model. While Empire-Cable solves its governing equations on a high-order, finite-element mesh with an an implicit-in-time formulation, Xyce must use a first order graph for the circuit and explicit-in-time approach to be compatible with non-linear electrical device models. Thus, given the different solution methodologies in Xyce as compared to Empire-Cable, the convergence rates are expected to be different but the overall quality of the solution should be the same. The original four canonical problems studied in the Empire-Cable verification report are replicated here running in Xyce using transmission line modeling parameters from the verification report. Overall, agreement between the codes is excellent with Xyce’s convergence rates being limited mostly to first order due to the circuit network approximation of a transmission line being a first order approximation.

42 ENGINEERING↗

DOE OSTI2024

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: Gray diffusion

Radiation hydrodynamics are a challenging multiscale and multiphysics set of equations. To capture the relevant physics of interest, one typically must time step on the hydrodynamics timescale, making explicit integration the obvious choice. On the other hand, the coupled radiation equations have a scaling such that implicit integration is effectively necessary in non-relativistic regimes. A first-order Lie-Trotter-like operator split is the most common time integration scheme used in practice, alternating between an explicit hydrodynamics step and an implicit radiation solve and energy deposition step. However, such a scheme is limited to first-order accuracy, and nonlinear coupling between the radiation and hydrodynamics equations makes a more general additive partitioning of the equations non-trivial. Here, we develop a new formulation and partitioning of radiation hydrodynamics with gray diffusion that allows us to apply (linearly) implicit-explicit Runge-Kutta time integration schemes. In conclusion, we prove conservation of total energy in the new framework, and demonstrate 2nd-order convergence in time on multiple radiative shock problems, achieving error 3–5 orders of magnitude smaller than the first-order Lie-Trotter operator split at the hydrodynamic CFL, even when Lie-Trotter applies a 3rd-order TVD Runge-Kutta scheme to the hydrodynamics equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Local Truncation Error-Informed Code Verification

The method of manufactured solutions (MMS) has become increasingly popular in conducting code verification studies on predictive codes, such as nuclear power system codes and computational fluid dynamic codes. The reason for the popularity of this approach is that it can be used when an analytical solution is not available. Using MMS, code developers are able to verify that their code is free of coding errors that impact the observed order of accuracy. While MMS is still an excellent tool for code verification, it does not identify coding errors that are of the same order as the numerical method. This paper presents a method that combines MMS with modified equation analysis (MEA), which calculates the local truncation error (LTE) to identify coding error up to and including the order of the numerical method. This method is referred to as modified equation analysis methd of manufactured solutions (MEAMMS). MEAMMS is then applied to a custom-built code, which solves the shallow water equations, to test the performance of the code verification method. MEAMMS is able to detect all coding errors that impact the implementation of the numerical scheme. To show how MEAMMS is different than MMS, they are both applied to the same first-order numerical method test problem with a first-order coding error. When there are first-order coding errors, only MEAMMS is able to identify them. Finally, this shows that MEAMMS is able to identify a larger set of coding errors while still being able to identify the coding errors MMS is able to identify.

42 ENGINEERING↗

DOE OSTI2020

Slow manifold reduction for plasma science

The classical Chapman–Enskog procedure admits a substantial geometrical generalization known as slow manifold reduction. This generalization provides a paradigm for deriving and understanding most reduced models in plasma physics that are based on controlled approximations applied to problems with multiple timescales. In this Review we develop the theory of slow manifold reduction with a plasma physics audience in mind. In particular we illustrate (a) how the slow manifold concept may be used to understand breakdown of a reduced model over sufficiently-long time intervals, and (b) how a discrete-time analogue of slow manifold theory provides a useful framework for developing implicit integrators for temporally-stiff plasma models. For readers with more advanced mathematical training we also use slow manifold reduction to explain the phenomenon of inheritance of Hamiltonian structure in dissipation-free reduced plasma models. Various facets of the theory are illustrated in the context of the Abraham–Lorentz model of a single charged particle experiencing its own radiation drag. As a culminating example we derive the slow manifold underlying kinetic quasineutral plasma dynamics up to first-order in perturbation theory. This first-order result incorporates several physical effects associated with small deviations from exact charge neutrality that lead to slow drift away from predictions based on the leading-order approximation n e = Z i n i .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Near-room-temperature ferromagnetic ordering in the pressure-induced collapsed-tetragonal phase in SrC⁢o 2 ⁢P 2

We present high-pressure electrical transport, magnetization, and single-crystal x-ray diffraction data on SrC⁢o 2 ⁢P 2 single crystals, revealing a transition from the uncollapsed tetragonal (ucT) to the collapsed tetragonal (cT) phase >∼10 GPa. This structural transition is accompanied by ferromagnetic (FM) ordering with a Curie temperature up to 260 K. First-principles-based density functional theory calculations show that there is a first-order transition between ucT and cT phases, with an onset ∼10 GPa as well as the appearance of FM ordering in the cT phase. Above ∼30 GPa, the experimental signatures of the magnetic ordering vanish in a first-order-like manner, consistent with the theoretical calculation results, indicating that SrC⁢o 2 ⁢P 2 is another example of the avoidance of quantum criticality in FM intermetallic compounds with almost room-temperature critical temperature. SrC⁢o 2 ⁢P 2 provides clear evidence that the structural, electronic, and magnetic properties associated with the cT transition are strongly entangled and are not only qualitatively captured by our first-principles-based calculations but are quantitatively reproduced as well.

density functional theory↗

DOE OSTI2021

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

QCD-collapsed domain walls: QCD phase transition and gravitational wave spectroscopy

Abstract For a discrete symmetry that is anomalous under QCD, the domain walls produced in the early universe from its spontaneous breaking can naturally annihilate due to QCD instanton effects. The gravitational waves generated from wall annihilation have their amplitude and frequency determined by both the discrete symmetry breaking scale and the QCD scale. The evidence of stochastic gravitational waves at nanohertz observed by pulsar timing array experiments suggests that the discrete-symmetry-breaking scale is around 100 TeV, assuming the domain-wall explanation. The annihilation temperature is about 100 MeV, which could naturally be below the QCD phase transition temperature. We point out that the QCD phase transition within some domains with an effective large QCDθangle could be a first-order one. To derive the phase diagram inθand temperature, we adopt a phenomenological linear sigma model with three quark flavors. The domain-wall explanation for the NANOGrav, EPTA, PPTA and CPTA results hints at a first-order QCD phase transition, which predicts additional gravitational waves at higher frequencies. If the initial formation of domain walls is also a first-order process, this class of domain-wall models predicts an interesting gravitational wave spectroscopy with frequencies spanning more than ten orders of magnitude, from nanohertz to 100 Hz.

Physics↗

DOE OSTI2023

Metal–insulator transition in composition-tuned nickel oxide films

Abstract Thin films of the solid solution Nd 1 − x La x NiO 3 are grown in order to study the expected 0 K phase transitions at a specific composition. We experimentally map out the structural, electronic and magnetic properties as a function of x and a discontinuous, possibly first order, insulator–metal transition is observed at low temperature when x = 0.2. Raman spectroscopy and scanning transmission electron microscopy show that this is not associated with a correspondingly discontinuous global structural change. On the other hand, results from density functional theory (DFT) and combined DFT and dynamical mean field theory calculations produce a 0 K first order transition at around this composition. We further estimate the temperature-dependence of the transition from thermodynamic considerations and find that a discontinuous insulator–metal transition can be reproduced theoretically and implies a narrow insulator–metal phase coexistence with x . Finally, muon spin rotation ( µ SR) measurements suggest that there are non-static magnetic moments in the system that may be understood in the context of the first order nature of the 0 K transition and its associated phase coexistence regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Mechanical analysis of heterogeneous materials with higher-order parameters

Abstract Even though heterogeneous porous materials are widely used in a variety of engineering and scientific fields, such as aerospace, energy-storage technology, and bio-engineering, the relationship between effective material properties of porous materials and their underlying morphology is still not fully understood. To contribute to this knowledge gap, this paper adopts a higher-order asymptotic homogenization method to numerically investigate the effect of complex micropore morphology on the effective mechanical properties of a porous system. Specifically, we use the second-order scheme that is an extension of the first-order computational homogenization framework, where a generalized continuum enables us to introduce length scale into the material constitutive law and capture both pore size and pore distribution. Through several numerical case studies with different combinations of porosity, pore shapes, and distributions, we systematically studied the relationship between the underlying morphology and effective mechanical properties. The results highlight the necessity of higher-order homogenization in understanding the mechanical properties and reveal that higher-order parameters are required to capture the role of realistic pore morphologies on effective mechanical properties. Furthermore, for specific pore shapes, higher-order parameters exhibit dominant influence over the first-order continuum.

42 ENGINEERING↗

DOE OSTI2025

Structural Distortions and Uniaxial Negative Thermal Expansion in the Polar Dion–Jacobson Oxide RbNdTa 2 O 7

We provide deeper insight into the crystal structures, sequential structural phase transitions (I2cm → Cmce → I4/mcm → P4/mmm), thermal expansion, and electronic properties of the n = 2 Dion–Jacobson polar oxide RbNdTa 2 O 7 , through X-ray powder diffraction, neutron powder diffraction, Raman studies, and density functional theory calculations. We observed a uniaxial negative thermal expansion (NTE) across the first-order transition, I2cm → Cmce, where the unit cell contracts along the c-axis, which is driven by a contraction of the NdTa 2 O 6 layer. Here, this NTE occurs within the temperature range of the first-order phase transition and contrasts with the corkscrew mechanism typically observed in Ruddlesden–Popper phases. In RbNdTa 2 O 7 , the I2cm (hybrid improper ferroelectric) → Cmce (antipolar) transition involves crucial changes in the bond lengths of Nd and Ta polyhedra, coupled with polar to antipolar displacement of the Nd ions, leading to a net contraction in the NdTa 2 O 6 layer along the c-axis, while preserving the overall octahedral tilting magnitude. This transition highlights the intricate interplay between the Nd and Ta coordination and the associated TaO 6 distortions. Temperature-dependent Raman spectra analysis further confirms the first-order structural transition and associated NTE, providing evidence for increased bond stiffness across this transition. Additionally, using neutron powder diffraction, we have determined that the transition I4/mcm → P4/mmm occurs at approximately 1150 K. Finally, we have calculated from DFT + U, the partial density of states, the energy bandgaps, and effective masses of the charge carriers of the polar ground structure.

Chemical structure↗

DOE OSTI2024

LaFeSi–LaFe 13-x Si x composites: Modulating magnetic and magnetocaloric properties through inherent stress manipulation

We examine structural and magnetic properties of a series of La–Fe–Si alloys in the region of concentrations where they naturally form two-phase LaFeSi–LaFe 13-x Si x composites with variable content and connectivity of LaFe 13-x Si x grains distributed within the LaFeSi matrix. Theoretical calculations confirm that the LaFeSi constituent is magnetically and structurally inert below room temperature and at pressures between -10 and 10 GPa. The LaFe 13-x Si x constituent, on the other hand, is magnetically and structurally active: it exhibits first-order magnetostructural transformations that, in addition to x Si , can be controlled with temperature, magnetic field, and pressure. In composites where the concentration of the inactive constituent is ~70 wt. % or greater, the standard, single-step, LaFe 13-x Si x first-order phase transformation proceeds in two steps separated by over 30 K in a zero magnetic field. Increasing the magnetic field recouples the two steps and restores the single-step phase transformation pathway. We analyze the roles of stresses caused by both thermal expansion mismatch and the first-order magnetic phase transition in LaFe 13-x Si x to rationalize the observed physical behaviors that emerge as the temperature or/and magnetic field vary.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Time-reversed Young’s experiment: Deterministic, diffractionless second-order interference effect

The classic Young’s double-slit experiment exhibits first-order interference, producing alternating bright and dark fringes shaped by the diffraction effect of the slits. In contrast, here we demonstrate that its time-reversed configuration generates an ideal, deterministic second-order ‘ghost’ interference pattern free from diffraction and first-order effects with the use of only a position-fixed detector. The pattern’s size is governed by the dimensions of the ‘effectively extended light source’ formed by point emitter(s). Beyond highlighting the nonreciprocal physics between the two configurations, this system unlocks a range of novel phenomena inaccessible in traditional double-slit experiments. These include fully programmable, digitized interference fringe formations and the ability to align the pattern plane with the source plane on the same side of the setup. Remarkably, our proposed experiment achieves these outcomes without relying on nonclassical correlations or quantum entanglement. Here, by restoring time-reversal symmetry and eliminating diffraction, our theoretical analysis reveals that this approach offers exciting potential for advancing optical imaging and sensing technologies with improved transverse and longitudinal phase-shift sensitivity and resolution limit beyond current limitations.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2020

Nonequilibrium steady state and heat transport in nonlinear open quantum systems: Stochastic influence action and functional perturbative analysis

In this paper, we show that a nonequilibrium steady state (NESS) exists at late times in open quantum systems with weak nonlinearity by following its nonequilibrium dynamics with a perturbative analysis. We consider an oscillator chain containing three-types of anharmonicity: cubic α- and quartic β-type Fermi–Pasta–Ulam–Tsingou (FPUT) nearest-oscillator interactions and the on-site (pinned) Klein–Gordon (KG) quartic self-interaction. Assuming weak nonlinearity, we introduce a stochastic influence action approach to the problem and obtain the energy flows in different junctures across the chain. The formal results obtained here can be used for quantum transport problems in weakly nonlinear quantum systems. For α-type anharmonicity, we observe that the first-order corrections do not play any role in the thermal transport in the NESS of the configuration we considered. For KG and β-types anharmonicity, we work out explicitly the case of two weakly nonlinearly coupled oscillators, with results scalable to any number of oscillators. We examine the late-time energy flows from one thermal bath to the other via the coupled oscillators, and show that both the zeroth- and the first-order contributions of the energy flows become constant in time at late times, signaling the existence of a late-time NESS to first order in nonlinearity. Our perturbative calculations provide a measure of the strength of nonlinearity for nonlinear open quantum systems, which may help control the mesoscopic heat transport distinct from or close to linear transport. Furthermore, our results also give a benchmark for the numerical challenge of simulating heat transport. Our setup and predictions can be implemented and verified by investigating heat flow in an array of Josephson junctions in the limit of large Josephson energy with the platform of circuit QED.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Phenomenological theory of the magnetic 90° helical state

We explore a phenomenological phase diagram for the magnetic helical state with 90° turn angle between neighboring spins in the external magnetic field. Such a state is formed by the Eu spin layers in the superconducting iron arsenide RbEuFe 4 As 4 . The peculiarity of this spin configuration is that it is not realized in the standard Heisenberg model with bilinear exchange interactions. A minimum model allowing for such a state requires the biquadratic nearest-neighbor interaction term. In addition, in tetragonal materials, the 90° helix state may be stabilized by the in-plane fourfold anisotropy term, which also fixes helix orientation with respect to the crystal lattice. Such a system has a very rich behavior in the external magnetic field. The magnetic field induces the metamagnetic transition to the double-periodic state with the moment angles (α α -α,-α) with respect to the field for the four subsequent spins. The transition field to this state from the deformed helix is determined by the strength of biquadratic interaction. The transition is second-order for small biquadratic coupling and becomes first-order when this coupling exceeds the critical value. On the other hand, the aligned state at high magnetic field becomes unstable with respect to the formation of an incommensurate fan state, which transforms into the double-periodic state with decreasing magnetic field. The range of this incommensurate state near the saturation field is proportional to square of the biquadratic coupling. In addition, when the magnetic field is applied along one of four the equilibrium moment directions, the deformed helix state experience the first-order rotation transition at the field determined by the fourfold anisotropy.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Catalytic Site Requirements for N2O Decomposition on Cu-, Co-, and Fe-SSZ-13 Zeolites

N2O decomposition is investigated on Cu, Co and Fe-exchanged SSZ-13 zeolite catalysts at relatively low metal loadings. The catalysts are synthesized by solution ion exchange, and subjected to X-ray diffraction (XRD), temperature-programed-reduction by H2 (H2-TPR), temperature-programed-reaction of N2O (N2O-TPR) coupled with in-situ transmission FTIR, and finally steady-state flow reaction tests. At low N2O pressures (< 0.05 kPa), all catalysts display pseudo first-order kinetics. From Arrhenius analysis, Cu and Fe-SSZ-13 display very different apparent activation energies but similar pre-exponential factors, suggesting their similar reaction mechanisms. N2O decomposition follows a dual-site mechanism, occurring on dimeric M-O-M sites in these catalysts, and O2 is formed by the combination of two O ad-atoms from two vicinal metal sites. Under low N2O pressure (0.05 kPa) and first-order kinetic regime, the reaction is limited by N-O cleavage on bare metal active sites. In comparison to Cu-SSZ-13, the much higher N2O decomposition rate over Fe-SSZ-13 is attributed to the much lower activation barriers for the N-O cleavage step. N2O decomposition occurs on isolated Co2+ ions in Co-SSZ-13. The rate-limiting step is N-O cleavage on an O-occupied Co site in the low-pressure first order kinetic regime. This single-site mechanism leads to much higher pre-exponential factors as compared to the dual-site mechanism. This beneficial factor for reaction rate enhancement, however, is compromised by the much higher activation barriers over this catalyst.

Lin, Fan↗