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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1995

The Method of Space-time Conservation Element and Solution Element: Development of a New Implicit Solver

The method of space-time conservation element and solution element is a nontraditional numerical method designed from a physicist's perspective, i.e., its development is based more on physics than numerics. It uses only the simplest approximation techniques and yet is capable of generating nearly perfect solutions for a 2-D shock reflection problem used by Helen Yee and others. In addition to providing an overall view of the new method, we introduce a new concept in the design of implicit schemes, and use it to construct a highly accurate solver for a convection-diffusion equation. It is shown that, in the inviscid case, this new scheme becomes explicit and its amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, its principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme.

Chang, S. C.↗

NASA NTRS1982

Radiation effects studies for the high-resolution spectrograph

The generation and collection of charge carriers created during the passage of energetic protons through a silicon photodiode array are modeled. Pulse height distributions of noise charge collected during exposure of a digicon type diode array to 21 and 75 MeV protons were obtained. The magnitude of charge collected by a diode from each proton event is determined not only by diffusion, but by statistical considerations involving the ionization process itself. Utilizing analytical solutions to the diffusion equation for transport of minority carriers, together with the Vavilov theory of energy loss fluctuations in thin absorbers, simulations of the pulse height spectra which follow the experimental distributions fairly well are presented and an estimate for the minority carrier diffusion length L sub d is provided.

Smith, L. C.↗

NASA NTRS1984

Effect of growth parameters on compositional variations in directionally solidified HgCdTe alloys

A series of Hg(1-x)Cd(x)Te alloy crystals was grown by directional solidification with compositions ranging from x = 0.2 to x = 0.4. The measured axial compositional profiles were interpreted in terms of an exact numerical solution of the appropriate diffusion equation that takes into account both the variations of the segregation coefficient and solidification rate with composition. The solutions for all growth rates agree generally well with the experimental data, and confirm earlier observations that employed approximate analytical solutions. The effective mass diffusion coefficients showed no significant composition or melt temperature dependence. Slightly higher diffusion coefficients were obtained, however, for the highest growth rates.

Szofran, F. R.↗

NASA NTRS1975

Behaviour of nitric oxide trails deposited in the mesosphere and stratosphere

The transient behavior of a nitric oxide trail deposited at approximately 60 km altitude is studied by the solution of the appropriate multidimensional diffusion equation which includes terms representing the effects of wind shear. Similar analysis is then carried out for the situation in the stratosphere. Trail behavior is found to be relatively independent of altitude and background ozone, but strongly dependent on the magnitude of eddy diffusity and the initial nitric oxide concentration. The nitric oxide trail reacts with ambient ozone to form nitrogen dioxide. For a trail 100 m initial radius, an ozone hole will form to a maximum size in 4 to 6 hours and then decay. The overall recovery time of the atmosphere following the creation of the trail is less than 12 hours.

Eberstein, I. J.↗

NASA NTRS1980

The interplanetary transport of solar cosmic rays

Numerical solutions are presented for the propagation of solar cosmic rays in interplanetary space, including the effects of pitch-angle scattering and adiabatic focusing. The intensity-time profiles can be well fitted by a simple radial spatial diffusion equation with scattering mean-free path lambda(fit). The radial mean-free path so obtained is significantly larger than the true scattering mean-free path for low-rigidity particles due to both adiabatic focusing and the inapplicability of the diffusive approximation early in the event. The well-known discrepancy between lambda(fit) and the theoretical predictions may be resolved by these calculations.

Gombosi, T. I.↗

NASA NTRS2019

Permeation Rate Equations for Hydrogen and Deuterium in a Palladium-Silver Alloy

Mass transfer of a gas through a selective, solid membrane is an effective method for separation of desired species. This selective permeability is evident in the flow of hydrogen and the isotope deuterium through palladium-silver metal alloy media. In this study, based upon Sieverts’s law and the Arrhenius diffusion equation, an empirical correlation was developed to determine the steady-state permeation rate R, dependence on media temperature T, gas supply pressure p(sub S), and gas backpressure p(sub B) on the lower pressure side. Because of an extensive range of experimental conditions and complete reporting of raw data, the research by Ackerman and Koskinas was used as a source for data allowing empirical equation fitting. Unfortunately, those authors reported best-fit equations that poorly represented their own results. To improve the modeling of the original data and demonstrate the quality of the measurements, the current study develops improved hydrogen and deuterium permeation rate equations: P = 4.22×10(exp −6) A[exp(−704/T)](sq. root p(sub S) − sq. root p(sub B))/t for hydrogen P = 2.12×10(exp −6) A[exp(−468/T )](sq. root p(sub S) − sq. root p(sub B))/t for deuterium for values of cross-sectional area A (sq.cm), medium thickness t (cm), pressure p (psia), and temperature T (Kelvin), giving a permeation rate in mole/minute. These equations model permeation rate data more closely than do several other existing literature sources.

