Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Nonlinear mode coupling and energetics of driven magnetized shear-flow turbulence

To comprehensively understand the saturation of two-dimensional (2D) magnetized Kelvin–Helmholtz-instability-driven turbulence, energy transfer analysis is extended from the traditional interaction between scales to include eigenmode interactions, by using the nonlinear couplings of linear eigenmodes of the ideal instability. While both kinetic and magnetic energies cascade to small scales, a significant fraction of turbulent energy deposited by unstable modes in the fluctuation spectrum is shown to be re-routed to the conjugate-stable modes at the instability scale. They remove energy from the forward cascade at its inception. The remaining cascading energy flux is shown to attenuate exponentially at a small scale, dictated by the large-scale stable modes. Guided by a widely used instability-saturation assumption, a general quasi-linear model of instability is tested by retaining all nonlinear interactions except those that couple to the large-scale stable modes. These complex interactions are analytically removed from the magnetohydrodynamic equations using a novel technique. Observations are an explosive large-scale vortex separation instead of the well-known merger of 2D, a dramatic enhancement in turbulence level and spectral energy fluxes, and a reduced small-scale dissipation length scale. These show the critical role of the stable modes in instability saturation. Possible reduced-order turbulence models are proposed for fusion and astrophysical plasmas, based on eigenmode-expanded energy transfer analyses.

Physics↗
DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

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

Constraints on Large-Scale White Noise in the Cosmic Density Field

We present observational constraints on large-scale white noise (LSWN) in the cosmic density field, a phenomenon predicted to arise from non-linear mode coupling during cosmological evolution. Building on the theoretical framework of Paper I, where we demonstrated that non-linearities inevitably redistribute power from small to large scales through mode mixing, we confront these predictions with current cosmological data. We modify the CLASS Boltzmann code to incorporate a white noise component $k_\mathrm{BH}/k$ in the primordial power spectrum and perform parameter estimation using current cosmological data. The non-detection of excess power on the largest observable scales places stringent upper bounds: $k_\mathrm{BH} \leq 1.80 \times 10^{-13}~\mathrm{Mpc}^{-1}$ at 99% confidence. These constraints imply the primordial power spectrum must deviate from perfect scale invariance on small scales, either through a cutoff at $k_{\mathrm{cut}} \lesssim 3~\mathrm{pc}^{-1}$ or through running of the spectral index with $α_s \lesssim -0.015$. Our results demonstrate that LSWN provides a powerful probe of the primordial spectrum at scales orders of magnitude smaller than directly observable, offering unique constraints on early-universe physics.

