Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equation”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Computations of turbulent evaporating sprays

A computational study of turbulent evaporating sprays is reported. The major focus is to examine the sensitivity of the vaporization behavior of turbulent sprays to the transient liquid-phase processes. Three models considered to represent these processes are the thin skin, infinite diffusion, and diffusion limit models. Favre-averaged equations with k-epsilon-g turbulence model are employed for the gas phase. The Lagrangian approach with a stochastic separated flow method is used for the liquid phase where the effects of gas turbulence on droplet trajectories and interphase transport rates are considered using random-walk computations. Also the variable-property effects are considered in detail. Results indicate that, depending upon the boiling temperature and heat of vaporization of the fuel considered, the vaporization behavior of turbulent sprays may be quite sensitive to the modeling of transient liquid-phase processes. Thus, it is important that for most hydrocarbon fuels these processes be adequately represented in any comprehensive spray computations. The present results also provide further support to the conclusions of earlier studies which have been based on simplified spray configurations.

Aggarwal, S. K.↗
NASA NTRS1991

A simplified treatment of SiB's land surface albedo parameterization

The earlier presented surface albedo parameterization is simplified by assuming that the reflectance of direct solar radiation is a simple function of solar zenith angle. The function chosen contains three parameters that vary with vegetation type, greenness, and leaf area index. Tables of parameter values are presented. Using these tables, SiB's (Simple Biosphere model) absorbances of direct solar radiation can be reproduced with an average relative error of less than 0.5 percent. Finally, the direct reflectance function is integrated over zenith angle to produce an equation for the surface reflectance of diffuse radiation.

Koster, Randal D.↗
NASA NTRS1991

Sinc-Galerkin estimation of diffusivity in parabolic problems

A fully Sinc-Galerkin method for the numerical recovery of spatially varying diffusion coefficients in linear partial differential equations is presented. Because the parameter recovery problems are inherently ill-posed, an output error criterion in conjunction with Tikhonov regularization is used to formulate them as infinite-dimensional minimization problems. The forward problems are discretized with a sinc basis in both the spatial and temporal domains thus yielding an approximate solution which displays an exponential convergence rate and is valid on the infinite time interval. The minimization problems are then solved via a quasi-Newton/trust region algorithm. The L-curve technique for determining an approximate value of the regularization parameter is briefly discussed, and numerical examples are given which show the applicability of the method both for problems with noise-free data as well as for those whose data contains white noise.

Smith, Ralph C.↗
NASA NTRS1991

Convective instability in superposed fluid and anisotropic porous layers

The onset of thermal convection due to heating from below in a system consisting of a fluid layer overlying a porous layer with anisotropic permeability and thermal diffusivity is studied. The linear perturbation equations are solved numerically. It is found that the effects of anisotropy on the onset of thermal convection are most profound for small values of the depth ratio xi (ratio of fluid layer thickness to porous layer thickness), since in that case, the onset of convection corresponds to significant motion in both layers. For fixed vlues of the vertical permeability in the porous medium, decreasing the value of zeta (ratio of horizontal to vertical permeability) leads to stabilization of the superposed layer configuration because of increased resistance to motion in the porous medium. For larger values of xi, the onset of motion is increasingly confined to the fluid layer, with the transport of heat through the porous layer occurring primarily by conduction. Accordingly, the influence of zeta on the stability characteristics for larger xi is less significant than the effects of an anisotropic thermal conductivity.

Chen, Falin↗
NASA NTRS1992

Magnetic field decay in isolated neutron stars

Three mechanisms that promote the loss of magnetic flux from an isolated neutron star - Ohmic decay, ambipolar diffusion, and Hall drift - are investigated. Equations of motions are solved for charged particles in the presence of a magnetic field and a fixed background of neutrons, while allowing for the creation and destruction of particles by weak interactions. Although these equations apply to normal neutrons and protons, the present interpretations of their solutions are extended to cover cases of neutron superfluidity and proton superconductivity. The equations are manipulated to prove that, in the presence of a magnetic force, the charged particles cannot be simultaneously in magnetostatic equilibrium and chemical equilibrium with the neutrons. The application of the results to real neutron stars is discussed.

Goldreich, Peter↗
NASA NTRS1995

Solar structure inversion with LOWL data

