Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electron Fluid Model”

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 109 records · Page 6

Some physical implications of recent solar wind measurements

The physical implications of the existence at about 1 AU of a quiet solar wind particle flux about 90 percent larger than that suggested in the past is investigated within the framework of the two-fluid solar wind model equations. During the spherically symmetric radial expansion of the quiet solar wind, the particle flux is conserved quantity. It is found that a pure collisional two-fluid model provides good particle density and streaming velocity at 1 AU, but predicts too large an electron temperature and too small a proton temperature. When noncollisional contributions to the transport coefficients are incorporated in the model equations, a complete satisfactory agreement with the available observations is obtained. Upper limits to the effective coupling between electrons and protons, as well as to the effective proton thermal conductivity, and both upper and lower limits to the effective electron thermal conductivity in the quiet solar wind, required to provide agreement with observations, are given.

Cuperman, S.↗

Hybrid simulations of the saturated Farley-Buneman instability in the ionosphere

Numerical simulations of the Farley-Buneman instability in 2-1/2 dimensions using particle ions and fluid electrons show the growth, saturation and nonlinear behavior of two-stream waves. This hybrid technique models the saturated state of the instability for a much longer period of time than the pure particle codes that preceded it. While focusing principally on modeling the topside E region equatorial electrojet, many of these results apply to the auroral two-stream instability as well. The following features are seen in all our hybrid simulations: (1) wave growth at an angle offset from the electron drift direction where the angle depends on the strength of the driving electric field, (2) nonlinear coupling to waves traveling perpendicular to the propagation direction of the principal two-stream waves, (3) a saturated wave phase velocity at or above the sound speed but well below the velocity predicted by linear theory and (4) phase velocities which remain almost constant as a simulated radar sweeps from a horizontal direction to nearly vertical. The nonlinear electron motion dominates the behavior of these waves. Further, these simulations indicate that ion kinetic effects are not essential for the saturation of the instability and that electron temperature effects have a minor impact on the final saturated state.

Oppenheim, Meers↗

Finite Element Models for Electron Beam Freeform Fabrication Process

Electron beam freeform fabrication (EBF3) is a member of an emerging class of direct manufacturing processes known as solid freeform fabrication (SFF); another member of the class is the laser deposition process. Successful application of the EBF3 process requires precise control of a number of process parameters such as the EB power, speed, and metal feed rate in order to ensure thermal management; good fusion between the substrate and the first layer and between successive layers; minimize part distortion and residual stresses; and control the microstructure of the finished product. This is the only effort thus far that has addressed computer simulation of the EBF3 process. The models developed in this effort can assist in reducing the number of trials in the laboratory or on the shop floor while making high-quality parts. With some modifications, their use can be further extended to the simulation of laser, TIG (tungsten inert gas), and other deposition processes. A solid mechanics-based finite element code, ABAQUS, was chosen as the primary engine in developing these models whereas a computational fluid dynamics (CFD) code, Fluent, was used in a support role. Several innovative concepts were developed, some of which are highlighted below. These concepts were implemented in a number of new computer models either in the form of stand-alone programs or as user subroutines for ABAQUS and Fluent codes. A database of thermo-physical, mechanical, fluid, and metallurgical properties of stainless steel 304 was developed. Computing models for Gaussian and raster modes of the electron beam heat input were developed. Also, new schemes were devised to account for the heat sink effect during the deposition process. These innovations, and others, lead to improved models for thermal management and prediction of transient/residual stresses and distortions. Two approaches for the prediction of microstructure were pursued. The first was an empirical approach involving the computation of thermal gradient, solidification rate, and velocity (G,R,V) coupled with the use of a solidification map that should be known a priori. The second approach relies completely on computer simulation. For this purpose a criterion for the prediction of morphology was proposed, which was combined with three alternative models for the prediction of microstructure; one based on solidification kinetics, the second on phase diagram, and the third on differential scanning calorimetry data. The last was found to be the simplest and the most versatile; it can be used with multicomponent alloys and rapid solidification without any additional difficulty. For the purpose of (limited) experimental validation, finite element models developed in this effort were applied to three different shapes made of stainless steel 304 material, designed expressly for this effort with an increasing level of complexity. These finite element models require large computation time, especially when applied to deposits with multiple adjacent beads and layers. This problem can be overcome, to some extent, by the use of fast, multi-core computers. Also, due to their numerical nature coupled with the fact that solid mechanics- based models are being used to represent the material behavior in liquid and vapor phases as well, the models have some inherent approximations that become more pronounced when dealing with multi-bead and multi-layer deposits.

Chandra, Umesh↗

First-Principles Modeling of IAT-Driven Anomalous Resistivity in Hollow Cathode Discharges I: Theory

