Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOLTZMANN DISTRIBUTION”

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 91 records · Page 5

Ion extraction from a plasma

The physical processes governing ion extraction from a plasma have been examined experimentally. The screen hole plasma sheath (the transition region wherein significant ion acceleration and complete electron retardation occurs) has been defined by equipotential plots for a variety of ion accelerator system geometries and operating conditions. It was found that the screen hole plasma sheath extends over a large distance, and influences ion and electron trajectories at least 15 Debye lengths within the discharge chamber. The electron density variation within the screen hole plasma sheath satisfied a Maxwell-Boltzmann density distribution at an effective electron temperature dependent on the discharge plasma primary-to-Maxwellian electron density ratio. Plasma ion flow up to and through the sheath was predominantly one-dimensional, and the ions entered the sheath region with a modified Bohm velocity. Low values of the screen grid thickness to screen hole diameter ratio were found to give good ion focusing and high extracted ion currents because of the effect of screen webbing on ion focusing.

Aston, G.↗

Theory of molecular formation by radiative association in interstellar clouds

A theory of molecular formation by radiative association is presented which is consistent with angular momentum being conserved during the encounter and which incorporates explicitly the long-range attraction between the reactants. It is pointed out that the activated complex would not have a Boltzmann energy distribution should the rotational and kinetic temperatures of the reactants differ, and it is shown how allowance for this may be made. Account is taken of the inaccessibility of a considerable fraction of the nuclear spin states of the complex. Methods are given for treating the effect which the finiteness of the dissociation frequency may have on the association rate. Calculations on some reactions of interest are performed. A very simple semiempirical formula for the rate coefficient for radiative association is also given.

Bates, D. R.↗

Boost-phase discrimination research

The final report describes the combined work of the Computational Chemistry and Aerothermodynamics branches within the Thermosciences Division at NASA Ames Research Center directed at understanding the signatures of shock-heated air. Considerable progress was made in determining accurate transition probabilities for the important band systems of NO that account for much of the emission in the ultraviolet region. Research carried out under this project showed that in order to reproduce the observed radiation from the bow shock region of missiles in their boost phase it is necessary to include the Burnett terms in the constituent equation, account for the non-Boltzmann energy distribution, correctly model the NO formation and rotational excitation process, and use accurate transition probabilities for the NO band systems. This work resulted in significant improvements in the computer code NEQAIR that models both the radiation and fluid dynamics in the shock region.

Langhoff, Stephen R.↗

Theoretical studies of a molecular beam generator

We consider the problem of modeling a high temperature diatomic gas N2 flowing through a converging-diverging high expansion nozzle. The problem is 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 population distribution which means that each degree of freedom has its own temperature and consequently each system state can be characterized by these temperatures. This suggests the formulation of a second model with a vibrational degree of freedom as having its own temperature along with a rotational-translational degree of freedom with its own temperature. Initially the vibrational degree of freedom is excited by heating the gas to a high temperature. As the high temperature gas expands through the nozzle throat there is a sudden drop in the rotational-translational temperature along with a finite relaxation time for the vibrational degree of freedom to achieve equilibrium with the rotational-translational degree of freedom. That is, the temperature change that occurs when the N2 gas passes through the nozzle throat is so great that the changes in the vibrational degree of freedom lags behind the rotational-translational energy changes. The resulting relaxation time is finite. It is in this context that the term nonequilibrium is used. That is, the term nonequilibrium denotes the fact that the energy content of the various degrees of freedom are characterized by two temperatures. We neglect any chemical reactions resulting from the high temperatures which could also add nonequilibrium effects. We develop the basic equations for the two models in various forms in order to check the derivations with other sources. The final form which is solved numerically consists of the scaled equations in a conservative dimensionless form.

Heinbockel, John H.↗

Kinetics of the Cl(2)P(J)) + CH4 Reaction: Effects of Secondary Chemistry Below 300 K

