Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boltzmann 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

A kinetic equation with kinetic entropy functions for scalar conservation laws

A nonlinear kinetic equation is constructed and proved to be well-adapted to describe general multidimensional scalar conservation laws. In particular, it is proved to be well-posed uniformly in epsilon - the microscopic scale. It is also shown that the proposed kinetic equation is equipped with a family of kinetic entropy functions - analogous to Boltzmann's microscopic H-function, such that they recover Krushkov-type entropy inequality on the macroscopic scale. Finally, it is proved by both - BV compactness arguments in the one-dimensional case, that the local density of kinetic particles admits a continuum limit, as it converges strongly with epsilon below 0 to the unique entropy solution of the corresponding conservation law.

Perthame, Benoit↗
NASA NTRS2001

Electron Dynamics in the Magnetotail

The goal of this research has been to study the effects of electrons on magnetotail dynamics and current sheet structure. The approach is to follow ion trajectories in a global model of the magnetotail, use a Boltzmann approximation to include electrons, and then to update the field model according to the currents that are generated by the cross-tail electric field and/or induced fields. Parallel (and perpendicular) electric fields that form are included through the Boltzmann relation. Transverse electron currents are to be included through adiabatic drift equations.

Schriver, David↗
NASA NTRS1976

Formation of the 0.511 MeV line in solar flares

The slowing down and annihilation of positrons and the formation of positronium in a solar flare plasma are investigated to determine how the width of the 0.511-MeV line and its strength relative to the three-photon continuum from positronium decay depend on the temperature and density of the medium in which the positron comes to rest. The calculations are limited to the cases of annihilation in a completely ionized plasma, in a partially ionized plasma with an electron/neutron density ratio of 1.0 or 0.1, and in an atomic gas with a very small ion density. Thermally averaged rate coefficients are obtained for the free annihilation of positrons and for positronium formation through radiative recombination in the fully ionized plasma. Positronium formation rates and the resultant energy distributions of the positronium atoms in the partially ionized medium are determined by numerically solving the Fokker-Planck equation in a medium where the ambient free electrons have a Maxwell-Boltzmann distribution of finite temperature but the density of the medium is sufficiently low that positronium atoms decay without further collisions following their formation. A Monte Carlo calculation is performed for the positron energy loss, positronium formation through charge exchange, and positronium breakup in the weakly ionized medium. The energy distributions of decaying positronium atoms and the relative number of triplet to singlet positronium decays are evaluated for ion concentrations and densities characteristic of the solar photosphere.

Crannell, C. J.↗
NASA NTRS1989

Dynamic stability of unidirectional fiber-reinforced viscoelastic composite plates

This paper deals with a dynamic stability analysis of unidirectional fiber-reinforced composite viscoelastic plates subjected to compressive edge loads. The integrodifferential equations governing the stability problem are obtained by using, in conjunction with a Boltzmann hereditary constitutive law for a three-dimensional viscoelastic medium, a higher-order shear deformation theory of orthotropic plates. Such a theory incorporates transverse shear deformation, transverse normal stress, and rotatory inertia effects. The solution of the stability problem as considered within this paper concerns the determination of the critical in-plane edge loads yielding the asymptotic instability. Numerical applications, based on material properties derived within the framework of Aboudi's micromechanical model, are presented and pertinent conclusions concerning the nature of the loss of stability and the influence of various parameters are outlined.

Chandiramani, N. K.↗
NASA NTRS1991

Spectroscopic investigations of microwave generated plasmas

The study deals with the plasma behavior as applied to spacecraft propulsion from the perspective of obtaining better design and modeling capabilities. The general theory of spectroscopy is reviewed, and existing methods for converting emission-line intensities into such quantities as temperatures and densities are outlined. Attention is focused on the single-atomic-line and two-line radiance ratio methods, atomic Boltzmann plot, and species concentration. Electronic temperatures for a helium plasma are determined as a function of pressure and a gas-flow rate using these methods, and the concentrations of ions and electrons are predicted from the Saha-Eggert equations using the sets of temperatures obtained as a function of the gas-flow rate. It is observed that the atomic Boltzmann method produces more reliable results for the electronic temperature, while the results obtained from the single-line method reflect the electron temperatures accurately.

