SEARCH · Engineering Papers

Results for “Second-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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Study of fully coupled three-dimensional envelope instability using automatic differentiation

Automatic differentiation is a powerful tool for computing derivatives of simulation results with respect to given parameters. In this Letter, we have applied this tool to investigate the instability of a dynamical system governed by 21 ordinary differential equations. This second-order instability (named envelope instability) is driven by space-charge effects and has a significant impact on the operational regimes of particle accelerators. Our study delves into the three-dimensional envelope instability, incorporating both transverse and longitudinal coupling. Conventionally, analyzing this complex system would necessitate solving 441 ordinary differential equations, which is computationally intractable. However, by employing automatic differentiation, we were able to track only 21 equations. This approach allowed us to uncover an additional instability stopband, which arises from space-charge-induced coupling and has not been reported in previous studies. This research highlights the significant advantages of automatic differentiation in analyzing complicated dynamical systems involving a large number of ordinary differential equations.

Qiang, Ji [Lawrence Berkeley National Laboratory (↗

DOE OSTI2026

Higher-order Zeno sequences

The quantum Zeno effect typically refers to freezing the dynamics of a quantum system through frequent observations. In general, quantum Zeno dynamics is obtained with an error of order 𝒪⁢(1/𝑁), where 𝑁 is the number of projective measurements performed within a fixed evolution time. In this work, we develop higher-order Zeno sequences that achieve faster convergence to Zeno dynamics, yielding an improved error scaling of 𝒪⁢(1/𝑁 2⁢𝑘 ), where 𝑘 describes the order of the Zeno sequence. This is achieved by relating higher-order Zeno sequences to higher-order Trotter formulas that achieve similar convergence behavior. We leverage this relation to develop higher-order Zeno sequences for different manifestations of the quantum Zeno effect, including frequent projective measurements and unitary kicks. We go on to discuss achieving quantum Zeno dynamics through periodic control fields of high frequency. We explicitly develop control fields that yield a second-order type improvement in the Zeno error scaling and present shorter Zeno sequences. Finally, we discuss the connection to randomized and Uhrig dynamical decoupling to develop more efficient implementations in the weak-coupling regime.

Quantum Zeno dynamics↗

DOE OSTI2026

One-dimensional neutron diffraction from layered graphite: Reciprocal space structure and grating behavior

In this work we report observations of one-dimensional neutron diffraction from highly oriented pyrolytic graphite (HOPG), where the scattering angle varies continuously with incident angle following classical grating-like behavior. The 2D polycrystalline structure of HOPG—with highly aligned layers along the 𝑐 axis but random in-plane rotations—creates planes of scattering intensity in reciprocal space at 𝑄 𝑐 =𝑛⁢(2⁢𝜋/𝑑) where 𝑑=3.35Å is the interlayer spacing. As the Ewald sphere sweeps through reciprocal space during sample rotation, it continuously intersects these planes, producing the observed angular dispersion. We observe both first-order (𝑛=1) and second-order (𝑛=2) diffraction at conventional scattering angles (25°–70°), with peak positions that remain temperature-independent between 10 K and 294 K and follow quantitative agreement with momentum conservation 𝑄 𝑐 =𝑘⁢[sin⁡𝜓−sin⁡𝜓 𝑓 ]=𝑛⁢(2⁢𝜋/𝑑). X-ray diffraction under similar conditions shows no comparable behavior, confirming that sharp nuclear-vacuum contrast is essential. While diffraction intensities are weak (∼10 −6 of Bragg peaks), the observations demonstrate how the interplay of atomic-scale periodicity, nuclear contrast, and structural disorder enables observation of continuous diffraction curves at thermal neutron wavelengths, illustrating how HOPG's unique microstructure determines its scattering properties beyond conventional Bragg diffraction.

2-dimensional systems↗

DOE OSTI2025

Dissipative Phase Transition in the Two-Photon Dicke Model

We explore the dissipative phase transition of the two-photon Dicke model, a topic that has garnered significant attention recently. Our analysis reveals that while single-photon loss does not stabilize the intrinsic instability in the model, the inclusion of two-photon loss restores stability, leading to the emergence of superradiant states, which coexist with the normal vacuum states. Using a second-order cumulant expansion for the photons, we derive an analytical description of the system in the thermodynamic limit, which agrees well with the exact calculation results. Additionally, we present the Wigner function for the system, shedding light on the breaking of the 𝑍4 symmetry inherent in the model. These findings offer valuable insights into stabilization mechanisms in open quantum systems and pave the way for exploring complex nonlinear dynamics in two-photon Dicke models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Tuning of altermagnetism by strain

For all collinear altermagnets, we sort out piezomagnetic free-energy invariants allowed in the nonrelativistic limit and relativistic piezomagnetic invariants bilinear in the Néel vector $\mathbf{L}$ and magnetization $\mathbf{M}$, which include strain-induced Dzyaloshinskii-Moriya interaction. The symmetry-allowed responses are fully determined by the nonrelativistic spin Laue group. In the nonrelativistic limit, two distinct mechanisms are discussed: the band-filling mechanism, which exists in metals and is illustrated using the simple two-dimensional Lieb lattice model, and the temperature-dependent exchange-driven mechanism, which is illustrated using first-principles calculations for transition-metal fluorides. The leading second-order nonrelativistic term in the strain-induced magnetization is also obtained for CrSb. Piezomagnetism due to the strain-induced Dzyaloshinskii-Moriya interaction is calculated from first principles for transition-metal fluorides, MnTe, and CrSb. Finally, we discuss triplet superconducting correlations supported by altermagnets and protected by inversion rather than time-reversal symmetry. We apply the nonrelativistic classification of Cooper pairs to describe the interplay between strain and superconductivity in the two-dimensional Lieb lattice and in bulk rutile structures. Here, we show that triplet superconductivity is, on average, unitary in an unstrained altermagnet, but becomes non-unitary under piezomagnetically active strain.

FOS: Physical sciences↗

DOE OSTI2026

High-order cumulants and correlation functions near the critical point from molecular dynamics

We present a systematic investigation of particle-number fluctuations in the crossover region near the critical end point of a first-order phase transition using molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our prior studies to third- and fourth-order cumulants in both coordinate- and momentum-space acceptances and integrated correlation functions (factorial cumulants). We find that, even near the critical point, non-Gaussian cumulants equilibrate on timescales comparable to those of the second-order cumulants, but show stronger finite-size effects. The presence of interactions and of the critical point leads to strong deviations of the cumulants from the ideal-gas baseline in coordinate space; these deviations are expected to persist in momentum space in the presence of collective expansion. In particular, the kurtosis becomes strongly negative, κσ 2 ≪ − 1 , on the crossover side of the critical point. However, this signal is significantly diluted once an efficiency cut used to distinguish protons from baryons is applied, leading to |κσ 2 | ≲ 1 even in the presence of the critical point. We discuss our results in the context of ongoing measurements of proton-number cumulants in heavy-ion collisions in RHIC-BES-II.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Correlation-aware binning for small-angle neutron scattering via Gaussian-process inference

Binning in small-angle neutron scattering (SANS) is typically performed empirically, with fixed parameters chosen for convenience rather than statistical optimality. Such practices often fail to balance statistical precision and spatial resolution, leading to inconsistencies across instruments and datasets. Here we establish a correlation-aware framework that determines the optimal bin width from first principles by extending the classical Freedman–Diaconis (FD) rule to account for inter-bin correlations with a Gaussian process. In this formulation, the scattering intensity is treated as a smooth stochastic field whose statistical coherence is described by a covariance matrix. Analytical expressions of errors derived from this model yield closed-form criteria that separate the total deviation into contributions from counting noise, aliasing distortion and curvature-dependent correlation effects. Expressed in reduced variables, the resulting dimensionless error surface reveals a continuous transition from the uncorrelated FD regime to the correlation-dominated limit, providing a unified description of noise suppression and resolution control. Because the formulation depends only on the profile characteristics of scattering intensity I(Q), specifically its average intensity and first- and second-order derivatives, it applies generally to any SANS measurement regardless of sample, instrument or geometry. Experimental validation using small- and ultra-small-angle neutron scattering data confirms the predicted scaling behavior, demonstrating that correlation-aware inference systematically reduces mean-squared error and enables information-efficient reproducible data reduction across materials and instruments.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗

DOE OSTI2025

Temperature and Strain Sensing Characteristics of a 128° YX-Cut LiNbO 3 Rayleigh-Mode SAW Sensor From Room to Cryogenic Temperatures

Accurate, passive, and wireless monitoring of cryogenic hardware is essential for high-energy physics, space propulsion, and biomedical instrumentation. Here, this study quantifies the coupled temperature-strain behavior of Rayleigh-mode surface acoustic-wave (SAW) delay-line sensors fabricated on 128° YX-cut LiNbO 3 . A nonlinear finite element (FE) model incorporating Varshni-based elastic constants, higher-order thermal expansion, and temperature-dependent piezo- and dielectric coefficients was developed and validated experimentally between 280K and 80K. Free-standing (first test condition) and bonded/wired (second test condition) devices exhibited indistinguishable thermal responses; the average temperature coefficient of delay (TCD) in the critical cryogenic range from 130K down to 80K differed by only 0.15 ppm/K (0.32%), confirming that bonding-induced stress is negligible. Over 280-80K the measured TCD was 61.77 ppm/K, while the FE model predicted an equivalent temperature coefficient of frequency (TCF) of −62.74 ppm/K with an overall coefficient of determination R2 = 0.998. In the critical cryogenic interval 130-80K the TCD fell to 47.66 ppm/K, indicating improved thermal stability at low temperature. Controlled loading (0-300 με) revealed a strain coefficient of delay (SCD) that rises from 0.53 ± 0.02 ppm/με at 300K to 1.05 ± 0.02 ppm/με at 80K. This modest sensitivity confirms that, for temperature sensing, strain is a second-order perturbation above 135K but must be compensated at deeper cryogenic levels. Overall, this work establishes a predictive multiphysics model together with repeatable wired measurements that confirm the suitability of SAW sensors for temperature and strain monitoring in extreme cryogenic environments, while also providing a baseline for future wireless implementations.

25 ENERGY STORAGE↗

DOE OSTI2026

Deep learning-assisted modeling for χ (2) nonlinear optics

Modeling second-order (χ(2)) nonlinear optical processes remains computationally expensive due to the need to resolve fast field oscillations and simulate wave propagation using methods such as the split-step Fourier method (SSFM). This can become a bottleneck in real-time applications, such as high-repetition-rate laser systems requiring rapid feedback and control. We present a long short-term memory-based surrogate model trained on SSFM simulations generated from a start-to-end model of the photocathode drive laser at SLAC National Accelerator Laboratory’s Linac Coherent Light Source II. The model achieves over 250× speedup while maintaining high fidelity, enabling future real-time optimization and laying the foundation for data-integrated modeling frameworks and digital twins of laser systems.

Accelerator Physics (physics.acc-ph)↗

DOE OSTI2026

Explicit Monotone Stable Super-Time-stepping Methods for Finite Time Singularities

We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$⁡(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

A Unidirectional Two-Compartment Neuron Circuit with On-chip STDP learning

Most neuromorphic chips implement the single-compartment point neuron model where synapse circuits connect directly to a leaky integrate and fire (LIF) soma circuit. However, when using a biologically plausible soma circuit (e.g., Hodgkin-Huxley neuron model), an interface circuitry, such as a current conveyor circuit, is needed to transmit synaptic current to the soma circuit. This is especially true for ultra-low power neuron circuits, where membrane capacitance is on the order of 20 fF. This need for an interface circuit arises because the parasitic capacitance and leakage current caused by fabrication mismatch and second-order effects of the output transistors in the synapse circuits can disturb the spiking dynamics of the soma circuit if connected without an interface. Using an interface circuit to isolate the soma’s membrane capacitor from synapses resolves this issue. We propose to use a unidirectional resistor (a transconductance circuit) to connect the synapse and soma circuits instead of conventional current conveyor circuits. Using a biologically plausible spike pattern detection model, we show that the on-chip spike-timing-dependent plasticity (STDP) learning performance of the proposed unidirectional two-compartment neuron circuit is similar to a single-compartment circuit (with a current conveyor as an interface) and additionally, it is more power-efficient and biologically plausible. The chip is fabricated in the Taiwan Semiconductor Manufacturing Company (TSMC) 250 nm technology node and comprises a single neuron circuit.

Gautam, Ashish [ORNL]↗

DOE OSTI2026

Did You Win the GPU Cloud Lottery? Benchmarking from TFLOPS to Tokens/$

Cloud GPUs are commonly assumed to deliver consistent performance for a given GPU model. This assumption does not always hold: cloud providers employ diverse system configurations and virtualization mechanisms, and GPUs themselves exhibit non-negligible manufacturing variability (the silicon lottery). In this work, we present a large-scale measurement study of GPU performance variability across 11 cloud providers, covering over 3,500 physical GPUs and 6,800 benchmark runs. Our hierarchical analysis shows that while execution-level variation stays below 9%, performance varies by up to 38% across devices and providers for the same GPU model. Regression analysis indicates that driver- and OS-related software factors contribute less than 1% of the variance; instead, silicon lottery effects dominate observed performance variation, and cloud providers further amplify them through persistent, systematic second-order effects.

Slynko, Platon [Silicon Data, New York, USA] (ORCI↗

DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗

DOE OSTI2026

Emergent pair density wave order across a Lifshitz transition

We numerically investigate the telltale signs of pair-density-wave order (PDW) in the Kondo-Heisenberg chain by focusing on the momentum resolved spectrum in different parameter regimes. Density matrix renormalization group calculations reveal that this phase is characterized by a dispersion with two minima and four Fermi points, indicating the emergence of an effective next-nearest-neighbor hopping that arises as a second-order effect to avoid magnetic frustration. The pairs appear in the spectrum as in-gap bound states with weight concentrated in the hole pockets. The low-energy physics can be understood by means of a generalized t−J model with next-nearest-neighbor hopping. Our results offer a guide for searching for experimental signatures, and for other models that can realize PDW physics.

Yang, Luhang [University of Tennessee Knoxville]↗

DOE OSTI2025

Hybrid HEFA-HDCJ Process for the Production of Jet Fuel Blendstocks

09 BIOMASS FUELS↗

DOE OSTI2025

Using quantum noise correlation analysis to measure ion temperature

Quantum noise correlation analysis was tested at the proof-of-concept level as a technique to measure ion temperature in a plasma. If eventually successful, this technique could enable a compact, inexpensive, and robust ion temperature diagnostic suitable for a burning plasma environment. Ion temperature is a key parameter determining the fusion performance of a burning plasma, as the fusion cross-section has a strong dependence on ion temperature. This ion temperature diagnostic would require only a small optical view of the plasma through a port to passively record impurity line emission. The instrumentation would be remote from the reactor behind the neutron and bio-shielding. The technique relies solely on quantum correlations of the photons emitted by a plasma impurity to measure ion temperature; there is no grating dispersion of an emission line-width or pulse-height analysis of photon energy. This measurement innovation was tested with instrumentation consisting of two single-photon detectors with high timing resolution, a time-tagging unit, and simple light collection optics. This instrumentation measures the second-order correlation between the light intensity falling on the two detectors. The next steps beyond the proof-of-concept level will be development of diagnostic designs for application of this technique to high-temperature and burning plasmas. Arrays of single-photon avalanche detectors to multiplex measurements of photon correlation will be required to reduce signal integration time to an acceptable duration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Hybrid HEFA-HDCJ Process for the Production of Jet Fuel Blendstocks

The hydrotreatment of bio-oil derived from the pyrolysis and biocrude from hydrothermal liquefaction of lignocellulosic materials to produce hydrocarbons faces significant technological challenges, mainly due to the high reactivity and poor thermal stability of bio-oil, resulting in the formation of large quantities of coke. This problem has been addressed by existing PNNL patents with a two-step hydrotreatment technology in which the bio-oil is first stabilized with a noble hydrogenation metal (often Pt or Ru). Then, in the second step, the bio-oil is deoxygenated with a Ni-Mo or Co-Mo sulfide catalyst. The main problem with this approach is that the Pt/Ru catalysts deactivate easily in the presence of S or other impurities, which are commonly present in pyrolysis oils. In this project, we explored technological solutions to mitigate coke formation, avoiding the use of Pt/Ru catalysts. Our strategy is based on three actions: (1) Bio-oil stabilization in the presence of alcohols. In this project, we studied the stabilization with butanol. (2) the use of a cosolvent to solubilize the bio-oil. Because coke formation reactions are second-order reactions a reduction in the concentration of reactive bio-oil molecules. In this case, we used yellow greases as a co-solvent. (3) Separation of bio-oil reactive fractions. In this project, we studied the removal of water-soluble fractions. Our batch co-hydrotreatment studies confirmed that the addition of butanol and methanol and the blend with lipids effectively contributed to mitigating coke formation (reducing coke yield to about 1 wt.% %). The removal of sugars did not have a noticeable effect on the overall coke yield, suggesting that coke precursors are present in all bio-oil fractions. Our analytical work suggests that they may be concentrated in the water-insoluble/CH 2 Cl 2 insoluble fractions of pyrolysis oils. Although the technological strategies tested resulted in significant coke reductions, the levels achieved were not sufficiently low to ensure a reliable operation in continuous, fixed-bed trickle-bed reactors. Long runs of more than 100 hours (maximum: 255 h) of co-processing time on stream were achieved in a continuous 40 mL reactor. When the same test was conducted in a larger 400 mL reactor, pressure drop increases associated with coke formation were observed. This increase in coke formation could be due to larger temperature gradients in the bed. Hydrodeoxygenation tests in moving bed reactors and using more active hydrogenation catalysts (for example Ni) could lead to more reliable operations. Unfortunately, our team did not have access to such experimental setups. The technoeconomic analysis suggests that although alcohol use is an effective means to reduce coke formation, the use of alcohol increases production cost. Thus, its use needs to be minimized. A delicate balance needs to be found between the use of technological solutions that allow the reliable operation of the system (stabilization with Ni catalysts, use of small quantities of solvents, processing in moving bed reactors) with a tolerable level of coke formation for the hydrodeoxygenation reactor used and that result in minimum production costs.

09 BIOMASS FUELS↗

DOE OSTI2025

Block-Structured Operator Inference for Coupled Multiphysics Model Reduction

This work presents a block-structured formulation of Operator Inference as a way to learn structured reduced-order models for multiphysics systems. The approach specifies the governing equation structure for each physics component and the structure of the coupling terms. Once the multiphysics structure is specified, the reduced-order model is learned from snapshot data following the nonintrusive Operator Inference methodology. In addition to preserving physical system structure, which in turn permits preservation of system properties such as stability and second-order structure, the block-structured approach has the advantages of reducing the overall dimensionality of the learning problem and admitting tailored regularization for each physics component. The numerical advantages of the block-structured formulation over a monolithic Operator Inference formulation are demonstrated for aeroelastic analysis, which couples aerodynamic and structural models. For the benchmark test case of the AGARD 445.6 wing, block-structured Operator Inference provides an average 20% online prediction speedup over monolithic Operator Inference across subsonic and supersonic flow conditions in both the stable and fluttering parameter regimes while preserving the accuracy achieved with monolithic Operator Inference.

42 ENGINEERING↗