Search NASA⌕ Search

SEARCH · Search NASA

Results for “Error Fields”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Stochastic tensor contraction for quantum chemistry

Many computational methods in ab initio quantum chemistry are formulated in terms of high-order tensor contractions, whose cost determines the size of system that can be studied. We introduce stochastic tensor contraction to perform such operations with greatly reduced cost, and present its application to the gold-standard quantum chemistry method, coupled cluster theory with up to perturbative triples. For total energy errors more stringent than chemical accuracy, we reduce the computational scaling to that of mean-field theory, while starting to approach the mean-field absolute cost, thereby challenging the existing cost-to-accuracy landscape. Benchmarks against state-of-the-art local correlation approximations further show that we achieve an order-of-magnitude improvement in both total computation time and error, with significantly reduced sensitivity to system dimensionality and electron delocalization. We conclude that stochastic tensor contraction is a powerful computational primitive to accelerate a wide range of quantum chemistry.

Chemical Physics (physics.chem-ph)↗
DOE OSTI2022

Fast cycling HTS based superconducting accelerator magnets: Feasibility study and readiness demonstration program driven by neutrino physics and muon collider needs

Development of energy-efficient fast cycling accelerator magnets is critical for the next generation of proton rapid cycling synchrotrons (RCS) for neutrino research and booster accelerators of future muon colliders. We see a unique opportunity for having such magnets to be built on base of High Temperature Superconductors (HTS). Besides being superconducting at relatively high temperatures, rare-earth HTS tapes have shown very small AC losses compared to low-temperature NbTi superconductor cables. Recent tests of the HTS-based 0.5 m long two-bore superconducting accelerator magnet have shown record high dB/dt ramping rates of about 300 T/s at 10 Hz repetition rate and 0.5 T B-field span. No temperature rise in 6 K cooling He was observed within the 0.003 K error setting the upper limit on the cryogenic power loss in the magnet conductor coil to be less than 0.2 W/m. Based on this result we outline a possible upgrade of this test magnet design to 2 T B-field in the 10 mm beam gap with the dB/dt ramping rates up to 1000 T/s. The power test results of this short sample magnet will be used to project both cryogenic and electrical power losses as a function of the magnet B-field and the dB/dt ramping rates. Then these projections will be scaled to the range of expected accelerator magnet beam gaps and B-fields for the proton and muon RCS accelerators. We invite collaborators to join these studies and call for support of the R&D program aimed at comprehensive demonstration of this approach that includes design, construction, and power tests of a long prototype of the HTS-based fast-cycling accelerator magnet by 2028.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Optimization of algorithmic errors in analog quantum simulations

Due to rapidly improving quantum computing hardware, Hamiltonian simulations of relativistic lattice field theories have seen a resurgence of attention. Furthermore, this computational tool requires turning the formally infinite-dimensional Hilbert space of the full theory into a finite-dimensional one. For gauge theories, a widely used basis for the Hilbert space relies on the representations induced by the underlying gauge group, with a truncation that keeps only a set of the lowest dimensional representations. This works well at large bare gauge coupling, but becomes less efficient at small coupling, which is required for the continuum limit of the lattice theory. In this work, we develop a new basis suitable for the simulation of an SU(2) lattice gauge theory in the maximal tree gauge. In particular, we show how to perform a Hamiltonian truncation so that the eigenvalues of both the magnetic and electric gauge-fixed Hamiltonian are mostly preserved, which allows for this basis to be used at all values of the coupling. Little prior knowledge is assumed, so this may also be used as an introduction to the subject of Hamiltonian formulations of lattice gauge theories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Density weighted angular redshift fluctuations: a new cosmological observable

We propose the use of angular fluctuations in the galaxy redshift field as a new way to extract cosmological information in the Universe. This new probe $\delta z (\hat{\mathbf {n}})$ consists of the statistics of sky maps built by projecting redshifts under a Gaussian window of width σz centred upon a redshift zobs, and weighted by the galaxy density field. We compute the angular power spectrum of the $\delta z (\hat{\mathbf {n}})$ field in both numerical simulations and in linear perturbation theory. From these, we find that the $\delta z (\hat{\mathbf {n}})$ field (i) is sensitive to the underlying density and peculiar velocity fields; (ii) is highly correlated, at the $\gtrsim 60\, {{\ \rm per\ cent}}$ level, to the line-of-sight projected peculiar velocity field; (iii) for narrow windows (σz < 0.03), it is almost completely uncorrelated to the projected galaxy angular density field under the same redshift window; and (iv) it is largely unaffected by multiplicative and additive systematic errors on the observed number of galaxies that are redshift-independent over ~σz. We conclude that $\delta z (\hat{\mathbf {n}})$ is a simple and robust tomographic measure of the cosmic density and velocity fields, complementary to angular clustering, that will contribute to more complete exploitations of current and upcoming galaxy redshift surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Mesh-based super-resolution of fluid flows with multiscale graph neural networks

