Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient calculation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Lattice parameters and electronic structure of BeMgZnO quaternary solid solutions: Experiment and theory

Be{sub x}Mg{sub y}Zn{sub 1−x−y}O semiconductor solid solutions are attractive for UV optoelectronics and electronic devices owing to their wide bandgap and capability of lattice-matching to ZnO. In this work, a combined experimental and theoretical study of lattice parameters, bandgaps, and underlying electronic properties, such as changes in band edge wavefunctions in Be{sub x}Mg{sub y}Zn{sub 1−x−y}O thin films, is carried out. Theoretical ab initio calculations predicting structural and electronic properties for the whole compositional range of materials are compared with experimental measurements from samples grown by plasma assisted molecular beam epitaxy on (0001) sapphire substrates. The measured a and c lattice parameters for the quaternary alloys Be{sub x}Mg{sub y}Zn{sub 1−x} with x = 0−0.19 and y = 0–0.52 are within 1%–2% of those calculated using generalized gradient approximation to the density functional theory. Additionally, composition independent ternary BeZnO and MgZnO bowing parameters were determined for a and c lattice parameters and the bandgap. The electronic properties were calculated using exchange tuned Heyd-Scuseria-Ernzerhof hybrid functional. The measured optical bandgaps of the quaternary alloys are in good agreement with those predicted by the theory. Strong localization of band edge wavefunctions near oxygen atoms for BeMgZnO alloy in comparison to the bulk ZnO is consistent with large Be-related bandgap bowing of BeZnO and BeMgZnO (6.94 eV). The results in aggregate show that precise control over lattice parameters by tuning the quaternary composition would allow strain control in Be{sub x}Mg{sub y}Zn{sub 1−x−y}O/ZnO heterostructures with possibility to achieve both compressive and tensile strain, where the latter supports formation of two-dimensional electron gas at the interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1976

The dc arc characteristics in subsonic orifice nozzle flow

The cold air flow field for a 1/27 cm orifice nozzle was determined for subsonic flow velocities. In addition, d.c. arc voltage, current, and diameter measurements were made for a range of velocities and arc gaps. Average voltage gradients increased rapidly as an arc extinguishing velocity was approached. Measured values of current and diameter were used as an input for relaxation solution of an energy balance equation to compute radial temperature profiles. Calculated arc voltage gradients compare favorably with measured average voltage gradients.

Nagamatsu, H. T.↗
DOE OSTI2015

SU-E-J-273: Skin Temperature Recovery Rate as a Potential Predictor for Radiation-Induced Skin Reactions

Purpose: To assess the potential of dynamic infrared imaging to evaluate early skin reactions during radiation therapy in cancer patients. Methods: Thermal images were captured by our home-built system consisting of two flash lamps and an infrared (IR) camera. The surface temperature of the skin was first raised by ∼ 6 °C from ∼1 ms short flashes; the camera then captured a series of IR images for 10 seconds. For each image series, a basal temperature was recorded for 0.5 seconds before flash was triggered. The temperature gradients (ε) were calculated between a reference point (immediately after the flash) and at a time point of 2sec, 4sec and 9sec after that. A 1.0 cm region of interest (ROI) on the skin was drawn; the mean and standard deviations of the ROIs were calculated. The standard ε values for normal human skins were evaluated by imaging 3 healthy subjects with different skin colors. All of them were imaged on 3 separate days for consistency checks. Results: The temperature gradient, which is the temperature recovery rate, depends on the thermal properties of underlying tissue, i.e. thermal conductivity. The average ε for three volunteers averaged over 3 measurements were 0.64±0.1, 0.72±0.2 and 0.80±0.3 at 2sec, 4sec and 9sec respectively. The standard deviations were within 1.5%–3.2%. One of the volunteers had a prior small skin burn on the left wrist and the ε values for the burned site were around 9% (at 4sec) and 13% (at 9sec) lower than that from the nearby normal skin. Conclusion: The temperature gradients from the healthy subjects were reproducible within 1.5%–3.2 % and that from a burned skin showed a significant difference (9%–13%) from the normal skin. We have an IRB approved protocol to image head and neck patients scheduled for radiation therapy.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2023

