Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic equations”

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 289 records · Page 16

Ion heating by kinetic cross-field streaming instability due to reflected ions at a quasiperpendicular shock

The present paper shows that the reflected ion at a supercritical quasi-perpendicular shock wave can excite a purely growing mode propagating parallel to the ambient magnetic field. To discuss the ion heating by such an unstable mode, the self-consistent quasi-linear kinetic equation is solved with the assumption that the present purely growing mode is the dominant unstable mode in the system. In the quasi-linear analysis of the instability, two particular cases are considered: the case of low initial ion beta and that of a high initial ion beta.

Yoon, Peter H.↗

Bounds on internal state variables in viscoplasticity

A typical viscoplastic model will introduce up to three types of internal state variables in order to properly describe transient material behavior; they are as follows: the back stress, the yield stress, and the drag strength. Different models employ different combinations of these internal variables--their selection and description of evolution being largely dependent on application and material selection. Under steady-state conditions, the internal variables cease to evolve and therefore become related to the external variables (stress and temperature) through simple functional relationships. A physically motivated hypothesis is presented that links the kinetic equation of viscoplasticity with that of creep under steady-state conditions. From this hypothesis one determines how the internal variables relate to one another at steady state, but most importantly, one obtains bounds on the magnitudes of stress and back stress, and on the yield stress and drag strength.

Freed, Alan D.↗

Stochastic acceleration of electrons. I - Effects of collisions in solar flares

Stochastic acceleration of thermal electrons to nonrelativistic energies is studied under solar flare conditions. We show that, in turbulent regions, electron-whistler wave interactions can result in the acceleration of electrons in times comparable to or shorter than the Coulomb collision time. The kinetic equation describing the evolution of the electron energy distribution including stochastic acceleration by whistlers and energy loss via Coulomb interactions is solved for an initial thermal electron energy spectrum. In general, the shape of the resulting electron distributions are characterized by the energy E(c) where systematic energy gain by turbulence equals energy loss due to Coulomb collisions. For energies less than E(c), the spectra are steep (quasi-thermal) whereas above E(c), the spectra are power laws. We find that hard X-ray spectra computed using the electron distributions obtained from our numerical simulations are able to explain the complex spectral shapes and variations observed in impulsive hard X-ray bursts. In particular, we show that the gradual steepening observed by Lin et al. (1981) could be due to a systematic increase in the density of the plasma (due to evaporation) and the increasing importance of collisions instead of the appearance of a superhot thermal component.

Hamilton, Russell J.↗

Analysis of the ionosphere-plasmasphere transport of superthermal electrons. I - Transport in the plasmasphere

Analytical solutions are developed for the kinetic equation which describes the transport of superthermal electrons in the terrestrial plasmasphere, together with a relationship which makes it possible to calculate the transparency of the plasmasphere to these electrons. In addition, analytic expressions are presented for the heating rate of the thermal plasma due to the passage of these superthermal electrons through the plasmasphere.

Khazanov, George V.↗

Numerical simulation of satellite-ring interactions - Resonances and satellite-ring torques

The Krook kinetic equation for planetary rings is numerically solved in two spatial dimensions and in time, with (1) interparticle collisions and (2) satellite-forcing, but (3) without self-gravity, for the case of a flattened planetary ring that undergoes gravitational perturbation by a nearby satellite. It is noted that the amplitude of wakes is limited by purely kinematic effects, even in the absence of collisions. Attention is given to the results of a simulation of an inner Lindblad-resonance location, as the distribution approaches steady state; these simulations do not show an increase in velocity dispersion in the resonance zone, obviating a net torque.

Brophy, Thomas G.↗

Nova dust nucleation - Kinetics and photodissociation

Dust is observed to form in nova ejecta. The grain temperature is determined by the diluted nova radiation field rather than the gas kinetic temperature, making classical nucleation theory inapplicable. We used kinetic equations to calculate the growth of carbon nuclei in these ejecta. For expected values of the parameters too many clusters grew, despite the small sticking probability of atoms to small clusters, and the clusters only reached radii of about 100 A when the carbon vapor was depleted. We then included the effects of cluster photodissociation by ultraviolet radiation from the nova. This suppresses nucleation, but too well, and no grains form at all. Finally we suggest that a few growing carbon nuclei may be protected from photodissociation by a sacrificial surface layer of hydrogen.

Johnson, D. J.↗

A kinetic description of electron beam ejection from spacecraft

A model based on the Boltzmann kinetic equation describing electron beam ejection from spacecraft is presented. Results are shown for steady-state, beam-atmosphere interaction (BAI) and for beam-plasma interaction (BPI). The BPI considered is that of elastic scattering of beam electrons with strong, short-scale Langmuir turbulence. This mechanism has been suggested for prompt echoes of beam electrons observed in sounding rocket experiments. It is shown that the interaction is sufficiently strong to account for observations of back-scattered electrons at the beam accelerator energy. However, it is clear from observations that particles undergo diffusion in energy at a rate exceeding levels expected for BAI. Inelastic scattering in BPI must therefore also be important.