Absolute rate data for the Cl(2)P(J) + CH4 yields HCl + CH3 reaction have been obtained from 218 to 298 K by using the discharge-flow resonance fluorescence technique at I Torr total pressure. The result at 298 K is (10.1 +/- 0.6) x 10(exp -14) cu cm/molecule/s. The temperature dependence in Arrhenius form is (6.5 +/- 0.9 ) x 10(exp -12) exp[(-1235 +/- 34 )/T]. The errors given are one standard deviation; overall experimental error is estimated at +/- 15%. Because of the relatively large disagreement among earlier measurements at low temperatures, the results were examined for possible effects of non-Boltzmann spin distribution and vibrational excitation of CH4, secondary chemistry of CH3 radicals, and impurities in the CH4 source. There was no significant change in the observed rate constant when an efficient spin quencher, CF4, was added and estimates indicate that vibrational partitioning in CH4 should be at the ambient reactor temperature before the start of the reaction. The results were also independent of the source of Cl atoms (microwave discharge or thermal decomposition of Cl2) and whether CH4 was purified in-situ. However, the observed rate constant did depend on initial Cl atom concentrations and to a lesser extent on CH4 concentrations. Numerical simulations were used to assess the importance of secondary chemistry over a range of reactant concentrations

Wang, J. J.↗

Kinetics of the CI(2PJ) + CH4 Reaction: Effects of Secondary Chemistry Below 300k

Because of the relatively large disagreement among earlier measurements at low temperatures, the results were examined for possible effects of non-Boltzmann spin distribution and vibrational excitation of CH4, secondary chemistry of CH3 radicals,and impuities in the CH4 source.

atmospheric chemistry kinetics low temperature↗

Gas Flux and Density Surrounding a Cylindrical Aperture in the Free Molecular Flow Regime

The equations for rigorously calculating the particle flux and density surrounding a cylindrical aperture in the free molecular flow regime are developed and presented. The fundamental equations for particle flux and density from a reservoir and a diffusely reflecting surface will initially be developed. Assumptions will include a Maxwell-Boltzmann speed distribution, equal particle and wall temperatures, and a linear flux distribution along the cylindrical aperture walls. With this information, the equations for axial flux and density surrounding a cylindrical aperture will be developed. The cylindrical aperture will be divided into multiple volumes and regions to rigorously determine the surrounding axial flux and density, and appropriate limits of integration will be determined. The results of these equations will then be evaluated. The linear wall flux distribution assumption will be assessed. The axial flux and density surrounding a cylindrical aperture with a thickness-to-radius ratio of 1.25 will be presented. Finally, the equations determined in this study will be verified using multiple methods.

Soulas, George C.↗

Advancement of Entry System Modeling to Support Exploration of Giant Planets

This paper describes NASA’s efforts to advance entry system modeling and simulation capabilities to support future exploration of Giant planets. The Giant planets are key destinations of interest to the planetary science community for their potential to provide insight into the formation and evolution of our Solar System, as well as extrasolar planetary systems. To date, the Galileo atmospheric probe is the only purpose-built entry probe to a Giant planet. Post-flight analysis of Galileo’s performance showed that there was significant recession of the thermal protection system (TPS), well beyond what was anticipated on the flank, and this was due in part to insufficiently accurate capability for estimating the flight environment and TPS response. While Galileo ultimately survived its flight, the example serves to highlight the great challenge of designing successful missions for environments that are poorly understood or where models have not yet been validated. An important means to reduce mission risks is the incorporation of physics-based modeling with well-quantified uncertainties. The emphasis on physics-based modeling – in contrast to empirically-driven models – is motivated by the fact that it is impossible to completely replicate entry environments through ground tests and, therefore, extrapolation to the flight environment is required. Basing analysis in fundamental physics removes the bias of ground test limitations, though one must then be careful to properly characterize model inputs, simplifying assumptions, and the limits wherein the model is valid. NASA’s Entry Systems Modeling (ESM) Project is tasked with investigating such considerations for planetary science missions across the Solar System, and in recent years has begun to do so for Giant planets. The most distinctive features of the Giant planets, from an entry system perspective, are the atmospheres composed primarily of hydrogen and helium. The entry velocities of proposed missions are generally very large and can therefore be expected to result in significant convective and radiative heating generated by the vehicle’s shock layer. Yet thermochemical behavior of the hydrogen-helium system is not well understood under such conditions. The ESM project is leading efforts to develop accurate thermochemical databases based on state-of-the-art measurements in the Electric Arc Shock Tube and detailed computational chemistry. The large heat fluxes anticipated by missions has driven interest in new TPS materials, in particular woven materials, which may be enabling but have never been flown before. Consequently, multiscale models are in development to describe properties and performance of the materials from micro- to system-scale. The goal is to not only provide accurate thermal response but also to inform thermostructural reliability predictions for extreme entries. Additionally, new computational models have been developed to evaluate performance of non-destructive evaluation techniques which are vital to establishing acceptance of systems to be free of manufacturing faults like material cracking, voids, and debonding. Finally, in the area of guidance and control, aerocapture has been shown conceptually to provide a number of mission benefits, including reducing transit time and increasing payload fraction. The ESM project is building a launch-to-landing trajectory simulation capability to enable detailed studies of aerocapture maneuvers in the context of Giant planets missions. The final presentation and paper will describe each of these topics in detail, including discussion of specific gaps and the technical approach to solving them. In addition, the final paper will briefly discuss ongoing coordination between ESM project work and an ESA-funded technology development activity comprised of validation testing in the Oxford T6, IRS PWK and IST ESTHER tunnels, as well as state-to-state modeling of the shock layer to better represent non-Boltzmann energy distributions leading to non-equilibrium radiation.

