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 325 records · Page 18

Plasma Diffusion in Self-Consistent Fluctuations

The problem of particle diffusion in position space, as a consequence ofeleclromagnetic fluctuations is addressed. Numerical results obtained with a self-consistent hybrid code are presented, and a method to calculate diffusion coefficient in the direction perpendicular to the mean magnetic field is proposed. The diffusion is estimated for two different types of fluctuations. The first type (resuiting from an agyrotropic in itiai setting)is stationary, wide band white noise, and associated to Gaussian probability distribution function for the magnetic fluctuations. The second type (result ing from a Kelvin-Helmholtz instability) is non-stationary, with a power-law spectrum, and a non-Gaussian probabi lity distribution function. The results of the study allow revisiting the question of loading particles of solar wind origin in the Earth magnetosphere.

Smets, R.↗

Effects of free stream turbulence intensity and integral length scale on heat transfer from a circular cylinder in crossflow

A model has been proposed for the momentum eddy diffusivity induced by free stream turbulence intensity and integral length scale. The eddy diffusivity model is applied to a cylinder situated in a uniform crossflow in the presence of free stream turbulence. Numerical solution of the governing momentum and energy equations with the proposed eddy diffusivity model yielded results for the skin friction coefficient and the Nusselt number. The numerical predictions of the present work are compared with experimental data and the agreement between the two is seen to be very good.

Gorla, R. S. R.↗

A numerical model of gravity wave breaking and stress in the mesosphere

The goal of the study is to calculate numerically the deceleration and heating caused by breaking gravity waves. The effect of the radiative dissipation of the wave is included as vertical-wavelength-dependent Newtonian cooling. The parameterization for zonal deceleration is extended by breaking gravity waves (Lindzen, 1981) to include the turbulent diffusion of heat and momentum. After describing the numerical model, the numerical results are presented and compared with the parameterizations in a noninteractive model of the mean zonal wind. Attention is then given to the transport of constituents by gravity waves and the attendant turbulent zone. It is noted that if gravity wave breaking were not an intermittent process, gravity wave stresses would produce an adiabatic mesosphere with a zonal mean velocity close to the phase speed of the breaking wave.

Schoeberl, M. R.↗

Analytic approximations in the study of the solar modulation of electrons

Numerical solutions to the transport equation of galactic cosmic rays in the interplanetary medium have been used to investigate the applicability of commonly used approximate analytic solutions for electrons (10 MeV to 10 GeV). We find that for a given cosmic-ray diffusion coefficient the force-field approximation is in reasonable agreement with the numerical solution at energies above 200 MeV, but deviates significantly at lower energies, depending on the shape of the interstellar electron spectrum. The diffusion-convection approximation agrees generally with the numerical solution within a factor of 2 over the entire energy range.

Cummings, A. C.↗

Finite-element analysis of turbulent flow in annular exhaust diffusers of gas turbine engines

A modified version of the Petrov-Galerkin weighted residual method coupled with a biquadratic finite element of the Lagrangian type was used to develop a finite-element model of the turbulent flow field in the annular exhaust diffuser of a gas turbine engine. The swirling flow field was analyzed with emphasis on the diffuser off-design operation. A comparison of the numerical results with experimental data shows that the model is applicable to moderately separating flows of the kind that are typically associated with the off-design performance of diffusing passages in gas turbines.

Baskharone, E. A.↗

Momentum deposition on Wolf-Rayet winds: Nonisotropic diffusion with effective gray opacity

We derive the velocity and mass-loss rate of a steady state Wolf-Rayet (WR) wind, using a nonisotropic diffusion approximation applied to the transfer between strongly overlapping spectral lines. Following the approach of Friend & Castor (1983), the line list is assumed to approximate a statistically parameterized Poisson distribution in frequency, so that photon transport is controlled by an angle-dependent, effectively gray opacity. We show the nonisotropic diffusion approximation yields good agreement with more accurate numerical treatments of the radiative transfer, while providing analytic insight into wind driving by multiple scattering. We illustrate, in particular, that multiple radiative momentum deposition does not require that potons be repeatedly reflected across substantial distances within the spherical envelope, but indeed is greatest when photons undergo a nearly local diffusion, e.g., through scattering by many lines closely spaced in frequency. Our results reiterate the view that the so-called 'momentum problem' of Wolf-Rayet winds is better characterized as an 'opacity problem' of simply identfying enough lines. One way of increasing the number of thick lines in Wolf-Rayet winds is to transfer opacity from saturated to unsaturated lines, yielding a steeper opacity distribution than that found in OB winds. We discuss the implications of this perspective for extending our approach to W-R wind models that incorporate a more fundamental treatment of the ionization and excitation processes that determine the line opacity. In particular, we argue that developing statistical descriptions of the lines to allow an improved effective opacity for the line ensemble would offer several advantages for deriving such more fundamental W-R wind models.

