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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Vibrationally excited nitrogen in stable auroral red arcs and its effect on ionospheric recombination

The time-dependent continuity equations, including diffusion, were solved for the first six energy levels of molecular nitrogen for conditions in the thermosphere corresponding to stable auroral red (SAR) arcs. The results show that molecular nitrogen is excited vibrationally to the degree that the rate constant for the ionospheric loss process, O(+) + N2 yields NO(+) + N, is increased by as much as a factor of 7.6 at F2 region altitudes. It was found that deviations from the energetically equivalent Boltzmann distribution were large, causing the rate constant to be as much as 1.6 times the rate constant calculated for the Boltzmann distribution. These results indicate that SAR arc intensities as small as 58 R can produce noticeable increases in the ionosphere ion-atom interchange reaction rate and hence in the rate of loss of ionospheric electrons. It is suggested that the observed decrease of electron density in the F2 region in SAR arcs can probably be explained by enhanced reaction rates for ion-atom interchange between O(+) and N2 caused by vibrational excitation of molecular nitrogen by electron impact.

Newton, G. P.↗
NASA NTRS1977

Systems improved numerical differencing analyzer

Analyzer is software system that possesses capabilities which make it well suited for solving lumped-parameter representations of physical problems governed by diffusion-type equations, such as Fourier, Poisson, or La Place.

Fink, L. C.↗
NASA NTRS1980

Systems improved numerical differencing analyzer

Program solves physical problems governed by diffusion-type equations, provided that equations can be modeled by lumped-parameter representation. Program is used for thermal analysis, and could be adapted to solve Fourier, Poisson, and Laplace differential equations. Program is in FORTRAN IV and Assembler for execution on UNIVAC 1100-series or CYBER 175.

Source record↗
NASA NTRS1981

Parameterization of linear wave chemical transport in planetary atmospheres by eddy diffusion

Two-dimensional transport by linear waves is parameterized by wave velocity and temperature correlations and nonconservative properties of the atmosphere in a residual Eulerian framework. Explicit expressions derived for the eddy tracer flux vector are written in eddy diffusion tensor form, but it is demonstrated that strong coupling of chemistry and temperature through the rate constants, as is the case for ozone, requires additional terms proportional to the mixing ratio and not to its gradient. A variety of limiting cases are discussed. In the residual Eulerian framework the K(yz), K(zy) components are nonzero for a nonconservative atmosphere and/or a reactive tracer. Tracers have common transport coefficients in the limit where their chemistry is negligible in comparison to transience. It is shown that the usual one-dimensional vertical eddy diffusion continuity equation is rigorously valid only when the temperature dependence of the rate constants is negligible and the tracer has strong vertical stratification relative to its horizontal variability.

Strobel, D. F.↗
NASA NTRS1982

Jacchia-Lineberry upper atmosphere density model

A series of increasingly accurate models which are a careful blend of empirical and theoretical formulae were developed. The exospheric temperature is assumed to be a function of: (1) the average and daily variations in the solar flux, (2) the average and three hourly variations in the geomagnetic index, (3) the angle between the position vector and the axis of the unsymmetric atmospheric bulge, and (4) the angle between the position vector and the geomagnetic pole. The exospheric temperature is related to the density by the solution of the diffusion equilibrium equations for the different constituents of the atmosphere as a function of altitude. Other variations are modeled directly as changes in the density. They are: (1) changes due to the semiannual effect, and (2) changes due to the seasonal latitudinal effect. The causes for these variations are not exactly known but may be modeled sufficiently by empirical formulae. The Jacchia model is assumed to be valid over the altitude range of 90 to 2500 km. The residuals between the observed density from satellite drag observations and the computed densities show the mean relative error to be generally less than 10 percent with occasional peak errors near 50 percent.

Mueller, A. C.↗
NASA NTRS1983

Acceleration and transport processes - Verification and observations

The general problem of diffusive transport and acceleration of energetic charged particles is considered. The transport of solar-flare particles, solar modulation of galactic cosmic rays and shock acceleration processes on the solar wind are examined and observational tests are summarized. It is concluded that the basic diffusive transport equation is a useful approximation in situations like the solar wind, where turbulent scattering by magnetic irregularities is sufficient to maintain near isotropy. The application of this equation to the interstellar medium andd other, more distant astrophysical regimes is then discussed and implications for gamma-ray astrophysics are outlined. Finally the evidence for interstellar turbulence is reviewed and its consequences briefly discussed.

