Search NASA⌕ Search

SEARCH · Search NASA

Results for “PARTICLE TRAJECTORY”

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 217 records · Page 12

Extended atmospheres of comets and outer planet-satellite systems

For the hydrogen coma of comet P/Halley, both a Lyman-alpha image and extensive Lyman-alpha scan data obtained by the Pioneer Venus Orbiter Ultraviolet spectrometer as well as H-alpha ground-based spectral observations obtained by the University of Wisconsin Space Physics Group were successfully interpreted and analyzed with our Monte Carlo particle trajectory model. The excellent fit of the model and the Halley data and the water production rate determined near perihelion (9 Feb. 1986) from 13 Dec. 1985 to 13 Jan. 1986 and from 1 Feb. to 7 Mar. 1986 are discussed. Studies for the circumplanetary distribution of atomic hydrogen in the Saturn and Neptune systems were undertaken for escape of H atoms from Titan and Triton, respectively. The discovery of a new mechanism which can dramatically change the normal cylindrically symmetric distribution of hydrogen about the planet is discussed. The implications for the Titan-Saturn and Triton-Neptune are summarized.

Smyth, William H.↗

Laboratory study of the temporal evolution of the current-voltage characteristic of a probe in the wake of an object immersed in a pulsed flowing plasma

Measurements of the current-voltage characteristics of a Langmuir probe in the near wake of a disk immersed in a pulsed flowing plasma were made. A 1 cm diameter biasable sphere was placed in the ion-free near wake region of a 10 cm diameter disk immersed in a Mach 8 pulsed flowing plasma. The current-voltage characteristic of the sphere was observed as a function of time as the sphere bias was scanned from -5000 V to +1000 V. The collected current is found to be monotonically increasing with increasing positive bias voltage but exhibits a threshold voltage for current collection as the bias voltage becomes more negative. Potential measurements in the wake region were made for a sphere bias voltages below, at, and above the current collection threshold for a number of times during the wake formation period. The time evolution of the potential profile is shown to change as the sheath around the biased sphere is established. Predictions from the particle trajectory code SIMION are compared with data, showing excellent agreement in the prediction of the current collection threshold.

Meassick, S.↗

A two-phase restricted equilibrium model for combustion of metalized solid propellants

An Eulerian-Lagrangian two-phase approach was adopted to model the multi-phase reacting internal flow in a solid rocket with a metalized propellant. An Eulerian description was used to analyze the motion of the continuous phase which includes the gas as well as the small (micron-sized) particulates, while a Lagrangian description is used for the analysis of the discrete phase which consists of the larger particulates in the motor chamber. The particulates consist of Al and Al2O3 such that the particulate composition is 100 percent Al at injection from the propellant surface with Al2O3 fraction increasing due to combustion along the particle trajectory. An empirical model is used to compute the combustion rate for agglomerates while the continuous phase chemistry is treated using chemical equilibrium. The computer code was used to simulate the reacting flow in a solid rocket motor with an AP/HTPB/Al propellant. The computed results show the existence of an extended combustion zone in the chamber rather than a thin reaction region. The presence of the extended combustion zone results in the chamber flow field and chemical being far from isothermal (as would be predicted by a surface combustion assumption). The temperature in the chamber increases from about 2600 K at the propellant surface to about 3350 K in the core. Similarly the chemical composition and the density of the propellant gas also show spatially non-uniform distribution in the chamber. The analysis developed under the present effort provides a more sophisticated tool for solid rocket internal flow predictions than is presently available, and can be useful in studying apparent anomalies and improving the simple correlations currently in use. The code can be used in the analysis of combustion efficiency, thermal load in the internal insulation, plume radiation, etc.

Sabnis, J. S.↗

Ion radial diffusion in an electrostatic impulse model for stormtime ring current formation

