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 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Aerothermal modeling program, phase 2

The main objective of the NASA sponsored Aerothermal Modeling Program, Phase 2--Element A, is to develop an improved numerical scheme for predicting combustor flow fields. This effort consists of the following three technical tasks. Task 1 involves the selection and evaluation of various candidate numerical techniques. Task 2 involves an in-depth evaluation of the selected numerical schemes. Task 3 involves the convection-diffusion scheme and the direct solver that will be incorporated in the NASA 3-D elliptic code (COM3S).

Karki, K. C.↗

NASA NTRS1991

A study of hydrogen diffusion flames using PDF turbulence model

The application of probability density function (pdf) turbulence models is addressed. For the purpose of accurate prediction of turbulent combustion, an algorithm that combines a conventional computational fluid dynamic (CFD) flow solver with the Monte Carlo simulation of the pdf evolution equation was developed. The algorithm was validated using experimental data for a heated turbulent plane jet. The study of H2-F2 diffusion flames was carried out using this algorithm. Numerical results compared favorably with experimental data. The computations show that the flame center shifts as the equivalence ratio changes, and that for the same equivalence ratio, similarity solutions for flames exist.

Hsu, Andrew T.↗

NASA NTRS1991

Modelling of the dynamics of HgCdTe growth by the vertical Bridgman method

The analysis of Kim and Brown (1989), which described the pseudo-steady state of the HgCdTe melt/crystal system grown by the vertical Bridgman method, is extended to the study of the transients during the crystal growth and, in particular, at the beginning of the growth. The time-dependent numerical simulations developed show the development of a diffusion layer in the melt and the direct relationship between this layer and the intensity and the pattern of convection in the melt. It is shown that the spread of the solute layer occurs on a diffusive time scale, with the initial transient to the steady-state growth rate lasting several hours for the slow growth rates that are typical of the HgCdTe growth. The simulations demonstrate all the features observed experimentally by Szofran et al. (1984).

Kim, Do H.↗

NASA NTRS1992

A pressure correction method for the calculation of compressible chemical reacting flows

A recently developed noniterative method for the solution of the transient fluid flow equations at all speed is extended to handle chemical reacting flows. The species conservation equations are loosely coupled into the predictor/multicorrector sequence of the solution procedure. A split-operator method separates the chemical kinetics terms from the fluid-dynamical terms, as well as an implicit differencing method enhance the numerical stability. The method was applied for turbulent diffusion flame calculations and for the analyses of high pressure, axisymmetric turbulent hypersonic nozzle flows. The diffusion flame results were compared with a similar pressure method for fast chemistry integration scheme without operator-splitting. Simulations of the nozzle flow indicated that the nonideal intermolecular effects must be included in the analysis and design of high pressure hypersonic nozzle.

Chen, Z. J.↗

NASA NTRS1991

Mechanisms of voids formation during cooldown and freezing of lithium in SP-100 type systems

The mechanisms of void formation during the cooldown and freezing of lithium coolant within the primary loop of SP-100 type systems are investigated. These mechanisms are: (a) homogeneous nucleation, (b) heterogeneous nucleation, (c) normal segregation of helium gas dissolved in liquid lithium, and (d) shrinkage of lithium during freezing. To evaluate the void formation potential due to segregation, a numerical scheme that couples the freezing and mass diffusion processes in both the solid and liquid regions is developed. The results indicated that the formation of He bubbles is unlikely by either homogeneous or heterogeneous nucleation during the cooldown process. However, homogeneous nucleation of He bubbles following the segregation of dissolved He in liquid Lithium ahead of the solid-liquid interface is likely to occur. Results also show that total volume of He void is insignificant when compared to that of shrinkage voids.

Yang, Jae Y.↗

NASA NTRS2002

Discontinuous Galerkin Methods for Turbulence Simulation

