Modeling snowflake divertors in MAST-U tokamak using UEDGE code
Not Available
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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.
Final report of activities for the LDRD-express project #223796 titled: “Fluid models of charged species transport: numerical methods with mathematically guaranteed properties”, PI: Ignacio Tomas, Co-PI: John Shadid
This report summarizes the findings and outcomes of the LDRD-express project with title “Fluid models of charged species transport: numerical methods with mathematically guaranteed properties”. The primary motivation of this project was the computational/mathematical exploration of the ideas advanced aiming to improve the state-of-the-art on numerical methods for the one-fluid Euler-Poisson models and gain some understanding on the Euler-Maxwell model. Euler-Poisson and Euler-Maxwell, by themselves are not the most technically relevant PDE plasma-models. However, both of them are elementary building blocks of PDE-models used in actual technical applications and include most (if not all) of their mathematical difficulties. Outside the classical ideal MHD models, rigorous mathematical and numerical understanding of one-fluid models is still a quite undeveloped research area, and the treatment/understanding of boundary conditions is minimal (borderline non-existent) at this point in time. This report focuses primarily on bulk-behaviour of Euler-Poisson’s model, touching boundary conditions only tangentially.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
Fluid and gyrokinetic turbulence in open field-line, helical plasmas
Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently, it has become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas [T. N. Bernard et al., Phys. Plasmas 26, 042301 (2019)]. Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs that can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, and improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
3D Modelling of the Lithium Injection Experiment in H-Mode Plasma of EAST
Investigations of lithium transport in H-mode plasma of EAST tokamak have been performed by the three-dimensional fluid code EMC3-EIRENE. The simulation results show that the toroidally-localized lithium injection at the low field side (LFS) leads to a toroidally non-axisymmetric distribution of the Li 1+ and Li 2+ ions. The field line tracing approach is utilized to carry out a detailed study on the distributions of the Li 1+ and Li 2+ ions. It is found that Li 1+ ions mainly reside at the LFS while Li 2+ ions can transport from LFS to the high field side. The relationship between mean free path and connection length for Li 1+ and Li 2+ ions has been discussed. Here, the 3D line-integration images of deuterium and lithium emission spectra modelled by EMC3-EIRENE are in qualitative consistency with the experiment results.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Validation study of RWM stability in DIII-D high- β N plasmas
The n = 1 (n is the toroidal mode number) resistive wall mode (RWM) stability is numerically investigated for two DIII-D high-β N discharges 176440 and 172461, utilizing the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) codes. Systematic validation efforts are attempted, for the first time, for discharges with very slow or vanishing toroidal flow for a large fraction of the plasma volume. While gaining physics insights in accessing stable operation regime at β N exceeding the Troyon no-wall limit in these slow-rotation experiments, the predictive capability of fluid and non-perturbative magnetohydrodynamic-kinetic hybrid models for the RWM is further confirmed. The MARS-F fluid model, with a strong but numerically tunable viscosity mimicking ion Landau damping of parallel sound waves, finds complete stabilization of the n = 1 RWM in the considered DIII-D plasmas under the experimental flow conditions. Similarly, either full stabilization (for discharge 176440) or marginal stability (for discharge 172461) of the mode is computed by the MARS-K hybrid model, which is first-principle based without free model parameters. In particular, all drift kinetic resonances, including those of thermal and energetic particles, are found to synergistically act to marginally stabilize the RWM in discharge 172461. These MARS-F/K modeling results explain the experimentally observed stable operational regime in DIII-D, as far as the RWM stability is concerned. Extensive numerical sensitivity studies, with respect to the plasma toroidal flow speed as well as the radial location of the resistive wall, are also carried out to further support the validation study.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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.↗
Implications of a two-fluid heat diffusion model for RF condensation effects in magnetic islands
Theoretical calculations have predicted an “RF condensation” effect through which the temperature perturbation and current drive may be nonlinearly enhanced in electron cyclotron (EC) heated magnetic islands. Previously, two-fluid effects had been neglected in the OCCAMI code used to model RF condensation. For realistic experimental conditions, we show that two-fluid effects in an EC heated plasma may have a major impact on RF condensation. Less collisional plasmas, where electron and ion temperatures are decoupled, exhibit an enhanced linear temperature perturbation in the island, and an even stronger nonlinear enhancement of the perturbed temperature. This can occur even when nonlinear effects are small in the single-fluid limit. We also calculate corresponding hysteresis effects. Finally, we demonstrate that other physical effects, including the rotation of the magnetic island and the Gaussian spread of the beam from a gyrotron, can significantly decrease the predicted observable RF condensation effect.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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.↗
Parameter space for magnetization effects in high-energy-density plasmas
Magnetic fields are well known to affect the evolution of fluids via the J × B force, where J is the current density and B is the magnetic field. This force leads to the influence of magnetic fields on hydrodynamics (magnetohydrodynamics). Magnetic fields are often neglected in modeling of high-energy-density plasmas, since J × B is very small compared with the plasma pressure gradients. However, many experiments lie in a separate part of parameter space where the plasma is indirectly affected via magnetization of the heat flux and charged particle transport. This is true even for initially unmagnetized plasmas, since misaligned density and temperature gradients can self-generate magnetic fields. By comparing terms in the induction equation, we go on to estimate the regions of parameter space where these self-generated fields are strong enough to affect the hydrodynamics.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas
JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Scale Separation Effects on Simulations of Plasma Turbulence
Abstract Understanding plasma turbulence requires a synthesis of experiments, observations, theory, and simulations. In the case of kinetic plasmas such as the solar wind, the lack of collisions renders the fluid closures such as viscosity meaningless and one needs to resort to higher-order fluid models or kinetic models. Typically, the computational expense in such models is managed by simulating artificial values of certain parameters such as the ratio of the Alfvén speed to the speed of light ( v A / c ) or the relative mass ratio of ions and electrons ( m i / m e ). Although, typically care is taken to use values as close as possible to realistic values within the computational constraints, these artificial values could potentially introduce unphysical effects. These unphysical effects could be significant at sub-ion scales, where kinetic effects are the most important. In this paper, we use the 10-moment fluid model in the Gkeyll framework to perform controlled numerical experiments, systematically varying the ion–electron mass ratio from a small value down to the realistic proton–electron mass ratio. We show that the unphysical mass ratio has a significant effect on the kinetic range dynamics as well as the heating of both plasma species. The dissipative process for both ions and electrons becomes more compressive in nature, although the ions remain nearly incompressible in all cases. The electrons move from being dominated by incompressive viscous-like heating/dissipation to very compressive heating/dissipation dominated by compressions/rarefactions. While the heating change is significant for the electrons, a mass ratio of m i / m e ∼ 250 captures the asymptotic behavior of electron heating.
79 ASTRONOMY AND ASTROPHYSICS↗
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.↗