SEARCH · Search NASA

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

All sources NASA NTRS DOE OSTI

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↗

NASA NTRS1997

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

An investigation is conducted on several numerical schemes for use in the computation of two-dimensional, spatially evolving, laminar variable-density compressible shear layers. Schemes with various temporal accuracies and arbitrary spatial accuracy for both inviscid and viscous terms are presented and analyzed. All integration schemes use explicit or compact finite-difference derivative operators. Three classes of schemes are considered: an extension of MacCormack's original second-order temporally accurate method, a new third-order variant of the schemes proposed by Rusanov and by Kutier, Lomax, and Warming (RKLW), and third- and fourth-order Runge-Kutta schemes. In each scheme, stability and formal accuracy are considered for the interior operators on the convection-diffusion equation U(sub t) + aU(sub x) = alpha U(sub xx). Accuracy is also verified on the nonlinear problem, U(sub t) + F(sub x) = 0. Numerical treatments of various orders of accuracy are chosen and evaluated for asymptotic stability. Formally accurate boundary conditions are derived for several sixth- and eighth-order central-difference schemes. Damping of high wave-number data is accomplished with explicit filters of arbitrary order. Several schemes are used to compute variable-density compressible shear layers, where regions of large gradients exist.

Kennedy, Christopher A.↗

NASA NTRS1997

Physiological Response of Plants Grown on Porous Ceramic Tubes

This research involves the manipulation of the root-zone water potential for the purposes of discriminating the rate limiting step in the inorganic nutrient uptake mechanism utilized by higher plants. This reaction sequence includes the pathways controlled by the root-zone conditions such as water tension and gradient concentrations. Furthermore, plant based control mechanisms dictated by various protein productions are differentiated as well. For the nutrients limited by the environmental availability, the kinetics were modeled using convection and diffusion equations. Alternatively, for the nutrients dependent upon enzyme manipulations, the uptakes are modeled using Michaelis-Menten kinetics. In order to differentiate between these various mechanistic steps, an experimental apparatus known as the Porous Ceramic Tube - Nutrient Delivery System (PCT-NDS) was used. Manipulation of the applied suction pressure circulating a nutrient solution through this system imposes a change in the matric component of the water potential. This compensates for the different osmotic components of water potential dictated by nutrient concentration. By maintaining this control over the root-zone conditions, the rate limiting steps in the uptake of the essential nutrients into tomato plants (Lycopersicon esculentum cv. Cherry Elite) were differentiated. Results showed that the uptake of some nutrients were mass transfer limited while others were limited by the enzyme kinetics. Each of these were adequately modeled with calculations and discussions of the parameter estimations provided.

Tsao, David↗

NASA NTRS1998

A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy

Stable and accurate interface conditions are derived for the linear advection-diffusion equation. The conditions are functionally independent of the spatial order of accuracy and rely only on the form of the discrete operator. We focus on high-order finite-difference operators that satisfy the summation-by-parts (SBP) property. We prove that stability is a natural consequence of the SBP operators used in conjunction with the new boundary conditions. In addition, we show that the interface treatments are conservative. New finite-difference operators of spatial accuracy up to sixth order are constructed: these operators satisfy the SBP property. Finite-difference operators are shown to admit design accuracy (p(sup th)-order global accuracy) when (p - 1)(sup th)-order stencil closures are used near the boundaries if the physical boundary conditions are implemented to at least p(sup th)-order accuracy. Stability and accuracy are demonstrated on the nonlinear Burgers' equation for an twelve-subdomain problem with randomly distributed interfaces.

Carpenter, Mark H.↗

NASA NTRS1997

A Thermal and Electrical Analysis of Power Semiconductor Devices

