Search NASASearch

SEARCH · Search NASA

Results for “Directional derivatives”

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 19 records

Augmenting subspace optimization methods with linear bandits

In this work, we consider the framework of methods for unconstrained minimization that are, in each iteration, restricted to a model that is only a valid approximation to the objective function on some affine subspace containing an incumbent point. These methods are of practical interest in computational settings where derivative information is either expensive or impossible to obtain. Recent attention has been paid in the literature to employing randomized matrix sketching for generating the affine subspaces within this framework. We consider a relatively straightforward, deterministic augmentation of such a generic subspace optimization method. In particular, we consider a sequential optimization framework where actions consist of one-dimensional linear subspaces and rewards consist of (approximations to) the magnitudes of directional derivatives computed in the direction of the action subspace. Reward maximization in this context is consistent with maximizing lower bounds on descent guaranteed by first-order Taylor models. This sequential optimization problem can be analysed through the lens of dynamic regret. We modify an existing linear upper confidence bound (UCB) bandit method and prove sublinear dynamic regret in the subspace optimization setting. We demonstrate the efficacy of employing this linear UCB method in a setting where forward-mode algorithmic differentiation can provide directional derivatives in arbitrary directions and in a derivative-free setting. For the derivative-free setting, we propose SS-POUNDers, an extension of the derivative-free optimization method POUNDers that employs the linear UCB mechanism to identify promising subspaces. Our numerical experiments suggest a preference, in either computational setting, for employing a linear UCB mechanism within a subspace optimization method.

97 MATHEMATICS AND COMPUTING

An Accurate SUPG-stabilized Continuous Galerkin Discretization for Anisotropic Heat Flux in Magnetic Confinement Fusion

We present a novel spatial discretization for the anisotropic heat conduction equation, aimed at improved accuracy at the high levels of anisotropy seen in a magnetized plasma, for example, for magnetic confinement fusion. The new discretization is based on a mixed formulation, introducing a form of the directional derivative along the magnetic field as an auxiliary variable and discretizing both the temperature and auxiliary fields in a continuous Galerkin (CG) space. Both the temperature and auxiliary variable equations are stabilized using the streamline upwind Petrov–Galerkin (SUPG) method, ensuring a better representation of the directional derivatives and therefore an overall more accurate solution. This approach can be seen as the CG-based version of our previous work (Wimmer, Southworth, Gregory, Tang, 2024), where we considered a mixed discontinuous Galerkin (DG) spatial discretization including DG-upwind stabilization. We prove consistency of the novel discretization, and demonstrate its improved accuracy over existing CG-based methods in test cases relevant to magnetic confinement fusion. This includes a long-run tokamak equilibrium sustainment scenario, demonstrating a 35% and 32% spurious heat loss for existing primal and mixed CG-based formulations versus 4% for our novel SUPG-stabilized discretization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St

Pareto-optimal target definition for multi-axis random vibration testing

In random vibration testing with multiple control channels, existing control laws require specification of a complete spectral density matrix at each control frequency. Spectral density matrices include autospectral densities on the diagonal and cross-spectral densities on the off-diagonal. In practice, the off-diagonal terms are often unknown, and recent vibration testing research has focused on fixing the diagonal and specifying the off-diagonal to minimize the required control energy, subject to a constraint that the target matrix is positive semidefinite. This paper shows that, even with a fixed diagonal, off-diagonal terms strongly affect control residuals. This overlooked effect occurs in both square and rectangular systems. By jointly considering input energy and control residuals, open-loop inputs are derived directly from the diagonal without specifying the off-diagonal terms. Vibration targets that can be used in closed-loop control are then derived using the optimal inputs, with positive semidefinite constraints applied during the derivation. The result is a set of Pareto-optimal control solutions. For each solution in the set, any other possible solution produces greater control error, greater input energy, or both. A balanced solution is selected automatically, though others can be chosen based on test needs. Simulations and experiments show that the proposed method outperforms state-of-the-art energy-minimizing approaches, achieving significant reductions in both control error and input energy.

Autospectral density

Reconstructing Quasar Spectra and Measuring the Lyα Forest with SpenderQ