Khazanov, George V.↗

Calculations of condensation and chemistry in an aircraft contrail

The flow field, chemistry, and condensation nucleation behind a transport airplane are calculated in two regimes using two separate reacting flow codes: first the axisymmetric plume, then the three dimensional vortex wake. The included chemical kinetics equations follow the evolution of the NO(y) and SO(x) chemical families. In the plume regime, the chemistry is coupled with the binary homogeneous formation of sulfate condensation nuclei, where the calculated nucleation rates predict that copious quantities of H2SO4/H2O nuclei are produced in subnanometer sizes. These sulfate aerosols could play a major role in the subsequent condensation of water vapor and the formation of contrails under favorable atmospheric conditions.

Miake-Lye, Richard C.↗

Cross-effects in microgravity flows

In microgravity materials growth (i.e., crystal growth via vapor deposition) experiments, it is of interest ot understand and control the flows that arise from the molecular rather than the mere continuum nature of the gases and the vapors. The project research is a theoretical and experimental investigation of the flow of gas/vapor mixtures in realistic geometries and environments, as well as the application of new fundamental understandings to simulating flows in the ampoules. Towards this goal, the research tasks are: to obtain a theoretical description of the flow by solving appropriate kinetics equations; to verify the results by acquiring experimental data in a diffusion cell; and to explore applications of the results above to simulations of flows in the ampoules. The paper provides a description of the cross-phenomena and the progress realized to-date by the project personnel.

Loyalka, Sudarshan K.↗

Toward a theory of interstellar turbulence. 1: Weak Alfvenic turbulence

We study weak Alfvenic turbulence of an incompressible, magnetized fluid in some detail, with a view to developing a firm theoretical basis for the dynamics of small-scale turbulence in the interstellar medium. We prove that resonant 3-wave interactions are absent. We also show that the Iroshnikov-Kraichnan theory of incompressible, magnetohydrodynamic turbulence -- which is widely accepted -- describes weak 3-wave turbulence; consequently, it is incorrect. Physical arguments, as well as detailed calculations of the coupling coefficients are used to demonstrate that these interactions are empty. We then examine resonant 4-wave interactions, and show that the resonance relations forbid energy transport to small spatial scales along the direction of the mean magnetic field, for both the shear Alfven wave and the pseudo Alfven wave. The three-dimensional inertial-range energy spectrum of 4-wave shear Alfven turbulence guessed from physical arguments reads E(k(sub z), k(sub perpendicular)) approximately V(sub A)v(sub L)L(exp -1/3)k(sub perpendicular)(exp -10/3), where V(sub A) is the Alfven speed, and v(sub L) is the velocity difference across the outer scale L. Given this spectrum, the velocity difference across lambda(sub perpendicular) approximately k(sub perpendicular exp -1) is v(sub lambda (sub perpendicular)) is approximately v(sub L)(lambda(sub perpendicular)/L)(exp 2/3). We derive a kinetic equation, and prove that this energy spectrum is a stationary solution and that it implies a positive flux of energy in k-space, along directions perpendicular to the mean magnetic field. Using this energy spectrum, we deduce that 4-wave interactions strengthen as the energy cascades to small, perpendicular spatial scales; beyond an upper bound in perpendicular wavenumber, k(sub perpendicular)L is approximately (V(sub A)/v(sub L))(exp 3/2), weak turbulence theory ceases to be valid. Energy excitation amplitudes must be very small for the 4-wave inertial-range to be substantial. When the excitation is strong, the width of the 4-wave inertial-range shrinks to zero. This seems likely to be the case in the interstellar medium.

Sridhar, S.↗

Photon-splitting limits to the hardness of emission in strongly magnetized soft gamma repeaters

Soft gamma repeaters are characterized by recurrent activity consisting of short-duration outbursts of high-energy emission that is typically of temperature less than 40 keV. One recent model of repeaters is that they originate in the environs of neutron stars with superstrong magnetic fields, perhaps greater than 10(exp 14) G. In such fields, the exotic process of magnetic photon splitting gamma yields gamma gamma acts very effectively to reprocess gamma-ray radiation down to hard X-ray energies. In this Letter, the action of photon splitting is considered in some detail, via the solution of photon kinetic equations, determining how it limits the hardness of emission in strongly magnetized repeaters, and thereby obtaining observational constraints to the field in SGR 1806-20.

Baring, Matthew G.↗

Resonant Compton scattering and gamma-ray burst continuum spectra

The Thomson limit of resonant inverse Compton scattering in the strong magnetic fields of neutron stars is considered as a mechanism for producing gamma-ray burst continuum spectra. Photon production spectra and electron cooling rates are presented using the full magnetic Thomson cross-section. Model emission spectra are obtained as self-consistent solutions of a set of photon and electron kinetic equations, displaying spectral breaks and other structure at gamma-ray energies.

Baring, M. G.↗

Non-steady-state transport of superthermal electrons in the plasmasphere

