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 559 records · Page 31

Global, Collisional Model of High-Energy Photoelectrons

A previously-developed colissional, interhemispheric flux tube model for photoelectrons (PE) has been extended to three dimensions by including transport due to vector E x vector B and magnetic gradient-curvature drifts. Using this model, initial calculations of the high-energy (greater then 50 eV) PE distribution as a function of time, energy, pitch angle, and spatial location in the equatorial plane, are reported for conditions of low geomagnetic activity. To explore both the dynamic and steady behaviors of the model, the simulation starts with the abrupt onset of photoelectron excitation, and is followed to steady state conditions. The results illustrate several features of the interaction of photoelectrons with typical magnetospheric plasmas and fields, including collisional diffusion of photoelectrons in pitch angle with flux tube filling, diurnal intensity and pitch angle asymmetries introduced by directional sunlight, and energization of the photoelectron distribution in the evening sector. Cross-field drift is shown to have a long time scale, taking 12 to 24 hours to reach a steady state distribution. Future applications of the model are briefly outlined.

Khazanov, G. V.↗

MAP stability, design, and analysis

The Microwave Anisotropy Probe (MAP) is a follow-on to the Differential Microwave Radiometer (DMR) instrument on the Cosmic Background Explorer (COBE) spacecraft. The design and analysis of the MAP attitude control system (ACS) have been refined since work previously reported. The full spacecraft and instrument flexible model was developed in NASTRAN, and the resulting flexible modes were plotted and reduced with the Modal Significance Analysis Package (MSAP). The reduced-order model was used to perform the linear stability analysis for each control mode, the results of which are presented in this paper. Although MAP is going to a relatively disturbance-free Lissajous orbit around the Earth-Sun L(2) Lagrange point, a detailed disturbance-torque analysis is required because there are only a small number of opportunities for momentum unloading each year. Environmental torques, including solar pressure at L(2), aerodynamic and gravity gradient during phasing-loop orbits, were calculated and simulated. Thruster plume impingement torques that could affect the performance of the thruster modes were estimated and simulated, and a simple model of fuel slosh was derived to model its effect on the motion of the spacecraft. In addition, a thruster mode linear impulse controller was developed to meet the accuracy requirements of the phasing loop burns. A dynamic attitude error limiter was added to improve the performance of the ACS during large attitude slews. The result of this analysis is a stable ACS subsystem that meets all of the mission's requirements.

Ericsson-Jackson, A. J.↗

MAP Stability, Design and Analysis

The Microwave Anisotropy Probe (MAP) is a follow-on to the Differential Microwave Radiometer (DMR) instrument on the Cosmic Background Explorer (COBE) spacecraft. The design and analysis of the MAP attitude control system (ACS) have been refined since work previously reported. The full spacecraft and instrument flexible model was developed in NASTRAN, and the resulting flexible modes were plotted and reduced with the Modal Significance Analysis Package (MSAP). The reduced-order model was used to perform the linear stability analysis for each control mode, the results of which are presented in this paper. Although MAP is going to a relatively disturbance-free Lissajous orbit around the Earth-Sun L2 Lagrange point, a detailed disturbance-torque analysis is required because there are only a small number of opportunities for momentum unloading each year. Environmental torques, including solar pressure at L2, and aerodynamic and gravity gradient during phasing-loop orbits, were calculated and simulated. A simple model of fuel slosh was derived to model its effect on the motion of the spacecraft. In addition, a thruster mode linear impulse controller was developed to meet the accuracy requirements of the phasing loop burns. A dynamic attitude error limiter was added to improve the performance of the ACS during large attitude slews. The result of this analysis is a stable ACS subsystem that meets all of the mission's requirements.

Ericsson -Jackson, A.J.↗

Closed Loop Software Control of the MIDEX Power System

The Microwave Anisotropy Probe (MAP) is a follow-on to the Differential Microwave Radiometer (DMR) instrument on the Cosmic Background Explorer (COBE) spacecraft. The design and analysis of the MAP attitude control system (ACS) have been refined since work previously reported. The full spacecraft and instrument flexible model was developed in NASTRAN, and the resulting flexible modes were plotted and reduced with the Modal Significance Analysis Package (MSAP). The reduced-order model was used to perform the linear stability analysis for each control mode, the results of which are presented in this paper. Although MAP is going to a relatively disturbance-free Lissajous orbit around the Earth-Sun L2 Lagrange point, a detailed disturbance-torque analysis is required because there are only a small number of opportunities for momentum unloading each year. Environmental torques, including solar pressure at L2, and aerodynamic and gravity gradient during phasing-loop orbits, were calculated and simulated. A simple model of fuel slosh was derived to model its effect on the motion of the spacecraft. In addition, a thruster mode linear impulse controller was developed to meet the accuracy requirements of the phasing loop burns. A dynamic attitude error limiter was added to improve the performance of the ACS during large attitude slews. The result of this analysis is a stable ACS subsystem that meets all of the mission's requirements.

