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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Velocity Profiling of Multiphase Flows Using Capacitive Sensor Sensitivity Gradient

Velocity profiling of a flow involves the task of determining the velocity vector at every point in a given flow volume. A new method is proposed for velocity profiling of multiphase flows based on electrical capacitance volume tomography (ECVT) sensors. The proposed method utilizes a mapping between the change in measured capacitances and the displacement of flow that is effected by the spatial gradient of the sensitivity distribution. Here, this novel mapping not only avoids the need for costly image cross-correlations but also is fully compatible with existing ECVT sensor and image reconstruction algorithms. Simulation and measurement results are provided to demonstrate the proposed method.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Inverse Normalization Method for Cross-Sectional Imaging and Velocimetry of Two-Phase Flows Based on Electrical Capacitance Tomography

Proper normalization of the measured capacitance data is a prerequisite in electrical capacitance tomography (ECT). Conventional methods, such as parallel and series normalization, assume a higher permittivity perturbation in a homogeneous (lower) permittivity background. Although this is applicable to a number of cases, some two-phase flows are better modeled instead as a lower permittivity perturbation in a homogeneous (higher) permittivity background. A different normalization method is proposed for such cases, which provides similar image quality as conventional, but is particularly beneficial for velocimetry applications based on sensitivity gradient of ECT sensor. Furthermore, simulation and experimental results are provided to illustrate the advantages for the proposed normalization technique.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Nanoscale electric field imaging with an ambient scanning quantum sensor microscope

Nitrogen-vacancy (NV) center in diamond is a promising quantum sensor with remarkably versatile sensing capabilities. While scanning NV magnetometry is well-established, NV electrometry has been so far limited to bulk diamonds. Here we demonstrate imaging external alternating (AC) and direct (DC) electric fields with a single NV at the apex of a diamond scanning tip under ambient conditions. A strong electric field screening effect is observed at low frequencies. We quantitatively measure its frequency dependence and overcome this screening by mechanically oscillating the tip for imaging DC fields. Our scanning NV electrometry achieved an AC E-field sensitivity of 26 mV μm -1 Hz -1/2 , a DC E-field gradient sensitivity of 2 V μm -2 Hz -1/2 , and sub-100 nm resolution limited by the NV-sample distance. Our work represents an important step toward building a scanning-probe-based multimodal quantum sensing platform.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Parameterization retrieval of trace gas volume mixing ratios from Airborne MAX-DOAS

We present a parameterization retrieval of volume mixing ratios (VMRs) from differential slant column density (dSCD) measurements by Airborne Multi-AXis Differential Optical Absorption Spectroscopy (AMAX-DOAS). The method makes use of the fact that horizontally recorded limb spectra (elevation angle 0°) are strongly sensitive to the atmospheric layer at instrument altitude. These limb spectra are analyzed using reference spectra that largely cancel out column contributions from above and below the instrument, so that the resulting limb dSCDs, i.e., the column integrated concentration with respect to a reference spectrum, are almost exclusively sensitive to the atmospheric layers around instrument altitude. The conversion of limb dSCDs into VMRs is then realized by calculating box air mass factors (Box-AMFs) for a Rayleigh atmosphere and applying a scaling factor constrained by O 4 dSCDs to account for aerosol extinction. An iterative VMR retrieval scheme corrects for trace gas profile shape effects. Benefits of this method are (1) a fast conversion that only requires the computation of Box-AMFs in a Rayleigh atmosphere; (2) neither local aerosol extinction nor the slant column density in the DOAS reference (SCD ref ) needs to be known; and (3) VMRs can be retrieved for every measurement point along a flight track, thus increasing statistics and adding flexibility to capture concentration gradients. Sensitivity studies are performed for bromine monoxide (BrO), iodine monoxide (IO) and nitrogen dioxide (NO 2 ), using (1) simulated dSCD data for different trace gas and aerosol profiles and (2) field measurements from the Tropical Ocean tRoposphere Exchange of Reactive halogen species and Oxygenated VOC (TORERO) field experiment. For simulated data in a Rayleigh atmosphere, the agreement between the VMR from the parameterization method (VMR para ) and the true VMR (VMR true ) is excellent for all trace gases. Offsets, slopes and R 2 values for the linear fit of VMR para over VMR true are, respectively (0.008 ± 0.001) pptv, 0.988 ± 0.001, 0.987 for BrO; (-0.0066 ± 0.0001) pptv, 1.0021 ± 0.0003, 0.9979 for IO; (-0.17 ± 0.03) pptv, 1.0036 ± 0.0001, 0.9997 for NO 2 . The agreement for atmospheres with aerosol shows comparable R 2 values to the Rayleigh case, but slopes deviate a bit more from one: (0.093 ± 0.002) pptv, 0.933 ± 0.002, 0.907 for BrO; (0.0021 ± 0.0004) pptv, 0.887 ± 0.001, 0.973 for IO; (8.5 ± 0.1) pptv, 0.8302 ± 0.0006, 0.9923 for NO 2 . VMR para from field data are further compared with optimal estimation retrievals (VMR OE ). Least orthogonal distance fit of the data give the following equations: BrO para = (0.1 ± 0.2) pptv + (0.95 ± 0.14) × BrO OE ; IO para = (0.01 ± 0.02) pptv + (1.00 ± 0.12) × IO OE ; NO 2para = (3.9 ± 2.5) pptv + (0.87 ± 0.15) × NO 2OE . Overall, we conclude that the parameterization retrieval is accurate with an uncertainty of 20 % for IO, 30 % for BrO and NO 2 , but not better than 0.05 pptv IO, 0.5 pptv BrO and 10 pptv NO 2 . Finally, the retrieval is applicable over a wide range of atmospheric conditions and measurement geometries and not limited to the interpretation of vertical profile measurements in the remote troposphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An optimization method for chaotic turbulent flow

