Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density gradient function”

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 145 records · Page 8

An analysis of the Voyager 2 Ultraviolet Spectrometer occultation data at Uranus - Inferring heat sources and model atmospheres

Heat source information is derived here from the Voyager 2 Ultraviolet Spectrometer occultation data of Uranus. Analytic functions for the local heat dependence on altitude are used to obtain a temperature profile by solving the heat equation. The stellar entrance and exit occultation and a solar occultation are used to infer the thermal and density structure of the atmosphere. The least squares fit solution to the solar occultation data gives one source located at 1.8 x 10 exp -5 microbar with a strength of 0.056 +/- 0.01 erg/sq cm/s. Latitudinal temperature gradients are obtained.

Stevens, Michael H.↗

Detailed Studies on the Structure and Dynamics of Reacting Dusty Flows at Normal and Microgravity

Dusty reacting flows are of particular interest for a wide range of applications. Inert particles can alter the flammability and extinction limits of a combustible mixture. Reacting particles can release substantial amount of heat and can be used either for power generation or propulsion. Accumulation of combustible particles in air can result in explosions which, for example, can occur in grain elevators, during lumber milling and in mine galleries. Furthermore, inert particles are used as flow velocity markers in reacting flows, and their velocity is measured by non-intrusive laser diagnostic techniques. Despite their importance, dusty reacting flows have been less studied and understood compared to gas phase as well as sprays. The addition of solid particles in a flowing gas stream can lead to strong couplings between the two phases, which can be of dynamic, thermal, and chemical nature. The dynamic coupling between the two phases is caused by the inertia that causes the phases to move with different velocities. Furthermore, gravitational, thermophoretic, photophoretic, electrophoretic, diffusiophoretic, centrifugal, and magnetic forces can be exerted on the particles. In general, magnetic, electrophoretic, centrifugal, photophoretic, and diffusiophoretic can be neglected. On the other hand, thermophoretic forces, caused by steep temperature gradients, can be important. The gravitational forces are almost always present and can affect the dynamic response of large particles. Understanding and quantifying the chemical coupling between two phases is a challenging task. However, all reacting particles begin this process as inert particles, and they must be heated before they participate in the combustion process. Thus, one must first understand the interactions of inert particles in a combustion environment. The in-detail understanding of the dynamics and structure of dusty flows can be only advanced by considering simple flow geometries such as the opposed-jet, stagnation-type. In such configurations the imposed strain rate is well characterized, and the in-depth understanding of the details of the physico-chemical processes can be systematically obtained. A number of computational and experimental studies on spray and particle flows have been conducted in stagnation-type configurations. Numerically, the need for a hybrid Eulerian-Lagrangian approach has been identified by Continillo and Sirignano, and the use of such approach has allowed for the prediction of the phenomenon of droplet flow reversal. Gomez and Rosner have conducted a detailed study on the particle response in the opposed-jet configuration, and the particle thermophoretic diffusivities were determined experimentally. Sung, Law and co-workers have conducted numerical studies on the effect of strain rate and temperature gradients on the dynamics of inert particles, as a way of understanding potential errors in experimental LDV data that may arise from thermophoretic forces. This investigation is a combined experimental and numerical study on the details of reacting dusty flows. The specific tasks are: (1) Experimental determination of laminar flame speeds, and extinction strain rates of dusty flows at normal- and micro-gravity as functions of the particle type, particle initial diameter, particle initial number density, and gas phase chemical composition; (2) Detailed numerical simulation of the experiments. Results are compared with experiments and the adequacy of theoretical models is assessed; and (3) Provision of enhanced insight into the thermo-chemical coupling between the two phases.

Andac, M. Gurhan↗

Nano-Transistor Modeling: Two Dimensional Green's Function Method

Two quantum mechanical effects that impact the operation of nanoscale transistors are inversion layer energy quantization and ballistic transport. While the qualitative effects of these features are reasonably understood, a comprehensive study of device physics in two dimensions is lacking. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL (Drain Induced Barrier Lowering), and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density-gradient and quantum-corrected MEDICI).

Svizhenko, Alexei↗

Supersonic free turbulent mixing layers

