Search NASA⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

Conditional statistics in a turbulent premixed flame derived from direct numerical simulation

The objective of this paper is to briefly introduce conditional moment closure (CMC) methods for premixed systems and to derive the transport equation for the conditional species mass fraction conditioned on the progress variable based on the enthalpy. Our statistical analysis will be based on the 3-D DNS database of Trouve and Poinsot available at the Center for Turbulence Research. The initial conditions and characteristics (turbulence, thermo-diffusive properties) as well as the numerical method utilized in the DNS of Trouve and Poinsot are presented, and some details concerning our statistical analysis are also given. From the analysis of DNS results, the effects of the position in the flame brush, of the Damkoehler and Lewis numbers on the conditional mean scalar dissipation, and conditional mean velocity are presented and discussed. Information concerning unconditional turbulent fluxes are also presented. The anomaly found in previous studies of counter-gradient diffusion for the turbulent flux of the progress variable is investigated.

Mantel, Thierry↗

Energy densities and temperature associated with the diffuse X-ray background, galaxies and large scale structures

Model independent evidence indicating numerical coincidences between energy densities and temperatures which indicate possible connections between the diffuse X-ray background and the formation of galaxies and large scale structures is presented. A consideration of characteristic temperatures associated with galaxies, peculiar velocities of galaxies, and voids suggests characteristic energies of emission. Barring complications induced by magnetic fields, these temperatures may also be associated with components of the X-ray background. Thermal bremsstrahlung emmission associated with the formation of large scale structures may be the dominant contributor to the high energy X-ray background, 5 to 200 keV. Thermal bremsstrahlung emission produced by galaxies at the epoch of their formation may contribute significantly to the low energy X-ray background, 0.1 to 3 keV. It is concluded that an important, if not dominant fraction of the X-ray background, was probably emitted by hot gas associated with the formation of large scale structures and galaxies.

Daly, Ruth A.↗

Numerical methods of solving a system of multi-dimensional nonlinear equations of the diffusion type

The principles of conservation and stability of difference schemes achieved using the iteration control method were examined. For the schemes obtained of the predictor-corrector type, the conversion was proved for the control sequences of approximate solutions to the precise solutions in the Sobolev metrics. Algorithms were developed for reducing the differential problem to integral relationships, whose solution methods are known, were designed. The algorithms for the problem solution are classified depending on the non-linearity of the diffusion coefficients, and practical recommendations for their effective use are given.

Agapov, A. V.↗

Analysis of turbulent free-jet hydrogen-air diffusion flames with finite chemical reaction rates

A numerical analysis is presented of the nonequilibrium flow field resulting from the turbulent mixing and combustion of an axisymmetric hydrogen jet in a supersonic parallel ambient air stream. The effective turbulent transport properties are determined by means of a two-equation model of turbulence. The finite-rate chemistry model considers eight elementary reactions among six chemical species: H, O, H2O, OH, O2 and H2. The governing set of nonlinear partial differential equations was solved by using an implicit finite-difference procedure. Radial distributions were obtained at two downstream locations for some important variables affecting the flow development, such as the turbulent kinetic energy and its dissipation rate. The results show that these variables attain their peak values on the axis of symmetry. The computed distribution of velocity, temperature, and mass fractions of the chemical species gives a complete description of the flow field. The numerical predictions were compared with two sets of experimental data. Good qualitative agreement was obtained.

Sislian, J. P.↗

Erosional dynamics, flexural isostasy, and long-lived escarpments: A numerical modeling study

Erosional escarpments common features of high-elevation rifted continets. Fission track data suffest that these escarpments form by base level lowering and/or marginal uplift during rifting, followed by lateral retreat of an erosion front across tens to hundreds of kioometers. Previous modeling studies have shown that this characteristic pattern of denudation can have a profound impact upon marginal isostatic uplift and the evolution of offshore sedimentary basins. Yet at present there is only a rudimentary understanding of the geomorphic mechanisms capable of driving such prolonged escarpment retreat. In this study we present a nonlinear, two-dimensional landscape evolution model tha tis used to asses the necessary and sufficient conditions for long-term retreat of a rift-generated escarpment. The model represents topography as a grid of cells, with drainage networkds evolving as water flows across the grid in the direction of steepest descent. The model accounts for sediment production by weathering, fluvial sediment transport, bedrock channel erosion, and hillslope sediment transport by diffusive mechanisms and by mass failure. Numerical experiments presented explore the effects of different combinations of erosion processes and of dynamic coupling between denudation and flexural isostatic uplift. Model results suggest that the necessary and sufficient conditions for long-term escarpment retreat are (1) incising bedrock channels in which the erosion rate increases with increasing drainage area, so that the channels steepen and propagate headward; (2) a low rate of sediment production relative to sediment transport efficiency, which promotes relief-generating processes over diffusive ones; (3) high continental elevation, which allows greater freedom for fluvial dissection; and (4) any process, including flexural isostatic uplift, that helps to maintain a drainage divide near an escarpment crest. Flexural isostatic uplift also facilitates escarpment, thereby increasing channel gradients and accelerating erosion which in turn generates additional isostatic uplift. Of all the above conditions, high continental elevation is common ot most rift margin escarpments and may ultimately be the most important factor.

