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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Numerical simulation of unsteady flows generated by dissociating nitrogen diffusion

An explicit time accurate predictor-corrector scheme, similar to MacCormack's, has been used to simulate unsteady nonequilibrium laminar chemically-reacting diffusion-reaction generated flows. A zonal approach, along with nonstructured boundaries, allows computation of internal and external flow with the same computer code. Results for equilibrium flows generated by dissociating, vibrationally relaxing nitrogen in a rectangular chamber are presented. The unsteady flows were generated strictly by the chemical reactions and concentration gradients.

Marraffa, Lionel↗

NASA NTRS2004

A Numerical Investigation of the Extinction of Low Strain Rate Diffusion Flames by an Agent in Microgravity

Our goal has been to investigate the influence of both dilution and radiation on the extinction process of nonpremixed flames at low strain rates. Simulations have been performed by using a counterflow code and three radiation models have been included in it, namely, the optically thin, the narrowband, and discrete ordinate models. The counterflow flame code OPPDIFF was modified to account for heat transfer losses by radiation from the hot gases. The discrete ordinate method (DOM) approximation was first suggested by Chandrasekhar for solving problems in interstellar atmospheres. Carlson and Lathrop developed the method for solving multi-dimensional problem in neutron transport. Only recently has the method received attention in the field of heat transfer. Due to the applicability of the discrete ordinate method for thermal radiation problems involving flames, the narrowband code RADCAL was modified to calculate the radiative properties of the gases. A non-premixed counterflow flame was simulated with the discrete ordinate method for radiative emissions. In comparison with two other models, it was found that the heat losses were comparable with the optically thin and simple narrowband model. The optically thin model had the highest heat losses followed by the DOM model and the narrow-band model.

Puri, Ishwar K.↗

NASA NTRS1988

Does the resistive tearing instability nonlinearly evolve to a fast reconnection mode

A fundamental problem in applying linear tearing instability theory to the rapid processes (particle acceleration, heating) in flares was the characteristically slow rate of reconnection. This problem can be at least partially overcome if the tearin mode nonlinearly evolves to a regime in which the reconnection rate is substantially enhanced, such as that for the Petschek configuration. This possibility was often suggested, and some numerical simulations appear to provide support for such a view. Numerical simulation are used to study the nonlinear evolution of the tearing stability and show that a fast Petschek-like regime may not be achieved. This conclusion follows when there are sufficient grid points within the diffusion region to completely resolve the nonlinear dynamic interactions in the diffusion layer. When the numerical resolution is not adequate, the solution does appear to approach a Petschek configuration. The resolved solution contains reverse flow vortices and current sheets, terminated with a current reversal, similar to those obtained by Syrovatsky (JEPT, 33, 933, 1971).

Steinolfson, R. S.↗

NASA NTRS1993

New developments in the method of space-time conservation element and solution element: Applications to the Euler and Navier-Stokes equations

A new numerical framework for solving conservation laws is being developed. This new approach differs substantially in both concept and methodology from the well-established methods--i.e., finite difference, finite volume, finite element, and spectral methods. It is conceptually simple and designed to avoid several key limitations to the above traditional methods. An explicit model scheme for solving a simple 1-D unsteady convection-diffusion equation is constructed and used to illuminate major differences between the current method and those mentioned above. Unexpectedly, its amplification factors for the pure convection and pure diffusion cases are identical to those of the Leapfrog and the DuFort-Frankel schemes, respectively. Also, this explicit scheme and its Navier-Stokes extension have the unusual property that their stabilities are limited only by the CFL condition. Moreover, despite the fact that it does not use any flux-limiter or slope-limiter, the Navier-Stokes solver is capable of generating highly accurate shock tube solutions with shock discontinuities being resolved within one mesh interval. An accurate Euler solver also is constructed through another extension. It has many unusual properties, e.g., numerical diffusion at all mesh points can be controlled by a set of local parameters.

Chang, Sin-Chung↗

NASA NTRS2001

Modelling of Trapped Radiation Near Jupiter

