Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle distributions”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Density-conserving shape factors for particle simulations in cylindrical and spherical coordinates

It is established that a particle-to-grid weighting which is based on both conventional particle-in-cell (PIC) and cloud-in-cell (CIC) shape factors will yield nonuniform grid densities, even in the case of uniform particle distributions, in both cylindrical and spherical coordinates. In view of this, alternative density-conserving weighting schemes are discussed which encompass modified PIC and CIC weighting. These alternative shape factors are mathematically no more complicated than the conventional nondensity-conserving expressions.

Ruyten, Wilhelmus M.↗
NASA NTRS1993

Shaping of the magnetotail from the mantle - Global and local structuring

Kinetic modeling of the properties of magnetotail formation from a plasma mantle source are discussed, and a unified view of the structure of the central part of the magnetotail plasma sheet as well as the structure of its boundary layer is developed. The Tsyganenko magnetic field model for quiet periods of magnetospheric activity is used to trace trajectories of mantle protons in the presence of a uniform dawn-dusk electric field. The present model shows the evolution of the global structure of the tail with a clearly defined central plasma sheet and plasma sheet boundary layer developing from its beginnings as a plasma 'nucleus' in the distant tail current sheet. This large-scale restructuring is accompanied by the creation of small-scale features in the particle distribution functions. It is also shown that complementary ring distribution structures should also exist in the central plasma sheet. The ion distribution functions in the central plasma sheet should take a rather specific form in velocity space with loss regions oriented predominantly orthogonal to the magnetic field.

Ashour-Abdalla, Maha↗
NASA NTRS1992

The evolution of large-amplitude MHD waves near quasi-parallel shocks in the solar wind

We discuss attempts to explain the steepening of magnetohydrodynamic waves upstream of the Earth's bow shock and the formation of short-wavelength 'shocklets' in terms of simple analytic models. The effort involves use of the Derivative Nonlinear Schroedinger Equation (DNLS) to describe the wave evolution. We review shortcomings of previous attempts to model foreshock wave phenomena in terms of wave packet evolution according to the DNLS. It is pointed out that: (1) the oblique propagation of the waves and (2) their growth by unstable particle distribution are conducive to the processes of wave steepening and shocklet formation. Numerical solutions of the DNLS with initial conditions corresponding to envelope-modulated, obliquely propagating wave packets yield evolved wave packets with many properties similar to the observed waves.

Spangler, S. R.↗
NASA NTRS1993

A detailed study of sheared Alfven waves and their possible application to gamma-ray bursts

Strong particle acceleration in neutron-star environments may result in the upscattering of soft photons to gamma-ray energies. Sheared Alfven waves, which may be generated by several plausible disturbances in the stellar crust, can energize the charges to Lorentz factors in excess of 10 exp 5. However, a full understanding of the physics of these MHD waves is lacking due to the complications arising from the charge separation within the sheared region. Here, we present a detailed study of this mechanism, with particular attention to a comparison of the nonrelativistic and relativistic domains. We find that the charge separation is consistent with our previous classical MHD solutions and that it provides a natural description of the sheared current via a relativistic flow. These waves can apparently produce the particle distributions that are necessary to produce the gamma-ray fluxes observed in gamma-ray bursts, regardless of whether the sources are Galactic or extragalactic.

Fatuzzo, Marco↗
NASA NTRS1993

Study of supersonic intersection flowfield at modified wing-body junctions

The problem of supersonic flow control using fillets and sweep for a wing-body junction has been investigated numerically using a three-dimensional Navier-Stokes code, which employs the MacCormack's time-split finite volume technique. An elliptic grid generation technique with direct control over spacing has been developed for constructing the grid at a filleted wing-body junction. The computed results for pressure distribution, particle paths, and limiting streamlines on the flat plate and fin surface for a swept fin show a decrease in the peak pressure on the fin leading edge and in the extent of the separated flow region. Moreover, the results for filleted juncture clearly show that the flow streamline patterns lose much of their vortical character with proper filleting. It has been demonstrated that fillets with a radius of three-and-one-half times the fin leading-edge diameter are required to weaken the vorticity in the horseshoe vortex by a factor of three for the Mach number and Reynolds number considered in the present study.

Lakshmanan, B.↗
NASA NTRS1993

Particle tracing in the magnetosphere - New algorithms and results

