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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

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

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

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

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

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↗
NASA NTRS1993

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

New Developments in the Method of Space-Time Conservation Element and Solution Element-Applications to Two-Dimensional Time-Marching Problems

A new numerical discretization method for solving conservation laws is being developed. This new approach differs substantially in both concept and methodology from the well-established methods, i.e., finite difference, finite volume, finite element, and spectral methods. It is motivated by several important physical/numerical considerations and designed to avoid several key limitations of the above traditional methods. As a result of the above considerations, a set of key principles for the design of numerical schemes was put forth in a previous report. These principles were used to construct several numerical schemes that model a 1-D time-dependent convection-diffusion equation. These schemes were then extended to solve the time-dependent Euler and Navier-Stokes equations of a perfect gas. It was shown that the above schemes compared favorably with the traditional schemes in simplicity, generality, and accuracy. In this report, the 2-D versions of the above schemes, except the Navier-Stokes solver, are constructed using the same set of design principles. Their constructions are simplified greatly by the use of a nontraditional space-time mesh. Its use results in the simplest stencil possible, i.e., a tetrahedron in a 3-D space-time with a vertex at the upper time level and other three at the lower time level. Because of the similarity in their design, each of the present 2-D solvers virtually shares with its 1-D counterpart the same fundamental characteristics. Moreover, it is shown that the present Euler solver is capable of generating highly accurate solutions for a famous 2-D shock reflection problem. Specifically, both the incident and the reflected shocks can be resolved by a single data point without the presence of numerical oscillations near the discontinuity.

Chang, Sin-Chung↗
NASA NTRS1995

Extensional flow convecting a reactant undergoing a first order homogeneous reaction and diffusional mass transfer from a sphere at low to intermediate Peclet and Damkohler numbers

Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 less than or equal to Pe less than or equal to 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 less than or equal to Da less than or equal to 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r, Theta) = Sum of c(n)(r)P(n)(cos Theta) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions Pn(cos Theta), the variable coefficient PDE can be reduced to a system of N + 1 second order ODEs for the radial functions c(sub n)(r), n = 0, 1, 2,..., N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N + 1) first order differential equations as a two-point boundary value problem on 1 less than or equal to r less than or equal to r(sub infinity). Convergence of the expansion for a specific value of N can thus be established and provides 'spectral' behavior as well as the full concentration field c(r, Theta).

Shah, N. Y.↗
NASA NTRS1995

Application of the space-time conservation element and solution element method to two-dimensional advection-diffusion problems

The existing 2-D alpha-mu scheme and alpha-epsilon scheme based on the method of space-time conservation element and solution element, which were constructed for solving the linear 2-D unsteady advection-diffusion equation and unsteady advection equation, respectively, are tested. Also, the alpha-epsilon scheme is modified to become the V-E scheme for solving the nonlinear 2-D inviscid Burgers equation. Numerical solutions of six test problems are presented in comparison with their exact solutions or numerical solutions obtained by traditional finite-difference or finite-element methods. It is demonstrated that the 2-D alpha-mu, alpha-epsilon, and nu-epsilon schemes can be used to obtain numerical results which are more accurate than those based on some of the traditional methods but without using any artificial tuning in the computation. Similar to the previous 1-D test problems, the high accuracy and simplicity features of the space-time conservation element and solution element method have been revealed again in the present 2-D test results.

Wang, Xiao-Yen↗
NASA NTRS1994

Perpendicular electron heating by absorption of auroral kilometric radiation

We investigate the possibility of perpendicular heating of electrons and the generation of '90 deg -electron conics' by particle diffusion in velocity space due to wave-particle interaction with intense auroral kilometric radiation. This interaction is made possible by the downward shift in the R-X cutoff below the electron cyclotron frequency that occurs in the presence of warm plasma. We stimulate this condition and solve the diffusion equation using a finite difference algorithm. The results show strong perpendicular electron heating and indicate that the main characteristics of an electron conic distribution can be reproduced under these conditions.

Morgan, D. D.↗
NASA NTRS1994

Nonaxisymmetric evolution in protostellar disks

We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.

Laughlin, Gregory↗
NASA NTRS1995

One-dimensional hybrid satellite track model for the Dynamics Explorer 2 (DE 2) satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. (1985) satellite track code but also includes elements of Roble et al. (1987b) global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the Mass Spectrometer Incoherent Scatter 1990 (MSIS-90) model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule heating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗
NASA NTRS1994

The effect of water on thermal stresses in polymer composites

