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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Features of a reattaching turbulent shear layer subject to an adverse pressure gradient

Experimental data have been obtained in an incompressible turbulent flow over a rearward-facing step with superimposed adverse pressure gradient. Mean velocities, Reynolds stresses and triple-products measured by a laser Doppler velocimeter are presented for two cases of adverse pressure gradient. Mixing lengths, eddy viscosities, production, convection, turbulent diffusion, and dissipation terms are extracted from the data. These data are compared with various mixing length and eddy-viscosity turbulence models. Numerical calculations incorporating the k-epsilon and the algebraic-stress turbulence models are compared with the data. When determining quantities of engineering interest, the modified algebraic-stress model (ASM) is a significant improvement over the unmodified ASM and the unmodified k-epsilon model

Driver, D. M.↗
NASA NTRS1985

The anisotropic radiative transfer problem in optically thick, strongly magnetized plasma - A comparison of results

The accuracy of four different methods for solving the coherent radiative transfer problem in a strongly magnetized plasma is compared. Specific attention is given to the case of a semi-infinite homogeneous plasma with a magnetic field perpendicular to the surface. The four approaches evaluated were: diffusion approximation schemes; the generalized H-function method of Silant'ev (1982) the Feautrier method of Nagel (1981) and the integral equation approach of Meszares and Bonazzola (1981). It is found that the coupled diffusion approximation schemes of Nagel (1981) and Kaminker et al. (1982) provided satisfactory results and were more efficient in many cases than direct numerical methods.

Pavlov, G. G.↗
NASA NTRS1984

Boundary layer effects on particle impaction and capture

The inertial impaction and deposition of small particles on larger bodies with viscous boundary layers are considered theoretically, in a detailed comment on a paper by Menguturk et al. (1983). Topics addressed include cushion effects, the dimensionless groups corresponding to the diameter range (3-6 microns) examined by Menguturk et al. in a numerical example, analogous effects of particle-gas energy and mass exchange in boundary layers, and the combined effects of particle inertia and diffusion. It is argued that the inertial effects can be characterized in terms of a body, boundary-layer, or sublayer Stokes number. In a reply by Menguturk et al., the focus is on the application of the theoretical model to the erosion of blade surfaces in large gas turbines; the Stokes number is found to be of limited practical value in these cases, because the particle motion is not primarily normal to the blade surfaces.

Rosner, D. E.↗
NASA NTRS1993

Hydrodynamical simulations of star-gas interactions in the interstellar medium with an external gravitational potential

We have calculated 2D models for the global structure of the interstellar medium in a disk galaxy where stars and gas are coupled through star formation, mass loss, stellar heating of the gas, and optically thin radiative cooling. In one set of models, the two dimensions are in the plane of a disk, and for one case, the total size of the region is 1 kpc square. Using a more accurate numerical approach, we confirm the results of Chiang and Bregman (1988) that H I forms into filaments or sheets separated by hotter, more diffuse gas. Although the filaments often have denser subregions, individual clouds of neutral gas are rare. We calculated the vertical distribution of the gas in another set of models where one dimension was in the plane and the other was perpendicular to it. A filament and hot gas network still occurs, although dense neutral complexes are more common and occur near the midplane. The hot gas often cools as it rises, leading to the highest velocity material having both positive and negative velocities. We find other features that are qualitatively similar to observations; for instance, H I 'worms' are produced.

Rosen, Alexander↗
NASA NTRS1997

Structure of Microgravity Transitional and Pulsed Jet Diffusion Flames

This paper describes results obtained in a study of pulsed gas jet diffusion flames to better characterize the recently observed vortex/flame interactions in microgravity transitional and turbulent diffusion flames, and to improve the understanding of large-scale structures in corresponding normal-gravity flames. In preparation for a space experiment, tests were conducted in the 5.18-Second Zero-Gravity Facility of the NASA Lewis Research Center. Both unpulsed and pulsed laminar flames were studied and numerical modeling of these flames was carried out for data comparison and model validation. In addition, complementary tests for a series of unpulsed flames were conducted on-board the NASA KC-135 research aircraft. The microgravity transitional and turbulent gas-jet diffusion flames have been observed to be dominated by large-scale disturbances, or structures. These structures first appear intermittently in the flame at Reynolds numbers (based on the cold jet injection properties) of about 2100. With increase in injection Reynolds number, the rate of intermittent disturbances increases until the generation becomes continuous at Reynolds numbers of 3000 and higher. The behavior of these structures depends upon the velocity and temperature characteristics of the jet/flame shear layer. These characteristics are different in normal gravity and microgravity.