A graph neural network (GNN) approach is introduced in this work which enables mesh-based three-dimensional super-resolution of fluid flows. In this framework, the GNN is designed to operate not on the full mesh-based field at once, but on localized meshes of elements (or cells) directly. To facilitate mesh-based GNN representations in a manner similar to spectral (or finite) element discretizations, a baseline GNN layer (termed a message passing layer, which updates local node properties) is modified to account for synchronization of coincident graph nodes, rendering compatibility with commonly used element-based mesh connectivities. Furthermore, the architecture is multiscale in nature, and is comprised of a combination of coarse-scale and fine-scale message passing layer sequences (termed processors) separated by a graph unpooling layer. The coarse-scale processor embeds a query element (alongside a set number of neighboring coarse elements) into a single latent graph representation using coarse-scale synchronized message passing over the element neighborhood, and the fine-scale processor leverages additional message passing operations on this latent graph to correct for interpolation errors. Demonstration studies are performed using hexahedral mesh-based data from Taylor–Green Vortex and backward-facing step flow simulations at Reynolds numbers of 1600 and 3200. Through analysis of both global and local errors, the results ultimately show how the GNN is able to produce accurate super-resolved fields compared to targets in both coarse-scale and multiscale model configurations. Reconstruction errors for fixed architectures were found to increase in proportion to the Reynolds number. Geometry extrapolation studies on a separate cavity flow configuration show promising cross-mesh capabilities of the super-resolution strategy.

Backward-facing step↗
DOE OSTI2025

High order interpolation of magnetic fields with vector potential reconstruction for particle simulations

We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary differential equation (ODE) integrators, such as the Dormand-Prince method. The method provides C(m) continuity and achieves high-order accuracy, making it particularly suited for particle trajectory integration and Poincaré section analysis under optimal integration order and timestep adjustments. Through numerical experiments, we demonstrate that the Hermite interpolation method preserves volume and continuity, which are critical for conserving toroidal canonical momentum and magnetic moment in guiding center simulations, especially over long-term trajectory integration. Furthermore, we analyze the impact of insufficient derivative continuity on Runge-Kutta schemes and show how it degrades accuracy at low error tolerances, introducing discontinuity-induced truncation errors. Lastly, we demonstrate performant Poincaré section analysis in two relevant settings of field data collocated from finite element meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

TopoSZ: Preserving Topology in Error-Bounded Lossy Compression

Existing error-bounded lossy compression techniques control the pointwise error during compression to guarantee the integrity of the decompressed data. However, they typically do not explicitly preserve the topological features in data. When performing post hoc analysis with decompressed data using topological methods, preserving topology in the compression process to obtain topologically consistent and correct scientific insights is desirable. In this paper, we introduce TopoSZ, an error-bounded lossy compression method that preserves the topological features in 2D and 3D scalar fields. Specifically, we aim to preserve the types and locations of local extrema as well as the level set relations among critical points captured by contour trees in the decompressed data. The main idea is to derive topological constraints from contour-tree-induced segmentation from the data domain, and incorporate such constraints with a customized error-controlled quantization strategy from the SZ compressor (version 1.4). In conclusion, our method allows users to control the pointwise error and the loss of topological features during the compression process with a global error bound and a persistence threshold.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

High-Accuracy Mesh-Free Quadrature for Trimmed Parametric Surfaces and Volumes

This work presents a high-accuracy, mesh-free, generalized Stokes theorem-based numerical quadrature scheme for integrating functions over trimmed parametric surfaces and volumes. The algorithm relies on two fundamental steps: (1) We iteratively reduce the dimensionality of integration using the generalized Stokes theorem to line integrals over trimming curves, and (2) we employ numerical antidifferentiation in the generalized Stokes theorem using high-order quadrature rules. Additionally, the scheme achieves exponential convergence up to trimming curve approximation error and has applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order curvilinear meshes, and initialization of multi-material simulations. We compare the quadrature scheme to commonly-used quadrature schemes in the literature and show that our scheme is much more efficient in terms of number of quadrature points used. We provide an open-source implementation of the scheme in MATLAB as part of QuaHOG, a software package for Quadrature of High-Order Geometries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Three-dimensional traction microscopy accounting for cell-induced matrix degradation

