Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Scaling Studies”

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

Effects of artificial collisions, filtering, and nonlocal closure approaches on Hermite-based Vlasov–Poisson simulations

Kinetic simulations of collisionless plasmas are computationally challenging due to phase-space mixing and filamentation, resulting in fine-scale velocity structures. This study compares three methods developed to reduce artifacts related to limited velocity resolution in Hermite-based Vlasov–Poisson simulations: artificial collisions, filtering, and nonlocal closure approaches. We evaluate each method's performance in approximating the linear kinetic response function and suppressing recurrence in linear and nonlinear regimes. Numerical simulations of Landau damping demonstrate that artificial collisions, particularly higher orders of the Lenard-Bernstein collisional operator, most effectively recover the correct damping rate across a range of wavenumbers. Moreover, Hou-Li filtering and nonlocal closures underdamp high wavenumber modes in linear simulations, and the Lenard-Bernstein collisional operator overdamps low wavenumber modes in both linear and nonlinear simulations. This study demonstrates that hypercollisions offer a robust approach to kinetic simulations, accurately capturing collisionless dynamics with limited velocity resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Enabling accurate chemical modeling of shocked energetic materials using a machine learning interatomic potential

Understanding the complex chemistry of organic materials under dynamic compression is important for many applications, but it is challenging due to the large number of reactions occurring at various time scales. Here, in this study, we develop a machine learning potential based on Chebyshev polynomials to study the insensitive energetic material 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under detonation. We discuss a strategy for constructing diverse training data needed to capture the complex chemistry of TATB. Our potential demonstrates strong transferability across a wide range of thermodynamic conditions and other explosives, enabling accurate and reliable chemical modeling of organic materials under extreme conditions. The efficiency of our approach allows for simulations over several nanoseconds and for large system sizes, providing detailed insights into the chemistry of shocked TATB. The model accurately reproduces experimental Hugoniot equation of state data, and our simulations reveal the rapid formation of nitrogen-rich carbon clusters following shock. The methods and datasets developed here offer a robust framework for accurate chemical modeling of other shocked organic energetic materials.

Chemistry↗
DOE OSTI2024

Machine learning without a processor: Emergent learning in a nonlinear analog network

Standard deep learning algorithms require differentiating large nonlinear networks, a process that is slow and power-hungry. Electronic contrastive local learning networks (CLLNs) offer potentially fast, efficient, and fault-tolerant hardware for analog machine learning, but existing implementations are linear, severely limiting their capabilities. These systems differ significantly from artificial neural networks as well as the brain, so the feasibility and utility of incorporating nonlinear elements have not been explored. Here, we introduce a nonlinear CLLN—an analog electronic network made of self-adjusting nonlinear resistive elements based on transistors. We demonstrate that the system learns tasks unachievable in linear systems, including XOR (exclusive or) and nonlinear regression, without a computer. We find our decentralized system reduces modes of training error in order (mean, slope, curvature), similar to spectral bias in artificial neural networks. The circuitry is robust to damage, retrainable in seconds, and performs learned tasks in microseconds while dissipating only picojoules of energy across each transistor. This suggests enormous potential for fast, low-power computing in edge systems like sensors, robotic controllers, and medical devices, as well as manufacturability at scale for performing and studying emergent learning.

Science & Technology - Other Topics↗
DOE OSTI2025

Importance of $\delta B_{\|}$ on ETG stability, turbulence, and transport in NSTX

This study employs electron-scale gyrokinetic simulations to investigate the electron temperature gradient (ETG) driven instabilities, turbulence, and transport in the pedestal region of the National Spherical Torus Experiment, comparing non-lithiated (narrow pedestal) and lithiated (wide pedestal) scenarios. Our findings reveal that, in the non-lithiated case, a branch of strongly unstable ETG modes exhibiting finite parallel magnetic field fluctuations ($\delta B_{\parallel} \neq 0$) emerges at the pedestal top and upper density pedestal region. This branch is uncovered only when $\delta B_{\parallel}$ is retained in the simulations and is associated with substantial electrostatic electron heat flux. This region of strong ETG transport corresponds to the only region in the plasma where the pressure gradient is far below the critical gradient for kinetic ballooning modes. We investigated the origin of this finite $\delta B_{\parallel}$ ETG branch by analyzing the gyrokinetic field equations. Nonlinear saturation is also analyzed and contrasted for simulations with and without $\delta B_{\parallel}$. In contrast with the nonlithiated case, ETG modes in the lithiated case produce substantial transport in the steep gradient region, but are negligible at the pedestal top.

ETG↗
DOE OSTI2024

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗
DOE OSTI2025

Physics-guided dual implicit neural representations for source separation