Identification of Sea Surface Temperature and Sea Surface Salinity Fronts along the California Coast: Application Using Saildrone and Satellite Derived Products

Coastal upwelling regions are one of the most dynamic areas of the world’s oceans. The California and Baja California Coasts are impacted by both coastal upwelling and the California Current, leading to frontal activity that is captured by gradients in both Sea Surface Temperature (SST) and Sea Surface Salinity (SSS). Satellite data are a great source of spatial data to study fronts. However, biases near coastal areas and coarse resolutions can impair its usefulness in upwelling areas. In this work gradients in SST from NASA Multi-Scale Ultra-High Resolution (MUR) and in two SSS products derived from the Soil Moisture Active Passive (SMAP) NASA mission are compared directly with gradients derived from the Saildrone uncrewed vehicles to validate the gradients as well as to assess their ability to detect known frontal features. The three remotely sensed data sets (MURSST/JPL, SMAP SSS/RSS, SMAP SSS) were co-located with the Saildrone data prior to the calculation of the gradients. Wavelet analysis is used to determine how well the satellite derived SST and SSS products are reproducing the Saildrone derived gradients. Overall results indicate the remote sensing products are reproducing features of known areas of coastal upwelling. Differences between the SST and SSS gradients are mainly associated with the limitations of the microwave derived SSS coverage near land and its reduced spatial resolution. The results are promising for using remote sensing data sets to monitor frontal structure along the California Coast and the application to long term changes in coastal upwelling and dynamics.

Sea Surface Salinity↗
NASA NTRS1978

Paleohydraulics and hydrodynamics of Scabland floods

The last major episode of scabland flooding (approx. 18,000-13,000 years B.P.) left considerable high-water mark evidence in the form of: (1) eroded channel margins; (2) depositional features; (3) ice-rafter erratics; and (4) divide crossings. These were used to reconstruct maximum flood stages and water-surface gradients. Engineering hydraulic calculation procedures allowed the analyses of flood discharges and mean velocities from these data. Secondary flow phenomena, including various forms of vortices and flow separations, are considered to have been the principal erosive processes. The intense pressure and velocity gradients of vortices along the irregular channel boundaries produced the plucking-type erosion.

Baker, V. R.↗
DOE OSTI2016

The structure and electronic properties of hexagonal Fe{sub 2}Si

On the basis of first principle calculations, we show that a hexagonal structure of Fe{sub 2}Si is a ferromagnetic crystal. The result of the phonon spectra indicates that it is a stable structure. Such material exhibits a spin-polarized and half-metal-like band structure. From the calculations of generalized gradient approximation, metallic and semiconducting behaviors are observed with a direct and nearly 0 eV band gap in various spin channels. The densities of states in the vicinity of the Fermi level is mainly contributed from the d-electrons of Fe. We calculate the reflection spectrum of Fe{sub 2}Si, which has minima at 275 nm and 3300 nm with reflectance of 0.27 and 0.49, respectively. Such results may provide a reference for the search of hexagonal Fe{sub 2}Si in experiments. With this band characteristic, the material may be applied in the field of novel spintronics devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Training material models using gradient descent algorithms

High temperature design requires accurate constitutive models to describe material inelastic deformation and failure behavior. Oftentimes, calibrating accurate models devolves into the problem of fitting the model parameters against experimental test data. Here, we present the pyopmat package, an open source framework for calibrating constitutive models against experiment data subjected to various loading conditions using machine learning techniques. The package calculates the exact gradient of the model response with respect to the parameters using a combination of automatic differentiation and the adjoint method. Given this exact gradient, we compare the performance of several gradient-based optimization techniques in fitting realistic constitutive models against data. Here, we demonstrate the efficiency and accuracy of our package through example problems using both synthetic data, generated using known parameter sets, under monotonic and cyclic loading conditions and also with an example applying the techniques developed here to actual high temperature creep-fatigue test data.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Efficient Molecular Dynamics Simulations of Multiple Radical Center Systems Based on the Fragment Molecular Orbital Method