A mean turbulent closure method with a length scale function is developed for predicting the reduced spreading rate in high Mach number, high Reynolds number free shear layers. The new feature of the method is the inclusion of effects of the pressure-dilation term which appears in the compressible form of the turbulence energy equation. The model indicates that the pressure-dilation correlation may act as a turbulence energy sink or source depending on whether the transverse gradients of mean Mach number and mean density have the same or opposite signs.

Oh, Y. H.↗

Density Irregularities Inside the Plasmasphere: Cluster Observations

The electron density profiles derived from the EFW and WHISPER instruments onboard the four Cluster spacecraft reveal small-scale density irregularities inside the plasmasphere and at its outer boundary, the plasmapause. We review statistics of the plasmapause position and thickness, as well as statistics of these density structures. We focus on a particular plasmasphere crossing on 11 April 2002, with several density irregularities, as well as two plasma tails, observed by Cluster on the legs of the inbound and outbound passes. We derive the density gradient vectors from simultaneous density measurements by the four spacecraft. We determine also the normal velocity of the surface of these irregularities, assuming they are planar boundaries, from the time delays between density structures in the four individual density profiles. These new observations yield novel insights about (1) the dimensions of plasma irregularities across and parallel to magnetic field lines, (2) the dynamics of these small-scale structures, (3) their bulk velocities, and (4) their position and distribution as a function of the magnetic local time (MLT) and the geomagnetic conditions (as determined by Kp).

Darrouzet, F.↗

Direct numerical simulations of stably-stratified sheared turbulence: Implications for oceanic mixing

Direct numerical simulations of the time evolution of homogeneous stably stratified turbulent shear flows have been performed for several Richardson numbers Ri and Reynolds numbers R(sub lambda) in earlier works. The results show excellent agreement with length scale models developed from laboratory experiments to characterize oceanic turbulence. When the Richardson number Ri is less than the stationary value Ri(sub s), the turbulence intensity grows at all scales, and the growth rate appears to be a function of Ri. The size of the vertical density inversions also increases. On the other hand, when Ri is greater than or equal to Ri(sub s) the largest turbulent eddies become vertically constrained by buoyancy when the Ellison (turbulence) scale L(sub E) and the Ozmidov (buoyancy) scale L(sub O) are equal. At this point, the mixing efficiency is maximal and corresponds to a flux Richardson number R(sub f) = 0.20. The vertical mass flux becomes counter-gradient when epsilon = 19(nu)N(exp 2) and vertical density overturns are suppressed in less than half a Brunt-Vaisala period. The results of the simulations were also recast in terms of the Hydrodynamic Phase Diagram introduced in fossil turbulence models. The so-called point of fossilization occurs when epsilon = 4DCN(exp 2); Gibson proposed 13DCN(exp 2). This value is in agreement with indirect laboratory observations and field observations. Finally, the validity of the steady-state models to estimate vertical eddy diffusivities in the oceanic thermocline is discussed.

Itsweire, E. C.↗

Investigations of turbulent scalar fields using probability density function approach

Scalar fields undergoing random advection have attracted much attention from researchers in both the theoretical and practical sectors. Research interest spans from the study of the small scale structures of turbulent scalar fields to the modeling and simulations of turbulent reacting flows. The probability density function (PDF) method is an effective tool in the study of turbulent scalar fields, especially for those which involve chemical reactions. It has been argued that a one-point, joint PDF approach is the one to choose from among many simulation and closure methods for turbulent combustion and chemically reacting flows based on its practical feasibility in the foreseeable future for multiple reactants. Instead of the multi-point PDF, the joint PDF of a scalar and its gradient which represents the roles of both scalar and scalar diffusion is introduced. A proper closure model for the molecular diffusion term in the PDF equation is investigated. Another direction in this research is to study the mapping closure method that has been recently proposed to deal with the PDF's in turbulent fields. This method seems to have captured the physics correctly when applied to diffusion problems. However, if the turbulent stretching is included, the amplitude mapping has to be supplemented by either adjusting the parameters representing turbulent stretching at each time step or by introducing the coordinate mapping. This technique is still under development and seems to be quite promising. The final objective of this project is to understand some fundamental properties of the turbulent scalar fields and to develop practical numerical schemes that are capable of handling turbulent reacting flows.

Gao, Feng↗

