Search

The evolution of the surfaces of bodies unprotected by either strong magnetic fields or thick atmospheres in the Solar System is caused by various processes, induced by photons, energetic ions and micrometeoroids. Among these processes, the continuous bombardment of the solar wind or energetic magnetospheric ions onto the bodies may significantly affect their surfaces, with implications for their evolution. Ion precipitation produces neutral atom releases into the exosphere through ion sputtering, with velocity distribution extending well above the particle escape limits. We refer to this component of the surface ejecta as sputtered high-energy atoms (SHEA). The use of ion sputtering emission for studying the interaction of exposed bodies (EB) with ion environments is described here. Remote sensing in SHEA in the vicinity of EB can provide mapping of the bodies exposed to ion sputtering action with temporal and mass resolution. This paper speculates on the possibility of performing remote sensing of exposed bodies using SHEA The evolution of the surfaces of bodies unprotected by either strong magnetic fields or thick atmospheres in the Solar System is caused by various processes, induced by photons, energetic ions and micrometeoroids. Among these processes, the continuous bombardment of the solar wind or energetic magnetospheric ions onto the bodies may significantly affect their surfaces, with implications for their evolution. Ion precipitation produces neutral atom releases into the exosphere through ion sputtering, with velocity distribution extending well above the particle escape limits. We refer to this component of the surface ejecta as sputtered high-energy atoms (SHEA). The use of ion sputtering emission for studying the interaction of exposed bodies (EB) with ion environments is described here. Remote sensing in SHEA in the vicinity of EB can provide mapping of the bodies exposed to ion sputtering action with temporal and mass resolution. This paper speculates on the possibility of performing remote sensing of exposed bodies using SHEA and suggests the need for quantitative results from laboratory simulations and molecular physic modeling in order to understand SHEA data from planetary missions. In the Appendix, referenced computer simulations using existing sputtering data are reviewed.

The need for a thin foil to serve both as an extreme ultraviolet (EUV) radiation attenuator and as a secondary electron generator for time-of-flight (TOF) instruments and energetic neutral atom (ENA) imagers has led us to test the scattering and transmission characteristics of a layered foil of nominal aerial thicknesses of 5 microns/sq cm Si, 4 microns/sq.cm Lexan ((C16)H(14)O3), and 1 micron/sq.cm C. We report the scattering and transmission results of subjecting this foil to H(sup +), He(sup +), and O(sup +) beams of energy less than 120 keV. The scattering of these particles, defined by the half-angular-width at half-maximum (HWHM), decreases with increasing energy and is found to be 10 deg for 10 keV H(sup +), 22 keV for He(sup +), and 100 keV for O(sup +). The 50% transmission threshold for these particles was found to be 3.9 keV for H(sup +), 8.6 keV for He(sup +), and 49.5 keV for O(sup +). Additionally, ion transmission shows evidence of resonance scattering from a rectangular potential independent of incident ion mass. Particularly evident is a local transmission minimum at 70 keV possessing a Gaussian-width of 20 keV. For incident energies below 50 keV, computer simulations of H+-foil interactions using SRIM96 agree well with the H(sup +) data, while yielding less scattering and higher transmission than indicated by the O(sup +) data.

The need for a thin foil to serve both as an extreme ultraviolet (EUV) radiation attenuator and as a secondary electron generator for time-of-flight (TOF) instruments and energetic neutral atom (ENA) imagers has led us to test the scattering and transmission characteristics of a layered foil of nominal aerial thicknesses of 5 micrograms (raised dot) per square centimeter Si, 4 micrograms (raised dot) per square centimeter Lexan (C16H14O3), and 1 microgram (raised dot) per square centimeter C. We report the scattering/transmission results of subjecting this foil to H(+), He(+), and O(+) beams of energies less than or equal to 120 keV. The scattering of these particles, defined by the half-angular-width at half-maximum (HWHM), decreases with increasing energy and found to be 10 degrees for 10 keV H(+), 22 keV for He(+), and 100 keV for O(+). The 50% transmission threshold for these particles was found to be 3.9 keV for H(+), 8.6 keV for He(+), and 49.5 keV for O(+). Computer simulations of H(+)-foil interactions using SRIM96 agree well with the H(+) data, while yielding less scattering and higher transmission than indicated by the O(+) data.

