SEARCH · Engineering Papers

Results for “Diffusion equation”

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,423 records · Page 79

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Viscosity and the Chew-Goldberger-Low equations in the solar corona

A general discussion of the dominant terms in the stress tensor in a magnetized plasma such as the solar corona is presented. The importance of dissipative terms such as electrical resistivity, heat conduction, and interspecies collisions is assessed. For average coronal conditions, the proton stress tensor is found to reduce to the dominant terms in the classical expression for the viscous stress. The classical expression can fail in the transition region, however. In the diffusion region of reconnection, classical viscosity will be appropriate if the resistivity is very large, so that the diffusion region is broad, but in that case the viscous heating is small compared to the resistive heating. On the other hand, the more general expression for the stress tensor is required if the diffusion region is thin; the stress tensor will be important in this case. The electron stress tensor is also considered, and it is shown how the classical expression for electron viscosity can fail in the transition region and lower corona.

Hollweg, J. V.↗

NASA NTRS1992

Segmented multigrid domain decomposition procedure for incompressible viscous flows

In this paper, the Navier-Stokes (NS) equations are approximated with a reduced NS system, that represents the lowest-order terms in an asymptotic Re expansion. This system allows for simplified boundary conditions, more generality in the location of the outflow boundary, and insures mass conservation in all subdomain grid interfaces, and also at the outflow boundary. The higher-order NS diffusion terms are included through a deferred corrector, in chosen subdomains, when required.

Srinivasan, Kumar↗

NASA NTRS1976

Modeling of short channel MOS transistors

Higher frequency response in MOS technology can be obtained by shortening the channel length. One approach for doing this involves an employment of higher resolution lithography technology. A second approach makes use of a double-diffused MOS transistor (DMOS). It is pointed out that the ordinary method of modeling the transistors used in both approaches is not accurate. An investigation is conducted of the questions which have to be considered for DMOS modeling. The modeling of a short channel MOS transistor is discussed, taking into account the derivation of the threshold voltage equation. Excellent agreement between theoretical and experimental data shows the accuracy of the described modeling approach.

Lin, H. C.↗

NASA NTRS1982

Design principles for high efficiency small-grain polysilicon solar cells, with supporting experimental studies

Design principles suggested here aim toward high conversion efficiency (greater than 15 percent) in polysilicon cells. The principles seek to decrease the liabilities of both intragranular and grain-boundary-surface defects. The advantages of a phosphorus atom concentration gradient in a thin (less than 50 microns) base of a p(+)/n(x)/n(+) drift-field solar cell, which produces favorable gradients in chemical potential, minority-carrier mobility and diffusivity, and recombination lifetime (via phosphorus gettering) are suggested. The degrading effects of grain boundaries are reduced by these three gradients and by substituting atoms (P, H, F or Li) for vacancies on the grain-boundary surface. From recent experiments comes support for the benefits of P diffusion down grain boundaries and, for quasi-grain-boundary-free and related structures. New analytic solutions for the n(x)-base include the effect of a power-law dependence between P concentration and lifetime. These provide an upper-bound estimate on the open circuit voltage. Finite-difference numerical solutions of the six Shockley equations furnish complete information about all solar-cell parameters and add insight concerning design.

Lindholm, F. A.↗

NASA NTRS1988

The high temperature creep deformation of Si3N4-6Y2O3-2Al2O3

The creep properties of silicon nitride containing 6 wt percent yttria and 2 wt percent alumina have been determined in the temperature range 1573 to 1673 K. The stress exponent, n, in the equation epsilon dot varies as sigma sup n, was determined to be 2.00 + or - 0.15 and the true activation energy was found to be 692 + or - 25 kJ/mol. Transmission electron microscopy studies showed that deformation occurred in the grain boundary glassy phase accompanied by microcrack formation and cavitation. The steady state creep results are consistent with a diffusion controlled creep mechanism involving nitrogen diffusion through the grain boundary glassy phase.

Todd, J. A.↗

