Diffusion of heavy ions in the solar corona
Diffusion of heavy ions in solar corona studied with one-fluid and two-fluid dynamic models
DYNAMIC MODEL↗
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.
Microphysics of Waves and Instabilities in the Solar Wind and Their Macro Manifestations in the Corona and Interplanetary Space
Investigations of the physical processes responsible for coronal heating and the acceleration of the solar wind were pursued with the use of our recently developed 2D MHD solar wind code and our 1D multifluid code. In particular, we explored: (1) the role of proton temperature anisotropy in the expansion of the solar (2) the role of plasma parameters at the coronal base in the formation of high (3) a three-fluid model of the slow solar wind (4) the heating of coronal loops (5) a newly developed hybrid code for the study of ion cyclotron resonance in wind, speed solar wind streams at mid-latitudes, the solar wind.
Habbal, Shadia Rifai↗
Oblique High Mach Number Heliospheric Shocks: The Role of α Particles
Spacecraft observations of heliospheric shocks often find oscillations in the magnetic field and density both, upstream and downstream. The downstream magnetic oscillations of oblique collisionless shocks were detected by Wind with 10.9 samples/s and DSCOVR spacecraft with high temporal resolution of 50 samples/s. The density oscillations associated with the shocks are also evident in proton and α particle density by Wind (with much lower temporal resolution). Recently, we have investigated low Mach number low-β oblique shock oscillations using satellite data and 2.5D hybrid modeling with electrons modeled as fluid and ions modeled as particles and found that α particles—an important component of heliospheric plasma—may affect considerably the downstream oscillations and the shock structure. The objective of the present study is to investigate the effects of α particles on high Mach number heliospheric shocks dynamics, oscillations, nonstationarity, and shock-front rippling. We extend the study to high Mach number shocks (M > 3), investigate several α particle typical densities, and compare the results for the various shock parameters. We model the effects of α particles on the shock ramp, wake, and downstream oscillation structure and the kinetic properties of proton and α particle velocity distributions at various locations downstream of the shocks. Using the 2.5D hybrid model we found that the modeled high Mach number quasi-perpendicular shock magnetic and density structures are significantly affected by α particles with typical solar wind relative abundances, suggesting that the observed high Mach number shocks are similarly affected by α particles.
Leon Ofman↗
Simulation of Gas Mixtures in RF Discharges for Nitride Deposition
Tungsten and titanium nitride films have long been grown by chemical vapor deposition (CVD) methods. However, there has been recent interest in low temperature growth using plasma enhanced CVD. For the present work, we focus on the radio frequency (BE) discharge characteristics of gas mixtures used in nitride deposition (for example, WF6 and ammonia). Because the radial variations for a standard 200 mm, parallel plate reactor are limited to a small zone near the edges of the electrodes, a 1-D (one-dimensional) analysis is considered. This model consists of a self-consistent, 3-D (three-dimensional) moment fluid simulation that solves the continuity, momentum, and energy equations for neutral and charged species. The results in terms of plasma structure, radical concentrations, and local deposition rate will be presented. We will also compare the 1-D results with those obtained from a 2-D hybrid plasma equipment model (HPEM) developed at the University of Illinois.
Hwang, Helen H.↗
Calculation of solar wind flows about terrestrial planets
A computational model was developed for the determination of the plasma and magnetic field properties of the global interaction of the solar wind with terrestrial planetary magneto/ionospheres. The theoretical method is based on an established single fluid, steady, dissipationless, magnetohydrodynamic continuum model, and is appropriate for the calculation of supersonic, super Alfvenic solar wind flow past terrestrial planets. A summary is provided of the important research results.
Stahara, S. S.↗
A mathematical model of the structure and evolution of small scale discrete auroral arcs
A three dimensional fluid model which includes the dispersive effect of electron inertia is used to study the nonlinear macroscopic plasma dynamics of small scale discrete auroral arcs within the auroral acceleration zone and ionosphere. The motion of the Alfven wave source relative to the magnetospheric and ionospheric plasma forms an oblique Alfven wave which is reflected from the topside ionosphere by the negative density gradient. The superposition of the incident and reflected wave can be described by a steady state analytical solution of the model equations with the appropriate boundary conditions. This two dimensional discrete auroral arc equilibrium provides a simple explanation of auroral acceleration associated with the parallel electric field. Three dimensional fully nonlinear numerical simulations indicate that the equilibrium arc configuration evolves three dimensionally through collisionless tearing and reconnection of the current layer. The interaction of the perturbed flow and the transverse magnetic field produces complex transverse structure that may be the origin of the folds and curls observed to be associated with small scale discrete arcs.
Seyler, C. E.↗
Titan's Plasma Environment : 3D Hybrid Kinetic Modeling of the TA Flyby and Comparison with CAPS-ELS and RPWS LP Observations
In this report we discuss the global plasma environment of the TA flyby from the perspective of 3D hybrid modeling. In our model the background, pickup, and ionospheric ions are considered as particles, whereas the electrons are described as a fluid. In homogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. Our modeling shows that mass loading of the background plasma (H(+), O(+)) by pick up ions H2(+), CH4(+) and N2(+) differs from theT9 encounter simulations when O(+) ions are not introduced into the background plasma. In our hybrid-modeling we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M=28 amu ions that were generated inside the ionosphere. Titan's interior is considered as a weakly conducting body. Special attention has been paid to comparing the simulated pickup ion density distribution with CAPS-ELS and with RPWS LP observations by the Cassini-Huygens spacecraft along theta trajectory. Our modeling shows an asymmetry of the ion density distribution and the magnetic field, including the formation of Alfven wing-like structures.
induced magnetospheres↗
Continuous blood densitometry - Fluid shifts after graded hemorrhage in animals
Rapid fluid shifts in four pigs and two dogs subjected to graded hemorrhage are investigated. Arterial blood density (BD), mean arterial pressure (MAP), central venous pressure (CVP), arterial plasma density (PD), hematocrit (Hct) and erythrocyte density were measured. The apparatus and mechancial oscillator technique for measuring density are described. Fluid shifts between red blood cells and blood plasma and alterations in the whole-body-to-large vessel Hct, F(cell) are studied using two models. The bases of the model calculations are discussed. A decrease in MAP, CVP, and BP is detected at the beginning of hemorrhaging; continued bleeding results in further BD decrease correlating with volume displacement. The data reveal that at 15 ml/kg blood loss the mean PD and BD dropped by 0.99 + or - 0.15 and 2.42 + or 0.26 g/liter, respectively, and the Hct dropped by 2.40 + or 0.47 units. The data reveal that inward-shifted fluid has a higher density than normal ultrafiltrate and/or there is a rise in the F(cell) ratio. It is noted that rapid fluid replacement ranged from 5.8 + or - 0.8 to 10.6 + or - 2.0 percent of the initial plasma volume.
Hinghofer-Szalkay, H.↗
Modeling the effects of fast shocks on solar winds ions
Observations show that, when alpha particles and other minor ions in the solar wind plasma encounter fast shocks, they are heated more than protons and their bulk motion is decelerated less than protons. These effects have been studied using a three-fluid model, and the model predictions have been compared with observations. The comparison indicates that, for supercritical fast shocks, the three-fluid model can explain cross-shock minor ion heating which is significantly greater than that of protons. When the ratio of specific heats for minor ions, gamma (alpha), equals 2, both the lesser cross-shock deceleration and the greater heating of minor ions than of protons can be predicted by the model; thus, the minor ion heating through the shock transition region is consistent with the involvement of two degrees of freedom. Because the analysis is formulated in the de Heffmann-Teller frame of reference, the method is not valid for perpendicular shocks or when the angle is large. These results agree with the few extant observations and might be confirmed by further observations at the earth's bow shock.
Zhao, Xuepu↗
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↗
Development of a computational model for predicting solar wind flows past nonmagnetic terrestrial planets
A computational model for the determination of the detailed plasma and magnetic field properties of the global interaction of the solar wind with nonmagnetic terrestrial planetary obstacles is described. The theoretical method is based on an established single fluid, steady, dissipationless, magnetohydrodynamic continuum model, and is appropriate for the calculation of supersonic, super-Alfvenic solar wind flow past terrestrial ionospheres.
Stahara, S. S.↗
Predicting the magnetospheric plasma of weather
The prediction of the plasma environment in time, the plasma weather, is discussed. It is important to be able to predict when large magnetic storms will produce auroras, which will affect the space station operating in low orbit, and what precautions to take both for personnel and sensitive control (computer) equipment onboard. It is also important to start to establish a set of plasma weather records and a record of the ability to predict this weather. A successful forecasting system requires a set of satellite weather stations to provide data from which predictions can be made and a set of plasma weather codes capable of accurately forecasting the status of the Earth's magnetosphere. A numerical magnetohydrodynamic fluid model which is used to model the flow in the magnetosphere, the currents flowing into and out of the auroral regions, the magnetopause, the bow shock location and the magnetotail of the Earth is discussed.
Dawson, John M.↗
Evidence for the presence of quasi-two-dimensional nearly incompressible fluctuations in the solar wind
Assuming that the slab and isotropic models of solar wind turbulence need modification (largely due to the observed anisotropy of the interplanetary fluctuations and the results of laboratory plasma experiments), this paper proposes a model of the solar wind. The solar wind is seen as a fluid which contains both classical transverse Alfvenic fluctuations and a population of quasi-transverse fluctuations. In quasi-two-dimensional turbulence, the pitch angle scattering by resonant wave-particle interactions is suppressed, and the direction of minimum variance of interplanetary fluctuations is parallel to the mean magnetic field. The assumed incompressibility is consistent with the fact that the density fluctuations are small and anticorrelated, and that the total pressure at small scales is nearly constant.
Matthaeus, William H.↗
Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model
The global dynamics of the ionized and neutral components in the environment of Io plays an important role in the interaction of Jupiter's corotating magnetospheric plasma with Io. The stationary simulation of this problem was done in the MHD and the electrodynamics approaches. One of the main significant results from the simplified two-fluid model simulations was a production of the structure of the double-peak in the magnetic field signature of the I0 flyby that could not be explained by standard MHD models. In this paper, we develop a method of kinetic ion simulation. This method employs the fluid description for electrons and neutrals whereas for ions multilevel, drift-kinetic and particle, approaches are used. We also take into account charge-exchange and photoionization processes. Our model provides much more accurate description for ion dynamics and allows us to take into account the realistic anisotropic ion distribution that cannot be done in fluid simulations. The first results of such simulation of the dynamics of ions in the Io's environment are discussed in this paper.
Lipatov, Alexander S.↗
Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model
The global dynamics of the ionized and neutral gases in the environment of Io plays an important role in the interaction of Jupiter s corotating magnetospheric plasma with Io. Stationary simulations of this problem have already been done using the magnetohydrodynamics (MHD) and the electrodynamics approaches. One of the major results of recent simplified two-fluid model simulations [Saur, J., Neubauer, F.M., Strobel, D.F., Summers, M.E., 2002. J. Geophys. Res. 107 (SMP5), 1-18] was the production of the structure of the double-peak in the magnetic field signature of the Io flyby. These could not be explained before by standard MHD models. In this paper, we present a hybrid simulation for Io with kinetic ions and fluid electrons. This method employs a fluid description for electrons and neutrals, whereas for ions a particle approach is used. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. Our model may provide a much more accurate description for the ion dynamics than previous approaches and allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. The first results of such a simulation of the dynamics of ions in Io s environment are discussed in this paper. Comparison with the Galileo IO flyby results shows that this approach provides an accurate physical basis for the interaction and can therefore naturally reproduce all the observed salient features.
Lipatov, Alexander S.↗
Application of a global solar wind/planetary obstacle interaction computational model: Earth, Venus, Mars, Jupiter and Saturn studies
This work relates to a series of collaborative investigations involving the application of a computational model for the determination of the detailed plasma and magnetic field properties associated with the global interaction of the solar wind with various planetary obstacles throughout the solar system. The theoretical method is based on an established single fluid, steady, dissipationless, magnetohydrodynamic continuum model, and is appropriate for the calculation of supersonic, super-Alfvenic solar wind flow past planetary obstacles. The investigations undertaken relate to studies of various solar wind interaction phenomena with Venus, Earth, Jupiter, and Saturn.
Stahara, S. S.↗
Decay of electrostatic hydrogen cyclotron waves into ion acoustic modes in auroral field lines
The coherent three-wave decay of a linearly unstable electrostatic hydrogen cyclotron (EHC) wave into stable EHC and ion acoustic modes is considered. The general problem of the three weakly interacting electrostatic normal modes in a Maxwellian plasma is discussed. EHC is examined in a fluid description, and the results are used to guide a similar study in a Vlasov plasma system intended to model the aurora acceleration region parameters. The time dependence of the decay in a simple three-wave interaction is presented in order to show how wave saturation can arise.
Bergmann, R.↗