Significant challenges exist in efficient data analysis of most advanced experimental and observational techniques because the collected signals often include unwanted contributions, such as background and signal distortions, that can obscure the physically relevant information of interest. To address this, we have developed a self-supervised machine-learning approach for source separation using a dual implicit neural representation framework that jointly trains two neural networks: one for approximating distortions of the physical signal of interest and the other for learning the effective background contribution. Our method learns directly from the raw data by minimizing a reconstruction-based loss function without requiring labeled data or pre-defined dictionaries. We demonstrate the effectiveness of our framework by considering a challenging case study involving large-scale simulated, as well as experimental, momentum-energy-dependent inelastic neutron scattering data in a four-dimensional parameter space, characterized by heterogeneous background contributions and unknown distortions to the target signal. The method is found to successfully separate physically meaningful signals from a complex or structured background even when the signal characteristics vary across all four dimensions of the parameter space. An analytical approach that informs the choice of the regularization parameter is presented. Our method offers a versatile framework for addressing source separation problems across diverse domains, ranging from superimposed signals in astronomical measurements to structural features in biomedical image reconstructions.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Binary operations on neuromorphic hardware with application to linear algebraic operations and stochastic equations

Abstract Non-von Neumann computational hardware, based on neuron-inspired, non-linear elements connected via linear, weighted synapses—so-called neuromorphic systems—is a viable computational substrate. Since neuromorphic systems have been shown to use less power than CPUs for many applications, they are of potential use in autonomous systems such as robots, drones, and satellites, for which power resources are at a premium. The power used by neuromorphic systems is approximately proportional to the number of spiking events produced by neurons on-chip. However, typical information encoding on these chips is in the form of firing rates that unarily encode information. That is, the number of spikes generated by a neuron is meant to be proportional to an encoded value used in a computation or algorithm. Unary encoding is less efficient (produces more spikes) than binary encoding. For this reason, here we present neuromorphic computational mechanisms for implementing binary two’s complement operations. We use the mechanisms to construct a neuromorphic, binary matrix multiplication algorithm that may be used as a primitive for linear differential equation integration, deep networks, and other standard calculations. We also construct a random walk circuit and apply it in Brownian motion simulations. We study how both algorithms scale in circuit size and iteration time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Scanning Transmission Electron Microscopy–Atom Probe Tomography Correlative Analysis for the Characterization of Solute-defect Interactions

Atom probe tomography (APT) and (scanning) transmission electron microscopy ((S)TEM) are complementary techniques that provide spatially resolved chemical and structural information at the atomic scale. Here, in this study, we employ two different STEM/APT correlative analysis methods to investigate Cr segregation at dislocation loops in ultra-high purity Fe–Cr alloys. APT needles for the correlative analysis were extracted either from bulk material or from thinned TEM lamellae. STEM analysis was used to determine the Burgers vectors of ion-irradiation-induced dislocation loops, while APT reconstruction of the same region revealed the Cr segregation to these loops. We extended the g•b = 0 invisibility criterion of dislocation loops from TEM mode in a lamella to STEM mode in a needle-shaped specimen. STEM and APT analysis on the same needle provide straightforward correlative analysis, although it is limited by a small observation volume. In contrast, iterative STEM analysis of TEM lamellae, followed by the selective extraction of specific regions of interest for APT analysis, expands the observation area by up to 100 times but requires additional time-consuming steps for APT needle extraction from the lamellae.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Hemispherical variance anomaly and reionization optical depth

ABSTRACT Cosmic microwave background (CMB) full-sky temperature data show a hemispherical asymmetry in power nearly aligned with the Ecliptic, with the Northern hemisphere displaying an anomalously low variance, while the Southern hemisphere appears consistent with expectations from the best-fitting theory, Lambda Cold Dark Matter (ΛCDM). The low signal-to-noise ratio in current polarization data prevents a similar comparison. Polarization realizations constrained by temperature data show that in ΛCDM the lack of variance is not expected to be present in polarization data. Therefore, a natural way of testing whether the temperature result is a fluke is to measure the variance of CMB polarization components. In anticipation of future CMB experiments that will allow for high-precision large-scale polarization measurements, we study how the variance of polarization depends on ΛCDM-parameter uncertainties by forecasting polarization maps with Planck’s Markov chain Monte Carlo chains. We show that polarization variance is sensitive to present uncertainties in cosmological parameters, mainly due to current poor constraints on the reionization optical depth τ, which drives variance at low multipoles. We demonstrate how the improvement in the τ measurement seen between Planck’s two latest data releases results in a tighter constraint on polarization variance expectations. Finally, we consider even smaller uncertainties on τ and how more precise measurements of τ can drive the expectation for polarization variance in a hemisphere close to that of the cosmic-variance-limited distribution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A differentiable perturbation-based weak lensing shear estimator

