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 811 records · Page 45

Determination of thermospheric quantities from simple ionospheric observations using numerical simulation

Measured ionospheric electron content and peak electron concentration data are introduced into a numerical simulation of the ionosphere to yield values of induced plasma drifts and exospheric neutral temperatures consistent with the observations. Data collected on 23-24 March 1970 on the East Coast of the U.S.A. are analyzed and the results are in agreement with incoherent radar measurements at Millstone Hill, Massachusetts. Neutral winds and meridional exospheric temperature gradients that give rise to the computed plasma drifts are calculated through the use of a dynamic model of the thermosphere.

Antoniadis, D. A.↗

The calculated and observed ionospheric properties during Atmospheric Explorer-C satellite crossings over Millstone Hill

Results are presented for calculated and observed incoherent scatter radar data on ionospheric and atmospheric properties during the crossings of the Atmospheric Explorer-C satellite over the Millstone Hill, Massachusetts incoherent scatter radar station. The incoherent scatter radar measured electron density vertical profiles and electron and ion temperatures. These values yielded the exospheric temperature and the vertical distribution of neutral gas temperature. The satellite measured electron and ion temperatures and densities, neutral gas composition, and photoelectron spectra. Measurements were also made of solar UV and EUV fluxes, the vertical profile of the NO number density, and the vertical profile of the 6300 A airglow-volume emission rate. The results are compared to those of a time-dependent coupled model of the ionospheric E- and F-regions. Good agreement is found between satellite results, incoherent scatter radar measurements, and model calculations. Along the orbital path satellite measurements show significant latitudinal gradients in the measured properties.

Roble, R. G.↗

Airfoil design by numerical optimization using a minicomputer

A computer program developed for the automated design of low speed airfoils utilizes a generalized Joukowski method for aerodynamic analysis coupled with a conjugate gradient, penalty function, numerical optimization algorithm to give an efficient calculation technique for use with minicomputers. The program designs airfoils with a prescribed pressure distribution as well as those which minimize or maximize some aerodynamic force coefficient. At present the method is restricted to inviscid, incompressible flow. A typical design problem will execute in 4.5 hr on an HP 9830 minicomputer.

Hicks, R. M.↗

A three-dimensional model of corotating streams in the solar wind. II - Hydrodynamic streams

Theoretical aspects of corotating solar wind dynamics on a global scale are examined by means of numerical simulations executed with a nonlinear, inviscid, adiabatic, single fluid, three-dimensional hydrodynamic formulation. The factors most strongly affecting stream evolution, namely, the intrinsic correlations among density, temperature, and velocity existing near the source, the amplitude, longitudinal breadth and latitudinal breadth of the stream, and the heliographic latitude of the centroid of the stream, are considered. Calculations indicate that even in the presence of substantial meridional gradients, the latitudinal transport of mass, energy, and momentum by streams appears to have little impact on the overall evolution of the structure.

Pizzo, V. J.↗

The absence of convective stability in incongruent crystal growth on Earth

It is often assumed that convection can be excluded in materials processing by avoiding horizontal temperature gradients, and employing top-heating and horizontal interfaces. Model calculations for vapor transport processes show that this supposition holds only for strictly congruent vapor-solid transitions. This result can be generalized to all fluid-solid phase transitions. Due to the omnipresence of impurity segregation, it is concluded that convective stability is practically unattainable in materials processing on Earth.

Rosenberger, F. E.↗

Advanced high temperature static strain sensor development

An examination was made into various techniques to be used to measure static strain in gas turbine liners at temperatures up to 1150 K (1600 F). The methods evaluated included thin film and wire resistive devices, optical fibers, surface acoustic waves, the laser speckle technique with a heterodyne readout, optical surface image and reflective approaches and capacitive devices. A preliminary experimental program to develop a thin film capacitive device was dropped because calculations showed that it would be too sensitive to thermal gradients. In a final evaluation program, the laser speckle technique appeared to work well up to 1150 K when it was used through a relatively stagnant air path. The surface guided acoustic wave approach appeared to be interesting but to require too much development effort for the funds available. Efforts to develop a FeCrAl resistive strain gage system were only partially successful and this part of the effort was finally reduced to a characterization study of the properties of the 25 micron diameter FeCrAl (Kanthal A-1) wire. It was concluded that this particular alloy was not suitable for use as the resistive element in a strain gage above about 1000 K.

Hulse, C. O.↗

The effect of anisotropic crystal-melt surface tension on grain boundary groove morphology

The shape of a stationary solid-liquid interface in a temperature gradient near a grain boundary in a pure material is calculated for anisotropic crystal-melt surface tension and equal thermal conductivities of crystal and melt. Results are compared with those for the well-known problem of the two-dimensional equilibrium shape of a crystal. For small anisotropy, the resulting interface shapes have continuously turning tangents but differ in detail from the grain boundary groove shapes that have been calculated for isotropic surface tension. For larger anisotropy, the interface shapes have discontinuities in slope as a result of missing orientations.