The fully analytic energy gradient has been developed and implemented for the restricted open-shell Hartree–Fock (ROHF) method based on the fragment molecular orbital (FMO) theory for systems that have multiple open-shell molecules. The accuracy of the analytic ROHF energy gradient is compared with the corresponding numerical gradient, illustrating the accuracy of the analytic gradient. The ROHF analytic gradient is used to perform molecular dynamics simulations of an unusual open-shell system, liquid oxygen, and mixtures of oxygen and nitrogen. These molecular dynamics simulations provide some insight about how triplet oxygen molecules interact with each other. Timings reveal that the method can calculate the energy gradient for a system containing 4000 atoms in only 6 h. Therefore, it is concluded that the FMO-ROHF method will be useful for investigating systems with multiple open shells.

36 MATERIALS SCIENCE↗
DOE OSTI2015

An atomic orbital-based formulation of analytical gradients and nonadiabatic coupling vector elements for the state-averaged complete active space self-consistent field method on graphical processing units

We recently presented an algorithm for state-averaged complete active space self-consistent field (SA-CASSCF) orbital optimization that capitalizes on sparsity in the atomic orbital basis set to reduce the scaling of computational effort with respect to molecular size. Here, we extend those algorithms to calculate the analytic gradient and nonadiabatic coupling vectors for SA-CASSCF. Combining the low computational scaling with acceleration from graphical processing units allows us to perform SA-CASSCF geometry optimizations for molecules with more than 1000 atoms. The new approach will make minimal energy conical intersection searches and nonadiabatic dynamics routine for molecular systems with O(10{sup 2}) atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Infrared fixed point of the SU(3) gauge theory with 𝑁 𝑓 = 10 flavors

We use lattice simulations and the continuous renormalization group method, based on the gradient flow, to calculate the 𝛽 function and anomalous dimensions of the SU(3) gauge theory with 𝑁 𝑓 = 10 flavors of fermions in the fundamental representation. We employ several improvements to extend the range of available renormalized couplings, including the addition of heavy Pauli-Villars bosons to reduce cutoff effects and the combination of a range of gradient flow transformations. While in the weak coupling regime our result is consistent with those of earlier studies, our techniques allow us to study the system at much stronger couplings than previously possible. We find that the renormalization group 𝛽 function develops a zero, corresponding to an infrared-stable fixed point, at gradient-flow coupling 𝑔 2 = 15.0⁢(5). We also determine the mass and tensor anomalous dimensions: At the fixed point we find 𝛾 𝑚 ≃0.6, suggesting that this system might be deep inside the conformal window.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Electronic transport, thermal transport, thermal expansion, and magnetization in the strongly correlated metal LaNi⁢O 3

Perovskite structured LaNiO 3 is a strongly correlated metal with intriguing thermal and magnetic properties. The volume dependence of calculated and measured physical properties can add additional critical information to develop a more in-depth understanding of this strongly correlated phenomenon. Taking advantage of recent single crystal LaNiO 3 growth using the floating-zone method, we have measured the thermal expansion, the magnetostriction, and the pressure dependence of the magnetic susceptibility, which then allows derivation of the Grüneisen parameters γ e = $\frac{dlnN(E_F)}{d lnV}$, γ χ = $\frac{dlnχ}{d lnV}$, as well as of electric and thermal transport properties. We simulate the volume dependence of structural and magnetic properties using Density Functional Theory calculations at the Generalized Gradient Approximation level. A large discrepancy between experimental values and calculated ones suggests that strong correlations are likely to be dynamic in nature. This study also provides a side-by-side comparison of measurements in single crystal and polycrystalline samples of LaNiO 3 to elucidate intrinsic materials properties. A broad hump at high temperatures in the temperature dependence of magnetization found in the single crystal sample of LaNiO 3 has been rationalized by a model that includes the influence of electron correlations on the Landau diamagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1995

Optical diagnostics of solution crystal growth