The linearized governing equations for an ideal fluid are presented for numerical analysis for the stability of free capillary surfaces in rectangular containers against unfavorable disturbances (accelerations,i.e. Rayleigh-Taylor instability). The equations are solved for the case of the right square cylinder. The results are expressed graphically in term of a critical Bond number as a function of system contact angle. A critical wetting phenomena in the corners is shown to significantly alter the region of stability for such containers in contrast to simpler geometries such as the right circular cylinder or the infinite rectangular slot. Such computational results provide additional constraints for the design of fluids systems for space-based applications.

Using solar wind particle data from the charge, element and isotope analysis system (CELIAS) experiment on the SOHO mission, densities of the elements O, Ne, Mg, Si, S, Ca, and Fe are derived, and their abundance is analyzed before, during and after the 6 Janaury 1997 coronal mass ejection event (CME). In the interstream and coronal hole regions before and after this event, typical solar wind abundances for the elements investigated were found. However, during the passage of the coronal mass ejection and during the passage of the erupted filament, the elemental composition differed markedly from typical solar wind. For the passage of the CME and for the passage of the erupted filament, a mass-dependent enhancement of the elements was found, with a monotonic increase towards heavier elements. Si/O and Fe/O ratios of the order of one during these time periods were observed.

The use of energetic neutral atoms (ENAs) to study the outer heliosphere is discussed. The detection of ENAs under stringent observational conditions is addressed. Imaging ENA instruments are examined.

An overview of recent advances in theoretical analyses of supercritical droplet vaporization and combustion is conducted. Both hydrocarbon and cryogenic liquid droplets over a wide range of thermodynamic states are considered. Various important high-pressure effects on droplet behavior, such as thermodynamic non-ideality, transport anomaly, and property variation, are reviewed. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influence on fluid transport, gas-liquid interfacial thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibit a considerable change in the combustion mechanism at the criticl pressure, mainly as a result of reduced mass diffusivity and latent heat of vaporization with increased pressure. The influence of droplet size on the burning characteristics is also noted.

The idea of using energetic neutral atoms (ENAs), produced by charge exchange between energetic ions and ambient neutral atoms, as a diagnostic tool to investigate planetary magnetospheres from a distance has been extended to the investigation of the heliosphere. The paper explores what one can reasonably expect of the heliospheric ENA (HSENA) and what criteria would be imposed on HSENA instruments by concentrating on 10-10 exp 3 keV protons in quiet-time interplanetary space, solar-flare events, corotating interaction regions, and populations have distinctive signatures and that the detection of these particles can reveal energy spatial and propagation of ions in 3D interplanetary space, including the solar-wind termination shock. Such breadth of information could not be gained by in situ means.

An investigation of the Runge-Kutta time-stepping, combined with compact difference schemes to solve the unsteady Euler equations, is presented. Initially, a generalized form of a N-step Runge-Kutta technique is derived. By comparing this generalized form with its Taylor's series counterpart, the criteria for the three-step and four-step schemes to be of third- and fourth-order accurate are obtained.

The solar-wind termination shock is inaccessible for repeated in situ investigation. We examine, therefore, the possibility of remote sensing the entire heliopause from Earth's orbit using the energetic neutral atoms (ENA) produced by charge exchange between energetic ions and the neutral atoms of the interstellar medium at and beyond the termination shock. We estimate the ENA fluxes at Earth's orbit coming from the thermalized solar-wind ions and the shock-accelerated anomalous cosmic rays (ACR) at the heliospheric boundary.

Studies have shown that Earth passages of fast coronal mass ejections (CMEs) trigger geomagnetic storms. Early identification of fast earth-directed CME can help provide storm warnings, but detection of such by coronagraphs is extremely difficult. We suggest that energetic hydrogen atoms (EHA) between 2 and 10 keV produced during the transit phase of an Earth-directed CME by recombination between protons and electrons in the CME can travel ahead of the CME and act as harbingers of a magnetic storm. This forecasting scheme should work if enough EHA are produced, because while CMEs decelerate continuously after their ejection, the EHA fluxes produced in the initial phase of fast CMEs propagate at their initial high speeds. Model simulations support this proposed mechanism.