Self-consistent expressions are derived for the impact of ion acoustic turbulence (IAT) on the fluid properties of electrons and ions in the plume of a hollow cathode. This instability is believed to have a dominant role in the plasmadynamics in this region. The random phase approximation is applied to the perturbed Boltzmann equation to determine the effective drag and heating of the plasma species induced by the IAT. These driving terms are then incorporated into momentum and energy fluid equations for the ions and electrons. A fluid-like relation governing the evolution of the total IAT wave energy density is also presented. The physical significance of the IAT-based terms is discussed in the context of energy and momentum conservation, and a path forward is outlined for implementing the effects of IAT in the fluid-based model in Part II.

Jorns, Benjamin A.↗

Self-Consistent Implementation of a Zero-Equation Transport Model Into a Predictive Model for a Hall Effect Thruster

The performance of an axisymmetric multi-fluid Hall thruster code that incorporates a self-consistent, data-driven closure model for the anomalous electron transport is investigated. Five different operating conditions of the H9 magnetically shielded Hall thruster are simulated with the Jet Propulsion Laboratory’s Hall2De. In order to capture the inherent uncertainty associated with the closure model, O(100) simulations are run for each condition, each of using a coefficient set sampled randomly from a probability distribution. The results of these simulations provide probabilistic predictions of thruster performance quantities including thrust, and discharge current, as well as several component efficiencies and centerline plasma properties. The model is found to yield converged solutions at all conditions, with large 10 kHzrange oscillations and performance trends with voltage and flow rate similar to experiment. The model under-predicts the thrust by 15-25% and over-predicts the discharge current by 20% on average compared to experiments at the same discharge voltage and mass flow rate. This performance discrepancy is due to lower beam utilization, mass utilization, and divergence efficiency than experiment, resulting from high Hall parameters in the acceleration region, which lead to a protracted ion acceleration region. The physical processes underlying this result are discussed in the context of future data-driven modeling efforts.

Jorns, Benjamin A.↗

Improved Chemistry and Attenuation Models for Communication Blackout Simulation During Mars 2020 Entry

As a blunt body enters a planetary atmosphere, a plasma forms in the hypersonic shock layer and attenuates radio communication causing signal blackout for some duration of the entry sequence. In our previous work,1 computational fluid dynamics (CFD) was applied to model the entry flow around the Mars 2020 spacecraft, including ionization and electron density throughout the flow field, and predict ultra-high frequency (UHF) radio wave attenuation due to electrons. In total, 17 chemical species and their spatial profiles are modelled around the Mars 2020 spacecraft at 11 different points in time during entry. Although the simulation predicted the onset of attenuation well, the timing of the end of the predicted blackout window significantly preceded the end time observed during the 2021 landing. The present work seeks to improve the attenuation model by accounting for the fact that electrons undergo collisions with heavier species in the flow, which is an effect that was neglected in previous analyses. It is determined that including electron collisions increases the overall magnitude of attenuation predicted especially towards the end of the measured attenuation period, improving qualitative agreement between predicted and measured attenuation to both spacecraft receiving the signal from Mars 2020. To explore the remaining uncertainty in signal attenuation predictions further, a sensitivity study is performed to investigate the impact of associative ionization and electron-impact ionization rate coefficients on the electron density predicted by CFD and on the resulting attenuation predictions. These coefficients are believed to contain up to order-of-magnitude uncertainty, and therefore may significantly affect the number density of electrons throughout the flow field. Variations in associative ionization coefficients demonstrate significant impact on the magnitude of attenuation due to variation in the electron density coming from associative ionization. However, the start and end times of the predicted signal attenuation period are only slightly impacted by said variation.

Eve Papajak↗

Improved Chemistry and Attenuation Models for Communication Blackout Simulation During Mars 2020 Entry

As a blunt body enters a planetary atmosphere, a plasma forms in the hypersonic shock layer and attenuates radio communication causing signal blackout for some duration of the entry sequence. In our previous work,1 computational fluid dynamics (CFD) was applied to model the entry flow around the Mars 2020 spacecraft, including ionization and electron density throughout the flow field, and predict ultra-high frequency (UHF) radio wave attenuation due to electrons. In total, 17 chemical species and their spatial profiles are modelled around the Mars 2020 spacecraft at 11 different points in time during entry. Although the simulation predicted the onset of attenuation well, the timing of the end of the predicted blackout window significantly preceded the end time observed during the 2021 landing. The present work seeks to improve the attenuation model by accounting for the fact that electrons undergo collisions with heavier species in the flow, which is an effect that was neglected in previous analyses. It is determined that including electron collisions increases the overall magnitude of attenuation predicted especially towards the end of the measured attenuation period, improving qualitative agreement between predicted and measured attenuation to both spacecraft receiving the signal from Mars 2020. To explore the remaining uncertainty in signal attenuation predictions further, a sensitivity study is performed to investigate the impact of associative ionization and electron-impact ionization rate coefficients on the electron density predicted by CFD and on the resulting attenuation predictions. These coefficients are believed to contain up to order-of-magnitude uncertainty, and therefore may significantly affect the number density of electrons throughout the flow field. Variations in associative ionization coefficients demonstrate significant impact on the magnitude of attenuation due to variation in the electron density coming from associative ionization. However, the start and end times of the predicted signal attenuation period are only slightly impacted by said variation.

