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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2010

Ablation, Thermal Response, and Chemistry Program for Analysis of Thermal Protection Systems

In previous work, the authors documented the Multicomponent Ablation Thermochemistry (MAT) and Fully Implicit Ablation and Thermal response (FIAT) programs. In this work, key features from MAT and FIAT were combined to create the new Fully Implicit Ablation, Thermal response, and Chemistry (FIATC) program. FIATC is fully compatible with FIAT (version 2.5) but has expanded capabilities to compute the multispecies surface chemistry and ablation rate as part of the surface energy balance. This new methodology eliminates B' tables, provides blown species fractions as a function of time, and enables calculations that would otherwise be impractical (e.g. 4+ dimensional tables) such as pyrolysis and ablation with kinetic rates or unequal diffusion coefficients. Equations and solution procedures are presented, then representative calculations of equilibrium and finite-rate ablation in flight and ground-test environments are discussed.

Milos, Frank S.↗
NASA NTRS1976

The partially averaged field approach to cosmic ray diffusion

The kinetic equation for particles interacting with turbulent fluctuations is derived by a new nonlinear technique which successfully corrects the difficulties associated with quasilinear theory. In this new method the effects of the fluctuations are evaluated along particle orbits which themselves include the effects of a statistically averaged subset of the possible configurations of the turbulence. The new method is illustrated by calculating the pitch angle diffusion coefficient D sub Mu Mu for particles interacting with slab model magnetic turbulence, i.e., magnetic fluctuations linearly polarized transverse to a mean magnetic field. Results are compared with those of quasilinear theory and also with those of Monte Carlo calculations. The major effect of the nonlinear treatment in this illustration is the determination of D sub Mu Mu in the vicinity of 90 deg pitch angles where quasilinear theory breaks down. The spatial diffusion coefficient parallel to a mean magnetic field is evaluated using D sub Mu Mu as calculated by this technique. It is argued that the partially averaged field method is not limited to small amplitude fluctuating fields and is hence not a perturbation theory.

Jones, F. C.↗
NASA NTRS1985

Long-time asymptotics of a system for plasma diffusion

A system of parabolic nonlinear equations that describe the diffusion of a fully collisional plasma across a strong magnetic field is discussed. It is demonstrated that the solution to this system tends to a time asymptotic state which is of space-time separable form, theta(t)f(x). Furthermore, f(x) is independent of the initial conditions and theta(t) depends slightly on the initial conditions. The rate of decay of the temporal part is governed by a nonlinear eigenvalue problem. Since the equations are considered in a bounded domain we are able to analyze the effect of boundary conditions on the evolution of the system. Additional effects on radiation, heating, and particle injection can also be accounted for. Essential differences between the behavior of a fully-coupled system and a scalar equation are observed.

Rosenau, P.↗
NASA NTRS1984

Formation and dynamics of large-scale magnetic structures in the ionosphere of Venus

The formation and dynamics of large-scale magnetic structures in the ionosphere of Venus are examined. It is shown that such structures must be the result of steady state convection of interplanetary field lines into the ionosphere by the small amount of solar wind plasma (less than or approximately equal to 1-5 percent) absorbed by the planetary atmosphere below the ionopause, rather than isolated remnants of large fields persisting for long periods without connection to the solar wind induced current and convection pattern. In particular, it is demonstrated that the magnetic diffusion of such structures would result in their dissipation with time scales of 1-10 min, if they were not steady state structures in convective and diffusive equilibriuim. It is shown that the equations governing the diffusion of these magnetic structures are similar to those governing diffusion of a gas out of an enclosed chamber with a porous wall, and a simple analog is illustrated. The application of these results to magnetic fields of astrophysical plasmas is discussed.

Cloutier, P. A.↗
NASA NTRS1972

Relation of the diffuse reflectance remission function to the fundamental optical parameters.

The Kubelka-Munk equations describing the diffuse reflectance of a powdered sample were compared to equations obtained using a uniformly-sized rough-surfaced spherical particle model. The comparison resulted in equations relating the remission function and the Kubelka-Munk constants to the index of refraction, the absorption coefficient, and the average particle diameter of a powdered sample. Published experimental results were used to test the equation relating to the remission function to the fundamental optical parameters.