The fundamentals of the thermodynamic theory of mixtures and continuum thermochemistry are reviewed for a mixture of condensed water and polymer. A specific mixture which is mechanically elastic with temperature and water concentration gradients present is considered. An expression for the partial pressure of water in the mixture is obtained based on certain assumptions regarding the thermodynamic state of the water in the mixture. Along with a simple diffusion equation, this partial pressure expression may be used to simulate the thermostructural behavior of polymer composite materials due to water in the free volumes of the polymer. These equations are applied to a specific polymer composite material during isothermal heating conditions. The thermal stresses obtained by the application of the theory are compared to measured results to verify the accuracy of the approach.

Sullivan, Roy M.↗
NASA NTRS1996

New Interpretation of the Wigner Function

I define a two-sided or forward-backward propagator for the pseudo-diffusion equation of the 'squeezed' Q function. This propagator leads to squeezing in one of the phase-space variables and anti-squeezing in the other. By noting that the Q function is related to the Wigner function by a special case of the above propagator, I am led to a new interpretation of the Wigner function.

Daboul, Jamil↗
NASA NTRS1995

Numerical Modeling of High-Temperature Corrosion Processes

Numerical modeling of the diffusional transport associated with high-temperature corrosion processes is reviewed. These corrosion processes include external scale formation and internal subscale formation during oxidation, coating degradation by oxidation and substrate interdiffusion, carburization, sulfidation and nitridation. The studies that are reviewed cover such complexities as concentration-dependent diffusivities, cross-term effects in ternary alloys, and internal precipitation where several compounds of the same element form (e.g., carbides of Cr) or several compounds exist simultaneously (e.g., carbides containing varying amounts of Ni, Cr, Fe or Mo). In addition, the studies involve a variety of boundary conditions that vary with time and temperature. Finite-difference (F-D) techniques have been applied almost exclusively to model either the solute or corrodant transport in each of these studies. Hence, the paper first reviews the use of F-D techniques to develop solutions to the diffusion equations with various boundary conditions appropriate to high-temperature corrosion processes. The bulk of the paper then reviews various F-D modeling studies of diffusional transport associated with high-temperature corrosion.

Nesbitt, James A.↗
NASA NTRS1996

A Parallel Multilevel Spectral Element Scheme

A parallel multilevel strategy is developed using spectral (p) finite elements. Hierarchic bases are particularly well suited since the element matrices and vectors are nested and the multilevel projections easiliy performed. Since the basis degree is used to specify the multigrid level, an EBE strategy is natural br the multilevel technique. Results are presented for two candidate nonlinear elliptic transport problems: the augmented drift-diffusion equations of semiconductor device modeling and the stream function-vorticity equations of incompressible fluid dynamics.

Davis, M. B.↗
NASA NTRS1994

A Coupled Soil-Atmosphere Model of H2O2 on Mars

The Viking Gas Chromatograph Mass Spectrometer failed to detect organic compounds on Mars, and both the Viking Labeled Release and the Viking Gas Exchange experiments indicated a reactive soil surface. These results have led to the widespread belief that there are oxidants in the martian soil. Since H2O2 is produced by photochemical processes in the atmosphere of Mars, and has been shown in the laboratory to reproduce closely the Viking LR results, it is a likely candidate for a martian soil oxidant. Here, we report on the results of a coupled soil/atmosphere transport model for H202 on Mars. Upon diffusing into the soil, its concentration is determined by the extent to which it is adsorbed and by the rate at which it is catalytically destroyed. An analytical model for calculating the distribution of H202 in the martian atmosphere and soil is developed. The concentration of H202 in the soil is shown to go to zero at a finite depth, a consequence of the nonlinear soil diffusion equation. The model is parameterized in terms of an unknown quantity, the lifetime of H202 against heterogeneous catalytic destruction in the soil. Calculated concentrations are compared with a H202 concentration of 30 nmoles/cu cm, inferred from the Viking Labeled Release experiment. A significant result of this model is that for a wide range of H202 lifetimes (up to 105 years), the extinction depth was found to be less than 3 m. The maximum possible concentration in the top 4 cm is calculated to be approx. 240 nmoles/cu cm, achieved with lifetimes of greater than 1000 years. Concentrations higher than 30 nmoles/cu cm require lifetimes of greater than 4.3 terrestrial years. For a wide range of H202 lifetimes, it was found that the atmospheric concentration is only weakly coupled with soil loss processes. Losses to the soil become significant only when lifetimes are less than a few hours. If there are depths below which H202 is not transported, it is plausible that organic compounds, protected from an oxidizing environment, may still exist. They would have been deposited by meteors, or be the organic remains of past life.

Bullock, Mark A.↗
NASA NTRS1995

One-Dimensional Hybrid Satellite Track Model for the Dynamics Explorer 2 (DE 2) Satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. satellite track code but also includes elements of Roble et al. global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the MSIS-90 model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule beating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