A high enthalpy plasma generator for entry heating simulation
High enthalpy plasma generator for atmospheric entry heating simulation
SEARCH · Engineering Papers
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.
High enthalpy plasma generator for atmospheric entry heating simulation
HPHall software simulates and tracks the time-dependent evolution of the plasma and erosion processes in the discharge chamber and near-field plume of Hall thrusters. HPHall is an axisymmetric solver that employs a hybrid fluid/particle-in-cell (Hybrid-PIC) numerical approach. HPHall, originally developed by MIT in 1998, was upgraded to HPHall-2 by the Polytechnic University of Madrid in 2006. The Jet Propulsion Laboratory has continued the development of HPHall-2 through upgrades to the physical models employed in the code, and the addition of entirely new ones. Primary among these are the inclusion of a three-region electron mobility model that more accurately depicts the cross-field electron transport, and the development of an erosion sub-model that allows for the tracking of the erosion of the discharge chamber wall. The code is being developed to provide NASA science missions with a predictive tool of Hall thruster performance and lifetime that can be used to validate Hall thrusters for missions.
Attention is given to the philosophical as well as physical principles that are essential to the establishment of MHD simulation studies for solar plasma research, assuming the capabilities of state-of-the-art computers and emphasizing the importance of 'local' MHD simulation. Solar-terrestrial plasma space is divided into several elementary regions where a macroscopic elementary energy conversion process could conceivably occur; the local MHD simulation is defined as self-contained in each of the regions. The importance of, and the difficulties associated with, the boundary condition are discussed in detail. The roles of diagnostics and of the finite difference method are noted.
We present laboratory spectra of dielectronic recombination (DR) satellite transitions attached to the He-like and H-like iron resonance lines obtained with the NASA Goddard Space Flight Center X-ray calorimeter and produced by a thermal plasma simu1ation technique on the EBIT-I electron beam ion trap at the Lawrence Livermore National Laboratory. We demonstrate that the calorimeter has sufficient spectral resolution in the 6-9 keV range to provide reliable measurements not only of standard DR satellite to resonance line intensities but also of DR satellite to DR satellite ratios that can be used to diagnose nonthermal electron distributions. Electron temperatures derived from the measured line intensities are consistent with the temperature of the simulated plasma. Temperature measurements based on DR satellite transitions have significant advantages over those based on collisional ionization equilibrium or continuum shape. Thus, successful demonstration of this method with the X-ray calorimeter is an important step fur its application in X-ray astronomy.
A three-dimensional electrostatic particle-in-cell (PIC) plasma simulation code has been developed on coarse-grain distributed-memory massively parallel computers with message passing communications. Our implementation is the generalization to three-dimensions of the general concurrent particle-in-cell (GCPIC) algorithm. In the GCPIC algorithm, the particle computation is divided among the processors using a domain decomposition of the simulation domain. In a three-dimensional simulation, the domain can be partitioned into one-, two-, or three-dimensional subdomains ("slabs," "rods," or "cubes") and we investigate the efficiency of the parallel implementation of the push for all three choices. The present implementation runs on the Intel Touchstone Delta machine at Caltech; a multiple-instruction-multiple-data (MIMD) parallel computer with 512 nodes. We find that the parallel efficiency of the push is very high, with the ratio of communication to computation time in the range 0.3%-10.0%. The highest efficiency (> 99%) occurs for a large, scaled problem with 64(sup 3) particles per processing node (approximately 134 million particles of 512 nodes) which has a push time of about 250 ns per particle per time step. We have also developed expressions for the timing of the code which are a function of both code parameters (number of grid points, particles, etc.) and machine-dependent parameters (effective FLOP rate, and the effective interprocessor bandwidths for the communication of particles and grid points). These expressions can be used to estimate the performance of scaled problems--including those with inhomogeneous plasmas--to other parallel machines once the machine-dependent parameters are known.
In view of the recent remarkable advancement of computer technology and simulation software, simulation studies are one of the most powerful academic tools for establishment of quantitative space physics and modelling of our space environment. The complex nature encountered in space plasma physics has motivated considerable development in computer simulations, which have played an essential role in the development of space plasma theory. This report describes research undertaken to understand physical processes involved in plasma waves observed in the magnetospheric plasmas, and associated nonlinear phenomena such as heating, diffusion, and acceleration of particles due to excited waves. The research explains and clarifies the observational data both qualitatively and quantitatively.
Particle free path estimates for axisymmetric high specific impulse plasma accelerators, discussing collisionless plasma simulated by computer
A closed pore ceramic foam insulation material (CPI) has been developed that offers possibilities for use as a reusable external heat shield for the NASA manned space shuttle. The outstanding characteristics of CPI are: (1) negligible water absorption due to a noninterconnecting network of cells; (2) high emittance at room and elevated temperature; (3) ability to survive at least 10 simulated reentry cycles to 1500 K using radiant heat lamps to simulate the reentry heat fluxes; (4) ability to survive, with no change in properties or appearance, at least 10 simulated plasma arc jet cycles to 1500 K (with the exception of some stress cracks induced either by the unduly severe nature of the initial arc splash heating pulse or by improper mechanical holding of the specimen in the test fixture); (5) strength (flexure); and (6) a low thermal conductivity throughout the temperature range of interest for the space shuttle.
Preliminary measurements have been made in a large laboratory discharge device configured to simulate plasma wake phenomena. A large relative electron drift impinging on a nonconducting disc gives rise to a perturbed density struture downstream, as well as a reflected beam of electrons upstream. Velocity distribution measurements were made using a novel energy analyzer with angular resolution.
Microscopic dissipation processes in quasi-perpendicular shocks are studied by two-dimensional plasma simulations in which electrons and ions are treated as particles moving in self-consistent electric and magnetic fields. Cross-field currents induce substantial turbulence at the shock front reducing the reflected ion fraction, increasing the bulk ion temperature behind the shock, doubling the average magnetic ramp thickness, and enhancing the upstream field aligned electron heat flow. The short scale length magnetic fluctuations observed in the bow shock are probably associated with this turbulence.
Plasma contactors could be used to ground satellites to space plasma to acquire a flow of electrons to propel or power the satellites. A tether would cut across geomagnetic field lines, producing a potential difference between the ends of the tether. Closing the connection between the ends would form a circuit in which an electrical load could be inserted. Design constraints of the circuit are low impedance and a fully reversible high current. The contactor would generate a neutral plasma to connect to the ionospheric plasma. The surface area of the connection would have to be kept large enough for the current density to be equal to the random electron current density in the unperturbed space plasmas. The other contactor would feed electrons and draw ions from the space plasma. Experimental results from spaceborne and ground-based space plasma simulator tests of hollow cathodes that have shown that multiampere currents can be collected are described.
Observations of a 100-mA, 1-keV electron beam were made with a low light level TV system and other instruments aboard the space shuttle STS 3 in March 1982. When fired at a moderate pitch angle with respect to the magnetic field, the beam electrons initially moved outward along a classical helical trajectory. Within a short distance, however, the beam electrons were deflected away from the helical trajectory into paths parallel to magnetic field lines threading the primary helix. This resulted in the formation of a thin-walled, cylindrical electron beam structure having a radius equal to the primary beam gyroradius. This structure, which is consistent with earlier laboratory vacuum chamber observations and recent two-dimensional computer electrostatic plasma simulations, has important implications for determining the extent to which pulsed electron beams can radiate electromagnetic energy in the ionosphere.
The electron heating required if protons scatter elastically in a parallel, collisionless shock is calculated. Near-elastic proton scattering off large amplitude background magnetic field fluctuations might be expected if the waves responsible for the shock dissipation are generated by the fire hose instability. The effects of an electrostatic potential jump in the shock layer are included by assuming that the energy lost by protons in traversing the potential jump is converted into electron thermal pressure. It is found that the electron temperature increase is a strong function of the potential jump. Comparison is made to the parallel shock plasma simulation of Quest (1987).
The present evaluation of the basic techniques and illustrative results of charged particle-modeling numerical codes suitable for particle acceleration at oblique, fast-mode collisionless shocks emphasizes the treatment of ions as test particles, calculating particle dynamics through numerical integration along exact phase-space orbits. Attention is given to the acceleration of particles at planar, infinitessimally thin shocks, as well as to plasma simulations in which low-energy ions are injected and accelerated at quasi-perpendicular shocks with internal structure.
The exotic phenomenon of energetic ion-conic and counterstreaming electron formation by lower hybrid waves along discrete auroral field lines in the earth magnetosphere is considered. Mean-particle calculations, plasma simulations, and analytical treatments of the acceleration processes are described. It is shown that, in the primary auroral electron-beam region, lower hybrid waves could be an efficient mechanism for the transverse heating of H (+) and O(+) ions of ionospheric origin, as well as for the field-aligned heating of the ambient electrons leading to coincident counterstreaming electron distributions. For O(+) ions to be energized by such a wave-particle interaction process, however, some sort of preheating mechanism is required.
Experiments have shown that plasma environments can induce discharges in solar arrays. These plasmas simulate the environments found in low earth orbits where current plans call for operation of very large power systems. The discharges could be large enough to couple into the power system and possibly disrupt operations. Here, the general concepts of the discharge mechanism and the techniques of coupling are discussed. Data from both ground and flight experiments are reviewed to obtain an expected basis for the interactions. These concepts were applied to the Space Station solar array and distribution system as an example of the large space power system. The effect of discharges was found to be a function of the discharge site. For most sites in the array discharges would not seriously impact performance. One location at the negative end of the array was identified as a position where discharges could couple to charge stored in system capacitors. This latter case could impact performance.
A parametric study of the evolution within, and signatures at, 1 AU of high-speed streams is performed with the use of a MHD two-and-a-half-dimensional time-dependent model. This study is an extension of an earlier one by Smith and Dryer (1990) who examined the ecliptic plane consequences of relatively short-duration, energetic solar disturbances. The present study examines both the erupting and corotating parts of long-duration, high-speed streams characteristic of coronal hole flows. By examining the variation of the simulated plasma velocity, density, temperature, and magnetic field at 1 AU, as well as the location of the solar coronal hole sources relative to the observer at 1 AU, it was possible to provide some insight into the identification of the solar sources of interplanetary disturbances. Two definitions for angle locating the solar source of interplanetary disturbances at 1 AU are presented and discussed. The results are applied to the suggestion by Hewish (1988) that low-latitude coronal holes are suitably positioned to be the sources of major geomagnetic storms when the holes are in the eastern half of the solar hemisphere at the time of the commencement of the storm. The results indicate that, for these cases, the streams emanating from within the hole must be very fast, greater than 1000 km/s, or very wide, greater than 60 deg, at the inner boundary of 18 solar radii.
Linear and nonlinear electron damping of the whistler precursor wave train to low Mach number quasi-perpendicular oblique shocks is studied using a one-dimensional electromagnetic plasma simulation code with particle electrons and ions. In some parameter regimes, electrons are observed to trap along the magnetic field lines in the potential of the whistler precursor wave train. This trapping can lead to significant electron heating in front of the shock for low beta(e). Use of a 64-processor hypercube concurrent computer has enabled long runs using realistic mass ratios in the full particle in-cell code and thus simulate shock parameter regimes and phenomena not previously studied numerically.