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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

NASA NTRS1991

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

NASA NTRS1991

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

NASA NTRS1992

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

NASA NTRS1992

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

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↗