A discontinuous Galerkin (DG) method is formulated, implemented, and tested for simulation of compressible turbulent flows. The method is applied to turbulent channel flow at low Reynolds number, where it is found to successfully predict low-order statistics with fewer degrees of freedom than traditional numerical methods. This reduction is achieved by utilizing local hp-refinement such that the computational grid is refined simultaneously in all three spatial coordinates with decreasing distance from the wall. Another advantage of DG is that Dirichlet boundary conditions can be enforced weakly through integrals of the numerical fluxes. Both for a model advection-diffusion problem and for turbulent channel flow, weak enforcement of wall boundaries is found to improve results at low resolution. Such weak boundary conditions may play a pivotal role in wall modeling for large-eddy simulation.

Collis, S. Scott↗

NASA NTRS1998

On the Origin of Chaos in the Asteroid Belt

We consider the effect of gravitational perturbations from Jupiter on the dynamics of asteroids, when Jupiter is itself perturbed by Saturn. The presence of Saturn introduces a number of additional frequencies into Jupiters orbit. These frequencies in turn produce chaos in narrow regions on either side of the chaotic zones associated with the mean motion resonances between the asteroids and Jupiter. The resonant arguments of these three-body resonances contain the longitudes of Jupiter and the asteroid together with either the secular frequency 9-6, or the longitude of Saturn. Resonances involving the longitude of Saturn are analogs of the Laplace resonance in the Jovian satellite system. We show that many three-body resonances involving the longitude of Saturn are chaotic. We give simple expressions for the width of the chaotic region and the associated Lyapunov time. In some cases the chaos can produce a diffusive growth in the 4 eccentricity of the asteroid that leads to ejection of the asteroid on times shorter than the age of the solar system. We give simple estimates for the diffusion time. Finally, we present the results of numerical integrations testing the theory.

Murray, N.↗

NASA NTRS2004

Flow field topology of transient mixing driven by buoyancy

Transient mixing driven by buoyancy occurs through the birth of a symmetric Rayleigh-Taylor morphology (RTM) structure for large length scales. Beyond its critical bifurcation the RTM structure exhibits self-similarity and occurs on smaller and smaller length scales. The dynamics of the RTM structure, its nonlinear growth and internal collision, show that its genesis occurs from an explosive bifurcation which leads to the overlap of resonance regions in phase space. This event shows the coexistence of regular and chaotic regions in phase space which is corroborated with the existence of horseshoe maps. A measure of local chaos given by the topological entropy indicates that as the system evolves there is growth of uncertainty. Breakdown of the dissipative RTM structure occurs during the transition from explosive to catastrophic bifurcation; this event gives rise to annihilation of the separatrices which drives overlap of resonance regions. The global bifurcation of explosive and catastrophic events in phase space for the large length scale of the RTM structure serves as a template for which mixing occurs on smaller and smaller length scales. Copyright 2004 American Institute of Physics.

Nonlinear Dynamics↗

NASA NTRS2015

Monolithic Solid-State Lasers for Spaceflight

A new solution for building high power, solid state lasers for space flight is to fabricate the whole laser resonator in a single (monolithic) structure or alternatively to build a contiguous diffusion bonded or welded structure. Monolithic lasers provide numerous advantages for space flight solid state lasers by minimizing misalignment concerns. The closed cavity is immune to contamination. The number of components is minimized thus increasing reliability. Bragg mirrors serve as the high reflector and output coupler thus minimizing optical coatings and coating damage. The Bragg mirrors also provide spectral and spatial mode selection for high fidelity. The monolithic structure allows short cavities resulting in short pulses. Passive saturable absorber Q-switches provide soft aerturing fro spatial mode filtering and improved pointing stability. We will review our recent commercial and in-house developments toward fully monolithic solid state lasers.

Monolthic laser↗

NASA NTRS1977

Cosmic diffuse gamma rays from 2 to 25 MeV

Results are reported for measurements of cosmic diffuse gamma rays performed with a double Compton scatter gamma-ray telescope flows aboard a balloon at 3.5-g/sq cm residual atmosphere. The operation and calibration of the telescope are briefly reviewed, unwanted backgrounds are identified, and the only significant background for the telescope used is shown to arise from the interaction of albedo neutrons with hydrogen and carbon in one of the liquid scintillators. The fluxes of cosmic diffuse gamma rays with energies of 2 to 3, 3 to 5, 5 to 7, 7.5 to 10 MeV are listed, and upper limits are provided for energies of 10 to 15 and 15 to 25 MeV. These data are compared with numerous previous measurements, and no indication is found that the diffuse radiation is statistically different from isotropic. Upper limits are placed on the gamma-ray fluxes (mainly 3 to 25 MeV) from various localized sources.