A fast and efficient method is employed to trace charged particles through realistic magnetospheric electric and magnetic fields, greatly reducing computer simulation times. The method is applicable for particles having arbitrary charge, energy, or pitch angle but which conserve the first two adiabatic invariants. An efficient method is used to classify drift orbits, which greatly facilitates identification of the last closed drift path or other drift boundaries. The time-independent evolution of the bounce-averaged phase space density along convective drift orbits is calculated. These three tools can be used to quantitatively describe convective evolution of the particle distribution from the tail, an essential step in understanding the production of unstable distributions in the magnetosphere.

Sheldon, R. B.↗
NASA NTRS1994

High Altitude Plasma Instrument (HAPI) data analysis

The objectives of the Dynamics Explorer mission are to investigate the coupling of energy, mass, and momentum among the earth's magnetosphere, ionosphere, and upper atmosphere. At launch, on August 3, 1981, DE-1 was placed into an elliptical polar orbit having an apogee of 23,130 km to allow global auroral imaging and crossings of auroral field lines at altitudes of several thousand kilometers. At the same time DE-2 was placed into a polar orbit, coplanar with that of DE-1 but with a perigee altitude low enough (309 km) for neutral measurements and an apogee altitude of 1012 km. The DE-1 High Altitude Plasma Instrument (HAPI) provided data on low and medium energy electrons and ions from August 13, 1981 until December 1, 1981, when a high-voltage failure occured. Analysis of HAPI data for the time period of this contract has produced new results on the source mechanisms for electron conical distributions, particle acceleration phenomena in auroral acceleration regions, Birkeland currents throughout the nightside auroral regions, the source region for auroral kilometric radiation (AKR), and plasma injection phenomena in the polar cusp.

Burch, J. L.↗
NASA NTRS1994

Spectral evolution of microwaves and hard X-rays in the 1989 March 18 flare and its interpretation

We analyze the time variation of microwave spectra and hard X-ray spectra of 1989 March 18, which are obtained from the Solar Array at the Owens Valley Radio Observatory (OVRO) and the Hard X-Ray Burst Spectrometer (HXRBS) on the Solar Maximum Mission (SMM), respectively. From this observation, it is noted that the hard X-ray spectra gradually soften over 50 - 200 keV on-and-after the maximum phase while the microwaves at 1 - 15 GHz show neither a change in spectral shape nor as rapid a decay as hard X-rays. This leads to decoupling of hard X-rays from the microwaves in the decay phase away from their good correlation seen in the initial rise phase. To interpret this observation, we adopt a view that microwave-emitting particles and hard X-ray particles are physically separated in an inhomogeneous magnetic loop, but linked via interactions with the Whistler waves generated during flares. From this viewpoint, it is argued that the observed decoupling of microwaves from hard X-rays may be due to the different ability of each source region to maintain high energy electrons in response to the Whistler waves passing through the entire loop. To demonstrate this possibility, we solve a Fokker-Planck equation that describes evolution of electrons interacting with the Whistler waves, taking into account the variation of Fokker-Planck coefficients with physical quantities of the background medium. The numerical Fokker-Planck solutions are then used to calculate microwave spectra and hard X-ray spectra for agreement with observations. Our model results are as follows: in a sronger field region, the energy loss by electron escape due to scattering by the waves is greatly enhanced resulting in steep particle distributions that reproduce the observed hard X-ray spectra. In a region with weaker fields and lower density, this loss term is reduced allowing high energy electrons to survive longer so that microwaves can be emitted there in excess of hard X-rays during the decay phase of the flare. Our results based on spectral fitting of a flare event are discussed in comparison with previous studies of microwaves and hard X-rays based on either temporal or spatial information.

Lee, Jeongwoo W.↗
NASA NTRS1994

An analysis of whistler waves at interplanetary shocks

We present an analysis of whistler wave magnetic and electric field amplitude ratios from which we compute wave propagation angles and energies of electrons in resonance with the waves. To do this analysis, we compute the theoretical dependence of ratios of wave components on the whistler wave propagation angle Theta for various combinations of orthogonal wave components. Ratios of wave components that would be observed by a spinning spacecraft are determined, and the effects of arbitrary inclinations of the spacecraft to the ambient magnetic field and to the whistler wave vector are studied. This analysis clearly demonstrates that B/E, the ratio of magnetic to electric field amplitudes, cannot be assumed to be the wave index of refraction, contrary to assumptions of some earlier studies. Therefore previous interpretations of whistler wave observations based on this assumption must be reinvestigated. B/E ratios derived using three orthogonal wave components can be used to unambiguously determine Theta. Using spin plane observations alone, a significant uncertainty occurs in the determination of Theta. Nevertheless, for whistler waves observed downstream of several interplanetary shocks by the Ulysses plasma wave experiment we find that Theta is highly oblique. We suggest that the analysis of wave amplitude ratios used in conjunction with traditional stability analyses provide a promising tool for determining which particle distributions and resonances are likely to be dominant contributors to wave growth.

