Search NASA⌕ Search

SEARCH · Search NASA

Results for “INTERPOLATION”

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

ReMU: regional minimal updating for model-based derivative-free optimization

Derivative-free optimization (DFO) problems are optimization problems where derivative information is unavailable or extremely difficult to obtain. Model-based DFO solvers have been applied extensively in scientific computing. Powell's NEWUOA (2004) [Powell, The NEWUOA software for unconstrained optimization without derivatives, in Large-Scale Nonlinear Optimization, Nonconvex Optimization and its Applications Vol. 83, G. Di Pillo and M. Roma, eds., Springer, 2006, pp. 255–297] and Wild's POUNDerS (2014) [Wild, Solving derivative-free nonlinear least squares problems with POUNDERS, in Advances and Trends in Optimization with Engineering Applications, T. Terlaky, M.F. Anjos, and S. Ahmed, eds., SIAM, 2017, pp. 529–540] explore the numerical power of the minimal norm Hessian (MNH) model for DFO and contributed to the open discussion on building better models with fewer data to achieve faster numerical convergence. Another decade later, we propose the regional minimal updating (ReMU) models, and extend the previous models into a broader class, including the H 2 norm models [Xie and Yuan, Least H 2 norm updating of quadratic interpolation models for derivative-free trust-region algorithms, IMA J. Numer. Anal. 46 (2025), pp. 21–50]. This paper shows motivation behind ReMU models, computational details, theoretical and numerical results on particular extreme points and the barycentre of ReMU's weight coefficient region, and the associated KKT matrix error and distance. Novel metrics, such as the truncated Newton step error, are proposed to numerically understand the new models' properties. A new algorithmic strategy, based on iteratively adjusting the ReMU model type, is also proposed, and shows numerical advantages by combining and switching between the barycentric model and the classic least Frobenius norm model in an online fashion.

derivative-free trust-region methods↗
NASA NTRSDate not supplied

Parahydrogen Properties Version 05 Database Release for Nuclear Thermal Propulsion Applications

Consistent modeling assumptions across any large project are crucial to minimize errors between different approaches. Use of consistent material and fluid properties across a large project supports consistent interpretation and application within modeling and simulation results as well as their relevancy to operational systems. NASA’s Space Nuclear Propulsion program dedicates extensive resources towards establishing consistent and, to the extent possible, accurate property databases for its internal staff and all external partners. This work highlights the extensive research performed to modernize the fluid property database of hydrogen which is the leading propellant option for in-space nuclear propelled spacecraft. Specifically, the parahydrogen spin state is of interest since the propellant is stored in a near normal boiling point liquid state which results in it consisting almost entirely of the parahydrogen spin isomer. This database tool has taken recent NASA work and modernized it into a python-based package for easy usage across all modeling entities. The package allows users to provide temperature and pressure pairs along with their desired output properties to yield results accounting for both real-gas and equilibrium dissociation effects while also sharing the default thermodynamic reference state provided by the National Institute of Standards and Technology (NIST) Standard Database 23. The suite also includes advanced capabilities to increase usability, such as on-the-fly interpolation and multidimensional plotting.

Nuclear Thermal Propulsion↗
NASA NTRSDate not supplied

Parahydrogen Properties Version 05 Database Release for NTP Applications

Consistent modeling assumptions across any large project are crucial to minimize errors between different approaches. Use of consistent material and fluid properties across a large project supports consistent interpretation and application within modeling and simulation results as well as their relevancy to operational systems. NASA’s Space Nuclear Propulsion program dedicates extensive resources towards establishing consistent and, to the extent possible, accurate property databases for its internal staff and all external partners. This work highlights the extensive research performed to modernize the fluid property database of hydrogen which is the leading propellant option for in-space nuclear propelled spacecraft. Specifically, the parahydrogen spin state is of interest since the propellant is stored in a near normal boiling point liquid state which results in it consisting almost entirely of the parahydrogen spin isomer. This database tool has taken recent NASA work and modernized it into a python-based package for easy usage across all modeling entities. The package allows users to provide temperature and pressure pairs along with their desired output properties to yield results accounting for both real-gas and equilibrium dissociation effects while also sharing the default thermodynamic reference state provided by the National Institute of Standards and Technology (NIST) Standard Database 23. The suite also includes advanced capabilities to increase usability, such as on-the-fly interpolation and multidimensional plotting.