Evidence indicates that quantities-of-interest in some turbulent flows can be controlled despite the overall chaotic dynamics. It is typically thought that this is via relatively deterministic, larger-scale components of the turbulence. However, finding such controls, if they exist, is challenging because chaos causes sensitivity gradients to explode and the search space to become intractably non-convex. This challenge is analyzed, and a penalty method is introduced to cope with it. In the new approach, the time domain is broken into segments approximately matching the chaos time scales, so that the solution within each segment is both physical and relatively deterministic. The initial condition of each segment is included in an adjoint-based gradient optimization, which temporarily introduces artificial Δq discontinuities in the overall solution. The optimization then proceeds in stages with increasing penalization of Δq. The method is developed and illustrated for a logistic map, the Lorenz Equation, and an advection augmented Kuramoto–Sivashinsky Equation. These examples show how the Δq temporarily increases the search scale prior to the strong Δq → 0 penalization that recovers a physical solution. It is then applied to turbulent Kolmogorov flow, for which it also far outperforms a standard adjoint-based gradient search. Finally, the utility of such an optimized chaotic solution is discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Incorporating Elevation in Traffic-Vehicle CO-Simulation: Issues, Impacts, and Solutions

Traffic-vehicle co-simulation couples microscopic traffic simulation with full-body vehicle dynamics to assess system-level impacts on mobility, energy, and safety with greater realism. Incorporating elevation is critical for accurately modeling vehicle behavior and energy use, especially for gradient-sensitive vehicles such as electric and heavy-duty trucks. However, raw elevation data often contain noise, discontinuities, and inconsistencies. While such issues may be negligible in traditional traffic simulations, they significantly affect traffic-vehicle co-simulations where vehicle dynamics are sensitive to road grade variations. This paper investigates the impact of unprocessed elevation data on vehicle behavior and energy consumption using a 42-mile simulation along Interstate 81. We propose an elevation processing workflow that can mitigate the effects stem from elevation data issues, improving the realism and stability of traffic-vehicle co-simulation. Results show that the method effectively removes noise and abrupt elevation transitions while preserving roadway geometry.

Xu, Guanhao [ORNL] (ORCID:0000000214326357)↗
DOE OSTI2023

Development of compact tokamak fusion reactor use cases to inform future transport studies

