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.

SourceResearch connectionFrom yearThrough year

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

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Understanding the large-scale structure of magnetic fields in the corona and heliosphere

To understand the non-radial nature of the southwestern streamer in the 3 November 1994 eclipse and the near absence of a latitudinal gradient of the interplanetary magnetic field strength seen by Ulysses over the southern pole, we calculate the coronal and inter planetary magnetic field from WSO photospheric field data using both the horizontal current-current sheet model (Zhao and Hoeksema, 1994) and the current sheet-source surface model (Zhao and Hoeksema, 1995). These models account for the effects of large-currents and current sheets in the corona. Comparison of the calculated results with white light, X-ray, and IMF observations show that non-radial coronal streamers and a near-zero latitudinal gradient of the radial component of the IMF may be explained by a large-scale current sheet. To best match the observations, the height of the cusp point (above which current sheets form) differs from one streamer to another. While heliospheric volume currents flowing above the source surface certainly affect the interplanetary field, they have little or no influence on the structure of the corona.

Zhao, X. P.↗

NASA NTRS1990

Spectra of circularly polarized radiation from astrophysical OH masers

A striking feature of astrophysical masers is the tendency for either one or the other of the circular polarizations to dominate in the radiation from the strong, widely observed masing transitions of OH at 18 cm. Spectral line profiles are calculated for polarized maser radiation due to the combined effects of a velocity gradient and, as is indicated for these transitions, a Zeeman splitting that is at least comparable with the thermal contributions to the breadths of the spectral lines. The resulting spectral features are similar in appearance, including the presence of large net circular polarization and narrow line breadths, to the commonly observed spectra of OH masers in molecular clouds. The calculations presented here are performed as a function of frequency without making the approximations of a large velocity gradient. Rapid cross relaxation, which has been advocated by others for the OH masers, is assumed.

Nedoluha, Gerald E.↗

NASA NTRS1973

Unsteady turbulent boundary layers in two-dimensional, incompressible flow.

The paper presents the results of some calculations aimed at studying the characteristics of turbulent boundary layers in unsteady flow. The calculations were performed using the method of Singleton and Nash, which is based on the turbulent kinetic energy equation. Emphasis is placed on the approach to zero wall stress, and the results show the extent to which onset of this condition can be delayed by the effects of time dependence. The delay is not associated merely with alleviation of the pressure gradients, and quasi-steady methods cannot readily be used to calculate this type of flow.

Nash, J. F.↗

DOE OSTI2021

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Thermal energy and charge currents in multi-terminal nanorings

We study in experiment and theory thermal energy and charge transfer close to the quantum limit in a ballistic nanodevice, consisting of multiply connected one-dimensional electron waveguides. The fabricated device is based on an AlGaAs/GaAs heterostructure and is covered by a global top-gate to steer the thermal energy and charge transfer in the presence of a temperature gradient, which is established by a heating current. The estimate of the heat transfer by means of thermal noise measurements shows the device acting as a switch for charge and thermal energy transfer. The wave-packet simulations are based on the multi-terminal Landauer-Büttiker approach and confirm the experimental finding of a mode-dependent redistribution of the thermal energy current, if a scatterer breaks the device symmetry.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

NASA NTRS1998

Operational Data Reduction Procedure for Determining Density and Vertical Structure of the Martian Upper Atmosphere from Mars Global Surveyor Accelerometer Measurements

The success of aerobraking by the Mars Global Surveyor (MGS) spacecraft was partly due to the analysis of MGS accelerometer data. Accelerometer data was used to determine the effect of the atmosphere on each orbit, to characterize the nature of the atmosphere, and to predict the atmosphere for future orbits. To interpret the accelerometer data, a data reduction procedure was developed to produce density estimations utilizing inputs from the spacecraft, the Navigation Team, and pre-mission aerothermodynamic studies. This data reduction procedure was based on the calculation of aerodynamic forces from the accelerometer data by considering acceleration due to gravity gradient, solar pressure, angular motion of the MGS, instrument bias, thruster activity, and a vibration component due to the motion of the damaged solar array. Methods were developed to calculate all of the acceleration components including a 4 degree of freedom dynamics model used to gain a greater understanding of the damaged solar array. The total error inherent to the data reduction procedure was calculated as a function of altitude and density considering contributions from ephemeris errors, errors in force coefficient, and instrument errors due to bias and digitization. Comparing the results from this procedure to the data of other MGS Teams has demonstrated that this procedure can quickly and accurately describe the density and vertical structure of the Martian upper atmosphere.