Smith, Phillip J.↗

NASA NTRS1984

On the acceleration of energetic ions in Jupiter's magnetosphere

Several aspects of the problem of high-energy ions in the Jovian magnetosphere are addressed. Voyager observations pertaining to the problem of high-energy ions in the magnetosphere are summarized, and the charge exchange emission of fast neutral sulfur and oxygen atoms and their subsequent recapture by electron impact, charge exchange, and photoionization is considered. Solutions are given to the diffusion equation assuming a source of ions injected with a gyroenergy corresponding to pickup in the middle and outer magnetosphere. It is concluded that no reasonable model parameters exist to produce the required steep spectra of the particle observations with only pickup and adiabatic radial diffusion included. A local acceleration mechanism based on nonadiabatic wave-particle interactions is needed. The assumptions and model predictions of stochastic acceleration by MHD turbulence for the Jovian magnetosphere are described. The model makes a specific correspondence between MHD wave spectrum properties and particle spectrum properties at energies above the Alfven energy.

Barbosa, D. D.↗

NASA NTRS1984

A spectral element method for fluid dynamics - Laminar flow in a channel expansion

A spectral element method that combines the generality of the finite element method with the accuracy of spectral techniques is proposed for the numerical solution of the incompressible Navier-Stokes equations. In the spectral element discretization, the computational domain is broken into a series of elements, and the velocity in each element is represented as a high-order Lagrangian interpolant through Chebyshev collocation points. The hyperbolic piece of the governing equations is then treated with an explicit collocation scheme, while the pressure and viscous contributions are treated implicitly with a projection operator derived from a variational principle. The implementation of the technique is demonstrated on a one-dimensional inflow-outflow advection-diffusion equation, and the method is then applied to laminar two-dimensional (separated) flow in a channel expansion. Comparisons are made with experiment and previous numerical work.

Patera, A. T.↗

NASA NTRS1995

The method of space-time and conservation element and solution element: A new approach for solving the Navier-Stokes and Euler equations

A new numerical framework for solving conservation laws is being developed. This new framework 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 conceptually simple and designed to overcome several key limitations of the above traditional methods. A two-level scheme for solving the convection-diffusion equation is constructed and used to illuminate the major differences between the present method and those previously mentioned. This explicit scheme, referred to as the a-mu scheme, has two independent marching variables.

Chang, Sin-Chung↗

NASA NTRS1976

Lossy radial diffusion of relativistic Jovian electrons

The radial diffusion equation with synchrotron losses is solved by the Laplace-transform method for near equatorially mirroring relativistic electrons. The evolution of a power-law distribution function is found, and the characteristics of synchrotron burnoff are stated in terms of explicit parameters for an arbitrary diffusion coefficient of a specific form. The peaking of the 10.4-cm volume emissivity from Jupiter at an L shell of about 1.8 provides an estimate of the diffusion coefficient in the radiation belts; one value is suggested as the appropriate modification, for an equatorial field strength of 4.2 G, of the Birmingham et al. (1974) result. Nonsynchrotron losses are included phenomenologically; from the phase-space densities reported by McIlwain and Fillius (1975), the particle lifetime is estimated. Asymptotic forms for the distribution in the strong synchrotron loss regime are provided.

Barbosa, D. D.↗

NASA NTRS1976

Lossy radial diffusion of relativistic Jovian electrons

The radial diffusion equation with synchrotron losses was solved by the Laplace transform method for near-equatorially mirroring relativistic electrons. The evolution of a power law distribution function was found and the characteristics of synchrotron burn-off are stated in terms of explicit parameters for an arbitrary diffusion coefficient. Emissivity from the radiation belts of Jupiter was studied. Asymptotic forms for the distribution in the strong synchrotron loss regime are provided.

Barbosa, D. D.↗

NASA NTRS1993

Recent developments in high order K-exact reconstruction on unstructured meshes

This paper presents recent improvements in high-order K-exact reconstruction on unstructured meshes. The new reconstruction procedures are incorporated into a basic upwind finite-volume scheme suitable for solving scalar advection-diffusion equations as well as the Euler and Navier-Stokes equations. Numerical calculations are performed comparing the present method with lower order accurate reconstruction procedures (piecewise constant and piecewise linear) and various competing technologies such as the fluctuation splitting method of Roe (1987) and Deconinck et al. (1992) and a system-variant of the streamline diffusion Petrov-Galerkin method developed by Hansbo (1991) and Hansbo and Johnson (1991). Five test problems are used in the numerical comparisons: scalar circular advection, transonic and supersonic Euler flow, laminar boundary-layer flow, and general compressible Navier-Stokes flow.

