SEARCH · Search NASA

Results for “Numerical diffusion”

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 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Computer simulation of impurity diffusion in silicon, part 1

The elementary classical models for idealized diffusion conditions are described, and the principles are then used in developing more realistic models. The practical models require some type of numerical analysis. The numerical techniques are outlined and details concerning their implementation are given. Some results are presented which were obtained with the computer programs implementing the numerical techniques with implicit and explicit methods. Special problems of impurity-rich interlayers forming between an oxide and silicon are considered. A set of computed curves for sheet resistance, junction depth, and oxide thickness for different diffusion schedules is included.

Gassaway, J. D.↗

NASA NTRS1996

Effect of chemistry and turbulence on NO formation in oxygen-natural gas flames

The effects of chemistry and turbulence on NO formation in oxygen-natural turbulent diffusion flames gas flames have been investigated. The chemistry of nitric oxides has been studied numerically in the counterflow configuration. Systematic calculations with the GRI 2.11 mechanism for combustion of methane and NO chemistry were conducted to provide a base case. It was shown that the 'simple' Zeldovich mechanism accounts for more than 75% of N2 consumption in the flame in a range of strain-rates varying between 10 and 1000 s-l. The main shortcomings of this mechanism are: 1) overestimation (15%) of the NO production rate at low strain-rates because it does not capture the reburn due to the hydrocarbon chemistry, and 2) underestimation (25%) of the NO production rate at high strainrates because it ignores NO production through the prompt mechanism. Reburn through the Zeldovich mechanism alone proves to be significant at low strain-rates. A one-step model based on the Zeldovich mechanism and including reburn has been developed. It shows good agreement with the GRI mechanism at low strain-rates but underestimates significantly N2 consumption (about 50%) at high strain-rates. The role of turbulence has been assessed by using an existing 3-D DNS data base of a diffusion flame in decaying turbulence. Two PDF closure models used in practical industrial codes for turbulent NO formation have been tested. A simpler version of the global one-step chemical scheme for NO compared to that developed in this study was used to test the closure assumptions of the PDF models, because the data base could not provide all the necessary ingredients. Despite this simplification, it was possible to demonstrate that the current PDF models for NO overestimate significantly the NO production rate due to the fact that they neglect the correlations between the fluctuations in oxygen concentration and temperature. A single scalar PDF model for temperature that accounts for such correlations based on laminar flame considerations has been developed and showed excellent agreement with the values given by the DNS.

Samaniego, J. -M.↗

NASA NTRS1987

Ion and electron velocity distributions within flux transfer events

The detailed nature of the thermal and suprathermal ion and electron distributions within magnetic flux transfer events (FTEs) is examined. Examples of both magnetosheath FTEs and magnetospheric FTEs are discussed. The detailed distributions confirm that FTEs contain a mixture of magnetosheath and magnetospheric plasmas. To lowest order, the distributions are consistent with a simple superposition of the two interpenetrating populations, with no strong interactions between them. To first order, some interesting differences appear, especially in the electron distributions, suggesting that considerable pitch angle scattering and some electron energy diffusion are also occurring. These observations should provide a useful test of analytical and numerical studies of interpenetrating plasmas.

Thomsen, M. F.↗

NASA NTRS1993

A study of the compressible flow through a diffusing S-duct

Benchmark aerodynamic data are presented for compressible flow through a representative S-duct configuration. A numerical prediction of the S-duct flow field, obtained from a subsonic parabolized Navier-Stokes algorithm, is also shown. The experimental and numerical results are compared. Measurements of the three-dimensional velocity field, total pressures, and static pressures were obtained at five cross-sectional planes. Aerodynamic data were gathered with calibrated pneumatic probes. Surface static pressure and surface flow visualization data were also acquired. All reported tests were conducted with an inlet centerline Mach number of 0.6. The Reynolds number, based on the inlet centerline velocity and duct inlet diameter, was 2.6 x 10(exp 6). Thin inlet turbulent boundary layers existed. The collected data should be beneficial to aircraft inlet designers and the measurements are suitable for the validation of computational codes. The results show that a region of streamwise flow separation occurred within the duct. Details about the separated flow region, including mechanisms which drive this complicated flow phenomenon, are discussed. Results also indicate that the duct curvature induces strong pressure driven secondary flows. The cross flows evolve into counter-rotating vortices. These vortices convect low momentum fluid of the boundary layer toward the center of the duct, degrading both the uniformity and magnitude of the total pressure profile.

Wellborn, Steven R.↗

NASA NTRS1992

Laboratory scaled simulation of lidar cloud sounding experiments

