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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Vapor Transport Within the Thermal Diffusion Cloud Chamber

A review of the equations used to determine the 1-D vapor transport in the thermal diffusion cloud chamber (TDCC) is presented. These equations closely follow those of the classical Stefan tube problem in which there is transport of a volatile species through a noncondensible, carrier gas. In both cases, the very plausible assumption is made that the background gas is stagnant. Unfortunately, this assumption results in a convective flux which is inconsistent with the momentum and continuity equations for both systems. The approximation permits derivation of an analytical solution for the concentration profile in the Stefan tube, but there is no computational advantage in the case of the TDCC. Furthermore, the degree of supersaturation is a sensitive function of the concentration profile in the TD CC and the stagnant background gas approximation can make a dramatic difference in the calculated supersaturation. In this work, the equations typically used with a TDCC are compared with very general transport equations describing the 1-D diffusion of the volatile species. Whereas no pressure dependence is predicted with the typical equations, a strong pressure dependence is present with the more general equations given in this work. The predicted behavior is consistent with observations in diffusion cloud experiments. It appears that the new equations may account for much of the pressure dependence noted in TDCC experiments, but a comparison between the new equations and previously obtained experimental data are needed for verification.

Ferguson, Frank T.↗
NASA NTRS1988

A closure model of diffusion transport of the Reynolds-stress equations and its application to a turbulent step flow

A study is made for the development and computations of the separating and reattaching shear flows. The highlight of the study is an attempt to predict the third-moments of turbulent velocity which is responsible for the diffusion transport of the Reynolds stresses. The present computations show that the third-moments obtained by employing the low-Reynolds number model of transport equations improve the prediction of the third-moments.

Amano, R. S.↗
NASA NTRS1983

Modeling of Fluctuating Mass Flux in Variable Density Flows

The approach solves for both Reynolds and Favre averaged quantities and calculates the scalar pdf. Turbulent models used to close the governing equations are formulated to account for complex mixing and variable density effects. In addition, turbulent mass diffusivities are not assumed to be in constant proportion to turbulent momentum diffusivities. The governing equations are solved by a combination of finite-difference technique and Monte-Carlo simulation. Some preliminary results on simple variable density shear flows are presented. The differences between these results and those obtained using conventional models are discussed.

So, R. M. C.↗
NASA NTRS1994

Transportation of volatile elements in thermally evolving planetesimals: An important role of metallic iron

Ordinary chondrites are considered to have experienced thermal metamorphism in small bodies. We are interested in behaviors of volatile elements in such a kind of thermally evolving planetesimals. Volatile elements generally have high vapor pressures at high temperature. In porous bodies, with a high gas permeability, volatile elements are transported efficiently over a long range. Behavior of volatile elements transported by permeable gas flow can be handled by an equation whose form is similar to that of the equation of thermal diffusion. We can follow transportation of heats and volatile elements in planetesimals, when parameters in these equations, initial conditions and chemical behavior of volatile elements are given. Recently, we discovered that nitrogen in equilibrated H-chondrites is mainly trapped in taenite (f.c.c. Fe-Ni), probably dissolved in interstitial sites. Fegley suggests that metallic iron cannot trap nitrogen in the solar nebula gas due to its very low nitrogen partial pressure. Approximately 1 bar of nitrogen pressure is required to explain the nitrogen content in taenite. We may expect high nitrogen gas partial pressure (possibly produced by vaporization of nitrogen-bearing solids such as organic materials) at the interior of thermally evolving planetesimals. Kinetic behavior of nitrogen in taenite suggests that it can easily be equilibrated with the ambient nitrogen gas at temperatures of approximately 500 C or higher. We consider that nitrogen is trapped in taenite through a nitrogen redistribution process occurred during the thermal metamorphic event.

Hashizume, K.↗
NASA NTRS1992

New explicit equations for the accurate calculation of the growth and evaporation of hydrometeors by the diffusion of water vapor