Guiding-center simulations of stormtime transport of ring-current and radiation-belt ions having first adiabatic invariants mu is approximately greater than 15 MeV/G (E is approximately greater than 165 keV at L is approximately 3) are surprisingly well described (typically within a factor of approximately less than 4) by the quasilinear theory of radial diffusion. This holds even for the case of an individual model storm characterized by substorm-associated impulses in the convection electric field, provided that the actual spectrum of the electric field is incorporated in the quasilinear theory. Correction of the quasilinear diffusion coefficient D(sub LL)(sup ql) for drift-resonance broadening (so as to define D(sub LL)(sup ql)) reduced the typical discrepancy with the diffusion coefficients D(sub LL)(sup sim) deduced from guiding-center simulations of representative-particle trajectories to a factor of approximately 3. The typical discrepancy was reduced to a factor of approximately 1.4 by averaging D(sub LL)(sup sim), D(sub LL)(sup ql), and D(sub LL)(sup rb) over an ensemble of model storms characterized by different (but statistically equivalent) sets of substorm-onset times.

Chen, Margaret W.↗

Analysis of IUE observations of hydrogen in comets

The large body of hydrogen Lyman-alpha observations of cometary comae obtained with the International Ultraviolet Explorer satellite has gone generally unanalyzed because of two main modeling complications. First, the inner comae of many bright (gas productive) comets are often optically thick to solar Lyman-alpha radiation. Second, even in the case of a small comet (low gas production) the large IUE aperture is quite small as compared with the immense size of the hydrogen coma, so an accurate model which properly accounts for the spatial distribution of the coma is required to invert the inferred brightnesses to column densities and finally to H atom production rates. Our Monte Carlo particle trajectory model (MPTM), which for the first time provides the realistic full phase space distribution of H atoms throughout the coma was used as the basis for the analysis of IUE observations of the inner coma. The MCPTM includes the effects of the vectorial ejection of the H atoms upon dissociation of their parent species (H2O and OH) and of their partial collisional thermalization. Both of these effects are crucial to characterize the velocity distribution of the H atoms. A new spherical radiative transfer calculation based on our MCPTM was developed to analyze IUE observations of optically thick H comae. The models were applied to observations of comets P/Giacobini-Zinner and P/Halley.

Combi, Michael R.↗

Plasma sheath effects on ion collection by a pinhole

This work presents tables to assist in the evaluation of pinhole collection effects on spacecraft. These tables summarize results of a computer model which tracks particle trajectories through a simplified electric field in the plasma sheath. A technique is proposed to account for plasma sheath effects in the application of these results and scaling rules are proposed to apply the calculations to specific situations. This model is compared to ion current measurements obtained by another worker, and the agreement is very good.

Herr, Joel L.↗

Low-frequency mobility response functions for the central plasma sheet with application to tearing modes

Consideration is given to the effect of constant cross-tail magnetic field By on the collisionless conductivity produced by chaotic scattering and stochastic diffusion of particles in the current sheet for a parabolic geometry. It is shown that the correlation time scales as (By/Bz)-squared, and from this strong By scaling a strong tendency toward stabilization of the linear tearing modes with increasing values of By is inferred. This effect of increased dawn-dusk mobility is particularly dramatic when electrons are introduced in the calculation, and is in agreement with the results of kinetic particle simulations. The collisionless conductivity is expressed in terms of the ensemble-averaged power spectrum of the single particle trajectories, which makes it possible to calculate directly the linear conductivity instead of deriving it from the calculation of the irreversible heating rates.

Hernandez, J.↗

The OH distribution in cometary atmospheres - A collisional Monte Carlo model for heavy species

The study presents an extension of the cometary atmosphere Monte Carlo particle trajectory model formalism which makes it both physically correct for heavy species and yet computationally reasonable. The derivation accounts for the collision path and scattering redirection of a heavy radical traveling through a fluid coma with a given radial distribution in outflow speed and temperature. The revised model verifies that the earlier fast-H atom approximations used in earlier work are valid, and it is applied to a case where the heavy radical formalism is necessary: the OH distribution. It is found that a steeper variation of water production rate with heliocentric distance is required for a water coma which is consistent with the velocity-resolved observations of Comet P/Halley.