NASA NTRS1989

The high temperature creep deformation of Si3N4-6Y2O3-2Al2O3

Todd, J. A.↗

NASA NTRS2001

Observations of O VI Emission from the Diffuse Interstellar Medium

We report the first Far Ultraviolet Spectroscopic Explorer (FUSE) measurements of diffuse O(VI) (lambda lambda 1032,1038) emission from the general diffuse interstellar medium outside of supernova remnants or superbubbles. We observed a 30 arcsec x 30 arcsec region of the sky centered at l = 315.0 deg and b = -41.3 deg. From the observed intensities (2930 +/- 290 (random) +/- 410 (systematic) and 1790 +/- 260 (random) +/- 250 (systematic) photons/sq cm/s/sr in the 1032 and 1038 angstrom emission lines, respectively), derived equations, and assumptions about the source location, we calculate the intrinsic intensity, electron density, thermal pressure, and emitting depth. The intensities are too large for the emission to originate solely in the Local Bubble. Thus, we conclude that the Galactic thick disk and lower halo also contribute. High velocity clouds are ruled out because there are none near the pointing direction. The calculated emitting depth is small, indicating that the O(VI)-bearing gas fills a small volume. The observations can also be used to estimate the cooling rate of the hot interstellar medium and constrain models. The data also yield the first intensity measurement of the C(II) 3s 2S(1/2) to 2p 2P(3/2) emission line at 1037 angstrom and place upper limits on the intensities of ultraviolet line emission from C(I), C(III), Si(II), S(III), S(IV), S(VI), and Fe(III).

Shelton, R. L.↗

NASA NTRS2003

Large Eddy Simulation Of Gravitational Effects In Transitional And Turbulent Gas-Jet Diffusion Flames

The influence of gravity on the spatial and the compositional structures of transitional and turbulent hydrocarbon diffusion flames are studies via large eddy simulation (LES) and direct numerical simulation (DNS) of round and planar jets. The subgrid-scale (SGS) closures in LES are based on the filtered mass density function (FMDF) methodology. The FMDF represents the joint probability density function (PDF) of the SGS scalars, and is obtained by solving its transport equation. The fundamental advantage of LES/FMDF is that it accounts for the effects of chemical reaction and buoyancy exactly. The methodology is employed for capturing some of the fundamental influences of gravity in equilibrium flames via realistic chemical kinetic schemes. Some preliminary investigation of the gravity effects in non-equilibrium flames is also conducted, but with idealized chemical kinetics models.

Jaberi, Farhad A.↗

NASA NTRS1981

An objective method for determining the generalized transport tensor for two-dimensional Eulerian models

An objective method for deriving the components of a generalized transport tensor for a two-dimensional model is presented. Representative meridional and vertical velocities and thermodynamic scalars at a uniform grid are used to reduce the problem to the solution of two flux equations for two unknowns. One unknown is the stream-function, coefficient of an antisymmetric tensor, which corrects the Eulerian mean motions for Stokes drift; the other is a time constant, which converts the deviatory velocity tensor to a symmetric transport tensor. The complete asymmetric tensor, called a transport tensor, has a divergence which yields both advection and diffusion by the deviatory velocities. Advantages and disadvantages of Lagrangian and Eulerian averages are discussed, and meridional-vertical velocity correlations are provided.

Danielsen, E. F.↗

NASA NTRS1990

Numerical prediction of turbulent flow over airfoil sections with a new nonequilibrium turbulence model

A new nonequilibrium turbulence closure model has been developed for computing wall bounded two-dimensional turbulent flows. This two-layer eddy viscosity model was motivated by the success of the Johnson-King model in separated flow regions. The influence of history effects are described by an ordinary differential equation developed from the turbulent kinetic energy equation. The performance of the present model has been evaluated by solving the flow around three airfoils using the Reynolds time-averaged Navier-Stokes equations. Excellent results were obtained for both attached and separated turbulent flows about the NACA 0012 airfoil, the RAE 2822 airfoil, and the Integrated Technology A 153W airfoil. Based on the comparison of the numerical solutions with the available experimental data, it is concluded that the new nonequilibrium turbulence model accurately captures the history effects of convection and diffusion on turbulence.

