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 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2013

Temperature-driven motion of liquid crystals confined in a microvolume

The theoretical description of the reorientational dynamics in microsized liquid crystal (LC) cell, where the nematic sample is confined by two horizontal and two lateral surfaces, under the influence of a temperature gradient ∇T, caused by a laser beam focused on the bounding surface with and without the orientational defect, whereas the rest of the bounding surfaces of the LC cell are kept at constant temperature, has been presented. Our calculations, based on the appropriate nonlinear extension of the classical Ericksen-Leslie theory, show that due to interaction between ∇T and the gradient of the director field ∇n-circumflex in the LC sample, a thermally excited vortical fluid flow is maintained in the vicinity of the orientational defect, with the motion in the positive sense (clockwise) around the middle part of that defect. In the case of the same LC cell, but without the orientational defect on the lower hotter boundary, the heating regime can also produce the vortical flow in the vicinity of the lower boundary, but with the motion in the negative sense (anti-clockwise) around the middle part of that boundary. At that, the second vortex is characterized by a much slower speed than the vortical flow in the first case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Effect of mechanical strain on electronic properties of bulk MoS{sub 2}

Ab-initio density functional theory based calculations of electronic properties of bulk and monolayer Molybdenum di-Sulfide (MoS{sub 2}) have been performed using all electron Full Potential Linearised Augmentad Plane Wave (FPLAPW) method using Elk code. We have used Generalised Gradient Approximation (GGA) for exchange and correlation functionals and performed calculaitons of Lattice parameters, Density Of States (DOS) and Band Structure (BS). Band structure calculations revealed that bulk MoS{sub 2} has indirect band gap of 0.97 eV and mono-layer MoS{sub 2} has direct band gap which has increased to 1.71 eV. These are in better agreement with experimental values as compared with the other calculations using pseudo-potential code. The effect of mechanical strain on the electronic properties of bulk MoS{sub 2} has also been studied. For the different values of compressive strain (varying from 2% to 8% in steps of 2%) along the c-axis, the corresponding DOS and BS are obtained. We observed that the band gap decreases by about 15% for every 2% increase in strain along the c-axis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1983

Stress studies in EFG

Stress distributions were calculated for a creep law to predict a rate of plastic deformation. The expected reduction in stresses is obtained. Improved schemes for calculating growth system temperature distributions were evaluated. Temperature field modeling examined the possibility of using horizontal temperature gradients to influence stress distribution in ribbon. The defect structure of 10 cm wide ribbon grown in the cartridge system was examined. A new feature is identified from an examination of cross sectional micrographs. It consists of high density dislocation bands extending through the ribbon thickness. A four point bending apparatus was constructed for high temperature study of the creep response of silicon, to be used to generate defects for comparison with as grown defects in ribbon. The feasibility of laser interferometric techniques for sheet residual stress distribution measurement is examined. The mathematical formalism for calculating residual stress from changes in surface topology caused by an applied stress in a rectangular specimen was developed, and the system for laser interferometric measurement to obtain surface topology data was tested on CZ silicon.

Source record↗

NASA NTRS1990

Ecosystem approach to a global nitrous oxide budget

An approach using relationships among soil fertility, nitrogen cycling, and nitrous oxide production is presented in order to estimate nitrous oxide flux from humid tropical forests. The effects of human disturbance are also investigated. It is noted that recent estimates suggest that 6-10 million ha. of tropical forest are cleared permanently each year. The results of measurements of nitrous oxide flux from cleared and burned sites in Costa Rica and Brazil, in conjunction with satellite-based classifications of ecosystem types, are used to calculate average hourly fluxes for a given region. In the NASA intensive site near Manaus, Brazil, it was found that although pasture areas cover only 11 percent of the region they account for 40 percent of the nitrous oxide flux. This analysis of the fluxes of trace gases through consideration of the gradients of factors that control both fluxes and ecosystem properties and processes will prove useful for extrapolating fluxes and for calculating budgets of nitrous oxide, methane, nonmethane hydrocarbons, and carbon dioxide.

Matson, P. A.↗

DOE OSTI2013

Combined short scale roughness and surface dielectric function gradient effects on the determination of tip-sample force in atomic force microscopy