Castell, Karen↗

Computation of Turbulent Heat Transfer on the Walls of a 180 Degree Turn Channel With a Low Reynolds Number Reynolds Stress Model

The Low Reynolds number version of the Stress-omega model and the two equation k-omega model of Wilcox were used for the calculation of turbulent heat transfer in a 180 degree turn simulating an internal coolant passage. The Stress-omega model was chosen for its robustness. The turbulent thermal fluxes were calculated by modifying and using the Generalized Gradient Diffusion Hypothesis. The results showed that using this Reynolds Stress model allowed better prediction of heat transfer compared to the k-omega two equation model. This improvement however required a finer grid and commensurately more CPU time.

Ameri, A. A.↗

Tracer Lamination in the Stratosphere: A Global Climatology

Vertical soundings of stratospheric ozone often exhibit laminated tracer structures characterized by strong vertical tracer gradients. The change in time of these gradients is used to define a tracer lamination rate. It is shown that this quantity can be calculated by the cross product of the horizontal temperature and horizontal tracer gradients. A climatology based on UARS satellite-borne ozone data and on ozone-like pseudotracer data is presented. Three stratospheric regions with high lamination rates were found: the part of the stratospheric overworld which is influenced by the polar vortex, the part of the lowermost stratosphere which is influenced by the tropopause and a third region in the subtropical lower stratosphere mainly characterized with strong vertical shear. High lamination rates in the stratospheric overworld were absent during summer, whereas in the lowermost stratosphere high lamination rates were found year-round. This is consistent with the occurrence and seasonal variation of the horizontal tracer gradient and vertical shear necessary for tilting the tracer surfaces. During winter, high lamination rates associated with the stratospheric polar vortex are present down to approximately 100 hPa. Several features of the derived climatology are roughly consistent with earlier balloon-borne studies. The patterns in the southern and northern hemisphere are comparable, but details differ as anticipated from a less disturbed and more symmetric southern polar vortex.

Appenzeller, Christof↗

A higher-order Lagrangian discontinuous Galerkin hydrodynamic method for solid dynamics

We present a new multidimensional high-order Lagrangian discontinuous Galerkin (DG) hydrodynamic method that supports hypoelastic and hyperelastic strength models for simulating solid dynamics with higher-order elements. We also present new one-dimensional test problems that have an analytic solution corresponding to a hyperelastic–plastic wave. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated high-order Taylor series polynomials. The stress fields are represented using nodal quantities. The constitutive models used to calculate the deviatoric stress are either a hypoelastic–plastic, infinitesimal strain hyperelastic–plastic, or finite strain hyperelastic–plastic model. These constitutive models require new methods for calculating high-order polynomials for the velocity gradient and deformation gradient in an element. The plasticity associated with the strength model is determined using a radial return method to with a J 2 yield criterion and perfect plasticity. The temporal evolution of the governing equations is achieved with the total variation diminishing Runge–Kutta (TVD RK) time integration method. A diverse suite of 1D and 2D test problems are calculated. Here, the new 1D piston test problems, which have analytic solutions for each elastic–plastic model, are presented and calculated to demonstrate the stability and formal accuracy of the various models with the new Lagrangian DG method. 2D test problems are calculated to demonstrate the stability and robustness of the new Lagrangian DG method on multidimensional problems with high-order elements, which have faces that can bend.

42 ENGINEERING↗

Effects of La Doping and Zn or O Vacancy on the Magnetic Property of ZnO