The results of lidar measurements carried out on laboratory scale models of clouds are presented. Measurements on laboratory scale models are important since one has the knowledge of the relevant parameters of the diffusing medium, such as: scattering and absorption coefficients, phase function, homogeneity, shape, etc. Knowledge of these parameters enables one to use the results to test the reliability of theoretical and numerical investigations. To obtain a laboratory scaled model of a lidar system sounding a cloud, it is necessary to scale down all the geometrical quantities by the same factor to reduce distances of the order of kilometers to the order of meters, keeping the size and the optical depth of the diffusers unchanged. If a time resolution of the order of nanoseconds is necessary for a lidar sounding actual clouds, the corresponding time resolution for the laboratory model should be of the order of picoseconds. It is possible to obtain this resolution by using picosecond laser systems and fast electrooptical detectors like the streak camera. The results of the laboratory measurements showed that the multiple scattering effect strongly depends on the size of the diffusers, as well as on the concentration. The experimental results were compared with the numerical results of a Monte Carlo code. A generally good agreement was obtained.

Zaccanti, G.↗

NASA NTRS1980

Magnetic pumping of particles in the outer Jovian magnetosphere

The mechanism of magnetic pumping consists of two processes, the adiabatic motion of charged particles in a time varying magnetic field and their pitch-angle diffusion. The result is a systematic increase in the energy of charged particles trapped in mirror (and particularly, magnetospheric) magnetic fields. A numerical model of the mechanism is constructed, compared with analytic theory where possible, and, through elementary exercises, is used to predict the consequences of the process for cases that are not tractable by analytical means. For energy dependent pitch angle diffusion rates, characteristic 'two temperature' distributions are produced. Application of the model to the outer Jovian magnetosphere shows that beyond 20 Jupiter radii in the outer magnetosphere, particles may be magnetically pumped to energies of the order of 1 - 2 MeV. Two temperature distribution functions with "break points" at 1 - 4 KeV for electrons and 8 - 35 KeV for ions are predicted.

Borovsky, J. E.↗

NASA NTRS1990

A turbulent two-phase flow model for nebula flows

A new and very efficient turbulent two-phase flow numericaly model is described to analyze the environment of a protoplanetary nebula at a stage prior to the formation of planets. Focus is on settling processes of dust particles in flattened gaseous nebulae. The model employs a perturbation technique to improve the accuracy of the numerical simulations of such flows where small variations of physical quantities occur over large distance ranges. The particles are allowed to be diffused by gas turbulence in addition to settling under gravity. Their diffusion coefficients is related to the gas turbulent viscosity by the non-dimensional Schmidt number. The gas turbulent viscosity is determined by the means of the eddy viscosity hypothesis that assumes the Reynolds stress tensor proportional to the mean strain rate tensor. Zero- and two-equation turbulence models are employed. Modeling assumptions are detailed and discussed. The numerical model is shown to reproduce an existing analytical solution for the settling process of particles in an inviscid nebula. Results of nebula flows are presented taking into account turbulence effects of nebula flows. Diffusion processes are found to control the settling of particles.

Champney, Joelle M.↗

NASA NTRS2017

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles.

Alexander Radkevich↗

NASA NTRS1988

SINDA - Systems Improved Numerical Differencing Analyzer

Capability to model fluid flow added. SINDA is software system developed to solve physical problems governed by diffusion-type equations mathematically modeled by lumped-parameter representations. Used as general thermal analyzer with resistor-and-capacitor-network representations. Adapted to wide range of problems represented by Fourier, Possion, Laplace, or other differential equations. Written in FORTRAN and Assembler.

Damico, Steve J.↗

NASA NTRS1995

An Experimental and Theoretical Study of Radiative Extinction of Diffusion Flames

The objective of this research was to experimentally and theoretically investigate the radiation-induced extinction of gaseous diffusion flames in microgravity. The microgravity conditions were required because radiation-induced extinction is generally not possible in 1-g but is highly likely in microgravity. In 1-g, the flame-generated particulates (e.g. soot) and gaseous combustion products that are responsible for flame radiation, are swept away from the high temperature reaction zone by the buoyancy-induced flow and a steady state is developed. In microgravity, however, the absence of buoyancy-induced flow which transports the fuel and the oxidizer to the combustion zone and removes the hot combustion products from it enhances the flame radiation due to: (1) transient build-up of the combustion products in the flame zone which increases the gas radiation, and (2) longer residence time makes conditions appropriate for substantial amounts of soot to form which is usually responsible for most of the radiative heat loss. Numerical calculations conducted during the course of this work show that even non-radiative flames continue to become "weaker" (diminished burning rate per unit flame area) due to reduced rates of convective and diffusive transport. Thus, it was anticipated that radiative heat loss may eventually extinguish the already "weak" microgravity diffusion flame. While this hypothesis appears convincing and our numerical calculations support it, experiments for a long enough microgravity time could not be conducted during the course of this research to provide an experimental proof. Space shuttle experiments on candle flames show that in an infinite ambient atmosphere, the hemispherical candle flame in microgravity will burn indefinitely. It was hoped that radiative extinction can be experimentally shown by the aerodynamically stabilized gaseous diffusion flames where the fuel supply rate was externally controlled. While substantial progress toward this goal was made during this project, identifying the experimental conditions for which radiative extinction occurs for various fuels requires further study. Details concerning this research which are discussed in published articles are included in the appendices.

