Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion equations”

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 613 records · Page 34

Computation of the unsteady facilitated transport of oxygen in hemoglobin

The transport of a reacting permeant diffusing through a thin membrane is extended to more realistic dissociation models. A new nonlinear analysis of the reaction-diffusion equations, using implicit finite-difference methods and direct block solvers, is used to study the limits of linearized and equilibrium theories. Computed curves of molecular oxygen permeating through hemoglobin solution are used to illustrate higher-order reaction models, the effect of concentration boundary layers at the membrane interfaces, and the transient buildup of oxygen flux.

Davis, Sanford↗

Fourier analysis of finite element preconditioned collocation schemes

The spectrum of the iteration operator of some finite element preconditioned Fourier collocation schemes is investigated. The first part of the paper analyses one-dimensional elliptic and hyperbolic model problems and the advection-diffusion equation. Analytical expressions of the eigenvalues are obtained with use of symbolic computation. The second part of the paper considers the set of one-dimensional differential equations resulting from Fourier analysis (in the tranverse direction) of the 2-D Stokes problem. All results agree with previous conclusions on the numerical efficiency of finite element preconditioning schemes.

Deville, Michel O.↗

Electric currents in E-like planetary ionospheres

In this paper an MHD approach is used to consider the conduction of electric current in a lightly ionized gas, taking into account the gradients of pressure in the ion and electron gases, in addition to the electric field. The coefficients of electrical conductivity are found for each driver of current. New expressions for the components of heat dissipation associated with each driver of current are developed, which are fully consistent with kinetic theory. The relationship of the results to those obtained by kinetic theory is discussed. New components of currents associated with planetary equatorial electrojets are found. A new diffusion equation for magnetic induction is found, applicable in E-like regions of planetary ionospheres, and stellar photospheres.

Cole, K. D.↗

Electron heating in quasi-perpendicular shocks - A Monte Carlo simulation

To study the problem of electron heating in quasi-perpendicular shocks, under the combined effects of 'reversible' motion, in the shock electric potential and magnetic field, and wave-particle interactions a diffusion equation was derived, in the drift (adiabatic) approximation and it was solved by using a Monte Carlo method. The results show that most of the observations can be explained within this framework. The simulation has also definitively shown that the electron parallel temperature is determined by the dc electromagnetic field and not by any wave particle induced heating. Wave-particle interactions are effective in smoothing out the large gradients in phase space produced by the 'reversible' motion of the electrons, thus producing a 'cooling' of the electrons. Some constraints on the wave-particle interaction process may be obtained from a detailed comparison between the simulation and observations. In particular, it appears that the adiabatic approximation must be violated in order to explain the observed evolution of the perpendicular temperature.

Veltri, Pierluigi↗

Pitch angle diffusion of newborn ions due to intrinsic turbulence in the solar wind

The objective of the present study is to understand the interaction of the solar wind with newborn ions in far upstream regions of a comet where the level of intrinsic turbulence is moderately low. Based on the assumption that quasi-linear theory is adequate and applicable, the pitch angle diffusion process and the time evolution of the newborn ion distribution function are investigated. Numerical solutions to the quasi-linear diffusion equation, including the effect of resonance broadening and that of continuous creation of newborn ions due to the ionization process, are obtained under several assumptions and approximations. It is found that theoretical results are consistent with the Giotto observations recently reported by Neugebauer et al. (1989).

Ziebell, L. F.↗

Modeling of chemical vapor infiltration for ceramic composites reinforced with layered, woven fabrics

A homogeneous model is developed for the chemical vapor infiltration by one-dimensional diffusion into a system of layered plies consisting of woven tows containing bundles of filaments. The model predictions of the amount of deposition and the porosity of the sample as a function of time are compared with the predictions of a recent nonhomogeneous model with aligned holes formed by the weave. The nonhomogeneous model allows for diffusion through the aligned holes, into the spaces between plies, and into the gaps around filaments; i.e., three diffusion equations apply. Relative to the nonhomogeneous results, the homogeneous model underestimates the amount of deposition, since the absence of holes and spaces allows earlier occlusion of gaps around filaments and restricts the vapor infiltration.

Chung, Gui-Yung↗

Flux-limited diffusion in a scattering medium

