Search NASA⌕ Search

SEARCH · Search NASA

Results for “sensitivity gradient”

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 55 records · Page 3

Post-optimality analysis in aerospace vehicle design

This analysis pertains to the applicability of optimal sensitivity information to aerospace vehicle design. The present analysis demonstrates that post-optimality information generated through first-order computations can be used to accurately predict file effect of constraint and parameter perturbations on the optimal solution. This assessment is based on the solution of an aircraft design problem in which the post-optimality estimates are shown to be within a few percent of the true solution over the practical range of constraint and parameter variations. Through solution of a reusable, single-stage-to-orbit, launch vehicle design problem, this optimal sensitivity information is also shown to improve the efficiency of the design process. For a hierarchically decomposed problem, this computational efficiency is realizable by estimating the main-problem objective gradient through optimal sensitivity calculations. By reducing the need for finite differentiation of a re-optimized subproblem, a significant decrease in the number of objective function evaluations required to reach the optimal solution is obtained.

Braun, Robert D.↗

Application of Adjoint Methodology in Various Aspects of Sonic Boom Design

One of the advances in computational design has been the development of adjoint methods allowing efficient calculation of sensitivities in gradient-based shape optimization. This paper discusses two new applications of adjoint methodology that have been developed to aid in sonic boom mitigation exercises. In the first, equivalent area targets are generated using adjoint sensitivities of selected boom metrics. These targets may then be used to drive the vehicle shape during optimization. The second application is the computation of adjoint sensitivities of boom metrics on the ground with respect to parameters such as flight conditions, propagation sampling rate, and selected inputs to the propagation algorithms. These sensitivities enable the designer to make more informed selections of flight conditions at which the chosen cost functionals are less sensitive.

Rallabhandi, Sriram K.↗

Dissipation and Bathymetric Sensitivities in an Unstructured Mesh Global Tidal Model

Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M 2 tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.

54 ENVIRONMENTAL SCIENCES↗

Large carbon cycle sensitivities to climate across a permafrost thaw gradient in subarctic Sweden

Permafrost peatlands store large amounts of carbon potentially vulnerable to decomposition. However, the fate of that carbon in a changing climate remains uncertain in models due to complex interactions among hydrological, biogeochemical, microbial, and plant processes. In this study, we estimated effects of climate forcing biases present in global climate reanalysis products on carbon cycle predictions at a thawing permafrost peatland in subarctic Sweden. The analysis was conducted with a comprehensive biogeochemical model ( ecosys ) across a permafrost thaw gradient encompassing intact permafrost palsa with an ice core and a shallow active layer, partly thawed bog with a deeper active layer and a variable water table, and fen with a water table close to the surface, each with distinct vegetation and microbiota. Using in situ observations to correct local cold and wet biases found in the Global Soil Wetness Project Phase 3 (GSWP3) climate reanalysis forcing, we demonstrate good model performance by comparing predicted and observed carbon dioxide (CO 2 ) and methane (CH 4 ) exchanges, thaw depth, and water table depth. The simulations driven by the bias-corrected climate suggest that the three peatland types currently accumulate carbon from the atmosphere, although the bog and fen sites can have annual positive radiative forcing impacts due to their higher CH 4 emissions. Our simulations indicate that projected precipitation increases could accelerate CH 4 emissions from the palsa area, even without further degradation of palsa permafrost. The GSWP3 cold and wet biases for this site significantly alter simulation results and lead to erroneous active layer depth (ALD) and carbon budget estimates. Biases in simulated CO 2 and CH 4 exchanges from biased climate forcing are as large as those among the thaw stages themselves at a landscape scale across the examined permafrost thaw gradient. Future studies should thus not only focus on changes in carbon budget associated with morphological changes in thawing permafrost, but also recognize the effects of climate forcing uncertainty on carbon cycling.

54 ENVIRONMENTAL SCIENCES↗

Yet another parameter-free shape optimization method

The use of node coordinates as design variables in shape optimization offers a larger design space than computer-aided design (CAD)-based shape parameterizations. It also allows for the optimization of legacy designs, i.e., a finite element mesh from an existing design can be readily optimized to meet new performance requirements without involving a CAD model. However, it is well known that the node coordinate parameterization method is fraught with numerical difficulties, which makes it impractical to use. This has led to several of “parameter-free” shape optimization methods that seek the advantages and avoid the pitfalls of the naïve node coordinate parameterization method. These methods come in two main varieties: sensitivity filtering (or gradient smoothing) and consistent filtering. The latter is analogous to the density filter method used in topology optimization (TO). In this work, we use the PDE filter from TO and energy-based filters to implement consistent shape optimization filtering schemes easily and efficiently. Numerical experiments demonstrate that consistent methods are more robust than sensitivity filtering methods.