Lengyel-Frey, D.↗
NASA NTRS1994

Ulysses field and plasma observations of magnetic holes in the solar wind and their relation to mirror-mode structures

The term 'magnetic hole' has been used to denote isolated intervals when the magnitude of the interplanetary magnetic field drops to a few tenths, or less, of its ambient value for a time that corresponds to a linear dimension of tens to a few hundreds of proton gyro-radii. Data obtained by the Ulysses magnetometer and solar wind anlayzer have been combined to study the properties of such magnetic holes in the solar wind between 1 AU and 5.4 AU and to 23 deg south latitude. In order to avoid confusion with decreases in field strength at interplanetary discontinuities, the study has focused on linear holes across which the field direction changed by less than 5 deg. The holes occurred preferentially, but not without exception, in the interaction regions on the leading edges of high-speed solar wind streams. Although the plasma surrounding the holes was generally stable against the mirror instability, there are indications that the holes may have been remnants of mirror-mode structures created upstream of the points of observation. Those indications include the following: (1) For the few holes for which proton of alpha-particle pressure could be measured inside the hole, the ion thermal pressure was always greater than in the plasma adjacent to the holes. (2) The plasma surrounding many of the holes was marginally stable for the mirror mode, while the plasma environment of all holes was significantly closer to mirror instability than was the average solar wind. (3) The plasma containing trains of closely spaced holes was closer to mirror instability than was the plasma containing isolated holes. (4) The near-hole plasma had much higher ion beta (ratio of thermal to magnetic pressure) than did the average solar wind. (5) Near the holes, T(sub perp)/T(sub parallel) tended to be either greater than 1 or larger than in the average wind. (6) The proton and alpha-particle distribution functions measured inside the holes occasionally exhibited the flattened phase-space-density contoures in nu(sub perp)/nu(sub parallel) space found in some numerical simulations of the mirror instability.

Winterhalter, Daniel↗
NASA NTRS1995

VLF waves in the foreshock

Plasma waves observed in the VLF range upstream of planetary bow shocks not only modify the particle distributions, but also provide important information about the acceleration processes that occur at the bow shock. Electron plasma oscillations observed near the tangent field line in the electron foreshock are generated by electrons reflected at the bow shock through a process that has been referred to as Fast Fermi acceleration. Fast Fermi acceleration is the same as shock-drift acceleration, which is one of the mechanisms by which ions are energized at the shock. We have generated maps of the VLF emissions upstream of the Venus bow shock, using these maps to infer properties of the shock energization processes. We find that the plasma oscillations extend along the field line up to a distance that appears to be controlled by the shock scale size, implying that shock curvature restricsts the flux and energy of reflected electrons. We also find that the ion acoustic waves are observed in the ion foreshock, but at Venus these emissions are not detected near the ULF forshock boundary. Through analogy with terrestrial ion observations, this implies that the ion acoustic waves are not generated by ion beams, but are instead generated by diffuse ion distributions found deep within the ion foreshock. However, since the shock is much smaller at Venus, and there is no magnetosphere, we might expect ion distributions within the ion foreshock to be different than at the Earth. Mapping studies of the terrestrial foreshock similar to those carried out at Venus appear to be necessary to determine if the inferences drawn from Venus data are applicable to other foreshocks.

Strangeway, R. J.↗
NASA NTRS1994

Kinetic and thermodynamic properties of a convecting plasma in a two-dimensional dipole field