Hawley, Martin C.↗
NASA NTRS1984

Acoustic-radiation stress in solids. I - Theory

The general case of acoustic-radiation stress associated with quasi-compressional and quasi-shear waves propagating in infinite and semiinfinite lossless solids of arbitrary crystalline symmetry is studied. The Boussinesq radiation stress is defined and found to depend directly on an acoustic nonlinearity parameter which characterizes the radiation-induced static strain, a stress-generalized nonlinearity parameter which characterizes the stress nonlinearity, and the energy density of the propagating wave. Application of the Boltzmann-Ehrenfest principle of adiabatic invariance to a self-constrained system described by the nonlinear equations of motion allows the acoustic-radiation-induced static strain to be identified with a self-constrained variation in the time-averaged product of the internal energy density and displacement gradient. The time-averaged product is scaled by the acoustic nonlinearity parameter and represents the first-order nonlinearity in the virial theorem. Finally, the relationship between the Boussinesq and the Cauchy radiation stress is obtained in a closed three-dimensional form.

Cantrell, J. H., Jr.↗
NASA NTRS1993

A nonequilibrium model for a moderate pressure hydrogen microwave discharge plasma

This document describes a simple nonequilibrium energy exchange and chemical reaction model to be used in a computational fluid dynamics calculation for a hydrogen plasma excited by microwaves. The model takes into account the exchange between the electrons and excited states of molecular and atomic hydrogen. Specifically, electron-translation, electron-vibration, translation-vibration, ionization, and dissociation are included. The model assumes three temperatures, translational/rotational, vibrational, and electron, each describing a Boltzmann distribution for its respective energy mode. The energy from the microwave source is coupled to the energy equation via a source term that depends on an effective electric field which must be calculated outside the present model. This electric field must be found by coupling the results of the fluid dynamics and kinetics solution with a solution to Maxwell's equations that includes the effects of the plasma permittivity. The solution to Maxwell's equations is not within the scope of this present paper.

Scott, Carl D.↗
NASA NTRS1992

Study of basic physical processes in liquid rocket engines

Inconsistencies between analytical results and measurements for liquid rocket thrust chamber performance, which escape suitable explanations, have motivated the examination of the basic phys ical modeling formulations as to their unlimited application. The publication of Prof. D. Straub's book, 'Thermofluid-dynamics of Optimized Rocket Propulsions,' further stimulated the interest of understanding the gas dynamic relationships in chemically reacting mixtures. A review of other concepts proposed by Falk-Ruppel (Gibbsian Thermodynamics), Straub (Alternative Theory, AT), Prigogine (Non-Equilibrium Thermodynamics), Boltzmann (Kinetic Theory), and Truesdell (Rational Mechanism) has been made to obtain a better understanding of the Navier-Stokes equation, which is now used extensively for chemically reacting flow treatment in combustion chambers. In addition to the study of the different concepts, two workshops were conducted to clarify some of the issues. The first workshop centered on Falk-Ruppel's new 'dynamics' concept, while the second one concentrated on Straub's AT. In this report brief summaries of the reviewed philosophies are presented and compared with the classical Navier-Stokes formulation in a tabular arrangement. Also the highlights of both workshops are addressed.

Wu, S. T.↗
NASA NTRS1996

A viscoelastic higher-order beam finite element

A viscoelastic internal variable constitutive theory is applied to a higher-order elastic beam theory and finite element formulation. The behavior of the viscous material in the beam is approximately modeled as a Maxwell solid. The finite element formulation requires additional sets of nodal variables for each relaxation time constant needed by the Maxwell solid. Recent developments in modeling viscoelastic material behavior with strain variables that are conjugate to the elastic strain measures are combined with advances in modeling through-the-thickness stresses and strains in thick beams. The result is a viscous thick-beam finite element that possesses superior characteristics for transient analysis since its nodal viscous forces are not linearly dependent an the nodal velocities, which is the case when damping matrices are used. Instead, the nodal viscous forces are directly dependent on the material's relaxation spectrum and the history of the nodal variables through a differential form of the constitutive law for a Maxwell solid. The thick beam quasistatic analysis is explored herein as a first step towards developing more complex viscoelastic models for thick plates and shells, and for dynamic analyses. The internal variable constitutive theory is derived directly from the Boltzmann superposition theorem. The mechanical strains and the conjugate internal strains are shown to be related through a system of first-order, ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. Equations of motion for the solid are derived from the virtual work principle using the total time-dependent stress. Numerical examples for the problems of relaxation, creep, and cyclic creep are carried out for a beam made from an orthotropic Maxwell solid.