Quasar spectra carry the imprint of foreground intergalactic medium (IGM) through absorption features. In particular, absorption caused by neutral hydrogen gas, the "Lyα forest," is a key spectroscopic tracer for cosmological analyses used to measure cosmic expansion and test physics beyond the standard model. Despite their importance, current methods for measuring Lyα absorption cannot directly derive the intrinsic quasar continuum and make strong assumptions on its shape, thus distorting the measured Lyα clustering. We present SpenderQ , a ML-based approach for directly reconstructing the intrinsic quasar spectra and measuring the Lyα forest from observations. SpenderQ uses the Spender spectrum autoencoder to learn a compact and redshift-invariant latent encoding of quasar spectra, combined with an iterative procedure to identify and mask absorption regions. To demonstrate its performance, we apply SpenderQ to 400,000 synthetic quasar spectra created to validate the Dark Energy Spectroscopic Instrument Year 1 Lyα cosmological analyses. SpenderQ accurately reconstructs the true intrinsic quasar spectra, including the broad Lyβ, Lyα, SiIV, CIV, and CIII emission lines. Redward of Lyα, SpenderQ provides percent-level reconstructions of the true quasar spectra. Blueward of Lyα, SpenderQ reconstructs the true spectra to < 5%. SpenderQ reproduces the shapes of individual quasar spectra more robustly than the current state-of-the-art. We, thus, expect it will significantly reduce biases in Lyα clustering measurements and enable studies of quasars and their physical properties. SpenderQ also provides informative latent variable encodings that can be used to, e.g., classify quasars with Broad Absorption Lines. Overall, SpenderQ provides a new data-driven approach for unbiased Lyα forest measurements in cosmological, quasar, and IGM studies.