Cancro, George J.↗

DOE OSTI2019

Computation of magnetic fields from field components on a plane grid

An algorithm is presented to calculate the field components of the magnetic field [B x (x, y, z), B y (x, y, z), B z (x, y, z)] at a point (x, y, z) in space, from the knowledge of the components [B x (x, y) = 0, z), B y (x, y = 0, z), B z (x, y = 0,z)] on a “reference plane”, which is normal to the y-axis at . The algorithm, which is a general one and is not restricted to fields with mid-plane symmetry is based on the Maclaurin series expansion of the magnetic field components at any point in space in terms of the distance (y) of the point from the reference plane. The coefficients of the Maclaurin series expansion are expressed in terms of the on-plane field components and their partial derivatives with respect to spatial coordinates (). The field components are usually generated from magnetic field measurements on a rectangular grid on the plane. This algorithm was employed in 1986 in the RAYTRACE computer code to help calculate the optical properties of magnets and of the Alternating Gradient Synchrotron (AGS) at the Brookhaven National Laboratory (BNL). Here, a general mathematical formulation of this algorithm based on the differential algebraic method was presented by Makino in 2011. This paper presents the step by step derivation of the algorithm and provides the necessary formulas to be introduced by the reader in any computer code which requires the field components generated by magnetic devices. In addition provides an example of the use of the algorithm and its limitations as applied to a Halbach type magnet with or without median plane symmetry.

43 PARTICLE ACCELERATORS↗

NASA NTRS1990

Application of implicit numerical techniques to the solution of the three-dimensional diffusion equation

Implicit techniques for calculating three-dimensional, time-dependent heat diffusion in a cube are tested with emphasis on storage efficiency, accuracy, and speed of calculation. For this purpose, a tensor product technique with both Chebyshev collocation and finite differences and a generalized conjugate gradient technique with finite differences are used in conjunction with Crank-Nicolson discretization. An Euler explicit finite difference calculation is performed for use as a benchmark. The implicit techniques are found to be competitive with the Euler explicit method in terms of storage efficiency and speed of calculation and offer advantages both in accuracy and stability. Mesh stretching in the finite difference calculations is shown to markedly improve the accuracy of the solution.

Peltier, Leonard Joel↗

DOE OSTI2018

Calculation of Permeability for Porous Media with Platy Grains - 18095

The permeability is calculated numerically for a porous medium with granules of platy shape when pressure gradient is imposed externally across the medium by considering slow flow of a viscous fluid in the pores. Under the assumption of periodicity a unit cell which contains solid phase of platy granules is selected on the microscale and the flow field is calculated by solving Stokes problem in which inertial effects are ignored. A Boundary-Vale Problem is obtained in the course of applying the theory of homogenization to fluid flow through porous media. It is shown that the permeability is strongly influenced by the flow direction along which the pressure gradient on the macro-scale is imposed and the particle arrangement pattern in the unit cell used in the numerical computation. (authors)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Analytical nonadiabatic coupling and state-specific energy gradient for the crystal field Hamiltonian describing lanthanide single-ion magnets

Paramagnetic molecules with a metal ion as an electron spin center are promising building blocks for molecular qubits and high-density memory arrays. However, fast spin relaxation and decoherence in these molecules lead to a rapid loss of magnetization and quantum information. Nonadiabatic coupling (NAC), closely related to spin-vibrational coupling, is the main source of spin relaxation and decoherence in paramagnetic molecules at higher temperatures. Predicting these couplings using numerical differentiation requires a large number of computationally intensive ab initio or crystal field electronic structure calculations. To reduce computational cost and improve accuracy, we derive and implement analytical NAC and state-specific energy gradient for the ab initio parametrized crystal field Hamiltonian describing single-ion molecular magnets. Our implementation requires only a single crystal field calculation. In addition, the accurate NACs and state-specific energy gradients can be used to model spin relaxation using sophisticated nonadiabatic molecular dynamics, which avoids the harmonic approximation for molecular vibrations. To test our implementation, we calculate the NAC values for three lanthanide complexes. Finally, the predicted values support the relaxation mechanisms reported in previous studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1988