The OMFIT STEP (Meneghini et al. , Nucl. Fusion , vol. 10, 2020, p. 1088) workflow has been used to develop inductive and steady-state H-mode core plasma scenario use cases for a $B_0 = 8 \, {\rm T}$ , $R_0 = 4 \, {\rm m}$ machine to help guide and inform future higher-fidelity studies of core transport and confinement in compact tokamak reactors. Both use cases are designed to produce 200 MW or more of net electric power in an up-down symmetric plasma with minor radius $a = 1.4 \, {\rm m}$ , elongation $\kappa = 2.0$ , triangularity $\delta = 0.5$ and effective charge $Z_{{\rm eff}} \simeq 2$ . Additional considerations based on the need for compatibility of the core with reactor-relevant power exhaust solutions and external actuators were used to guide and constrain the use case development. An extensive characterization of core transport in both scenarios is presented, the most important feature of which is the extreme sensitivity of the results to the quantitative stiffness level of the transport model used as well as the predicted critical gradients. This sensitivity is shown to arise from different levels of transport stiffness exhibited by the models, combined with the gyroBohm-normalized fluxes of the predictions being an order of magnitude larger than other H-mode plasmas. Additionally, it is shown that although heating in both plasmas is predominantly to the electrons and collisionality is low, the plasmas remain sufficiently well coupled for the ions to carry a significant fraction of the thermal transport. As neoclassical transport is negligible in these conditions, this situation inherently requires long-wavelength ion gyroradius-scale turbulence to be the dominant transport mechanism in both plasmas. These results are combined with other basic considerations to propose a simple heuristic model of transport in reactor-relevant plasmas, along with simple metrics to quantify coupling and core transport properties across burning and non-burning plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Liquid Crystal Orientation and Shape Optimization for the Active Response of Liquid Crystal Elastomers

Liquid crystal elastomers (LCEs) are responsive materials that can undergo large reversible deformations upon exposure to external stimuli, such as electrical and thermal fields. Controlling the alignment of their liquid crystals mesogens to achieve desired shape changes unlocks a new design paradigm that is unavailable when using traditional materials. While experimental measurements can provide valuable insights into their behavior, computational analysis is essential to exploit their full potential. Accurate simulation is not, however, the end goal; rather, it is the means to achieve their optimal design. Such design optimization problems are best solved with algorithms that require gradients, i.e., sensitivities, of the cost and constraint functions with respect to the design parameters, to efficiently traverse the design space. In this work, a nonlinear LCE model and adjoint sensitivity analysis are implemented in a scalable and flexible finite element-based open source framework and integrated into a gradient-based design optimization tool. To display the versatility of the computational framework, LCE design problems that optimize both the material, i.e., liquid crystal orientation, and structural shape to reach a target actuated shapes or maximize energy absorption are solved. Multiple parameterizations, customized to address fabrication limitations, are investigated in both 2D and 3D. The case studies are followed by a discussion on the simulation and design optimization hurdles, as well as potential avenues for improving the robustness of similar computational frameworks for applications of interest.

42 ENGINEERING↗
DOE OSTI2016

Three axis vector atomic magnetometer utilizing polarimetric technique

The three axis vector magnetic field measurement based on the interaction of a single elliptically polarized light beam with an atomic system is described. The magnetic field direction dependent atomic responses are extracted by the polarimetric detection in combination with laser frequency modulation and magnetic field modulation techniques. The magnetometer geometry offers additional critical requirements like compact size and large dynamic range for space application. Further, the three axis magnetic field is measured using only the reflected signal (one polarization component) from the polarimeter and thus can be easily expanded to make spatial array of detectors and/or high sensitivity field gradient measurement as required for biomedical application.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

A comparative study of calibration techniques for finite strain elastoplasticity: Numerically-exact sensitivities for FEMU and VFM

Accurate identification of material parameters is crucial for predictive modeling in computational mechanics. Here, the two primary approaches in the experimental mechanics community for calibration from full-field digital image correlation data are known as finite element model updating (FEMU) and the virtual fields method (VFM). In VFM, the objective function is a squared mismatch between internal and external virtual work or power. In FEMU, the objective function quantifies the weighted mismatch between model predictions and corresponding experimentally measured quantities of interest. It is minimized by iteratively updating the parameters of an FE model. While FEMU is seen as more flexible, VFM is commonly used instead of FEMU due to its considerably greater computational expense. However, comparisons between the two methods usually involve approximations of gradients or sensitivities with finite difference schemes, thereby making direct assessments difficult. Hence, in this study, we compare VFM and FEMU in the context of numerically-exact sensitivities obtained through local sensitivity analyses and the application of automatic differentiation software. To this end, we conduct a series of test cases to assess both methods under practical challenges using a finite strain elastoplasticity model.