A diffusion equation (FDT) is presented with a coefficient that reduces to the appropriate limiting form in the streaming and near thermodynamic limits for a moving fluid in which the dominant source of opacity is Thomson scattering. The present results are compared to those obtained with the corresponding equations for an absorptive medium. It is found that FDT for a scattering medium is accurate to better than less than about 17 percent over the range of optical depths of tau in the range of about 0 to 3.

Melia, Fulvio↗

Layer tracking, asymptotics, and domain decomposition

A preliminary report is presented on the work on the tracking of internal layers in a singularly-perturbed convection-diffusion equation. It is shown why such tracking may be desirable, and it is also shown how to do it using domain decomposition based on asymptotic analysis.

Brown, D. L.↗

Characterization of directionally solidified Hg(1-x)Zn(x)Se semiconducting alloys

Hg(1-x)Zn(x)Se alloys with composition between x = 0.08 and 0.115 were synthesized from elemental constituents and were resolidified using a Bridgman-Stockbarger growth technique. By performing precision mass-density measurements on selected wafers cut perpendicular to the growth axis, it was shown that the axial compositional profiles fit a numerical solution to the 1D diffusion equation which takes into account the variation of interface velocity with time. Infrared transmission-edge measurements performed on selected transverse slices from each ingot showed that the relative radial variations in composition decreased with decreasing growth rate. Van der Pauw measurements on selected wafers showed that the 10 exp 18/cu cm electron concentration, typical of as-grown crystals, could be reduced by approximately an order of magnitude by annealing in Se vapor.

Cobb, S. D.↗

Discrete observability and numerical quadrature

The authors consider the problem of approximate observability of a one-dimensional diffusion equation on a finite spatial domain with spatial point measurements. The problem of the optimal selection of the measurement points is considered under three conditions: (1) no preassigned measurement nodes; (2) one preassigned node and; (3) two preassigned nodes. The main observation is that the optimal choice is related to three classical procedures in numerical analysis: (1) Gaussian quadrature; (2) Radau quadrature and; (3) Lobatto quadrature. It is shown that the existence of the Radau and Lobatto quadrature is closely related to classical root locus theory.

Martin, Clyde F.↗

Subsurface emissions from Mercury - VLA radio observations at 2 and 6 centimeters

Radio observations of Mercury made with the VLA; once in 1986, and on two dates in February of 1988 are presented. These observations are the first to spatially map both hot regions associated with the theoretical hot poles. These 'hot poles' are separated by 180 deg and are a result of the unusual diurnal heating from Mercury's 3/2 spin-orbit resonance and eccentric orbit. The highest resolution data maps areas of the planet as small as 330 km. Maps of total intensity, brightness temperature, polarized intensity, fractional polarization, depolarization, and spectral index are included. It is found that the subsurface thermal emissions from Mercury are characteristic of blackbody reradiation from the solar insolation over a diurnal cycle. These observations to produce full-disk thermophysical models are used. The one-dimensional, time-dependent heat-diffusion equation for all observed disk elements at each epoch in order to constrain thermophsyical parameters and properties of the subsurface material are solved. Using typical lunar values for several of the parameters, it is possible to reproduce the temperature morphology and most of the observed temperature values. It is found that the best-fit models require a substantial contribution of the heat transport in the subsurface to be radiative in nature. The primary difficulty in the models is in predicting the observed temperature differences as a function of frequency.

Ledlow, Michael J.↗

Estimation of unknown variable parameters in moving boundary problems

The problem of estimating unknown variable parameters appearing in moving boundary problems is considered; these are specifically nonlinear diffusion equations defined on a moving spatial domain. A spline-based approximation method that results in a sequence of computationally tractable approximate parameter estimation problems has been developed. A convergence result is proved for a certain class of these moving boundary problems. The paper is concluded with a set of representative numerical examples.

Murphy, K. A.↗

Physics of pitch angle scattering and velocity diffusion. I - Theory

A general theory for the pitch angle scattering and velocity diffusion of particles in the field of a spectrum of waves in a magnetized plasma is presented. The test particle theory is used to analyze the particle motion. The form of diffusion surfaces is examined, and analytical expressions are given for the resonance width and bounce frequency. The resonance widths are found to vary strongly as a function of harmonic number. The resulting diffusion can be quite asymmetric with respect to pitch angle of 90 deg. The conditions for the onset of pitch angle scattering and energy diffusion are explained in detail. Some of the known shortcomings of the standard quasi-linear theory are also addressed, and ways to overcome them are shown. In particular, the often stated quasi-linear gap at 90 deg is found to exist only under very special cases. For instance, oblique wave propagation can easily remove the gap. The conditions for the existence of the gap are described in great detail. A new diffusion equation which takes into account the finite resonance widths is also discussed. The differences between this new theory and the standard resonance broadening theory is explained.

