Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution function”

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 415 records · Page 23

A scaling law for energy transfer by inelastic electron-molecule collisions in mixtures

The equation governing the electron energy distribution in the presence of a spatially uniform electric field in a weakly ionized gas was reformulated into an integral equation for the logarithmic slope of the distribution function. For gas mixtures in which the dominant electron energy loss mechanism is by vibrational excitation of the molecules, this equation is suitable for approximate analysis and exact numerical solution by iteration. Superelastic collisions are easily included in this formulation, and do not seriously effect the convergence of the numerical scheme. The approximate analytical results are only qualitatively correct, but suggest appropriate parameters which correlate the exact numerical results very well. The distribution function as well as certain gross properties such as net energy transfer into vibration, mean energy, and drift velocity depend primarily on a single nondimensional parameter involving only E/N and the cross sections.

Bienkowski, G. K.↗

A family of models for spherical stellar systems

We describe a one-parameter family of models of stable sperical stellar systems in which the phase-space distribution function depends only on energy. The models have similar density profiles in their outer parts (rho propotional to r(exp -4)) and central power-law density cusps, rho proportional to r(exp 3-eta), 0 less than eta less than or = 3. The family contains the Jaffe (1983) and Hernquist (1990) models as special cases. We evaluate the surface brightness profile, the line-of-sight velocity dispersion profile, and the distribution function, and discuss analogs of King's core-fitting formula for determining mass-to-light ratio. We also generalize the models to a two-parameter family, in which the galaxy contains a central black hole; the second parameter is the mass of the black hole. Our models can be used to estimate the detectability of central black holes and the velocity-dispersion profiles of galaxies that contain central cusps, with or without a central black hole.

Tremaine, Scott↗

Electron Energization and Mixing Observed by MMS in the Vicinity of an Electron Diffusion Region During Magnetopause Reconnection

Measurements from the Magnetospheric Multiscale (MMS) mission are reported to show distinct features of electron energization and mixing in the diffusion region of the terrestrial magnetopause reconnection. At the ion jet and magnetic field reversals, distribution functions exhibiting signatures of accelerated meandering electrons are observed at an electron out-of-plane flow peak. The meandering signatures manifested as triangular and crescent structures are established features of the electron diffusion region (EDR). Effects of meandering electrons on the electric field normal to the reconnection layer are detected. Parallel acceleration and mixing of the inflowing electrons with exhaust electrons shape the exhaust flow pattern. In the EDR vicinity, the measured distribution functions indicate that locally, the electron energization and mixing physics is captured by two-dimensional reconnection, yet to account for the simultaneous four-point measurements, translational invariant in the third dimension must be violated on the ion-skin-depth scale.

Chen, Li-Jen↗

Wall boundary equations with slip and catalysis for multicomponent, nonequilibrium gas flows

Boundary equations obtained for a low Reynolds number, high enthalpy gas flow in regions of velocity slip and temperature jump are presented. The formulation treats, through a first-order distribution function used to include multicomponent diffusion, a multicomponent gas mixture that may be in nonequilibrium with finite-rate catalytic recombination occurring on the wall. In the boundary equations, which are obtained for use in flow-field calculations applicable to low-density flow regimes, a simplified gas/wall interaction is assumed wherein individual atoms or molecules either reflect specularly off the wall or stick and are fully accommodated. Fluxes in terms of evaluated integrals over the distribution function and integrals necessary for determining fluxes are given.

Scott, C. D.↗

Electron Bulk Acceleration and Thermalization at Earth's Quasiperpendicular Bow Shock

Electron heating at Earth's quasiperpendicular bow shock has been surmised to be due to the combined effects of a quasistatic electric potential and scattering through wave-particle interaction. Here we report the observation of electron distribution functions indicating a new electron heating process occurring at the leading edge of the shock front. Incident solar wind electrons are accelerated parallel to the magnetic field toward downstream, reaching an electron-ion relative drift speed exceeding the electron thermal speed. The bulk acceleration is associated with an electric field pulse embedded in a whistler-mode wave. The high electron-ion relative drift is relaxed primarily through a nonlinear current-driven instability. The relaxed distributions contain a beam traveling toward the shock as a remnant of the accelerated electrons. Similar distribution functions prevail throughout the shock transition layer, suggesting that the observedacceleration and thermalization is essential to the cross-shock electron heating.

Chen, L.-J.↗

Numerical investigations of low-density nozzle flow by solving the Boltzmann equation