Eve Papajak↗

The spiral field inhibition of thermal conduction in two-fluid solar wind models

The paper reports on two-field models which include the inhibition of thermal conduction by the spiraling interplanetary field to determine whether any of the major conclusions obtained by Nerney and Barnes (1977) needs to be modified. Comparisons with straight field line models reveal that for most base conditions, the primary effect of the inhibition of thermal conduction is the bottling-up of heat in the electrons as well as the quite different temperature profiles at a large heliocentric radius. The spiral field solutions show that coronal hole boundary conditions do not correspond to states of high-speed streams as observed at 1 AU. The two-fluid models suggest that the spiral field inhibition of thermal conduction in the equatorial plane will generate higher gas pressures in comparison with flows along the solar rotation axis (between 1 and 10 AU). In particular, massive outflows of stellar winds, such as outflow from T Tauri stars, cannot be driven by thermal conduction. The conclusions of Nerney and Barnes remain essentially unchanged.

Nerney, S.↗

Hybrid simulation techniques applied to the earth's bow shock

The application of a hybrid simulation model, in which the ions are treated as discrete particles and the electrons as a massless charge-neutralizing fluid, to the study of the earth's bow shock is discussed. The essentials of the numerical methods are described in detail; movement of the ions, solution of the electromagnetic fields and electron fluid equations, and imposition of appropriate boundary and initial conditions. Examples of results of calculations for perpendicular shocks are presented which demonstrate the need for a kinetic treatment of the ions to reproduce the correct ion dynamics and the corresponding shock structure. Results for oblique shocks are also presented to show how the magnetic field and ion motion differ from the perpendicular case.

Winske, D.↗

Three-Fluid Magnetohydrodynamic Modeling of the Solar Wind in the Outer Heliosphere

We have developed a three-fluid, fully three-dimensional magnetohydrodynamic model of the solar wind plasma in the outer heliosphere as a co-moving system of solar wind protons, electrons, and interstellar pickup protons, with separate energy equations for each species. Our approach takes into account the effects of electron heat conduction and dissipation of Alfvenic turbulence on the spatial evolution of the solar wind plasma and interplanetary magnetic fields. The turbulence transport model is based on the Reynolds decomposition of physical variables into mean and fluctuating components and uses the turbulent phenomenologies that describe the conversion of fluctuation energy into heat due to a turbulent cascade. We solve the coupled set of the three-fluid equations for the mean-field solar wind and the turbulence equations for the turbulence energy, cross helicity, and correlation length. The equations are written in the rotating frame of reference and include heating by turbulent dissipation, energy transfer from interstellar pickup protons to solar wind protons, and solar wind deceleration due to the interaction with the interstellar hydrogen. The numerical solution is constructed by the time relaxation method in the region from 0.3 to 100 AU. Initial results from the novel model are presented.

Usmanov, Arcadi V.↗

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations↗

Interaction of the solar wind with interstellar neutral hydrogen - Three-fluid model

It is commonly assumed in models of the solar wind-interstellar neutral hydrogen interaction that the ionized interstellar particles are quickly assimilated into the solar wind proton population and 'become indistinguishable' from the original solar wind. This assumption leads to the prediction that the solar wind proton temperature should increase with radius in the outer heliosphere. This temperature increase has not been observed. It is pointed out that assimilation of the interstellar particles to the point of indistinguishability takes place on the very long Coulomb collision time scale, and is not expected to occur within the heliosphere. Results are presented of a three-fluid model of the solar wind which consists of comoving thermal populations of protons of solar origin, protons produced by ionization of interstellar hydrogen, and electrons. The steady-state results yield a solar wind with a 'core' proton distribution which cools adiabatically, and a 'halo' of interstellar pickup protons which is maintained near 10 to the 7th K by the energy input of continued ionization and pickup. Such a distribution will not be observed to manifest the temperature increase at large heliocentric distances which is predicted from a one-fluid analysis. Further time-dependent calculations show a strong correlation between the densities of the solar wind and the interstellar pickup protons. It is suggested that the interstellar pickup population may be observable by the Voyager plasma instruments in low resolution mode during periods of high solar wind density and low solar wind temperature.