Tractions exerted by cells on the extracellular matrix (ECM) are critical in many important physiological and pathological processes such as embryonic morphogenesis, wound healing, and cancer metastasis. Three-dimensional Traction Microscopy (3DTM) is a tool to quantify cellular tractions by first measuring the displacement field in the ECM in response to these tractions, and then using this measurement to infer tractions. Most applications of 3DTM have assumed that the ECM has spatially-uniform mechanical properties, but cells secrete enzymes that can locally degrade the ECM. In this work, a novel computational method is developed to quantify both cellular tractions and ECM degradation. In particular, the ECM is modeled as a hyperelastic, Neo-Hookean solid, whose material parameters are corrupted by a single degradation parameter. The feasibility of determining both the traction and the degradation parameter is first demonstrated by showing the existence and uniqueness of the solution. An inverse problem is then formulated to determine the nodal values of the traction vector and the degradation parameter, with the objective of minimizing the difference between a predicted and measured displacement field, under the constraint that the predicted displacement field satisfies the equation of equilibrium. The inverse problem is solved by means of a gradient-based optimization approach, and the gradient is computed efficiently using appropriately derived adjoint fields. The computational method is validated in-silico using a geometrically realistic neuronal cell model and synthetic traction and degradation fields. It is found that the method accurately recovers both the traction and degradation fields. Moreover, it is found that neglecting ECM degradation can yield significant errors in traction measurements. Finally, our method can extend the range of context where tractions can be appropriately measured.

42 ENGINEERING↗
DOE OSTI2026

Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems

The Material Point Method (MPM) has been increasingly used to simulate large strain deformations. Linear interpolation functions are commonly used to perform the spatial integration. It is well-known that the discontinuities in the interpolation function derivatives induce shock-like artifacts known as ‘cell-crossing’ error. These errors compound with volumetric locking errors when used with hydromechanical formulations for porous media, where different velocity fields are used for each phase. The capabilities of higher-order MPM frameworks have not been explored for real-scale geotechnical problems. As such, this paper aims to assess, validate, and further discuss a higher-order B-spline MPM (BS-MPM) framework. First, the BS-MPM framework is verified against the large-strain oedometer consolidation problem. Second, the framework is validated against a real-scale slope failure experiment triggered by pore water pressure recharge. Landslide features that are captured using the higher-order framework are specifically highlighted, and results (e.g., pore water pressure and deformation) are validated with field measurements. A generally convergent numerical solution is observed when using cubic interpolation functions. Third, a footing penetration problem is simulated using the multi-patch BS-MPM. Trends are examined with respect to penetration velocity and variation in hydraulic conductivity. The BS-MPM framework ultimately presents a stabilized numerical solution that captures plausible hydromechanical interaction trends important in geotechnical engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A perspective on single-crystal layered oxide cathodes for lithium-ion batteries

As the demand for lithium-ion batteries grows exponentially to feed the nascent electric-vehicle and grid-storage markets, the need for higher energy density and longer cycle life becomes more apparent. Increasing the nickel content in the layered-oxide cathodes has been a dominant strategy to increase energy density, but this has exacerbated the surface reactivity concerns. Furthermore, the high states of charge associated with increased capacity lead to cracking of the polycrystalline cathode particles, which exposes fresh surfaces and accelerates capacity fade. Single-crystal cathode particles with low susceptibility to cracking have recently taken on intense interest due to their remarkable cyclability. Many different approaches to single-crystal synthesis have been explored, but there has been no systematic analysis of the different techniques to date. Additionally, many reports have taken on a trial-and-error approach, ignoring well established literature on grain growth outside the battery field. This perspective reviews the literature on single-crystal cathode synthesis through the lens of classic grain growth theories, providing insight and hopefully accelerating the adoption of the promising single-crystal cathode morphology.

25 ENERGY STORAGE↗
DOE OSTI2020

Assessing Collision Cross Section Calibration Strategies for Traveling Wave-Based Ion Mobility Separations in Structures for Lossless Ion Manipulations

The collision cross section (CCS) is an important property that aids in the structural characterization of molecules. Here, we investigated the CCS calibration accuracy with traveling waves ion mobility spectrometry (TWIMS) separations in structures for lossless ion manipulations (SLIM) using three sets of calibrants. A series of singly negatively-charged phospholipids and bile acids were calibrated in nitrogen buffer gas using two different TW waveform profiles (square and sine) and amplitudes (20 V0-p, 25 V0-p, and 30 V 0-p ). The calibration errors for the three calibrant sets Agilent tuning mixture and polyalanine, and one assembled in-house showed negligible differences using a sine-shaped TW waveform. Furthermore, calibration errors were all within 1-2 % of the DTIMS measurements, with lower errors for sine waveforms, presumably due to the lower average and maximum fields experienced by ions. Finally, ultrahigh resolution multi-pass (long path length) SLIM TWIMS demonstrated improved CCS calibration for phospholipid and bile acid isomers.

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