Hahn, ChangHoon [Arizona U., Astron. Dept. - Stewa

Dual C–H functionalization of polyolefins for covalent adaptable network formation via cooperative electrolysis

The increasing demand for carbon fibre-reinforced polymer composites with polyolefin thermoset matrices will lead to the generation of large volumes of low molecular weight waste oligomers during the fibre recovery process. To address this challenge, we present a cooperative electrolytic dual C–H bond functionalization strategy for the one-step installation of two key functional groups essential for dynamic linkages of these deconstructed oligomers. Through direct competition studies, we reveal the predominance of tertiary allylic C–H activation along the complex, branched oligomer backbone. The resulting difunctionalized oligomers form covalently adaptable networks with exceptional circularity. By establishing a sustainable pathway for upcycling deconstructed oligomers from carbon fibre-reinforced polymer fibre recovery, this strategy introduces thermoset materials derived directly from waste. Here, it advances sustainable manufacturing practices, contributing to a net-zero waste synthetic ecosystem, and highlights the prospects of mediated electrolysis in macromolecular transformations.

Electrocatalysis

Experimental determination of energy-level structures of diamond silicon-vacancy centers in off-axis magnetic fields

Here, we report the development of an experimental approach to efficiently determine the energy level structure of an individual silicon vacancy (SiV) center in a magnetic field along an arbitrary direction. This approach uses two coupling rates (one each for the ground and the excited states) to characterize the combined effects of static strain and dynamic Jahn-Teller coupling and exploits the fact that orbital Zeeman effects vanish when the magnetic field is normal to the SiV axis. With an analytical expression for the energy level structure of the SiV under a transverse magnetic field, the two coupling rates can be directly derived from two measurements: one on the frequency separation between two spin-conserved transitions and the other on the coherent population trapping resonance of the SiV ground spin states. A detailed comparison between the numerical calculation and the experimental result on the dependence of the spin-conserved splitting on both the amplitude and direction of the magnetic field further reveals unequal orbital magnetic coupling for the ground and excited states, indicating that unequal orbital quenching factors are needed for an accurate description of the SiV energy level structure in a magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

CUDO: closed-form universal dwell-time optimization for computer-controlled optical surfacing

Precision optical figuring demands fast and accurate dwell time optimization to reach nanometer- and sub-nanometer-level accuracy in next-generation optical systems. We introduce CUDO (closed-form universal dwell-time optimization), the first, to the best of our knowledge, unified closed-form analytical framework that supports both function-form and matrix-form dwell time models in computer-controlled optical surfacing (CCOS). In contrast to traditional methods, which rely on iterative optimization and hyperparameter tuning, our framework derives direct analytical solutions with no adjustable parameters. This approach unifies the solution principles of existing methods within a single mathematical model, delivering three key advantages: (1) accuracy on par with, or superior to, iterative solvers, (2) substantial reduction in computation time, and (3) numerical robustness. Comparative studies with prior art confirm that closed-form solutions achieve equivalent residual error while removing runtime bottlenecks. By simplifying the implementation and enabling real-time, scalable deployment, CUDO establishes a practical foundation for future deterministic fabrication of large-aperture and high-performance optics.

36 MATERIALS SCIENCE

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Derivation and verification of the direct-sampling method for simulating Monte Carlo flight paths in tetrahedral meshes with linear finite-element cross sections

This paper provides a derivation of a direct-sampling approach for modeling continuously varying cross sections in tetrahedral-mesh-based Monte Carlo codes. Specifically, cross sections are spatially approximated using linear nodal finite elements. A linearization strategy is provided for non-linearly varying cross sections. The method is verified against seven analytical pure-absorber test problems. These test problems also highlight the benefit of using linear finite elements over element-wise-constant cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

U-net architected deep material network training with microstructure local field information

The Deep Material Network (DMN) has recently emerged as a powerful reduced-order modeling framework for simulating the mechanical response of heterogeneous materials such as composites. Unlike most data-driven approaches that directly learn a material’s response under prescribed loading, the DMN acts as a homogenization operator, learning the kinematic constraints and mechanical interactions of the underlying microstructure. However, traditional DMN training relies exclusively on homogenized effective properties derived from Direct Numerical Simulations (DNS), discarding the rich local field data that govern microstructural interactions. In this work, we extend the DMN framework to incorporate such local field information into the offline training process. Utilizing a U-Net architecture, we augment the DMN training objective to include the first and second statistical moments of the local stress fields obtained from linear DNS. This ensures that the learned network topology not only fits the effective stiffness but also accurately reflects the internal local stress and strain partitioning of the microstructure. The results confirm that supervising the localization process during training yields a superior surrogate model, reducing local prediction errors by an order of magnitude and significantly improving generalization to unseen nonlinear constitutive behaviors compared to traditional DMNs.

36 MATERIALS SCIENCE

A Review of Direct Ink Writing of Polymer Derived Ceramics

With the growing demand for materials capable of withstanding extreme temperatures and pressures, ceramic components with exceptional corrosion resistance and reliable mechanical properties have experienced a significant surge in demand. However, traditional ceramic forming methods involve high-temperatures and energy-intensive processes that often struggle to produce complex parts or composites efficiently. Polymer-Derived Ceramics (PDCs) offer a transformative solution by using polymeric precursors that can be converted into a wide variety of silicon-based and non-silicon-based ceramics through heat treatment. The polymeric nature of PDC precursors enables the fabrication of geometrically intricate components using conventional polymer-forming techniques at significantly lower processing temperatures. Furthermore, PDCs are well-suited for additive manufacturing (AM), allowing the production of complex structural and functional components through cost-effective, low-temperature processes. By leveraging the diverse properties of PDC materials – each with unique advantages and limitations – manufacturers can optimise performance for specific applications. This review provides an overview of the types of PDCs developed to date and their broad range of applications. Specifically, it delves into the Direct Ink Writing (DIW) process, exploring its rheological requirements and the critical role of fillers in tailoring the rheological properties of polymeric precursors to meet the specific demands of DIW.

Victoria Bishop

Energy shift of Fe-K fluorescence lines due to low ionization demonstrated with XRISM in Centaurus X-3

The Fe K$\alpha$ fluorescence line at 6.4 keV is a powerful probe of cold matter surrounding X-ray sources and has been widely used in various astrophysical contexts. The X-ray microcalorimeter spectrometer onboard XRISM can measure line shifts with unprecedented precision of ${\sim }$0.2 eV, equivalent to a line-of-sight velocity of ${\sim }$10 km s$^{-1}$. At this level of accuracy, however, several factors that influence the line energy must be carefully considered prior to astrophysical interpretation. One such important factor is the ionization degree, Fe$^{q+}$. The K$\alpha$ line shifts redward by ${\sim }$4 eV as q increases from 0 (neutral) to 8 (Ar-like). Additionally, the accompanying Fe K$\beta$ line at 7.06 keV shifts blueward by ${\sim }$30 eV from $q=0$ to 8. We demonstrate that this effect is actually observable in the XRISM data of the high-mass X-ray binary Centaurus X-3 (Cen X-3). We advocate that the differential energy shift between the K$\alpha$ and K$\beta$ line provides a robust estimate of q by decoupling from other effects that shift the two lines in the same direction. We derived $q \sim 5$ (Sc-like) for the fluorescing matter by comparing the observation with atomic structure calculations of our own and in the literature. By accounting for the derived charge state and the corresponding shift in the rest-frame line energy, we made corrections for this effect and reached a consistent residual shift among the K$\alpha$, K$\beta$, and the optical measurement attributable to the systemic velocity of the system. Consequently, we obtained a new constraint on the location of the cold matter. This ionization effect needs to be assessed in all use cases of the Fe K$\alpha$ line shift beyond Cen X-3, and the proposed metric is generally applicable to all of them.

Astronomy and AstroPhysics

SPRUCE Ground Observations of Phenology in Experimental Plots, 2024

This data set consists of one comma separated (*.csv) file containing phenological transition dates, as derived from direct observations of vegetative and reproductive phenology recorded by a human observer, from the SPRUCE experiment during 2024 (2025-03-06 to 2025-11-21), the ninth full year of whole-ecosystem warming (Hanson et al. 2017). Both spring and autumn phenological events are included. Since April 2016, human observers have been directly tracking the phenology of both woody and herbaceous species on a weekly schedule within the SPRUCE experimental chambers, these data are reported in annual ground observations data sets (see Related Data Sets). The observed date reported here is the first survey date in 2024 on which an event/phenophase was definitively observed. This data set also contains a companion file in HTML (*.html) containing figures showing the relationship between the day of year and temperature treatment for different phenological phases by species for 2024.

54 ENVIRONMENTAL SCIENCES

Direct Synthesis of Polyester from Biomass Derived 1,6-hexanediol using a Copper-Ceria Catalyst

Biomass derived polyesters are gaining attention due to their increased mechanical and thermal properties, and also because of their biodegradability. These polyesters are typically synthesized from polycondensation reactions between diols with dicarboxylic acids using catalysts. Other methods such as alcoholysis and acidolysis of low molecular weight esters and alcoholysis of acyl chlorides are also well known. Direct conversion of diol to ester without using dicarboxylic acid is not well known, and to our knowledge has not been reported over a solid catalyst. In this study, we report a direct synthesis of polyester from 1,6-hexanediol using a copper-ceria catalyst without using a dicarboxylic acid monomer. The copper-ceria catalysts were synthesized by incipient wetness impregnation method and characterized using X-ray photoelectron spectroscopy (XPS), X-ray diffraction analysis (XRD), Brunauer-Emmett-Teller (BET), and infrared spectroscopy (IR). Our initial studies show that 1,6-hexanediol can be converted to corresponding polyester at higher temperatures under atmospheric pressures of CO2 or N2. The reaction products were analyzed using nuclear magnetic resonance spectroscopy (NMR), gel permeation chromatography (GPC), high performance liquid chromatography (HPLC), and IR spectroscopy.

Biomass derived diols

JetGP: A derivative enhanced Gaussian process library

Derivative enhanced Gaussian Processes (DEGPs) can significantly improve surrogate model accuracy over standard Gaussian Process (GP) formulations by incorporating derivative information. However, standard implementations scale poorly with dimension, limiting their use in high dimensional engineering problems. JetGP is a Python framework that unifies existing derivative enhanced GP methodologies into a single library and extends them to support arbitrary order derivative information. The library implements four complementary formulations: standard derivative enhanced Gaussian Processes (DEGP), directional DEGP (DDEGP), generalized directional DEGP (GDDEGP), and weighted DEGP (WDEGP). By unifying these approaches in a consistent interface with robust numerical implementations, JetGP enables practitioners to balance predictive accuracy and computational efficiency for high dimensional optimization, uncertainty quantification, and sensitivity analysis in engineering design.

Derivative enhanced Gaussian process

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

Excited-State Densities from Time-Dependent Density Functional Response Theory

While the variational principle for excited-state energies leads to a route to obtaining excited-state densities from time-dependent density functional theory, relatively little attention has been paid to the quality of the resulting densities in real space obtained with different exchange-correlation functional approximations or how nonadiabatic approximations developed for energies of states of double-excitation character perform for their densities. Here we derive an expression directly in real space for the excited-state density, which includes the case of nonadiabatic kernels and consequently is able, for the first time, to yield densities of states of double-excitation character. Under some well-defined simplifications, we compare the performance of the local-density approximation and exact-exchange approximation, which are in a sense at the opposite extremes of the fundamental functional approximations, on local and charge-transfer excitations in one-dimensional model systems and show that the dressed Time-Dependent Density Functional Theory (TDDFT) approach gives good densities of double excitations.

approximation