Combi, Michael R.↗

Polar wind ion dynamics in the magnetotail

The circulation of polar wind ions from the high-latitude ionosphere to the plasma sheet is investigated. Single-particle trajectory simulations for the geomagnetic tail show, in addition to the breaking of adiabaticity, the possible creation of new high-altitude mirror points. This trajectory feature results from an intense parallel deceleration imparted by the magnetic field rotation during fast ExB transport. This centrifugal deceleration yields a critical parallel velocity which depends on the magnitude of the convection electric field and below which ions remain trapped inside the neutral sheet.

Delcourt, D. C.↗

The internal dynamics of slowly rotating biological systems

The structure and the dynamics of biological systems are complex. Steady gravitational forces that act on organisms cause hydrostatic pressure gradients, stress in solid components, and ordering of movable subsystems according to density. Rotation induces internal motion; it also stresses and or deforms regions of attachment and containment. The disrupted gravitationally ordered layers of movable entities are replaced by their orbital movements. New ordering geometries may arise also, especially if fluids of various densities occur. One novel result obtained concerns the application of scheduled variation of clinostat rotation rates to the management of intracellular particle trajectories. Rotation and its consequences are discussed in terms of scaling factors for parameters such as time, derived from mathematical models for simple rotating mechanical systems.

Kessler, John O.↗

Analysis of hydrogen H-alpha observations of the coma of Comet P/Halley

The Monte Carlo Particle Trajectory Model of Combi and Smyth (1988) is used here to analyze observations of the H-alpha coma of Comet Halley. The solar excitation mechanism for the H-alpha emissions line is described. The H2O production rates derived for the H-alpha brightness measurements are shown to be very consistent with the H2O production rates determined from other Comet Halley observations of the H, O, and OH comae. Revised H2O production rates determined from 6300 A brightness measurements are presented.

Smyth, William H.↗

Observations and analysis of O(1D) and NH2 line profiles for the coma of comet P/Halley

A set of high-resolution Fabry-Perot measurements of the coma of comet P/Halley was acquired in the (O I) 6300 A and NH2 6298.62 A emission lines. These high-resolution measurements provide the first optical observations capable of studying directly the photochemical kinetics and dynamic outflow of the coma. The observations were analyzed by a Monte Carlo Particle Trajectory Model. The agreement of the model and observed line profiles was excellent and verified the underlying dynamics, exothermic photodissociative chemistry, and collisional thermalization in the coma. The somewhat wider intrinsic line profile width for the O(1D) emission in 1986 January compared to 1986 May, is, for example, produced by the larger outflow speeds and gas temperatures nearer perihelion in January. The January O(1D) profile, which is wider than the January NH2 profile, is indicative of the photochemical kinetics in the dissociation of the parent molecules H2O and OH in the coma. The absolute calibration of the observations in 1986 January allowed the production rates for H2O and the NH2-parent molecules to be determined. The average daily water production rates derived from the O(1D) emission data for January 16 and 17 are presented. These very large water production rates are consistent with the extrapolated (and 7.6 day time variable) water production rates determined from the analysis of lower spectral resolution observations for O(1D) and H-alpha emissions that covered the time period up to January 13. The large production rates on January 16 and 17 establish that the maximum water production rate for comet Halley accurred pre-perihelion in January. Implications drawn from comparison with 18 cm radio emission data in January suggest that the peak water production rate was even larger. The average production rate for NH3 determined from the NH2 emission data for January 17 was (1.48 +/- 0.10) x 10(exp 28) molecules/s, yielding an NH3/H2O production rate ratio of 0.55%.

Smyth, William H.↗

Analysis of hydrogen Lyman-alpha observations of the coma of Comet P/Halley near the perihelion