Voorhees, P. W.↗

High resolution flow field prediction for tail rotor aeroacoustics

The prediction of tail rotor noise due to the impingement of the main rotor wake poses a significant challenge to current analysis methods in rotorcraft aeroacoustics. This paper describes the development of a new treatment of the tail rotor aerodynamic environment that permits highly accurate resolution of the incident flow field with modest computational effort relative to alternative models. The new approach incorporates an advanced full-span free wake model of the main rotor in a scheme which reconstructs high-resolution flow solutions from preliminary, computationally inexpensive simulations with coarse resolution. The heart of the approach is a novel method for using local velocity correction terms to capture the steep velocity gradients characteristic of the vortex-dominated incident flow. Sample calculations have been undertaken to examine the principal types of interactions between the tail rotor and the main rotor wake and to examine the performance of the new method. The results of these sample problems confirm the success of this approach in capturing the high-resolution flows necessary for analysis of rotor-wake/rotor interactions with dramatically reduced computational cost. Computations of radiated sound are also carried out that explore the role of various portions of the main rotor wake in generating tail rotor noise.

Quackenbush, Todd R.↗

Estimation of tropical forest canopy temperatures, thermal response numbers, and evapotranspiration using an aircraft-based thermal sensor

Thermal infrared Multispectral Scanner (TIMS) data were collected at a resolution of 5 to 10 m from a tropical rain forest over an elevation gradient from 35 to 2700 m in the Braulio Carrillo National Park in Costa Rica. Flight lines were repeated with a 15 to 30 minute time difference for measurement of forest canopy thermal response over time. Concurrent radiosonde measurements of atmospheric profiles of air temperature and moisture provided inputs to LOWTRAN6 for atmospheric radiance corrections of the TIMS data. Techniques for using calibrated aircraft-based thermal scanner data to examine tropical forest canopy thermal properties are described. Forest canopy temperature changes over time assessed between repeated, duplicated flight lines were combined with estimates of surface radiative energy measurements from towers above the forest canopy to determine temperature spatial variability, calculate Thermal Response Numbers (TRN), and estimate evapotranspiration along the elevation gradient from selected one hectare forest inventory plots.

Luvall, Jeffrey C.↗

Investigations on entropy layer along hypersonic hyperboloids using a defect boundary layer

A defect approach coupled with matched asymptotic expansions is used to derive a new set of boundary layer equations. This method ensures a smooth matching of the boundary layer with the inviscid solution. These equations are solved to calculate boundary layers over hypersonic blunt bodies involving the entropy gradient effect. Systematic comparisons are made for both axisymmetric and plane flows in several cases with different Mach and Reynolds numbers. After a brief survey of the entropy layer characteristics, the defect boundary layer results are compared with standard boundary layer and full Navier-Stokes solutions. The entropy gradient effects are found to be more important in the axisymmetric case than in the plane one. The wall temperature has a great influence on the results through the displacement effect. Good predictions can be obtained with the defect approach over a cold wall in the nose region, with a first order solution. However, the defect approach gives less accurate results far from the nose on axisymmetric bodies because of the thinning of the entropy layer.

Brazier, J. P.↗

Line profiles variations from atmospheric eclipses: Constraints on the wind structure in Wolf-Rayet stars

Binary systems in which one of the components has a stellar wind may present a phenomenon known as 'wind' or 'atmospheric eclipse', in which that wind occults the luminous disk of the companion. The enhanced absorption profile, relative to the spectrum at uneclipsed orbital phases, can be be modeled to yield constraints on the spatial structure of the eclipsing wind. A new, very efficient approach to the radiative transfer problem, which makes no requirements with respect to monotonicity of the velocity gradient or size of that gradient, is presented. The technique recovers both the comoving frame calculation and the Sobolev approximation in the appropiate limits. Sample computer simulations of the line profile variations induced by wind eclipses are presented. It is shown that the location of the wind absorption features in frequency is a diagnostic tool for identifying the size of the wind acceleration region. Comparison of the model profile variations with the observed variations in the Wolf-Rayet (W-R)+6 binary system V444 Cyg illustrate how the method can be used to derive information on the structure of the wind of the W-R star constrain the size of the W-R core radius.

Auer, L. H.↗

Neural Network and Regression Approximations in High Speed Civil Transport Aircraft Design Optimization

