Search NASA⌕ Search

SEARCH · Search NASA

Results for “Charged particle motion”

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.

250 records · Page 14

FOSSIL—Finding Our Cosmic Roots

FOSSIL (Fragments from the Origins of the Solar System and our Interstellar Locale) is a concept to explore the largest solar system object visible to the unaided eye—the zodiacal cloud. The cloud’s many dust particles are each a tiny time capsule from a comet or asteroid. The FOSSIL concept is for the in situ compositional analysis of interplanetary dust particles, and also of the particles passing through the solar system from interstellar space. By measuring the zodiacal and interstellar particles’ velocity vectors and compositions, the approach resolves fundamental questions about the solar system’s origins. The FOSSIL concept is to: (1) Discover whether today’s local interstellar dust matches the composition of the feedstock from which the solar system formed; (2) Determine whether comets’ fine-grained component preserves unprocessed pre-solar dust or shows signs of processing in the early solar system; and (3) Learn whether comets’ and asteroids’ organic material share a common source or formed from distinct reservoirs. The FOSSIL concept is based on the use of a Dust Telescope (DT) with the capability to measure the composition and the velocity vector of dust particles to unambiguously separate interstellar from interplanetary particles, and identify a subset of zodiacal particles that are exclusively cometary. In the DT, particles pass through the trajectory sensor unharmed and impact a target plate, where they are vaporized and partially ionized. The ions are electrostatically focused onto a detector where the time-of-flight mass spectrum is recorded, enabling measurement of the composition of the interstellar solids, detection of cometary minerals altered by high temperatures or exposure to liquid water in the early solar system, and the characterization of the organic materials from the comets and asteroids that are still being delivered to the Earth today. FOSSIL’s objectives crosscut several disciplines with planetary science. Astrophysics interests lie in understanding interstellar solar matter and the only debris disk accessible in situ. Heliophysics interest is in the verification of large-scale heliospheric magnetic field models by the measured effect on the motion of electrically charged interstellar dust. Reporting the makeup of interplanetary dust that ablates in our atmosphere is valuable to Earth sciences. FOSSIL’s DT is scalable and can be accommodated to a variety of mission opportunities, without restriction on launch dates, and for a large number of possible orbits.

Turner, Neal↗

Controlling and monitoring the space-station plasma interaction: A baseline for performing plasma experiments and using advanced technology

The size, complexity, and motion of space station through the Earth's environmental plasma means that there will be a large, complicated interaction region, involving a sheath, wake, charging of surfaces, induced electric fields, secondary emission, outgassing with ionization, etc. This interaction will necessarily be a factor in carrying out and interpreting plasma experiments and in the use of certain technologies. Attention should be given ahead of time to: (1) monitoring this interaction so that it is well described; (2) implifying the interaction by appropriate design and construction of the spacecraft and by appropriate planning of technology use; and (3) controlling the interaction by both active and passive means. Plasma emitters for modifying and controlling the spacecraft charge should be placed in several locations. Portable electrostatic shields could be deployed around noisy sections of the spacecraft in order to carry out sensitive experiments. A particle umbrella could be raised to deflect the ram ions and neutrals in order to provide a controlled environment. These interactions are briefly discussed.

Whipple, Elden C.↗

The ring arcs of Neptune

