Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma fluid modeling”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Structure of perpendicular shocks in collisionless plasma

The macrostructure of perpendicular shocks in the supercritical regime is investigated theoretically. The coupling between the shock precursor (foot region) associated with reflected ions and the magnetic ramp associated with electron ohmic heating is analyzed using a model in which the ions are treated as a multifluid and the electrons as a massless fluid. The results are found to be in good agreement with a number of features exhibited by recent numerical simulations of perpendicular shocks, including the morphology of the reflected-gyrating ion stream, the presence of potential and magnetic field overshoots, and the Mach number dependence of the shock structure.

Leroy, M. M.↗
NASA NTRS1992

Numerical simulation of geometric scale effects in cylindrical self-field MPD thrusters

A 2-D, two-temperature, single fluid magnetohydrodynamic code which incorporates classical plasma transport coefficients and Hall effects was developed to predict steady-state, self-field magnetoplasmadynamic (MPD) thruster performance. The governing equations and numerical methods of solution are outlined and discussed. Experimental comparisons are used to validate model predictions. The model accurately predicts thrust and reproduces trends in the discharge voltage for discharge currents below experimentally measured onset values. However, because the model does not include electrode effects, the calculated voltage drops are significantly lower than experimentally measured values. Predictions of thrust and flow efficiency are made for a matrix of fifteen cylindrical thruster geometries, assuming a fully ionized argon propellant. A maximum predicted specific impulse of 1680 s is obtained for a thruster with an anode radius of 2.5 cm, a cathode radius of 0.5 cm, and equal electrode lengths of 2.5 cm. A scaling relation is developed to predict, within limits, the onset of cylindrical, self-field thruster instability as a function of geometry and operating condition.

Lapointe, M.↗
NASA NTRS1992

Numerical simulation of geometric scale effects in cylindrical self-field MPD thrusters

A 2D, two-temperature, single fluid MHD code which incorporates classical plasma transport coefficients and Hall effects has been developed to predict steady-state, self-field MPD thruster performance. The governing equations and numerical methods of solution are outlined and discussed. Experimental comparisons are used to validate model predictions. The model accurately predicts thrust and reproduces trends in the discharge voltage for discharge currents below experimentally measured onset values. However, because the model does not include electrode effects the calculated voltage drops are significantly lower than experimentally measured values. Predictions of thrust and flow efficiency are made for a matrix of fifteen cylindrical thruster geometries assuming a fully ionized argon propellant.

Lapointe, M.↗
NASA NTRS1980

Strongly turbulent stabilization of electron beam-plasma interactions

The stabilization of electron beam interactions due to strongly turbulent nonlinearities is studied analytically and numerically for a wide range of plasma parameters. A fluid mode coupling code is described in which the effects of electron and ion Landau damping and linear growth due to the energetic electron beam are included in a phenomenological manner. Stabilization of the instability is found to occur when the amplitudes of the unstable modes exceed the threshold of the oscillating two-stream instability. The coordinate space structure of the turbulent spectrum which results clearly shows that soliton-like structures are formed by this process. Phenomenological models of both the initial stabilization and the asymptotic states are developed. Scaling laws between the beam-plasma growth rate and the fluctuations in the fields and plasma density are found in both cases, and shown to be in good agreement with the results of the simulation.

Freund, H. P.↗
NASA NTRS2014

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗
NASA NTRS1976

Electrostatic shocks as nonlinear ion acoustic waves

A cold fluid approach is presented which yields exact double-layer solutions that can be used to model electrostatic S-type shocks. Solutions derived from the two-temperature electron model are presented in order to show the range of amplitudes and scale lengths possible for this plasma model and to examine how shock properties depend on orientation of the shock. The dependence of various double-layer quantities on plasma composition is then considered for the cold electron beam model. Double-layer solutions pertinent to describing electrostatic shocks are pointed out.

Witt, E.↗
NASA NTRS1993

A comparison of arcjet plume properties to model predictions

This paper describes an experimental study of the plasma plume properties of a 1 kW class hydrogen arcjet thruster and the comparison of measured temperature and velocity field to model predictions. The experiments are based on laser-induced fluorescence excitation of the Balmer-alpha transition. The model is based on a single-fluid magnetohydrodynamic description of the flow originally developed to predict arcjet thruster performance. Excellent agreement between model predictions and experimental velocity is found, despite the complex nature of the flow. Measured and predicted exit plane temperatures are in disagreement by as much as 2000K over a range of operating conditions. The possible sources for this discrepancy are discussed.

Cappelli, M. A.↗
NASA NTRS1993