A two-dimensional finite-difference code to solve the BGK-Boltzmann equation has been developed. The solution procedure consists of three steps: (1) transforming the BGK-Boltzmann equation into two simultaneous partial differential equations by taking moments of the distribution function with respect to the molecular velocity u(sub z), with weighting factors 1 and u(sub z)(sup 2); (2) solving the transformed equations in the physical space based on the time-marching technique and the four-stage Runge-Kutta time integration, for a given discrete-ordinate. The Roe's second-order upwind difference scheme is used to discretize the convective terms and the collision terms are treated as source terms; and (3) using the newly calculated distribution functions at each point in the physical space to calculate the macroscopic flow parameters by the modified Gaussian quadrature formula. Repeating steps 2 and 3, the time-marching procedure stops when the convergent criteria is reached. A low-density nozzle flow field has been calculated by this newly developed code. The BGK Boltzmann solution and experimental data show excellent agreement. It demonstrated that numerical solutions of the BGK-Boltzmann equation are ready to be experimentally validated.

Deng, Zheng-Tao↗

Negative-ion formation in the explosives RDX, PETN, and TNT using the Reversal Electron Attachment Detection (READ) technique

In the search for high sensitivity and direct atmospheric sampling of trace species, techniques have been developed such as atmospheric-sampling, glow-discharge ionization (ASGDI), corona discharge, atmospheric pressure ionization (API), electron-capture detection (ECD), and negative-ion chemical ionization (NICI) that are capable of detecting parts-per-billion to parts-per-trillion concentrations of trace species. These techniques are based on positive- or negative-ion formation via charge-transfer to the target, or electron capture under multiple-collision conditions in a Maxwellian distribution of electron energies at the source temperature. One drawback of the high-pressure, corona- or glow-discharge devices is that they are susceptible to interferences either through indistinguishable product masses, or through undesired ion-molecule reactions. The ASGDI technique is relatively immune from such interferences, since at target concentrations of less than 1 ppm the majority of negative ions arises via electron capture rather than through ion-molecule chemistry. A drawback of the conventional ECD, and possibly of the ASGDI, is that they exhibit vanishingly small densities of electrons with energies in the range 0-10 millielectron volts (meV), as can be seen from a typical Maxwellian electron energy distribution function at T = 300 K. Slowing the electrons to these subthermal (less than 10 meV) energies is crucial, since the cross section for attachment of several large classes of molecules is known to increase to values larger than 10(exp -12) sq cm at near-zero electron energies. In the limit of zero energy these cross sections are predicted to diverge as epsilon(exp -1/2), where epsilon is the electron energy. In order to provide a better 'match' between the electron energy distribution function and attachment cross section, a new concept of attachment in an electrostatic mirror was developed. In this scheme, electrons are brought to a momentary halt by reversing their direction with electrostatic fields. At this turning point the electrons have zero or near-zero energy. A beam of target molecules is introduced, and the resultant negative ions extracted. This basic idea has been recently improved to allow for better reversal geometry, higher electron currents, lower backgrounds, and increased negative-ion extraction efficiency. We present herein application of the so-called reversal electron attachment detector (READ) to the study of negative-ion formation in the explosives molecules RDX, PETN, and TNT under single-collision conditions.

Chutijian, Ara↗

Test particle study of minor ions in solar wind turbulence

We perform a parameter study of the temporal evolution of a test particle distribution function in MHD turbulence. The turbulent fields are calculated using a pseudo-spectral method and periodic boundary conditions on a regular grid of 180(exp 3) points, appropriate for incompressible, homogeneous and isotropic turbulence. Initially, the kinetic and the magnetic energy are equal on the average. Both, deterministic and random initial conditions are used, in the former case with zeros of the magnetic field located at grid points, in the latter case located by interpolation between grid points. The evolution of the minor ion distribution function is studied in detail as these turbulent fields evolve, developing strong current and vorticity sheets. Using the full collisionless equation of motion for the test particles, the efficiency of nonlinear interactions can be studied. The results are compared to theoretical predictions and are then discussed in connection with the observations of the dynamical properties of solar wind minor ions derived from in situ observations.

Zurbuchen, Th.↗

Prediction of crosslink density of solid propellant binders

A quantitative theory is outlined which allows calculation of crosslink density of solid propellant binders from a small number of predetermined parameters such as the binder composition, the functionality distributions of the ingredients, and the extent of the curing reaction. The parameter which is partly dependent on process conditions is the extent of reaction. The proposed theoretical model is verified by independent measurement of effective chain concentration and sol and gel fractions in simple compositions prepared from model compounds. The model is shown to correlate tensile data with composition in the case of urethane-cured polyether and certain solid propellants. A formula for the branching coefficient is provided according to which if one knows the functionality distributions of the ingredients and the corresponding equivalent weights and can measure or predict the extent of reaction, he can calculate the branching coefficient of such a system for any desired composition.

Marsh, H. E., Jr.↗

Solar wind proton velocity distributions - Comparison of the bi-Maxwellian based 16-moment expansion with observations