Solution crystal growth monitoring of LAP/TGS crystals by various optical diagnostics systems, such as conventional and Mach-Zehnder (M-Z) interferometers, optical heterodyne technique, and ellipsometry, is under development. The study of the dynamics of the crystal growth process requires a detailed knowledge of crystal growth rate and the concentration gradient near growing crystals in aqueous solution. Crystal growth rate can be measured using conventional interferometry. Laser beam reflections from the crystal front as well as the back surface interfere with each other, and the fringe shift due to the growing crystal yields information about the growth rate. Our preliminary results indicate a growth rate of 6 A/sec for LAP crystals grown from solution. Single wavelength M-Z interferometry is in use to calculate the concentration gradient near the crystal. Preliminary investigation is in progress using an M-Z interferometer with 2 cm beam diameter to cover the front region of the growing crystal. In the optical heterodyne technique, phase difference between two rf signals (250 KHZ) is measured of which one is a reference signal, and the other growth signal, whose phase changes due to a change in path length as the material grows. From the phase difference the growth rate can also be calculated. Our preliminary results indicate a growth rate of 1.5 A/sec. the seed and solution temperatures were 26.46 C and 27.92 C respectively, and the solution was saturated at 29.0 C. an ellipsometer to measure the growth rate and interface layer is on order from JOBIN YVON, France. All these systems are arranged in such a manner that measurements can be made either sequentially or simultaneously. These techniques will be adapted for flight experiment.

Kim, Yongkee↗
DOE OSTI2015

First-principles study on the bulk and (1 1 1) surface half-metallicity of KS and RbS in CsCl structure

Graphical abstract: Spin-polarized total and atomic DOS at S-(1 1 1) terminated slab and bulk in CsCl-type RbS. - Highlights: • The half metallic properties of CsCl-type RbS and KS have been studied. • The RbS's and KS's (1 1 1) slabs have been investigated. • Surface energy of RbS's and KS's (1 1 1) slabs are calculated. - Abstract: The electronic and magnetic properties of RbS and KS in CsCl structure have been investigated by using the full-potential local-orbital minimum-basis method. Calculating the relation between the total energies and lattice parameters for RbS and KS, we find out that the equilibrium lattice parameters are 4.02 Å and 3.84 Å for RbS and KS, respectively. According to our calculations in generalized gradient approximation approximation, both RbS and KS are half-metallic ferromagnets with the magnetic moments of 1 μ{sub B} per formula unit, and band gap of 4.287 eV for RbS and 4.395 eV for KS. We also have studied the electronic and magnetic properties of (1 1 1) surfaces of RbS and KS, and have found out that the half-metallicity of their bulk is preserved in all of those surfaces. Finally, through the calculations of formation energy of RbS and KS, it is found that their thin films are stable in the equilibrium conditions, and the Rb-terminated (1 1 1) slab of RbS and the K-terminated (1 1 1) slab of KS are more stable than their S-terminated (1 1 1) slabs. All of the above properties lead the compounds of RbS and KS in CsCl structure to be promising candidates for spintronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Classical-Quantum Algorithm for Solving Stochastic Programs

Stochastic programming provides a rigorous mathematical framework for making decisions under uncertainty in a risk-aware manner. Two-stage stochastic programming is, perhaps, the simplest form of this framework. Here the first-stage variables represent decisions that must be made "here and now" in the face of uncertainty, while the second-stage variables are decisions made after uncertain events. However, the broad adoption of stochastic programming has been hindered by computational challenges caused by the two-stage stochastic programming formulation which requires solving an ensemble of optimization problems. Using quantum amplitude estimation (QAE), quantum computers have shown the theoretic ability to compute expectations with Monte-Carlo methods with quadratically fewer samples than classical methods. In this work, we present a quantum algorithm for computing the expectation term using QAE for given first-stage decisions. Further, we detail methods of computing gradient information from the quantum calculation enabling the application of classical gradient-based optimization techniques. The result is a classical-quantum hybrid method of solving two-stage stochastic programs. These techniques are demonstrated with computational experiments based an engineering optimization problem.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Parallel algorithm for dominant points correspondences in robot binocular stereo vision

