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 865 records · Page 48

On the dynamics of agglomerating ensembles of clouds

It is shown that for dynamical time scales longer than 10 million years the ensemble of interstellar cold diffuse clouds may be described by equations of mass and momentum conservation where the pressure has the form p/v-squared/ where /v-squared/ is the mass averaged square velocity dispersion. Equilibrium cloud spectra are derived and an equation of state obtained. The equation of state undergoes a phase transition at a critical density close to the local interstellar gas density. At higher densities, more massive clouds are formed and the cloud fluid softens. It is suggested that this mechanism may naturally regulate the interstellar gas density if star formation is very efficient in the cooler high-density systems.

Cowie, L. L.↗

Analytical description of charged particle transport along arbitrary guiding-field configurations

Focused particle transport is analyzed in terms of partial differential equations which are applicable to arbitrary spatial variations of the mean free path and the focusing length. A distribution function for transport is written in two forms of the Boltzmann equation and a coefficient of spatial diffusion is defined. Changes from a nearly coherent propagation of flare particles in regions of strong focusing near the sun to ordinary diffusion in the outer solar system are described as supercoherent transitions. A solution to the steady state problem is developed and results are presented for a monopolar guiding field configuration with isotropic diffusion gradually giving way to anisotropy in a region suggestive of a collimated beam of particles streaming away from its source.

Earl, J. A.↗

Ion temperature variations in the daytime high-latitude F region

The Schunk and Sojka (1981a, b) high latitude ionospheric model is improved through the inclusion of thermal conduction and diffusion terms in the ion energy equation, permitting the study of daytime, high latitude F layer temperature variations in a region poleward of the auroral oval. It is found that ion temperature variation with solar cycle, season and geomagnetic activity closely follows the neutral atomic oxygen variation, and that meridional electric fields of more than 40 mV/m can cause larger ion temperature changes than those due to solar cycle, seasonal or geomagnetic activity variations. In the presence of meridional electric fields, there is an upward flow of heat from the lower ionosphere that also acts to raise ion temperatures at high altitudes. Zonal electric fields affect ion temperature indirectly, through electron density changes.

Schunk, R. W.↗

The diffusivity of hydrogen in Nb stabilized stainless steel

The evolution of hydrogen from 347 stainless steel has been studied by using a real time dynamic technique under ultrahigh vacuum conditions. Auger electron spectroscopy was used to determine the surface composition as a function of time and temperature. The surface film on the electropolished samples was found to be approximately 15 A thick and consisted of a carbon-oxygen complex and a metal oxide (FexOy). Upon heating to 400 C, the carbon-oxygen complex desorbed as CO and the remaining oxygen and carbon began to incorporate. Also at this temperature sulfur began to diffuse out of the bulk to the surface and at approximately 800 C formed a complete monolayer. At 900 C, carbon and oxygen virtually disappeared, leaving the monolayer of sulfur as the only surface contaminant. The hydrogen diffusivity was found to follow closely the equation D = 7.01 x 10 to the -7th exp(-48.0/RT) sq m per second over the entire temperature range studied, thus indicating that hydrogen evolution is not significantly affected by the changing surface composition. The somewhat higher value of the diffusivity obtained in this work compared to past measurements in austenitic stainless steels may indicate the importance of sample preprocessing and ultrahigh vacuum conditions in minimizing the effects of surface layers.

Outlaw, R. A.↗

Losses in Channels with Increased External Turbulence

An approximate method for determining the effect of the level of turbulence on the aerodynamic characteristics of convergent and diffuser channels is examined. A momentum equation for the boundary layer is in the method, introducing external flow turbulence on the basis of experimental values of the coefficient of friction and the form factor. It is found that at significant levels of external turbulence, losses must be considered not only in the boundary layer but also in the central region of the channel.

Zaryankin, A. Y.↗

Group-kinetic theory of turbulence

The two phases are governed by two coupled systems of Navier-Stokes equations. The couplings are nonlinear. These equations describe the microdynamical state of turbulence, and are transformed into a master equation. By scaling, a kinetic hierarchy is generated in the form of groups, representing the spectral evolution, the diffusivity and the relaxation. The loss of memory in formulating the relaxation yields the closure. The network of sub-distributions that participates in the relaxation is simulated by a self-consistent porous medium, so that the average effect on the diffusivity is to make it approach equilibrium. The kinetic equation of turbulence is derived. The method of moments reverts it to the continuum. The equation of spectral evolution is obtained and the transport properties are calculated. In inertia turbulence, the Kolmogoroff law for weak coupling and the spectrum for the strong coupling are found. As the fluid analog, the nonlinear Schrodinger equation has a driving force in the form of emission of solitons by velocity fluctuations, and is used to describe the microdynamical state of turbulence. In order for the emission together with the modulation to participate in the transport processes, the non-homogeneous Schrodinger equation is transformed into a homogeneous master equation. By group-scaling, the master equation is decomposed into a system of transport equations, replacing the Bogoliubov system of equations of many-particle distributions. It is in the relaxation that the memory is lost when the ensemble of higher-order distributions is simulated by an effective porous medium. The closure is thus found. The kinetic equation is derived and transformed into the equation of spectral flow.

Tchen, C. M.↗

Finite element solution techniques for large-scale problems in computational fluid dynamics

Element-by-element approximate factorization, implicit-explicit and adaptive implicit-explicit approximation procedures are presented for the finite-element formulations of large-scale fluid dynamics problems. The element-by-element approximation scheme totally eliminates the need for formation, storage and inversion of large global matrices. Implicit-explicit schemes, which are approximations to implicit schemes, substantially reduce the computational burden associated with large global matrices. In the adaptive implicit-explicit scheme, the implicit elements are selected dynamically based on element level stability and accuracy considerations. This scheme provides implicit refinement where it is needed. The methods are applied to various problems governed by the convection-diffusion and incompressible Navier-Stokes equations. In all cases studied, the results obtained are indistinguishable from those obtained by the implicit formulations.

Liou, J.↗

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.↗

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.↗

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.↗

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↗

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↗

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↗

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.↗

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.↗

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.↗

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.↗

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.↗