Upcoming imaging surveys will use weak gravitational lensing to study the large-scale structure of the Universe, demanding sub-per cent accuracy for precise cosmic shear measurements. We present a new differentiable implementation of our perturbation-based shear estimator (fpfs), using jax, which is publicly available as part of a new suite of analytic shear algorithms called anacal. This code can analytically calibrate the shear response of any non-linear observable constructed with the fpfs shapelets and detection modes utilizing autodifferentiation (ad), generalizing the formalism to include a family of shear estimators with corrections for detection and selection biases. Using the ad capability of jax, it calculates the full Hessian matrix of the non-linear observables, which improves the previously presented second-order noise bias correction in the shear estimation. As an illustration of the power of the new anacal framework, we optimize the effective galaxy number density in the space of the generalized shear estimators using an LSST-like galaxy image simulation for the 10 yr LSST. For the generic shear estimator, the magnitude of the multiplicative bias |m| is below 3 × 10 –3 (99.7 per cent confidence interval), and the effective galaxy number density is improved by 5 per cent. We also discuss some planned future additions to the anacal software suite to extend its applicability beyond the fpfs measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Phenotyping stomatal closure by thermal imaging for GWAS and TWAS of water use efficiency-related genes

Abstract Stomata allow CO2 uptake by leaves for photosynthetic assimilation at the cost of water vapor loss to the atmosphere. The opening and closing of stomata in response to fluctuations in light intensity regulate CO2 and water fluxes and are essential for maintaining water-use efficiency (WUE). However, a little is known about the genetic basis for natural variation in stomatal movement, especially in C4 crops. This is partly because the stomatal response to a change in light intensity is difficult to measure at the scale required for association studies. Here, we used high-throughput thermal imaging to bypass the phenotyping bottleneck and assess 10 traits describing stomatal conductance (gs) before, during and after a stepwise decrease in light intensity for a diversity panel of 659 sorghum (Sorghum bicolor) accessions. Results from thermal imaging significantly correlated with photosynthetic gas exchange measurements. gs traits varied substantially across the population and were moderately heritable (h2 up to 0.72). An integrated genome-wide and transcriptome-wide association study identified candidate genes putatively driving variation in stomatal conductance traits. Of the 239 unique candidate genes identified with the greatest confidence, 77 were putative orthologs of Arabidopsis (Arabidopsis thaliana) genes related to functions implicated in WUE, including stomatal opening/closing (24 genes), stomatal/epidermal cell development (35 genes), leaf/vasculature development (12 genes), or chlorophyll metabolism/photosynthesis (8 genes). These findings demonstrate an approach to finding genotype-to-phenotype relationships for a challenging trait as well as candidate genes for further investigation of the genetic basis of WUE in a model C4 grass for bioenergy, food, and forage production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Three-dimensional lattice modulations in the charge density wave system Lu 2 Ir 3 Si 5

Using total and resonant x-ray scattering coupled to large-scale computer modeling, we study the lattice modulations in the complex charge density wave (CDW) material Lu 2 ⁢Ir 3 ⁢Si 5 . Here, we find that it is a unique quantum system where periodic lattice modulations related to emergent CDW order occur in three orthogonal atomic planes of the crystal lattice, leading to the emergence of an unusual three-dimensional (3D) pattern of short and long Ir-Ir and Lu-Lu bonds. The 3D character of observed lattice modulations explains the largely isotropic character of the changes in the electronic properties occurring when the CDW order sets in, demonstrating the strong electron-lattice coupling in Lu 2 ⁢Ir 3 ⁢Si 5 . The result is supported by DFT calculations based on the experimental structure data. Altogether, our work provides strong evidence for the presence of a relationship between the dimensionality of emergent lattice distortions and that of concurrent changes in the electronic properties of CDW materials. The relationship may need to be accounted for when these materials are explored for practical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Emergent surface resonance from charge density wave symmetry breaking in TiSe 2

Surface confined electronic states provide a fertile ground for discovering emergent phenomena that have no counterpart in the bulk, offering new routes to manipulate correlations, symmetry breaking, and dimensionality at the atomic scale. Here, in this study, we show that charge density wave (CDW) symmetry breaking can yield surface states in 1⁢T−TiSe 2 . Micro–angle-resolved photoemission spectroscopy (µ-ARPES) resolves a sharp, two-dimensional surface resonant state (SRS) that emerges within the CDW reconstructed low energy spectrum. The SRS exhibits notable temperature dependence and its spectral weight collapses around ∼ 160 K, while CDW transition temperature 𝑇 CDW is commonly reported as ≈ 202 K. Slab DFT + 𝑈 calculations reproduce a surface localized resonance when CDW folding brings valence and conduction states into near degeneracy, suggesting a correlation tuned, surface selective origin. These results point to a form of correlation-tuned surface resonance in a layered CDW compound and suggest a framework for engineering low-dimensional quantum states in van der Waals materials via symmetry breaking and electronic structure tuning.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Low-rank Sachdev-Ye-Kitaev models