The purpose of this paper is to study the possible types of velocity distributions that can be obtained from the bi-Maxwellian based 16-moment expansion of the distribution function, assuming macroscopic parameter values characteristic of the range of solar wind conditions. While previous studies also took heat flow into account, the theoretical expansions for f and the definitions of the physical moments adopted in these studies were different from those used in this paper. The choice of the 16-moment expansion and corresponding moment definitions was motivated by the fact that this is the correct generalization of the widely-used Maxwellian-based 13-moment expansion to the case where the zeroth-order distribution is a bi-Maxwellian. It is found that most of the features characteristic of solar wind proton distributions can be reproduced with the 16-moment distribution, including the appearance of secondary peaks. It is also shown how each of the physically significant velocity moments affects the shape of the distribution function.

Demars, H. G.↗

Heat flux in a non-Maxwellian plasma

A hybrid numerical scheme is applied to solve the Landau equation for the electron distribution function over all velocity space. Evidence is presented for the first time of the degree and character of the failure of the classical Spitzer-Haerm heat flux approximation in a realistic solar coronal loop structure. In the loop model used, the failure is so severe at some points that the role of the heat flux in the plasma's energy balance is completely misinterpreted. In the lower corona the Spitzer-Haerm approximation predicts that the heat flux should act as an energy source, whereas the more accurate distribution functions calculated here show this to be an energy sink.

Ljepojevic, N. N.↗

Onion-shell model of cosmic ray acceleration in supernova remnants

A method is devised to approximate the spatially averaged momentum distribution function for the accelerated particles at the end of the active lifetime of a supernova remnant. The analysis is confined to the test particle approximation and adiabatic losses are oversimplified, but unsteady shock motion, evolving shock strength, and non-uniform gas flow effects on the accelerated particle spectrum are included. Monoenergetic protons are injected at the shock front. It is found that the dominant effect on the resultant accelerated particle spectrum is a changing spectral index with shock strength. High energy particles are produced in early phases, and the resultant distribution function is a slowly varying power law over several orders of magnitude, independent of the specific details of the supernova remnant.

Bogdan, T. J.↗

Feature Profile Evolution of SiO2 Trenches In Fluorocarbon Plasmas

Etching of silicon microstructures for semiconductor manufacturing in chlorine plasmas has been well characterized. The etching proceeds in a two-part process, where the chlorine neutrals passivate the Si surface and then the ions etch away SiClx. However, etching in more complicated gas mixtures and materials, such as etching of SiO2 in Ar/C4F8, requires knowledge of the ion and neutral distribution functions as a function of angle and velocity, in addition to modeling the gas surface reactions. In order to address these needs, we have developed and integrated a suite of models to simulate the etching process from the plasma reactor level to the feature profile evolution level. This arrangement allows for a better understanding, control, and prediction of the influence of equipment level process parameters on feature profile evolution. We are currently using the HPEM (Hybrid Plasma Equipment Model) and PCMCM (Plasma Chemistry Monte Carlo Model) to generate plasma properties and ion and neutral distribution functions for argon/fluorocarbon discharges in a GEC Reference Cell. These quantities are then input to the feature scale model, Simulation of Profile Evolution by Level Sets (SPELS). A surface chemistry model is used to determine the interaction of the incoming species with the substrate material and simulate the evolution of the trench profile. The impact of change of gas pressure and inductive power on the relative flux of CFx and F to the wafer, the etch and polymerization rates, and feature profiles will be examined. Comparisons to experimental profiles will also be presented.

Hwang, Helen↗

Test particle propagation in magnetostatic turbulence. 1. Failure of the diffusion approximation

The equation which governs the quasi-linear approximation to the ensemble and gyro-phase averaged one-body probability distribution function is constructed from first principles. This derived equation is subjected to a thorough investigation in order to calculate the possible limitations of the quasi-linear approximation. It is shown that the reduction of this equation to a standard diffusion equation in the Markovian limit can be accomplished through the application of the adiabatic approximation. A numerical solution of the standard diffusion equation in the Markovian limit is obtained for the narrow parallel beam injection. Comparison of the diabatic and adiabatic results explicitly demonstrates the failure of the Markovian description of the probability distribution function. Through the use of a linear time-scale extension the failure of the adiabatic approximation, which leads to the Markovian limit, is shown to be due to mixing of the relaxation and interaction time scales in the presence of the strong mean field.

Klimas, A. J.↗

Decay of equatorial ring current ions and associated aeronomical consequences