Research on the magnetic origin and mechanism of La-doped ZnO system has been controversial. To solve this problem, we studied the effects of La doping and Zn or O vacancy on the magnetic property of ZnO using geometry optimization and energy calculation on the basis of first-principle generalized gradient approximation method of density functional theory. The electronic structures and magnetic properties of undoped ZnO, La-doped ZnO system, and systems with coexisting La doping and Zn or O vacancies were calculated. Results showed that the systems of La-doped ZnO and the coexistence of La doping and O vacancy are nonmagnetic. In addition, the system of La replacing Zn and one Zn vacancy exhibits long-range orderly ferromagnetism, and the Curie temperature of the doping system can achieve room temperature. The magnetism source of the systems with La doping and Zn vacancy coexisting in ZnO demonstrates strong hybrid coupling electron exchange effects existing among Zn-4s, Zn-3p, O-2p, and La-5s orbits that are nearest to the Zn vacancy. The next nearest distance between doping and Zn vacancy leads to the lowest formation energy and highest stability under the same doping condition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Geometries, stabilities and electronic properties of photo sensitized (CdS) 2n (n=5–27) nanoclusters

As advanced materials for energy storage or conversion applications, metal sulfides with hollow structures are as widely used as supercapacitors and sensors. The geometries and electronic properties of (CdS) 2n (n = 5–27) nanocages are calculated here by density functional theory with the gradient corrected Perdew-Burke-Ernzerh of exchange-correlation functional and including the relativistic effect. The obtained results showed that the inner radius and charge-transfers, surface stabilities along with ionic semiconductor and metallic properties of the hollowed (CdS) 2n are increased as cage size being increased. The distinct metallic features in large-size (CdS) 2n are revealed. Particularly, the calculated energy gaps for (CdS) 2n are generally decreased as the size of (CdS) 2n being increased, together with absorption spectra moving from short to long wavelength. The photon to current conversion efficiency related energy conversion is discussed. Interestingly, the encapsulated double-layer CdS nanocages have bigger capacity for energy conversion than the monolayer nanocages.

36 MATERIALS SCIENCE↗

Linear simulation of kinetic electromagnetic instabilities in a tokamak plasma with weak magnetic shear

Gyrokinetic simulation and eigenvalue calculation of electromagnetic instabilities are carried out for an experimentally observed low-n mode in weak magnetic shear discharge. With different magnetic shear values, the Gyrokinetic Toroidal Code simulation of the ion temperature gradient mode is consistent with the eigenvalue code (HD7) calculation. Due to the sensitivity of the kinetic ballooning mode (KBM) to global equilibrium, the simulation of the KBM deviates from the eigenvalue results, for the ballooning representation used in HD7 satisfies its spatial scale separations. Under a flat safety factor profile, the KBM is more unstable and its mode structure tends to move with the peak of the ion temperature drive. Further simulation of the KBM in an HL-2A-like equilibrium shows that the β excitation threshold of the mode is lower than 0.2% and the dominant toroidal mode number is n = 4, which is consistent with the measured experimental spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Polarization of massive fermions in a vortical fluid

For this study, fermions become polarized in a vortical fluid due to spin-vorticity coupling. Such a polarization can be calculated from the Wigner function in a quantum kinetic approach. By extending previous results for chiral fermions, we derive the Wigner function for massive fermions up to next-to-leading order in spatial gradient expansion. The polarization density of fermions can be calculated from the axial vector component of the Wigner function and is found to be proportional to the local vorticity ω. The polarizations per particle for fermions and antifermions decrease with the chemical potential and increase with energy (mass). Both quantities approach the asymptotic value ℏω/4 in the large energy (mass) limit. The polarization per particle for fermions is always smaller than that for antifermions, whose ratio of fermions to antifermions also decreases with the chemical potential. The polarization per particle on the Cooper-Frye freeze-out hypersurface can also be formulated and is consistent with the previous result of Becattini et al.

43 PARTICLE ACCELERATORS↗

Application of a planetary wave breaking parameterization to stratospheric circulation statistics

The planetary wave parameterization scheme developed recently by Garcia is applied to statospheric circulation statistics derived from 12 years of National Meteorological Center operational stratospheric analyses. From the data a planetary wave breaking criterion (based on the ratio of the eddy to zonal mean meridional potential vorticity (PV) gradients), a wave damping rate, and a meridional diffusion coefficient are calculated. The equatorward flank of the polar night jet during winter is identified as a wave breaking region from the observed PV gradients; the region moves poleward with season, covering all high latitudes in spring. Derived damping rates maximize in the subtropical upper stratosphere (the 'surf zone'), with damping time scales of 3-4 days. Maximum diffusion coefficients follow the spatial patterns of the wave breaking criterion, with magnitudes comparable to prior published estimates. Overall, the observed results agree well with the parameterized calculations of Garcia.

Randel, William J.↗

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development of an efficient procedure for calculating the aerodynamic effects of planform variation