Nutrition and human physiological adaptations to space flight

Space flight provides a model for the study of healthy individuals undergoing unique stresses. This review focuses on how physiological adaptations to weightlessness may affect nutrient and food requirements in space. These adaptations include reductions in body water and plasma volume, which affect the renal and cardiovascular systems and thereby fluid and electrolyte requirements. Changes in muscle mass and function may affect requirements for energy, protein and amino acids. Changes in bone mass lead to increased urinary calcium concentrations, which may increase the risk of forming renal stones. Space motion sickness may influence putative changes in gastro-intestinal-hepatic function; neurosensory alterations may affect smell and taste. Some or all of these effects may be ameliorated through the use of specially designed dietary countermeasures.

Review↗
NASA NTRS1994

The structure of mass-loading shocks. 1: Comets

A new multifluid model to describe the solar wind-cometary ion plasma in the outer cometary coma is derived. This model is distinguished from previous multifluid models in that wave-particle effects are included explicitly. By considering hydrodynamic timescales, self-consistent dissipative terms (which correspond to the spatial diffusion of cometary ions) are derived, thereby allowing the structure of a cometary bow shock to be resolved. In the hypersonic limit, our diffusive multifluid model reveals an attractive connection to the one-fluid model formulated by Biermann et al. (1967), a connection which serves as an important guide when investigating our more complicated nonhypersonic diffusive model. The detailed model consists of cometary ions diffusing in a warm solar wind which is heated by the dissipation of Alfvenic turbulence (in situ and generated by the initial ring beam distribution). The quasi-parallel cometary shock is found to be at least an order of magnitude thicker than a quasi-perpendicular cometary shock. It is furthermore found that for undisturbed solar wind (sonic) Mach numbers in excess of approximately 6, the cometary bow shock (for both perpendicular and parallel cases) is smoothed completely by the cometary ions and no 'proton' or thermal subshock is necessary. Finally, a detailed comparison of the model with the plasma observations made at Halley by Giotto is undertaken. It is found that the observed shock structure, thickness, location, and plasma parameters compare well with the predictions of the theoretical model, particularly in the case of the quasi-parallel shock.

Zank, G. P.↗
NASA NTRS1977

Thermodynamics of Thomas-Fermi screened Coulomb systems

We obtain in closed analytic form, estimates for the thermodynamic properties of classical fluids with pair potentials of Yukawa type, with special reference to dense fully ionized plasmas with Thomas-Fermi or Debye-Hueckel screening. We further generalize the hard-sphere perturbative approach used for similarly screened two-component mixtures, and demonstrate phase separation in this simple model of a liquid mixture of metallic helium and hydrogen.

Firey, B.↗
NASA NTRS1984

Statistical mechanics of light elements at high pressure. VI - Liquid-state calculations with Thomas-Fermi-Dirac theory

A model free energy is developed for hydrogen-helium mixtures based on solid-state Thomas-Fermi-Dirac calculations at pressures relevant to the interiors of giant planets. Using a model potential similar to that for a two-component plasma, effective charges for the nuclei (which are in general smaller than the actual charges because of screening effects) are parameterized, being constrained by calculations at a number of densities, compositions, and lattice structures. These model potentials are then used to compute the equilibrium properties of H-He fluids using a charged hard-sphere model. The results find critical temperatures of about 0 K, 500 K, and 1500 K, for pressures of 10, 100, and 1000 Mbar, respectively. These phase separation temperatures are considerably lower (approximately 6,000-10,000 K) than those found from calculations using free electron perturbation theory, and suggest that H-He solutions should be stable against phase separation in the metallic zones of Jupiter and Saturn.

Macfarlane, J. J.↗
NASA NTRS1986

Simulations of high Mach number perpendicular shocks with resistive electrons

A simulation code which models the ions as microparticles and the electrons as a resistive massless fluid is employed to study the structure of high Mach number perpendicular shocks. It is found that stable stationary shock solutions can be obtained for Alfven Mach numbers (M sub A) between 5 and 60 for upstream plasmas where the ratio of the plasma pressure to the magnetic pressure is 1, providing that the upstream resistive diffusion length is much smaller than the ion inertial length. For much larger resistive diffusion lengths, the magnetic field overshoot is damped, and the imbalance in the electron momentum equation results in a periodic fluctuation of the fraction of reflected ions. In the limit of M sub A of less than 10, the magnetic overshoot and the fraction of reflected ions increase with increasing M sub A, while at higher Mach numbers the fraction of reflected ions peaks at about 40 percent and the magnetic field overshoot increases at a much slower rate. Electron inertial effects are also considered.