Ahmed, S.↗

NASA NTRS1993

Evaporation of LOX under supercritical and subcritical conditions

The evaporation of LOX under supercritical and subcritical conditions was studied experimentally and theoretically. In experiments, the evaporation rate and surface temperature were measured for LOX strand vaporizing in helium environments at pressures ranging from 5 to 68 atmospheres. Gas sampling and chromatography analysis were also employed to profile the gas composition above the LOX surface for the purpose of model validation. A comprehensive theoretical model was formulated and solved numerically to simulate the evaporation process of LOX at high pressures. The model was based on the conservation equations of mass, momentum, energy, and species concentrations for a multicomponent system, with consideration of gravitational body force, solubility of ambient gases in liquid, and variable thermophysical properties. Good agreement between predictions and measured oxygen mole fraction profiles was obtained. The effect of pressure on the distribution of the Lewis number, as well as the effect of variable diffusion coefficient, were further examined to elucidate the high-pressure transport behavior exhibited in the LOX vaporization process.

Yang, A. S.↗

NASA NTRS1993

Ion acceleration in impulsive solar flares

Nonrelativistic spectra of protons and ions accelerated in impulsive solar flares are derived using more realistic turbulence power spectra. The calculation is based on a particle transport equation extracted from a second step acceleration model containing stochastic acceleration. The turbulence model is generalized to waves with a small angle to the magnetic field vector and to turbulence power spectra with spectral indices s smaller than 2. Due to the occurrence of impulsive flares at low coronal heights, Coulomb losses at the dense coronal plasma and diffusive particle escape are taken into account. The ion spectra show deviations from long-duration spectra near the Coulomb barrier, where the losses become maximal. The Z-squared/A-dependence of the Coulomb losses leads to spectral variations for different ions. We present a method to estimate the turbulence parameters and injection conditions of the flare particles using ion ratios like Fe/O of impulsive flares.

Steinacker, Jurgen↗

NASA NTRS2010

Ares I Reaction Control System Propellant Feedline Decontamination Modeling

The objective of the work presented here is to quantify the effects of purge gas temperature, pressure, and mass flow rate on Hydrazine (Hz) decontamination rates of the Ares I Roll Control System and Reaction Control System. A survey of experts in this field revealed the absence of any decontamination rate prediction models. Three basic decontamination methods were identified for analysis and modeling. These include low pressure eduction, high flow rate purge, and pulse purge. For each method, an approach to predict the Hz mass transfer rate, as a function of system pressure, temperature, and purge gas mass flow rate, is developed based on the applicable physics. The models show that low pressure eduction is two orders of magnitude more effective than the high velocity purge, which in turn is two orders of magnitude more effective than the pure diffusion component of pulse purging of deadheads. Eduction subjects the system to low pressure conditions that promote the extraction of Hz vapors. At 120 F, Hz is saturated at approximately 1 psia. At lower pressures and 120 F, Hz will boil, which is an extremely efficient means to remove liquid Hz. The Hz boiling rate is predicted by equating the rate at which energy is added to the saturated liquid Hz through heaters at the tube outer wall with the energy removed from the liquid through evaporation. Boil-off fluxes were predicted by iterating through the range of local pressures with limits set by the minimum allowed pressure of 0.2 psia and maximum allowed wall temperature of 120 F established by the heaters, which gives a saturation pressure of approximately 1.0 psia. Figure 1 shows the resulting boil-off fluxes as a function of local eduction pressure. As depicted in figure 1, the flux is a strong inverse function of eduction pressure, and that minimizing the eduction pressure maximizes the boil-off flux. Also, higher outer wall temperatures lead to higher boil-off fluxes and allow for boil-off over a greater range of eduction pressures.

