SEARCH · Search NASA

Results for “Gradient estimation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS2020

The bioclimatic extent and pattern of the cold edge of the boreal forest: the circumpolar taiga-tundra ecotone

Current configurations of forest structure at the cold edge of the boreal may help understand the future of ecosystem functioning in high northern latitudes. The circumpolar biome boundary at the boreal (taiga) forest and tundra interface is an ecological transition zone (taiga-tundra ecotone; TTE) experiencing changes that affect its forest structure. We accounted for the TTE's horizontal forest structure with an estimate of its extent and pattern as represented by tree canopy cover (TCC). We quantified TCC patterns with an algorithm that describes its spatial gradient, and summarized landscape patterns of structure to represent heterogeneity, capturing abrupt, diffuse, and uniform forest at mesoscales. We used these landscape patterns to constrain the spatial extent of sparse and open canopy forest, and non-forest (forest-adjacent) edge that defines the TTE extent. The resulting map of the TTE extent is based on forest structure spatial patterns resolved at 30 m, highlights structural variability across landscapes, and helps distinguish tundra from boreal domains. We classified 14 594 landscapes as those associated with the TTE within a circumpolar bioclimatic envelope (11.575 million km2), where 44.83% of the area of these landscapes were forest and non-forest edge, yet 36.43% contributed to the TTE extent. We report the overall extent of the TTE (3.032 million km2) across North America and Greenland (53%), and Eurasia (47%), where 0.697 million km2 is non-forest edge, 0.549 million km2 is sparse forest, and 1.787 million km2 is open canopy forest. Diffuse forest landscapes dominate the TTE (79%), and abrupt landscapes (~19%) indicate portions of the TTE where sparse forest and non-forest edge are the prevailing structural patterns. This account of the TTE quantifies the area of the cold edge of the boreal forest where previous global estimates show high discrepancies, and can help target monitoring and prediction of circumpolar dynamics.

Paul M Montesano↗

NASA NTRS1983

Seismic velocity structure of the lunar mantle

The recently completed set of seismic arrival times from the Apollo lunar seismic network are inverted to estimate the average seismic velocities in three sections of the lunar mantle: two for the upper mantle and one for the middle mantle. The method used is a variation of the linearized least squares inversion where the inversion is accomplished in steps. The estimated average velocities in the upper mantle decrease from Vp = 7.74 km/sec and Vs = 4.49 km/sec in the section above 270-km depth to Vp = 7.46 km/sec and Vs = 4.25 km/sec in the section between 270- and 500-km depth, confirming the earlier finding of negative gradients based on seismic amplitude variations. The average velocities in the middle mantle between the depths of 500 km and 1000 km of Vp = 8.26 km/sec and Vs = 4.65 km/sec are significantly higher than those in the upper mantle, contradicting earlier estimates based on more limited data. The higher velocities may suggest initial melting of the moon down to at least 1000-km depth.

Nakamura, Y.↗

DOE OSTI2025

Robustness of Deep Learning Classification to Adversarial Input on GPUs: Asynchronous Parallel Accumulation Is a Source of Vulnerability

The ability of machine learning (ML) classification models to resist small, targeted input perturbations—known as adversarial attacks—is a key measure of their safety and reliability. We show that floating-point non associativity (FPNA) coupled with asynchronous parallel programming on GPUs is sufficient to result in misclassification, without any perturbation to the input. Additionally, we show that this misclassification is particularly significant for inputs close to the decision boundary and that standard adversarial robustness results may be overestimated up to 4.6 when not considering machine-level details. We first study a linear classifier, before focusing on standard Graph Neural Network (GNN) architectures and datasets used in robustness assessments. We develop a novel black-box attack using Bayesian optimization to discover external workloads that can change the instruction scheduling which bias the output of reductions on GPUs and reliably lead to misclassification. Motivated by these results, we present a new learnable permutation (LP) gradient-based approach to learning floating-point operation orderings that lead to misclassifications. The LP approach provides a worst-case estimate in a computationally efficient manner, avoiding the need to run identical experiments tens of thousands of times over a potentially large set of possible GPU states or architectures. Finally, using instrumentation-based testing, we investigate parallel reduction ordering across different GPU architectures under external background workloads, when utilizing multi-GPU virtualization, and when applying power capping. Our results demonstrate that parallel reduction ordering varies significantly across architectures under the first two conditions, substantially increasing the search space required to fully test the effects of this parallel scheduler-based vulnerability. These results and the methods developed here can help to include machine-level considerations into adversarial robustness assessments, which can make a difference in safety and mission critical applications.