Nonlinear mathematical-programming-based design optimization can be an elegant method. However, the calculations required to generate the merit function, constraints, and their gradients, which are frequently required, can make the process computational intensive. The computational burden can be greatly reduced by using approximating analyzers derived from an original analyzer utilizing neural networks and linear regression methods. The experience gained from using both of these approximation methods in the design optimization of a high speed civil transport aircraft is the subject of this paper. The Langley Research Center's Flight Optimization System was selected for the aircraft analysis. This software was exercised to generate a set of training data with which a neural network and a regression method were trained, thereby producing the two approximating analyzers. The derived analyzers were coupled to the Lewis Research Center's CometBoards test bed to provide the optimization capability. With the combined software, both approximation methods were examined for use in aircraft design optimization, and both performed satisfactorily. The CPU time for solution of the problem, which had been measured in hours, was reduced to minutes with the neural network approximation and to seconds with the regression method. Instability encountered in the aircraft analysis software at certain design points was also eliminated. On the other hand, there were costs and difficulties associated with training the approximating analyzers. The CPU time required to generate the input-output pairs and to train the approximating analyzers was seven times that required for solution of the problem.

Patniak, Surya N.↗

Jupiter's Atmospheric Temperatures: From Voyager IRIS to Cassini CIRS

Retrievals run on Cassini Composite Infrared Spectrometer data obtained during the distant Jupiter flyby have been used to generate global temperature maps of the planet in the troposphere and stratosphere. Similar retrievals were performed on Voyager 1 IRIS data and have provided the first detailed IRIS map of the stratosphere. In both data sets, high latitude troposphere temperatures are presented for the first time, and the meridional gradients indicate the presence of circumpolar jets. Thermal winds were calculated for each data set and show strong vertical shears in the zonal winds at low latitudes. The temperatures retrieved from the two spacecraft were also compared with yearly ground-based data obtained over the intervening two decades. Tropospheric temperatures reveal gradual changes at low latitudes, with little obvious seasonal or short-term variation (Orton et al. 1994). Stratospheric temperatures show much more complicated behavior over short timescales, consistent with quasi-quadrennial oscillations at low latitudes, as suggested in prior analyses of shorter intervals of ground- based data (Orton et al. 1991, Friedson 1999). A scaling analysis indicates that meridional motions, mechanically forced by wave or eddy convergence, play an important role in modulating the temperatures and winds in the upper troposphere and stratosphere on seasonal and shorter time scales. At latitudes away from the equator, the mechanical forcing can be derived simply from a temporal record of temperature and its vertical derivative. Ground-based observations with improved vertical resolution and/or long-term monitoring from spacecraft are required for this purpose.

Simon-Miller, Amy A.↗

Energetic and spatial parameters for gating of the bacterial large conductance mechanosensitive channel, MscL

MscL is multimeric protein that forms a large conductance mechanosensitive channel in the inner membrane of Escherichia coli. Since MscL is gated by tension transmitted through the lipid bilayer, we have been able to measure its gating parameters as a function of absolute tension. Using purified MscL reconstituted in liposomes, we recorded single channel currents and varied the pressure gradient (P) to vary the tension (T). The tension was calculated from P and the radius of curvature was obtained using video microscopy of the patch. The probability of being open (Po) has a steep sigmoidal dependence on T, with a midpoint (T1/2) of 11.8 dyn/cm. The maximal slope sensitivity of Po/Pc was 0.63 dyn/cm per e-fold. Assuming a Boltzmann distribution, the energy difference between the closed and fully open states in the unstressed membrane was DeltaE = 18.6 kBT. If the mechanosensitivity arises from tension acting on a change of in-plane area (DeltaA), the free energy, TDeltaA, would correspond to DeltaA = 6.5 nm2. MscL is not a binary channel, but has four conducting states and a closed state. Most transition rates are independent of tension, but the rate-limiting step to opening is the transition between the closed state and the lowest conductance substate. This transition thus involves the greatest DeltaA. When summed over all transitions, the in-plane area change from closed to fully open was 6 nm2, agreeing with the value obtained in the two-state analysis. Assuming a cylindrical channel, the dimensions of the (fully open) pore were comparable to DeltaA. Thus, the tension dependence of channel gating is primarily one of increasing the external channel area to accommodate the pore of the smallest conducting state. The higher conducting states appear to involve conformational changes internal to the channel that don't involve changes in area.

Non-NASA Center↗

TU-H-206-03: Characterizing B1 Inhomogeneities in DCE MRI