The contribution of tip roughness to the van der Waals force between an atomic force microscopy probe tip and the sample is calculated using the multilayer effective medium model, which allows us to consider the relevant case of roughness characterized by correlation length and amplitude in the nanometer scale. The effect of the surface dielectric function gradient is incorporated in the tip-sample force model. It is concluded that for rms roughness in the few nanometers range the effect of short scale tip roughness is quite significant.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2013

Thermal transpiration through single walled carbon nanotubes and graphene channels

Thermal transpiration through carbon nanotubes (CNTs) and graphene channels is studied using molecular dynamics (MD) simulations. The system consists of two reservoirs connected by a CNT. It is observed that a flow is developed inside the CNT from the low temperature reservoir to the high temperature reservoir when the two reservoirs are maintained at different temperatures. The influence of channel size and temperature gradient on the mean velocity is analysed by varying the CNT diameter and the temperature of one of the reservoirs. Larger flow rate is observed in the smaller diameter CNTs showing an increase in the mean velocity with increase in the temperature gradient. For the flow developed inside the CNTs, slip boundaries occur and the slip length is calculated using the velocity profile. We examine the effect of fluid-wall interaction strength (ε{sub fw}), diffusivity (D), and viscosity of the fluid (μ) on the temperature induced fluid transport through the CNTs. Similar investigations are also carried out by replacing the CNT with a graphene channel. Results show that the mean velocity of the fluid atoms in the graphene channel is lower than that through the CNTs. This can be attributed to the higher degree of confinement observed in the CNTs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1994

Near-surface thermal gradients and their effects on mid-infrared emission spectra of planetary surfaces

We model the heat transfer by radiation and conduction in the top few millimeters of a planetary surface to determine the magnitude of near-surface (approximately 100 micrometers) thermal gradients and their effects on mid-infrared emission spectra for a number of planetary environments. The model is one-dimensional and uses a finite difference scheme for approximately 10 micrometers layers. Calculations are peformed for samples heated at the base and from above by sunlight. Our results indicate that near-surface radiative cooling creates significant thermal gradients in the top few hundred microns of surfaces in which radiation is an importamnt heat transfer mechanism. The effect is maximized in evacuated, underdense particulate media with sufficiently high temperatures. Near-surface thermal gradients will be significant in fine-grained particulate surfaces on the Moon (40-60 K/100 micrometers) and Mercury (approximately 80 K/100 micrometers), increasing spectral contrast and creating emission maxima in the transparent regions of the spectra. They will be of lesser importance on the surface of Mars, with a maximum value of around 5 k/100 micrometers in areas of low thermal inertia, and will be negligible on planets with more substantial atmospheres (less than 1 K/100 micrometers). We conclude that the effects that thermal gradients have on mid-IR emission spectra are predictable and do not negate the utility of emission spectroscopy for remote determination of planetary surface composition.

Henderson, B. G.↗

NASA NTRS1971

Proposed techniques for launching instrumented balloons into tornadoes

A method is proposed to introduce instrumented balloons into tornadoes by means of the radial pressure gradient, which supplies a buoyancy force driving to the center. Presented are analytical expressions, verified by computer calculations, which show the possibility of introducing instrumented balloons into tornadoes at or below the cloud base. The times required to reach the center are small enough that a large fraction of tornadoes are suitable for the technique. An experimental procedure is outlined in which a research airplane puts an instrumented, self-inflating balloon on the track ahead of the tornado. The uninflated balloon waits until the tornado closes to, typically, 750 meters; then it quickly inflates and spirals up and into the core, taking roughly 3 minutes. Since the drive to the center is automatically produced by the radial pressure gradient, a proper launch radius is the only guidance requirement.

Grant, F. C.↗

DOE OSTI2016

Multilevel summation with B-spline interpolation for pairwise interactions in molecular dynamics simulations