Entry systems↗

An upwind, kinetic flux-vector splitting method for flows in chemical and thermal non-equilibrium

We have developed new upwind kinetic difference schemes for flows with non-equilibrium thermodynamics and chemistry. These schemes are derived from the Boltzmann equation with the resulting Euler schemes developed as moments of the discretized Boltzmann scheme with a locally Maxwellian velocity distribution. Splitting the velocity distribution at the Boltzmann level is seen to result in a flux-split Euler scheme and is called Kinetic Flux Vector Splitting (KFVS). Extensions to flows with finite-rate chemistry and vibrational relaxation is accomplished utilizing nonequilibrium kinetic theory. Computational examples are presented comparing KFVS with the schemes of Van Leer and Roe for a quasi-one-dimensional flow through a supersonic diffuser, inviscid flow through two-dimensional inlet, and viscous flow over a cone at zero angle-of-attack. Calculations are also shown for the transonic flow over a bump in a channel and the transonic flow over an NACA 0012 airfoil. The results show that even though the KFVS scheme is a Riemann solver at the kinetic level, its behavior at the Euler level is more similar to the existing flux-vector splitting algorithms than to the flux-difference splitting scheme of Roe.

Eppard, W. M.↗

The solution of the relaxation problem for the Boltzmann equation by the integral iteration method

The Boltzmann equation is considered in terms of the problem of relaxation of some initial distribution function which depends only on velocities, to Maxwell's distribution function. The Boltzmann equation is given for the relaxation problem in which the distribution function f(t, u, v) is time dependent and is also dependent on two other variables u and v (the velocities of rigid spherical molecules). An iteration process is discussed in which the velocity space u, v is subdivided into squares, the distribution function in each square being approximated by the second-order surface from the values of the distribution function at nine points. The set of all of these points forms a network of u, v values at the nodes of which the distribution function can be found.

Limar, Y. F.↗

Monte Carlo modeling of exospheric bodies - Mercury

In order to study the interaction with the surface, a Monte Carlo program is developed to determine the distribution with altitude as well as the global distribution of density at the surface in a single operation. The analysis presented shows that the appropriate source distribution should be Maxwell-Boltzmann flux if the particles in the distribution are to be treated as components of flux. Monte Carlo calculations with a Maxwell-Boltzmann flux source are compared with Mariner 10 UV spectrometer data. Results indicate that the presently operating models are not capable of fitting the observed Mercury exosphere. It is suggested that an atmosphere calculated with a barometric source distribution is suitable for more realistic future exospheric models.

Smith, G. R.↗

A spectroscopic study of CW chemical lasers.