On the radial structure of planetary rings

Viscous shear stress of a ring of particles in orbit about a planet due to the radial gradient of orbital velocity is discussed. This stress tends to spread the ring with time. At low optical depth (t) and at high t, the shear stress is an increasing function of t. In the intermediate range stress may decrease with increasing t, leading to a diffusive instability which tends to break an initially uniform ring into ringlets of high and low optical depths. According to the shepherd satellite model of narrow ring confinement, the viscous shear stress is opposed by tidal torques from neighboring shepherds, so that radial spreading does not occur. By requiring that the gradient of the viscous shear stress (i.e., angular momentum deposited) be equal to the shepherd satellite torque density at all radii in the rings, equilibrium radial density profiles of a confined ring were constructed. The most noteworthy feature of these profiles is a narrow, high density ringlet surrounded by a low optical density halo of greater width.

Harris, A. W.↗

Semidiscrete Galerkin modelling of compressible viscous flow past a circular cone at incidence

A numerical study of the laminar and compressible boundary layer, about a circular cone in a supersonic free stream, is presented. It is thought that if accurate and efficient numerical schemes can be produced to solve the boundary layer equations, they can be joined to numerical codes that solve the inviscid outer flow. The combination of these numerical codes is competitive with the accurate, but computationally expensive, Navier-Stokes schemes. The primary goal is to develop a finite element method for the calculation of 3-D compressible laminar boundary layer about a yawed cone. The proposed method can, in principle, be extended to apply to the 3-D boundary layer of pointed bodies of arbitrary cross section. The 3-D boundary layer equations governing supersonic free stream flow about a cone are examined. The 3-D partial differential equations are reduced to 2-D integral equations by applying the Howarth, Mangler, Crocco transformations, a linear relation between viscosity, and a Blasius-type of similarity variable. This is equivalent to a Dorodnitsyn-type formulation. The reduced equations are independent of density and curvature effects, and resemble the weak form of the 2-D incompressible boundary layer equations in Cartesian coordinates. In addition the coordinate normal to the wall has been stretched, which reduces the gradients across the layer and provides high resolution near the surface. Utilizing the parabolic nature of the boundary layer equations, a finite element method is applied to the Dorodnitsyn formulation. The formulation is presented in a Petrov-Galerkin finite element form and discretized across the layer using linear interpolation functions. The finite element discretization yields a system of ordinary differential equations in the circumferential direction. The circumferential derivatives are solved by an implicit and noniterative finite difference marching scheme. Solutions are presented for a 15 deg half angle cone at angles of attack of 5 and 10 deg. The numerical solutions assume a laminar boundary layer with free stream Mach number of 7. Results include circumferential distribution of skin friction and surface heat transfer, and cross flow velocity distributions across the layer.

Meade, Andrew James, Jr.↗

Radial basis function network learns ceramic processing and predicts related strength and density

Radial basis function (RBF) neural networks were trained using the data from 273 Si3N4 modulus of rupture (MOR) bars which were tested at room temperature and 135 MOR bars which were tested at 1370 C. Milling time, sintering time, and sintering gas pressure were the processing parameters used as the input features. Flexural strength and density were the outputs by which the RBF networks were assessed. The 'nodes-at-data-points' method was used to set the hidden layer centers and output layer training used the gradient descent method. The RBF network predicted strength with an average error of less than 12 percent and density with an average error of less than 2 percent. Further, the RBF network demonstrated a potential for optimizing and accelerating the development and processing of ceramic materials.

Cios, Krzysztof J.↗

Large infrared and optical color gradients in the Cartwheel ring galaxy - Evidence for the first epoch of star formation in the wake of an expanding ring

Substantial near-infrared and optical color gradients have been observed within the disk of the classical ring galaxy A0035-324 (the 'Cartwheel'). The slope of the radial B-V, V-K color gradient coincides with starburst color evolution predictions, indicating that the age of the disk population is a function of radius. Observations indicate that the youngest stars are on the leading edge of the outer ring, and the stellar population grows increasingly older toward the nucleus. The radial color gradient, taken together with other evidence, strongly support the collisional picture for the formation of ring galaxies. Strong IR emission observed in the outer ring appears to be directly associated with the young star-forming regions rather than an evolved stellar population. We conclude that star formation is occurring in a mainly gaseous expanding density wave and that the Cartwheel is undergoing its first epoch of star formation since its formation.