After the corotation resonance with an exterior satellite proved inapplicable to the Neptune ring arc confinement, a search for other mechanisms settled on the possible influence of Neptune's magnetic field. The areas of greater optical depth around the ring are much dustier than the low optical depth regions. These particles reside in a plasma; therefore, they must carry some charge. The components of Neptune's magnetic field on the equator at the radius of the ring arcs as a function of Neptunian longitude are shown. The components are those of an offset tilted dipole model. Although the dipole model is probably not a good approximation so close to the planet, the magnitude of the field that is given is probably close to the actual value. The possible importance of the magnetic field on the smallest particles in the ring is indicated by the ratio of the magnetic field on the smallest particles in the ring is indicated by the ratio of the magnetic force to the central gravitation attraction with the field strength of B = 0.01 gauss at the ring distance. A preferred position in the orbit for magnetically perturbed particles seems to require a commensurability between the rotation of the planet and the motion of the particle in the orbit. The period of rotation is assumed to be that of the radio bursts at 16.11 hours. However, without a model for the radio emission, one cannot be absolutely sure. Jupiter's decametric radiation depends on Io's orbital position as well as the rotation, so a synodic periodicity might be appropriate. But the latter radiation is highly directed, whereas Neptune's was seen all along the spacecraft trajectory on the 16.11 hour schedule, i.e., with no shifts in phase relative to a fixed longitude on the planet. The ring orbital period is 10.536 hours which is not commensurate with the rotation period. If the 16.11 hours is interpreted as a synodic period between the rotation and a satellite motion, the closest rotation periods to 16 hours are 15.9 hours if the satellite is 1989N4 and 18.2 hours if the satellite is Triton. The former is near a 3:2 resonance with the ring particle motion. The problem deserves some more thought before a possible herding of small particles by the magnetic field is abandoned.

Peale, S. J.↗

Radial transport of storm time ring current ions

Radial transport of energetic ions for the development of the main phase of geomagnetic storms is investigated with data from the medium energy particle analyzer (MEPA) on the Charge Composition Explorer spacecraft, which monitored protons, helium ions, and the carbon-nitrogen-oxygen group, which is mostly dominated by oxygen ions. From a study of four geomagnetic storms, we show that the flux increase of these ions in the inner ring current region can be accounted for by an inward displacement of the ring current population by 0.5 to 3.5 R(E). There is a general trend that a larger inward displacement occurs at higher L shells than at lower ones. These results are in agreement with previous findings. The radially injected population consists of the prestorm population modified by substorm injections which occur on a much shorter time scale than that for a storm main phase. It is also found that the inward displacement is relatively independent of ion mass and energy, suggesting that the radial transport of these energetic ions is effected primarily by convective motion from a large electric field or by diffusion resulting from magnetic field fluctuations.

Lui, A. T. Y.↗

Airborne observations of electric fields around growing and decaying cumulus clouds

Airborne electric field data were gathered in an atmospheric electrification study near Cape Canaveral, FL. A Learjet 36A was instrumented with eight electric field meters (mills) and five different particle probes. The local electric field enhancements at each field mill site were determined under lab conditions and verified using in-flight data. The overdetermined system of eight equations (one for each field mill) was solved using a weighted least squares algorithm to compute the magnitude and direction of the ambient electric field. The signal processing system allowed the measured data to be expressed in terms of earth coordinates, regardless of the attitude of the aircraft. Thus, it was possible to take maximum advantage of the Learjet's speed and maneuverability in studying the electric field structure in the vicinity of the clouds. Data gathered while circling just outside the boundary of a growing cumulus cloud show a nonsymmetric pattern of electric field strength. Field intensity grew rapidly over a period of less than 10 minutes. The observed direction of the ambient electric field vector can be explained by an ascending motion of the charge centers of a classic tripole model of a thunderstorm.

Giori, K. L.↗

Particle acceleration very near an x-line in a collisionless plasma

In a previous paper, we applied a simplified model for particle motion in the vicinity of a magnetic X-line that had been introduced by Dungey. We used the model to quantitatively show that an electric force along an X-line can be balanced by the gyroviscous force associated with the off-diagonal elements of the pressure tensor. Distribution functions near the X-line were shown to be skewed in azimuth about the magnetic field and to include particles accelerated to very high energies. In the present paper, we apply the previous model and use the distribution functions to evaluate the energization that results from particle interactions with the X-line. We find that, in general, this interaction gives a spectrum of energized particles that can be represented by a Maxwellian distribution. A power-law, high-energy tail does not develop. The thermal energy, K, of the Maxwellian can be expressed simply in terms of the field parameters and particle mass and charge. It is independent of the thermal energy, K(sub i), of the particle distribution incident upon the region of the X-line, provided that K(sub i) is less than K. Significant energization is not found for K(sub i) is greater than K.