The traditional explicit growth equation has been widely used to calculate the growth and evaporation of hydrometeors by the diffusion of water vapor. This paper reexamines the assumptions underlying the traditional equation and shows that large errors (10-30 percent in some cases) result if it is used carelessly. More accurate explicit equations are derived by approximating the saturation vapor-density difference as a quadratic rather than a linear function of the temperature difference between the particle and ambient air. These new equations, which reduce the error to less than a few percent, merit inclusion in a broad range of atmospheric models.

Srivastava, R. C.↗
NASA NTRS1983

Accuracy of approximations to the Navier-Stokes equations

The systems of truncated differential equations that have been proposed to reduce the complexity and large computational costs of solutions to the full Navier-Stokes equations are considered. These systems are computationally efficient and capture all the physically relevant behavior. The systems follow a certain hierarchy: (1) the classical boundary-layer equations with specified edge properties (usually the streamwise pressure distribution); (2) the coupled boundary-layer/inviscid equations; (3) the so-called thin-layer equations that discard streamwise diffusion; and (4) the Navier-Stokes equations. Consideration is given to each of these approximations applied to an incompressible, laminar-separating flow at low and moderate Reynolds numbers. It is pointed out that for any flow or region of flow for which viscous-inviscid interaction effects are small, classical boundary-layer equations will provide a satisfactory description of the viscous flow at a fraction of the computational cost of any higher approximations.

Murphy, J. D.↗
NASA NTRS1994

The Prediction of Nozzle Performance and Heat Transfer in Hydrogen/Oxygen Rocket Engines with Transpiration Cooling, Film Cooling, and High Area Ratios

An advanced engineering computational model has been developed to aid in the analysis of chemical rocket engines. The complete multispecies, chemically reacting and diffusing Navier-Stokes equations are modelled, including the Soret thermal diffusion and Dufour energy transfer terms. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film-cooled nozzle, and a transpiration-cooled plug-and-spool rocket engine. The results indicate that the thrust coefficient predictions of the 1030:1 nozzle and the film-cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. It is demonstrated that thermal diffusion has a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle and was shown to represent a significant fraction of the diffusion fluxes occurring in the transpiration-cooled rocket engine.

Kacynski, Kenneth J.↗
NASA NTRS1994

On Bi-Grid Local Mode Analysis of Solution Techniques for 3-D Euler and Navier-Stokes Equations

A procedure is presented for utilizing a bi-grid stability analysis as a practical tool for predicting multigrid performance in a range of numerical methods for solving Euler and Navier-Stokes equations. Model problems based on the convection, diffusion and Burger's equation are used to illustrate the superiority of the bi-grid analysis as a predictive tool for multigrid performance in comparison to the smoothing factor derived from conventional von Neumann analysis. For the Euler equations, bi-grid analysis is presented for three upwind difference based factorizations, namely Spatial, Eigenvalue and Combination splits, and two central difference based factorizations, namely LU and ADI methods. In the former, both the Steger-Warming and van Leer flux-vector splitting methods are considered. For the Navier-Stokes equations, only the Beam-Warming (ADI) central difference scheme is considered. In each case, estimates of multigrid convergence rates from the bi-grid analysis are compared to smoothing factors obtained from single-grid stability analysis. Effects of grid aspect ratio and flow skewness are examined. Both predictions are compared with practical multigrid convergence rates for 2-D Euler and Navier-Stokes solutions based on the Beam-Warming central scheme.

Ibraheem, S. O.↗
NASA NTRS1995

Concentration Dependence of Solution Shear Viscosity and Solute Mass Diffusivity in Crystal Growth from Solutions

The physical properties of a supersaturated binary solution such as its density rho, shear viscosity eta, and solute mass diffusivity D are dependent on the solute concentration c: rho = rho(c), eta = eta(c), and D = D(c). The diffusion boundary layer equations related to crystal growth from solution are derived for the case of natural convection with a solution density, a shear viscosity, and a solute diffusivity that are all depen- dent on solute concentration. The solution of these equations has demonstrated the following. (1) At the vicinity of the saturation concentration c(sub s) the solution shear viscosity eta depends on rho as eta(sub s) = eta(rho(sub s))varies as square root of rho(c(sub s)). This theoretically derived result has been verified in experiments with several aqueous solutions of inorganic and organic salts. (2) The maximum solute mass transfer towards the growing crystal surface can be achieved for values of c where the ratio of d ln(D(c)/dc) to d ln(eta(c)/dc) is a maximum.