Simmons, E. L.↗
NASA NTRS2018

Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. I: Approach Based on the Far-Field Foldy Equations

In this paper, we revisit, with further enhancements and clarifications, the self-consistent first-principles approach developed previously for deriving the vector radiative transfer theory for a discrete random medium with a sparse concentration of particles. We specifically consider the case of a plane-parallel particulate layer embedded in an otherwise homogeneous unbounded medium. The solution method is based on the far-field Foldy equations, an order-of-scattering expansion for the total field derived under the Twersky approximation, the computation of the coherent field by assuming that the positions of the particles are uncorrelated, and the ladder approximation for the coherency dyadic. The latter yields an integral equation for the diffuse specific coherency dyadic, defined through an angular spectrum repre- sentation for the coherency dyadic, which in turn, gives the vector radiative transfer equation for the diffuse specific intensity column vector. We analyze specifically the computation of the coherent field for inhomogeneous particulate media and multiple species of particles, the continuous extension of the far-field representation to the near field, the Foldy approximation, and the Foldy integral equation for the coherent field. Finally, we discuss the transition from the vector to the scalar radiative transfer equation.

Doicu, Adrian↗
NASA NTRS1998

Multi-Component Diffusion with Application To Computational Aerothermodynamics

The accuracy and complexity of solving multicomponent gaseous diffusion using the detailed multicomponent equations, the Stefan-Maxwell equations, and two commonly used approximate equations have been examined in a two part study. Part I examined the equations in a basic study with specified inputs in which the results are applicable for many applications. Part II addressed the application of the equations in the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) computational code for high-speed entries in Earth's atmosphere. The results showed that the presented iterative scheme for solving the Stefan-Maxwell equations is an accurate and effective method as compared with solutions of the detailed equations. In general, good accuracy with the approximate equations cannot be guaranteed for a species or all species in a multi-component mixture. 'Corrected' forms of the approximate equations that ensured the diffusion mass fluxes sum to zero, as required, were more accurate than the uncorrected forms. Good accuracy, as compared with the Stefan- Maxwell results, were obtained with the 'corrected' approximate equations in defining the heating rates for the three Earth entries considered in Part II.

Sutton, Kenneth↗
NASA NTRS1977

The relative diffusion of strong magnetic fields and tenuous gases

General equations are given for treating the diffusion of a magnetic field in a fluid with electrical conductivity and a resistive diffusion coefficient. Consequences of nonnegligible gas motion are investigated by considering the one-dimensional problem where a tenuous isothermal gas with pressure expressed as a function of x-coordinate and time is enmeshed in a magnetic field in the y-direction. Approximate equations are obtained for the case of weak fluctuations in gas or field pressure, and the properties of a nonlinear diffusion equation are illustrated with the example of an initial thin slab of gas confined by a uniform field on either side. The advance of a gas front into a nonvanishing gas density is illustrated by using the progressive wave solution. Application of the results to the evolution of streamers and filaments observed in the solar corona is briefly discussed

Parker, E. N.↗
NASA NTRS2011

A Simple, Analytical Model of Collisionless Magnetic Reconnection in a Pair Plasma

A set of conservation equations is utilized to derive balance equations in the reconnection diffusion region of a symmetric pair plasma. The reconnection electric field is assumed to have the function to maintain the current density in the diffusion region, and to impart thermal energy to the plasma by means of quasi-viscous dissipation. Using these assumptions it is possible to derive a simple set of equations for diffusion region parameters in dependence on inflow conditions and on plasma compressibility. These equations are solved by means of a simple, iterative, procedure. The solutions show expected features such as dominance of enthalpy flux in the reconnection outflow, as well as combination of adiabatic and quasi-viscous heating. Furthermore, the model predicts a maximum reconnection electric field of E(sup *)=0.4, normalized to the parameters at the inflow edge of the diffusion region.

Hesse, Michael↗
NASA NTRS1975