Lyons, L. R.↗

The striated dust tail of Comet West 1976 VI as a particle fragmentation phenomenon

The motions of 16 striae in the dust tail of Comet West are studied. It is found that all 16 striae have originated from particle ejections, and that the particles responsible for the formation of a discrete striae must be emitted simultaneously, be subjected to the same repulsive acceleration in the tail, and break up simultaneously. A strong correlation is found between the ejection times and the times of known explosive events. The repulsive accelerations of particle fragments in the striae vary from 0.6 to 2.7 times the solar attraction, indicating submicron grains of strongly absorbing materials. The repulsive acceleration of parent particles range from 0.55 to 1.10, which is only slightly smaller than those of their fragments and suggests highly nonspherical shapes of parents. The mass of dust in an average striae is estimated to be one-billion grams. Rotational bursting is discussed as a possible fragmentation mechanism. The absence of measurable effects of the Lorentz force indicate an upper limit of a few volts for the electric charge of the fragments.

Sekanina, Z.↗

Three-Dimensional Electromagnetic Monte Carlo Particle-in-Cell Simulations of Critical Ionization Velocity Experiments in Space

Although the existence of the critical ionization velocity (CIV) is known from laboratory experiments, no agreement has been reached as to whether CIV exists in the natural space environment. In this paper we move towards more realistic models of CIV and present the first fully three-dimensional, electromagnetic particle-in-cell Monte-Carlo collision (PIC-MCC) simulations of typical space-based CIV experiments. In our model, the released neutral gas is taken to be a spherical cloud traveling across a magnetized ambient plasma. Simulations are performed for neutral clouds with various sizes and densities. The effects of the cloud parameters on ionization yield, wave energy growth, electron heating, momentum coupling, and the three-dimensional structure of the newly ionized plasma are discussed. The simulations suggest that the quantitative characteristics of momentum transfers among the ion beam, neutral cloud, and plasma waves is the key indicator of whether CIV can occur in space. The missing factors in space-based CIV experiments may be the conditions necessary for a continuous enhancement of the beam ion momentum. For a typical shaped charge release experiment, favorable CIV conditions may exist only in a very narrow, intermediate spatial region some distance from the release point due to the effects of the cloud density and size. When CIV does occur, the newly ionized plasma from the cloud forms a very complex structure due to the combined forces from the geomagnetic field, the motion induced emf, and the polarization. Hence the detection of CIV also critically depends on the sensor location.

Wang, J.↗

Particle Acceleration in Active Galactic Nuclei

The high efficiency of energy generation inferred from radio observations of quasars and X-ray observations of Seyfert active galactic nuclei (AGNs) is apparently achieved only by the gravitational conversion of the rest mass energy of accreting matter onto supermassive black holes. Evidence for the acceleration of particles to high energies by a central engine is also inferred from observations of apparent superluminal motion in flat spectrum, core-dominated radio sources. This phenomenon is widely attributed to the ejection of relativistic bulk plasma from the nuclei of active galaxies, and accounts for the existence of large scale radio jets and lobes at large distances from the central regions of radio galaxies. Reports of radio jets and superluminal motion from galactic black hole candidate X-ray sources indicate that similar processes are operating in these sources. Observations of luminous, rapidly variable high-energy radiation from active galactic nuclei (AGNs) with the Compton Gamma Ray Observatory show directly that particles are accelerated to high energies in a compact environment. The mechanisms which transform the gravitational potential energy of the infalling matter into nonthermal particle energy in galactic black hole candidates and AGNs are not conclusively identified, although several have been proposed. These include direct acceleration by static electric fields (resulting from, for example, magnetic reconnection), shock acceleration, and energy extraction from the rotational energy of Kerr black holes. The dominant acceleration mechanism(s) operating in the black hole environment can only be determined, of course, by a comparison of model predictions with observations. The purpose of the work proposed for this grant was to investigate stochastic particle acceleration through resonant interactions with plasma waves that populate the magnetosphere surrounding an accreting black hole. Stochastic acceleration has been successfully applied to the problem of ion and electron energization in solar flares, and is capable of accounting for a wide range of both neutral and charged particle emissions. It is also a component in diffusive shock acceleration, since pitch-angle scattering (which is necessary for multiple shock crossings) is accompanied by diffusion in momentum space, which in turn yields a net systematic energy gain; however, stochastic energization will dominate the first-order shock process only in certain parameter regimes. Although stochastic acceleration has been applied to particle energization in the lobes of radio galaxies, its application to the central regions of AGNs has only recently been considered, but not in detail. We proposed to systematically investigate the plasma processes responsible for stochastic particle acceleration in black hole magnetospheres along with the energy-loss processes which impede particle energization. To this end we calculated acceleration rates and escape time scales for protons and electrons resonating with Alfven waves, and for electrons resonating with whistlers. Assuming either a Kolmogorov or Kraichnan wave spectrum, accretion at the Eddington limit, magnetic field strengths near equipartition, and turbulence energy densities approx. 10% of the total magnetic field energy density, we find that Alfven waves accelerate protons to Lorentz factors approx, equals 10(exp 4) - 10(exp 6) before they escape from the system. Acceleration of electrons by fast mode and whistler waves can produce a nonthermal population of relativistic electrons whose maximum energy is determined by a competition with radiation losses.

