Engineering Papers⌕ Search

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.

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

At least 901 records · Page 50

Test-particle method in kinetic theory of a plasma.

The introduction of a test particle into a system is considered. The system may be described by the Born-Bogoliubov-Green-Kirkwood-Yvon hierarchy. The field particles form a cloud which surrounds the test particle. The cloud is described by a conditional probability function which satisfies a certain equation. A generalization of the superposition principle reported by Rostoker (1964) to higher order correlation functions is discussed. Kinetic equations with the generalized Lenard-Balescu term are obtained, taking into account also diffusion by waves. The characteristics regarding the absorption or emission of waves by particles can be calculated.

Matsuda, K.↗

Design and calibration of a rocket-borne electron spectrometer for investigation of particle ionization in the nighttime midlatitude E region

An explanation was developed for the formation, near midnight at midlatitudes, of a broad electron density layer extending approximately from 120 to 180 km and usually referred to as the intermediate E layer. The responsible mechanism is believed to be the converging vertical ion drifts resulting from winds of the solar semidiurnal tide. Numerical solutions of the continuity equation appropriate to the intermediate layer is described for particular models of ion drift, diffusion coefficents, and ionization production. Analysis of rocket observations of the layer show that the ionization rate is highly correlated with the planetary geomagnetic index, K sub p. Particle flux measurements support the idea that energetic electrons are the principal source of this ionization. A semiconductor spectrometer experiment for investigation of the particle flux, spectrum, and angular properties was designed and successfully flown on a Nike Apache rocket. A detailed description of the theory, design, and calibration of the experiment and some preliminary results presented.

Voss, H. D.↗

Two-stream theory of spectral reflectance of snow

Spectral reflectance of snow under diffuse illumination is studied using the two-stream approximation of the radiative transfer equation. The scattering and absorption within the snowcover due to the randomly distributed ice grains are characterized by the single scattering albedo and anisotropic phase function. Geometric optics calculations are used to relate the scattering and absorption parameters to grain size and density of snow. Analytical expressions for the intensity within the snowpack and the asymptotic flux extinction coefficient are also obtained. Good agreement is shown between the theory and available experimental data on visible and near-infrared reflectance and asymptotic flux extinction coefficient. The theory also may be used to explain the observed effect of aging on the snow reflectance.

Choudhury, B. J.↗

Two-stream theory of reflectance of snow

Spectral reflectance of snow under diffuse illumination is studied using the two-stream approximation of the radiative transfer equation. The scattering and absorption parameters of the radiative transfer equation - the single scattering albedo, the optical depth, and the integrated phase function are obtained from the grain size and density of snow. Analytical expressions for the intensity within the snowpack, the reflectance, and the asymptotic flux extinction coefficient, are given. Good agreement is shown between the theory and available experimental data on visible and near-infrared reflectance, and the asymptotic flux extinction coefficient. The theory may also be used to explain the observed effect of aging on the snow reflectance.

Choudhury, B. J.↗

Viscous, resistive MHD stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one dimensional magnetohydrodynamic (MHD) quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds like numbers. Marginal stability curves, growth rates versus Reynolds like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result which appears general is that instability was found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

Simulations of magnetic-flux transport in solar active regions

The evolution of several observed solar active regions is simulated by solving a transport equation for magnetic flux at the photosphere. The rates of rotation, meridional flow, and diffusion of the flux are determined self-consistently in the calculations. The findings are in good quantitative agreement with previous measures of the rotation rate and diffusion constant associated with photospheric magnetic fields. Although the meridional velocities are consistent in direction and magnitude with recently reported poleward flows, relatively large uncertainties in the velocity determinations make this result inconclusive.

Devore, C. R.↗

A method for calculating turbulent boundary layers and losses in the flow channels of turbomachines

An interactive inviscid core flow-boundary layer method is presented for the calculation of turbomachine channel flows. For this method, a one-dimensional inviscid core flow is assumed. The end-wall and blade surface boundary layers are calculated using an integral entrainment method. The boundary layers are assumed to be collateral and thus are two-dimensional. The boundary layer equations are written in a streamline coordinate system. The streamwise velocity profiles are approximated by power law profiles. Compressibility is accounted for in the streamwise direction but not in the normal direction. Equations are derived for the special cases of conical and two-dimensional rectangular diffusers. For these cases, the assumptions of a one-dimensional core flow and collateral boundary layers are valid. Results using the method are compared with experiment and good quantitative agreement is obtained.

Schumann, Lawrence F.↗

A method for calculating turbulent boundary layers and losses in the flow channels of turbomachines