Inversion results for the radial hydrostatic structure of the Sun, using six months of oscillation data obtained with the LOWL instrument, are presented. Both low and intermediate degree modes are used, thus avoiding the systematic errors that might have occurred in previous inversions by merging more than one data set. Using modes of between 0 deg and 90 deg and frequencies of between 1.5 mHz and 3.5 mHz, the variations with depth of the speed of sound, the density and the pressure were inferred for radii of between 0.05 and 0.85 stellar radius. It was found that in this region, the sound speed was within 0.15% of that of a model constructed using an equation of state that incorporated helium diffusion. The density difference between the Sun and the model was less than 0.8%. Given the small error bars on the inversion results, these differences are considered as being significant.

Basu, Sarbani↗
NASA NTRS1985

Diffusion and solubility of oxygen in silver

The diffusion and solubility of oxygen in Ag in the temperature range between 412 and 862 C was determined. The following interpolation formula was found for the solubility: L = 8.19.1/100.exp(-11 860/RT)Mol O2/g.At.Ag.at 1/.5. The process obeys the Sieverts square root law within the limits of error. The dissolution of oxygen in Ag may be accompanied by the dissociation of the oxygen molecules into atoms. The tests on Ag-foils reveal that below a temperature of about 500 C a higher solubility is simulated by the adsorption of oxygen. The diffusion coefficient of oxygen in silver obeys the following equation: D = 2.72.1/100.exp(-11 000/RT)sq cm/s. The relatively low activation energy of 11 kcal/g.At suggests that the diffusion of oxygen takes places over interstitial sites.

Eichenauer, W.↗
NASA NTRS1987

Downstream boundary effects on the frequency of self-excited oscillations in transonic diffuser flows

An investigation of downstream boundary effects on the frequency of self-excited oscillations in two-dimensional, separated transonic diffuser flows has been conducted numerically by solving the compressible, Reynolds-averaged, thin-layer Navier-Stokes equation with a two-equation turbulence model. It was found that the unsteady diffuser flowfields are very sensitive to the location of the downstream boundary. Extension of the diffuser downstream boundary significantly reduces the frequency and amplitude of oscillations for pressure, velocity and shock. Computational results suggest that the mechanism causing the self-excited oscillation changes from viscous convective wave dominated oscillations to inviscid acoustic wave dominated oscillations when the location of downstream boundary varies from 8.66 to 134.7 throat height. The existence of a suction slot in the experimental setup obscures the physical downstream boundary and, therefore, presents a difficulty for quantitative comparisons between computation and experiment.

Hsieh, T.↗
NASA NTRS1984

Investigation of Aeronomic Processes on the Basis of Ionospheric Sporadic E

A method by means of which the characteristics of atmospheric turbulence can be determined is discussed. Turbulent diffusivity deduced from ionospheric sporadic E in different altitudes above the ionospheric stations at Bekescsaba and Juliusrah are given. Seasonal variations in different altitudes above Bekescsaba are shown. An attempt was made to determine the variation of the turbulent diffusivity during and after geomagnetic disturbances by computing the horizontal wind velocity from the thermal wind equation. It is argued that the change of the turbulent diffusivity below the turbopause is opposed to that above the turbopause.

Bencze, P.↗
NASA NTRS1974

Turbulent diffusion with memories and intrinsic shear

The first part of the present theory is devoted to the derivation of a Fokker-Planck equation. The eddies smaller than the hydrodynamic scale of the diffusion cloud form a diffusivity, while the inhomogeneous, bigger eddies give rise to a nonuniform migratory drift. This introduces an eddy-induced shear which reflects on the large-scale diffusion. The eddy-induced shear does not require the presence of a permanent wind shear and is intrinsic to the diffusion. Secondly, a transport theory of diffusivity is developed by the method of repeated-cascade and is based upon a relaxation of a chain of memories with decreasing information. The full range of diffusion consists of inertia, composite, and shear subranges, for which variance and eddy diffusivities are predicted. The coefficients are evaluated. Comparison with experiments in the upper atmosphere and oceans is made.

Tchen, C. M.↗
NASA NTRS1978

An integral solution to a nonlinear diffusion problem

The generalized Galerkin method (or the method of integral relations) is applied to the type of problem described by quasilinear parabolic equations. As an example the problem of nonlinear transient slab diffusion with a general reservoir boundary condition is worked out. The integral relations are given for an arbitrary number of strips, and solutions using up to seven strips have been obtained in order to investigate the convergence of the method.

Hussaini, M. Y.↗
NASA NTRS1989

Interpretation and modeling of the averaged equations for a fluid-solid flow