The state-of-art power semiconductor devices require a thorough understanding of the thermal behavior for these devices. Traditional thermal analysis have (1) failed to account for the thermo-electrical interaction which is significant for power semiconductor devices operating at high temperature, and (2) failed to account for the thermal interactions among all the levels involved in, from the entire device to the gate micro-structure. Furthermore there is a lack of quantitative studies of the thermal breakdown phenomenon which is one of the major failure mechanisms for power electronics. This research work is directed towards addressing. Using a coupled thermal and electrical simulation, in which the drift-diffusion equations for the semiconductor and the energy equation for temperature are solved simultaneously, the thermo-electrical interactions at the micron scale of various junction structures are thoroughly investigated. The optimization of gate structure designs and doping designs is then addressed. An iterative numerical procedure which incorporates the thermal analysis at the device, chip and junction levels of the power device is proposed for the first time and utilized in a BJT power semiconductor device. In this procedure, interactions of different levels are fully considered. The thermal stability issue is studied both analytically and numerically in this research work in order to understand the mechanism for thermal breakdown.

Vafai, Kambiz↗

NASA NTRS1996

Frequency Response of Pressure Sensitive Paints

An experimental method for measuring the frequency response of Pressure Sensitive Paints (PSP) is presented. These results lead to the development of a dynamic correction technique for PSP measurements which is of great importance to the advancement of PSP as a measurement technique. The ability to design such a dynamic corrector is most easily formed from the frequency response of the given system. An example of this correction technique is shown. In addition to the experimental data, an analytical model for the frequency response is developed from the one dimensional mass diffusion equation.

Winslow, Neal A.↗

NASA NTRS1999

Application of the Space-Time Conservation Element and Solution Element Method to One-Dimensional Advection-Diffusion Problems

Test problems are used to examine the performance of several one-dimensional numerical schemes based on the space-time conservation and solution element (CE/SE) method. Investigated in this paper are the CE/SE schemes constructed previously for solving the linear unsteady advection-diffusion equation and the schemes derived here for solving the nonlinear viscous and inviscid Burgers equations. In comparison with the numerical solutions obtained using several traditional finite-difference schemes with similar accuracy, the CE/SE solutions display much lower numerical dissipation and dispersion errors.

Wang, Xiao-Yen↗

NASA NTRS1998

A Parallel Non-Overlapping Domain-Decomposition Algorithm for Compressible Fluid Flow Problems on Triangulated Domains

This paper considers an algebraic preconditioning algorithm for hyperbolic-elliptic fluid flow problems. The algorithm is based on a parallel non-overlapping Schur complement domain-decomposition technique for triangulated domains. In the Schur complement technique, the triangulation is first partitioned into a number of non-overlapping subdomains and interfaces. This suggests a reordering of triangulation vertices which separates subdomain and interface solution unknowns. The reordering induces a natural 2 x 2 block partitioning of the discretization matrix. Exact LU factorization of this block system yields a Schur complement matrix which couples subdomains and the interface together. The remaining sections of this paper present a family of approximate techniques for both constructing and applying the Schur complement as a domain-decomposition preconditioner. The approximate Schur complement serves as an algebraic coarse space operator, thus avoiding the known difficulties associated with the direct formation of a coarse space discretization. In developing Schur complement approximations, particular attention has been given to improving sequential and parallel efficiency of implementations without significantly degrading the quality of the preconditioner. A computer code based on these developments has been tested on the IBM SP2 using MPI message passing protocol. A number of 2-D calculations are presented for both scalar advection-diffusion equations as well as the Euler equations governing compressible fluid flow to demonstrate performance of the preconditioning algorithm.

Barth, Timothy J.↗

NASA NTRS1998

Re-examination of a Classic Experiment to Measure the Positronium-Helium Cross Section