Numerical procedures to compute gradients in aerodynamic loading due to planform shape changes using panel method codes were studied. Two procedures were investigated: one computed the aerodynamic perturbation directly; the other computed the aerodynamic loading on the perturbed planform and on the base planform and then differenced these values to obtain the perturbation in loading. It is indicated that computing the perturbed values directly can not be done satisfactorily without proper aerodynamic representation of the pressure singularity at the leading edge of a thin wing. For the alternative procedure, a technique was developed which saves most of the time-consuming computations from a panel method calculation for the base planform. Using this procedure the perturbed loading can be calculated in about one-tenth the time of that for the base solution.

Mercer, J. E.↗

Planning fuel-conservative descents in an airline environmental using a small programmable calculator: Algorithm development and flight test results

A simple, airborne, flight-management descent algorithm was developed and programmed into a small programmable calculator. The algorithm may be operated in either a time mode or speed mode. The time mode was designed to aid the pilot in planning and executing a fuel-conservative descent to arrive at a metering fix at a time designated by the air traffic control system. The speed model was designed for planning fuel-conservative descents when time is not a consideration. The descent path for both modes was calculated for a constant with considerations given for the descent Mach/airspeed schedule, gross weight, wind, wind gradient, and nonstandard temperature effects. Flight tests, using the algorithm on the programmable calculator, showed that the open-loop guidance could be useful to airline flight crews for planning and executing fuel-conservative descents.

Knox, C. E.↗

Cooling flows and the stability of radio jets

This paper considers what happens to radio jets which, propagating out from the nucleus of an elliptical galaxy, encounters a point at which the ambient gas pressure increases very rapidly. The pressure variation within the sonic radius of a cooling flow is calculated, and a semianalytical calculation of the stability of a jet which encounters a rapidly increasing pressure gradient is given. The results of numerical hydrodynamic simulations which test the results of the analytic calculations are shown. The case where the jet encounters the sonic region at an oblique angle is considered, and it is suggested that this may provide an explanation for jets which are bent on small scales. These results are applied to radio sources associated with central dominant cluster galaxies and with more normal elliptical galaxies.

Soker, Noam↗

Distributed gradient descent for solving optimal power flow in radial networks

Node controllers and power distribution networks in accordance with embodiments of the invention enable distributed power control. One embodiment includes a node controller comprising a memory containing: a plurality of node operating parameters; and a plurality of node operating parameters describing operating parameters for a set of at least one node selected from the group consisting of at least one downstream node and at least one upstream node; wherein the processor is configured by the node controller application to: receive and store in memory a plurality of coordinator parameters describing operating parameters of a node coordinator by the network interface; and calculate updated node operating parameters using an iterative gradient projection process to determine updated node parameters using node operating parameters that describe operating parameters of node and operating parameters of the set of at least one node, where each iteration is determined by the coordinator parameters.

Gan, Lingwen↗

Stabilization of electron-scale turbulence by electron density gradient in national spherical torus experiment

Theory and experiments have shown that electron temperature gradient (ETG) turbulence on the electron gyro-scale, k{sub ⊥}ρ{sub e} ≲ 1, can be responsible for anomalous electron thermal transport in NSTX. Electron scale (high-k) turbulence is diagnosed in NSTX with a high-k microwave scattering system [D. R. Smith et al., Rev. Sci. Instrum. 79, 123501 (2008)]. Here we report on stabilization effects of the electron density gradient on electron-scale density fluctuations in a set of neutral beam injection heated H-mode plasmas. We found that the absence of high-k density fluctuations from measurements is correlated with large equilibrium density gradient, which is shown to be consistent with linear stabilization of ETG modes due to the density gradient using the analytical ETG linear threshold in F. Jenko et al. [Phys. Plasmas 8, 4096 (2001)] and linear gyrokinetic simulations with GS2 [M. Kotschenreuther et al., Comput. Phys. Commun. 88, 128 (1995)]. We also found that the observed power of electron-scale turbulence (when it exists) is anti-correlated with the equilibrium density gradient, suggesting density gradient as a nonlinear stabilizing mechanism. Higher density gradients give rise to lower values of the plasma frame frequency, calculated based on the Doppler shift of the measured density fluctuations. Linear gyrokinetic simulations show that higher values of the electron density gradient reduce the value of the real frequency, in agreement with experimental observation. Nonlinear electron-scale gyrokinetic simulations show that high electron density gradient reduces electron heat flux and stiffness, and increases the ETG nonlinear threshold, consistent with experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