Marcum, P. M.↗

Electrodeposition of metals and metal/cermet composites in low gravity

Electrodeposition experiments were carried out on the bench and a KC-135 aircraft at 0.01 g in anticipation of microgravity flights on NASA's Space Transportation System Shuttle. Experimental results obtained by interferometry compare concentration gradients as a function of time in the vicinity of a reducing electrode (cathode) for Cu(+2) and Co(+2) electrodeposition cells. No difference was found between bench and 0.01 g produced gradients for a .1M CuSO4 cell, but a significant difference was noted between the gradients in a 1M CoSO4 cell even though the bench cells were operated in a nonconvecting shielded (cathode over anode) mode. The gradient for Co(+2) depletion produced at 0.01 g was greater and the entire layer was thicker than found on the bench. Neutral buoyancy/matched density codeposition experiments were performed on the bench in an attempt to physically duplicate the results of metal/cermet codepositions in microgravity. Polystyrene spheres with average diameter 11.8 microns and density approximately matching that of 1M CoSO4 were utilized to emulate nonsedimenting cermets in microgravity. The cells were operated in a shielded convectionless mode. Comparison with literature data on codeposition with stirred cells indicate significant improvement in volume percent neutral occluded in the depositing metal matrix. A multicell electrodeposition flight apparatus that has been designed, constructed and is undergoing testing is discussed.

Riley, Clyde↗

Similarity considerations and conservation laws for magneto-static atmospheres

