Electron Physics in 3‐D Two‐Fluid 10‐Moment Modeling of Ganymede's Magnetosphere
Explore the source record for details and available documents.
Liang Wang↗
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.
Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)
Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A cancellation problem in hybrid particle-in-cell schemes due to finite particle size
The quasi-neutral hybrid particle-in-cell algorithm with kinetic ions and fluid electrons is a popular model to study multi-scale problems in laboratory, space, and astrophysical plasmas. Here, it is shown that the different spatial discretizations of ions as finite-spatial-size particles and electrons as a grid-based fluid can lead to significant numerical wave dispersion errors in the long wavelength limit (kd i «1, where k is the wavenumber and di is the ion skin-depth). The problem occurs when high-order particle-grid interpolations, or grid-based smoothing, spreads the electric field experienced by the ions across multiple spatial cells and leads to inexact cancellation of electric field terms in the total (ion + electron) momentum equation. Practical requirements on the mesh spacing Δx/d i are suggested to bound these errors from above. The accuracy impact of not respecting these resolution constraints is shown for a non-linear shock problem.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Computational fluid dynamic modeling to determine the indoor environment of an electron-ion collider service building
The design for the Electron-Ion Collider (EIC) calls for several service buildings that house various power supplies and control electronics for the collider ring itself. In order to operate within specified conditions, the ambient air entering the power supplies needs to be within a certain temperature range while dissipating the heat from losses. Proper cooling is therefore a necessity in the service buildings to ensure that every aspect of the EIC works as intended. Since the EIC is in the design stage, we are evaluating the indoor environment of the service building using the current design specifications. We have researched multiple forms of literature and performed the necessary calculations to compile a list of boundary conditions that accurately represent the situation at hand. We are using computational fluid dynamics modeling to solve the conservation equations for mass, momentum, and energy (Navier-Stokes). This allows us to perform a finite element analysis which will give us the flow distribution in the room as well as temperature profiles throughout the building. We have obtained a simulation result giving us the temperature profiles for the building and it shows that the placement of the racks and supply vents are essential to obtaining an even temperature distribution. This model will provide a basis for design decisions which will affect the overall cooling of the service building without extending the schedule and avoiding a costly reworking of the cooling system.
43 PARTICLE ACCELERATORS↗
Multifluid Simulations of the Global Solar Wind Including Pickup Ions and Turbulence Modeling
I will describe a three-dimensional magnetohydrodynamic model of the solar wind that takes into account turbulent heating of the wind by velocity and magnetic fluctuations as well as a variety of effects produced by interstellar pickup protons. The interstellar pickup protons are treated in the model as one fluid and the protons and electrons are treated together as a second fluid. The model equations include a Reynolds decomposition of the plasma velocity and magnetic field into mean and fluctuating quantities, as well as energy transfer from interstellar pickup protons to solar wind protons that results in the deceleration of the solar wind. The model is used to simulate the global steady-state structure of the solar wind in the region from 0.3 to 100 AU. The simulation assumes that the background magnetic field on the Sun is either a dipole (aligned or tilted with respect to the solar rotation axis) or one that is deduced from solar magnetograms.
Goldstein, Melvyn L.↗
Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D
The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Inertial and anisotropic pressure effects on cross-field electron transport in low-temperature magnetized plasmas
In this paper, a one-dimensional (1D) particle-in-cell Monte Carlo collision (PIC-MCC) model is developed to investigate the effects of anisotropic pressure and inertial terms due to non-Maxwellian velocity distribution functions on cross-field electron transport. The conservation of momentum is evaluated by taking the moments of the first-principles gas-kinetic equation. A steady-state discharge is obtained without any low-frequency ionization oscillations by considering an anomalous electron scattering profile. Further, the results obtained from the 1D PIC-MCC model are compared with fluid models, including the quasi-neutral drift-diffusion (DD), non-neutral DD, and full fluid moment models. The discharge current obtained from the PIC-MCC model is in good agreement with the fluid models. The cross-field electron transport due to the inertial terms, i.e. the gradient of axial and azimuthal drift, is evaluated. Moreover, PIC-MCC simulation results show non-zero, anisotropic, off-diagonal pressure tensor terms due to asymmetric non-Maxwellian electron velocity distribution function, potentially contributing to cross-field electron transport.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Transient behavior of drift and ionization in atmospheric pressure nitrogen discharge
The fluid models are frequently used to describe a non-thermal plasma such as a streamer discharge. The required electron transport data and rate coefficients for the fluid model are parametrized using the local field approximation (LFA) in first order models and the local-mean-energy approximation (LMEA) in second order models. We performed Monte Carlo simulations in Nitrogen gas with step changes in the E/N (reduced electric field) to study the behavior of the transport properties in the transient phase. During the transient phase of the simulation, we extract the instantaneous electron mean energy, which is different from the steady state mean electron energy, and the corresponding transport parameters and rate coefficients. Our results indicate that the mean electron energy is not a suitable parameter for mobility/drift of electrons due to big difference in momentum relaxation and energy relaxation. Furthermore, the high energy threshold rates such as ionization show a strong correlation to mean electron energy. In second order models where the energy-balance equation is solved, we suggest that it would rather be appropriate to use the local electric field to find electron drift velocity in gases such as Nitrogen and the local mean electron energy to determine the ionization and excitation rates.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A Two-Fluid, MHD Coronal Model
We describe first results from a numerical two-fluid MHD model of the global structure of the solar corona. The model is two-fluid in the sense that it accounts for the collisional energy exchange between protons and electrons. As in our single-fluid model, volumetric heat and momentum sources are required to produce high speed wind from coronal holes, low speed wind above streamers, and mass fluxes similar to the empirical solar wind. By specifying different proton and electron heating functions we obtain a high proton temperature in the coronal hole and a relatively low proton temperature in the streamer (in comparison with the electron temperature). This is consistent with inferences from SOHO/UVCS, and with the Ulysses/SWOOPS proton and electron temperature measurements which we show from the fast latitude scan. The density in the coronal hole between 2 solar radii and 5 solar radii (2RS and 5RS) is similar to the density reported from SPARTAN 201-01 measurements by Fisher and Guhathakurta. The proton mass flux scaled to 1 AU is 2.4 x 10(exp 8)/sq cm s, which is consistent with Ulysses observations. Inside the closed field region, the density is sufficiently high so that the simulation gives equal proton and electron temperatures due to the high collision rate. In open field regions (in the coronal hole and above the streamer) the proton and electron temperatures differ by varying amounts. In the streamer, the temperature and density are similar to those reported empirically by Li et al and the plasma beta is larger than unity everywhere above approx. 1.5 R(sub s), as it is in all other MHD coronal streamer models.
Suess, Steven T.↗
A Two-Fluid, MHD Coronal Model
We describe first results from a numerical two-fluid MHD model of the global structure of the solar Corona. The model is two-fluid in the sense that it accounts for the collisional energy exchange between protons and electrons. As in our single-fluid model, volumetric heat and Momentum sources are required to produce high speed wind from Corona] holes, low speed wind above streamers, and mass fluxes similar to the empirical solar wind. By specifying different proton and electron heating functions we obtain a high proton temperature in the coronal hole and a relatively low proton temperature above the streamer (in comparison with the electron temperature). This is consistent with inferences from SOHO/UltraViolet Coronagraph Spectrometer instrument (UVCS), and with the Ulysses/Solar Wind Observations Over the Poles of the Sun instrument (SWOOPS) proton and electron temperature measurements which we show from the fast latitude scan. The density in the coronal hole between 2 and 5 solar radii (2 and 5 R(sub S)) is similar to the density reported from SPARTAN 201.-01 measurements by Fisher and Guhathakurta [19941. The proton mass flux scaled to 1 AU is 2.4 x 10(exp 8)/sq cm s, which is consistent with Ulysses observations. Inside the closed field region, the density is sufficiently high so that the simulation gives equal proton and electron temperatures due to the high collision rate. In open field regions (in the coronal hole and above the streamer) the proton and electron temperatures differ by varying amounts. In the streamer the temperature and density are similar to those reported empirically by Li et al. [1998], and the plasma beta is larger than unity everywhere above approx. 1.5 R(sub S), as it is in all other MHD coronal streamer models [e.g., Steinolfson et al., 1982; also G. A. Gary and D. Alexander, Constructing the coronal magnetic field, submitted to Solar Physics, 1998].
Suess, S. T.↗
Thermal Management System of an Outer-Rotor-Motor-Based Traction Drive With Integrated Power Electronics in its Central Cavity
Increasing power densities of electric vehicle traction drive systems necessitates combining the electric motor and the power electronics into one unit. A compact, integrated traction drive unit with fewer components also drives production costs down, enabling wider adoption of electric vehicles. However, the integration of power electronics in the electric machine is associated with challenges of designing an effective thermal management solution for the combined traction drive system. This paper focuses on the thermal management approach selected for Oak Ridge National Laboratory’s outer-rotor-motor-based integrated traction drive and evaluates its potential performance. The outer-rotor-motor configuration provided an opportunity for integration of the six-phase inverter in the available space in the central cavity of the internal stator. A cylindrical inverter enclosure with integrated coolant (water-ethylene glycol) channels in its walls was designed to enable heat removal from the power electronics. Numerical thermal-fluid modeling and initial channel/fin optimization results for the cylindrical heat sink are presented here. As permanent magnets are integrated in a high-speed (20 000 RPM) outer rotor, forced air convection provides cooling for the magnets. The magnets were segmented axially to mitigate the eddy current losses. Heat generated in the stator windings and laminations is removed by water-ethylene glycol coolant circulating in interpolar T-shape ceramic heat exchangers inserted between windings. General design concepts and numerical thermal-fluid simulations illustrating the electric motor thermal management solution are also presented.
Kekelia, Bidzina↗
2003 NASA Faculty Fellowship Program at Glenn Research Center
The Office of Education at NASA Headquarters provides overall policy and direction for the NASA Faculty Fellowship Program (NFFP). The American Society for Engineering Education (ASEE) and the Universities Space Research Association (USRA) have joined in partnership to recruit participants, accept applications from a broad range of participants, and provide overall evaluation of the NFFP. The NASA Centers, through their University Affairs Officers, develop and operate the experiential part of the program. In concert with co-directing universities and the Centers, Fellows are selected and provided the actual research experiences. This report summarizes the 2003 session conducted at the Glenn Research Center (GRC).Research topics covered a variety of areas including, but not limited to, biological sensors, modeling of biological fluid systems, electronic circuits, ceramics and coatings, unsteady probablistic analysis and aerodynamics, gas turbines, environmental monitoring systems for water quality, air quality, gaseous and particulate emissions, bearings for flywheel energy storage, shape memory alloys,photonic interrogation and nanoprocesses,carbon nanotubes, polymer synthesis for fuel cells, aviation communications, algorithm development and RESPlan Database.
Prahl, Joseph M.↗
Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales
Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Understanding the Variability of Helium Abundance in the Solar Corona Using Three Fluid Modeling and Ultraviolet Observations
The variability of helium abundance in the solar corona and the solar wind is an important signature of solar activity, solar cycle, and solar wind sources, as well as coronal heating processes. Motivated by recently reported remote-sensing UV imaging observations by Helium Resonance Scattering in the Corona and Heliosphere payload sounding rocket of helium abundance in the inner corona on 2009 September 14 near solar minimum, we present the results of the first three-dimensional three-fluid (electrons, protons, and alpha particles) model of tilted coronal streamer belt and slow solar wind that illustrates the various processes leading to helium abundance differentiation and variability. We find good qualitative agreement between the three-fluid model and the coronal helium abundance variability deduced from UV observations of streamers, providing insight on the effects of the physical processes, such as heating, gravitational settling, and interspecies Coulomb friction in the outflowing solar wind that produce the observed features. The study impacts our understanding of the origins of the slow solar wind.
Leon Ofman↗
Saturn's Magnetosphere Interaction with Titan for T9 Encounter: 3D Hybrid Modeling and Comparison with CAPS Observations
Global dynamics of ionized and neutral gases in the environment of Titan plays an important role in the interaction of Saturn s magnetosphere with Titan. Several hybrid simulations of this problem have already been done (Brecht et al., 2000; Kallio et al., 2004; Modolo et al., 2007a; Simon et al., 2007a, 2007b; Modolo and Chanteur, 2008). Observational data from CAPS for the T9 encounter (Sittler et al., 2009) indicates an absence of O(+) heavy ions in the upstream that change the models of interaction which were discussed in current publications (Kallio et al., 2004; Modolo et al., 2007a; Simon et al., 2007a, 2007b; Ma et al., 2007; Szego et al., 2007). Further analysis of the CAPS data shows very low density or even an absence of H(+) ions in upstream. In this paper we discuss two models of the interaction of Saturn s magnetosphere with Titan: (A) high density of H(+) ions in the upstream flow (0.1/cu cm), and (B) low density of H(+) ions in the upstream flow (0.02/cu cm). The hybrid model employs a fluid description for electrons and neutrals, whereas a particle approach is used for ions. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. The model atmosphere includes exospheric H(+), H(2+), N(2+)and CH(4+) pickup ion production as well as an immobile background ionosphere and a shell distribution for active ionospheric ions (M(sub i)=28 amu). The hybrid model allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of Alfven wing-like structures. The results of the ion dynamics in Titan s environment are compared with Cassini T9 encounter data (CAPS).
Lipatov, A. S.↗
The formation of ion acoustic shocks
Recent experiments performed in the double plasma (DP) device have verified the existence of electrostatic ion acoustic laminar shocks. The influence of the piston on the shock structure is investigated by modeling the DP device and by numerically solving the temporal and spatial evolution of the shock. In order to isolate piston effects, as opposed to kinetic theory effects such as reflected ions and trapped electrons, the DP plasma is modeled as a cold ion fluid with isothermal Boltzmann electrons. It is shown that laminar shock transitions with structure agreeing with DP shock experiments can be excited.
White, R. B.↗
Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile
New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.
Kinsey, J. E. (ORCID:0000000193347473)↗
Nonlocal effects on thermal transport in hydrodynamic simulations of unmagnetized MagLIF-relevant gaspipes on NIF
We present simulations of heat flow relevant to gaspipe experiments on the National Ignition Facility to investigate kinetic effects on transport phenomena. D 2 and neopentane (C 5 H 12 ) filled targets are used to study the laser preheat stage of a MagLIF scheme where an axial magnetic field is sometimes applied to the target. Simulations were done with the radiation-MHD code HYDRA with a collision-dominated fluid model and the SNB nonlocal electron thermal conduction model. Using the SNB model to evolve the electron temperature increased the heat front propagation of neopentane gas targets compared to a local model by limiting radial heat flow. This increases electron temperature near the axis, which decreases laser absorption. We find that the effect of heat flow models on temperature profiles and laser propagation is modest. Beyond the SNB model, we utilize HYDRA to initialize plasma conditions for the Vlasov–Fokker–Planck K2 code. We run K2 until a quasi-steady state is reached and examine the impact of kinetic effects on heat transport. Although axial heat flow is well predicted by fluid models, the fluid model consistently overpredicts radial heat flow up to 150% in regions with the largest temperature gradient of D 2 filled gaspipes. On the other hand, the SNB nonlocal electron conduction model is found to be adequate for capturing kinetic heat flow in gaspipes.
Lau, Ryan Y. [Univ. of Colorado, Boulder, CO (Unit↗