Johnson, Arthur R.↗
NASA NTRS1992

Electron-impact vibrational relaxation in high-temperature nitrogen

Vibrational relaxation process of N2 molecules by electron-impact is examined for the future planetary entry environments. Multiple-quantum transitions from excited states to higher/lower states are considered for the electronic ground state of the nitrogen molecule N2 (X 1Sigma-g(+)). Vibrational excitation and deexcitation rate coefficients obtained by computational quantum chemistry are incorporated into the 'diffusion model' to evaluate the time variations of vibrational number densities of each energy state and total vibrational energy. Results show a non-Boltzmann distribution of number densities at the earlier stage of relaxation, which in turn suppresses the equilibrium process but affects little the time variation of total vibrational energy. An approximate rate equation and a corresponding relaxation time from the excited states, compatible with the system of flow conservation equations, are derived. The relaxation time from the excited states indicates the weak dependency of the initial vibrational temperature. The empirical curve-fit formula for the improved e-V relaxation time is obtained.

Lee, Jong-Hun↗
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 NTRS1976

The broadening of the calcium resonance line in a high-temperature helium atmosphere

The paper describes experimental measurements of the shape of the Ca I resonance line at 4227 A in a high-temperature helium atmosphere. A ballistic piston compressor was used to produce hot helium in the temperature range from 3000 to 7000 K and the number-density range from 2 to 4 by 10 to the 20th power per cu cm, which conditions approximate those postulated for the atmospheres of certain cool white-dwarf stars. The Boltzmann temperature of the calcium atoms was measured by the brightness-emissivity method, the absorption line shape was measured with a twelve-channel polychromator, and the helium kinetic temperature was calculated from the equation of an ideal gas. The observed deviation from thermodynamic equilibrium is illustrated. The value of the ratio of the damping constant to the helium number density at 5000 K is found to be 45 (+ or - 4) by 10 to the -22nd power A cu cm, or 4.7 (+ or - 0.4) by 10 to the -9th power rad/s cu cm. It is noted that no existing line-broadening calculation for the calcium-helium system is consistent with the present results. Recent studies of the 4227-A line in spectra of the sun and the white dwarf van Maanen 2 are discussed.

Driver, R. D.↗
NASA NTRS1983

Ion-acoustic solitary waves in a magnetized plasma with arbitrary electron equation of state

The oblique propagation of fully nonlinear, slow ion-acoustic solitary waves in a collisionless, low-beta, magnetized plasma is examined. The analysis includes the effects of a finite ion pressure, electron trapping, and multicomponent particle populations. The existence of both compressive and rarefactive modes propagating obliquely to the magnetic field in a plasma with two distinct Boltzmann electron populations and cold ions is demonstrated. It is shown that paired electrostatic shocks observed over the earth's auroral zone may be closely related to the rarefactive modes. As a measure of the collisionless dissipation encountered by the solitary waves, the linear response of the plasma to slow ion-acoustic waves is also examined.

Witt, E.↗
NASA NTRS1993

Ionospheric plasma flow about a system of electrically biased flat plates