Automatic differentiation↗
DOE OSTI2022

Approximating the linear response of physical chaos

Abstract Parametric derivatives of statistics are highly desired quantities in prediction, design optimization and uncertainty quantification. In the presence of chaos, the rigorous computation of these quantities is certainly possible, but mathematically complicated and computationally expensive. Based on Ruelle’s formalism, this paper shows that the sophisticated linear response algorithm can be dramatically simplified in higher-dimensional systems featuring statistical homogeneity in the physical space. We argue that the contribution of the SRB (Sinai–Ruelle–Bowen) measure gradient, which is an integral yet the most cumbersome part of the full algorithm, is negligible if the objective function is appropriately aligned with unstable manifolds. This abstract condition could potentially be satisfied by a vast family of real-world chaotic systems, regardless of the physical meaning and mathematical form of the objective function and perturbed parameter. We demonstrate several numerical examples that support these conclusions and that present the use and performance of a simplified linear response algorithm. In the numerical experiments, we consider physical models described by differential equations, including Lorenz 96 and Kuramoto–Sivashinsky.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Plasma characterization using ultraviolet Thomson scattering from ion-acoustic and electron plasma waves (invited)

Collective Thomson scattering is a technique for measuring the plasma conditions in laser-plasma experiments. Simultaneous measurements of ion-acoustic and electron plasma-wave spectra were obtained using a 263.25-nm Thomson-scattering probe beam. A fully reflective collection system was used to record light scattered from electron plasma waves at electron densities greater than 10{sup 21} cm{sup −3}, which produced scattering peaks near 200 nm. An accurate analysis of the experimental Thomson-scattering spectra required accounting for plasma gradients, instrument sensitivity, optical effects, and background radiation. Practical techniques for including these effects when fitting Thomson-scattering spectra are presented and applied to the measured spectra to show the improvements in plasma characterization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Phase retrieval for refraction-enhanced x-ray radiography using a deep neural network

X-ray refraction-enhanced radiography (RER) or phase contrast imaging is widely used to study internal discontinuities within materials. The resulting radiograph captures both the decrease in intensity caused by material absorption along the x-ray path, as well as the phase shift, which is highly sensitive to gradients in density. A significant challenge lies in effectively analyzing the radiographs to decouple the intensity and phase information and accurately ascertain the density profile. Conventional algorithms often yield ambiguous and unrealistic results due to difficulties in including physical constraints and other relevant information. We have developed an algorithm that uses a deep neural network to address these issues and applied it to extract the detailed density profile from an experimental RER. To generalize the applicability of our algorithm, we have developed a technique that quantitatively evaluates the complexity of the phase retrieval process based on the characteristics of the sample and the configuration of the experiment. Accordingly, this evaluation aids in the selection of the neural network architecture for each specific case. Beyond RER, the model has potential applications for other diagnostics where phase retrieval analysis is required.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Toward a Holistic Performance Evaluation of Large Language Models Across Diverse AI Accelerators

Artificial intelligence (AI) methods have become critical in scientific applications to help accelerate scientific discovery. Large language models (LLMs) are being considered a promising approach to address some challenging problems because of their superior generalization capabilities across domains. The effectiveness of the models and the accuracy of the applications are contingent upon their efficient execution on the underlying hardware infrastructure. Specialized Al accelerator hardware systems have recently become available for accelerating Al applications. However, the comparative performance of these AI accelerators on large language models has not been previously studied. In this paper, we systematically study LLMs on multiple AI accelerators and GPUs and evaluate their performance characteristics for these models. We evaluate these systems with (i) a micro-benchmark using a core transformer block, (ii) a GPT-2 model, and (iii) an 1,I,M-driven science use case, GenSLM. We present our findings and analyses of the models' performance to better understand the intrinsic capabilities of AI accelerators. Furthermore, our analysis takes into account key factors such as sequence lengths, scaling behavior, and sensitivity to gradient accumulation steps.