A self-consistent derivation of the conservation laws is given for flows of a fluid-solid mixture. A unified analytical framework for obtaining constitutive relations is provided. This analysis uses a control volume/control surface approach that is widely used in fluid mechanics. All terms in the governing equations and the constitutive relations are written in terms of the mass-weighted averages except solid concentration. It is believed that the mass-weighted average is the natural bridge between micromechanics and constitutive relations. The derived momentum equations contain terms that differ from all existing models except that of Prosperetti and Jones (1984). However, their assumptions are not needed here. Special attention is given to the solid phase pressure. The physical basis of the previously assumed form for this pressure (Givler 1987) becomes clear. A number of related phenomena are also discussed. These include the anti-diffusion and anisotropic normal stresses. The energy equations are also different from existing models.

Shen, Hayley H.↗
NASA NTRS1976

The ionization of a low-density intercloud medium by a single O star

The ionization structure of hydrogen in a cloud-containing low-density gas surrounding an O star is calculated, emphasizing how UV photons cause the ionization of hydrogen in regions which are totally obscured from the star by neutral clouds. Particular consideration is given to the ionization of such obscured regions by diffuse Lyman continuum radiation produced when H(+) recombines directly to the ground state. The assumed physical structure of the cloud-containing neighborhood of an O star is discussed, a method is outlined for expressing the emission and absorption properties of the neutral clouds and their shadows in terms of effective volume emissivities and absorption coefficients on a line of sight, and equations of radiative transfer are derived for diffuse and stellar radiation in a cloud-containing H II region. Reradiation of Lyman continuum photons from the ionized boundary layer of a neutral cloud or shadow is analyzed. The expected emission measure across a large cloud-containing H II region is determined and found to compare well with that obtained from observations of H-alpha emission around runaway O stars. It is shown that the hydrogen in cloud shadows may be completely ionized by diffuse radiation within approximately half the Stroemgren radius, depending on the extent of that radius, the average frequency of the Lyman continuum radiation from the star, and the average cloud size.

Elmegreen, B. G.↗
NASA NTRS1996

A Systematic Solution Approach for Neutron Transport Problems in Diffuse Regimes

A systematic solution approach for the neutron transport equation, based on a least-squares finite-element discretization, is presented. This approach includes the theory for the existence and uniqueness of the analytical as well as of the discrete solution, bounds for the discretization error, and guidance for the development of an efficient multigrid solver for the resulting discrete problem. To guarantee the accuracy of the discrete solution for diffusive regimes, a scaling transformation is applied to the transport operator prior to the discretization. The key result is the proof of the V-ellipticity and continuity of the scaled least-squares bilinear form with constants that are independent of the total cross section and the absorption cross section. For a variety of least-squares finite-element discretizations this leads to error bounds that remain valid in diffusive regimes. Moreover, for problems in slab geometry a full multigrid solver is presented with V(1, 1)-cycle convergence rates approximately equal to 0.1, independent of the size of the total cross section and the absorption cross section.

Manteuffel, T. A.↗
NASA NTRS1969

Diffusion of rays in random media.

Fokker-Planck equation applied to probability density for finding rays larger than correlation distance during propagation in random medium

Liu, C. H.↗
NASA NTRS1972

An improved design method and experimental performance of two dimensional curved wall diffusers

A computer design program was developed to incorporate the suction slots in solving the potential flow equations with prescribed boundary conditions. Using the contour generated from this program two Griffith diffusers were tested having area ratios AR = 3 and 4. The inlet Reynolds number ranged from 600,000 to 7 million. It was found that the slot suction required for metastable operation depends on the sidewall suction applied. Values of slot suction of 8% of the inlet flow rate was required for AR = 4 with metastable condition, provided that enough sidewall suction was applied. For AR = 3, the values of slot suction was about 25% lower than those required for AR = 4. For nearly all unseparated test runs, the effectiveness was 100% and the exit flow was uniform. In addition to the Griffith diffusers, dump and cusp diffusers of comparable area ratios were built and tested. The results obtained from these diffusers were compared with those of the Griffith diffusers. Flow separation occurred in all test runs with the dump and cusp diffusers.

Yang, T.↗
NASA NTRS1988

Numerical simulation of axisymmetric turbulent flow in combustors and diffusors

A method for predicting turbulent flow in combustors and diffusers is developed. The Navier-Stokes equations, incorporating a turbulence kappa-epsilon model equation, were solved in a nonorthogonal curvilinear coordinate system. The solution applied the finite volume method to discretize the differential equations and utilized the SIMPLE algorithm iteratively to solve the differenced equations. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass, momentum and energy fluxes. The numerical accuracy was assessed using different finite differencing methods, i.e., hybrid, quadratic upwind and skew upwind, to represent the convection terms. Flows in different geometries of combustors and diffusers were simulated and results compared with experimental data and good agreement was obtained.

Yung, Chain Nan↗