Tucker, Gregory E.↗

Gravity wave propagation in a diffusively separated atmosphere with height-dependent collision frequencies

Numerical calculations of gravity wave propagation through a two-component thermosphere with vertically-varying inter-species collision frequencies are performed using a direct integration method. Reflection of upgoing gravity wave energy into downgoing gravity waves appears to be small though nonnegligible at typical thermospheric periods and wavelengths. Coupling into diffusion waves is completely insignificant in most cases and may be determined in part by the relative abundance of the minor species. Temperature perturbation differences between species can be greater than 10% above 300 km altitude even though coupling into diffusion waves is small. Diffusive dissipation of gravity wave energy can be significant in the lower thermosphere and may be comparable to dissipation by viscosity and heat conduction below about 225 km altitude. Heating by diffusive dissipation may also be comparable to direct heating by solar radiation in the lower thermosphere. As a result of dissipative filtering of long-period and short-wavelength waves and the sensitivity of He-N2 density perturbation phase differences to diffusion, the AE-C satellite wave observations can only be fit by internal gravity waves with periods no greater than about 25 min and horizontal wavelengths no less than about 200 km.

Del Genio, A. D.↗

Numerical calculations of turbulent reacting flow in a gas-turbine combustor

A numerical study for confined, axisymmetrical, turbulent diffusion flames is presented. Local mean gas properties are predicted by solving the appropriate conservation equations in the finite difference form with the corresponding boundary conditions. The k-epsilon two-equation turbulence model is employed to describe the turbulent nature of the flow. A two-step kinetic model is assumed to govern the reaction mechanism. The finite reaction rate is the smaller of an Arrhenius type of reaction rate and a modified version of eddy-breakup model. Reasonable agreement is observed between calculations and measurements, but to obtain better agreement, more work is needed on improvements of the above mathematical models. However, the present numerical study offers an improvement in the analysis and design of the gas turbine combustors.

Lin, Chin-Shun↗

Effects of detailed droplet heating models on turbulent sprays vaporization behavior

The effects of three different droplet heating models on the vaporization history and internal structure of turbulent liquid fuel sprays were investigated. The models considered were the infinite-diffusion, diffusion-limit, and effective-conductivity models. A numerical solution for the models was developed and implemented in the KIVA-II computer code. Low temperature and relatively high temperature numerical studies were conducted. The low temperature computations were compared with existing experimental data. The comparisons showed that while the infinite-diffusion and diffusion-limit models respectively overpredicted and underpredicted the fuel vapor peak concentration and distribution in the combustor, the effective-conductivity model gave results that were in better agreement with measurements. A limited study for the high temperature case was performed due to lack of experimental data and predictions using the three models were compared with each other.

Mawid, M. A.↗

Triple flame structure and diffusion flame stabilization

The stabilization of diffusion flames is studied using asymptotic techniques and numerical tools. The configuration studied corresponds to parallel streams of cold oxidizer and fuel initially separated by a splitter plate. It is shown that stabilization of a diffusion flame may only occur in this situation by two processes. First, the flame may be stabilized behind the flame holder in the wake of the splitter plate. For this case, numerical simulations confirm scalings previously predicted by asymptotic analysis. Second, the flame may be lifted. In this case a triple flame is found at longer distances downstream of the flame holder. The structure and propagation speed of this flame are studied by using an actively controlled numerical technique in which the triple flame is tracked in its own reference frame. It is then possible to investigate the triple flame structure and velocity. It is shown, as suggested from asymptotic analysis, that heat release may induce displacement speeds of the triple flame larger than the laminar flame speed corresponding to the stoichiometric conditions prevailing in the mixture approaching the triple flame. In addition to studying the characteristics of triple flames in a uniform flow, their resistance to turbulence is investigated by subjecting triple flames to different vortical configurations.

Veynante, D.↗

Algorithms for the Fractional Calculus: A Selection of Numerical Methods

Many recently developed models in areas like viscoelasticity, electrochemistry, diffusion processes, etc. are formulated in terms of derivatives (and integrals) of fractional (non-integer) order. In this paper we present a collection of numerical algorithms for the solution of the various problems arising in this context. We believe that this will give the engineer the necessary tools required to work with fractional models in an efficient way.

Diethelm, K.↗

Theoretical approach to oxygen atom degradation of silver