Izmailov, Alexander F.↗
NASA NTRS1972

Modulation of low energy cosmic rays

The power spectrum of the interplanetary magnetic field were tested, as well as the power spectrum-cosmic diffusion ray coefficient relation at low energies. A first order perturbation solution of the Fokker-Planck equation governing the diffusion, convection, and adiabatic deceleration of galactic cosmic rays in the solar medium was found to relate intensity fluctuations of low energy cosmic rays to local changes in the propagation parameters. Diffusion coefficients and their day to day variations were calculated from interplanetary magnetic field data obtained by Pioneer 6 in 1965/1966. These are compared to simultaneous observations by IMP 3 of the proton flux in three energy channels (20 to 40, 40 to 60, 60 to 80 MeV).

Sari, J. W.↗
NASA NTRS1990

Compact scheme for systems of equations applied to fundamental problems of mechanics of continua

Compact scheme formulation was used in the treatment of boundary conditions for a system of coupled diffusion and Poisson equations. Models and practical solutions of specific engineering problems arising in solid mechanics, chemical engineering, heat transfer and fuid mechanics are described and analyzed for efficiency and accuracy. Only 2-D cases are discussed and a new method of numerical treatment of boundary conditions common in the fundamental problems of mechanics of continua is presented.

Klimkowski, Jerzy Z.↗
NASA NTRS1981

A computer program for the simulation of heat and moisture flow in soils

A computer program that simulates the flow of heat and moisture in soils is described. The space-time dependence of temperature and moisture content is described by a set of diffusion-type partial differential equations. The simulator uses a predictor/corrector to numerically integrate them, giving wetness and temperature profiles as a function of time. The simulator was used to generate solutions to diffusion-type partial differential equations for which analytical solutions are known. These equations include both constant and variable diffusivities, and both flux and constant concentration boundary conditions. In all cases, the simulated and analytic solutions agreed to within the error bounds which were imposed on the integrator. Simulations of heat and moisture flow under actual field conditions were also performed. Ground truth data were used for the boundary conditions and soil transport properties. The qualitative agreement between simulated and measured profiles is an indication that the model equations are reasonably accurate representations of the physical processes involved.

Camillo, P.↗
NASA NTRS1990

Radial diffusion of low-energy plasma ions in Saturn's magnetosphere

Radial diffusion of low-energy plasma ions in Saturn's magnetosphere is investigated using a comprehensive set of equations for radial diffusion that incorporate distributed sources and sinks of ions. The results of calculations indicate that the radial-diffusion transport of low-energy O(+) ions with a source in the neutral H2O cloud of the satellites Dione and Tethys can account for Voyager observations of thermal heavy ions in Saturn's magnetosphere. The source rate was calculated to be about 10 to the 26th O(+) ions/sec, in good agreement with the sputtering calculations of Johnson et al. (1989). It is estimated that, due to fast radial diffusion, the residence time of O(+) ions in the Dione-Tethys torus is about 30 days, sufficiently short to account for the plasma density observed there. The densities of hot H(+) and N(+) resulting from the ionization and pickup of Titan's neutral clouds in the outer magnetosphere can also be accounted for within the framework of diffusive ion transport.

Barbosa, D. D.↗
NASA NTRS1973

The application of far IR lasers to the investigation of the optical properties of materials

The effect of diffusion on the saturation and absorption properties of molecular gases is considered. The analysis starts with a simple derivation of the multi-level rate equations. Incorporating the diffusion term, calculations are carried out both for the homogeneous and inhomogeneous broadening cases. Both graphical and analytical solutions are obtained for the absorption coefficient and the saturation intensity. It is shown that the effect of diffusion is characterized by a diffusion function, which varies inversely with the radius of the beam at low pressures and small beam sizes. The analysis is done in terms of two parameters k1a and k2a, which are simple functions of the intensity of the beam I, the effective relaxation rate w, and the effective diffusion constant D the a is designated as the radius of the beam.