Miller, James A.↗

Interaction between Convection and Heat Transfer in Crystal Growth

Crystals are integral components in some of our most sophisticated and rapidly developing industries. Single crystals are solids with the most uniform structures that can be obtained on an atomic scale. Because of their structural uniformity, crystals can transmit acoustic and electromagnetic waves and charged particles with essentially no scattering or interferences. This transparency, which can be selectively modified by controlled additions of impurities known as dopants, is the foundation of modern electronic industry. It has brought about widespread application of crystals in transistors, lasers, microwave devices, infrared detectors, magnetic memory devices, and many other magnets and electro-optic components. The performance of a crystal depends strongly on its compositional homogeneity. For instance, in modern microcircuitry, compositional variations of a few percent (down to a submicron length scale) can seriously jeopardize predicted yields. Since crystals are grown by carefully controlled phase transformations, the compositional adjustment in the solid is often made during growth from the nutrient. Hence, a detailed understanding of mass transfer in the nutrient is essential. Moreover, since mass transfer is often the slowest process during growth, it is usually the rate limiting mechanism. Crystal growth processes are usually classified according to the nature of the parent phase. Nevertheless, whether the growth occurs by solidification from a melt (melt growth), nucleation from a solution (solution growth), condensation from a vapor (physical vapor transport) or chemical reaction of gases (chemical vapor deposition), the parent phase is a fluid. As is with most non-equilibrium processes involving fluids, liquid or vapor, fluid motion plays an important role, affecting both the concentration and temperature gradients at the soli-liquid interface.

Source record↗

O(+)-O collision frequency in high-speed flows

Throughout much of the terrestrial thermosphere and ionosphere, the motions of the neutral and ionized constituents are closely coupled and relative velocities are small, of the order of 100 m/s or less. This is particularly true at midlatitudes to low latitudes where typical velocities in the neutral gas due to tidal forcing are only 20-50 m/s. However, the solar wind-magnetosphere interaction drives a large-scale convection pattern in the polar ionosphere. When the rapid adjustment of the plasma to changes in the solar wind is combined with the slower response of the more massive neutral gas, large relative velocities on the order of 1 km/s can exist for substantial lengths of time. This will be more common during periods of high geomagnetic activity, as a result of the greater number of magnetic substorms and other particle precipitation events. When a significant relative velocity is present, the calculation of interaction parameters of the two gases passing through each other, such as collision frequency, must include that velocity. These effects are usually neglected when interpreting wind and ion drift observations. We show how the collision frequency is affected by a directed velocity between any two gases interacting with a power law or exponential potential energy curved. The directed velocity increases the collision frequency at all temperatures for most ion-neutral interactions. For certain power law potentials, such as the charge quadrupole, the collision frequency is decreased. We present an analytic solution for the high-speed collision integral using the resonance charge exchange cross section.