Karimabadi, H.↗

Adjustment of footprint correction for airborne flux mapping over the FIFE site

A meaningful interpretation of airborne flux estimates over nonuniform terrain must consider local advection from nonhomogeneous source distributions (footprint correction). An empirical procedure for footprint correction of airborne flux data obtained in the First ISLSCP Field Experiment (FIFE) 1989 is presented, based on simple analytical solutions of the diffusion equation which have been adjusted to approximate more realistic solutions. Optimized estimates of surface flux 'maps' derived from airborne observations are then constructed for maximum correlation between flux estimates and independently observed surface characteristics. This paper also includes an analysis of resolution for the given procedure, in terms of amplitude recovery of hypothetical square-wave distributions of surface flux density, and discusses its implications for the interpretation of the optimized estimates of FIFE surface flux maps. Results show that corrected spatial sink distributions are obtained by the current empirical procedure at the cost of considerable reduction in amplitude recovery of small-scale variations in surface flux density. The high spatial correlation observed between corrected flux maps and surface characteristics such as greenness or surface temperature excess is ascribed to the approximately correct positioning of the main flux gradients across the site, which can be expected to have been reproduced with amplitude recovery equal to or greater than 60 percent.

Schuepp, P. H.↗

Pickup protons and water ions at Comet Halley - Comparisons with Giotto observations

The cometary ion pickup process along the sun-comet line at Comet Halley is investigated using a quasi-linear diffusion model including both pitch angle and energy diffusion, adiabatic compression, and convective motion with the solar wind flow. The model results are compared with energetic ion distributions observed by instruments on board the Giotto spacecraft. The observed power spectrum index of magnetic turbulence (gamma) is 2-2.5. The present simulation shows that when gamma was 2, the calculated proton distributions were much more isotropic than the observed ones. The numerical solutions of the quasi-linear diffusion equations show that the isotropization of the pickup ion distribution, particularly at the pickup velocity, is not complete even close to the bow shock. Given the observed turbulence level, quasi-linear theory yields pickup ion energy distributions that agree with the observed ones quite well and easily produces energetic ions with energies up to hundreds of keV.

Ye, G.↗

The pulsar planet production process

Most plausible scenarios for the formation of planets around pulsars end with a disk of gas around the pulsar. The supplicant author then points to the solar system to bolster faith in the miraculous transfiguration of gas into planets. We here investigate this process of transfiguration. We derive analytic sequences of quasi-static disks which give good approximations to exact solutions of the disk diffusion equation with realistic opacity tables. These allow quick and efficient surveys of parameter space. We discuss the outward transfer of mass in accretion disks and the resulting timescale constraints, the effects of illumination by the central source on the disk and dust within it, and the effects of the widely different elemental compositions of the disks in the various scenarios, and their extensions to globular clusters. We point out where significant uncertainties exist in the appropriate grain opacities, and in the effect of illumination and winds from the neutron star.

Phinney, E. S.↗

Extrapolation techniques applied to matrix methods in neutron diffusion problems

A general matrix method is developed for the solution of characteristic-value problems of the type arising in many physical applications. The scheme employed is essentially that of Gauss and Seidel with appropriate modifications needed to make it applicable to characteristic-value problems. An iterative procedure produces a sequence of estimates to the answer; and extrapolation techniques, based upon previous behavior of iterants, are utilized in speeding convergence. Theoretically sound limits are placed on the magnitude of the extrapolation that may be tolerated. This matrix method is applied to the problem of finding criticality and neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron-diffusion equations is treated. Results for this example are indicated.

Mccready, Robert R↗

The thickness of eucritic crust in the HED parent body

Because cumulate eucrites are generally thought to be located at the lower part in eucritic crust on a diogenitic mantle in the HED parent body, the burial depth of cumulate eucrites gives information on the thickness of the eucritic crust. We estimated the burial depth and cooling rate of cumulate eucrites, Serra de Mage and Moore County on the basis of the width of augite lamellae and compositional gradients of Ca in pyroxenes by numerically solving the diffusion equation. We obtained the burial depth of the eucrites of 7-8 km and cooling rate of 0.00016-0.0002 degrees C/yr.

Miyamoto, M.↗