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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Improvement of the second- and third-moment modeling of turbulence: A study of Reynolds-stress closure model

Four parts of the Reynolds-stress closure modeling are reported: (1) improvement of the k and epsilon equaitons; (2) development of the third-moment transport equation; (3) formulation of the diffusion coefficient of the momentum equation by using the algebraic-stress model of turbulence; and (4) the application of the Reynolds-stress model to a heat exchanger problem. It was demonstrated that the third-moment transport model improved the prediction of the triple-velocity products in the recirculating and reattaching flow regions in comparison with the existing algebraic models for the triple-velocity products. Optimum values for empirical coefficients are obtained for the prediction of the backward-facing step flows. A functional expression is derived for the coefficient of the momentum diffusion by employing the algebraic-stress model. The second-moment closure is applied to a heat transfer problem. The computations for the flow in a corrugated-wall channel show that the second-moment closure improves the prediction of the heat transfer rates by 30% over the k - epsilon model.

Amano, R. S.↗
NASA NTRS1972

Resonant diffusion in strongly turbulent plasmas.

The effect of turbulent fluctuations on plasma particles is considered, and equations are derived which describe the evolution of macroscopic properties such as temperature and flow speed of the turbulent plasma. Initially, a diffusion equation for a single-particle distribution function averaged over an ensemble of plasmas is derived for an unmagnetized plasma. For the resonant diffusion in strongly turbulent plasmas, an ensemble of three-dimensional plasmas is considered with an approximately homogeneous and stationary distribution of random electromagnetic fluctuations. For each realization, the single-particle distribution function satisfies the Vlasov equation.

Birmingham, T. J.↗
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.↗
DOE OSTI2020

Mixtures of octanol and an ionic liquid: Structure and transport

Ionic liquids (ILs) with long alkyl substituents are amphiphilic, which leads to a bicontinuous liquid structure. The strongly interacting anionic and cationic head groups form a long range charge network, with the hydrocarbon tails forming a nonpolar domain. Such nonpolar domains have been shown to dissolve a variety of neutral organic solvents. In mixtures of ILs with solvents the neutral organic molecules residing in the nonpolar domains experience different environments and friction from the charged cations and anions. Thus, the neutral molecules diffuse much faster than predicted by hydrodynamic scaling using the average viscosity of the mixture. Here, we report studies on the structure and transport properties of mixtures of 1-octanol with the IL trihexyltetradecylphosphonium bis(trifluoromethylsulfonyl)imide (P 6,6,6,14 + /NTf 2 – ). The majority of the atom fraction in the P 6,6,6,14 + cation comprises four hydrocarbon substituents. The unique amphiphilic nature of ILs with the P 6,6,6,14 + cation makes 1-octanol fully miscible with the IL at ambient temperatures. X-ray scattering experiments show that the IL structure persists in the mixtures for 1-octanol mole fractions as large as x oct = 0.90. The self-diffusion coefficients of the three molecular species in the mixtures were measured by NMR experiments. The self-diffusion of the P 6,6,6,14 + cation is well described by the Stokes–Einstein equation, while the diffusivity of the NTf 2 – anion is slightly lower than the hydrodynamic prediction. The measured diffusivities of octanol in these mixtures are 1.3–4 times higher than the hydrodynamic predictions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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.↗
DOE OSTI2023

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗
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.↗