Internal temperature gradient of alloy semiconductor melts from interrupted growth experiments

Interrupted growth experiments on Hg(0.8)Cd(0.2)Se alloys were performed in an attempt to estimate the actual internal temperature gradient in the melt during directional solidification. The results have been analyzed in terms of a theoretical model which assumes an axial composition profile characteristic of diffusion controlled solute redistribution during growth. A comparison of the calculated and measured values suggests a reduction in the applied (empty furnace) temperature gradient by about a factor of three. Values of the interface segregation coefficient (k) determined from the phase diagram, predict values of solute concentration build up in the solid during growth interruptions which are inconsistent with measured results. It appears that the k values used tend to underestimate the actual k values for higher alloy compositions and overestimate k for the lower alloy compositions.

Andrews, R. N.↗

NASA NTRS1976

Interplanetary stream magnetism - Kinematic effects

The particle density and the magnetic-field intensity and direction are calculated for volume elements of the solar wind as a function of the initial magnetic-field direction and the initial speed gradient. It is assumed that the velocity is constant and radial. These assumptions are approximately valid between about 0.1 and 1.0 AU for many streams. Time profiles of the particle density, field intensity, and velocity are calculated for corotating streams, neglecting effects of pressure gradients. The compression and rarefaction of the magnetic field depend sensitively on the initial field direction. By averaging over a typical stream, it is found that the average radial field intensity is inversely proportional to the square of the heliocentric distance, whereas the average intensity in the direction of the planets' motion does not vary in a simple way, consistent with deep space observations. Changes of field direction may be very large, depending on the initial angle; but when the initial angle at 0.1 AU is such that the base of the field line corotates with the sun, the spiral angle is the preferred direction at 1 AU. The theory is also applicable to nonstationary flows.

Burlaga, L. F.↗

NASA NTRS1993

Effects of thermal-solutal convection on temperature and solutal fields under various gravitational orientations

Semiconductor crystals such as Hg(1-x)CD(x)Te grown by unidirectional solidification Bridgmann method have shown compositional segregations in both the axial and radial directions. Due to the wide separation between the liquidus and the solidus of its pseudobinary phase diagram, there is a diffusion layer of higher HgTe content built up in the melt near the melt-solid interface which gives a solute concentration gradient in the axial direction. The value of effective diffusion coefficient calculated from fitting of the data to 1D model varies with Hg(1-x)Cd(x)Te growth conditions. This indicates that the growth condition of the Hg(1-x)Cd(x)Te is not purely diffusion controlled. Because of the higher thermal conductivity in the melt than that in the crystal in the growth system, there is a thermal leakage through the fused silica crucible wall near the melt-solid interface. This gives a thermal gradient in the radial direction. Hart, and Thorpe, Hutt and Soulsby have shown that under such conditon a fluid will become convectively unstable as a result of different diffusitivities of temperature and solute. It is quite important to understand the effects of this thermosolute convection on the compositonal segregation in both axial and radial directions in the unidirectionally solidified crystals under various gravitational directions. To reach this goal, we start with a simplified problem to study the effects of thermal-solutal convection on the temperature and solutal fields under various gravitional orientations. We begin by reviewing model governing equations.

Wang, Jai-Ching↗

DOE OSTI2020

Localization in the SCAN meta-generalized gradient approximation functional leading to broken symmetry ground states for graphene and benzene

Density functional theory calculations play a central role in understanding chemical and solid-state systems. Progress depends on density functionals that accurately reproduce both energies, for thermochemistry, and properly describe ground states and other properties that are of interest. The Cr dimer, benzene and graphene are especially important benchmark systems for quantum chemistry and condensed matter physics. The Strongly Constrained and Appropriately Normed (SCAN) functional, which is an advanced meta-generalized gradient approximation functional that significantly improves molecular energies is shown to perform poorly for the Cr dimer. This is connected with its poor performance for itinerant solid-state magnets and is a consequence of over localization of electrons, thus illustrating an analogy between the Cr dimer and itinerant magnets. The Cr dimer is a notoriously difficult system for density functionals. However, we additionally find that SCAN predicts an incorrect symmetry broken ground state for 2D graphene and for the benzene molecule, which is surprising considering that ground states of these are known to be well described even by the simplest local density approximation. We show that SCAN overly favors localized spin polarized states, which is a serious deficiency of this approach. Thus, the challenge of finding density functionals that accurately treat both localized and delocalized electronic systems remains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1999