Numerical solutions to the time-dependent kinetic equation, which describes the transport of superthermal electrons in the splasmasphere between the two conjugate ionospheres, are presented. The model calculates the distribution function as a function of time, field-aligned distance, energy, and pitch-angle. The processes of refilling, depleting, and establishing steady-state conditions of superthermal electrons in the plasmasphere are discussed.

Khazanov, George V.↗

Relativistic hadrons and the origin of relativistic outflows in active galactic nuclei

We examine the hydrodynamic origin of relativistic outflows in active galactic nuclei (AGN). Specifically, we propose that the presence of a population of relativistic hadrons in the AGN 'central engine' and the associated neutron production suffices to produce outflows which under rather general conditions could be relativistic. The main such condition is that the size of the neutron production region be larger than the neutron flight path tau(sub n) approximately 3 x 10(exp 13) cm. This condition guarantees that the mean energy per particle in the proton fluid, resulting from the decay of the neutrons outside their production region, be greater than the proton rest mass. The expansion of this fluid can then lead naturally to a relativistic outflow by conversion of its internal energy to directed motion. We follow the development of such flows by solving the mass, energy as well as the kinetic equation for the proton gas in steady state, taking into account the source terms due to compute accurately the adiabatic index of the expanding gas, and in conjunction with Bernoulli's equation the detailed evolution of the bulk Lorentz factor. We further examine the role of large-scale magnetic fields in confining these outflows to produce the jets observed at larger scales.

Contopoulos, John↗

A three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawn and dusk sides of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H+ fluxes at tens of keV, which are always over-estimated. A newly-invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm-time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, j(sub o)(1+Ay(exp n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗

Three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawnside and duskside of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H(+) fluxes at tens of keV, which are always overestimated. A newly invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, J(sub o)(1 + Ay(sup n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗

A Large-Particle Monte Carlo Code for Simulating Non-Linear High-Energy Processes Near Compact Objects

High-energy radiation processes in compact cosmic objects are often expected to have a strongly non-linear behavior. Such behavior is shown, for example, by electron-positron pair cascades and the time evolution of relativistic proton distributions in dense radiation fields. Three independent techniques have been developed to simulate these non-linear problems: the kinetic equation approach; the phase-space density (PSD) Monte Carlo method; and the large-particle (LP) Monte Carlo method. In this paper, we present the latest version of the LP method and compare it with the other methods. The efficiency of the method in treating geometrically complex problems is illustrated by showing results of simulations of 1D, 2D and 3D systems. The method is shown to be powerful enough to treat non-spherical geometries, including such effects as bulk motion of the background plasma, reflection of radiation from cold matter, and anisotropic distributions of radiating particles. It can therefore be applied to simulate high-energy processes in such astrophysical systems as accretion discs with coronae, relativistic jets, pulsar magnetospheres and gamma-ray bursts.

Stern, Boris E.↗

A Burst Mode, Ultrahigh Temperature UF4 Vapor Core Reactor Rankine Cycle Space Power System Concept

Static and dynamic neutronic analyses have been performed on an innovative burst mode (100's of MW output for a few thousand seconds) Ulvahigh Temperature Vapor Core Reactor (UTVR) space nuclear power system. The NVTR employs multiple, neutronically-coupled fissioning cores and operates on a direct, closed Rankine cycle using a disk Magnetohydrodynamic (MHD) generater for energy conversion. The UTVR includes two types of fissioning core regions: (1) the central Ultrahigh Temperature Vapor Core (UTVC) which contains a vapor mixture of highly enriched UF4 fuel and a metal fluoride working fluid and (2) the UF4 boiler column cores located in the BeO moderator/reflector region. The gaseous nature of the fuel the fact that the fuel is circulating, the multiple coupled fissioning cores, and the use of a two phase fissioning fuel lead to unique static and dynamic neutronic characteristics. Static neutronic analysis was conducted using two-dimensional S sub n, transport theory calculations and three-dimensional Monte Carlo transport theory calculations. Circulating-fuel, coupled-core point reactor kinetics equations were used for analyzing the dynamic behavior of the UTVR. In addition to including reactivity feedback phenomena associated with the individual fissioning cores, the effects of core-to-core neutronic and mass flow coupling between the UTVC and the surrounding boiler cores were also included in the dynamic model The dynamic analysis of the UTVR reveals the existence of some very effectlve inherent reactivity feedback effects that are capable of quickly stabilizing this system, within a few seconds, even when large positive reactivity insertions are imposed. If the UTVC vapor fuel density feedback is suppressed, the UTVR is still inherently stable because of the boiler core liquid-fuel volume feedback; in contrast, suppression of the vapor fuel density feedback in 'conventional" gas core cavity reactors causes them to become inherently unstable. Due to the strength of the negative reactivity feedback in the UTVR, it is found that external reactivity insertions alone are inadequate for bringing about significant power level changes during normal reactor operations. Additional methods of reactivity control such as variations in the gaseous fuel mass flow rate, are needed to achieve the desired power level oontrol.

Dugan, E. T.↗