Investigation B: Laminar confined coaxial entrance flow with heat generation

The results of a parametric study on the entrance flow region in a gas core nuclear reactor are presented. The physical system is modeled as laminar confined, coaxial flow with heat generation in the inner fluid. The governing equations include the boundary layer approximations and the assumptions of only radial radiative transport of energy represented as an energy diffusion term. The Von Mises transformation and a zeta transformation are used to transform the equations into nonlinear nonhomogeneous convective-diffusion equations. A unique combination of forward and backward difference equations which yields accurate results at moderate computational times, is used in the numerical method. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.

Bobba, G. K. M.↗
NASA NTRS1973

Vertical distribution of vibrational energy of molecular nitrogen in a stable auroral red arc and its effect on ionospheric electron densities

Previous solutions of the problem of the distribution of vibrationally excited molecular nitrogen in the thermosphere have either assumed a Boltzmann distribution and considered diffusion as one of the loss processes or solved for the energy level populations and neglected diffusion. Both of the previous approaches are combined by solving the time dependent continuity equations, including the diffusion process, for the first six energy levels of molecular nitrogen for conditions in the thermosphere corresponding to a stable auroral red arc. The primary source of molecular nitrogen excitation was subexcitation, and inelastic collisions between thermal electrons and molecular nitrogen. The reaction rates for this process were calculated from published cross section calculations. The loss processes for vibrational energy were electron and atomic oxygen quenching and vibrational energy exchange. The coupled sets of nonlinear, partial differential equations were solved numerically by employing finite difference equations.

Newton, G. P.↗
NASA NTRS1986

Migration of the separation point on a deforming cylinder

An iterative scheme of solving the Navier-Stokes equations for unsteady two-dimensional flow about a deforming and translating cylinder is given. In the kth iteration, the convected vorticity appears as the source term in an equation of transient diffusion of vorticity. A novel integral transform is used to reduce this transient vorticity equation to a k-dimensional heat equation. The bounded solution of this equation is obtained with a general method of superposition for problems involving a moving boundary. An equation describing the migration of the separation point on a deforming cylinder in unsteady cross flows is derived from the analytically obtained velocity field. The radial expansion of the cylinder surface is shown to hasten the separation time and to increase the separation angle. The results imply that in the steady flow past a body that is moving forward and sinking at constant rates, the locus of points at which the azimuthal component of skin friction changes sign originates on the leeward ray at a point downstream of the front tip.

Lin, S. P.↗
NASA NTRS1986

Improvement of the second- and third-moment modeling of turbulence: A study of Reynolds-stress closure model

Four parts of the Reynolds-stress closure modeling are reported: (1) improvement of the k and epsilon equaitons; (2) development of the third-moment transport equation; (3) formulation of the diffusion coefficient of the momentum equation by using the algebraic-stress model of turbulence; and (4) the application of the Reynolds-stress model to a heat exchanger problem. It was demonstrated that the third-moment transport model improved the prediction of the triple-velocity products in the recirculating and reattaching flow regions in comparison with the existing algebraic models for the triple-velocity products. Optimum values for empirical coefficients are obtained for the prediction of the backward-facing step flows. A functional expression is derived for the coefficient of the momentum diffusion by employing the algebraic-stress model. The second-moment closure is applied to a heat transfer problem. The computations for the flow in a corrugated-wall channel show that the second-moment closure improves the prediction of the heat transfer rates by 30% over the k - epsilon model.

Amano, R. S.↗
NASA NTRS1972

Resonant diffusion in strongly turbulent plasmas.

The effect of turbulent fluctuations on plasma particles is considered, and equations are derived which describe the evolution of macroscopic properties such as temperature and flow speed of the turbulent plasma. Initially, a diffusion equation for a single-particle distribution function averaged over an ensemble of plasmas is derived for an unmagnetized plasma. For the resonant diffusion in strongly turbulent plasmas, an ensemble of three-dimensional plasmas is considered with an approximately homogeneous and stationary distribution of random electromagnetic fluctuations. For each realization, the single-particle distribution function satisfies the Vlasov equation.

Birmingham, T. J.↗