A survey is conducted for state-of-the-art techniques for detecting energetic neutral atoms (ENAs) in the 100-300 keV range, in regions from the heliospheric boundary to the auroral zones where the solar wind plays a crucial role. While ENA spectrometry allows sampling of the mass and energy distributions of a distant plasma, ENA imaging gives a global view of the structures and dynamics of an extended plasma. The ENA instrument designs discussed share many components which exhibit excellent flight performance as elements in charged-particle analyzers for space missions.

Coded aperture sensors for photons or energetic neutral atoms (ENAs), which incorporate FOV limiters and subdivide the object field into a number of elements which is smaller than the number of detector pixels, are described. A least squares fit to the data is made in reconstructing the object field. To evaluate the optics and reconstruction algorithms, two 'breadboard' sensors have been constructed, one based on a film camera and the other upon a UV-light sensitive microchannel plate detector system. Results obtained thus far show that the concept is viable, and no special difficulties should be encountered in adapting the detector geometries to neutral particle imaging systems. Charged particle deflection plates could be incorporated into the region between the FOV limiter and the aperture, or installed ahead of the limiter.

A class of instruments designed for remote sensing of space plasmas by measuring energetic neutral atoms (ENA) uses a thin foil as both a signal generator and a light shield. An ENA imager must look directly at the ENA source region, which is also usually an intense source of H Ly-alpha (1216 A) photons. It is desirable to minimize the energy threshold for ENA detectors, at the same time maximizing the blocking of H Ly-alpha. Optimizing filter design to meet these two contrary requirements has led us to measure the transmittance of thin C, Si/C, and Al/C foils at H Ly-alpha. Our results indicate that (1) transmittance of less than 0.0007 can be achieved with 7 micro-g/sq cm Si on 1.7 micro-g/sq cm C; (2) an Si/C composite foil with a thin carbon layer is more effective in blocking UV radiation while having the lowest energy threshold of all the foils measured; and (3) transmittance of Si/C foils of known Si and C thicknesses cannot be accurately predicted, but must be measured.

The spatial distribution and expansion velocity of the Comets Wilson (1987 VII) and pre- and postperihelion P/Halley are derived on the bases of velocity-resolved H2O spectral line profiles, using a kinematic model which synthesizes line profiles for comparison with observed line shapes. The results thus obtained demonstrate that the spherically symmetric outflow at constant velocity is a poor characterization of cometary neutral-gas outflow. While the radial dependence of the H2O expansion velocity is noted to be consistent with theoretically envisioned trends, the high H2O outflow velocity observed in Comet Wilson resists reconciliation with any existing kinematic model.

The present overview of cometary spectral-line profiles obtainable through advancements in high-resolution spectroscopic studies, which allow the retrieval of coma kinematic properties from velocity-resolved spectral-line profiles, incorporates the most important gas dynamic processes into an outflow model which is tailored to the interpretation of spectroscopic observations of parent molecules. The model is then used to study the influence on parent-molecule spectral line profile formation of the field-of-view, the expansion velocity, the kinetic temperature, and the anisotropic outflow distributions.

The systems of equations governing chemically nonequilibrium flows are solved numerically by two lower-upper(LU) algorithms, namely, LU-SSOR (Yoon and Jameson, 1987) and LU-SW (Steger and Warming, 1981). The methods of flux-Jacobian splitting differentiate the two LU algorithms. The LU-SSOR gives a simpler implicit operator with less temporal damping, while the LU-SW converges faster but requires more computation for each iteration. Relative efficiencies and mutual merits of these two LU algorithms are compared. The results show that LU-SW with upwind-differenced right-hand-side gives the best convergence.

Compact electrostatic analyzers that employ a new geometry are described. They accept collimated beams of ions through relatively wide (non-slit-like) inlet apertures and focus them at various locations, depending on the ion energy per charge, along rectangular position-sensing detectors. Each analyzer has a dynamic range of 10 or more for a fixed deflection voltage. The focal line width at a given energy is a small fraction of the inlet aperture width, and the focal lines are approximately coplanar.