Energetic (62 to approx. 130 MeV) proton fluxes measured with the Galileo Probe inside Jupiter's main ring (radii 122,500 to 128,940 km) show a modest increase as Jupiter is approached. Solutions of a reduced (1-d) equatorial diffusion equation with-constant losses described by a lifetime tau match those data for tau greater than or equal to 10(exp 9) s for diffusion coefficient D(sub LL) = 10(exp -9) L(exp 4)/s, 10(exp -9) L(exp -3)/s, or 10(exp -10) L(exp 4)/s (this particle population may be undergoing nearly loss-free inward radial diffusion). Tau would increase as Jupiter is approached if its value were determined by pitch angle scattering due to EM wave-particle interactions. if the absolute amplitudes of the waves' magnetic fluctuations did not vary with radial distance. Exploration of numerical solutions of the diffusive transport equation for Jupiter's magnetospheric ions at die planetary ring and closer to Jupiter has been done. Explicit range-energy relationships for energy loss in ring matter (SiO2), modeled as a continuous disk, are incorporated. A spatial resolution of approx. 1000 km in the radial direction is used in the vicinity of the main ring; off-equatorial ion fluxes have been included (assuming the atmosphere is a perfect absorber). For protons the maximum effect of energy loss in the microscopic size range of ring particles is expected at energies approx. 0.1 Mev. Presumably there could be sufficiently large, unobserved bodies within the ring that would have to be modeled as discrete objects, as planetary satellites would be. Within limitations of the model and computational resources it's planned to characterize solutions for transport of magnetospheric ions past the planetary ring; the procedures that are used are fairly well known and may be applied generally to the simpler magnetospheric models.

Mihalov, John D.↗

NASA NTRS1985

Diffusion-convection function of cosmic rays

The fundamental properties and some numerical results of the solution of the diffusion equation of an impulsive cosmic-ray point source in an uniform, unbounded and spherically symmetrical moving medium is presented. The diffusion-convection(D-C) function is an elementary composite function of the solution of the D-C equation for the particles injected impulsively from a diffusive point source into the medium. It is the analytic solution derived by the dimensional method for the propagation equation of solar cosmic rays in the heliosphere, i.e. the interplanetary space. Because of the introduction of convection effect of solar wind, a nonhomogeneous term appears in the propagation equation, it is difficult to express its solution in terms of the ordinary special functions. The research made so far has led to a solution containing only the first order approximation of the convection effect.

Zhang, G.↗

NASA NTRS1998

Volume Diffusion Growth Kinetics and Step Geometry in Crystal Growth

The role of step geometry in two-dimensional stationary volume diff4sion process used in crystal growth kinetics models is investigated. Three different interface shapes: a) a planar interface, b) an equidistant hemispherical bumps train tAx interface, and c) a train of right angled steps, are used in this comparative study. The ratio of the super-saturation to the diffusive flux at the step position is used as a control parameter. The value of this parameter can vary as much as 50% for different geometries. An approximate analytical formula is derived for the right angled steps geometry. In addition to the kinetic models, this formula can be utilized in macrostep growth models. Finally, numerical modeling of the diffusive and convective transport for equidistant steps is conducted. In particular, the role of fluid flow resulting from the advancement of steps and its contribution to the transport of species to the steps is investigated.

Mazuruk, Konstantin↗

NASA NTRS1998

Temperature Dependence of Diffusivities in Liquid Elements (LMD)

This research was to advance the understanding of diffusion mechanisms in liquid metals and alloys through accurate diffusivity measurements over a wide range of temperatures, including the proximity of the materials melting points. Specifically, it was driven towards developing a methodology (and subsequent flight hardware) to enable several diffusion coefficient measurements (i.e., at several different temperatures) to be performed using a single sample. The Liquid Metal Diffusion (LMD) was funded as a Flight Definition Project in February 1993 in response to NRA 91-OSSA-20 (Microgravity Science and Applications Division). The Science Concept Review for LAID was held during April 1994. In January 1995 we were informed that we had failed this review and the project was change to ground-based activities only. A new proposal was submitted for the next NRA addressing the panels concerns. As part of NASA's Risk Mitigation program, a scaled-down version of the hardware was funded in July of 1995 for a flight opportunity utilizing experiment on the Microgravity Isolation Mount. This experiment was to determine the self-diffusivity of indium at 185 C. The LMD was transferred to the Mir Space Station in STS-81 and returned on STS-84 (January - May 1997). Three, out of five, self-diffusion data sets were returned. A description of this experiment/hardware is included below. This summary is only intended to give the reader an overview of the results obtained for the tasks outlined in the original proposal. Research that was not published is explained in more detail. At the end of this report is a list of refereed publications and invited talks that were given as a result of this work. The reader is directed to these for further details. Attachment: Real-time diffusivity measurements in liquids at several temperatures with one sample, On the insensitivity of liquid diffusivity measurements to deviations from 1D transport, and Numerical simulations of the convective contamination of diffusivity measurements in liquids.