Bahadori, M. Yousef↗
NASA NTRS1995

Numerical modeling of flame-balls in fuel-air mixtures

At low gravity, when buoyancy effects are small, flame-balls can be generated. These are stationary spherical structures whose existence appears to require a near-limit mixture, a small Lewis number and heat losses from radiation. It is our goal to combine computational modeling with existing experimental and theoretical studies (NASA) of these structures so that an improved understanding of flammability limits and near-limit phenomena will occur. The question of flammability limits is of fundamental importance and has long been examined. It is of great practical importance to predict, from first principles, a limit mixture strength that agrees with experimental values for the configuration at hand. Flame-balls provide an excellent configuration in which convective losses can be eliminated and the resulting stable solutions are produced from a diffusive, reactive and radiative balance. Although analytical modeling provides convincing evidence that the key physical ingredients of flame-balls have been identified, quantitative confirmation can only come from detailed numerical simulations. Our goal is to predict theoretically the mass fractions of the species and the temperature as functions of the independent coordinate r.

Smooke, Mitchell D.↗
NASA NTRS1977

Gravitational effects on chemically reacting laminar boundary layer flows over a horizontal flat plate

A theoretical study of the chemically reacting laminar boundary layer flow over a horizontal flat plate with gravitationally induced buoyant force is presented. A diffusion flame sheet model was used to describe the combustion process. The effects of gravity on the purely force convection flow can be characterized by a dimensionless coordinate quantity, which is involved in the generation of the governing equations. A numerical solution of the zero and first order governing equations subject to the appropriate physical boundary conditions was obtained. It is shown that the cross stream buoyancy induced body force acts effectively to produce a streamwise pressure gradient in the fluid adjacent to the plate surface. It is concluded that buoyancy plays an important role in boundary layer diffusion flames.

Lavid, M.↗
NASA NTRS1982

An iterative radiative transfer code for ocean-atmosphere systems

The details of an iterative radiative transfer code for computing the intensity and degree of polarization of diffuse radiation in models of the ocean-atmosphere system are described. The present code neglects the upwelling radiation from below the ocean surface and as such can be applied to the part of the spectrum where the absorption by water is strong. To establish the reliability of the numerical scheme and the computer code, the results are compared with those of Fraser and Walker (1968), Dave (1972), and Mullamaa (1964); they are found to be in excellent agreement. The computations also show that both the intensity and the degree of polarization of the upwelling diffuse radiation at the top of the atmosphere vary significantly when the rough ocean at the base of the atmosphere is replaced by a Lambertian surface that reflects the same energy as the rough ocean.

Ahmad, Z.↗
NASA NTRS2021

Modelling the Water and Carbon Dioxide Production Rates of Comet 67P/Churyumov-Gerasimenko

The European Space Agency Rosetta/Philae mission to Comet 67P/Churyumov–Gerasimenko in 2014–2016 is the most complete and diverse investigation of a comet carried out thus far. Yet, many physical and chemical properties of the comet remain uncertain or unknown, and cometary activity is still not a well–understood phenomenon. We here attempt to place constraints on the nucleus abundances and sublimation front depths of H 2 O and CO 2 ice, and to reconstruct how the nucleus evolved throughout the perihelion passage. We employ the thermophysical modelling code ‘Numerical Icy Minor Body evolUtion Simulator’, or NIMBUS, to search for conditions under which the observed H 2 O and CO 2 production rates are simultaneously reproduced before and after perihelion. We find that the refractories to water–ice mass ratio of relatively pristine nucleus material is μ ≈ 1, that airfall material has μ ≈ 2, and that the molar abundance of CO 2 relative H 2 O is near 30 per cent. The dust mantle thickness is typically ≲ 2 cm. The average CO 2 sublimation front depths near aphelion were ~ 3.8 m and ~ 1.9 m on the northern and southern hemispheres, respectively, but varied substantially with time. We propose that airfall material is subjected to substantial fragmentation and pulverisation due to thermal fatigue during the aphelion passage. Sub–surface compaction of material due to CO 2 activity near perihelion seems to have reduced the diffusivity in a measurable way.