Nuclear Thermal Propulsion↗
DOE OSTI2025

Recommendations for developing, documenting, and distributing data products derived from NEON data

The National Ecological Observatory Network (NEON) provides over 180 distinct data products from 81 sites (47 terrestrial and 34 freshwater aquatic sites) within the United States and Puerto Rico. These data products include both field and remote sensing data collected using standardized protocols and sampling schema, with centralized quality assurance and quality control (QA/QC) provided by NEON staff. Such breadth of data creates opportunities for the research community to extend basic and applied research while also extending the impact and reach of NEON data through the creation of derived data products—higher level data products derived by the user community from NEON data. Derived data products are curated, documented, reproducibly-generated datasets created by applying various processing steps to one or more lower level data products—including interpolation, extrapolation, integration, statistical analysis, modeling, or transformations. Derived data products directly benefit the research community and increase the impact of NEON data by broadening the size and diversity of the user base, decreasing the time and effort needed for working with NEON data, providing primary research foci through the development via the derivation process, and helping users address multidisciplinary questions. Creating derived data products also promotes personal career advancement to those involved through publications, citations, and future grant proposals. However, the creation of derived data products is a nontrivial task. Here we provide an overview of the process of creating derived data products while outlining the advantages, challenges, and major considerations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sensitivity of Fine‐Resolution Urban Heat Island Simulations to Soil Moisture Parameterization

ABSTRACT Urban areas experience the impact of natural disasters, such as heatwaves and flash floods, disparately in different neighbourhoods across a city. The demand for precise urban hydrometeorological and hydroclimatological modelling to examine this disparity, and the interacting challenges posed by climate change and urbanisation, has thus surged. The Weather Research and Forecasting (WRF) model has served such operational and research purposes for decades. Recent advancements in WRF, including enhanced numerical schemes and sophisticated urban atmospheric‐hydrological parameterizations, have empowered the simulation of urban geophysical processes at high resolution (~1 km), but even this resolution misses significant urban microclimate variability. This study applies the large‐eddy simulations (LES) mode within WRF, coupled with single‐layer urban canopy models (SLUCM), to enable even finer‐scale modelling (150 m) of the Urban Heat Island (UHI) effect in the Baltimore metropolitan area. We run nine scenarios to evaluate various methods of initializing soil moisture and various spinup lead times, and to assess the impact of WRF's Mosaic approach in depicting subgrid‐scale processes. We evaluate the scenarios by comparing the WRF simulated land surface temperature (LST) against Landsat LST and the WRF simulated hourly 2‐m air temperatures (AT) with observations from eight weather stations across the domain. Results underscore the paramount influence of the lead spinup time on the spatiotemporal distribution of simulated soil moisture, consequently shaping WRF's efficacy in predicting the UHI. Furthermore, interpolating soil moisture‐related parameters from the parent for child domain initialization yields a notable reduction in mean and root‐mean‐squared errors. This improvement was particularly evident in simulations with the longest spinup time, affirming the importance of carefully designing the initialization of soil moisture for improved urban temperature predictions.

Talebpour, Mahdad↗
DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
DOE OSTI2025

Central Compact Finite‐Difference Scheme With High Spectral Resolution for KdV Equation

This work presents a combination of cell‐node and cell‐centered compact finite difference scheme for the approximation of third derivatives involved in Korteweg–de Vries (KdV) equations. This approach employs a half‐shifted derivative construction at cell centers, avoiding the need for compact interpolation, thereby removing transfer errors; hence, it improves spectral resolution and maintains high‐order accuracy. Fourier analysis is performed to show the spectral properties of the proposed formulation, which provides higher spectral resolutions as compared to node‐based compact schemes. A filtering strategy is incorporated to suppress high‐frequency oscillations without compromising the accuracy of the numerical scheme, and the total variation diminishing Runge Kutta (TVDRK3) method is applied for time integration. Numerical experiments on linear, nonlinear, and coupled KdV systems are conducted, and a comparative analysis with cell‐node compact schemes confirms that the proposed scheme consistently reduces errors by up to an order of magnitude and achieves high spectral resolution properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Randomized Algorithms for Low-Rank Matrix and Tensor Decompositions