White, R. S.↗

NASA NTRS1990

A new class of solutions for interstellar magnetohydrodynamic shock waves

An analysis is presented of the equations of motion for steady MHD shock waves proopagating in interstellar clouds, for boundary conditions that preclude C shocks. In addition to J shocks, in which the neutral fluid component becomes subsonic at an adiabatic jump front, the equations admit a new class of solutions, called C-asterisk shocks, in which the transition to subsonic flow occurs continuously at a sonic point. Numerical methods are developed for computing the structure of J and C-asterisk shocks propagating in diffuse interstellar clouds. The effects of chemical, ionization, and recombination processes are included in this treatment. An alternative numerical method, which uses artificial viscosity to facilitate integration through sonic points, is analyzed and shown to be invalid. A set of exemplary solutions, computed for realistic shock parameters, shows that C-asterisk shocks occur for a broad range of conditions relevant to diffuse interstellar clouds.

Roberge, W. G.↗

NASA NTRS1988

Diffusional decay of striations

This paper investigates the amount of decay of growth striations in bulk crystals grown from melts and solutions, due to solid-state diffusion during the cooling of the crystal to room temperature. A model is developed based on the assumption of a sinusoidal variation in the solid composition at the interface of a semiinfinite ingot growing at constant velocity, and the damping behavior of this variation is followed as it moves away from the solid-liquid interface. Numerical results are presented showing that an increase in the diffusivity causes the striations to penetrate deeper into the solid and, consequently causes the distance for 99.9 percent decay to increase. Analytical results are also presented for the decay of striations in an isothermal ingot induced by a sinusoidal fluctuation in the rate of crystal growth.

Gray, R. T.↗

NASA NTRS1978

Galactic civilizations: Population dynamics and interstellar diffusion

The interstellar diffusion of galactic civilizations is reexamined by potential theory; both numerical and analytical solutions are derived for the nonlinear partial differential equations which specify a range of relevant models, drawn from blast wave physics, soil science, and, especially, population biology. An essential feature of these models is that, for all civilizations, population growth must be limited by the carrying capacity of the environment. Dispersal is fundamentally a diffusion process; a density-dependent diffusivity describes interstellar emigration. Two models are considered: the first describing zero population growth (ZPG), and the second which also includes local growth and saturation of a planetary population, and for which an asymptotic traveling wave solution is found.

Newman, W. I.↗

NASA NTRS1993

Recent developments in high order K-exact reconstruction on unstructured meshes

This paper presents recent improvements in high-order K-exact reconstruction on unstructured meshes. The new reconstruction procedures are incorporated into a basic upwind finite-volume scheme suitable for solving scalar advection-diffusion equations as well as the Euler and Navier-Stokes equations. Numerical calculations are performed comparing the present method with lower order accurate reconstruction procedures (piecewise constant and piecewise linear) and various competing technologies such as the fluctuation splitting method of Roe (1987) and Deconinck et al. (1992) and a system-variant of the streamline diffusion Petrov-Galerkin method developed by Hansbo (1991) and Hansbo and Johnson (1991). Five test problems are used in the numerical comparisons: scalar circular advection, transonic and supersonic Euler flow, laminar boundary-layer flow, and general compressible Navier-Stokes flow.

Barth, Timothy J.↗

NASA NTRS2003

Transient Evolution of a Planar Diffusion Flame Aft of a Translating Flat Plate

The high degree of spatial symmetry of a planar diffusion flame affords great simplifications for experimental and modeling studies of gaseous fuel combustion. Particularly, in a microgravity environment, where buoyancy effects are negligible, an effectively strain-rate-free, vigorous flame may be obtained. Such a flame can also provide long residence times and large length scales for practical probing of flame structures and soot processes. This 2-D numerical study explores the feasibility of establishing such a planar diffusion flame in an enclosed container utilizing a realistic test protocol for a microgravity experiment. Fuel and oxygen mixtures, initially segregated into two half-volumes of a squat rectangular container by a thin separator, are ignited as soon as a flammable mixture is formed in the wake of the separator withdrawn in the centerplane. A triple-flame ensues that propagates behind the trailing edge of the separator. The results of calculations show that the mechanically- and thermally-induced convection decays in about two seconds. The establishment of a planar diffusion flame after this period seems feasible in the central region of the container with sufficient quantities of reactants left over for subsequent studies. An analysis of the flame initiation and formation process suggests how the feasibility of creating such a flame can be further improved.