Emani, Murali↗
DOE OSTI2023

Lowering the computational barrier: Partially Bayesian neural networks for transparency in medical imaging AI

Deep Neural Networks (DNNs) can provide clinicians with fast and accurate predictions that are highly valuable for high-stakes medical decision-making, such as in brain tumor segmentation and treatment planning. However, these models largely lack transparency about the uncertainty in their predictions, potentially giving clinicians a false sense of reliability that may lead to grave consequences in patient care. Growing calls for Transparent and Responsible AI have promoted Uncertainty Quantification (UQ) to capture and communicate uncertainty in a systematic and principled manner. However, traditional Bayesian UQ methods remain prohibitively costly for large, million-dimensional tumor segmentation DNNs such as the U-Net. In this work, we discuss a computationally-efficient UQ approach via the partially Bayesian neural networks (pBNN). In pBNN, only a single layer, strategically selected based on gradient-based sensitivity analysis, is targeted for Bayesian inference. We illustrate the effectiveness of pBNN in capturing the full uncertainty for a 7.8-million parameter U-Net. We also demonstrate how practitioners and model developers can use the pBNN's predictions to better understand the model's capabilities and behavior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Patient-Specific Mathematical Model of the Clear Cell Renal Cell Carcinoma Microenvironment

The interactions between cells and molecules in the tumor microenvironment can give insight into the initiation and progression of tumors and their optimal treatment options. In this paper, we developed an ordinary differential equation (ODE) mathematical model of the interaction network of key players in the clear cell renal cell carcinoma (ccRCC) microenvironment. We then performed a global gradient-based sensitivity analysis to investigate the effects of the most sensitive parameters of the model on the number of cancer cells. The results indicate that parameters related to IL-6 have high a impact on cancer cell growth, such that decreasing the level of IL-6 can remarkably slow the tumor’s growth.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Sensitivity of magnetic islands in permanent magnet stellarators using the gradient and Hessian methods

Stellarator plasmas are known to be very sensitive to perturbations in the magnetic field. The permanent magnet stellarator was in part developed as a solution to high machining tolerances placed on the shape properties of electromagnetic coils in traditional stellarators. However, as a consequence of this high sensitivity to the field structure, sensitivities of permanent magnet stellarator plasmas to perturbations of permanent magnet properties must necessarily be well-understood. The gradient and Hessian matrix methods have been previously demonstrated to be useful sensitivity analysis methods for modular coils. We apply these two methods to the study of island width sensitivities in both the MUSE and PM4STELL permanent magnet stellarator projects. These sensitivity methods were used to determine the relative impacts of permanent magnet parameter perturbations on island widths in the vacuum field approximation of both stellarator equilibria. The square of resonant magnetic field perturbation is used here as a proxy for island width. In particular, gradients of magnetizations of individual magnets were examined in MUSE, as well as gradients of magnet group displacements informed by device design. Three different forms of permanent magnet magnetization perturbations are investigated for MUSE, and the flux surface response to perturbations is demonstrated. The Hessian matrix method is applied to PM4STELL, illustrating the sensitivity of dominant island widths to displacements of toroidal wedge structures. These methods allow for selective direction of experimental resources toward regions of heightened sensitivity, while constraints on less impactful permanent magnet parameters can be relaxed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria

The shape gradient quantifies the change in some figure of merit resulting from differential perturbations to a shape. Shape gradients can be applied to gradient-based optimization, sensitivity analysis and tolerance calculation. An efficient technique for computing the shape gradient for toroidal three-dimensional magnetohydrodynamic (MHD) equilibria is introduced. The method is based on the self-adjoint property of the equations for driven perturbations of MHD equilibria and is similar to the Onsager symmetry of transport coefficients. Two versions of the shape gradient are considered. One describes the change in a figure of merit due to an arbitrary displacement of the outer flux surface; the other describes the change in the figure of merit due to the displacement of a coil. The method is implemented for several example figures of merit and compared with direct calculation of the shape gradient. In these examples the adjoint method reduces the number of equilibrium computations by factors of $O(N)$ , where $N$ is the number of parameters used to describe the outer flux surface or coil shapes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