A subsonic transverse-flow chemical laser has been instrumented to permit measurements of vibration-rotation state populations and laser output spectra under flow conditions typical of present CW hydrogen halide laser operation. Measurements are presented for the F + HI, F + HBr, F + CH4, F + HCl, F + H2, Cl + HI, and Cl + HBr atom-exchange reaction systems. CW operation was found to be characterized by relatively modest departures from Maxwell-Boltzmann vibrational energy distributions at all stages of reaction in contrast to pulsed operation based on similar reactions.

Rosen, D. I.↗

Estimation of Boltzmann damping coefficients in beam models

A distributed parameter model of a flexible structure with Boltzmann type viscoelastic damping is discussed. A computational method for the estimation of the damping parameters is developed, and theoretical convergence results are given. An example is presented in which actual experimental data is used, demonstrating the efficacy of the computational method and the plausibility of the model for predicting response in damped structures.

Banks, H. T.↗

Interaction with the lower ionosphere of electromagnetic pulses from lightning - Heating, attachment, and ionization

A Boltzmann formulation of the electron distribution function and Maxwell's equations for the EM fields are used to simulate the interaction of lightning radiated EM pulses with the lower ionosphere. Ionization and dissociative attachment induced by the heated electrons cause significant changes in the local electron density, N(e). Due to 'slow' field changes of typical lightning EM pulses over time scales of tens of microsec, the distribution function follows the quasi-equilibrium solution of the Boltzmann equation in the altitude range of interest (70 to 100 km). The EM pulse is simulated as a planar 100 microsec long single period oscillation of a 10 kHz wave injected at 70 km. Under nighttime conditions, individual pulses of intensity 10-20 V/m (normalized to 100 km horizontal distance) produce changes in N(e) of 1-30 percent while a sequence of pulses leads to strong modification of N(e) at altitudes less than 95 km. The N(e) changes produce a 'sharpening' of the lower ionospheric boundary by causing a reduction in electron density at 75-85 km (due to attachment) and a substantial increase at 85-95 km (due to ionization) (e.g., the scale height decreases by a factor of about 2 at about 85 km for a single 20 V/m EM pulse). No substantial N(e) changes occur during daytime.

Taranenko, Y. N.↗

Fundamental equations of a mixture of gas and small spherical solid particles from simple kinetic theory.

The fundamental equations of a mixture of a gas and pseudofluid of small spherical solid particles are derived from the Boltzmann equation of two-fluid theory. The distribution function of the gas molecules is defined in the same manner as in the ordinary kinetic theory of gases, but the distribution function for the solid particles is different from that of the gas molecules, because it is necessary to take into account the different size and physical properties of solid particles. In the proposed simple kinetic theory, two additional parameters are introduced: one is the radius of the spheres and the other is the instantaneous temperature of the solid particles in the distribution of the solid particles. The Boltzmann equation for each species of the mixture is formally written, and the transfer equations of these Boltzmann equations are derived and compared to the well-known fundamental equations of the mixture of a gas and small solid particles from continuum theory. The equations obtained reveal some insight into various terms in the fundamental equations. For instance, the partial pressure of the pseudofluid of solid particles is not negligible if the volume fraction of solid particles is not negligible as in the case of lunar ash flow.

Pai, S. I.↗

Kinetic theory based new upwind methods for inviscid compressible flows

Two new upwind methods called the Kinetic Numerical Method (KNM) and the Kinetic Flux Vector Splitting (KFVS) method for the solution of the Euler equations have been presented. Both of these methods can be regarded as some suitable moments of an upwind scheme for the solution of the Boltzmann equation provided the distribution function is Maxwellian. This moment-method strategy leads to a unification of the Riemann approach and the pseudo-particle approach used earlier in the development of upwind methods for the Euler equations. A very important aspect of the moment-method strategy is that the new upwind methods satisfy the entropy condition because of the Boltzmann H-Theorem and suggest a possible way of extending the Total Variation Diminishing (TVD) principle within the framework of the H-Theorem. The ability of these methods in obtaining accurate wiggle-free solution is demonstrated by applying them to two test problems.

Deshpande, S. M.↗