Gokoglu, Suleyman A.↗

NASA NTRS2003

In General, the Total Voltammetric Current from a Mixture of Redox-Active Substances will Not be the Sum of the Currents that Each Substance would Produce Independently at the Same Concentration as in the Mixture

At the potential range where both decamethylferrocene (dMeFc) and ferrocene (Fc) are oxidized with rates controlled by linear diffusion, electrogenerated Fc(+) radicals diffusing outwards from the electrode react quantitatively (K23 C=5.8 x 10(exp 8) with dMeFc diffusing towards the electrode and produce Fc and dMeFc. That reaction replaces dMeFc with Fc, whose diffusion coefficient is higher than that of dMeFc(+), and the total mass-transfer limited current from the mixture is increased by approximately 10%. Analogous observations are made when mass-transfer is controlled by convective-diffusion as in RDE voltammetry. Similar results have been obtained with another, and for all practical purposes randomly selected pair of redox-active substances, [Co(bipy)3](2+) and N - methylphenothiazine (MePTZ); reaction of MePTZ(+) with [Co(bipy)3](2+) replaces the latter with MePTZ, which diffuses faster and the current increases by approximately 20%. The experimental voltammograms have been simulated numerically and the role of (a) the rate constant of the homogeneous reaction; (b) the relative concentrations; and, (c) the diffusion coefficients of all species involved have been studied in detail. Importantly, it was also identified that within any given redox system the dependence of the mass-transfer limited current on the bulk concentrations of the redox-active species is expected to be non-linear. These findings are discussed in terms of their electroanalytical implications.

Leventis, Nicholas↗

NASA NTRS1972

Effect of eddy diffusivity on wind-driven currents in a two-layer stratified lake

The steady state wind-driven circulation was numerically calculated in a rectangular stratified lake. The lake is composed of two layers having uniform but unequal densities and eddy diffusivities. The position in thermocline and the three-dimensional velocities in both layers calculated using shallow lake equations. The results show that, as the eddy diffusivity in the hypolimnion is increased, the thermocline tilt and hypolimnetic velocities increase. The effect of the other variables such as wind stress, density, basin length, and mean thermocline depth are also shown.

Gedney, R. T.↗

NASA NTRS1990

Axions and SN 1987A: Axion trapping

If an axion of mass between about 10(exp -3) and 10 eV exists, axion emission would have significantly affected the cooling of the nascent neutron star associated with SN 1987A. For an axion of mass greater than about 10(exp -2) eV axions would, like neutrinos, have a mean-free path that is smaller than the size of a neutron star, and thus would become trapped and radiated from an axion sphere. The trapping regime is treated by using numerical models of the initial cooling of a hot neutron star that incorporate a diffusion approximation for axion-energy transport. The axion opacity due to inverse nucleon-nucleon, axion bremsstrahlung is computed; and then the numerical models are used to calculate the integrated axion luminosity, the temperature of the axion sphere, and the effect of axion emission on the neutrino bursts detected by the Kamiokande II (KII) and Irvine-Michigan-Brookhaven (IMB) water-Cherenkov detectors. The larger the axion mass, the stronger the trapping and the smaller the axion luminosity. The estimate of the axion mass is confirmed above which trapping is so strong that axion emission does not significantly affect the neutrino burst. Based upon the neutrino-burst duration - the most sensitive barometer of axion cooling - it is concluded that for an axion mass greater than about 3 eV axion emission would not have had a significant effect on the neutrino bursts detected by KII and IMB. It is strongly suggested that an axion with mass in the interval 10(exp -3) to 3 eV is excluded by the observation of neutrinos from SN 1987A.

Burrows, Adam↗