Pesnell, William Dean↗

A comparison of Lorentz, planetary gravitational, and satellite gravitational resonances

We consider a charged dust grain whose orbital motion is dominated by a planet's point-source gravity, but perturbed by higher-order terms in the planet's gravity field as well as by the Lorentz force arising from an asymmetric planetary magnetic field. Perturbations to Keplerian orbits due to a nonspherical gravity field are expressed in the traditional way: in terms of a disturbing function which can be expanded in a series of spherical harmonics (W. M. Kaula, 1966). In order to calculate the electromagnetic perturbation, we first write the Lorentz force in terms of the orbital elements and then substitute it into Gauss' perturbation equations. We use our result to derive strengths of Lorentz resonances and elucidate their properties. In particular, we compare Lorentz resonances to two types of gravitational resonances: those arising from periodic tugs of a satellite and those due to the attraction of an arbitrarily shaped planet. We find that Lorentz resonances share numerous properties with their gravitational counterparts and show, using simple physical arguments, that several of these patterns are fundamental, applying not only to our expansions, but to all quantities expressed in terms of orbital elements. Some of these patterns have been previously called 'd'Alembert rules' for satellite resonances. Other similarities arise because, to first-order in the perturbing force, the three problems share an integral of the motion. Yet there are also differences; for example, first-order inclination resonances exist for perturbations arising from planetary gravity and from the Lorentz force, but not for those due to an orbiting satellite. Finally, we provide a heuristic treatment of a particle's orbital evolution under the influence of drag and resonant forces. Particles brought into mean-motion resonances experience either trapping or resonant 'jumps,' depending on the direction from which the resonance is approached. We show that this behavior does not depend on the details of the perturbing force but rather is fundamental to all mean-motion resonances.

Hamilton, Douglas P.↗

Bounce-Averaged Hamiltonian for Charged Particles in an Axisymmetric but Nondipolar Model Magnetosphere

In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M, J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. The L value of a field line is (by definition) inversely proportional to the flux enclosed by the corresponding magnetic shell of equatorial radius r(sub 0), and the L value at the neutral line (r(sub 0) = b) is denoted L*. The azimuthal coordinate phi measures magnetic local time. The best functional representation found for the normalized difference (L(exp 3)a(exp 3)/mu(sub E))(B(sub m) - B(sub 0)) between mirror-point field B(sub m) and equatorial field B(sub 0) along any field line is a 5-term expansion in powers (2/3 through 6/3) of the quantity X equivalent to (La/mu(sub E))(exp 1/2)K, where K equivalent to (J(exp 2)/8m(sub 0)M)(exp 1/2) is an adiabatically conserved quantity independent of particle energy, m(sub 0) is the rest mass of the particle, and a is the radius of the Earth. This functional form is motivated by results for limiting cases in which particles mirror very near and very far from the magnetic equator. Expansion coefficients corresponding to various powers of X are obtained from least squares fits to numerically computed results for X as a function of L and B(sub m). These are accurately expressible as fourth-order polynomials in (r(sub 0)/b)(exp 3), hence indirectly as functions of L/L* = 3La/2b. This representation, which leads (except for a manageably small region of parameter space) to better than 1% accuracy in the specification of B(sub m) as a function of K and L, allows bounce-averaged guiding center simulations to be performed without actually tracing the bounce motions of individual particles. Bounce-averaged drifts L' (meridional) and phi' (azimuthal) are proportional to derivatives of the Hamiltonian H (sum of kinetic and potential energies) with respect to phi and L, respectively. Our formulation thus provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Schulz, Michael↗

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the POLAR satellite is studied by using a 3-dimensional Particle-In-Cell (PIC) code. When the satellite passes through the region of low density plasma, the satellite charges to positive potentials as high as 4050Volts due to the photoelectrons emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, an ion-rich Xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at about 2Volts. Accordingly, in our 3-dimensional PIC simulation, we considered that the potential of the satellite is 2Volts as a fixed bias. Considering the relatively high density of the Xenon plasma in the sheath (approx. 10 - 10(exp 3)/cc), the ambient plasma of low density (less than 1/cc) is neglected. In the simulations, the electric fields and plasma dynamics are calculated self-consistently. We found that an "Apple" shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission, a high positive potential hill develops. Near the Thermal Ion Detection Experiment (TIDE) detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations, it takes about a hundred electron gyroperiods for the sheath to reach a quasi-steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. Using the steady state sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Leung, Wing C.↗