An interactive inviscid core flow-boundary layer method is presented for the calculation of turbomachine channel flows. For this method, a one-dimensional inviscid core flow is assumed. The end-wall and blade surface boundary layers are calculated using an integral entrainment method. The boundary layers are assumed to be collateral and thus are two-dimensional. The boundary layer equations are written in a streamline coordinate system. The streamwise velocity profiles are approximated by power law profiles. Compressibility is accounted for in the streamwise direction but not in the normal direction. Equations are derived for the special cases of conical and two-dimensional rectangular diffusers. For these cases, the assumptions of a one-dimensional core flow and collateral boundary layers are valid. Results using the method are compared with experiment and good quantitative agreement is obtained.

Schumann, Lawrence F.↗

The computation of non-equilibrium chemically-reacting flows

Numerical techniques for the analysis of flows with nonequilibrium thermodynamics and chemistry are developed and demonstrated. Particular attention is given to the formulation of the thermodynamic and chemical gas models, the derivation of the governing equations, and the solution procedures. Results for typical problems involving shock tubes and supersonic diffusers are presented in extensive graphs and briefly characterized.

Grossman, B.↗

Cosmic-ray streaming in the Born approximation

The present work invokes the Born approximation to derive a more accurate expression for the streaming of cosmic rays parallel to the mean magnetic field. While all prior results pertaining to the helicity dependence of the diffusion coefficient and convection speed can be recovered as special cases from this streaming equation, it is concluded that a new set of transport parameters presented here is more appropriate for the solar modulation of galactic cosmic rays. In addition, a new parameter related to time variability, which may be a dominant cause of charge sign-dependent transport of solar particles, is introduced.

Bieber, J. W.↗

A time-accurate implicit method for chemically reacting flows at all Mach numbers

The objective of this work is to develop a unified solution algorithm capable of treating time-accurate chemically reacting flows at all Mach numbers, ranging from molecular diffusion velocities to supersonic speeds. A rescaled pressure term is used in the momentum equation to circumvent the singular behavior of pressure at low Mach numbers. A dual time-stepping integration procedure is established. The system eigenvalues become well behaved and have the same order of magnitude, even in the very low Mach number regime. The computational efficiency for moderate and high speed flow is competitive with the conventional density-based scheme. The capabilities of the algorithm are demonstrated by applying it to selected model problems including nozzle flows and flame dynamics.

Withington, J. P.↗

Simulation of the present-day atmospheric ozone, odd nitrogen, chlorine and other species using a coupled 2-D model in isentropic coordinates

A zonally averaged two-dimensional model for chemical species in the atmosphere (Olaguer et al., 1990), with coupled dynamics, radiation transfer, and chemistry, is used to study the present-day atmospheric ozone, odd nitrogen, chlorine, other species. The model utilizes all zonally averaged physical equations of momentum, energy, and mass and determines self-consistently the advective and diffusive tranasport parameters from the temperature specific to the period of observation. A comparison of the model results with observations showed good agreement.

Yang, H.↗

The implementation and validation of improved landsurface hydrology in an atmospheric general circulation model

Landsurface hydrological parameterizations are implemented in the NASA Goddard Institute for Space Studies (GISS) General Circulation Model (GCM). These parameterizations are: (1) runoff and evapotranspiration functions that include the effects of subgrid scale spatial variability and use physically based equations of hydrologic flux at the soil surface, and (2) a realistic soil moisture diffusion scheme for the movement of water in the soil column. A one dimensional climate model with a complete hydrologic cycle is used to screen the basic sensitivities of the hydrological parameterizations before implementation into the full three dimensional GCM. Results of the final simulation with the GISS GCM and the new landsurface hydrology indicate that the runoff rate, especially in the tropics is significantly improved. As a result, the remaining components of the heat and moisture balance show comparable improvements when compared to observations. The validation of model results is carried from the large global (ocean and landsurface) scale, to the zonal, continental, and finally the finer river basin scales.

Johnson, Kevin D.↗

Mixing and mixing enhancement in supersonic reacting flowfields

Propulsion systems planned for use late in this century and beyond will require appropriate physical models for describing supersonic combustion and numerical techniques for solving the model governing equations. A computer program to study these flows is reported which considers the multicomponent diffusion and convection of important chemical species, the finite-rate reaction of these species, and the resulting interaction of the field mechanics and the chemistry. The application of the program to a spatially developing and reacting mixing layer, which serves an an excellent physical model for the mixing and reaction processes that take place in a scramjet combustor, is reported. Several techniques to enhance the fuel-air mixing and growth of that layer and improve its overall combustion efficiency are considered.

Drummond, J. P.↗

The implementation and validation of improved land-surface hydrology in an atmospheric general circulation model

New land-surface hydrologic parameterizations are implemented into the NASA Goddard Institute for Space Studies (GISS) General Circulation Model (GCM). These parameterizations are: 1) runoff and evapotranspiration functions that include the effects of subgrid-scale spatial variability and use physically based equations of hydrologic flux at the soil surface and 2) a realistic soil moisture diffusion scheme for the movement of water and root sink in the soil column. A one-dimensional climate model with a complete hydrologic cycle is used to screen the basic sensitivities of the hydrological parameterizations before implementation into the full three-dimensional GCM. Results of the final simulation with the GISS GCM and the new land-surface hydrology indicate that the runoff rate, especially in the tropics, is significantly improved. As a result, the remaining components of the heat and moisture balance show similar improvements when compared to observations. The validation of model results is carried from the large global (ocean and land-surface) scale to the zonal, continental, and finally the regional river basin scales.