The multilevel summation method for calculating electrostatic interactions in molecular dynamics simulations constructs an approximation to a pairwise interaction kernel and its gradient, which can be evaluated at a cost that scales linearly with the number of atoms. The method smoothly splits the kernel into a sum of partial kernels of increasing range and decreasing variability with the longer-range parts interpolated from grids of increasing coarseness. Multilevel summation is especially appropriate in the context of dynamics and minimization, because it can produce continuous gradients. This article explores the use of B-splines to increase the accuracy of the multilevel summation method (for nonperiodic boundaries) without incurring additional computation other than a preprocessing step (whose cost also scales linearly). To obtain accurate results efficiently involves technical difficulties, which are overcome by a novel preprocessing algorithm. Numerical experiments demonstrate that the resulting method offers substantial improvements in accuracy and that its performance is competitive with an implementation of the fast multipole method in general and markedly better for Hamiltonian formulations of molecular dynamics. The improvement is great enough to establish multilevel summation as a serious contender for calculating pairwise interactions in molecular dynamics simulations. In particular, the method appears to be uniquely capable for molecular dynamics in two situations, nonperiodic boundary conditions and massively parallel computation, where the fast Fourier transform employed in the particle–mesh Ewald method falls short.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2003

Measuring Diffusion of Liquids by Common-Path Interferometry

A method of observing the interdiffusion of a pair of miscible liquids is based on the use of a common-path interferometer (CPI) to measure the spatially varying gradient of the index refraction in the interfacial region in which the interdiffusion takes place. Assuming that the indices of refraction of the two liquids are different and that the gradient of the index of refraction of the liquid is proportional to the gradient in the relative concentrations of either liquid, the diffusivity of the pair of liquids can be calculated from the temporal variation of the spatial variation of the index of refraction. This method yields robust measurements and does not require precise knowledge of the indices of refraction of the pure liquids. Moreover, the CPI instrumentation is compact and is optomechanically robust by virtue of its common- path design. The two liquids are placed in a transparent rectangular parallelepiped test cell. Initially, the interface between the liquids is a horizontal plane, above which lies pure liquid 2 (the less-dense liquid) and below which lies pure liquid 1 (the denser liquid). The subsequent interdiffusion of the liquids gives rise to a gradient of concentration and a corresponding gradient of the index of refraction in a mixing layer. For the purpose of observing the interdiffusion, the test cell is placed in the test section of the CPI, in which a collimated, polarized beam of light from a low-power laser is projected horizontally through a region that contains the mixing layer.

Rashidnia, Nasser↗

DOE OSTI2018

Prediction of continuous porosity gradients in ceramics using ZnO as a model material

Materials with gradient microstructures have a wide range of applications such as cutting tools, armor, and electronic devices. However, it is difficult to predict and control the gradient microstructure during processing. In the present work, a continuous porosity gradient was successfully achieved in ZnO material via spark plasma sintering with a large induced thermal gradient. The porosity is overestimated if isothermal prediction is applied. The current work proposed a more accurate prediction of the porosity by considering the stress-shielding effect caused by the thermal gradient. The shielding effect results from different stress states in the sample due to differential sintering: the hotter side of the specimen experiences a higher strain rate and more shrinkage while the colder side experiences a lower strain rate and less shrinkage simultaneously. Therefore, the axial strains are varied throughout the sample thickness. Using the constituent equations in advanced sintering analysis, the shield stress was calculated to be approximately 13 MPa for the viscoelastic assumption of sintering. To improve the accuracy of predicting porosity gradient, it is necessary to add a load to overcome the shield stress when the materials are sintered with a thermal gradient.

36 MATERIALS SCIENCE↗

DOE OSTI2015

The complex mixed Wentzel–Kramers–Brillouin-full-wave approach and its application to the two dimensional mode structure analysis of ion temperature gradient/collisionless trapped electron mode drift waves

The complex mixed Wentzel–Kramers–Brillouin (WKB)-full-wave approach is applied to the 2D mode structure analysis of ion temperature gradient/collisionless trapped electron mode drift waves in tokamak plasmas. The parallel mode structure is calculated with the full-wave approach, while the radial envelope is calculated with the complex WKB method. The tilting of the global mode structure along radius is demonstrated analytically. The effects of the phase and amplitude variation of the radial envelope on the parallel mode structure are included in terms of a complex radial wave vector in the parallel mode equation. It is shown that the radial equilibrium non-uniformity leads to the asymmetry of the parallel mode structure not only in configuration space but also in spectrum space. The mixed approach provides a practical way to analyze the asymmetric component of the global mode structure due to radial equilibrium non-uniformity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Push it to the limit: comparing periodic and local approaches to density functional theory for intermolecular interactions