67P/Churyumov–Gerasimenko↗
NASA NTRS2007

Transient Numerical Modeling of Catalytic Channels

This paper presents a transient model of catalytic combustion suitable for isolated channels and monolith reactors. The model is a lumped two-phase (gas and solid) model where the gas phase is quasi-steady relative to the transient solid. Axial diffusion is neglected in the gas phase; lateral diffusion, however, is accounted for using transfer coefficients. The solid phase includes axial heat conduction and external heat loss due to convection and radiation. The combustion process utilizes detailed gas and surface reaction models. The gas-phase model becomes a system of stiff ordinary differential equations while the solid phase reduces, after discretization, into a system of stiff ordinary differential-algebraic equations. The time evolution of the system came from alternating integrations of the quasi-steady gas and transient solid. This work outlines the numerical model and presents some sensitivity studies on important parameters including internal transfer coefficients, catalytic surface site density, and external heat-loss (if applicable). The model is compared to two experiments using CO fuel: (1) steady-state conversion through an isothermal platinum (Pt) tube and (2) transient propagation of a catalytic reaction inside a small Pt tube. The model requires internal mass-transfer resistance to match the experiments at lower residence times. Under mass-transport limited conditions, the model reasonably predicted exit conversion using global mass-transfer coefficients. Near light-off, the model results did not match the experiment precisely even after adjustment of mass-transfer coefficients. Agreement improved for the first case after adjusting the surface kinetics such that the net rate of CO adsorption increased compared to O2. The CO / O2 surface mechanism came from a sub-set of reactions in a popular CH4 / O2 mechanism. For the second case, predictions improved for lean conditions with increased external heat loss or adjustment of the kinetics as in the first case. Finally, the results show that different initial surface-species distribution leads to different steady-states under certain conditions. These results demonstrate the utility of a lumped two-phase model of a transient catalytic combustor with detailed chemistry.

Struk, Peter M.↗
NASA NTRS1978

Molecular evolution of contracting clouds - Basic methods and initial results

The relationship between the dynamics of the interstellar gas and the thermal and chemical effects associated with interstellar molecules and dust is investigated. The evolution of a rather massive isolated initially diffuse cloud under self-gravity is studied, using the equations of hydrodynamics; only radial motions are considered, and the heat, chemical, and radiative-transfer equations are solved simultaneously with the hydrodynamic equations. The relevant chemistry is described along with the thermal model, the radiative-transfer process, and the numerical methods employed. Results for a contracting cloud are discussed in terms of the problem of initial conditions, the dynamical evolution of the cloud, its chemical and thermal evolution, time scales, and column densities. It is shown that the chemical evolution of a massive contracting diffuse cloud is sensitive to such physical properties as temperature and ion abundances, that warm and cool versions of a typical cloud evolve differently, and that the physical origin of this effect is the level of heating due to H2 formation on interstellar dust grains.

Gerola, H.↗
NASA NTRS1974

Calculation of a realistic gradient as a function of radial distance from the sun

Numerical calculations of radial gradients of galactic cosmic ray protons have been performed for a model of the interplanetary cavity which includes processes of diffusion along the spiral interplanetary magnetic field lines, convection, and adiabatic energy change. The increasing spiral nature of the field lines with increasing distance from the sun results in the bulk of the modulation region being beyond the earth. The radial gradients for this model are found to be small near the earth and out to a few AU, increasing as one approaches the boundary.

Lezniak, J. A.↗
NASA NTRS1980

Composition of the nighttime ionospheric F 1 region near the magnetic equator

The effects of vertical E x B transport on NO(+), O2(+) and O(+) densities in the nighttime equatorial ionospheric F 1 region are investigated. Ion densities are calculated as functions of altitude, latitude and local time by the numerical solution of coupled, time-dependent ion continuity equations, taking into account production, loss by charge exchange and dissociative recombination and transport by diffusion and E x B drift. The results of the calculations are compared with measurements of NO(+), O2(+) and O(+) ion densities obtained at low altitudes by a mass spectrometer on board the Atmospheric Explorer C satellite, and are found to be consistent with the observations, suggesting that in the equatorial region, vertical transport by E x B drift is primarily responsible for producing the observed NO(+), O2(+) and O(+) density profiles. In addition, the reaction of O2(+) with N(4S) is found to be an important sink for O2(+) and a source of NO(+) ions. Implications of the observed and calculated near constancy of electron and ion densities with altitude when NO(+) is the dominant ion on the growth of large-scale irregularities are also considered