Symmetrized Drude Oscillator Force Fields Improve Numerical Performance of Polarizable Molecular Dynamics

Drude oscillator potentials are a popular and computationally efficient class of polarizable models that represent each polarizable atom as a positively charged Drude core harmonically bound to a negatively charged Drude shell. Here we show that existing force fields that place all non-Coulomb forces on the Drude core and none on the shell inadvertently couple the dipole to non-Coulombic forces. This introduces errors where interactions with neutral particles can erroneously induce atomic polarization, leading to spurious polarizations in the absence of an electric field, exacerbating violations of equipartition in the employed Carr–Parinello scheme. A suitable symmetrization of the interaction potential that correctly splits the force between the Drude core and shell can correct this shortcoming, improving the stability and numerical performance of Drude oscillator-based simulations. The symmetrization procedure is straightforward and only requires the rescaling of a few force field parameters.

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

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

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

Evaluation of dynamic testing of full-scale wind turbine drivetrains with hardware-in-the-loop

Dynamometer testing of full-scale wind turbine drivetrains involves subjecting integrated components to extreme static loads, fatigue loads, or dynamic loads to replicate field conditions. This paper provides a detailed evaluation of the uncertainty and errors in applied loads and measured responses, and their magnitude relative to established repeatability bounds for a drivetrain subjected to static and dynamic loads on the 7.5 MW test bench at the Clemson University Wind Turbine Drivetrain Testing Facility. An ideal drivetrain simulation model is utilized to isolate the influence of load tracking errors on test article responses. It is shown that the test article response variations are mainly driven by run out and clearances inherent to the drivetrain. For the dynamic profiles, the load tracking errors are within acceptable limits. Furthermore, a frequency analysis is used to show that the test bench controller tracking performance is acceptable for profiles with frequency content up to 1.5 Hz.

17 WIND ENERGY↗
DOE OSTI2025

High-Spectral Resolution Lidar (HSRL) Instrument Handbook

High-spectral-resolution lidar (HSRL) systems provide vertical profiles of optical depth, backscatter cross-section, depolarization, and backscatter phase function. All HSRL measurements are calibrated by reference to molecular scattering, which is measured at each point in the lidar profile. Like the Raman lidar (RL), but unlike simple backscatter lidars such as the micropulse lidar (MPL), this enables the HSRL to measure backscatter cross-sections and optical depths without prior assumptions about the scattering properties of the atmosphere. The depolarization observations allow robust discrimination between ice and water clouds. Rigorous error estimates can be computed for all measurements. A very narrow angular field of view reduces multiple scattering contributions. The small field of view, coupled with a narrow optical bandwidth, nearly eliminates noise due to scattered sunlight. The laser transmitter is a diode-pumped, frequency-doubled Nd:YAG laser. Narrow-band, single-frequency operation is provided by injection seeding with a single-frequency, cw-diode-pumped diode laser. The main laser cavity is maintained in resonance with the seed laser by adjusting the cavity length to minimize the time between the Q-switch trigger and the emission of the laser pulse. The emission wavelength is tuned via temperature control of the seed laser crystal and is locked to line #1109 of the iodine absorption spectra. Locking is accomplished by minimizing the transmission through a 2-cm-long iodine absorption cell. Use of a high-repetition-rate laser and expansion of the transmitted beam through a 400-mm telescope reduces the transmitted energy density to eye-safe levels. It is possible to look directly into the output beam without hazard. The receiver and transmitter use the same afocal telescope, simplifying the maintenance of stable alignment of the transmitter and receiver although the angular FOV is only 100 μrad. The small FOV and the 4-kHz repetition rate also limit the near-field signal strength, making it possible to record continuous profiles that start at an altitude of ~100 m and extend to 30 km using photon counting detectors. The small FOV also suppresses multiple scattering contributions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Towards Practical Quantum Simulation for High Energy Physics (Final report for DE-SC0019452)

Quantum computing offers the potential of transformative improvements in our ability to simulate strongly interacting quantum systems. Each application area of quantum simulation requires careful study and new insights to determine the potential advantages of quantum computation. In particular, high energy physics and quantum field theory present specific challenges to existing approaches. In this project we studied quantum simulation in the light-front formulation of quantum field theory. In this formulation, quantum field theory more closely resembles other many-body quantum systems such as quantum chemistry, which are well studied from the perspective of quantum computation. We considered applications suitable for both near term and future error-corrected quantum computers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