A model for the production of energetic neutral atoms (ENAs) in the Saturn/Titan system, that takes into account existing data on the spatial distribution of the background neutral density, the energy and spatial distributions of the energetic ions, and the energy dependence of the charge-exchange cross sections, is presented. Hydrogen production via H(+), H is presented in detail. The extended Saturn/Titan system is 'imaged' in energetic neutral H at 66 R(s) and 25-deg latitude. From this vantage point, the energetic H flux is 6 + or - 1 H/sq cm per sec between 10 and 100 keV. Comparing the predicted H spectrum with that observed by Kirsch et al. (1981), the production of neutral O is estimated, and it is concluded that O(+), H charge exchange in the inner magnetosphere is sufficient to account for all the ENA O production.

This research is motivated by the improved ballistic performance of large-caliber guns using stick propellant charges. A comprehensive theoretical model for predicting the flame spreading, combustion, and grain deformation phenomena of long, unslotted stick propellants is presented. The formulation is based upon a combined Eulerian-Lagrangian approach to simulate special characteristics of the two phase combustion process in a cartridge loaded with a bundle of sticks. The model considers five separate regions consisting of the internal perforation, the solid phase, the external interstitial gas phase, and two lumped parameter regions at either end of the stick bundle. For the external gas phase region, a set of transient one-dimensional fluid-dynamic equations using the Eulerian approach is obtained; governing equations for the stick propellants are formulated using the Lagrangian approach. The motion of a representative stick is derived by considering the forces acting on the entire propellant stick. The instantaneous temperature and stress fields in the stick propellant are modeled by considering the transient axisymmetric heat conduction equation and dynamic structural analysis.

To unravel the liquid-jet breakup process in the nondilute region, a newly developed system of real-time X-ray radiography, an advanced digital image processor, and a high-speed video camera were used. Based upon recorded X-ray images, the inner structure of a liquid jet during breakup was observed. The jet divergence angle, jet breakup length, and fraction distributions along the axial and transverse directions of the liquid jets were determined in the near-injector region. Both wall- and free-jet tests were conducted to study the effect of wall friction on the jet breakup process.

Various computer programs based upon the SIMPLE or SIMPLER algorithm were studied and compared for numerical accuracy, efficiency, and grid dependency. Four two-dimensional and one three-dimensional code originally developed by a number of research groups were considered. In general, the accuracy and computational efficieny of these TEACH type programs were improved by modifying the differencing schemes and their solvers. A brief description of each program is given. Error reduction, spline flux and second upwind differencing programs are covered.

The neutral gas density profile of comet Halley measured by the Neutral Gas Experiment on Vega 1 showed an asymmetry between the inbound and the outbound legs during the fly-by on 6 March 1986. The implications of this asymmetry are discussed, and it is shown how the asymmetry detected by NGE on Vega 1 can be traced back to regions on or near the nucleus to obtain their relative gas production activities at specific times of emission.

Data from the Vega 1 permitted the determination of the total neutral gas density profile along the spacecraft trajectory. Discounting small fluctuations, the field ionization source instrument measured a density profile which varied approximately as the inverse radial distance squared. Data from the electron impact ionization instrument yielded a series of calibration points; e.g., the neutral density at 100,000 km is 10,000/cc. The combined data provide a calibrated total density profile, and imply a neutral production rate of 10 to the 30th power molecules/sec.

An ion source that utilizes a pair of thin-film field emission cathodes to provide current for electric impact ionization is studied. The field emission cathodes are composed of microscopically small needle tips and an extraction electrode; the cathodes generate a 200-eV, 0.1-mA electron beam. The multistep process for the manufacturing of the field emission cathodes is examined. The operation of the ion source, which ionizes atoms and molecules that enter the mass spectrometer with a ram velocity of about 80 km/sec, is described.

The formulation and numerical analysis of several problems related to the behavior of pneumatic tires are considered. These problems include the general rolling contact problem of a rubber-like viscoelastic cylinder undergoing finite deformations and the finite deformation of cord-reinforced rubber composites. New finite element models are developed for these problems. Numerical results obtained for several representative cases are presented.

By analyzing data on energetic particle fluxes measured simultaneously with detector systems on several earth satellites, signatures of different acceleration mechanisms for these particles were searched for. One of the samples is an event observed on ATS-6 and IMP-7. IMP-7 was in the dusk quarter at 38 earth radii while ATS-6 was located at local midnight at a distance of 6.6 earth radii. Although the flux variations as observed on the two spacecraft both showed 1.5 min periodicity, there was a 40-second time lag with IMP-7 behind. The results indicate that the particles are accelerated by magnetic field line annihilation, with the x-point located at about 10 earth radii.