Jokipii, J. R.↗
NASA NTRS1988

SINDA - Systems Improved Numerical Differencing Analyzer

Capability to model fluid flow added. SINDA is software system developed to solve physical problems governed by diffusion-type equations mathematically modeled by lumped-parameter representations. Used as general thermal analyzer with resistor-and-capacitor-network representations. Adapted to wide range of problems represented by Fourier, Possion, Laplace, or other differential equations. Written in FORTRAN and Assembler.

Damico, Steve J.↗
NASA NTRS1988

Finite element procedures for time-dependent convection-diffusion-reaction systems

New finite element procedures based on the streamline-upwind/Petrov-Galerkin formulations are developed for time-dependent convection-diffusion-reaction equations. These procedures minimize spurious oscillations for convection-dominated and reaction-dominated problems. The results obtained for representative numerical examples are accurate with minimal oscillations. As a special application problem, the single-well chemical tracer test (a procedure for measuring oil remaining in a depleted field) is simulated numerically. The results show the importance of temperature effects on the interpreted value of residual oil saturation from such tests.

Tezduyar, T. E.↗
NASA NTRS1990

SINDA '85/FLUINT, Version 2.2

SINDA, Systems Improved Numerical Differencing Analyzer, software system developed for solving physical problems governed by diffusion-type equations modeled by lumped-parameter representations. Used as general thermal analyzer with resistor-and-capacitor-network representations, adapted to wide range of problems represented by such differential equations as Fourier, Poisson, or LaPlace. Software system updated. FLUINT, FLUid INTegrator, is advanced, one-dimensional fluid-analysis program solving equations of arbitrary fluid-flow networks. Written in FORTRAN and ASSEMBLER.

Cullimore, Brent↗
NASA NTRS1989

Tau approximation techniques for identification of coefficients in parabolic PDE

A variant of the Tau method, called the weak Tau method, is developed on the basis of the weak form of the PDE for use in least-squares parameter estimation; also presented is a suitable abstract convergence framework. The emphasis is on the theoretical framework that allows treatment of the weak Tau method when it is applied to a wide class of inverse problems, including those for diffusion-advection equations, the Fokker-Planck model for population dynamics, and damped beam equations. Extensive numerical testing of the weak Tau method has demonstrated that it compares quite favorably with existing methods.

Banks, H. T.↗
NASA NTRS1991

A simple estimate of funneling-assisted charge collection

Funneling is qualitatively discussed in detail, and a quantitative analysis is given for the total (time-integrated) collected charge. It is shown that, for an n(+)/p junction, the total collected charge QT is given by QT = (1 + mu-n/mu-p) QD + 2Q(diff), where QD is the charge initially liberated in the depletion region and Qdiff is charge collected by diffusion. This equation does not apply to very short ion tracks or to devices having a thin epilayer. A numerical example is presented, and a comparison is made with PISCES results.

Edmonds, Larry D.↗
NASA NTRS1992

Evolution of a magnetic flux tube in two-dimensional penetrative convection

Highly supercritical compressible convection is simulated in a two-dimensional domain in which the upper half is unstable to convection while the lower half is stably stratified. This configuration is an idealization of the layers near the base of the solar convection zone. Once the turbulent flow is well developed, a toroidal magnetic field B sub tor is introduced to the stable layer. The field's evolution is governed by an advection-diffusion-type equation, and the Lorentz force does not significantly affect the flow. After many turnover times the field is stratified such that the absolute value of B sub tor/rho is approximately constant in the convective layer, where rho is density, while in the stable layer this ratio decreases linearly with depth. Consequently most of the magnetic flux is stored in the overshoot layer. The inclusion of rotation leads to travelling waves which transport magnetic flux latitudinally in a manner reminiscent of the migrations seen during the solar cycle.

Jennings, R. L.↗
NASA NTRS1993

Particle injection and the structure of energetic-particle-modified shocks