The pioneer Venus Orbiter Ultraviolet Spectrometer measurements of the Lyman-alpha intensity of atomic hydrogen excited by solar resonance scattering in the coma of Comet P/Halley acquired from December 28, 1985, to January 6, 1986, and from January 31, 1986, to March 6, 1986, are simulated with the Monte Carlo Particle Trajectory Model corrected for optical depth effects. Spatially detailed comparisons between data and model show excellent agreement and are used to infer that the highest cometary activity may not be at perihelion, but about 2 1/2 weeks before. An improved set of H2O production rates is presented for the period of time that the spectrometer was observing and is found to be consistent with the rates from other types of measurements. The apparent discrepancy between Stewart (1987) in early March and International Ultraviolet Explorer OH derived rates is resolved. The problem with the conversion of 18-cm OH radio brightness to H2O production rates is also discussed.

Smyth, William H.↗

Geomagnetic transmission of solar energetic protons during the geomagnetic disturbances of October 1989

Orbit-averaged geomagnetic transmission measurements during the large solar energetic particle events of October 1989 are presented using proton data from the NOAA-10 and GOES-7 satellies. The measurements are compared to geomagnetic transmission calculations determined by tracing particle trajectories through the combination of the International Geomagnetic Reference Field (IGRF) model and the 1989 Tsyganenko magnetospheric magnetic field model. The effective 'ring current' parameter in the 1989 Tsyganenko model based on the Dst data. Results are compared to calculations employing only the IGRF and to a parameterization of geomagnetically quiet-time cutoff rigidities derived from Cosmos/intercosmos observations. The 3-hour orbit-averaged results have approximately 15% accuracy during the October 1989 events.

Boberg, P. R.↗

Interactive data exploration and particle tracking for general circulation models

The SCENE environment for interactive visualization of complex data sets is discussed. This environment is used to create tools for graphical exploration of atmospheric flow models. These tools may be extended by the user in a seamless manner, so that no programming is required. A module for accurately tracing field lines and particle trajectories in SCENE is presented. This is used to examine the flowfield qualitatively with streamlines and pathlines and to identify critical points in the velocity field. The paper also describes a visualization tool for general circulation models on which the primary features of the environment are demonstrated.

Rosenbaum, R. I.↗

Turbulent Compressible Convection with Rotation: Mean Flows and Differential Rotation - 2

The effects of rotation on turbulent, compressible convection within stellar envelopes are studied through three-dimensional numerical simulations conducted within a local f-plane model. This work seeks to understand the types of differential rotation that can be established in convective envelopes of stars like the Sun, for which recent helioseismic observations suggest an angular velocity profile with depth and latitude at variance with many theoretical predictions. This paper analyzes the mechanisms that are responsible for the mean (horizontally averaged) zonal and meridional flows that are produced by convection influenced by Coriolis forces. The compressible convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers encompassing both laminar and turbulent flow conditions under weak and strong rotational constraints. When the nonlinearities are moderate, the effects of rotation on the resulting laminar cellular convection leads to distinctive tilts of the cell boundaries away from the vertical. These yield correlations between vertical and horizontal motions that generate Reynolds stresses that can drive mean flows, interpretable as differential rotation and meridional circulations. Under more vigorous forcing, the resulting turbulent convection involves complicated and contorted fluid particle trajectories, with few clear correlations between vertical and horizontal motions, punctuated by an evolving and intricate downflow network that can extend over much of the depth of the layer. Within such networks are some coherent structures of vortical downflow that tend to align with the rotation axis. These yield a novel turbulent alignment mechanism, distinct from the laminar tilting of cellular boundaries, that can provide the principal correlated motions and thus Reynolds stresses and subsequently mean flows. The emergence of such coherent structures that can persist amidst more random motions is a characteristic of turbulence with symmetries broken by rotation and stratification. Such structure is here found to play a crucial role in defining the mean zonal and meridional flows that coexist with the convection. Though they are subject to strong inertial oscillations, the strength and type of the mean flows are determined by a combination of the laminar tilting and the turbulent alignment mechanisms. Varying the parameters produces a wide range of mean motions. Among these, some turbulent solutions exhibit a mean zonal velocity profile that is nearly constant with depth, much as deduced by helioseismology at midlatitudes within the Sun. The solutions exhibit a definite handedness, with the direction of the persistent mean flows often prescribing a spiral with depth near the boundaries, also in accord with helioseismic deductions. The mean helicity has a profile that is positive in the upper portion of the domain and negative in the lower portion, a property bearing on magnetic dynamo processes that may be realized within such rotating layers of turbulent convection.