The equations of magnetohydrostatic equilibria for a plasma in a gravitational field are investigated analytically. For equilibria with one ignorable spatial coordinate, the equations reduce to a single nonlinear elliptic equation for the magnetic potential. Similarity solutions of the elliptic equation are obtained for the case of an isothermal atmosphere in a uniform gravitational field. The solutions are obtained from a consideration of the invariance group of the elliptic equation. The importance of symmetries of the elliptic equation also appears in the determination of conservation laws. It turns out that the elliptic equation can be written as a variational principle, and the symmetries of the variational functional lead (via Noether's theorem) to conservation laws for the equation. As an example of the application of the similarity solutions, a model magnetostatic atmosphere is constructed in which the current density J is proportional to the cube of the magnetic potential, and falls off exponentially with distance vertical to the base, with an 'e-folding' distance equal to the gravitational scale height. The solutions show the interplay between the gravitational force, the J x B force (B, magnetic field induction) and the gas pressure gradient.

Webb, G. M.↗

A hybrid simulation of contact discontinuity

Contact discontinuities in a collisionless plasma are studied by hybrid simulations, in which ions are treated as particles and electrons are considered as a fluid. It is demonstrated that contact discontinuity with a stable density ramp can exist in cases with a finite electron temperature. An electron pressure gradient is present across the contact discontinuity, leading to the presence of a parallel electric field and hence field-aligned potential increase (Delta Phi (sub parallel)) in the transition region. By reflecting ions at the discontinuity, this parallel electric potential peak reduces the interpenetration between hot and cold ions and maintains a stable density ramp across the contact discontinuity. The ratio of the field-aligned electric potential energy to ion thermal energy, e(Delta) Phi(sub parallel)/kT(sub i), is found to be an increasing function of T(sub e)/T(sub i), where T(sub e) and T(sub i) are respectively the electron and ion temperature.

Wu, B. H.↗

Adhesion of a bimetallic interface

The Hohenberg-Kohn and Kohn-Sham formalisms are used to examine binding (binding energy as a function of separation) for combinations of the simple metals Al(111), Zn(0001), Mg(0001), and Na(110) in contact. Similar metal contacts between Al, Zn, Mg, and Na are examined self-consistently in an ab initio calculation using the Kohn-Sham formalism. Crystallinity is included using the Aschroft pseudopotential via first order perturbation theory for the electron-ion interaction; and the ion-ion interaction is included exactly via a lattice sum. Binding energy was determined both in the local-density approximation and including gradient corrections to the exchange and correlation energy. Binding was found in all cases. In dissimilar metal contacts, interfacial bonding was greater than that in the weaker material predicting the possibility of metallic transfer. The nonzero position of the energy minimum in like metal contacts is explained in terms of consistency between the Ashcroft pseudopotential and the bulk charge density. Good agreement with experimental surface energies is obtained in the self-consistent calculation when nonlocal terms are included.

Ferrante, J.↗

A Simple Algebraic Grid Adaptation Scheme with Applications to Two- and Three-dimensional Flow Problems

An algebraic adaptive grid scheme based on the concept of arc equidistribution is presented. The scheme locally adjusts the grid density based on gradients of selected flow variables from either finite difference or finite volume calculations. A user-prescribed grid stretching can be specified such that control of the grid spacing can be maintained in areas of known flowfield behavior. For example, the grid can be clustered near a wall for boundary layer resolution and made coarse near the outer boundary of an external flow. A grid smoothing technique is incorporated into the adaptive grid routine, which is found to be more robust and efficient than the weight function filtering technique employed by other researchers. Since the present algebraic scheme requires no iteration or solution of differential equations, the computer time needed for grid adaptation is trivial, making the scheme useful for three-dimensional flow problems. Applications to two- and three-dimensional flow problems show that a considerable improvement in flowfield resolution can be achieved by using the proposed adaptive grid scheme. Although the scheme was developed with steady flow in mind, it is a good candidate for unsteady flow computations because of its efficiency.

Hsu, Andrew T.↗

Observations of neutral wind profiles between 115- and 175-km altitude in the dayside auroral oval

Lithium trail neutral wind measurements have been made in the morning dayside auroral oval as a function of altitude. During all four experiments the winds displayed a large altitude shear. In the most extreme case the wind at 135-km altitude was toward the south at 210 m/s, and at 185-km altitude toward the northwest at 170 m/s. A one-dimensional numerical solution to the Navier Stokes equation including an initial steady-state wind profile, viscosity, pressure gradient, Lorenz force and the Coriolis force was made in order to study the effect of auroral processes on the neutral atmosphere. Since electric fields and plasma density were measured, the Lorenz force could be modeled with reasonable accuracy. The pressure gradient was determined in an iterative way as a best fit to the observations. Within the limitations of a one-dimensional model the numerical solution is in reasonable agreement with one of the observations, which is applicable to such a modeling, and indicates that both Lorenz forces and pressure gradients due to auroral processes are required to explain the observation.

Pereira, E.↗

Potential Applications of Satellite Based Low Frequency Microwave Measurements of Snow and Ice

In this presentation we will examine some potential applications of low frequency microwave radar and radiometer data to remote sensing of snow and ice conditions over land. We focus on the following low frequency bands; L-band (1.28 GHz), S-band (2.69 GHz), and C-band (5.0 GHz). To predict the effect of snow cover on microwave emission, we need to characterize the dielectric properties of the snow in terms of important physical parameters such as snow density, temperature, and wetness. For the case of dry snow, we use an empirical mixing formula which gives the effective complex permittivity of snow as a function of density, temperature, and frequency. For wet snow, we use another mixing formula which gives the effective complex permittivity as a function of snow wetness. With dry snow, the loss at low microwave frequencies is very small so the corresponding penetration depths are very large (eg., 100 m). Clearly seasonal snow covers are far too thin to have a direct scattering or emission effect on these low frequency bands. There are, however, indirect effects introduced because of altered reflection at the snow/soil boundary, and an extra reflecting interface at the snow/air boundary. For C-band, layers of snow with different densities can have an impact if the number of layers grows sufficiently large. (eg., many meters of snow pack with cm-scale density layering.) Wet snow poses a more difficult problem. Liquid water is much more effective than ice at scattering and absorbing L-band radiation, so even a small amount of wetness will greatly reduce penetration. For example, with a snow wetness of 1 percent (volume fraction), and a density of 300 kg/cu. m, the penetration depths for L-band, S-band, and C-band are about 1.6 m, 0.3 m, and 0.15 m respectively. Very wet snow has more than 10 percent liquid water, and the corresponding L-band penetration is less than 0.2 m. Because the penetration distance is a strong function of frequency, it may be possible to identify a particular class of snow where the wetness is a few percent, and the depth is around one meter by looking at the frequency gradient.

West, Richard D.↗