Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear approximation”

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

Evaluation and calibration of MERRA-2 and CAMS reanalysis for PM 2.5 in a semi-urbanized area in the south of the Amazon

Air pollution has significant implications for the climate and poses irreversible risks to human health. The Amazon region of Brazil is severely affected by biomass burning (BB) emissions, yet air quality monitoring remains highly inadequate. Given the scarcity of surface-based observations, reanalysis models have become essential tools for assessing air pollution. Although MERRA-2 and CAMS PM 2.5 products are widely utilized, their validation and comprehensive evaluation for the Amazon Basin remain limited. Here, this study assesses the performance of these products in a semi-urbanized region in the southern Amazon. The calibrated time series was employed to analyze PM 2.5 concentrations from 2003 to 2023. Our results showed satisfactory performance of both products for the 24-h averages of PM 2.5 , with linear correlations above 0.76. However, it was found that both products overestimate surface concentrations. MERRA-2 performed better, with approximately 30% lower bias than CAMS. Time series analysis showed that the study area is strongly impacted by emissions BB in the dry period, mainly in August and September. Furthermore, our findings indicate a positive trend in increasing PM 2.5 concentrations, with a notable rise observed since 2014. The average PM 2.5 levels frequently exceed the daily air quality guidelines established by the WHO in 2021. It has been estimated that the population of this region is exposed to concentrations above 15 μg.m -3 , on average, more than 30 days per year. Our results contribute to the evaluation of MERRA-2 and CAMS products for Amazon and provide a corrected estimate for surface PM 2.5 . Recent concerns about air quality and the implementation of new surface monitoring networks may improve the evaluation of reanalysis products. In the short term, the need for this information makes our assessments indispensable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Uncertainty-Guided Prediction Horizon of Phase-Resolved Ocean Wave Forecasting Under Data Sparsity: Experimental and Numerical Evaluation

Accurate short-term wave forecasting is critical for the safe and efficient operation of marine structures that rely on real-time, phase-resolved ocean wave information for control and monitoring purposes (e.g., digital twins). These systems often depend on environmental sensors (e.g., waverider buoys, wave-sensing LIDAR). Challenges arise when upstream sensor data are missing, sparse, or phase-shifted due to drift. This study investigates the performance of two machine learning models, time-series dense encoder (TiDE) and long short-term memory (LSTM), for forecasting phase-resolved ocean surface elevations under varying degrees of data degradation. We introduce the τ-trimming algorithm, which adapts the prediction horizon based on uncertainty thresholds derived from historical forecasts. Numerical wave tank (NWT) and wave basin experiments are used to benchmark model performance under short- and long-term data masking, spatially coarse sensor grids, and upstream phase shifts. Results show under a 50% probability of upstream data loss, the τ-trimmed TiDE model achieves a 46% reduction in error at the most upstream target, compared to 22% for LSTM. Furthermore, phase misalignment in upstream data introduces a near-linear increase in forecast error. Under moderate model settings, a ±3 s misalignment increases the mean absolute error by approximately 0.5 m, while the same error is accumulated at ±4 s using the more conservative approach. These findings inform the design of resilient, uncertainty-aware wave forecasting systems suited for realistic offshore sensing environments.

42 ENGINEERING↗
DOE OSTI2025

Multiqubit entanglement generation with squeezed modes

We present a hybrid continuous variable–discrete variable entanglement generation protocol using linear optics and homodyne measurements, capable of producing multiple high-fidelity Bell pairs per protocol iteration, with an approximate 0.5 success probability. The effectiveness of the protocol is determined by the squeezing strength. To increase the number of Bell pairs, approximately 3 dB of extra squeezing is needed for each additional Bell pair. The protocol also generates single Bell pairs with an approximate 0.75 probability for squeezing strengths ⪅ 15 dB , achievable with current technology.