Shanmugavelu, Sanjif [Maxeler Technologies, a Groq↗

DOE OSTI2024

Concurrent Measurement of O 2 Production and Isoprene Emission During Photosynthesis: Pros, Cons and Metabolic Implications of Responses to Light, CO 2 and Temperature

Traditional leaf gas exchange experiments have focused on net CO 2 exchange (A net ). Here, using California poplar (Populus trichocarpa), we coupled measurements of net oxygen production (NOP), isoprene emissions and δ 18 O in O 2 to traditional CO 2 /H 2 O gas exchange with chlorophyll fluorescence, and measured light, CO 2 and temperature response curves. This allowed us to obtain a comprehensive picture of the photosynthetic redox budget including electron transport rate (ETR) and estimates of the mean assimilatory quotient (AQ = A net /NOP). We found that A net and NOP were linearly correlated across environmental gradients with similar observed AQ values during light (1.25 ± 0.05) and CO 2 responses (1.23 ± 0.07). In contrast, AQ was suppressed during leaf temperature responses in the light (0.87 ± 0.28), potentially due to the acceleration of alternative ETR sinks like lipid synthesis. A net and NOP had an optimum temperature (Topt) of 31°C, while ETR and δ 18 O in O2 (35°C) and isoprene emissions (39°C) had distinctly higher T opt . The results confirm a tight connection between water oxidation and ETR and support a view of light-dependent lipid synthesis primarily driven by photosynthetic ATP/NADPH not consumed by the Calvin–Benson cycle, as an important thermotolerance mechanism linked with high rates of (photo)respiration and CO 2 /O 2 recycling.

H218O labelling↗

NASA NTRS1981

Pinhole-coronograph tracking control

The pinhole occulter system, a space shuttle based experiment for the production of hard X-ray images taken primarily from the Sun is described. The system is a pinhole camera utilizing a deployable 50 m flexible boom for separating the pinhole from the recording devices located in the shuttle. The system is subjected to the disturbances of the shuttle. The worst of these are thruster firing for orbit correction, man motion, motion induced by other systems, and gravity gradient torques. The control system of the pointing mount senses both position and velocity of the mask tip and uses these to estimate accurately the flexible modes of the system.

Greene, M.↗

NASA NTRS1986

Ray tracing survey of Z mode emissions from source regions in the high-altitude auroral zone

The ray tracing computer program described by Menietti and Lin (1985) was used to investigate the propagation of Z mode emissions from sources in the auroral zone. The calculation results indicate that the propagation is primarily perpendicular to the magnetic field line (horizontal propagation) for almost all sources considered, except when the wave source is assumed to be unrestricted by the cyclotron maser mechanism. While the propagation is initially downward, refraction rapidly increases the wave normal angle, and only propagation at near-constant altitude occurs. The frequency bandwidth estimated from the ray paths was found to be less than that observed by the DE 1 satellite. Large horizontal magnetic field gradients associated with field-aligned currents produce downward propagation but are still not sufficient to produce the observed bandwidths. The results suggest that wave growth may be due to mechanisms in addition to cyclotron maser resonance.

Menietti, J. D.↗

NASA NTRS1988

Non-Maxwellian distribution functions in flaring coronal loops - Comparison of Landau-Fokker-Planck and BGK solutions

The high-velocity tail of the electron distribution has been calculated by solving the high-velocity form of the Landau equation for a thermal structure representative of a flaring coronal loop. These calculations show an enhancement of the tail population above Maxwellian for electrons moving down the temperature gradient. The results obtained are used to test the reliability of the BGK approximation. The comparison shows that the BGK technique can estimate contributions to the heat flux from the high-energy tail to within an order of magnitude.

Ljepojevic, N. N.↗

NASA NTRS1996

The Character of the Beaufort Sea Ice Cover From Active and Passive Microwave Observations

By definition, ice which survives the summer is classified as multiyear ice. Thus, multiyear ice concentration during the winter should be nearly equivalent to the ice concentration during the previous summer's minima. This condition provides a reasonable criterion for evaluation of ice concentration and ice type retrieval algorithms from remote sensing datasets. From SSM/I data, the NASA Team algorithm estimates the multiyear, first-year and total ice concentrations throughout the sea ice season using combinations of the polarization and spectral gradient ratios.

Beaufort Sea↗

NASA NTRS2012

An Empirical Temperature Variance Source Model in Heated Jets

An acoustic analogy approach is implemented that models the sources of jet noise in heated jets. The equivalent sources of turbulent mixing noise are recognized as the differences between the fluctuating and Favre-averaged Reynolds stresses and enthalpy fluxes. While in a conventional acoustic analogy only Reynolds stress components are scrutinized for their noise generation properties, it is now accepted that a comprehensive source model should include the additional entropy source term. Following Goldstein s generalized acoustic analogy, the set of Euler equations are divided into two sets of equations that govern a non-radiating base flow plus its residual components. When the base flow is considered as a locally parallel mean flow, the residual equations may be rearranged to form an inhomogeneous third-order wave equation. A general solution is written subsequently using a Green s function method while all non-linear terms are treated as the equivalent sources of aerodynamic sound and are modeled accordingly. In a previous study, a specialized Reynolds-averaged Navier-Stokes (RANS) solver was implemented to compute the variance of thermal fluctuations that determine the enthalpy flux source strength. The main objective here is to present an empirical model capable of providing a reasonable estimate of the stagnation temperature variance in a jet. Such a model is parameterized as a function of the mean stagnation temperature gradient in the jet, and is evaluated using commonly available RANS solvers. The ensuing thermal source distribution is compared with measurements as well as computational result from a dedicated RANS solver that employs an enthalpy variance and dissipation rate model. Turbulent mixing noise predictions are presented for a wide range of jet temperature ratios from 1.0 to 3.20.

Khavaran, Abbas↗

NASA NTRS1996

(abstract) Three Dimensional Ice-Flow Velocity Estimation Using Satellite Radar Interferometry

An understanding of the flow dynamics of an ice sheet's outlet glaciers and ice streams requires knowledge of their flow velocity and strain rates (i.e., velocity gradients). Prior to the recent advent of satellite radar interferometry, it was not possible to measure detailed ice-flow velocity over the vast featureless areas that comprise most of the ice sheets. Since the launch of ERS-1, the use of satellite radar interferometry data for making densly sampled ice-flow velocity measurements has been firmly established by several studies. We have combined data from nonparallel orbits with surface slope information to make vector ice-flow measurements for the Ryder Glacier, Greenland. Our results for the Ryder are promising and indicate that repeat-pass interferometric data can be used to make vector measurements of ice velocity.

ice sheets radar interferometry glaciers flow dyna↗

NASA NTRS2009

Conical-Domain Model for Estimating GPS Ionospheric Delays

The conical-domain model is a computational model, now undergoing development, for estimating ionospheric delays of Global Positioning System (GPS) signals. Relative to the standard ionospheric delay model described below, the conical-domain model offers improved accuracy. In the absence of selective availability, the ionosphere is the largest source of error for single-frequency users of GPS. Because ionospheric signal delays contribute to errors in GPS position and time measurements, satellite-based augmentation systems (SBASs) have been designed to estimate these delays and broadcast corrections. Several national and international SBASs are currently in various stages of development to enhance the integrity and accuracy of GPS measurements for airline navigation. In the Wide Area Augmentation System (WAAS) of the United States, slant ionospheric delay errors and confidence bounds are derived from estimates of vertical ionospheric delay modeled on a grid at regularly spaced intervals of latitude and longitude. The estimate of vertical delay at each ionospheric grid point (IGP) is calculated from a planar fit of neighboring slant delay measurements, projected to vertical using a standard, thin-shell model of the ionosphere. Interpolation on the WAAS grid enables estimation of the vertical delay at the ionospheric pierce point (IPP) corresponding to any arbitrary measurement of a user. (The IPP of a given user s measurement is the point where the GPS signal ray path intersects a reference ionospheric height.) The product of the interpolated value and the user s thin-shell obliquity factor provides an estimate of the user s ionospheric slant delay. Two types of error that restrict the accuracy of the thin-shell model are absent in the conical domain model: (1) error due to the implicit assumption that the electron density is independent of the azimuthal angle at the IPP and (2) error arising from the slant-to-vertical conversion. At low latitudes or at mid-latitudes under disturbed conditions, the accuracy of SBAS systems based upon the thin-shell model suffers due to the presence of complex ionospheric structure, high delay values, and large electron density gradients. Interpolation on the vertical delay grid serves as an additional source of delay error. The conical-domain model permits direct computation of the user s slant delay estimate without the intervening use of a vertical delay grid. The key is to restrict each fit of GPS measurements to a spatial domain encompassing signals from only one satellite. The conical domain model is so named because each fit involves a group of GPS receivers that all receive signals from the same GPS satellite (see figure); the receiver and satellite positions define a cone, the satellite position being the vertex. A user within a given cone evaluates the delay to the satellite directly, using (1) the IPP coordinates of the line of sight to the satellite and (2) broadcast fit parameters associated with the cone. The conical-domain model partly resembles the thin-shell model in that both models reduce an inherently four-dimensional problem to two dimensions. However, unlike the thin-shell model, the conical domain model does not involve any potentially erroneous simplifying assumptions about the structure of the ionosphere. In the conical domain model, the initially four-dimensional problem becomes truly two-dimensional in the sense that once a satellite location has been specified, any signal path emanating from a satellite can be identified by only two coordinates; for example, the IPP coordinates. As a consequence, a user s slant-delay estimate converges to the correct value in the limit that the receivers converge to the user s location (or, equivalently, in the limit that the measurement IPPs converge to the user s IPP).

Sparks, Lawrence↗

NASA NTRS1986

History matching by spline approximation and regularization in single-phase areal reservoirs

An automatic history matching algorithm is developed based on bi-cubic spline approximations of permeability and porosity distributions and on the theory of regularization to estimate permeability or porosity in a single-phase, two-dimensional real reservoir from well pressure data. The regularization feature of the algorithm is used to convert the ill-posed history matching problem into a well-posed problem. The algorithm employs the conjugate gradient method as its core minimization method. A number of numerical experiments are carried out to evaluate the performance of the algorithm. Comparisons with conventional (non-regularized) automatic history matching algorithms indicate the superiority of the new algorithm with respect to the parameter estimates obtained. A quasioptimal regularization parameter is determined without requiring a priori information on the statistical properties of the observations.

Lee, T. Y.↗

DOE OSTI2024

Measure this, not that: Optimizing the cost and model-based information content of measurements

Model-based design of experiments (MBDoE) is a powerful framework for selecting and calibrating science-based mathematical models from data. Here, this work extends popular MBDoE workflows by proposing a convex mixed integer (non)linear programming (MINLP) to optimize the selection of measurements. The solver MindtPy is modified to support calculating the D-optimality objective and its gradient via an external package, scipy, using the grey-box module in Pyomo. The new approach is demonstrated in two case studies: estimating highly correlated kinetics from a batch reactor and estimating transport parameters in a large-scale rotary packed bed for CO 2 capture. Both case studies show how examining the Pareto optimal trade-offs between information content measured by A- and D-optimality versus measurement budget offers practical guidance for selecting measurements for scientific experiments.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Mapping strain and structural heterogeneities around bubbles in amorphous ionically conductive Bi 2 O 3

While amorphous materials are often approximated to have a statistically homogeneous atomic structure, they frequently exhibit localized structural heterogeneity that challenges simplified models. This study uses 4D scanning transmission electron microscopy to investigate the strain and structural modifications around gas bubbles in amorphous Bi 2 O 3 induced by argon irradiation. We present a method for determining strain fields surrounding bubbles that can be used to measure the internal pressure of the gas. Compressive strain is observed around the cavities, with higher-order crystalline symmetries emerging near the cavity interfaces, suggesting paracrystalline ordering as a result of bubble coarsening. This ordering, along with a compressive strain gradient, indicates that gas bubbles induce significant localized changes in atomic packing. By analyzing strain fields with maximum compressive strains of 3%, we estimate a lower bound on the internal pressure of the bubbles at 2.5 GPa. These findings provide insight into the complex structural behavior of amorphous materials under stress, particularly in systems with gas inclusions, and offer new methods for probing the local atomic structure in disordered materials. Although considering structural heterogeneity in amorphous systems is non-trivial, these features have crucial impacts on material functionalities, such as mechanical strength, ionic conductivity, and electronic mobility.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Interpretation of high-harmonic fast-wave propagation in the scrape-off layer of NSTX-U as a geometrically bounded-waveguide mode

This paper revisits the onset of high-harmonic fast-wave (HHFW) propagation in the scrape-off layer (SOL) plasma of the NSTX/NSTX-U spherical tokamak, motivated by past HHFW heating and current drive experiments and modelling. Previously, the fast-wave propagation in the SOL was correlated with the opening (suppression) of the fast-wave right-hand cutoff layer in front of the antenna. In this work, the SOL propagation is interpreted as a geometric waveguide mode guided by the SOL geometry, whose radial width can be comparable to the wave’s perpendicular wavelength in the HHFW regime. A two-dimensional circular model is first employed using full-wave solvers to characterise the poloidal eigenmode structures and to clarify their relationship to annulus resonance. By progressively adding a tokamak-like magnetic-field configuration, starting from a uniform axial field, the effects of magnetic-field gradients and pitch on the SOL waveguide mode are characterised. The poloidal mode numbers supported by the SOL plasma agree well with analytic estimates, indicating that an anisotropic plasma in a bounded geometry selectively supports and amplifies the resonant poloidal mode number. Additionally, two-dimensional axisymmetric NSTX-U simulations demonstrate that the SOL eigenmode features identified in the circular model persist in the experimentally relevant configurations. A control approach based on lengthening the wave perpendicular wavelength relative to the SOL width, including a higher magnetic (B)-field operation, is discussed. The analysis here shows the key role of a bounded geometry in interpreting HHFW eigenmode coupling and propagation in the NSTX-U SOL plasma.

plasma heating↗

DOE OSTI2025

Influence of Soil Moisture on the Development of Organized Convective Systems in South America

Abstract This study investigates the role of soil moisture (SM) on the initiation and organization of convective systems using the convection‐permitting ICOsahedral Non‐hydrostatic (ICON) model. We conduct two sets of experiments: a Control experiment with interactive SM and a fixed SM experiment (FixedSM) with invariable SM conditions. We focus on two regions in South America: the Amazon and southeastern South America (SESA). Larger organized convective systems are associated with greater SM heterogeneity in both regions, though other large‐scale synoptic influences affect the robustness of this relationship in SESA. These results remain largely unaffected by disabling the effects of precipitation on SM in the FixedSM experiment, and complementary analyses using satellite‐based estimates of SM and precipitation support these findings. Spatial compositing of mesoscale environments in the Amazon shows the presence of well‐defined SM gradients, at a length scale of a few hundred kilometers, many hours before convective system detection. Larger SM gradients correspond to larger gradients in thermodynamic variables, particularly surface temperature and sensible heat flux, and are associated with larger convective systems. Overall, our findings suggest that surface heterogeneities such as SM gradients not only affect deep convection initiation, as previously suggested, but they can also encourage the growth and organization of convective systems into larger clusters, particularly in the absence of significant synoptic influences.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Lattice B -field correlators for heavy quarks

We analyze the color-magnetic (or “ B ”) field two-point function that encodes the finite-mass correction to the heavy-quark momentum-diffusion coefficient. The simulations are done on fine isotropic lattices in the quenched approximation at 1.5 T c , using a range of gradient flow times for noise suppression and operator renormalization. The continuum extrapolation is performed at fixed flow time followed by a second extrapolation to zero flow time. Perturbative calculations to next-to-leading order of this correlation function, matching gradient-flowed correlators to MS ¯ , are used to resolve nontrivial renormalization issues. We perform a spectral reconstruction based on perturbative model fits to estimate the coefficient κ B of the finite-mass correction to the heavy-quark momentum-diffusion coefficient. The approach we present here yields high-precision data for the correlator with all renormalization issues incorporated at next-to-leading order and is also applicable for actions with dynamical fermions. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1979

Covariance analyses of satellite-derived mesoscale wind fields

Statistical structure functions have been computed independently for nine satellite-derived mesoscale wind fields that were obtained on two different days. Small cumulus clouds were tracked at 5 min intervals, but since these clouds occurred primarily in the warm sectors of midlatitude cyclones the results cannot be considered representative of the circulations within cyclones in general. The field structure varied considerably with time and was especially affected if mesoscale features were observed. The wind fields on the 2 days studied were highly anisotropic with large gradients in structure occurring approximately normal to the mean flow. Structure function calculations for the combined set of satellite winds were used to estimate random error present in the fields. It is concluded for these data that the random error in vector winds derived from cumulus cloud tracking using high-frequency satellite data is less than 1.75 m/s. Spatial correlation functions were also computed for the nine data sets. Normalized correlation functions were considerably different for u and v components and decreased rapidly as data point separation increased for both components. The correlation functions for transverse and longitudinal components decreased less rapidly as data point separation increased.

Maddox, R. A.↗