This paper surveys randomized algorithms in numerical linear algebra for low-rank decompositions of matrices and tensors. The survey begins with a review of classical matrix algorithms that can be accelerated by randomized dimensionality reduction, such as the singular value decomposition (SVD) or interpolative (ID) and CUR decompositions. Recent advances in randomized dimensionality reduction are discussed, including new methods of fast matrix sketching and sampling techniques, which are incorporated into classical matrix algorithms for fast low-rank matrix approximations. The extension of randomized matrix algorithms to tensors is then explored for several low-rank tensor decompositions in the CP and Tucker formats, including the higher-order SVD, ID, and CUR decomposition.

Pearce, Katherine J. [The University of Texas at A↗
DOE OSTI2026

Graph theory inspired anomaly detection at the LHC

Designing model-independent anomaly detection algorithms for analyzing LHC data remains a central challenge in the search for new physics, due to the high dimensionality of collider events. In this work, we develop a graph autoencoder as an unsupervised, model-agnostic tool for anomaly detection, using the LHC Olympics dataset as a benchmark. By representing jet constituents as a graph, we introduce a method to systematically control the information available to the model through sparse graph constructions that serve as physically motivated inductive biases. Specifically, (1) we construct graph autoencoders based on locally rigid Laman graphs and globally rigid unique graphs, and (2) we explore the clustering of jet constituents into subjets to interpolate between high- and low-level input representations. We obtain the best performance, measured in terms of the Significance Improvement Characteristic curve for an intermediate level of subjet clustering and certain sparse unique graph constructions. We further investigate the role of graph connectivity in jet classification tasks. Our results demonstrate the potential of leveraging graph-theoretic insights to refine and increase the interpretability of machine learning tools for collider experiments.

Automation↗
DOE OSTI2025

Spinning black hole scattering at $ \mathcal{O} $(G 3 S 2 )(G 3 S 2 ): Casimir terms, radial action and hidden symmetry

We resolve subtleties in calculating the post-Minksowskian dynamics of binary systems, as a spin expansion, from massive scattering amplitudes of fixed finite spin. In particular, the apparently ambiguous spin Casimir terms can be fully determined from the gradient of the spin-diagonal part of the amplitudes with respect to S 2 = −s(s+1)ħ 2 , using an interpolation between massive amplitudes with different spin representations. From two-loop amplitudes of spin-0 and spin-1 particles minimally coupled to gravity, we extract the spin Casimir terms in the conservative scattering angle between a spinless and a spinning black hole at $ \mathcal{O} $(G 3 S 2 ), finding agreement with known results in the literature. This completes an earlier study [Phys. Rev. Lett. 130 (2023), 021601] that calculated the non-Casimir terms from amplitudes. We also illustrate our methods using a model of spinning bodies in electrodynamics, finding agreement between scattering amplitude predictions and classical predictions in a root-Kerr electromagnetic background up to $ \mathcal{O} $(α 3 S 2 ). For both gravity and electrodynamics, the finite part of the amplitude coincides with the two-body radial action in the aligned spin limit, generalizing the amplitude-action relation beyond the spinless case. Surprisingly, the two-loop amplitude displays a hidden spin-shift symmetry in the probe limit, which was previously observed at one loop. We conjecture that the symmetry holds to all orders in the coupling constant and is a consequence of integrability of Kerr orbits in the probe limit at the first few orders in spin.

Classical Theories of Gravity↗
DOE OSTI2025

Broken conformal window

We show that near the edges of the conformal window of supersymmetric SU(Nc) QCD, perturbed by Anomaly Mediated Supersymmetry Breaking (AMSB), chiral symmetry can be broken depending on the initial conditions of the RG flow. We do so by perturbatively expanding around Banks-Zaks fixed points and taking advantage of Seiberg duality. Interpolating between the edges of the conformal window, we predict that non-supersymmetric QCD breaks chiral symmetry up to N f ≤ 3N c − 1, while we cannot say anything definitive for N f ≥ 3N c at this moment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Thermal bootstrap of matrix quantum mechanics

We implement a bootstrap method that combines stationary state conditions, thermal inequalities, and semidefinite relaxations of matrix logarithm in the ungauged one-matrix quantum mechanics, at finite rank N as well as in the large N limit, and determine finite temperature observables that interpolate between available analytic results in the low and high temperature limits respectively. We also obtain bootstrap bounds on thermal phase transition as well as preliminary results in the ungauged two-matrix quantum mechanics.

1/N Expansion↗
DOE OSTI2024

Conservative scattering of Reissner-Nordström black holes at third post-Minkowskian order

Using a recently developed effective field theory formalism for extreme mass ratios arXiv:2308.14832, we present a calculation of charged black hole scattering at third post-Minkowskian order. The charges and masses are kept arbitrary, and the result interpolates from the scattering of Schwarzschild to extremal charged black holes, and beyond to charged particles in electrodynamics — agreeing with previously reported results in all such limits. The computation of the radial action is neatly organized in powers of the mass ratio. The probe (0SF) contributions are readily computed by direct integration of the radial momentum, and we use the effective field theory to compute the subleading (1SF) contributions via background-field Feynman rules supplemented by an operator encoding recoil of the background. Together these contributions completely determine the conservative physics at order O (G 3 ). .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Fuzzy spheres in stringy matrix models: quantifying chaos in a mixed phase space

We consider a truncation of the BMN matrix model to a configuration of two fuzzy spheres, described by two coupled non-linear oscillators dependent on the mass parameter μ. The classical phase diagram of the system generically (μ ≠ 0) contains three equilibrium points: two centers and a center-saddle; as μ → 0 the system exhibits a pitchfork bifurcation. We demonstrate that the system is exactly integrable in quadratures for μ = 0, while for very large values of μ, it approaches another integrable point characterized by two harmonic oscillators. The classical phase space is mixed, containing both integrable islands and chaotic regions, as evidenced by the classical Lyapunov spectrum. At the quantum level, we explore indicators of early and late time chaos. The eigenvalue spacing is best described by a Brody distribution, which interpolates between Poisson and Wigner distributions; it dovetails, at the quantum level, the classical results and reemphasizes the notion that the quantum system is mixed. We also study the spectral form factor and the quantum Lyapunov exponent, as defined by out-of-time-ordered correlators. These two indicators of quantum chaos exhibit weak correlations with the Brody distribution. We speculate that the behavior of the system as μ → 0 dominates the spectral form factor and the quantum Lyapunov exponent, making these indicators of quantum chaos less effective in the context of a mixed phase space.

AdS-CFT correspondence↗
DOE OSTI2025

Supersymmetry breaking in SYK and the black hole spectrum

The spectrum and dynamics of near-extremal black holes is strongly modified by quantum effects at low temperatures. When the extremal limit does not preserve any supersymmetry, the density of states goes to zero at extremality and no extremal black holes remain. However, when the extremal limit is supersymmetric, a large microscopic degeneracy survives and there is a gap to the first excited black hole visible from gravity. In this article we study large N quantum mechanical models where supersymmetry is explicitly broken, allowing us to interpolate between these two qualitatively different pictures. We propose and analyze deformations of N = 2 SYK models with such a pattern of (super) symmetry breaking which violate the U(1) R-symmetry. These models feature a lifting of the BPS degeneracy and a closing of the spectral gap, and we further show that the large N soft effective action is given by a modification of the N = 2 Schwarzian theory in which the U(1) R mode becomes massive.

1/N Expansion↗
DOE OSTI2024

DK/Dπ scattering and an exotic virtual bound state at the SU(3) flavour symmetric point from lattice QCD

Abstract ElasticS-wave scattering of a charm meson with a light pseudoscalar meson inJ P = 0 + is investigated in the flavour$$ \overline{\textbf{3}} $$ 3 ¯ ,6and$$ \overline{\textbf{15}} $$ 15 ¯ sectors at the SU(3) f flavour point using lattice QCD, working on three volumes withm π ≈ 700 MeV. Large bases of interpolating operators are employed to extract finite-volume spectra, which are subsequently used with the Lüscher method to provide constraints on infinite-volume scattering amplitudes. Examining the singularities of the amplitudes, theS-wave amplitude in the flavour$$ \overline{\textbf{3}} $$ 3 ¯ sector is found to contain a deeply bound state, strongly coupled to elastic threshold, corresponding to the$$ {J}^P={0}^{+}\ {D}_{s0}^{\ast }(2317) $$ J P = 0 + D s 0 ∗ 2317 . In the exotic flavour6sector a virtual bound state is found at$$ \sqrt{s_{\textrm{pole}}} $$ s pole = 2510 – 2610 MeV, roughly 40 – 140 MeV below threshold, whereas the$$ \overline{\textbf{15}} $$ 15 ¯ channel shows weak repulsion.

Physics↗