The decay of the major ion species which constitute the ring current is studied by solving the time evolution of their distribution functions during the recovery phase of a moderate geomagnetic storm. In this work, only equatorially mirroring particles are considered. Particles are assumed to move subject to E x B and gradient drifts. They also experience loses along their drift paths. Two loss mechanisms are considered: charge exchange with neutral hydrogen atoms and Coulomb collisions with thermal plasma in the plasmasphere. Thermal plasma densities are calculated with a plasmaspheric model employing a time-dependent convection electric field model. The drift-loss model successfully reproduces a number of important and observable features in the distribution function. Charge exchange is found to be the major loss mechanism for the ring current ions; however the important effects of Coulomb collisions on both the ring current and thermal populations are also presented. The model predicts the formation of a low-energy (less than 500 eV) ion population as a result of energy degradation caused by Coulomb collision of the ring current ions with the plasmaspheric electrons; this population may be one source of the low-energy ions observed during active and quiet periods in the inner magnetosphere. The energy transferred to plasmaspheric electrons through Coulomb collisions with ring current ions is believed to be the energy source for the electron temperature enhancement and the associated 6300 A (stable auroral red (SAR) arc) emission in the subauroral region. The calculated energy deposition rate is sufficient to produce a subauroral electron temperature enhancement and SAR arc emissions that are consistent with observations of these quantities during moderate magnetic activity levels.

Fok, M.-C.↗

A convected kappa distribution model for hot ions in the Jovian magnetodisc

Hot ion angular anisotropies measured by the Low Energy Charged Particle (LECP) instrument during the Voyager 2 encounter with the Jovian dayside outer magnetosphere (60-30 RJ) have been fitted to a 2 species convected kappa distribution function using a nonlinear least squares technique. The resulting parameters are well constrained by the data. The heavy ion species was assumed to be either sulfur or oxygen of unknown charge. The light species was assumed to be protons. The bulk flow speeds deduced from the model were found, contrary to some theories, to increase with increasing radial distance from Jupiter within the radial region addressed, remaining a substantial fraction (0.6) of the rigid corotation speed. Agreement with the averaged Voyager Plasma Science results was obtained near 30 RJ. The core Maxwellian temperature of the heavy ion distribution functions (30-100 keV) increased with increasing radial distance, following the trend anticipated from the corotation pickup of heavy ions. The proton temperature (20 keV) remained nearly constant.

Kane, M.↗

Electron Bernstein waves and narrow band plasma waves near the electron cyclotron frequency in the near-Sun solar wind

Context. Recent studies of the solar wind sunward of 0.25 AU reveal that it contains quiescent regions, with low-amplitude plasma and magnetic field fluctuations, and a magnetic field direction similar to the Parker spiral. The quiescent regions are thought to have a more direct magnetic connection to the solar corona than other types of solar wind, suggesting that waves or instabilities in the quiescent regions are indicative of the early evolution of the solar wind as it escapes the corona. The quiescent solar wind regions are highly unstable to the formation of plasma waves near the electron cyclotron frequency (fce). Aims. We examine high time resolution observations of these waves in an effort to understand their impact on electron distribution functions of the quiescent near-Sun solar wind. Methods. High time resolution waveform captures of near-fce waves were examined to determine variations of their amplitude and frequency in time as well as their polarization properties. Results. We demonstrate that the near-fce wave intervals contain several distinct wave types, including electron Bernstein waves and extremely narrowband waves that are highly sensitive to the ambient magnetic field orientation. Using the properties of these waves, we suggest possible plasma wave mode classifications and possible instabilities that generate these waves. The results of this analysis indicate that these waves may modify the cold core of the electron distribution functions in the quiescent near-Sun solar wind.

David M. Malaspina↗

Consequences of magnetotail ion dynamics

The trajectories of a large ensemble of particles are calculated in a modified Tsyganenko magnetic field model with a uniform cross-tail electric field. The model magnetotail can be divided into several distinct dynamical regimes of ion motion. Near Earth, where the field lines are dipolar the adiabatic formalism is adequate. In the mid-tail and distant tail, guiding-center theory breaks down and must be replaced by a quasi-adiabatic formalism. There is an important transition region between the adiabatic and quasi-adiabatic regions where ion trajectories become more complicated and no simple analytical description holds. This wall region is characterized by rapid ion acceleration and a major loss of particles to the dusk flank. The moments of the ion distribution function are constructued from the ion trajectories, including density, temperature, and pressure in the x-z and x-y planes. In the noon-midnight meridian plane, parameters are relatively constant except near the Earth, while the x-y plots show strong gradients across the magnetotail. Magnetotail plasma convects earthward, drifts toward dusk, and is squeezed out of the tail in the near-Earth region. A thin current sheet forms in the quasi-adiabatic region, and the pressure tensor has significant off-diagonal terms at its edges. These terms are the result of quasi-adiabatic ion trajectories which lead to azimuthally asymmetric distribution functions capable of maintaining approximate stress balance across the current sheet. Simplified analytical descriptions provide further physical insight into ion dynamics that are observed.

Ashour-Abdalla, Maha↗