A macroscopic 'self-consistent' nonlinear two-fluid model is developed for energetic-particle-modeled shocks. The model incorporates particle exchange between a thermal gas and an energetic population via a 'thermal leakage' mechanism. In the model, injection is regarded as a process whereby particles in the cooler thermal core of a total particle distribution are energized sufficiently to cross a 'momentum boundary' into the wings of the distribution. Particles in the wings are identified as 'energetic particles' or 'cosmic rays' and are assumed to propagate according to the diffusive transport equation of cosmic ray theory. Thermal particles are energized either as a result of adiabatic compression in a decelerating flow or directly via anomalous heating within a subshock. It is shown that the nature and model of particle injection directly affects the nature and dynamics of the shock, and hence the efficiency of particle acceleration, and that the shock itself regulates particle injection dynamically.

Zank, G. P.↗
NASA NTRS1994

On the dynamics of some grid adaption schemes

The dynamics of a one-parameter family of mesh equidistribution schemes coupled with finite difference discretisations of linear and nonlinear convection-diffusion model equations is studied numerically. It is shown that, when time marched to steady state, the grid adaption not only influences the stability and convergence rate of the overall scheme, but can also introduce spurious dynamics to the numerical solution procedure.

Sweby, Peter K.↗
NASA NTRS1993

Steady state and dynamical structure of a cosmic-ray-modified termination shock

A hydrodynamic model is developed for the structure of a cosmic-ray-modified termination shock. The model is based on the two-fluid equations of diffuse shock acceleration (Drury and Volk, 1981). Both the steady state structure of the shock and its interaction with outer heliospheric disturbances are considered. Under the assumption that the solar wind is decelerated by diffusing interstellar cosmic rates, it is shown that the natural state of the termination shock is a gradual deceleration and compression, followed by a discontinuous jump to a downstream state which is dominated by the pressure contribution of the cosmic rays. A representative model is calculated for the steady state which incorporates both interstellar cosmic ray mediation and diffusively accelerated anomalous ions through a proposed thermal leakage mechanism. The interaction of large-scale disturbances with the equilibrium termination shock model is shown to result in some unusual downstream structure, including transmitted shocks and cosmic-ray-modified contact discontinuities. The structure observed may be connected to the 2-kHz outer heliospheric radio emission (Cairns et al., 1992a, b). The time-dependent simulations also demonstrate that interaction with solar wind compressible turbulence (e.g., traveling interplanetary shocks, etc.) could induce the termination shock to continually fluctuate between cosmic-ray-dominated and gas-dynamic states. This fluctuation may represent a partial explanation of the galactic cosmic ray modulation effect and illustrates that the Pioneer and Voyager satellites will encounter an evolving shock whose structure and dynamic properties are strongly influence by the mediation of interstellar and anomalous cosmic rays.

Donohue, D. J.↗
NASA NTRS1996

Some Experiences with Nonoverlapping Schur Complement Parallel Preconditioning for CFD Calculations

In this work we consider solving matrices which arise from the discretization of advection-diffusion field equations on arbitrary triangulated domains using stabilized numerical methods. The talk will discuss several candidate matrix preconditioning algorithms based on the 2 x 2 block factorization induced by an apriori partitioning of the triangulated domain. Application of the 2 x 2 block preconditioner requires the formation and inversion of the Schur complement submatrix. We consider several strategies for simplifying this task: incomplete Schur complement factorizations, drop tolerance element filling, Schur complement probing, and localized Schur complement inversion. Numerical results will be shown comparing performance and efficiency of these approximations. The matrix preconditioner has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The remainder of the talk will show numerous examples in CFD to demonstrate the efficiency and robustness of the techniques.

Barth, Timothy J.↗
NASA NTRS1990

Thermal Analysis

Packaging engineers can avoid unnecessary and costly redesign of product through use of tools such as the SINDA '85/FLUINT software system for predicting thermal performance of electronic aviation equipment. SINDA '85/FLUINT is an acronym for Systems Improved Numerical Differencing Analyzer and Fluid Integrators. It is a program for solving the physical problems governed by flow and diffusion type equations and is most widely used as a general thermal analyzer.

Source record↗