Based on available Rutherford backscattering spectrometry (RBS), proton induced X-ray emission (PIXE) and ellipsometry data obtained on silver specimens subjected to atomic oxygen attack in low Earth orbit STS flight 41-G, a theory was developed to model the oxygen atom degradation of silver. The diffusion of atomic oxygen in a microscopically nonuniform medium is an essential constituent of the theory. The driving force for diffusion is the macroscopic electrochemical potential gradient developed between the specimen surface exposed to the ambient and the bulk of the silver specimen. The longitudinal electric effect developed parallel to the gradient is modified by space charge of the diffusing charged species. Lateral electric fields and concentration differences also exist due to the nonuniform nature of the medium. The lateral concentration differences are found to be more important than the lateral electric fields in modifying the diffusion rate. The model was evaluated numerically. Qualitative agreement exists between the kinetics predicted by the theory and kinetic data taken in ground-based experiments utilizing a plasma asher.

Fromhold, Albert T., Jr.↗

Computation of diffuse sky irradiance from multidirectional radiance measurements

Accurate determination of the diffuse solar spectral irradiance directly above the land surface is important in characterizing the reflectance properties of these surfaces, especially vegetation canopies. This determination is also needed to infer the net radiation budget of the earth-atmosphere system above these surfaces. An algorithm is developed here for the computation of hemispheric diffuse irradiance using the measurements from an instrument called PARABOLA, which rapidly measures upwelling and downwelling radiances in three selected wavelength bands. The validity of the algorithm is established from simulations. The standard reference data set of diffuse radiances of Dave (1978), obtained by solving the radiative transfer equation numerically for realistic atmospheric models, is used to simulate PARABOLA radiances. Hemispheric diffuse irradiance is estimated from a subset of simulated radiances by using the algorithm described. The algorithm is validated by comparing the estimated diffuse irradiance with the true diffuse irradiance of the standard data set. The validations include sensitivity studies for two wavelength bands (visible, 0.65-0.67 micron; near infrared, 0.81-0.84 micron), different atmospheric conditions, solar elevations, and surface reflectances. In most cases the hemispheric diffuse irradiance computed from simulated PARABOLA radiances and the true irradiance obtained from radiative transfer calculations agree within 1-2 percent. This technique can be applied to other sampling instruments designed to estimate hemispheric diffuse sky irradiance.

Ahmad, Suraiya P.↗

Turbulent diffusivity

A method for determining the turbulent diffusivity tensor is proposed which has application to problems such as the mixing in the interior of the sun for the solar neutrino problem, the chemical evolution of the galaxy, and grain sedimentation on the primitive solar nebula. A model for determining the turbulent energy spectral function, the eddy decorrelation rate, and the growth rate of the underlying instability is discussed. Good agreement is obtained between the present values for the coefficient of turbulent diffusivity and those obtained from both numerical simulation studies and laboratory data on shear turbulence.

Canuto, V. M.↗

Process modelling for Space Station experiments

Examined here is the sensitivity of a variety of space experiments to residual accelerations. In all the cases discussed the sensitivity is related to the dynamic response of a fluid. In some cases the sensitivity can be defined by the magnitude of the response of the velocity field. This response may involve motion of the fluid associated with internal density gradients, or the motion of a free liquid surface. For fluids with internal density gradients, the type of acceleration to which the experiment is sensitive will depend on whether buoyancy driven convection must be small in comparison to other types of fluid motion, or fluid motion must be suppressed or eliminated. In the latter case, the experiments are sensitive to steady and low frequency accelerations. For experiments such as the directional solidification of melts with two or more components, determination of the velocity response alone is insufficient to assess the sensitivity. The effect of the velocity on the composition and temperature field must be considered, particularly in the vicinity of the melt-crystal interface. As far as the response to transient disturbances is concerned, the sensitivity is determined by both the magnitude and frequency of the acceleration and the characteristic momentum and solute diffusion times. The microgravity environment, a numerical analysis of low gravity tolerance of the Bridgman-Stockbarger technique, and modeling crystal growth by physical vapor transport in closed ampoules are discussed.

Alexander, J. Iwan D.↗

The Effects of Flame Structure on Extinction of CH4-O2-N2 Diffusion Flames

The effects of flame structure on the extinction limits of CH4-O2-N2 counterflow diffusion flames were investigated experimentally and numerically by varying the stoichiometric mixture fraction Z(sub st), Z(sub st) was varied by varying free-stream concentrations, while the adiabatic flame temperature T(sub ad) was held fixed by maintaining a fixed amount of nitrogen at the flame. Z(sub st) was varied between 0.055 (methane-air flame) and 0.78 (diluted- methane-oxygen flame). The experimental results yielded an extinction strain rate K(sub ext) of 375/s for the methane-air flame, increasing monotonically to 1042/s for the diluted-methane-oxygen flame. Numerical results with a 58-step Cl mechanism yielded 494/s and 1488/s, respectively. The increase in K(sub ext) with Z(sub st) for a fixed T(sub ad) is explained by the shift in the O2 profile toward the region of maximum temperature and the subsequent increase in rates for chain-branching reactions. The flame temperature at extinction reached a minimum at Z(sub st) = 0.65, where it was 200 C lower than that of the methane-air flame. This significant increase in resistance to extinction is seen to correspond to the condition in which the OH and O production zones are centered on the location of maximum temperature.