Brummell, Nicholas H.↗

Analysis of IUE Observations of Hydrogen in Comets

The 15-years worth of hydrogen Lyman-alpha observations of cometary comae obtained with the International Ultraviolet Explorer (IUE) satellite had gone generally unanalyzed because of two main modeling complications. First, the inner comae of many bright (gas productive) comets are often optically thick to solar Lyman-alpha radiation. Second, even in the case of a small comet (low gas production) the large IUE aperture is quite small as compared with the immense size of the hydrogen coma, so an accurate model which properly accounts for the spatial distribution of the coma is required to invert the infrared brightnesses to column densities and finally to H atom production rates. Our Monte Carlo particle trajectory model (MCPTM), which for the first time provides the realistic full phase space distribution of H atoms throughout the coma has been used as the basis for the analysis of IUE observations of the inner coma. The MCPTM includes the effects of the vectorial ejection of the H atoms upon dissociation of their parent species (H2O and OH) and of their partial collisional thermalization. Both of these effects are crucial to characterize the velocity distribution of the H atoms. This combination of the MCPTM and spherical radiative transfer code had already been shown to be successful in understanding the moderately optically thick coma of comet P/Giacobini-Zinner and the coma of comet Halley that varied from being slightly to very optically thick. Both of these comets were observed during solar minimum conditions. Solar activity affects both the photochemistry of water and the solar Lyman-alpha radiation flux. The overall plan of this program here was to concentrate on comets observed by IUE at other time during the solar cycle, most importantly during the two solar maxima of 1980 and 1990. Described herein are the work performed and the results obtained.

Combi, Michael R.↗

Three-Dimensional Simulations of Electron Beams Focused by Periodic Permanent Magnets

A fully three-dimensional (3D) model of an electron beam focused by a periodic permanent magnet (PPM) stack has been developed. First, the simulation code MAFIA was used to model a PPM stack using the magnetostatic solver. The exact geometry of the magnetic focusing structure was modeled; thus, no approximations were made regarding the off-axis fields. The fields from the static solver were loaded into the 3D particle-in-cell (PIC) solver of MAFIA where fully 3D behavior of the beam was simulated in the magnetic focusing field. The PIC solver computes the time-integration of electromagnetic fields simultaneously with the time integration of the equations of motion of charged particles that move under the influence of those fields. Fields caused by those moving charges are also taken into account; thus, effects like space charge and magnetic forces between particles are fully simulated. The electron beam is simulated by a number of macro-particles. These macro-particles represent a given charge Q amounting to that of several million electrons in order to conserve computational time and memory. Particle motion is unrestricted, so particle trajectories can cross paths and move in three dimensions under the influence of 3D electric and magnetic fields. Correspondingly, there is no limit on the initial current density distribution of the electron beam, nor its density distribution at any time during the simulation. Simulation results including beam current density, percent ripple and percent transmission will be presented, and the effects current, magnetic focusing strength and thermal velocities have on beam behavior will be demonstrated using 3D movies showing the evolution of beam characteristics in time and space. Unlike typical beam optics models, this 3D model allows simulation of asymmetric designs such as non- circularly symmetric electrostatic or magnetic focusing as well as the inclusion of input/output couplers.

Kory, Carol L.↗