Gayley, Kenneth G.↗

An evaluation of computer codes for simulating the Galileo Probe aerothermal entry environment

The approaches of three computer flow field codes (HYVIS, COLTS, and RASLE), used to determine the Galileo Probe aerothermal environment and its effect on the design of the thermal protection system, are analyzed in order to resolve differences in their predicted results. All three codes account for the hypersonic, massively blown, radiation shock layers, characteristic of Jupiter entry. Significant differences, however, are evident in their solution procedures: the governing conservation equations, the numerical differencing methods, the governing physics (chemical, radiation, diffusion, and turbulence models), and the basic physical data (thermodynamic, transport, chemical, and spectral properties for atomic and molecular species). Solutions are compared for two near peak heating entry conditions for a Galileo Probe baseline configuration, having an initial mass of 242 kg and simulating entry into the Orton nominal atmosphere. The modern numerical methodology of COLTS and RASLE appear to provide an improved capability for coupled flow-field solutions.

Menees, G. P.↗

Prediction of the diffuse far-infrared flux from the galactic plane

A basic model and simple numerical relations useful for future far-infrared studies of the Galaxy are presented. Making use of recent CO and other galactic surveys, the diffuse far-infrared flux distribution from the galactic plane is predicted as a function of galactic longitude for longitudes between 4 and 90 deg; the far-infrared emissivity is predicted as a function of galactocentric distance. Future measurements of the galactic far-infrared flux would yield valuable information on the physical properties and distribution of dust and molecular clouds in the Galaxy, particularly the inner region.

Fazio, G. G.↗

Prediction of the diffuse far-infrared flux from the galactic plane

A basic model and simple numerical relations useful for future far-infrared studies of the Galaxy are presented. Making use of recent CO and other galactic surveys, the diffuse far-infrared flux distribution from the galactic planes is predicted as a function of galactic longitude for the region between 4 and 90 deg longitude; the far-infrared emissivity is predicted as a function of galactocentric distance. Future measurements of the galactic far-infrared flux would yield valuable information on the physical properties and distribution of dust and molecular clouds in the Galaxy, particularly the inner region.

Fazio, G. G.↗

Turbulent Flame Processes Via Diffusion Flame-Vortex Ring Interactions

Flame-vortex interactions are canonical configurations that can be used to study the underlying processes occurring in turbulent reacting flows. This configuration contains many of the fundamental aspects of the coupling between fluid dynamics and combustion that could be investigated with more controllable conditions than are possible under direct investigations of turbulent flames. Diffusion flame-vortex ring interaction contains many of the fundamental elements of flow, transport, combustion, and soot processes found in turbulent diffusion flames. Some of these elements include concentrated vorticity, entrainment and mixing, strain and nonequilibrium phenomena, diffusion and differential diffusion, partial premixing and diluent effects, soot formation and oxidation, and heat release effects. Such simplified flowfield allows the complex processes to be examined more closely and yet preserving the physical processes present in turbulent reacting flows. Furthermore, experimental results from the study of flame-vortex interactions are useful for the validation of numerical simulations and more importantly to deepen our understanding of the fundamental processes present in reacting flows. Experimental and numerical results obtained under microgravity conditions of the diffusion flame-vortex ring interaction are summarized in this paper. Results are obtained using techniques that include Flame Luminosity Imaging (FLI), Laser Soot-Mie Scattering (LSMS), Computational Fluid Dynamics and Combustion (CFDC), and Diode Laser Spectroscopy/Iterative Temperature with Assumed Chemistry (DLS/ITAC).

Dahm, Werner J. A.↗

Atlas of lunar radar maps at 70-cm wavelength

The intensity distribution of lunar radar echoes has been mapped for two-thirds of the earth-visible lunar surface at a wavelength of 70 cm. The depolarizing effects of the lunar surface were observed by simultaneously receiving the radar echoes in opposite polarizations. These echoes were mapped with areal resolutions of 25-100 sq km. Mapping with this resolution confirmed that the young craters have enhanced returns. A few craters were found to have enhanced echoes only from their rims. Backscattering differences were also observed between various areas within a mare, between different highland areas, and between maria and adjacent highlands. These scattering differences were interpreted with a simple model, which assumed that the surface backscattered with varying amounts of quasi-specular and diffuse power. Only an increase in the diffuse power was needed to give the numerical values of the enhancements.

Thompson, T. W.↗

A finite difference method of solving anisotropic scattering problems