Atreya, Arvind↗

NASA NTRS1999

Numerical Simulations of Acoustically Driven, Burning Droplets

This computational study focuses on understanding and quantifying the effects of external acoustical perturbations on droplet combustion. A one-dimensional, axisymmetric representation of the essential diffusion and reaction processes occurring in the vicinity of the droplet stagnation point is used here in order to isolate the effects of the imposed acoustic disturbance. The simulation is performed using a third order accurate, essentially non-oscillatory (ENO) numerical scheme with a full methanol-air reaction mechanism. Consistent with recent microgravity and normal gravity combustion experiments, focus is placed on conditions where the droplet is situated at a velocity antinode in order for the droplet to experience the greatest effects of fluid mechanical straining of flame structures. The effects of imposed sound pressure level and frequency are explored here, and conditions leading to maximum burning rates are identified.

Kim, H.-C.↗

NASA NTRS1959

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↗

NASA NTRS1972

Experimental determination of the turbulence in a liquid rocket combustion chamber

The intensity of turbulence and the Lagrangian correlation coefficient for a liquid rocket combustion chamber were determined experimentally using the tracer gas diffusion method. The results indicate that the turbulent diffusion process can be adequately modeled by the one-dimensional Taylor theory; however, the numerical values show significant disagreement with previously accepted values. The intensity of turbulence is higher by a factor of about two, while the Lagrangian correlation coefficient which was assumed to be unity in the past is much less than unity.

Hara, J.↗

NASA NTRS1975

Thermodynamics and kinetics of reactions in protective coating systems

A study of the aluminization of Ni from packs containing various percentages of unalloyed Al confirmed that the surface aluminum content of specimens aluminized tends to decrease with time and consequently a simple parabolic law for the weight-gain vs. time relationship is not obeyed. The diffusivity-composition relationship in NiAl was examined, and a set of curves is presented. A numerical method for the calculation of coating dissolution rates was developed and applied to NiAl-Ni3Al type of coatings.

Gupta, B.↗

NASA NTRS1979

On the acceleration of the auroral plasma of the Earth due to conic instability

The role of electrostatic conic instability in the diffusion of spilled particles and the acceleration of the auroral plasma of the Earth are numerically investigated by means of quasi-linear equations. Linear increments to the conic instability under magnetospheric conditions are introduced. The quasi-linear equation is averaged along tubes of force taking into account the drift across the magnetic field lines. It is shown that the existence of a conic instability leads to small spills and isotropic distribution of particles, but is also accompanied by significant acceleration of electrons.

Churayev, R. S.↗

NASA NTRS1979

A limit to the X-ray luminosity of nearby normal galaxies

Emission is studied at luminosities lower than those for which individual discrete sources can be studied. It is shown that normal galaxies do not appear to provide the numerous low luminosity X-ray sources which could make up the 2-60 keV diffuse background. Indeed, upper limits suggest luminosities comparable with, or a little less than, that of the galaxy. This is consistent with the fact that the average optical luminosity of the sample galaxies within approximately 20 Mpc is slightly lower than that of the galaxy. An upper limit of approximately 1% of the diffuse background from such sources is derived.

Worrall, D. M.↗

NASA NTRS1979

Effects of turbulence model selection on the prediction of complex aerodynamic flows

Numerical simulations of viscous transonic flow over a circular-arc airfoil and in a diffuser are described. The simulations are made with a new computer program designed to serve as a tool in the development of improved turbulence models for complex flows. The program incorporates zero-, one-, and two-equation eddy viscosity models and includes a variety of subsonic and supersonic boundary conditions. The airfoil flow contains a shock-separated boundary-layer interaction that has resisted previous attempts at simulation. The diffuser flow also contains a shock-boundary-layer interaction, which has not been simulated previously. Calculations using standard turbulence models, developed originally for incompressible unseparated flows, are described. Results indicate that although there are interesting differences in predictions between the various models, none of them predict the flows accurately. Suggestions for improved turbulence models are discussed.

Coakley, T. J.↗

NASA NTRS1983

An EBIC equation for solar cells

When an electron beam of a scanning electron microscope (SEM) impinges on an N-P junction, the generation of electron-hole pairs by impact ionization causes a characteristic short circuit current I(sc) to flow. The I(sc), i.e., EBIC (electron beam induced current) depends strongly on the configuration used to investigate the cell's response. In this paper the case where the plane of the junction is perpendicular to the surface is considered. An EBIC equation amenable to numerical computations is derived as a function of cell thickness, source depth, surface recombination velocity, diffusion length, and distance of the junction to the beam-cell interaction point for a cell with an ohmic contact at its back surface. It is shown that the EBIC equation presented here is more general and easier to use than those previously reported. The effects of source depth, ohmic contact, and diffusion length on the normalized EBIC characteristic are discussed.

Luke, K. L.↗