42 ENGINEERING↗

Sensitivity analysis for characterizing the impact of HNGD model on the prediction of hydrogen redistribution in Zircaloy cladding using BISON code

Hydrogen in zirconium cladding is able to precipitate into zirconium hydrides which impacts cladding integrity. The Hydride Nucleation-Growth-Dissolution (HNGD) model in the BISON code accounts for the precipitation and dissolution kinetics of hydride in Zircaloy material. This paper presents global sensitivity analyses of the HNGD model aiming to enhance our understanding of the hydride precipitation phenomena by quantifying the variance that key parameters have on the prediction of hydrogen behavior under various environmental conditions. Here, model predictions are compared to experimental data obtained under two different conditions: 1) with uniformly precharged specimens subjected to a linear thermal gradient, and 2) specimens precharged with a cathodically applied hydride rim at one end of the sample and subjected to an asymmetric thermal gradient. The Sobol sensitivity analysis identifies the key parameters in the HNGD model for both types of specimens. For linear temperature cases, the heat of transport dominates the accuracy of predictions when no precipitation occurs at the cold end, while Terminal Solid Solubility for Dissolution (TSSD) is the most important parameter when precipitation occurs. A large variation in the predicted hydrogen concentration profiles is found in the range of high TSSD due to the occurrence of precipitation. For asymmetric temperature cases, the solubility coefficient gives the largest impact on the predicted hydrogen distribution, as it determines the amount of solute hydrogen dissolved from the initially applied hydride rim. A large discrepancy in hydrogen distribution between simulations and experiments exists with the asymmetric specimens because BISON simulations fail to predict the precipitation of hydride at the cooler end. Comparative studies using former and updated models verifies the significant impact of the hydride growth mechanism on predicted hydrogen concentration profiles. In particular, when hydride initially exists, changes in TSSD generate a large variation in the predicted amount of precipitation by hydride growth, giving large uncertainty in predicting the hydrogen distribution over the sample length. The outputs characterize the significant impact of the hydride growth mechanism in the HNGD model on predicting hydrogen behavior, and improve the understanding of the precipitation of hydride in Zircaloy cladding within a range of expected environmental conditions. The analyses indicate work is still needed to improve the hydride solvus models in the BISON code to accurately predict experimentally observed hydride concentrations and distributions.

36 MATERIALS SCIENCE↗

Stochastic properties of ultralight scalar field gradients

Ultralight axion-like particles are well-motivated dark matter candidates that are the target of numerous direct detection efforts. In the vicinity of the Solar System, such particles can be treated as oscillating scalar fields. The velocity dispersion of the Milky Way determines a coherence time of about 10 6 oscillations, beyond which the amplitude of the axion field fluctuates stochastically. Any analysis of data from an axion direct detection experiment must carefully account for this stochastic behavior to properly interpret the results. This is especially true for experiments sensitive to the gradient of the axion field that are unable to collect data for many coherence times. Indeed, the direction, in addition to the amplitude, of the axion field gradient fluctuates stochastically. We present the first complete stochastic treatment for the gradient of the axion field, including multiple computationally efficient methods for performing likelihood-based data analysis, which can be applied to any axion signal, regardless of coherence time. Additionally, we demonstrate that ignoring the stochastic behavior of the gradient of the axion field can potentially result in failure to discover a true axion signal

79 ASTRONOMY AND ASTROPHYSICS↗

Aerodynamic design optimization using sensitivity analysis and computational fluid dynamics

A new and efficient method is presented for aerodynamic design optimization, which is based on a computational fluid dynamics (CFD)-sensitivity analysis algorithm. The method is applied to design a scramjet-afterbody configuration for an optimized axial thrust. The Euler equations are solved for the inviscid analysis of the flow, which in turn provides the objective function and the constraints. The CFD analysis is then coupled with the optimization procedure that uses a constrained minimization method. The sensitivity coefficients, i.e. gradients of the objective function and the constraints, needed for the optimization are obtained using a quasi-analytical method rather than the traditional brute force method of finite difference approximations. During the one-dimensional search of the optimization procedure, an approximate flow analysis (predicted flow) based on a first-order Taylor series expansion is used to reduce the computational cost. Finally, the sensitivity of the optimum objective function to various design parameters, which are kept constant during the optimization, is computed to predict new optimum solutions. The flow analysis of the demonstrative example are compared with the experimental data. It is shown that the method is more efficient than the traditional methods.