On an out-of-the-ecliptic mission, in situ measurements of densities and temperature of interstellar neutral H and He in the heliosphere should complement observations based on backscattered Lyman-alpha intensities. A means of performing the in situ measurements is briefly described.

The field ionization neutral detector FIND operates according to the following principle: Neutral atoms are attracted toward the field ionization tips since they are polarized by the electric field of the tips. The atoms are singly ionized and repelled by the positive potential of the tips toward the detector situated behind a grid at ground potential. The ions deposit in the detector their kinetic energies, typically 26 keV, corresponding to the potential difference between the ionization region and the detector. Laboratory results show that FIND can have the resolution, sensitivities and durability required to perform in situ measurements of neutral H and He fluxes in interplanetary space, cometary halos and exospheres.

Description of laboratory model of a neutral-particle mass spectrometer that could make localized direct measurements of the neutral densities at 1 AU. Using a magnet that has a circular pole face 2.54 cm in diameter and a gap flux density of 8.0 kG, the following sensitivities were obtained for the various species: 10 to the 10th power counts per second per Torr for H2, 3 times 10 to the eighth power counts per second per Torr for He, and 2 times 10 to the 9th power counts per second per Torr for N2. Overall sensitivity may be increased through the use of electrostatic focusing elements between the ion source and the magnetic analyzer.

On April 29, 1978 a low-energetic particle event originating from a 3B flare at N22 E38 was observed by HELIOS 1 at 0.31 AU and W 38 and by HELIOS 2 at 0.29 AU and W 67 simultaneously with an interplanetary shock wave. Using plasma and magnetic field data these shocks were classified; the time variations were examined along with the angular and spectral distributions in the ecliptic plane of the intensities of 80 or more keV ions.

The time of flight technique for particle identification was extended to 2 keV/nucleon and the size of the start-time detector was reduced considerably by the use of carbon foils of few micrograms/cm square in thickness combined with microchannel plates for detecting secondary electrons. Time of flight telescopes incorporating this start-time device were used to measure the stopping power of a number of low energy heavy ions in thin carbon foils and the charge states of these ions emerging from such foils. Applications for the detection and identification of low energy interplanetary and magnetospheric particles are suggested.

A neutral particle detector using field ionization and liquid-nitrogen-cooled solid state detectors has been developed to investigate the pulse height defect of low energy He(+) ions in silicon. The field ionization tip is fabricated from a short length of tungsten wire; a typical detector is of the surface-barrier Si type with a thin Au window and low capacitance (about 15 pF); the assembly is mounted in a cylindrical vacuum chamber which is pumped to about 10 to the -6th torr. As an example, pulse height distributions for He(+) and H(+) taken at tip bias of 25.7 kV are presented. The detected energy of the H and He ions are studied as a function of the actual energies of accelerated ions, for a detector having a gold front conducting surface layer of 32 micrograms/sq cm.

Differential energy spectrum of cosmic ray deuterium flux of galactic origin from IMP 5 satellite measurements

Interplanetary space He 3 differential energy spectrum in 10-30 MeV range from IMP 4 satellite observations, noting solar contributions

Solar flare relative abundance and energy spectra of He 3 and He 4, using charged particle telescope on IMP 4

Solar modulation of low energy galactic cosmic ray protons, proposing diffusion-convection model with time dependent diffusion coefficient

Cosmic ray origin in high energy H and He isotropic abundances and energy spectra observed by IMP 4 satellite during solar quiet times

Telescopic measurements of primary energy spectra and abundances of He 3 and He 4 in galactic cosmic radiation

A self-consistent model for cosmic-ray hydrogen and helium propagation from the sources in the Galaxy to the orbit of earth is obtained, taking into account experimental information now available on the isotopes H-1, H-2, He-3, and He-4. The only adjustable parameters include the shape of the energy spectra of H-1 and He-4 at the time of source injection, the distribution of particle path lengths in interstellar space, and the solar modulation parameters. It is found that the allowed form of the source differential spectra of the H-1 and He-4 nuclei is dominated by a power law in total energy.