Isenberg, P. A.↗

A model of Comet P/Giacobini-Zinner

A computer model of Comet P/Giacobini-Zinner is presented which contains photoprocesses, gas-phase chemical kinetics, energy balance, multifluid hydrodynamics with a transition to free molecular flow, and solar wind interaction. Recently the physics for electrons in the model has been improved by including electron impact ionization and dissociation and separately accounting for electron energetics. Electron heating and cooling mechanisms include photoprocesses, recombination processes, inelastic and elastic collisions with heavy molecules, and expansion cooling. The model incorporates an internally consistent interaction of the solar wind with the coma gas using the axisymmetric ideal fluid dynamic equations. The nuclear size and composition have been chosen to make the calculations relevant to the September 11, 1985 International Cometary Explorer (ICE) encounter with Comet P/Giacobini-Zinner. Model profiles of the temperature, velocity, and number density of the electrons are in good agreement with measurements along ICE's trajectory. The results indicate the probe passed through a region of the coma at the onset of the plasma tail.

Boice, D. C.↗

Pioneer 8 plasma-wave measurements at distant bow-shock crossings.

Evaluation of enhanced low-frequency plasma-wave levels detected near the Pioneer 8 multiple bow-shock crossings encountered beyond 120 earth radii (Bavassano et al., 1971) and of high-frequency plasma waves detected in the upstream region. It is suggested that the distant interaction of the solar wind and the magnetosheath produced nonthermal electrons of the type commonly found upstream from the subsolar shock. The bow-shock position is compared with fluid model predictions, and some distinctions between the standing bow shock and the propagating interplanetary shock are considered.

Scarf, F. L.↗

Flow properties of the solar wind obtained from white light data and a two-fluid model

The flow properties of the solar wind from 1 R(sub s) to 1 AU were obtained using a two fluid model constrained by density and scale height temperatures derived from white light observations, as well as knowledge of the electron temperature in coronal holes. The observations were obtained with the white light coronographs on SPARTAN 201-1 and at Mauna Loa (Hawaii), in a north polar coronal hole from 1.16 to 5.5 R(sub s) on 11 Apr. 1993. By specifying the density, temperature, Alfven wave velocity amplitude and heating function at the coronal base, it was found that the model parameters fit well the constraints of the empirical density profiles and temperatures. The optimal range of the input parameters was found to yield a higher proton temperature than electron temperature in the inner corona. The results indicate that no preferential heating of the protons at larger distances is needed to produce higher proton than electron temperatures at 1 AU, as observed in the high speed solar wind.

Habbal, Shadia Rifai↗

Nonthermal turbulent heating in the solar envelope.

It is shown that MHD pulses, in the form of fast magnetosonic waves or solitons, can produce a strong electron-ion coupling capable of maintaining electron-proton temperature equilibrium in the solar envelope. The mechanism producing the nonthermal heating is the fluid-like modified two-stream instability, which, since it is essentially independent of the electron-proton temperature ratio and the value of beta, becomes a prime candidate for the anomalous collisions required by the fluid models inside a distance less than 30 solar radii, in order to explain the dominant features of the solar-wind flow.

Papadopoulos, K.↗

Current status of the dissipative thermal model for solar hard X-ray bursts

The existing dissipative thermal models for hard X-ray bursts are briefly examined, and it is shown that the model with additional acceleration is the best candidate for explaining the whole hard X-ray burst. In both phases, but especially in the thermal phase, the plasma beta approaches unity, and two-dimensional modeling of the hydrodynamics is required. Following the accelerated electrons only, without taking into account the response of the bulk of the plasma, is inadequate. It is suggested that a useful approach might be a multifluid one using approximately 15 fluids with the possibility of transfer between fluids in one and two dimensions. It is concluded that, while the model is a promising one, many details remain to be worked out.

Smith, D. F.↗

Cosmic-ray acceleration during the impact of shocks on dense clouds

In order to elucidate the properties of diffusive shock acceleration in nonuniform environments, an extensive set of simulations of the dynamical interactions between plane nonradiative shocks and dense gas clouds was carried out initially in static equilibrium with their environments. These time-dependent calculations are based on the two-fluid model for diffusive cosmic ray transport, and include the dynamically active energetic proton component of the cosmic rays as well as passive electron and magnetic field components. Except when the incident shock is itself already dominated by cosmic ray pressure, it is found that the presence of the cloud adds little to the net acceleration efficiency of the original shock and can, in fact, reduce slightly the net amount of energy transferred to cosmic rays after a given time. It is found that, in 2D cloud simulations, the always-weak bow shock and the shock inside the cloud are less important to acceleration during the interaction than the tail shock.

Jones, T. W.↗