Quest, K. B.↗
NASA NTRS2001

Laser Powered Launch Vehicle Performance Analyses

The purpose of this study is to establish the technical ground for modeling the physics of laser powered pulse detonation phenomenon. Laser powered propulsion systems involve complex fluid dynamics, thermodynamics and radiative transfer processes. Successful predictions of the performance of laser powered launch vehicle concepts depend on the sophisticate models that reflects the underlying flow physics including the laser ray tracing the focusing, inverse Bremsstrahlung (IB) effects, finite-rate air chemistry, thermal non-equilibrium, plasma radiation and detonation wave propagation, etc. The proposed work will extend the base-line numerical model to an efficient design analysis tool. The proposed model is suitable for 3-D analysis using parallel computing methods.

Chen, Yen-Sen↗
NASA NTRS1999

On Europa's Magnetospheric Interaction: A MHD Simulation of the E4 Flyby

The global three-dimensional interaction of Europa with the Jovian magnetosphere is modeled by using a complete set of ideal magnetohydrodynamic (MHD) equations. The model accounts for exospheric mass loading, ion-neutral charge exchange, recombination, and a possible intrinsic dipole magnetic field of Europa. The single-fluid MHD equations are solved by using a modem, finite volume, higher-order, Godunov-type method on an adoptively refined unstructured grid, which allows detailed modeling of the region near Europa while still resolving both the upstream region and the satellite's wake. The magnetic field and plasma density measured during Galileo's EGA flyby of December 19, 1996, are reproduced reasonably well in the simulation. We find the agreement between the data and our model particularly convincing if we assume that the plasma velocity during the EGA flyby deviated from the nominal corotation direction by approximately 20 deg. Evidence from the Galileo energetic particle detector also supports this assumption. In this case, we can fit the data using a dipole with orientation close to that of an induced dipole arising from the interaction of a hypothetical conducting subsurface layer on Europa with the periodically changing magnetic field of Jupiter. However, the magnitude of the dipole in our model is somewhat smaller (70%) than that suggested by Khurana et al. The total mass loading and ion-neutral charge exchange rates are consistent with the estimates of Europa's atmosphere and ionosphere.

Kabin, K.↗
NASA NTRS2016

Preferential Heating of Oxygen 5+ Ions by Finite-Amplitude Oblique Alfven Waves

Minor ions in the fast solar wind are known to have higher temperatures and to flow faster than protons in the interplanetary space. In this study we combine previous research on parametric instability theory and 2.5D hybrid simulations to study the onset of preferential heating of Oxygen 5+ ions by large-scale finite-amplitude Alfven waves in the collisionless fast solar wind. We consider initially non-drifting isotropic multi-species plasma, consisting of isothermal massless fluid electrons, kinetic protons and kinetic Oxygen 5+ ions. The external energy source for the plasma heating and energization are oblique monochromatic Alfven-cyclotron waves. The waves have been created by rotating the direction of initial parallel pump, which is a solution of the multi-fluid plasma dispersion relation. We consider propagation angles theta less than or equal to 30 deg. The obliquely propagating Alfven pump waves lead to strong diffusion in the ion phase space, resulting in highly anisotropic heavy ion velocity distribution functions and proton beams. We discuss the application of the model to the problems of preferential heating of minor ions in the solar corona and the fast solar wind.

Maneva, Yana G.↗
NASA NTRS2022

Design of Hypervelocity Flow Generator (HFG) and Its Diagnostics