This paper presents an algorithm to find the correspondences of points representing dominant feature in robot stereo vision. The algorithm consists of two main steps: dominant point extraction and dominant point matching. In the feature extraction phase, the algorithm utilizes the widely used Moravec Interest Operator and two other operators: the Prewitt Operator and a new operator called Gradient Angle Variance Operator. The Interest Operator in the Moravec algorithm was used to exclude featureless areas and simple edges which are oriented in the vertical, horizontal, and two diagonals. It was incorrectly detecting points on edges which are not on the four main directions (vertical, horizontal, and two diagonals). The new algorithm uses the Prewitt operator to exclude featureless areas, so that the Interest Operator is applied only on the edges to exclude simple edges and to leave interesting points. This modification speeds-up the extraction process by approximately 5 times. The Gradient Angle Variance (GAV), an operator which calculates the variance of the gradient angle in a window around the point under concern, is then applied on the interesting points to exclude the redundant ones and leave the actual dominant ones. The matching phase is performed after the extraction of the dominant points in both stereo images. The matching starts with dominant points in the left image and does a local search, looking for corresponding dominant points in the right image. The search is geometrically constrained the epipolar line of the parallel-axes stereo geometry and the maximum disparity of the application environment. If one dominant point in the right image lies in the search areas, then it is the corresponding point of the reference dominant point in the left image. A parameter provided by the GAV is thresholded and used as a rough similarity measure to select the corresponding dominant point if there is more than one point the search area. The correlation is used as a final decision tool when there is still more than one point in the search area. If there is no dominant point in the search area of if the points in the search area are below a correlation threshold, then the dominant point in the reference image is occluded and can not be corresponded. The algorithm has been modeled, implemented and shown to be fast, robust and parallel. The parallelism is created from three main features: locality of the operators; a memory optimization scheme; and the ability to fully parallelize the extraction phase which is the most computational intensive task in the algorithm. The last feature is achieved by performing the extraction phase on the two images simultaneously.

Al-Tammami, A.↗
NASA NTRS1991

A 3-D upwind Euler solver for unstructured meshes

A three-dimensional finite-volume upwind Euler solver is developed for unstructured meshes. The finite-volume scheme solves for solution variables at vertices of the mesh and satisfies the integral conservation law on nonoverlapping polyhedral control volumes surrounding vertices of the mesh. The schene achieves improved solution accuracy by assuming a piecewise linear variation of the solution in each control volume. This improved spatial accuracy hinges heavily upon the calculation of the solution gradient in each control volume given pointwise values of the solution at vertices of the mesh. Several algorithms are discussed for obtaining these gradients. Details concerning implementation procedures and data structures are discussed. Sample calculations for inviscid Euler flow about isolated aircraft wings at subsonic and transonic speeds are compared with established Euler solvers as well as experiment.

Barth, Timothy J.↗
NASA NTRS2012

Adaptive Aft Signature Shaping of a Low-Boom Supersonic Aircraft Using Off-Body Pressures

The design and optimization of a low-boom supersonic aircraft using the state-of-the- art o -body aerodynamics and sonic boom analysis has long been a challenging problem. The focus of this paper is to demonstrate an e ective geometry parameterization scheme and a numerical optimization approach for the aft shaping of a low-boom supersonic aircraft using o -body pressure calculations. A gradient-based numerical optimization algorithm that models the objective and constraints as response surface equations is used to drive the aft ground signature toward a ramp shape. The design objective is the minimization of the variation between the ground signature and the target signature subject to several geometric and signature constraints. The target signature is computed by using a least-squares regression of the aft portion of the ground signature. The parameterization and the deformation of the geometry is performed with a NASA in- house shaping tool. The optimization algorithm uses the shaping tool to drive the geometric deformation of a horizontal tail with a parameterization scheme that consists of seven camber design variables and an additional design variable that describes the spanwise location of the midspan section. The demonstration cases show that numerical optimization using the state-of-the-art o -body aerodynamic calculations is not only feasible and repeatable but also allows the exploration of complex design spaces for which a knowledge-based design method becomes less effective.

Ordaz, Irian↗