Charged particle guiding center motion is considered in the magnetic field of a two-dimensional ('line') dipole on which is superimposed a small, static, perpendicular electric field. The parallel equation of motion is that of a simple harmonic oscillator for cos theta, the cosine of magnetic colatitude theta. Equations for the perpendicular electric and magnetic drifts are derived as well as their bounce-averaged forms. The latter are solved to yield a bounce-averaged guiding center trajectory, which is the same as that obtained from conversation of magnetic moment mu, longitudinal invariant J, and total (kinetic plus electrostatic) energy K. The algebraic simplicity of the trajectory equations is also manifest in the forms of the invariants. An interesting result is that guiding centers drift in such a way that they preserve the values of their equatorial pitch angles and (equivalently) mirror latitudes. The most general Maxwellian form of the equilibrium one-particle distribution function f is constructed from the invariants, and spatially varying density and pressure moments, parallel and perpendicular to the magnetic field, are identified. Much of the paper deals with the more restricted problem in which f is specified as a bi-Maxwellian over a straight line of finite length in the equatorial plane of the dipole and perpendicular to field lines. This might be thought of as specifying a cross-tail ion injection source; our formalism then describes the subsequent spatial development. The distribution away from the source is a scaled bi-Maxwellian but one that is cut off at large and small kinetic energies, which depend on position. Density and pressure components are reduced from the values they would have if the total content of individual flux tubes convected intact. The equatorial and meridional variations of density and pressure components are examined and compared systematically for the isotropic and highly anisotropic situations. There appears to be little qualitative difference due to anisotropy. An anisotropy measure is defined, and its spatial variation determined as a signature of possible MHD instability. Extreme values are found, larger than at the source, but the plasma beta in such regions is probably so low as to render the effect inconsequential energetically. Finally, the possible consequence of 'nonadia- batic' pressure profiles on electrostatic interchanges is considered, and a boundary delineating stabilizing and destabilizing regions determined.

Huang, T. S.↗
NASA NTRS1995

Thermal emission from particulate surfaces: A comparison of scattering models with measured spectra

Emissivity spectra of particulate mineral samples are highly dependent on particle size when that size is comparable to the wavelength of light emitted (5-50 micrometers for the midinfrared). Proper geologic interpretation of data from planetary infrared spectrometers will require that these particle size effects be well understood. To address this issue, samples of quartz powders were produced with narrow, well-characterized particle size distributions. Mean particle diameters in these samples ranged from 15 to 227 micrometers. Emission spectra of these powders allow the first detailed comparison of the complex spectral variations with particle size observed in laboratory data with the predictions of radiative transfer models. Four such models are considered here. Hapke's relectance theory (converted to emissivity via Kirchoff's law) is the first model tested. Hapke's more recently published emission theory is also employed. The third model, the 'Mie/Conel' model, uses Mie single scattering with a two-stream approximation for multiple scattering. This model, like the first, is a converted reflec- tance model. Mie scattering assumes particles are both spherical and well separated, which is not true for the quartz powders, but includes diffraction effects. The fourth model uses the Mie solution for single scattering by spheres and inputs those results into the multiple scattering formalism of Hapke's emission theory. The results of the four models are considered in relation to the values of the optical constants n and k. We have grouped these as class 1 (k large), class 2 (k moderate, n is approximately 2), class 3 (k small, n is approximately 2), and class 4 (k small, n is approximately 1). In general, the Mie/Hapke hybrid model does best at predicting variations with grain size. In particular, it predicts changes of the correct pattern, although incorrect magnitude, for class 1 bands, where large increases in emissivity with decreasing grain size are observed. This model also does an excellent job on moderate (class 2) and very weak and intraband (class 3) regions, and correctly predicts the emission maximum and its invariance with grain size near the Christiansen frequency (class 4). The Mie/Hapke hybrid model also has the fewest free parameters of the four models examined, while maintaining the most physical treatment of the radiative transfer. The Mie/Conel model performs as well as the Mie/Hapke hybrid model in strong bands (class 1) but does not accurately model the behavior of moderate (class 2) and very weak (class 3) bands.

Moersch, J. E.↗
NASA NTRS1995

Eight-moment approximation solar wind models

Heat conduction from the corona is important in the solar wind energy budget. Until now all hydrodynamic solar wind models have been using the collisionally dominated gas approximation for the heat conductive flux. Observations of the solar wind show particle distribution functions which deviate significantly from a Maxwellian, and it is clear that the solar wind plasma is far from collisionally dominated. We have developed a numerical model for the solar wind which solves the full equation for the heat conductive flux together with the conservation equations for mass, momentum, and energy. The equations are obtained by taking moments of the Boltzmann equation, using an 8-moment approximation for the distribution function. For low-density solar winds the 8-moment approximation models give results which differ significantly from the results obtained in models assuming the gas to be collisionally dominated. The two models give more or less the same results in high density solar winds.