Du, J.↗

Magnetic Reconnection: Theoretical and Observational Perspectives: Preface

Magnetic reconnection is a fundamental plasma-physical process by which energy stored in a magnetic field is converted, often explosively, into heat and the kinetic energy of the charged particles that constitute the plasma. It occurs in a variety of astrophysical settings, ranging from the solar corona to pulsar magnetospheres and winds, as well as in laboratory fusion experiments, where it is responsible for sawtooth crashes. First proposed by R.G. Giovanelli in the late I 940s as the mechanism responsible for solar flares, magnetic reconnection was invoked at the beginning of the space age to explain not just solar flares but also the transfer of energy, mass, and momentum from the solar wind to Earth's magnetosphere and the subsequent storage and release of the transferred energy in the magnetotai\. During the half century or so that has followed the seminal theoretical works by J.W. Dungey, P.A. Sweet, E.N. Parker, and H.E. Petschek, in-situ measurements by Earth-orbiting satellites and remote-sensing observations of the solar corona have provided a growing body of evidence for the occurrence of reconnection at the Sun, in the solar wind, and in the near-Earth space environment. The last thirty years have also seen the development of laboratory reconnection experiments at a number of institutions. In parallel with the efforts of experimentalists in both space and laboratory plasma physics, theorists have investigated, analytically and with the help of increasingly powerful MHD, hybrid, and kinetic numerical simulations, the structure of the diffusion region, the factors controlling the rate, onset, and cessation of reconnection, and the detailed physics that enables the demagnetization of the ions and electrons and the topological reconfiguration of the magnetic field. Moreover, the scope of theoretical reconnection studies has been extended well beyond solar system and laboratory plasmas to include more exotic astrophysical plasma systems whose strong (10(exp 14)-10(exp 15) G) magnetic fields require that models of reconnection in these systems incorporate quantum electrodynamical, special relativistic, and radiative effects. The papers collected in this topical issue of Space Science Reviews cover different aspects of recent theoretical and observational work on magnetic reconnection in solar and space physics, astrophysics, and laboratory plasma physics. They derive from presentations given at a workshop on magnetic reconnection held in the Yosemite National Park, February 8-12,2010. The intent of the workshop was to stimulate, through a combination of tutorial talks, shorter focused talks, and extensive informal discussions, an interdisciplinary dialogue among members of the different research communities working on the problem of magnetic reconnection. One of the motivating considerations for holding the workshop was its relevance to NASA's Magnetospheric Multiscale (MMS) mission, scheduled for launch in 2014. The four identically instrumented MMS spacecraft are designed to study reconnect ion in Earth's magnetosphere and, specifically, to probe the electron diffusion region in order to determine the microphysical processes that enable the change in the topology of the magnetic field. Building on the achievements of the multi spacecraft Cluster and THEMIS missions, MMS will use the magnetosphere as an astrophysical plasma laboratory in which to test, through in-situ measurement of the plasma, energetic particles, and electric and magnetic fields, various models and theories that have emerged during the past twenty years, a period of extraordinarily productive theoretical and observational work.

Lewis, W. S.↗

A numerical solution of the supersonic flow over a rearward facing step with transverse non-reacting hydrogen injection

The present investigation represents an application of computational fluid dynamics to a problem associated with the flow in the combustor region of a supersonic combustion ramjet engine (scramjet). The governing equations are considered, taking into account the Navier-Stokes equations, a molecular viscosity calculation, the molecular thermal conductivity, molecular diffusion, and a turbulence model. The employed numerical solution is patterned after the explicit, time-dependent, unsplit, predictor-corrector, finite-difference method given by MacCormack (1969). The calculation is concerned with the supersonic flow over a rearward-facing step with transverse H2 injection at conditions germane to the combustor region of a scramjet engine. The H2 jet acts as an effective body which essentially shields the primary flow from the rearward-facing step, thus substantially changing the wave pattern in the primary flow.

Berman, H. A.↗

Transient and diffusion analysis of Hg Cd Te

The solution of the one dimensional planar interface solidification problem is presented. Numerical solutions are applied to experimental solute profiles to determine the effective diffusion constant for Hg Co Te. Progress on the design and construction of a high gradient furnace is reported.

Clayton, J. C.↗