Banish, R. Michael↗

NASA NTRS1995

Putting quasilinear theory to the test: A comparison of analytical theory with results from a trajectory code

Parallel and perpendicular diffusion coefficients were computed numerically by following particle orbits in a simulated magnetic field. The simulated field was chosen to have delta B/B(sub o) small, so as to provide a test of quasilinear theory in a regime where the theory should be most accurate. The simulation space is large enough to contain many magnetic field correlation lengths, so that effects of field line random walk can be studied. After presenting results for parallel diffusion, we will focus on two controversial issues relating to perpendicular diffusion: (1) Do quasilinear descriptions of perpendicular diffusion retain any validity for particles whose Larmor radius is smaller than a correlation length? (2) Does field line random walk lead to particle diffusion in the usual sense, or does it produce 'compound' diffusion for which particles spread out proportionally to t(exp 1/4) instead of t(exp 1/2)?

Bieber, J. W.↗

NASA NTRS1992

ULTRA-SHARP solution of the Smith-Hutton problem

Highly convective scalar transport involving near-discontinuities and strong streamline curvature was addressed in a paper by Smith and Hutton in 1982, comparing several different convection schemes applied to a specially devised test problem. First order methods showed significant artificial diffusion, whereas higher order methods gave less smearing but had a tendency to overshoot and oscillate. Perhaps because unphysical oscillations are more obvious than unphysical smearing, the intervening period has seen a rise in popularity of low order artificially diffusive schemes, especially in the numerical heat transfer industry. The present paper describes an alternate strategy of using non-artificially diffusive high order methods, while maintaining strictly monotonic transitions through the use of simple flux limited constraints. Limited third order upwinding is usually found to be the most cost effective basic convection scheme. Tighter resolution of discontinuities can be obtained at little additional cost by using automatic adaptive stencil expansion to higher order in local regions, as needed.

Leonard, B. P.↗

NASA NTRS1975

Simulation of the velocity diffusion of charged particles in turbulent magnetic fields

Numerical simulation experiments have been performed in order to study the velocity diffusion of charged particles in a static turbulent magnetic field. By following orbits of particles moving in a large ensemble of random magnetic field realizations with suitably-chosen statistical properties, a pitch-angle diffusion coefficient is derived. Results are presented for a variety of particle rigidities and r.m.s. random field strengths and a power spectrum typical of conditions found in interplanetary space.

Kaiser, T. B.↗

NASA NTRS1986

Unsteady separated boundary layer in a transonic diffuser flow with self-excited oscillations

A numerical investigation of two-dimensional unsteady boundary layer in a transonic diffuser flow with self-excited oscillations and strong flow separation by solving the compressible, Reynolds-averaged, thin-layer Navier-Stokes equations with two-equations turbulence model is described. Three different meshes with constant streamwise mesh distribution and varying vertical mesh distribution were used. Results obtained indicate that a refinement of mesh studied here has minimal effect on the mean boundary layer flow but significantly increases the amplitude of oscillation of all flow variables. Comparisons of unsteady wall pressure, velocity profile, terminal shock, and separation pocket among computations and with experiment are presented.

Hsieh, T.↗

NASA NTRSDate not supplied

NASA HECC Geometry and Performance Review Part 3: A Numerical and Experimental Investigation of Tip Clearance Effects on the Vaneless Diffuser Configuration

Tip clearance effects in centrifugal compressors have been extensively investigated to understand the losses associated with the flow in the impeller tip clearance gap between the rotating blades and the stationary shroud. In Part 3 of this multipart investigation, experimental data and validated numerical simulations from Parts 1 and 2 were used to analyze the effects of the size of the tip clearance gap on the High Efficiency Centrifugal Compressor performance and aerodynamics. Four tip clearance gaps ranging from 0.012-in to 0.030-in (2.0% to 4.9% of the impeller exit blade height) were considered at both design and off-design operating conditions. The total pressure ratio and efficiency of the stage are found to decay linearly with increasing tip gaps. The sensitivity of the impeller performance to the tip gap was found to vary with the rotational speed of the compressor. Spanwise surveys of flow angle, total pressure, and total temperature collected at the impeller exit at design speed are used to validate numerical simulations at each experimental tip gap condition. Numerical simulations show that increased turbulence kinetic energy near the shroud at larger tip gaps leads to spanwise mixing of high entropy fluid near the impeller trailing edge which decreases the useful work input by the impeller. The data presented have been made available to the public in the HECC Data Archive located at https://storage.googleapis.com/hecc-data/NASA-HECC-Data-Archive.zip.