Stafsudd, O. M.↗
NASA NTRS1986

Vapor transport processes of Ge-chalcogenides under normal and reduced gravity conditions

Experimental and theoretical studies of the GeSe-GeI4 system revealed the multicomponent, multireaction nature of this system and the existence of a diffusion boundary layer; they also yielded equations for the prediction of diffusion limited mass transport rates. Vapor transport and crystal growth experiments of the GeSe-GeI4 and GeSe-Xe systems on earth and in a microgravity environment demonstrated the effects of convection and of microgravity on mass flux and crystal morphology. The surface and bulk morphology of GeSe crystals obtained for both systems in microgravity is considerably improved relative to ground-based specimens. Unexpected crystal growth phenomena were observed for the GeSe-Xenon system on recent Shuttle flights.

Wiedemeier, H.↗
NASA NTRS1989

Evolution of the sun's polar fields during sunspot cycle 21 - Poleward surges and long-term behavior

Longitudinally averaged observations of the photospheric field during 1976-1986 are analyzed using a flux transport model. The way in which source eruptions, supergranular diffusion, and meridional flow collaborate to produce strong, highly concentrated polar fields near sunspot minimum is clarified as follows: (1) widespread eruptions of individual bipolar magnetic regions, with their leading polarity flux equatorward of their trailing polarity flux, collectively establish a large-scale separation of polarities in latitude; (2) the low-latitude, leading polarity flux diffuses across the equator and merges with its opposite hemisphere counterpart; and (3) meridional flow carries the resulting surplus of trailing polarity flux to the poles, and concentrates it there against the spreading effect of diffusion. Episodic 'surges' of flux to the poles are induced by fluctuations in the source eruption rate. Simulations indicate that relatively weak, trailing polarity surges may occur even in a steady flow field. However, in order to account for the giant surges of alternating polarity and the resulting oscillations in the polar field strength observed during 1980-1982, both accelerated flow and enhanced eruption rates are required.

Wang, Y.-M.↗
NASA NTRS1986

Convective diffusion in protein crystal growth

A protein crystal modeled as a flat plate suspended in the parent solution, with the normal to the largest face perpendicular to gravity and the protein concentration in the solution adjacent to the plate taken to be the equilibrium solubility, is studied. The Navier-Stokes equation and the equation for convective diffusion in the boundary layer next to the plate are solved to calculate the flow velocity and the protein mass flux. The local rate of growth of the plate is shown to vary significantly with depth due to the convection. For an aqueous solution of lysozyme at a concentration of 40 mg/ml, the boundary layer at the top of a 1-mm-high crystal has a thickness of 80 microns at 1 g, and 2570 microns at 10 to the -6th g.

Baird, J. K.↗
NASA NTRS1986

The Martian climate and energy balance models with CO2/H2O atmospheres

The analysis begins with a seasonal energy balance model (EBM) for Mars. This is used to compute surface temperature versus x = sin(latitude) and time over the seasonal cycle. The core model also computes the evolving boundaries of the CO2 icecaps, net sublimational/condensation rates, and the resulting seasonal pressure wave. Model results are compared with surface temperature and pressure history data at Viking lander sites, indicating fairly good agreement when meridional heat transport is represented by a thermal diffusion coefficient D approx. 0.015 W/sq. m/K. Condensational wind distributions are also computed. An analytic model of Martian wind circulation is then proposed, as an extension of the EMB, which incorporates vertical wind profiles containing an x-dependent function evaluated by substitution in the equation defining the diffusion coefficient. This leads to a parameterization of D(x) and of the meridional circulation which recovers the high surface winds predicted by dynamic Mars atmosphere models (approx. 10 m/sec). Peak diffusion coefficients, D approx. 0.6 w/sq m/K, are found over strong Hadley zones - some 40 times larger than those of high-latitude baroclinic eddies. When the wind parameterization is used to find streamline patterns over Martian seasons, the resulting picture shows overturning hemispheric Hadley cells crossing the equator during solstices, and attaining peak intensities during the south summer dust storm season, while condensational winds are most important near the polar caps.

Hoffert, M. I.↗