Baysal, Oktay↗

Evaluation of the Effect of Burnup Acceleration on UO 2 Microstructure Evolution

Accelerated fuel qualification (AFQ) is a methodology by which new nuclear fuels are developed in an accelerated time frame compared with historical fuel qualification approaches. AFQ generally relies on high-fidelity physics-based modeling and simulation tools to adequately describe fuel performance as well as on revolutionary methods to accelerate burnup accumulation and collect relevant data more quickly. This report summarizes the use of advanced fuel modeling and simulation tools to evaluate microstructures from commercially irradiated fuel and microstructures from proposed MiniFuel irradiations, in which burnup accumulation is accelerated while prototypic temperature conditions are maintained. In this milestone, we used the mesoscale fuel performance code MARMOT to model the evolution of irradiated UO 2 microstructures and their potential restructuring at high burnup. The simulation conditions were informed by BISON models of both commercially irradiated fuel and MiniFuel. A first set of simulations investigated the recrystallization behavior of fully dense microstructures and showed full recrystallization at burnups as low as 52 MWd/kgU at 950°C. However, these simulations did not account for the presence of fission gas bubbles (predicted by BISON). Therefore, a second set of simulations including fission gas bubbles was performed and indicated that at the lowest temperature considered (650°C), the porous UO 2 microstructures have the highest total Gibbs free energies and are likely to recrystallize earlier than higher temperature cases (800 and 950°C), which agrees with high-burnup fuel characterization data. The results also showed that at lower temperature (650°C), the total free energies of the PWR fuel and MiniFuel microstructures are not significantly different. However, at the highest temperature (950°C), MiniFuel microstructures have a lower free energy than that of the PWR fuel microstructure. The competing effects between the temperature-dependent grain nucleation rate and the reduction of the free energy of the microstructure at higher temperature as a result of diffusion indicated that restructuring may occur at even higher temperatures than those considered in this study. In addition to the microstructure evolution modeling efforts, the burnup gradient across a single fuel specimen was also considered. This evaluation was for the VXF-15 position of the High Flux Isotope Reactor (HFIR) using the code suite HFIRCON, which was developed to automate the workflow for evaluating targets and fuel as they are irradiated in HFIR. The burnup gradient evaluation showed a dependence on both the axial and radial locations within the specimen, with a maximum difference of 1.7 between the inner and outermost radial layers. This relationship was further supported by considering the fission product speciation with respect to location within the specimen, which showed a higher concentration of 239 Pu, 240 Pu, and 241 Pu on the outer radial locations of the specimen than the center. The findings of the burnup and speciation evaluation show that some amount of self-shielding is occurring in the specimen when irradiated in the high-flux environment of HFIR; however, this impact is more pronounced for natural uranium when compared to 6% enrichment due to the higher ratio of 238 U in the specimen. Further analyses are required to understand the sensitivity of this gradient to spatial mesh and enrichment of the specimen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Bulk and Nanocrystalline Cesium Lead-Halide Perovskites as Seen by Halide Magnetic Resonance

Lead-halide perovskites increasingly mesmerize researchers because they exhibit a high degree of structural defects and dynamics yet nonetheless offer an outstanding (opto)electronic performance on par with the best examples of structurally stable and defect-free semiconductors. This highly unusual feature necessitates the adoption of an experimental and theoretical mindset and the reexamination of techniques that may be uniquely suited to understand these materials. Surprisingly, the suite of methods for the structural characterization of these materials does not commonly include nuclear magnetic resonance (NMR) spectroscopy. The present study showcases both the utility and versatility of halide NMR and NQR (nuclear quadrupole resonance) for probing the structure and structural dynamics of CsPbX 3 (X = Cl, Br, I), in both bulk and nanocrystalline forms. The strong quadrupole couplings, which originate from the interaction between the large quadrupole moments of, e.g., the 35 Cl, 79 Br, and 127 I nuclei, and the local electric-field gradients, are highly sensitive to subtle structural variations, both static and dynamic. The quadrupole interaction can resolve structural changes with accuracies commensurate with synchrotron X-ray diffraction and scattering. It is shown that space-averaged site-disorder is greatly enhanced in the nanocrystals compared to the bulk, while the dynamics of nuclear spin relaxation indicates enhanced structural dynamics in the nanocrystals. The findings from NMR and NQR were corroborated by ab initio molecular dynamics, which point to the role of the surface in causing the radial strain distribution and disorder. These findings showcase a great synergy between solid-state NMR or NQR and molecular dynamics simulations in shedding light on the structure of soft lead-halide semiconductors.