Anderson, D. N.↗
NASA NTRS1981

Mixed spectral/finite difference approximations for slightly viscous flows

A numerical scheme is presented, which employs a standard finite difference approximation for the viscous terms in high Reynolds number flows, and resorts to pseudo-spectral methods for the advection terms to greatly reduce the major source of numerical error without dramatically increasing computational cost. The spectral/finite difference (SFD) method evaluates the advection term and second-order differences to evaluate the diffusion term. The fully finite difference (FFD) method with second-order central differences on both terms is also used. The SFD method can handle strong shocks and can outperform the FFD method at moderate viscosity.

Zang, T. A.↗
NASA NTRS1989

A survey of upwind methods for flows with equilibrium and non-equilibrium chemistry and thermodynamics

Several versions of flux-vector split and flux-difference split algorithms were compared with regard to general applicability and complexity. Test computations were performed using curve-fit equilibrium air chemistry for an M = 5 high-temperature inviscid flow over a wedge, and an M = 24.5 inviscid flow over a blunt cylinder for test computations; for these cases, little difference in accuracy was found among the versions of the same flux-split algorithm. For flows with nonequilibrium chemistry, the effects of the thermodynamic model on the development of flux-vector split and flux-difference split algorithms were investigated using an equilibrium model, a general nonequilibrium model, and a simplified model based on vibrational relaxation. Several numerical examples are presented, including nonequilibrium air chemistry in a high-temperature shock tube and nonequilibrium hydrogen-air chemistry in a supersonic diffuser.

Grossman, B.↗
NASA NTRS1990

Numerical prediction of turbulent flow over airfoil sections with a new nonequilibrium turbulence model

A new nonequilibrium turbulence closure model has been developed for computing wall bounded two-dimensional turbulent flows. This two-layer eddy viscosity model was motivated by the success of the Johnson-King model in separated flow regions. The influence of history effects are described by an ordinary differential equation developed from the turbulent kinetic energy equation. The performance of the present model has been evaluated by solving the flow around three airfoils using the Reynolds time-averaged Navier-Stokes equations. Excellent results were obtained for both attached and separated turbulent flows about the NACA 0012 airfoil, the RAE 2822 airfoil, and the Integrated Technology A 153W airfoil. Based on the comparison of the numerical solutions with the available experimental data, it is concluded that the new nonequilibrium turbulence model accurately captures the history effects of convection and diffusion on turbulence.

Ahmed, S.↗
NASA NTRS1989

An embedded grid formulation applied to delta wings

An embedded grid algorithm for the Euler and/or Navier-Stokes equations is developed and applied to delta wings at high angles of attack in low speed flow. The Navier-Stokes code is an implicit, finite volume algorithm, using flux difference splitting for the convective and pressure terms and central differencing for the viscous and heat transfer terms. Calculations are compared with detailed experimental results over an angle of attack range up to and beyond the maximum lift coefficient, corresponding to vortex breakdown at the trailing edge, for a delta wing nominally of unit aspect ratio. The results indicate that the overall flowfield, including surface pressures, surface streamlines, and vortex trajectories, can be simulated accurately with the global grid version of the present algorithm. However, comparison of computed velocities and vorticity with experimentally measured off-body values at an angle of attack of 20.5 deg indicates the core region is substantially more diffuse in the computations than that measured with either a five-hole probe or a laser velocimeter. Embedded grids, used to improve the numerical discretization in the core region, are formulated within the framework of the implicit, upwind-biased multi-grid algorithm. Structured levels of local nested refinements are made. Three-dimensional results for both Euler and Navier-Stokes calculations are shown, with up to 3 levels of embedded refinement. The embedding procedure was effective in eliminating a crossflow secondary separation produced in the Euler solutions on coarse grids.

Thomas, James L.↗