The Kinematic and Microphysical Control of Storm Integrated Lightning Flash Extent

The objective of this preliminary study is to investigate the kinematic and microphysical control of lightning properties, particularly those that may govern the production of nitrogen oxides (NOx) in thunderstorms, such as flash rate, type and extent. The mixed-phase region is where the noninductive charging (NIC) process is thought to generate most storm electrification during rebounding collisions between ice particles in the presence of supercooled water. As a result, prior radar-based studies have demonstrated that lightning flash rate is well correlated to kinematic and microphysical properties in the mixed-phase region of thunderstorms such as updraft volume, graupel mass, or ice mass flux. There is also some evidence that lightning type is associated with the convective state. Intracloud (IC) lightning tends to dominate during the updraft accumulation of precipitation ice mass while cloud-to-ground (CG) lightning is more numerous during the downdraft-driven descent of radar echo associated with graupel and hail. More study is required to generalize these relationships, especially regarding lightning type, in a wide variety of storm modes and meteorological conditions. Less is known about the co-evolving relationship between storm kinematics, microphysics, morphology and three-dimensional flash extent, despite its importance for lightning NOx production. To address this conceptual gap, the NASA MSFC Lightning Nitrogen Oxides Model (LNOM) is applied to North Alabama Lightning Mapping Array (NALMA) and Vaisala National Lightning Detection NetworkTM (NLDN) observations following ordinary convective cells through their lifecycle. LNOM provides estimates of flash type, channel length distributions, lightning segment altitude distributions (SADs) and lightning NOx production profiles. For this study, LNOM is applied in a Lagrangian sense to well isolated convective cells on 3 April 2007 (single cell and multi-cell hailstorm, non-severe multicell) and 6 July 2007 (non-severe multi-cell) over Northern Alabama. The LNOM lightning characteristics are compared to the evolution of updraft and precipitation properties inferred from dual-Doppler and polarimetric radar analyses applied to observations from a nearby Doppler radar network, including the UA Huntsville Advanced Radar for Meteorological and Operational Research (ARMOR, C-band, polarimetric). The LNOM estimated SAD and lightning NOx production profiles are placed in the context of radar derived profiles of vertical motion, precipitation types and amounts. Finally, these analyses are used to determine if storm integrated flash channel extent is as well correlated to volumetric updraft and precipitation ice characteristics in the mixed phase region as flash rate for these individual convective cells.

Carey, Lawrence↗

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the Polar satellite is studied by using a three-dimensional particle-in-cell (PIC) code. When the satellite passes through the region of low-density plasma, the satellite charges to positive potentials as high as 40-50 V, owing to the photoelectron emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, ion-rich xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at approximately 2 V. Accordingly, in our three-dimensional PIC simulation we considered that the potential of the satellite is 2 V as a fixed bias. Considering the relatively high density of the xenon plasma in the sheath (10-10(exp 3)/cc), the ambient plasma of low density (<1/cc) is neglected. In the simulations the electric fields and plasma dynamics are calculated self-consistently. We found that an 'apple'-shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission a high positive potential hill develops. Near the Thermal Ion Dynamics Experiment detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations it takes only about a couple of tens of electron gyroperiods for the sheath to reach a quasi steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. After this time the expansion of the sheath in directions transverse to the ambient magnetic field slows down because the electrons are magnetized. Using the quasi steady sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Singh, N.↗