Johnson, Kevin D.↗

Forced Convection and Sedimentation Past a Flat Plate

The steady laminar flow of a well-mixed suspension of monodisperse solid spheres, convected steadily past a horizontal flat plate and sedimenting under the action of gravity, is examined. It is shown that, in the limit as Re approaches infinity and epsilon approaches 0, where Re is the bulk Reynolds number and epsilon is the ratio of the particle radius a to the characteristic length scale L, the analysis for determining the particle concentration profile has several aspects in common with that of obtaining the temperature profile in forced-convection heat transfer from a wall to a fluid stream moving at high Reynolds and Prandtl numbers. Specifically, it is found that the particle concentration remains uniform throughout the O(Re(exp -1/2)) thick Blasius boundary layer except for two O(epsilon(exp 2/3)) thin regions on either side of the plate, where the concentration profile becomes non-uniform owing to the presence of shear-induced particle diffusion which balances the particle flux due to convection and sedimentation. The system of equations within this concentration boundary layer admits a similarity solution near the leading edge of the plate, according to which the particle concentration along the top surface of the plate increases from its value in the free stream by an amount proportional to X(exp 5/6), with X measuring the distance along the plate, and decreases in a similar fashion along the underside. But, unlike the case of gravity settling on an inclined plate in the absence of a bulk flow at infinity considered earlier, here the concentration profile remains continuous everywhere. For values of X beyond the region near the leading edge, the particle concentration profile is obtained through the numerical solution of the relevant equations. It is found that, as predicted from the similarity solution, there exists a value of X at which the particle concentration along the top side of the plate attains its maximum value phi(sub m) and that, beyond this point, a stagnant sediment layer will form that grows steadily in time. This critical value of X is computed as a function of phi(sub s), the particle volume fraction in the free stream. In contrast, but again in conformity with the similarity solution, for values of X sufficiently far removed from the leading edge along the underside of the plate, a particle-free region is predicted to form adjacent to the plate. This model, with minor modifications, can be used to describe particle migration in other shear flows, as, for example, in the case of crossflow microfiltration.

Pelekasis, Nikolaos A.↗

A Numerical Method for Obtaining Monoenergetic Neutron Flux Distributions and Transmissions in Multiple-Region Slabs

This method is investigated for semi-infinite multiple-slab configurations of arbitrary width, composition, and source distribution. Isotropic scattering in the laboratory system is assumed. Isotropic scattering implies that the fraction of neutrons scattered in the i(sup th) volume element or subregion that will make their next collision in the j(sup th) volume element or subregion is the same for all collisions. These so-called "transfer probabilities" between subregions are calculated and used to obtain successive-collision densities from which the flux and transmission probabilities directly follow. For a thick slab with little or no absorption, a successive-collisions technique proves impractical because an unreasonably large number of collisions must be followed in order to obtain the flux. Here the appropriate integral equation is converted into a set of linear simultaneous algebraic equations that are solved for the average total flux in each subregion. When ordinary diffusion theory applies with satisfactory precision in a portion of the multiple-slab configuration, the problem is solved by ordinary diffusion theory, but the flux is plotted only in the region of validity. The angular distribution of neutrons entering the remaining portion is determined from the known diffusion flux and the remaining region is solved by higher order theory. Several procedures for applying the numerical method are presented and discussed. To illustrate the calculational procedure, a symmetrical slab ia vacuum is worked by the numerical, Monte Carlo, and P(sub 3) spherical harmonics methods. In addition, an unsymmetrical double-slab problem is solved by the numerical and Monte Carlo methods. The numerical approach proved faster and more accurate in these examples. Adaptation of the method to anisotropic scattering in slabs is indicated, although no example is included in this paper.

Schneider, Harold↗

Effects of non-unity Lewis numbers in diffusion flames

The purpose of this work is to carry out direct numerical simulations of diffusion controlled combustion with non-unity Lewis numbers for the reactants and products, thus accounting for the differential diffusion effects of the temperature and concentration fields. We use a formulation based on combining the conservation equations in a way to eliminate the reaction terms similar to the method used by Burke and Schumann (1928) for unity Lewis numbers. We present calculations for an axisymmetric fuel jet and for a planar, time evolving mixing layer, leaving out the effects of thermal expansion and variations of the transport coefficients due to the heat release. Our results show that the front of the flame shifts toward the fuel or oxygen sides owing to the effect of the differential diffusion and that the location of maximum temperature may not coincide with the flame. The dependence of the distribution of the reaction products on their Lewis number has been investigated.

Linan, A.↗