Macridin, Alexandru [Fermilab] (ORCID:000000022228↗
DOE OSTI2024

Matter power spectra in modified gravity: a comparative study of approximations and N -body simulations

ABSTRACT Testing gravity and the concordance model of cosmology, $\Lambda$CDM, at large scales is a key goal of this decade’s largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon, and K-mouflage. In particular, we compare the predictions from N-body simulations solving the full scalar field equation, two N-body codes with approximate time integration schemes, a parametrized modified N-body implementation, and the analytic halo model reaction approach. We find the modification to the $\Lambda$CDM spectrum is in 2 per cent agreement at $z\le 1$ and $k\le 1~h\,{\rm Mpc}^{-1}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation.

Bose, B. (ORCID:0000000319658614)↗
DOE OSTI2025

Randomized algorithms for accelerating linear algebraic computations

The project supported the development of new methodologies for performing matrix computations that form key building blocks in modern scientific computing, such as low rank approximation of matrices, and efficient representations of global operators that arise in simulations of physical phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Active learning for SNAP interatomic potentials via Bayesian predictive uncertainty

Bayesian inference with a simple Gaussian error model is used to efficiently compute prediction variances for energies, forces, and stresses in the linear SNAP interatomic potential. Here, the prediction variance is shown to have a strong correlation with the absolute error over approximately 24 orders of magnitude. Using this prediction variance, an active learning algorithm is constructed to iteratively train a potential by selecting the structures with the most uncertain properties from a pool of candidate structures. The relative importance of the energy, force, and stress errors in the objective function is shown to have a strong impact upon the trajectory of their respective net error metrics when running the active learning algorithm. Batched training of different batch sizes is also tested against singular structure updates, and it is found that batches can be used to significantly reduce the number of retraining steps required with only minor impact on the active learning trajectory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The TDHF code Sky3D version 1.2

The Sky3D code has been widely used to describe nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The approach is based on Skyrme forces or related energy density functionals. The static and dynamic equations are solved on a three-dimensional grid, and pairing is been implemented in the BCS approximation. This updated version of the code aims to facilitate the calculation of nuclear strength functions in the regime of linear response theory, while retaining all existing functionality and use cases. The strength functions are benchmarked against available RPA codes, and the user has the freedom of choice when selecting the nature of external excitation (from monopole to hexadecapole and more). Some utility programs are also provided that calculate the strength function from the time-dependent output of the dynamic calculations of the Sky3D code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Chromatic Transport Models for ImpactX

This note describes the basic levels of Hamiltonian approximation (for straight-axis elements) that are used in the symplectic beam dynamics modeling code ImpactX. These include: purely linear models, models based on a paraxial (chromatic) expansion of the Hamiltonian, and models based on the exact nonlinear Hamiltonian. The focus here is on paraxial (chromatic) models.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Electronic Properties of Ultra‐Wide Bandgap B x Al 1− x N Computed from First‐Principles Simulations

Abstract Ultra‐wide bandgap (UWBG) materials such as AlN and BN hold great promise for future power electronics due to their exceptional properties. They exhibit large bandgaps, high breakdown fields, high thermal conductivity, and high mechanical strengths. AlN and BN have been extensively researched, however, their alloys, B x Al 1− x N, are much less studied despite their ability to offer tunable properties by adjusting x . In this article, the electronic properties of 17 recently predicted ground states of B x Al 1− x N in the x = 0 − 1 range are predicted using first‐principles density functional theory and many‐body perturbation theory within GW approximation. All the B x Al 1− x N structures are found to be UWBG materials and have bandgaps that vary linearly from that of wurtzite‐phase ( w ) AlN (6.19 eV) to that of w ‐BN (7.47 eV). The bandstructures of B x Al 1− x N show that a direct‐to‐indirect bandgap crossover occurs near x = 0.25. Furthermore, it is found that B x Al 1− x N alloys have much larger dielectric constants than the constituent bulk materials (AlN = 9.3 ɛ 0 or BN = 7.3 ɛ 0 ), with values reaching as high as 12.1 ɛ 0 . These alloys are found to exhibit large dielectric breakdown fields in the range 9–35 MV cm −1 with a linear dependence on x . This work provides the much needed advancement in the understanding of the properties of B x Al 1− x N to aid their application in next‐generation devices.

Milne, Cody L.↗
DOE OSTI2025

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thickness scaling and ferroelectric switching in wurtzite structure zinc magnesium oxide thin films

This study reports the thickness scaling of sputtered ferroelectric Zn 0.61 Mg 0.39 O (ZMO) thin films down to 43 nm. Encapsulated IrO 2 /ZMO/Ir capacitors exhibited switchable polarizations exceeding 50 μC cm −2 and coercive fields that increased from 3.9 to 4.4 MV cm −1 as the thickness decreased. Switching kinetics are best described by the simultaneous non-linear nucleation and the growth model. Bimodal switching is prevalent at low thicknesses, with the fastest switching times measured to be approximately 400 ns. Device encapsulation made ZMO switching kinetics more abrupt, potentially due to changes in the concentration of atmosphere-induced defects such as hydroxides. These results demonstrate stable ferroelectricity and sub-microsecond switching in sub-50 nm wurtzite ZMO, highlighting its potential as a low-voltage ferroelectric for integrated nonvolatile memory applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

WEC-Sim Modeling of Laminar Scientific's Patented Seesaw Wave Energy Converter: Preprint

Laminar Scientific's patented seesaw wave energy converter was modeled in WEC-Sim to predict performance. The device operates by utilizing ocean surface waves to rotate a truss in pitch about a pivot. The pivot is located at the top of two pylons, which are embedded in the seafloor. The seesaw has a float on either end, and the buoyancy forces from each float cause the system to rise or fall with passing waves. Device performance relies upon seesaw length and ocean wavelength creating an antiphase effect. The seesaw truss has an adjustable length intended to achieve this effect. The operation method enforces a narrow band of wavelengths which induce the largest rotational motion from the device. The hydrodynamic analysis of the device was performed using Capytaine, and the results confirmed that the device operates best in a narrow frequency band. Four float-to-float spacing cases and three pylon radii were examined. The hydrodynamic results indicate a match between the model and the physical expectations for the device, and that varying the pylon radii by 0.1-m increments for three instances creates minimal changes in hydrodynamic properties. Power matrices for three float spacing cases of the device were simulated with Joint North Sea Wave Project spectra waves and optimal power take-off damping in WEC-Sim. The maximum average power production for the 15-m spacing case was 14.1 kW with a 5.0-s peak wave period and 4-m significant wave heights. Plots of capture-width ratios indicated that the device performance was linear and confirmed that the device is optimal in a narrow frequency band. The maximum percentage of the available wave power produced by the 15-m device was approximately 16%. Simulations of the device in regular waves were used to produce plots of average power compared to a ratio of float spacing to wavelength. These plots indicate that the power production is maximized at a ratio of 0.5, and further confirm that the device has a narrow frequency response. The device was simulated at an example field location, where the device produced an annual average power rating of 1.6 kW given an average omnidirectional wave climate of 10.3 kW m-1 and an optimal, linearized power take-off model. While the maximum predicted device performance is reliant upon a narrow band of wave frequencies, the conducted analysis provides an opportunity to improve device design prior to prototyping and testing. Modifying the design to respond to a broader frequency range would improve device performance.

numerical modeling↗
DOE OSTI2024

Linear gyrokinetic simulations of toroidal Alfvén eigenmodes in the Mega-Amp Spherical Tokamak

Linear gyrokinetic (GK) simulations using the Gyrokinetic Toroidal Code (GTC) have been performed to investigate Toroidicity-driven Alfvén Eigenmodes (TAEs) driven by the neutral beam injection (NBI) induced fast ions in the Mega-Amp Spherical Tokamak (MAST) to identify the non-perturbative and kinetic effects of thermal plasma. A specific TAE in MAST discharge 26887, with an on-axis NBI power of approximately 1.5 MW and plasma current around 800 kA, exhibited frequency chirping, and the tangential soft x-ray camera array resolved the radial mode structure peaked near |q|=1.5. Various excitation methods were used in the GTC linear simulations, illustrating this code's capability to realistically represent the mechanisms and behaviors of fast ion-driven TAEs in spherical tokamaks. The radial structures from these GK simulations closely match measurements and calculations performed using the NOVA ideal MHD code, though with the frequencies approximately 10 kHz lower, likely due to various kinetic and non-perturbative effects. The simulations measured the damping rates due to continuum damping, radiative damping, and ion Landau damping, revealing that ion Landau damping has the most significant contribution to the total damping rate of the TAE. A comparison of growth rates of TAEs excited by fast ion Maxwellian and slowing-down distributions shows that the TAEs excited by a fast ion anisotropic pitch distribution (as part of the slowing-down distributions) are more unstable compared to those excited by a Maxwellian distribution with an equivalent fast ion beta. This shows that the use of fast ion anisotropy alters the number of fast ions to be in shear Alfvén resonance, and hence, it can greatly affect the stability of TAEs. These tests can be performed with the GTC but impossible with ideal MHD simulations, highlighting the necessity of kinetic simulations such as the GTC for a precise prediction of the TAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

NoRA: A Tensor Network Ansatz for Volume-Law Entangled Equilibrium States of Highly Connected Hamiltonians

Motivated by the ground state structure of quantum models with all-to-all interactions such as mean-field quantum spin glass models and the Sachdev-Ye-Kitaev (SYK) model, we propose a tensor network architecture which can accomodate volume law entanglement and a large ground state degeneracy. We call this architecture the non-local renormalization ansatz (NoRA) because it can be viewed as a generalization of MERA, DMERA, and branching MERA networks with the constraints of spatial locality removed. We argue that the architecture is potentially expressive enough to capture the entanglement and complexity of the ground space of the SYK model, thus making it a suitable variational ansatz, but we leave a detailed study of SYK to future work. We further explore the architecture in the special case in which the tensors are random Clifford gates. Here the architecture can be viewed as the encoding map of a random stabilizer code. We introduce a family of codes inspired by the SYK model which can be chosen to have constant rate and linear distance at the cost of some high weight stabilizers. We also comment on potential similarities between this code family and the approximate code formed from the SYK ground space.

Physics↗
DOE OSTI2025

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗
DOE OSTI2024

Gradient Coding With Iterative Block Leverage Score Sampling

Gradient coding is a method for mitigating straggling servers in a centralized computing network that uses erasure-coding techniques to distributively carry out first-order optimization methods. Randomized numerical linear algebra uses randomization to develop improved algorithms for large-scale linear algebra computations. In this study, we propose a method for distributed optimization that combines gradient coding and randomized numerical linear algebra. The proposed method uses a randomized ℓ 2 -subspace embedding and a gradient coding technique to distribute blocks of data to the computational nodes of a centralized network, and at each iteration the central server only requires a small number of computations to obtain the steepest descent update. The novelty of our approach is that the data is replicated according to importance scores, called block leverage scores, in contrast to most gradient coding approaches that uniformly replicate the data blocks. Furthermore, we do not require a decoding step at each iteration, avoiding a bottleneck in previous gradient coding schemes. We show that our approach results in a valid ℓ 2 -subspace embedding, and that our resulting approximation converges to the optimal solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Updimensioning strategy derived from synthetic equiaxed grain structures for approximating 3D grain size distributions from 2D visualizations with 1D parameters

We generated synthetic equiaxed grain structures using computer graphics software to explore the relationship between various grain size determination methods and true three-dimensional (3D) grain diameters. Mirroring grain measurement techniques, the synthetic 3D grain structures are imaged as 2D micrographs which are measured to yield 1D grain size parameters. Synthetic grain structures provide data at a mass scale and permit exploration of both polished and fractured surface micrographs, revealing one-to-one correspondence between exposed 2D grain cross-sections and individual 3D grains. Analysis of this correspondence yielded a procedure to approximate 3D equiaxed grain size and volume distributions based on the mode of the 2D fractograph grain size distribution. The 3D approximation procedure is shown to be less susceptible to different imaging conditions that affect small, undiscernible grains compared to the standard planimetric and linear intercept methods, which by design also tend to underestimate the 3D grain diameter. The procedure requires larger sample sizes to lower variance and a deeper analysis which could become more practical with machine learning (ML) models for grain boundary segmentation, which synthetic grain structures can help train. This work lays the foundation for analyzing other grain distributions such as columnar and composite grains in similar depth.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Photoionization of the O 2 molecule

Ground state total and partial photoionization cross sections of the O 2 molecule have been produced with the configuration-average distorted-wave (CADW) approximation, and an energy-corrected pseudo independent-atom (PIA) method. The PIA approach is based on molecular orbitals expressed as a linear combination of atomic orbitals, and utilizes the corresponding atomic-orbital photoionization cross sections. Here we use an energy-corrected PIA approach with atomic photoionization cross sections computed using a CADW method. Our results are found to be in good agreement with photoionization measurements across a wide photon energy range, and in reasonably good agreement with the CADW calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