36 MATERIALS SCIENCE↗

A non-intrusive parallel-in-time approach for simultaneous optimization with unsteady PDEs

This paper presents a non-intrusive framework for integrating existing unsteady partial differential equation (PDE) solvers into a parallel-in-time simultaneous optimization algorithm. The time-parallelization is provided by the non-intrusive software library XBraid [Parallel multigrid in time, software available at http://llnl.gov/casc/xbraid], which applies an iterative multigrid reduction technique to the time domain of existing time-marching schemes for solving unsteady PDEs. Its general user-interface has been extended in [S. Günther, N.R. Gauger, and J.B. Schroder, A non-intrusive parallel-in-time adjoint solver with the XBraid library, Comput. Vis. Sci. 19 (2018), pp. 85–95. Available at arXiv, math.OC/1705.00663] for computing adjoint sensitivities such that gradients of output quantities with respect to design changes can be computed parallel-in-time alongside with the primal PDE solution. In this paper, the primal and adjoint XBraid iterations are embedded into a simultaneous optimization framework, namely the One-shot method. In this method, design updates towards optimality are employed after each state and adjoint update such that optimality and feasibility of the design and the PDE solution are reached simultaneously. Finally, the time-parallel optimization method is validated on an advection-dominated flow control problem which shows significant speedup over a classical time-serial optimization algorithm.

97 MATHEMATICS AND COMPUTING↗

Probing primordial power spectrum and non-Gaussianities with fast radio bursts

We use the precision measurements of the arrivaltime differences of the same fast radio burst (FRB) source along multiple sightlines to measure the primordial power spectrum and Non-Gaussianities. The anticipated experiment requires a sightline separation of 100 AU, achieved by sending three or more radio telescopes to the outer solar system.The Shapiro time delays, measured relatively between different telescopes, are sensitive to the gradient field of the gravitational potential between different sightlines. Since the arrival time difference is independent of when the transient signal is emitted from the source, every measurement of the detected FRB source can be correlated. With enough FRB sources discovered, we can map the gravitational potential across the sky. We further calculate the two-point and three-point correlation function of the arrival time difference between telescopes for different FRB sources in the sky. If 10$^{4}$ FRBs were to be detected, our results suggest that this technique can test the inflationary scale-invariant power spectrum down to ∼ 10$^{3}$ Mpc$^{-1}$ and primordial Non-Gaussianities at a level of f$_{NL}$ ∼ 1.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evaluating probabilistic deep learning methods for uncertainty quantification of temperature downscaling

Deep learning (DL) has emerged as a promising tool for downscaling coarse-resolution climate data to high-resolution outputs, enabling improved regional climate predictions. A critical aspect of DL-based downscaling is the incorporation of uncertainty quantification (UQ), which enhances the interpretability and reliability of predictions—key factors for climate risk assessment and decision-making. This study develops a DL model to downscale 2 m temperature across the contiguous United States using reanalysis datasets. We systematically evaluate three epistemic UQ methods—deep ensembles (DEns), Monte Carlo dropout (MCD), and Flipout—based on their probabilistic accuracy, downscaling performance, sensitivity to geographical features, and computational efficiency. Results indicate that MCD generally outperforms Flipout and DEns in terms of calibration and downscaling accuracy. However, DEns demonstrate lower calibration errors in coastal regions, indicating its higher confidence within these areas. Flipout, in contrast, is more sensitive to elevation gradients and exhibits higher calibration errors in mountainous regions. Hence, the choice of UQ method for this task depends on the specific requirements of the application. For applications that prioritize overall calibration, downscaling accuracy, and computational efficiency, MCD is a strong candidate. These findings highlight the importance of selecting UQ methods based on application-specific requirements, such as geographical context and computational constraints. By addressing the trade-offs between UQ methods, this study provides actionable insights for improving the reliability, scalability, and utility of DL-based downscaling in climate science.

Environmental sciences↗

Velocity Profiling of a Gas–Solid Fluidized Bed Using Electrical Capacitance Volume Tomography

In this study, a method of producing velocity profile maps from electrical capacitance volume tomography (ECVT) measurements by reconstructing displacement from measured changes in capacitance is developed and applied to fluidized bed systems. The mapping of the reconstruction leverages the gradient of the sensitivity distribution of the ECVT sensor to circumvent the need for image cross correlation techniques. Experimental data of both bubbling and slugging fluidized beds are collected in a cold flow model. Adaptation of the technique is discussed in detail, and velocity profiles are obtained for a range of gas flow rates. The produced velocity maps are compared against the established methods of cross correlation and against empirical correlations from the literature and are found to agree well in tracking slug and bubble velocity. The exception is when the tracked object is large relative to the ECVT sensor dimensions, a scenario that can be avoided through proper sensor design. The quantities of average velocity, momentum, and solid and gas volume fraction are derived from the image and velocity profiles. The results demonstrate and extend the power of ECVT as a measurement tool for the study and monitoring of gas–solid fluidized beds by providing a computationally cheaper alternative to 3-D cross correlation for deriving velocity profiles.

47 OTHER INSTRUMENTATION↗

Effect of Surface Roughness on Wind Turbine Performance

Wind farm operators observe production deficits as machines age. Quantifying deterioration on individual components is difficult, but one potential explanation is accumulation of blade surface roughness. Historically, wind turbine airfoils were designed for lift to be insensitive to roughness by simulating roughness with trip strips. However, roughness was still shown to negatively affect performance. Furthermore, experiments illustrated distributed roughness is not properly simulated by trip strips. To understand how real-world roughness affects performance, field measurements of turbine-blade roughness were made and simulated on a NACA 63 3 -418 airfoil in a wind tunnel. Insect roughness and paint chips were characterized and recreated as distributed roughness and a forward-facing step. Distributed roughness was tested in three heights and five density configurations. The model chord Reynolds number was varied between 0.8 to 4.8 x 10 6 . Measurements of lift, drag, pitching moment, and boundary-layer transition were completed. Results indicate minimal effect from paint-chip roughness. As distributed roughness height and density increase, lift-curve slope, maximum lift, and lift-to-drag ratio decrease. As Reynolds number increases, bypass transition occurs earlier. The critical roughness Reynolds number varies between 178 to 318, within the historical range. Little sensitivity to pressure gradient is observed. At a chord Reynolds number of 3.2 x 10 6 , the maximum lift-to-drag ratio decreases 40% for 140 m roughness, corresponding to a 2.3% loss in annual energy production. Simulated performance loss compares well to measured performance loss on an in-service wind turbine.

17 WIND ENERGY↗

Static-pressure probe for small geometries

Contoured pressure probe with static orifices located near tip is more effective than conventional probes in taking measurements in small high-Reynolds number geometries such as nozzles or flow inlets. Probe is less sensitive to pressure gradients and off axis variations in flow directions over short distances.

Pinckney, S. Z.↗

IUE observations of the 1982-84 eclipse of Epsilon Aurigae

The major characteristics in the ultraviolet of the 1982-84 eclipse of Eps Aur as observed with IUE by various workers are summarized. This star can be observed over the entire IUE wavelength range, from 1200 to 3200 A, in low dispersion, allowing eclipse light curves to be obtained in broadband regions, but due to its steep spectral gradient and the sensitivities of IUE cameras, high resolution exposures adequately cover only the regions from 1700-1900 and 2400-3200 A. In many ways, the UV data confirms or expands upon interpretations of the stem made from observations in other wavelength regions, but in other respects the system remains as enigmatic as before.

Ake, T. B.↗

Influence of Rayleigh-Doppler broadening on the selection of H2O dial system parameters

Computer simulations have enabled the performance of a H2O Differential Absorption Lidar (DIAL) system to be studied by spectrally analyzing the forward propagating and backscattered laser energy. The simulations were done for a high altitude (21 km) DIAL system operating in a nadir-viewing mode. The influence of Rayleigh Doppler broadening on DIAL measurement accuracies were evaluated and show that the Rayleigh broadening influence, which can be corrected to first order in regions free of large aerosol gradients, reduces the sensitivity of DIAL H2O measurement errors in the upper tropospheric region. The ability to correct the Rayleigh broadening and the selection of H2O DIAL parameters when all the systematic effects are combined, were discussed.

Ismail, S.↗