Double Diffusive Convection in Materials Processing

A great number of crystals grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity (g-jitter). As a specific example, past HgCdTe crystal growth space experiments by Lehoczky and co-workers indicate radial compositional asymmetry in the grown crystals. In the case of HgCdTe the rejected component into the melt upon solidification is HgTe which is denser than the melt. The space grown crystals indicate the presence of three dimensional flow with the heavier HgTe-rich material clearly aligned with the residual gravity (0.55-1.55 micro g) vector. This flow stems from double-diffusive convection, namely, thermal and solutal buoyancy driven flow in the melt. The study of double-diffusive convection is multi-faceted and rather vast. In our investigation, we seek to focus on one specific aspect of this discipline that is of direct relevance to materials processing especially crystal growth, namely, the side ways heating regime. This problem has been widely studied, both experimentally and numerically, in the context of solar ponds wherein the system is characterized by a linear salt (solutal) gradient with an imposed lateral temperature gradient. The induced flow instabilities arise from the wide disparity between the fluid thermal diffusivity and the solute diffusivity. The extension of the analysis to practical crystal growth applications has however not been rigorously made and understood. One subtle but important difference in crystal growth systems is the fact that die system solute gradient is non-linear (typically exponential). Besides, the crystal growth problem has the added complexities of solidification, both lateral and longitudinal thermal gradients and segregation phenomena in systems where binary and ternary compounds are being grown. This paper treats the side ways heating problem alone in a model fluid system. Results from detailed numerical calculations, mainly two dimensional are provided. The interactions between a non-linear solute gradient and an imposed transverse thermal gradient are investigated. The buoyancy effects are treated in the traditional Boussinesq approximation and also in a more complete density formulation to address recent concerns of the first approach especially in simulations of the system response in a reduced gravity environment. Detailed flow, temperature and solute field plots along with heat and mass transfer results are presented in the paper. Implications to practical crystal growth systems as discerned from the modeling results are also explored and reported.

Ramachandra, Narayanan↗

NASA NTRS1973

On empirical models of the upper atmosphere in the polar regions.

The expression for the exospheric temperature in Jacchia's static diffusion models of the upper atmosphere has a discontinuous gradient at the poles. Therefore it cannot describe the true state of the upper atmosphere in the polar regions. Furthermore, it cannot be used to calculate quantities that depend on the derivative of the exospheric temperature, or the density, like pressure gradients, horizontal forces or horizontal heat fluxes. A modified expression for the exospheric temperature is suggested. This modification yields variables of state of the upper atmosphere that deviate little from Jacchia's values, but it has continuous gradient at the poles and is therefore more suitable for treating dynamical problems like the global wind pattern.

Blum, P. W.↗

DOE OSTI2015

Pressure dependence of electronic properties of BaI{sub 2}

We present Density Functional Theoretical (DFT) calculations of the electronic properties of scintillator material BaI{sub 2} under pressure which were carried out using Perdew-Burke-Ernzerhof genralized gradient approximation. We found that BaI{sub 2} is a direct band gap material with band gap calculated as 3.35 eV. The pressure effect on BaI{sub 2} indicates a linear monotonously decreasing band gap and increasing valence band width with pressure. We have observed the shifting of band extrema from the Γ point with pressure. The pressure coefficient of band gap is found to be −0.047 eV/GPa. The interatomic ionicity factor of BaI{sub 2} is found to be 0.51. Trends in bonding and ionicity under pressure are also discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data

In this work, we present a machine learning approach to calculating electronic specific heat capacities for a variety of benchmark molecular systems. Our models are based on data from density matrix quantum Monte Carlo, which is a stochastic method that can calculate the electronic energy at finite temperature. As these energies typically have noise, numerical derivatives of the energy can be challenging to find reliably. In order to circumvent this problem, we use Gaussian process regression to model the energy and use analytical derivatives to produce the specific heat capacity. From there, we also calculate the entropy by numerical integration. We compare our results to cubic splines and finite differences in a variety of molecules in which Hamiltonians can be diagonalized exactly with full configuration interaction. We finally apply this method to look at larger molecules where exact diagonalization is not possible and make comparisons with more approximate ways to calculate the specific heat capacity and entropy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