Ground facilities for hypersonic research are a key element for successful development of aerodynamically proven hypervelocity vehicles. Design concepts and diagnostics of a hypervelocity flow generator (HFG) were made as a test platform for hypersonic/hypervelocity spacecraft models at the NASA Langley Research Center. The HFG is a hypersonic flow field generator using optically heated gas which is blown into an 80 cu m vacuum chamber. The vacuum chamber is kept at a stable vacuum pressure with a combination of three large vacuum pumps, while the HFG is in the test mode. The HFG provides a relatively small test section with approximately a 20 cm window. This facility was designed to generate 2.45 km/sec of flow speed, and potentially generate a continuous flow with the nozzle and vacuum system. The window in test section provides a direct view of the shock wave around a model in order to measure temperature, pressure, and density profile within the shock layer. One of the key test goals under this project is to understand why the emission spectra from the standing shockwave plasma predicted by Lora-Loran codes are significantly different from the measured emission spectra from the Flight Investigation Reentry Environment (FIRE II) Flight. The correct estimate of the thermal loading on the leading edge of hypersonic vehicles greatly affects the aerodynamic design, the material selection for the vehicle, and the cooling requirement and can be obtained by the precise modeling of emission spectra from shockwave plasma. However, the estimation of thermal loading is not an easy task due to complex non-equilibrium radiative process within high temperature shock layers that still falls into a category of cold plasma. Direct flight experiments are the most desirable, but not a cost-effective approach. Analysis by computational fluid dynamics (CFD) offers many test flexibilities. However, the CFD codes must be fully tested and validated with experimental data before the codes are effectively used for practical design. Large discrepancies between experiments and codes appear in hypersonic/hypervelocity flow regimes at high altitudes of 60 km ~ 90 km. This HFG facility offers some important parameters for CFD code validation, such as collision cross-sections, relaxation times, reaction rate coefficients and transportation coefficients. The HFG test facility is based on the ejection flow of high temperature gas heated over to 3500 K through a nozzle. The tungsten gas chamber of the HFG is heated up to a desired temperature by a 60 kW optical power beam source. This system consists of an optical power source, a thermal chamber, an expansion nozzle, a test section, and an 80 cu m vacuum tank. 60 kW optical input power is obtained from the 150 kW Vortek arc lamp system (by Vortek Industries, Vancouver, Canada). This optical beam is focused to heat the gas chamber within which a flow media is heated. The maximum achievable temperature of the flow medium reached approximately 3500 K or even higher but is limited by the melting point of the chamber material used. The exhaust velocity through the nozzle was determined by the stagnation temperature and the molecular weight of the working medium at the test section. To provide design parameters, a NASA Chemical Equilibrium with Application (CEA) computer program is used for the simulation of aerothermal data. This CEA program can calculate chemical equilibrium and properties of complex mixtures using shock tube parameters. For nitrogen gas at 2666 K stagnation temperature, the maximum achievable velocity at the test section is approximately 2.45 km/sec which is within the range of the thermal velocity of 8000 K shock layer. Based on the calculation through the CEA program, the design parameters of a HFG were determined and implemented for the test section that includes an expansion nozzle. The installed test section of the experimental facility can sustain a condition of a re-entry vehicle from the Space at an altitude of 60-90 km. After installation of the HFG, the system was fully tested and its operational parameters were measured. An 80 cu m vacuum chamber of HFG was set at 1 torr level to keep a stable downstream condition. The pumping time to reach the minimum vacuum pressure (~ 1 torr) at the test chamber from atmospheric pressure was approximately an hour using Kinney (MBV-14000/MB -1600/KT-300) and Stokes (Model 1772 and 412) pumps. Such a setup condition allowed a continued stable operation of the HFG experiment with flowrates through 1-, 3-, or 5-mm diameters of nozzle throat. The flow characteristics of the HFG for various operating conditions were performed using a focused Schlieren method. At 0.8 torr chamber pressure, a barrel shock was observed at the test section. With a sphere obstacle of 12 mm diameter in the flow, a bow shock (~ 2-mm thickness) was observed by a focused Schlieren visualization method.

Hypervelocity flow↗
NASA NTRS1975

A wave model of near wakes

A new model for the near-wake studies is presented. It pertains to a tenuous plasma medium with negligible magnetic field effects. The wake-generating body moves at a mesothermal speed, i.e., intermediate between the thermal speeds of ions and electrons of the ambient plasma which is bithermal (Ti less than Te) and has a Debye length much smaller than the body size. The model is constructed by using the concept that the cavity left by a rapidly moving body acts as a generator of finite-amplitude ion acoustic waves which are excited by an initial discontinuity of plasma density at the leading edge of the cavity. This ionic wave motion, represented by fluid equations, is represented approximately by the time dependent one-dimensional Korteweg-de Vries equation. Numerical and analytic solutions are both given.

Liu, V. C.↗
NASA NTRS1988

Electron heating in superhigh Mach number shocks

Fluid and MHD models, as well as direct extrapolation of the earth's bow shock measurements in the high Mach number (HMN) range to the superhigh Mach number (SHMN) range predict that the downstream electron pressure is only a negligible fraction of the Rankine-Hugoniot downstream pressure. Following Alfven, plasma physics experimental-theoretical data combined with magnetospheric observations were used to probe the physics of the SHMN shocks. It is shown that inclusion of proper plasma physics considerations in the interaction of the reflected and transmitted ions and the electrons at the 'foot' of the shock leads to the surprising result that electron heating can dominate in the SHNM range. A stationary model of the shock structure is derived and shown to be the result of extrapolation of the high Mach number shock physics wiht incorporation of collective interactions at the foot.

Papadopoulos, K.↗