Pasch, James J.↗

NASA NTRS1983

The nonlocal theory of periodic density drift instabilities

The nonlocal theory of electrostatic density drift instabilities is developed for an arbitrary, periodic variation in plasma density. Linear Vlasov theory is applied to Fourier series representations of the plasma density, and the resulting equations are solved numerically. The growth rates of the universal drift instability are compared by using local and nonlocal theories. It is demonstrated that the local theory is inadequate if the nonlocal eigenfunctions of the electrostatic modes have significant amplitude over spatial regions wider than the regions of large plasma density gradients. Nonlocal effects can enhance wave-particle diffusion of a triangular-shaped variation and reduce wave-particle dissipation of a square-shaped irregularity.

Bernhardt, P. A.↗

NASA NTRS1995

An Extended Lagrangian Method

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method,' is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. The present method and the Arbitrary Lagrangian-Eulerian (ALE) method have a similarity in spirit-eliminating the cross-streamline numerical diffusion. For this purpose, we suggest a simple grid constraint condition and utilize an accurate discretization procedure. This grid constraint is only applied to the transverse cell face parallel to the local stream velocity, and hence our method for the steady state problems naturally reduces to the streamline-curvature method, without explicitly solving the steady stream-coordinate equations formulated a priori. Unlike the Lagrangian method proposed by Loh and Hui which is valid only for steady supersonic flows, the present method is general and capable of treating subsonic flows and supersonic flows as well as unsteady flows, simply by invoking in the same code an appropriate grid constraint suggested in this paper. The approach is found to be robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

Liou, Meng-Sing↗

NASA NTRS2004

The Liquid-Water Oscillation in Modeling Boundary-Layer Cumuli with Third-Order Turbulence Closure Models

A hierarchy of third-order turbulence closure models are used to simulate boundary-layer cumuli in this study. An unrealistically strong liquid-water oscillation (LWO) is found in the fully prognostic model, which predicts all third moments. The LWO propagates from cloud base to cloud top with a speed of 1 m/s. The period of the oscillation is about 1000 s. Liquid-water buoyancy terms in the third-moment equations contribute to the LWO. The LWO mainly affects the vertical profiles of cloud fraction, mean liquid-water mixing ratio and the fluxes of liquid-water potential temperature and total water, but has less impact on the vertical profiles of other second-moments and third-moments. In order to minimize the LWO, a moderate large diffusion coefficient and a large turbulent dissipation at its originating level are needed. However, this approach distorts the vertical distributions of cloud fraction and liquid-water mixing ratio. A better approach is to parameterize liquid-water buoyancy more reasonably. A minimally prognostic model, which diagnoses all third moments except for vertical velocity, is shown to produce better results, compared to a fully prognostic model.

Cheng, A.↗

NASA NTRS1974

Observations and calculations of atmospheric helium 10830 emission

Observations and computer model calculations of helium excitation and transport were carried out. The following conclusions can be drawn: (1) the dawn helium density in the thermosphere is greater than the dusk helium density by 50 per cent or more, in contrast to the static diffusion models but in agreement with the dynamic models; (2) the computed He 10830 intensities are inconsistent with the observations when use of the new Penning ionization cross section for loss of atomic oxygen is used; (3) conjugate photoelectron precipitation is responsible for excitation of atmospheric helium at middle latitudes, but it is not evident in the observations at the equator or Brazil (L equals 1.1); (4) the computations predict greater conjugate effects than observed, leading to the possibility of strong interhemispheric absorption.

Christensen, A. B.↗

NASA NTRS1984

Turbulent solutions of the equations of fluid motion

Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. Initial three-dimensional cosine velocity fluctuations and periodic boundary conditions are used in most of the work considered. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations, the initially nonrandom flow develops into an apparently random turbulence. Thus randomness or turbulence can arise as a consequence of the structure of the Navier-Stokes equations. The cases considered include turbulence which is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A mean shear is present in some cases. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components, and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.

Deissler, R. G.↗