centrifugal compressor↗

NASA NTRSDate not supplied

NASA HECC Geometry and Performance Review Part 3: A Numerical and Experimental Investigation of Tip Clearance Effects on the Vaneless Diffuser Configuration

centrifugal compressor↗

NASA NTRS1983

Overlay coating degradation by simultaneous oxidation and coating/substrate interdiffusion

Degradation of NiCrAlZr overlay coatings on various NiCrAl substrates was examined after cyclic oxidation. Concentration/distance profiles were measured in the coating and substrate after various oxidation exposures at 1150 C. For each stubstrate, the Al content in the coating decreased rapidly. The concentration/distance profiles, and particularly that for Al, reflected the oxide spalling resistance of each coated substrate. A numerical model was developed to simulate diffusion associated with overlay-coating degradation by oxidation and coating/substrate interdiffusion. Input to the numerical model consisted of the Cr and Al content of the coating and substrate, ternary diffusivities, and various oxide spalling parameters. The model predicts the Cr and Al concentrations in the coating and substrate after any number of oxidation/thermal cycles. The numerical model also predicts coating failure based on the ability of the coating to supply sufficient Al to the oxide scale. The validity of the model was confirmed by comparison of the predicted and measured concentration/distance profiles. The model was subsequently used to identify the most critical system parameters affecting coating life.

Nesbitt, J. A.↗

NASA NTRS1976

Test particle propagation in magnetostatic turbulence. 1. Failure of the diffusion approximation

The equation which governs the quasi-linear approximation to the ensemble and gyro-phase averaged one-body probability distribution function is constructed from first principles. This derived equation is subjected to a thorough investigation in order to calculate the possible limitations of the quasi-linear approximation. It is shown that the reduction of this equation to a standard diffusion equation in the Markovian limit can be accomplished through the application of the adiabatic approximation. A numerical solution of the standard diffusion equation in the Markovian limit is obtained for the narrow parallel beam injection. Comparison of the diabatic and adiabatic results explicitly demonstrates the failure of the Markovian description of the probability distribution function. Through the use of a linear time-scale extension the failure of the adiabatic approximation, which leads to the Markovian limit, is shown to be due to mixing of the relaxation and interaction time scales in the presence of the strong mean field.

Klimas, A. J.↗

NASA NTRS1995

Diffusion at the magnetopause: Hybrid simulations

Electromagnetic wave generation and resulting cross-field diffusion of plasma are considered at a tangential discontinuity, which characterizes the magnetopause for northward interplanetary magnetic field. Two-dimensional hybrid (particle ions, massless fluid electrons) simulations, in which the tangential discontinuity is generated self-consistently via a stream-stream interaction, are used to show that wave growth occurs when the ambient magnetic field is predominantly perpendicular to the direction of the density gradient. Low-frequency (much less than ion gyrofrequency) waves, with amplitudes delta B/B less than or equal to 0.2 and anticorrelated density fluctuations delta n/n less than or equal to 0.6, are generated at the discontinuity, resulting in cross-field diffusion that is comparable to the Bohm rate. Both the fluctuation level and the lack of ion heating in the calculations are consistent with observations at the magnetopause. The magnitude of the diffusion is considered in the presence of numerical effects and in the context of the inferred diffusion rate at the magnetopause. The relation of the low-frequency waves and their consequences to faster growing, short-wavelength waves due to the lower hybrid drift instability is also addressed. The overall conclusion of this initial study is that diffusion due to low frequency waves is not likely to be a major effect at the magnetopause.

Winske, D.↗

NASA NTRS1992

Approximate Solutions Of Equations Of Steady Diffusion

Rigorous analysis yields reliable criteria for "best-fit" functions. Improved "curve-fitting" method yields approximate solutions to differential equations of steady-state diffusion. Method applies to problems in which rates of diffusion depend linearly or nonlinearly on concentrations of diffusants, approximate solutions analytic or numerical, and boundary conditions of Dirichlet type, of Neumann type, or mixture of both types. Applied to equations for diffusion of charge carriers in semiconductors in which mobilities and lifetimes of charge carriers depend on concentrations.

Edmonds, Larry D.↗