Barenboim, Gabriela [Valencia U., IFIC; Valencia U↗
DOE OSTI2026

Primary defect production from molecular dynamics simulations of high-energy displacement cascades in NbMoTaW alloys

In this work, we report on large-scale molecular-dynamics (MD) simulations of displacement cascades in equiatomic NbMoTaW alloys at PKA energies ranging from 0.15 to 150 keV. We find defect production to be strongly dependent on recoil energy, scaling sublinearly up to 10 keV, and linearly thereafter. We find the sublinear regime to be defined by low values of surviving Frenkel pairs, typically found as isolated point defects or small defect clusters, while at higher recoil energies dense cascades become more frequent, leading to splitting into subcascades and the production of relatively large prismatic-dislocation loops with ⟨111⟩ and ⟨001⟩ Burgers vectors. These loops immobilize large fractions of defects, leading to a rapid growth of the number of surviving defects in the linear regime. We also anneal post-cascade defect configurations using object-kinetic Monte Carlo (OKMC) simulations to account for intracascade recombination on time scales not accessible to MD simulations. Cascade annealing is strongly temperature dependent, with the OKMC simulations only showing significant recovery at 1000 K but not below. Our results are in general agreement with existing published data for refractory concentrated alloys.

Zhou, Xinran [Lawrence Berkeley National Laborator↗
DOE OSTI2024

Crustal permeability generated through microearthquakes is constrained by seismic moment

We link changes in crustal permeability to informative features of microearthquakes (MEQs) using two field hydraulic stimulation experiments where both MEQs and permeability evolution are recorded simultaneously. The Bidirectional Long Short-Term Memory (Bi-LSTM) model effectively predicts permeability evolution and ultimate permeability increase. Our findings confirm the form of key features linking the MEQs to permeability, offering mechanistically consistent interpretations of this association. Transfer learning correctly predicts permeability evolution of one experiment from a model trained on an alternate dataset and locale, which further reinforces the innate interdependency of permeability-to-seismicity. Models representing permeability evolution on reactivated fractures in both shear and tension suggest scaling relationships in which changes in permeability (Δk) are linearly related to the seismic moment (M) of individual MEQs as Δk ∝ M. This scaling relation rationalizes our observation of the permeability-to-seismicity linkage, contributes to its predictive robustness and accentuates its potential in characterizing crustal permeability evolution using MEQs.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Mixed-strategy approach to band-edge analysis and modeling in semiconductors

In semiconductor physics, the Urbach rule describes an exponential dependence of the UV-Vis absorption coefficient on photon energy at the band edge. The rule is so ubiquitously applied that a single common mechanism, variability in the local potential, is generally considered to be the cause of exponential band edges. Variability in the local potential may be caused by alloy inhomogeneity, structural defects, impurities, thermal excitation, or any other mechanism that interrupts ideal crystallinity and creates disorder. A crucial point of the Urbach rule is its development for band edges which approach linearity when plotted on a logarithmic scale. The rule was not developed in the context of nonlinear band edges with complex fine structure. Even so, this model has been applied extensively to band edges with such features. Using differential analysis and deconvolution, we develop an analytical technique capable of separating the components responsible for the fine structure observed. This approach reduces challenging signal-artifact convolution problems to a simple linear combination. Furthermore, we develop a model that can be applied to the nonlinear band-edge signals extracted. In this work, we show that in the limiting case of linearity, our model is consistent with and produces the same results as the linear Urbach model. Our model is capable of bypassing “tailing” as an indirect measure of the spread in the local potential and gives a direct measure of this quantity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Harmonizing tau positron emission tomography in Alzheimer's disease: The CenTauR scale and the joint propagation model

Abstract INTRODUCTION Tau‐positron emission tomography (PET) outcome data of patients with Alzheimer's disease (AD) cannot currently be meaningfully compared or combined when different tracers are used due to differences in tracer properties, instrumentation, and methods of analysis. METHODS Using head‐to‐head data from five cohorts with tau PET radiotracers designed to target tau deposition in AD, we tested a joint propagation model (JPM) to harmonize quantification (units termed “CenTauR” [CTR]). JPM is a statistical model that simultaneously models the relationships between head‐to‐head and anchor point data. JPM was compared to a linear regression approach analogous to the one used in the amyloid PET Centiloid scale. RESULTS A strong linear relationship was observed between CTR values across brain regions. Using the JPM approach, CTR estimates were similar to, but more accurate than, those derived using the linear regression approach. DISCUSSION Preliminary findings using the JPM support the development and adoption of a universal scale for tau‐PET quantification. Highlights Tested a novel joint propagation model (JPM) to harmonize quantification of tau PET. Units of common scale are termed “CenTauRs”. Tested a Centiloid‐like linear regression approach. Using five cohorts with head‐to‐head tau PET, JPM outperformed linearregressionbased approach. Strong linear relationship was observed between CenTauRs values across brain regions.

Neurosciences & Neurology↗
DOE OSTI2025

Non-linear enhancement of ultrafast X-ray diffraction through transient resonances

Diffraction-before-destruction imaging with ultrashort X-ray pulses can visualize non-equilibrium processes, such as chemical reactions, with sub-femtosecond precision in the native environment. Here, a nanospecimen diffracts a single X-ray flash before it disintegrates. The sample structure can be reconstructed from the coherent diffraction image (CDI). State-of-the-art X-ray snapshots lack high spatial resolution because of weak diffraction signal. Bleaching effects from photo-ionization significantly restrain image brightness scaling. We find that non-linear transient ion resonances can overcome this barrier if X-ray laser pulses are shorter than in most experiments. We compared snapshots from individual ≈ 100 nm Xe nanoparticles as a function of pulse duration and incoming X-ray fluence. Our experimental results and Monte Carlo simulations suggest that transient resonances can increase ionic scattering cross sections significantly beyond literature values. This provides a novel avenue towards substantial improvement of the spatial resolution in CDI in combination with sub-femtosecond temporal precision at the nanoscale.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

ORBIT-2: Scaling Exascale Vision Foundation Models for Weather and Climate Downscaling

Sparse observations and coarse-resolution climate models limit effective regional decision-making, underscoring the need for robust downscaling. However, existing AI methods struggle with generalization across variables and geographies and are constrained by the quadratic complexity of Vision Transformer (ViT) self-attention. We introduce ORBIT-2, a scalable foundation model for global, hyper-resolution climate downscaling. ORBIT-2 incorporates two key innovations: (1) Residual Slim ViT (Reslim), a lightweight architecture with residual learning and Bayesian regularization for efficient, robust prediction; and (2) TILES, a tile-wise sequence scaling algorithm that reduces self-attention complexity from quadratic to linear, enabling long-sequence processing and massive parallelism. ORBIT-2 scales to 10 billion parameters across 65,536 GPUs, achieving up to 4.1 ExaFLOPS sustained throughput and 74–98% strong scaling efficiency. It supports downscaling to 0.9 km global resolution and processes sequences up to 4.2 billion tokens. On 7 km resolution benchmarks, ORBIT-2 achieves high accuracy with R2 scores in range of 0.98–0.99 against observation data.

Wang, Xiao [ORNL] (ORCID:0000000165451943)↗
DOE OSTI2025

ORBIT-2 Weather and Climate Downscaling Software Repository

ORBIT-2 is a scalable foundation model for global, hyper-resolution climate and weather downscaling. ORBIT-2 incorporates two key innovations: (1) Residual Slim ViT (Reslim), a lightweight architecture with residual learning and Bayesian regularization for efficient, robust prediction; and (2) TILES, a tile-wise sequence scaling algorithm that reduces self-attention complexity from quadratic to linear, enabling long-sequence processing and massive parallelism. ORBIT-2 scales to 10 billion parameters across 65,536 GPUs, achieving up to 4.1 ExaFLOPS sustained throughput and 74–98% strong scaling efficiency. It supports downscaling to 0.9 km global resolution and processes sequences up to 4.2 billion tokens. On 7 km resolution benchmarks, ORBIT-2 achieves high accuracy with 𝑅2 scores in range of 0.98–0.99 against observation data.

Wang, Xiao [Oak Ridge National Laboratory]↗
DOE OSTI2024

Generalized Relationship Linking Water Balance and Vegetation Productivity across Site-to-Regional Scales

Evapotranspiration (ET) is a pivotal component in catchment-scale water balance and is essential for informed watershed management. Nevertheless, uncertainties in ET observation or modeling have been hindering effective water resources management. This study addresses this gap by establishing a robust, generalized linear relationship between ET and gross primary productivity (GPP) at the catchment scale. We test the linearity of the relationships between monthly GPP and ET data at 380 near-natural catchments across various climatic and landscape conditions in the contiguous U.S., yielding Pearson’s r ≥ 0.6 for 97% of the 380 catchments. We then develop a regionalization strategy to parameterize this GPP-ET relationship at the catchment scale by identifying and utilizing the linkages between the parameter values and extensively available hydroclimatic and landscape data. We demonstrate the efficacy of the proposed GPP-ET relationship and parameter regionalization strategy by their combined predictive capacity, where the predicted monthly GPP matches well with remote-sensing-based GPP product, achieving Kling-Gupta Efficient (KGE) values ≥ 0.5 for 92% of the catchments. In addition, we verify the relationship and its parameter regionalization at 35 AmeriFlux sites with KGE ≥ 0.5 for 25 sites, suggesting that the new relationship is transferable across the site, catchment, and regional scales. Furthermore, our findings are valuable for improving remote-sensing-based estimation of monthly ET and diagnosing coupled water–carbon simulations in land surface and Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Sustainable hydrogen manufacturing via renewable-integrated intensified process for refueling stations

The widescale consumer adoption of hydrogen fuel cell electric vehicles (HFCEVs) is currently hindered by the high cost of small-scale hydrogen generation and the lack of extensive hydrogen refueling infrastructure. Natural gas-based hydrogen is cheaper when produced in large volumes but is also associated with high CO 2 emissions. To counter these challenges, we propose a hybrid approach where both natural gas and renewables are integrated in a synergistic manner using a dynamic process intensification technology that can be deployed on-site for meeting local demands of refueling stations. The technology is based on sorption enhanced steam methane reforming (SE-SMR) that utilizes a combination of reaction with in-situ CO 2 adsorption for enhancing process modularity, productivity and efficiency thereby outperforming conventional SMR at small scale. We develop a mixed integer linear programming (MILP)-based optimization framework for simultaneous design and scheduling of the SE-SMR process. The simultaneous optimization provides a synergistic combination whereby the renewables allow sustainable hydrogen manufacturing and the dynamic SE-SMR allows optimal use of the intermittency of the renewables. The U.S. nationwide analysis indicates that for futuristic renewable prices and a hydrogen production capacity of 2 ton/day, hydrogen can be produced at 50% less cost compared to the current cost of small-scale hydrogen generation. Finally, the city-wise analysis with varying hydrogen demand shows that even with just 5% HFCEV market penetration level, hydrogen production cost less than $3/kg can be obtained at small scales across the United States with even cheaper hydrogen for large cities.

08 HYDROGEN↗
DOE OSTI2022

Beware of fake ν ’s: The effect of massive neutrinos on the nonlinear evolution of cosmic structure

Massive neutrinos suppress the growth of cosmic structure on small, nonlinear scales. It is thus often proposed that using statistics beyond the power spectrum can tighten constraints on the neutrino mass by extracting additional information from these nonlinear scales. We study the information content regarding neutrino mass at the field level, quantifying how much of this information arises from the difference in nonlinear evolution between a cosmology with one fluid [cold dark matter (CDM)] and two fluids ( CDM + neutrinos ). We do so by running two N -body simulations, one with and one without massive neutrinos, both with the same phases, and matching their linear power spectrum at a given low redshift. This effectively isolates the information encoded in the linear initial conditions from the nonlinear cosmic evolution. We demonstrate that, for k ≲ 1 h / Mpc , and for a single redshift, there is negligible difference in the real-space CDM field between the two simulations. This suggests that all the information regarding neutrino mass is in the linear power spectrum set by the initial conditions. Thus, any probe based on the CDM field alone will have negligible constraining power beyond that which exists at the linear level over the same range of scales. Consequently, any probe based on the halo field will contain little information beyond the linear power. We find similar results for the matter field responsible for weak lensing. We also demonstrate that there may be much information beyond the power spectrum in the 3D matter field; however, this is not observable in modern surveys via dark matter halos or weak lensing. Finally, we show that there is additional information to be found in redshift space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Scale-dependent halo bias and the squeezed limit bispectrum in the presence of radiation

Here, we investigate the gravitational effect of large-scale radiation perturbations on small-scale structure formation. In addition to making the growth of matter perturbations scale dependent, the free streaming of radiation also affects the coupling between structure formation at small and large scales. We study this using separate universe N -body simulations to compute the (isotropized) squeezed limit matter bispectrum and the linear halo bias. Our results show that the scale dependence in the growth of long-wavelength matter perturbations, caused by radiation, translates into these quantities acquiring a nontrivial scale dependence at k ≲ 0.05 Mpc -1 . In a universe with radiation composed of cosmic microwave background photons and three species of massless neutrinos, the bias of halos with b = 2 at high k will decrease by 0.29%, 0.45%, and 0.8% between k = 0.05 Mpc -1 and k = 0.0005 Mpc -1 at redshifts z = 0, 1, and 3, respectively. For objects with b >> 1 , these differences approach 0.43%, 0.68%, and 1.2%, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Gauge-invariant gravitational waves in matter beyond linearized gravity

Modeling the propagation of gravitational waves (GWs) in media other than vacuum is complicated by the gauge freedom of linearized gravity in that, once nonlinearities are taken into consideration, gauge artifacts can cause spurious acceleration of the matter. To eliminate these artifacts, we propose how to keep the theory of dispersive GWs gauge-invariant beyond the linear approximation and, in particular, obtain an unambiguous gauge-invariant expression for the energy–momentum of a GW in a dispersive medium. Using analytic tools from plasma physics, we propose an exactly gauge-invariant 'quasilinear' theory, in which GWs are governed by linear equations and also affect the background metric on scales large compared to their wavelength. As a corollary, the gauge-invariant geometrical optics of linear dispersive GWs in a general background is formulated. As an example, we show how the well-known properties of vacuum GWs are naturally and concisely yielded by our theory in a manifestly gauge-invariant form. We also show how the gauge invariance can be maintained within a given accuracy to an arbitrary order in the GW amplitude. These results are intended to form a physically meaningful framework for studying dispersive GWs in matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Application of linear prolongation to coarse mesh finite difference acceleration in CASMO5

The Coarse-Mesh Finite Difference (CMFD) method has been used for over a decade to accelerate the convergence of the Method of Characteristics (MOC) solution to the two- dimensional particle transport equation in CASMO5. Numerical testing, along with widespread use in production-level calculations, have shown that the current CMFD implementation provides stability and robustness for a wide range of realistic reactor physics problems. However, the recent development of linear prolongation has attracted attention from the community as a way to further improve the performance and stability of CMFD. Two interpolation methods for linear prolongation are presented in this work and implemented into a test version of CASMO5. The performance of the proposed interpolations, referred to as the bilinear and linear directional schemes, is evaluated in terms of runtime relative to the default constant or uniform scaling update. Numerical results indicate that the use of linear prolongation can reduce the transport solver runtime on average by approximately 10% when tested with two hundred randomly selected cases. The new directional linear interpolation, combined with default constant boundary updates, is found to provide the highest reduction in runtime for the cases analyzed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

A Tensor Network-Based Quantum Algorithm for the Nonlinear 1D Burgers' Equation

In this work, we implement a tensor network-based quantum algorithm to solve unsteady, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the compressible 1-dimensional (1D) Burgers' equation as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts to solve nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. Our framework is based on matrix product states (MPSs) and matrix product operators (MPOs). For example, the velocity field is represented by MPS, whereas the linear and nonlinear spatial differential terms of the velocity field are processed by MPOs. Our primary focus herein is to verify and validate the various tensor network components of the algorithm using solutions obtained by the classical algorithms on high performance computing (HPC) architectures. We use a classical time marching method to demonstrate the functionality of the tensor network operations to model the PDE and their robustness with the time evolution of the system. Our classical simulation results demonstrate the utility of tensor network-based operations in modeling nonlinear PDEs and highlight the necessity as well as potential advantages of using quantum simulations for these techniques.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