The steady state interaction of two electrically biased parallel plates immersed in a flowing plasma characteristic of low earth orbit is studied numerically. Fluid equations are developed to describe the motion of the cold positively charged plasma ions, and are solved using finite-differences in two dimensions on a Cartesian grid. The behavior of the plasma electrons is assumed to be described by the Maxwell-Boltzmann distribution. Results are compared to an analytical and a particle simulation technique for a simplified flow geometry consisting of a single semi-infinite negatively biased plate. Comparison of the extent of the electrical disturbance into the flowing plasma and the magnitude of the current collected by the plate is very good. The interaction of two equally biased parallel plates is studied as a function of applied potential. The separation distance at which the current collected by either plate decreases by five and twenty percent is determined as a function of applied potential. The percent decreases were based on a non-interacting case. The decrease in overall current is caused by a decrease in ionic density in the region between the plates. As the separation between the plates decreases, the plates collect the ions at a faster rate than they are supplied to the middle region by the oncoming plasma flow. The docking of spacecraft in orbit is simulated by moving two plates of unequal potential toward one another in a quasi-static manner. One plate is held at a large negative potential while the other floats electrically in the resulting potential field. It is found that the floating plate does not charge continuously negative as it approaches the other more negatively biased plate. Instead, it charges more and then less negative as ionic current decreases and then increases respectively upon approach. When the two plates come into contact, it is expected that the electrically floating plate will charge rapidly negative to a potential near that of the other plate.

Herr, Joel L.↗
NASA NTRS1995

Nozzle flow with vibrational nonequilibrium

This research concerns the modeling and numerical solutions of the coupled system of compressible Navier-Stokes equations in cylindrical coordinates under conditions of equilibrium and nonequilibrium thermodynamics. The problem considered was the modeling of a high temperature diatomic gas N2 flowing through a converging-diverging high expansion nozzle. The problem was modeled in two ways. The first model uses a single temperature with variable specific heats as functions of this temperature. For the second model we assume that the various degrees of freedom all have a Boltzmann distribution and that there is a continuous redistribution of energy among the various degrees of freedom as the gas passes through the nozzle. Each degree of freedom is assumed to have its own temperature and, consequently, each system state can be characterized by these temperatures. This suggests that formulation of a second model with a vibrational degree of freedom along with a rotational-translation degree of freedom, each degree of freedom having its own temperature. Initially the vibrational degree of freedom is excited by heating the gas to a high temperature. As the high temperature gas passes through the nozzle throat there is a sudden drop in temperature along with a relaxation time for the vibrational degree of freedom to achieve equilibrium with the rotational-translation degree of freedom. That is, we assume that the temperature change upon passing through the throat is so great that the changes in the vibrational degree of freedom occur at a much slower pace and consequently lags behind the rotational-translational energy changes. This lag results in a finite relaxation time. In this context the term nonequilibrium is used to denote the fact that the energy content of the various degrees of freedom are characterized by two temperatures. We neglect any chemical reactions which could also add nonequilibrium effects. We develop the energy equations for the nonequilibrium model from first principles. The resulting equations, which model the nozzle flow, can be expressed in various forms. In most forms the resulting equations are coupled systems of nonlinear partial differential equations subject to certain boundary conditions. To solve the resulting coupled system of nonlinear partial differential equations, several numerical techniques were investigated: (1) the explicit MacCormack method, (2) the explicit-implicit MacCormack method, (3) the method of operator splitting, (4) factorization schemes, and (5) the Steger-Warming scheme.

Heinbockel, J. H.↗
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 NTRS1985

Application of closed-form solutions to a mesh point field in silicon solar cells

A computer simulation method is discussed that provides for equivalent simulation accuracy, but that exhibits significantly lower CPU running time per bias point compared to other techniques. This new method is applied to a mesh point field as is customary in numerical integration (NI) techniques. The assumption of a linear approximation for the dependent variable, which is typically used in the finite difference and finite element NI methods, is not required. Instead, the set of device transport equations is applied to, and the closed-form solutions obtained for, each mesh point. The mesh point field is generated so that the coefficients in the set of transport equations exhibit small changes between adjacent mesh points. Application of this method to high-efficiency silicon solar cells is described; and the method by which Auger recombination, ambipolar considerations, built-in and induced electric fields, bandgap narrowing, carrier confinement, and carrier diffusivities are treated. Bandgap narrowing has been investigated using Fermi-Dirac statistics, and these results show that bandgap narrowing is more pronounced and that it is temperature-dependent in contrast to the results based on Boltzmann statistics.

Lamorte, M. F.↗