Barth, Timothy J.↗

NASA NTRS1999

Efficient Multi-Stage Time Marching for Viscous Flows via Local Preconditioning

A new method has been developed to accelerate the convergence of explicit time-marching, laminar, Navier-Stokes codes through the combination of local preconditioning and multi-stage time marching optimization. Local preconditioning is a technique to modify the time-dependent equations so that all information moves or decays at nearly the same rate, thus relieving the stiffness for a system of equations. Multi-stage time marching can be optimized by modifying its coefficients to account for the presence of viscous terms, allowing larger time steps. We show it is possible to optimize the time marching scheme for a wide range of cell Reynolds numbers for the scalar advection-diffusion equation, and local preconditioning allows this optimization to be applied to the Navier-Stokes equations. Convergence acceleration of the new method is demonstrated through numerical experiments with circular advection and laminar boundary-layer flow over a flat plate.

Kleb, William L.↗

NASA NTRS1999

Efficient Multi-Stage Time Marching for Viscous Flows via Local Preconditioning

Kleb, William L.↗

NASA NTRS1989

Radial diffusion and losses of energetic protons in the 5 to 12 R(S) region of Saturn's magnetosphere

Phase space density profiles for nearly equatorial mirroring energetic protons, derived from Pioneer 11 measurements during the 1979 Saturn encounter, were analyzed using solutions of the time-averaged radial diffusion equation in a dipolar planetary magnetic field. Then, several loss models were considered, ranging from a model that takes into account satellite absorption only to the one that accounts for satellite plus ring-E absorption plus added distributed losses, and the corresponding form of the time-averaged radial diffusion coefficient was determined for each case by a minimum-variance fit to the data-derived phase space densiy profiles. The model comparisons clearly indicated the need for distributed losses, in addition to satellite absorption, in order to explain the smooth decline of the experimental energetic proton phase space density profiles in the 5-12-Saturn radius region of Saturn's magnetosphere.

Hood, L. L.↗

NASA NTRS1985

The modulation features of the long-period cosmic ray variations in connection with the sign change of the general magnetic field of the Sun

On the basis of the model and experimental investigations, the spatial distribution of cosmic ray anisotropy for the different epochs of solar activity is studied. A solution is offered to the anisotropic diffusion equation with regard to the electromagnetic conditions for the periods of minimum and maximum solar activity and cosmic ray particle drift in a regular interplanetary magnetic field. It is shown that the long period changes amplitude and phase of the diurnal variations of cosmic rays is limited not only by convection and diffusion of particles but also by the drift effect before and after the sing change of the general magnetic field of the Sun. The calculated model is compared with the results obtained on the basis of an analysis of the experimental data from the neutron super monitor station, and it is shown that for the periods when the lines of magnetic force of the Sun come from the Northern Hemisphere the phase of the first harmonic diurnal variation is shifted forwards to an earlier time.

Iskra, K.↗

NASA NTRS1997

Adsorption and solvation of ethanol at the water liquid-vapor interface: a molecular dynamics study

The free energy profiles of methanol and ethanol at the water liquid-vapor interface at 310K were calculated using molecular dynamics computer simulations. Both alcohols exhibit a pronounced free energy minimum at the interface and, therefore, have positive adsorption at this interface. The surface excess was computed from the Gibbs adsorption isotherm and was found to be in good agreement with experimental results. Neither compound exhibits a free energy barrier between the bulk and the surface adsorbed state. Scattering calculations of ethanol molecules from a gas phase thermal distribution indicate that the mass accommodation coefficient is 0.98, and the molecules become thermalized within 10 ps of striking the interface. It was determined that the formation of the solvation structure around the ethanol molecule at the interface is not the rate-determining step in its uptake into water droplets. The motion of an ethanol molecule in a water lamella was followed for 30 ns. The time evolution of the probability distribution of finding an ethanol molecule that was initially located at the interface is very well described by the diffusion equation on the free energy surface.

NASA Discipline Exobiology↗

NASA NTRS1985

Upwind relaxation algorithms for the Navier Stokes equations

The development of upwind relaxation algorithms for obtaining efficient steady-state solutions to the compressible Navier-Stokes equations is described. The method is second-order accurate spatially and naturally disipative, using third-order flux splitting of the pressure and convective terms and second-order central differencing for shear and heat flux terms. A line Gauss-Seidel relaxation approach, shown to be unconditionally stable for model convection and diffusion equations, is used. The algorithm is demonstrated for several flows using the thin-layer form of the equations, including the problem of shock-induced separation over a flat plate.

Thomas, J. L.↗