Motivated by recent proposals of experimental realization of fast scramblers, we study a family of solvable variants of the ($\textit{q}$ = 4) Sachdev-Ye-Kitaev model in which the rank and eigenvalue distribution of the coupling matrix $J_{ ij,kl}$ are tuneable. When the rank is proportional to the number of fermions, the low temperature behavior is sensitive to the eigenvalue distribution. We obtain a complete classification of the possible non-Fermi liquid quantum phases. These include two previously studied phases whose fermion scaling dimension depends continuously on the rank; we show that they are maximally chaotic, but necessitate an extensively degenerate or negative semidefinite coupling matrix. In this work, more generic distributions give rise to “almost Fermi liquids” with a scaling dimension Δ = 1/2, but which differ from a genuine Fermi liquid in quasiparticle decay rate, quantum Lyapunov exponent, and/or specific heat.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Insights into negative differential resistance in MoS 2 Esaki diodes: A first-principles perspective

MoS 2 is a two-dimensional material with a band gap depending on the number of layers and tunable by an external electric field. The experimentally observed intralayer band-to-band tunneling and interlayer band-to-band tunneling in this material present an opportunity for new electronic applications in tunnel field-effect transistors. However, such a widely accepted concept has yet to be been supported by theoretical investigations based on first principles. In this paper, using density functional theory, in conjunction with nonequilibrium Green's function techniques and our electric field gating method, enabled by a large-scale computational approach, we study the relation between band alignment and transmission in planar and side-stack MoS 2 p–i–n junction configurations. Here, we demonstrate the presence of negative differential resistance for both in-plane and interlayer current, a staple characteristic of tunnel diode junctions, and analyze the physical origin of such an effect. Electrostatic potentials, the van der Waals barrier, and a complex band analysis are also examined for a thorough understanding of Esaki diodes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Optical magnons with dominant bond-directional exchange interactions in the honeycomb lattice iridate α — Li 2 IrO 3

We have used resonant inelastic x-ray scattering to reveal optical magnons in a honeycomb lattice iridate α - Li 2 IrO 3 . The spectrum in the energy region 20–25 meV exhibits momentum dependence, of which energy is highest at the location of the magnetic Bragg peak, ( h , k ) = ( ± 0.32 , 0 ) , and lowered toward ( 0 , 0 ) and ( ± 1 , 0 ) . We compare our data with a linear spin-wave theory based on a generic nearest-neighbor spin model. We find that a dominant bond-directional Kitaev interaction of order 20 meV is required to explain the energy scale observed in our study. The observed excitations are understood as stemming from optical magnon modes whose intensity is modulated by a structure factor, resulting in the apparent momentum dependence. We also observed diffuse magnetic scattering arising from the short-range magnetic correlation well above T N . In contrast to Na 2 IrO 3 , this diffuse scattering lacks the C 3 rotational symmetry of the honeycomb lattice, suggesting that the bond anisotropy is far from negligible in α - Li 2 IrO 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Isospin symmetry breaking in the mirror pair 73 Sr- 73 Br

The recent experimental observation of isospin symmetry breaking (ISB) in the ground states of the T=3/2 mirror pair 73 Sr- 73 Br is theoretically studied using large-scale shell-model calculations. The large valence space and the successful PFSDG-U effective interaction used for the nuclear part of the problem capture possible structural changes and provide a robust basis to treat the ISB effects of both electromagnetic and nonelectromagnetic origin. The calculated shifts and mirror-energy differences are consistent with the inversion of the I π =1/2 - ,5/2 - states between 73 Sr- 73 Br and suggest that the role played by the Coulomb interaction is dominant. Finally, an isospin breaking contribution of nuclear origin is estimated to be ≈25keV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Quantum Spin Probe of Single Charge Dynamics

Electronic defects in semiconductors form the basis for emerging quantum technologies, but many defect centers are difficult to access at the single-particle level. A method for probing optically inactive spin defects would reveal semiconductor physics at the atomic scale and advance the study of new quantum systems. Here, we exploit the intrinsic correlation between the charge and spin states of defect centers to measure the charge populations and dynamics of single substitutional nitrogen spin defects in diamond. By probing their steady-state spin population, read out at the single-defect level with a nearby nitrogen vacancy center, we directly measure the defect ionization-corroborated by first-principles calculations-an effect we do not have access to with traditional coherence-based quantum sensing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