Olsen, Espen Lyngdal↗
NASA NTRS1995

Heat flux and viscosity of ions in the collisionless solar wind

Between 1 and 2 solar radii, the Coulomb-collision mean free path for thermal ions exceeds the scale height of the solar atmosphere. The expanding solar plasma becomes collisionless and the kinetics of the solar wind are no longer dominated by thermalizing collisions. The usual Braginskii-type expressions for solar wind ion heat flux and viscosity are no longer valid. However, another microscale still exists in the solar wind, dictated by the gyro-radius of ions in the turbulent embedded solar wind magnetic field. Wave-particle interactions will act to isotropize (but not thermalize) particle distributions, and the relevant microscale for this process is the ion gyro-radius. The ion distribution can be modelled as undergoing isotropizing 'collisions,' with the relevant mean free path scaling with gyro-radius. Here, the author presents the heat flux and viscosity expected for solar wind protons which are relaxing to isotropy on a microscale that scales with gyro-radius. The collisionless viscosity and heat flux have a functional dependence different than their collisional analogs. The collisional expressions for ion viscosity and heat flux drastically overestimate the efficiency of diffusive energy and momentum transport actually operative in the solar wind.

Williams, L. L.↗
NASA NTRS1996

Studies of Interstellar Pickup Ions in the Solar Wind

The work under this grant involves studies of the interaction of interstellar pickup ions with the solar wind, with the goal of a comprehensive model of the particle distributions and wave intensities to be expected throughout the heliosphere, as well as the interactions of those distributions with the solar wind termination shock. In the past year, we have completed a number of projects, including observations and modeling of the effects of a large scattering mean free path on the pickup He(+) seen at AMPTE, an analytical model of anisotropic pickup tons in a steady radial magnetic field, and a derivation of a reduced solar wind Mach number due to increased estimates on the inflowing hydrogen density allowing for a weak termination shock. In the next year, we plan to investigate in more detail the correspondence between our models of anisotropic pickup ions and the data on spectra, variations, and proton-He(+) correlation provided by AMPTE, Ulysses, and our instrument on SOHO. We will model the time-dependent pickup ion density resulting from finite periods of radial magnetic field. We will also incorporate the effects of a large mean free path into our analysis of the He(+) focusing cone, leading to more accurate parameter values for the interstellar helium gas. This progress report also includes a discussion of our Space Physics Educational Outreach activities in the past year and plans for the next year.

Isenberg, Philip A.↗
NASA NTRS1998

A Gas-Kinetic Scheme for Reactive Flows

In this paper, the gas-kinetic BGK scheme for the compressible flow equations is extended to chemical reactive flow. The mass fraction of the unburnt gas is implemented into the gas kinetic equation by assigning a new internal degree of freedom to the particle distribution function. The new variable can be also used to describe fluid trajectory for the nonreactive flows. Due to the gas-kinetic BGK model, the current scheme basically solves the Navier-Stokes chemical reactive flow equations. Numerical tests validate the accuracy and robustness of the current kinetic method.

Lian,Youg-Sheng↗
NASA NTRS1999

Bidirectional Reflectance of Flat, Optically Thick Particulate Layers: An Efficient Radiative Transfer Solution and Applications to Snow and Soil Surfaces

We describe a simple and highly efficient and accurate radiative transfer technique for computing bidirectional reflectance of a macroscopically flat scattering layer composed of nonabsorbing or weakly absorbing, arbitrarily shaped, randomly oriented and randomly distributed particles. The layer is assumed to be homogeneous and optically semi-infinite, and the bidirectional reflection function (BRF) is found by a simple iterative solution of the Ambartsumian's nonlinear integral equation. As an exact Solution of the radiative transfer equation, the reflection function thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. Since this technique bypasses the computation of the internal radiation field, it is by far the fastest numerical approach available and can be used as an ideal input for Monte Carlo procedures calculating BRFs of scattering layers with macroscopically rough surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web at http://ww,,v.giss.nasa.gov/-crmim/brf.html and can be applied to a wide range of remote sensing, engineering, and biophysical problems. We also examine the potential effect of ice crystal shape on the bidirectional reflectance of flat snow surfaces and the applicability of the Henyey-Greenstein phase function and the 6-Eddington approximation in calculations for soil surfaces.

Mishchenko, Michael I.↗