With the aim of systematically comparing two popular approaches to density functional theory–all-electron calculations with local basis sets, and periodic calculations employing plane wave basis sets and norm-conserving pseudopotentials–we have computed complete-basis binding energies across the S22 set of intermolecular interactions, a dataset consisting of noncovalent interactions of small- and medium-sized molecules containing first- and second-row atoms, using the Troullier-Martins norm-conserving pseudopotentials with SPW92, a local spin-density approximation; and PBE, a generalised gradient approximation. Here, we have found that it is challenging to reach the basis set limit with these periodic calculations; for the methods and systems examined, a minimum vacuum distance of 30 Å between a system and its nearest images is necessary–unless some form of dipole correction is employed–as is a kinetic energy cutoff of at least 80 Ry. The trends in convergence with respect to vacuum size and kinetic energy cutoff are largely independent of the level of density functional approximation employed. A sense of the impact of each hyperparameter on basis set error provides a foundation for ensuring quality calculations in future studies and allows us to quantify the basis set errors incurred in existing studies on similar systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS2012

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗

NASA NTRS1988

Evaluation of ERBE scanner pointing accuracy based upon a coastline detection algorithm

Measurements from the Earth Radiation Budget Experiment (ERBE) scanning radiometers must be accurately located at the top of the Earth's atmosphere for proper interpretation of derived longwave and shortwave components of the Earth radiation field. The accuracy of the location calculations is affected by the orbit characteristics and attitude control systems of the three spacecraft, ERBS, NOAA-9, and NOAA-10. A technique was developed which makes use of the longwave scanner measurements to detect the thermal gradients at the boundaries of ocean and land masses. By analyzing many such coastline boundaries, estimates of the overall accuracy of the scanner measurement location calculations have been derived for the three satellites. The resulting measurement location errors are found to be smaller than the resolution of the radiometers, and this accuracy is well within the required location knowledge for useful science analysis of the ERBE scanner data.

Weaver, William L.↗

NASA NTRS1974

On the distribution of He+ in the plasmasphere from observations of resonantly scattered He II 304-A radiation

The intensity of the He II 304-A radiation in the night sky was measured from a sounding rocket launched from White Sands Missile Range on June 9, 1972. The data obtained on this flight are compared with predictions of theoretical models of the He(+) altitude distribution within the plasmasphere. The observations are not consistent with an isothermal equilibrium model of the helium ion in the plasmasphere, confirming earlier results. The data are best fit by a diffussive equilibrium model of exospheric He(+) having a plasma temperature gradient of 0.5 K per km and a number density at 1600 km of 200,000 He(+) ions per cu cm. In this model the ion densities are calculated without allowance for thermal diffusion. A similar model that allows for thermal diffusion and that has a He(+) number density at 1600 km of 1500 per cu cm and a temperature gradient of 1.2 K per km yields a less satisfactory fit to the observations but is more consistent with in situ satellite observations.

Paresce, F.↗

NASA NTRS1993

Numerical modeling of crystal growth on a centrifuge for unstable natural convection configurations

The fluid mechanics associated with crystal growth processes on centrifuges is modeled using 2D and 3D models. Two-dimensional calculations show that flow bifurcations exist in such crystal growth configurations where the ampoule is oriented in the same direction as the resultant gravity vector and a temperature gradient is imposed on the melt. A scaling analysis is formulated to predict the flow transition point from the natural convection dominated regime to the Coriolis force dominated regime. Results of 3D calculations are presented for two thermal configurations of the crystal growth cell: top heated and bottom heated with respect to the centrifugal acceleration. In the top heated configuration, a substantial reduction in the convection intensity within the melt can be attained by centrifuge operations, and close to steady diffusion-limited thermal conditions can be achieved over a narrow range of the imposed microgravity level. In the bottom heated configuration the Coriolis force has a stabilizing effect on fluid motion by delaying the onset of unsteady convection.

Ramachandran, N.↗