In 1975, before the advent of positronium beams, a clever experiment was carried out in an attempt to measure low-energy Ps-atom scattering cross-sections, especially that of Ps-He. A series of parallel thin plates was placed in the vessel containing the sample of gas, and positrons were allowed to enter the chamber from a radioactive source. The time spectrum of annihilation radiation was observed in the traditional way, and the pickoff annihilation due to thermalized Ps striking one of the plates was to be the special feature of the experiment. Using a diffusion analysis the authors derived cross- sections for several gases, and for helium the value obtained was sigma = 0.0166 pi(alpha)(sub o, sup 2). Even then this value was thought to be very much too small, while recent measurements and some calculations suggest a more reasonable value would be sigma = 3 to 13 pi(alpha)(sub o,sup 2). It has always been puzzling that an apparently well-designed experiment should give such an unsatisfactory result. We have now re-examined the analysis of the data with some interesting consequences which may explain the discrepancy. Two main observations support our re-analysis. First, we note that the mean free path for Ps-He scattering would be quite long if the cross-section were as small as that quoted above; the diffusion method itself would seem to be questionable. For the larger values, however, there would be no such problem. Second, it was assumed that when the annihilation rate had settled down and was following an exponential decay curve the asymptotic solution of the diffusion equation had been reached. We find, instead, that a superposition of exponentials can accurately represent the decay curve, consistent with the higher cross-section value.

Drachman, Richard J.↗

NASA NTRS2001

Ar Atmosphere: Implications for Structure and Composition of Mercury's Crust

We examine the possibilities of sustaining an argon atmosphere by diffusion from the upper 10 km of crust, and alternatively by effusion from a molten or previously molten area at great depth . Ar-40 in the atmospheres of the planets is a measure of potassium abundance in the interiors since Ar-40 is a product of radiogenic decay of K-40 by electron capture with the subsequent emission of a 1.46 eV gamma-ray. Although the Ar-40 in the earth's atmosphere is expected to have accumulated since the late bombardment, Ar-40 in surface-bounded exospheres is eroded quickly by photoionization and electron impact ionization. Thus, the argon content in the exospheres of the Moon, Mercury and probably Europa is representative of current effusion rather than accumulation over the lifetime of the body. Argon content will be a function of K content, temperature, grain size distribution, connected pore volume and possible seismic activity. Although Mercury and the Moon differ in many details, we can train the solutions to diffusion equations to predict the average lunar atmosphere. Then these parameters can be varied for Hermean conditions. Assuming a lunar crustal potassium abundance of 300 ppm, the observed argon atmosphere requires equilibrium between the argon production in the upper 9 Km of the moon (1.135 x 10(exp -3) cm(exp -3) s(exp -1)) and its loss. Hodges et al. conclude that this loss rate and the observed time variability requires argon release through seismic activity, tapping a deep argon source. An important observation is that the extreme surface of the Moon is enhanced in argon rather than depleted, as one would expect from outgassing of radiogenic argon. Manka and Michel concluded that ion implantation explains the surface enhancement of Ar-40. About half of the argon ions produced in the lunar atmosphere would return to the surface, where they would become embedded in the rocks. Similarly, at Mercury we expect the surface rocks to be enhanced in Ar-40 wherever the magnetosphere has been open over time. Thus the measurement of surface composition will reveal the long-term effects of solar wind-magnetosphere interaction. Additional information is contained in the original extended abstract.

Killen, R. M.↗

NASA NTRS1995

On the Calculation of Anisotropic Extinction Coefficients for Rigid Fibrous Ceramic Insulations

The specific anisotropic extinction coefficient e* couples the effective Radiative properties of a fibrous insulation into the radiation diffusion equation. This coefficient can be calculated using Mie scattering theory if fiber diameters, refractive indices and fiber orientation distributions are known. In general, fiber orientation distributions are not readily accessible and past calculations have considered fibers as either randomly distributed or normal to the heat flow direction. In certain rigid fibrous ceramic insulations neither of these cases apply well and a simple procedure is described for approximating e* from values calculated for the random and normal orientation cases. The intrinsic error associated with this scaling procedure is investigated. Numerical computations for several test structures and fiber materials show the average error to be less than 5% for net heat flux and radiation conductivity calculations.

Marshall, Jochen↗