Purpose: Dynamic Contrast Enhanced (DCE) MRI is a valuable technique for measuring perfusion and permeability characteristics of tumors. Exogenous contrast concentrations are calculated based on changes in T{sub 1} measured using fast 3D gradient echo (FLASH) sequences. However, the slab selective pulses used in 3D MRI may result in B{sub 1} inhomogeneities across the volume of interest that can lead to errors in T{sub 1} and thus the estimated gadolinium concentration. We compared three FLASH DCE sequences (GRE, TWIST, and VIBE) to determine their signal homogeneity across slices and the accuracy in calculating T{sub 1} using acquisitions with variable flip angles. Methods: The sequences were tested at 3 T on a Siemens mMR (VB20P) using a doped water phantom 3.75 g/L NiSO{sub 4} - 6H{sub 2}O + 5 g/L NaCl (T{sub 1} = 104 ms) and a 2% agar, 0.67% NaCl phantom (T{sub 1}= 1.71 s). 2D EPI B{sub 1} maps and inversion recovery T{sub 1}maps were acquired for ground truth. 3D MRI was acquired at different flip angles to generate a T{sub 1} map. Regions of interest were drawn to measure signal inside the phantoms as a function of slice position. The T{sub 1} for each slice ROI was fit to the FLASH steady-state model of magnetization. Results: Based on the data, GRE gave the most uniform signal homogeneity and T{sub 1} values in the middle slices of the 3D volume. The 3D VIBE sequence had the largest region of signal inhomogeneity compared to the 3D GRE and TWIST sequences. VIBE’s B{sub 1} inhomogeneity is inconsistent at low flip angles. However, VIBE resulted in more slices with T{sub 1} values similar to the ground truth. Conclusion: The central 1/3 of the slices yielded signals that result in T{sub 1} fits consistent with the ground truth. However, the remaining slices required some form of B{sub 1} inhomogeneity correction for quantitative DCE analysis. The research was supported in part by NIH NCI Grant R01CA159471.

60 APPLIED LIFE SCIENCES↗

Towards an Optimal Gradient-dependent Energy Functional of the PZ-SIC Form

Results of Perdew–Zunger self-interaction corrected (PZ-SIC) density functional theory calculations of the atomization energy of 35 molecules are compared to those of high-level quantum chemistry calculations. While the PBE functional, which is commonly used in calculations of condensed matter, is known to predict on average too high atomization energy (overbinding of the molecules), the application of PZ-SIC gives a large overcorrection and leads to significant underestimation of the atomization energy. The exchange enhancement factor that is optimal for the generalized gradient approximation within the Kohn-Sham (KS) approach may not be optimal for the self-interaction corrected functional. The PBEsol functional, where the exchange enhancement factor was optimized for solids, gives poor results for molecules in KS but turns out to work better than PBE in PZ-SIC calculations. The exchange enhancement is weaker in PBEsol and the functional is closer to the local density approximation. Furthermore, the drop in the exchange enhancement factor for increasing reduced gradient in the PW91 functional gives more accurate results than the plateaued enhancement in the PBE functional. A step towards an optimal exchange enhancement factor for a gradient dependent functional of the PZ-SIC form is taken by constructing an exchange enhancement factor that mimics PBEsol for small values of the reduced gradient, and PW91 for large values. The average atomization energy is then in closer agreement with the high-level quantum chemistry calculations, but the variance is still large, the F 2 molecule being a notable outlier.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Kohn-Sham Solver (KSSOLV) v2.0

KSSOLV is a MATLAB toolbox for solving Kohn-Sham density functional theory based electronic structure eigenvalue problems. It uses an object oriented features of MATLAB to represent atom, molecules, wavefunctions and Hamiltonians and their operations. It is designed to make it easier for users to prototype and test new algorithms for solving the Kohn-Sham problem. KSSOLV2.0 contains significant improvement over the original KSSOLV described in a paper published in ACM Transaction on Mathematical Software (attached). In addition to performing ground state calculation for small molecules, it can also perform geometry optimization for both molecules and solids. It uses standard pseudopotentials and implements local density approximation, generalized gradient approximation and hybrid functionals. Future releases will also include time-dependent DFT and post DFT calculations such as the GW quasi-particle energy calculation and Bethe-Salpeter equation solver for optical absorption.

Yang, Chao↗

Colossal magnetoresistance in a nonsymmorphic antiferromagnetic insulator

Here we investigate antiferromagnetic Eu 5 In 2 Sb 6 , a nonsymmorphic Zintl phase. Our electrical transport data show that Eu 5 In 2 Sb 6 is remarkably insulating and exhibits an exceptionally large negative magnetoresistance, which is consistent with the presence of magnetic polarons. From ab initio calculations, the paramagnetic state of Eu 5 In 2 Sb 6 is a topologically nontrivial semimetal within the generalized gradient approximation (GGA), whereas an insulating state with trivial topological indices is obtained using a modified Becke–Johnson potential. Notably, GGA + U calculations suggest that the antiferromagnetic phase of Eu 5 In 2 Sb 6 may host an axion insulating state. Our results provide important feedback for theories of topological classification and highlight the potential of realizing clean magnetic narrow-gap semiconductors in Zintl materials.

36 MATERIALS SCIENCE↗