A new method of solving radiative transfer problems is described including a comparison of its speed with that of the doubling method, and a discussion of its accuracy and suitability for computations involving variable optical properties. The method uses a discretization in angle to produce a coupled set of first-order differential equations which are integrated between discrete depth points to produce a set of recursion relations for symmetric and anti-symmetric angular sums of the radiation field at alternate depth points. The formulation given here includes depth-dependent anisotropic scattering, absorption, and internal sources, and allows arbitrary combinations of specular and non-Lambertian diffuse reflection at either or both boundaries. Numerical tests of the method show that it can return accurate emergent intensities even for large optical depths. The method is also shown to conserve flux to machine accuracy in conservative atmospheres

Barkstrom, B. R.↗

Solar particle propagation from 1 to 5 AU

The radial variation of the interplanetary propagation of solar cosmic rays is investigated on a grand scale by performing a statistical analysis of charged-particle data obtained by Pioneers 10 and 11 from launch to Jovian encounter (i.e., at distances of 1 to 5 AU). A numerical propagation model is adopted which includes diffusion, convection, and adiabatic deceleration, as well as a variable power-law diffusion coefficient. The study of solar particle events is carried out by analyzing individually each of five key parameters that are uniquely defined in each solar particle event: the time to maximum flux starting from the initial release at the sun, the anisotropy at the time to maximum flux, the flux at that time, the time width of the event at half the flux at the time to maximum, and the decay time. Combined results for the proton energy ranges from 3.4 to 5.2 MeV and 24 to 30 MeV are presented, the average radial diffusion coefficients between 1 and 5 AU are determined for each energy range, and it is found that the combination of a near-impusive injection with a temporally constant spectral index produces a reasonably good fit to the high-energy data but not as good a fit to the low-energy data. Implications of the results obtained are discussed for cosmic-ray propagation theory, solar modulation studies, and interplanetary acceleration of solar particles.

Zwickl, R. D.↗

Systems improved numerical differencing analyzer

Program solves physical problems governed by diffusion-type equations, provided that equations can be modeled by lumped-parameter representation. Program is used for thermal analysis, and could be adapted to solve Fourier, Poisson, and Laplace differential equations. Program is in FORTRAN IV and Assembler for execution on UNIVAC 1100-series or CYBER 175.

Source record↗

Application of a Galerkin finite element method to atmospheric transport problems

Numerical simulation of the movement of a contaminant within the atmosphere presents difficulties due to the multidimensionality of the problem and the fact that the horizontal transport is usually convection dominated, that the boundary conditions are mixed, and that both slow and fast atmospheric chemical reactions can be important. In this study, numerical experiments using a Crank-Nicolson Galerkin finite element method to solve the time-dependent partial differential equations demonstrate the applicability and accuracy of this method for the variety of conditions encountered in atmospheric pollutant modeling. The Crank-Nicolson Galerkin method using piecewise linear, piecewise cubic Hermite polynomials, and upwind finite elements is shown to accurately model the pure convection of initial wave forms. Numerical results studying the interactions of convection, diffusion, chemical reaction, pollutant removal, and the effects of contaminant emission source strength, source location and multiple sources are also presented.

Carmichael, G. R.↗

Solute partitioning under continuous cooling conditions as a cooling rate indicator

A model of solute partitioning in a finite body under conditions of continuous cooling is developed for the determination of cooling rates from concentration profile data, and applied to the partitioning of zirconium between ilmenite and ulvospinel in the Apollo 15 Elbow Crater rocks. Partitioning in a layered composite solid is described numerically in terms of concentration profiles and diffusion coefficients which are functions of time and temperature, respectively; a program based on the model can be used to calculate concentration profiles for various assumed cooling rates given the diffusion coefficients in the two phases and the equilibrium partitioning ratio over a range of temperatures. In the case of the Elbow Rock gabbros, the cooling rates are calculated from measured concentration ratios 10 microns from the interphase boundaries under the assumptions of uniform and equilibrium initial conditions at various starting temperatures. It is shown that the specimens could not have had uniform concentrations profiles at the previously suggested initial temperature of 1350 K. It is concluded that even under conditions where the initial temperature, grain sizes and solute diffusion coefficients are not well characterized, the model can be used to estimate the cooling rate of a grain assemblage to within an order of magnitude.

Onorato, P. I. K.↗

Conservative streamtube solution of steady-state Euler equations

This paper presents a new method for solving the steady state Euler equations. The method is similar to streamline curvature methods but has a conserative finite volume formulation which ensures correct shock capturing. Either wall position or wall pressure may be prescribed as boundary conditions, permitting both direct and inverse calculations. In supersonic applications the solution is obtained by space-marching while in subsonic and transonic applications iterative relaxation methods are used. Numerical results are given for: (1) supersonic diffuser with oblique shocks (direct calculation); (2) supersonic jet entering still reservoir (inverse calculation); (3) subsonic bump in a channel with 25 percent blockage (direct and inverse); (4) subsonic high-work turbine cascade (direct); and (5) transonic bump in a channel with 12 percent blockage (direct calculation).

Drela, M.↗