Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA SIMULATION”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A resistive MHD model and simulation on plasma flow evolution in the presence of resonant magnetic perturbation in a tokamak

Nonaxisymmetric magnetic fields such as the intrinsic error field and the externally applied resonant magnetic perturbation (RMP) in a tokamak are known to influence the plasma momentum transport and flow evolution through plasma response, which itself strongly depends on the plasma flow as well. The nonlinear interaction between plasma response and flow has been previously modeled in the conventional error field theory with the “no-slip” condition, which has been recently extended to allow the “free-slip” condition. In this work, we further target this specific process and numerically simulate the nonlinear plasma response and flow evolution in the presence of a single-helicity RMP in a circular-shaped model tokamak configuration, based on the full resistive MHD model in the initial-value code NIMROD. Time evolution of the parallel (to k) flow or “slip frequency” profile and its asymptotic steady state obtained from the NIMROD simulations are compared with both conventional and extended nonlinear response theories. Here, k is the wave vector of the propagating island. Good agreement with the extended theory with free-slip condition has been achieved for the parallel flow profile evolution in response to RMP in all resistive regimes, whereas the difference from the conventional theory with the no-slip condition tends to diminish as the plasma resistivity approaches zero.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

INTERMITTENT HEATING IN SOLAR WIND AND KINETIC SIMULATIONS

Low-density astrophysical plasmas may be described by magnetohydrodynamics at large scales, but require kinetic description at ion scales in order to include dissipative processes that terminate the cascade. Here kinetic plasma simulations and high-resolution spacecraft observations are compared to facilitate the interpretation of signatures of various dissipation mechanisms. Kurtosis of increments indicates that kinetic scale coherent structures are present, with some suggestion of incoherent activity near ion scales. Conditioned proton temperature distributions suggest heating associated with coherent structures. The results reinforce the association of intermittent turbulence, coherent structures, and plasma dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Impact of screening of resonant magnetic perturbations in three dimensional edge plasma transport simulations for DIII-D

The impact of resonant magnetic perturbations (RMPs) on the plasma edge can be analyzed in detail by three dimensional computer simulations, which take the underlying magnetic field structure as input. Previously, the 'vacuum approximation' has been used to calculate the magnetic field structure although plasma response effects may result in a screening (or even an amplification) of the external perturbations. Simulation results for an ITER similar shape plasma at the DIII-D tokamak are presented for the full vacuum perturbation field and an ad hoc screening case in comparison to the unperturbed configuration. It is shown that the RMP induced helical patterns in the plasma edge and on the divertor target shrink once screening is taken into account. However, a flat temperature profile is still found in the 'open field line domain' inside the separatrix, while the 'density pump out effect' found in the vacuum RMP case is considerably weakened.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Flux tube train model for local turbulence simulation of toroidal plasmas

A new simulation method for local turbulence in toroidal plasmas is developed by extending the conventional idea of the flux tube model. In the new approach, a train of flux tubes is employed, where flux tube simulation boxes are serially connected at each end along a field line so as to preserve a symmetry of the local gyrokinetic equations for image modes in an axisymmetric torus. Validity of the flux tube train model is confirmed against the toroidal ion temperature gradient turbulence for a case with a long parallel correlation of fluctuations, demonstrating numerical advantages over the conventional method in the time step size and the symmetry-preserving property.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulations of plasma current induced by toroidal field ramping down on tokamaks

In tokamak experiments, the BT ramping down can drive off-axis and parallel inductive current. This approach leads to a decrease in q95 and an increase in normalized beta, βN. The off-axis inductive current also broadens the current profile. Typically, the 1.5D transport code is used to simulate the time evolution of the plasma current profile, which is based on the flux-surface-averaged Faraday’s law. The ONETWO code is one of those transport codes. However, this code cannot simulate the situation of the evolving toroidal field BT. This article proposes a new time-dependent model to take into account the BT ramping down situation. A modified formula of flux-surface-averaged Faraday’s law was derived to consider the effect of the BT ramping rate on the current evolution, and it was implemented in the ONETWO code. Then, the modified ONETWO code was used to simulate the BT ramping down experiment on DIII-D. The simulation result of the plasma current evolution with BT ramping down shows a broader current profile with smaller ohmic current induced by the poloidal field, compared to that without BT ramping down.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simulation of the dynamics of gas mixtures during plasma processing in the C75 Cavity

Plasma processing using a mixture of noble gas and oxygen is a technique that is currently being used to reduce field emission and multipacting in accelerating cavities. Plasma is created inside the cavity when the gas mixture is exposed to an electromagnetic field that is generated by applying RF power through the fundamental power or higher-order mode couplers. Oxygen ions and atomic oxygen are created in the plasma which breaks down the hydrocarbons on the surface of the cavity and the residuals from this process are removed as part of the process gas flow. Removal of hydrocarbons from the surface increases the work function and reduces the secondary emission coefficient. This work describes the initial results of plasma simulation, which provides insight into the ignition process, distribution of different species, and interactions of free oxygen and oxygen ions with the cavity surfaces. The simulations have been done with an Ar/?2 plasma using COMSOL® multiphysics. These simulations help in understanding the dynamics and control of plasma inside the cavity and the exploration of different gas mixtures.

Dhakal, Pashupati↗
DOE OSTI2016

Finite time step and spatial grid effects in δf simulation of warm plasmas

This paper introduces a technique for analyzing time integration methods used with the particle weight equations in δf method particle-in-cell (PIC) schemes. The analysis applies to the simulation of warm, uniform, periodic or infinite plasmas in the linear regime and considers the collective behavior similar to the analysis performed by Langdon for full-f PIC schemes [1,2]. We perform both a time integration analysis and spatial grid analysis for a kinetic ion, adiabatic electron model of ion acoustic waves. An implicit time integration scheme is studied in detail for δf simulations using our weight equation analysis and for full-f simulations using the method of Langdon. It is found that the δf method exhibits a CFL-like stability condition for low temperature ions, which is independent of the parameter characterizing the implicitness of the scheme. The accuracy of the real frequency and damping rate due to the discrete time and spatial schemes is also derived using a perturbative method. The theoretical analysis of numerical error presented here may be useful for the verification of simulations and for providing intuition for the design of new implicit time integration schemes for the δf method, as well as understanding differences between δf and full-f approaches to plasma simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

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↗
DOE OSTI2024

Impact of Variable Perpendicular Transport Coefficients in WEST Simulations Using SolEdge-HDG

Plasma–wall interaction is one of the key research topics on the way to controlled fusion. To study the best operational designs with reduced heat and particle fluxes onto tokamak plasma facing components (PFCs) comprehensive plasma simulations are required. A recent implementation of a hybridized discontinuous Galerkin scheme into a new version of SolEdge code has the advantage of using magnetic equilibrium-free mesh. This allows us to conduct pioneering 2-D transport simulations of a full discharge in the WEST tokamak. In this work, we implemented plasma transport coefficients as functions of coordinate in the poloidal plane and neutral diffusion as a function of neutral mean free path. Moreover, the perpendicular convection flux terms were added to the code. Using the new features, a few test cases were investigated. Finally, the influence of nonconstant transport coefficients on the simulated particle and heat fluxes onto the WEST tokamak PFCs are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Extreme Scale Plasma Turbulence Simulations on Top Supercomputers Worldwide

The goal of the extreme scale plasma turbulence studies described in this paper is to expedite the delivery of reliable predictions on confinement physics in large magnetic fusion systems by using world-class supercomputers to carry out simulations with unprecedented resolution and temporal duration. This has involved architecture-dependent optimizations of performance scaling and addressing code portability and energy issues, with the metrics for multi-platform comparisons being 'time-to-solution' and 'energy-to-solution'. Realistic results addressing how confinement losses caused by plasma turbulence scale from present-day devices to the much larger $25 billion international ITER fusion facility have been enabled by innovative advances in the GTC-P code including (i) implementation of one-sided communication from MPI 3.0 standard; (ii) creative optimization techniques on Xeon Phi processors; and (iii) development of a novel performance model for the key kernels of the PIC code. Our results show that modeling data movement is sufficient to predict performance on modern supercomputer platforms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Simulations of plasma behavior during pellet injection in ITER

Plasma behavior during pellet injection in ITER is investigated using a 1.5D BALDUR integrated predictive modeling code. In these simulations, the pellet ablation is described using the neutral gas shielding (NGS) model developed by Parks and Turnbull [Phys. Fluids 21, 1735 (1978)]. The NGS pellet ablation model that includes the {nabla}B drift effect is coupled with a plasma core transport model, which is a combination of an MMM95 anomalous transport model and an NCLASS neoclassical transport model. The combination of core transport models, together with pellet model, is used to simulate the time evolution of plasma current, ion and electron temperatures, and density profiles for ITER standard type-I ELMy H-mode discharges during the pellet injection. It is found that the injection of pellet can result in either enhancement or degradation of plasma performance. The {nabla}B drift effect on the pellet deposition is very strong in ITER. The plasma density with high field side pellets, which favorable with the {nabla}B drift effect, is much higher and pellet can penetrate much deeper than that with low field side pellets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulation of plasma and neutral transport in PISCES-RF using SOLPS-ITER

In this research, the fluid plasma transport code SOLPS-ITER is applied and validated against experimental data from the plasma interaction surface component experimental station (PISCES)-RF linear plasma device to establish a physics basis for plasma and neutral transport in its two magnetic field (B-field) geometry setups-(1)the cusp and (2) non-cusp or linear B-field. The main focus of this study is to understand (1) radial plasma transport (2) heat and particle loads on the upstream dump and downstream target plate, and (3) the physics of plasma-neutral interactions in PISCES-RF. The simulation setup adheres to typical PISCES-RF experimental conditions, with a 2D helicon power deposition profile as an input heating source. SOLPS-ITER simulations reproduce experimental conditions with the Bohm diffusion model for both B-field configurations of the PISCES-RF experiment. Major energy loss channels include neutral radiation and power deposited on the wall and dump plate, with only 1% of the input power reaching the target. The ionization front is well confined near the dump plate due to the heating and puffing regions. Additionally, SOLPS-ITER simulation results are also found to be in very good agreement with the particle-in-cell calculations using the code-PICOS++ which supports the validity of SOLPS in low collisionality regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Simulation of neoclassical transport with the continuum gyrokinetic code COGENT

The development of the continuum gyrokinetic code COGENT for edge plasma simulations is reported. The present version of the code models a nonlinear axisymmetric 4D (R, v∥, μ) gyrokinetic equation coupled to the long-wavelength limit of the gyro-Poisson equation. Here, R is the particle gyrocenter coordinate in the poloidal plane, and v∥ and μ are the guiding center velocity parallel to the magnetic field and the magnetic moment, respectively. The COGENT code utilizes a fourth-order finite-volume (conservative) discretization combined with arbitrary mapped multiblock grid technology (nearly field-aligned on blocks) to handle the complexity of tokamak divertor geometry with high accuracy. Furthermore, topics presented are the implementation of increasingly detailed model collision operators, and the results of neoclassical transport simulations including the effects of a strong radial electric field characteristic of a tokamak pedestal under H-mode conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

NUMERICAL SIMULATIONS OF PLASMA DYNAMICS IN THE VICINITY OF A RETRACTING FLUX TUBE

In a previous paper, we presented an analytical, zero- β model for supra-arcade downflows in which a retracting flux tube deforms the surrounding magnetic field, constricting the flow of plasma along affected field lines and, in some cases, forcing the plasma to exhibit collimated shocks. Here we present a numerical simulation based on the same model construction—a retracting flux tube is treated as a rigid boundary around which the plasma is forced to flow and the magnetic field and plasma evolve according to the governing equations of magnetohydrodynamics. We find that the collimated shocks described in our previous study are recovered for plasma β in the range of 0 ≤ β ≲ 1, while for 1 ≲ β the behavior is similar to the simpler hydrodynamic case, with classical bow shocks forming when the acoustic Mach number approaches or exceeds unity. Furthermore, we find that while the plasma β is important for identifying the various types of behaviors, more important still is the Alfvén Mach number, which, if large, implies that the bulk kinetic energy of the fluid exceeds the internal energy of the magnetic field, thereby leading to the formation of unconfined, fast-mode magnetosonic shocks, even in the limit of small β .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Fast low-to-high confinement mode bifurcation dynamics in a tokamak edge plasma gyrokinetic simulation

Transport barrier formation and its relation to sheared flows in fluids and plasmas are of fundamental interest in various natural and laboratory observations and of critical importance in achieving an economical energy production in a magnetic fusion device. Here we report the first observation of an edge transport barrier formation event in an electrostatic gyrokinetic simulation carried out in a realistic diverted tokamak edge geometry under strong forcing by a high rate of heat deposition. Here, the results show that turbulent Reynolds-stress-driven sheared E x B flows act in concert with neoclassical orbit loss to quench turbulent transport and form a transport barrier just inside the last closed magnetic flux surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Three dimensional simulations of plasma filaments in the scrape off layer: A comparison with models of reduced dimensionality

This paper presents simulations of isolated 3D filaments in a slab geometry obtained using a newly developed 3D reduced fluid code, written using the BOUT++ framework. First, systematic scans were performed to investigate how the dynamics of a filament are affected by its amplitude, perpendicular size, and parallel extent. The perpendicular size of the filament was found to have a strong influence on its motions, as it determined the relative importance of parallel currents to polarization and viscous currents, whilst drift-wave instabilities were observed if the initial amplitude of the blob was increased sufficiently. Next, the 3D simulations were compared to 2D simulations using different parallel closures; namely, the sheath dissipation closure, which neglects parallel gradients, and the vorticity advection closure, which neglects the influence of parallel currents. The vorticity advection closure was found to not replicate the 3D perpendicular dynamics and overestimated the initial radial acceleration of all the filaments studied. In contrast, a more satisfactory comparison with the sheath dissipation closure was obtained, even in the presence of significant parallel gradients, where the closure is no longer valid. Specifically, it captured the contrasting dynamics of filaments with different perpendicular sizes that were observed in the 3D simulations which the vorticity advection closure failed to replicate. However, neither closure successfully replicated the Boltzmann spinning effects nor the associated poloidal drift of the blob that was observed in the 3D simulations. Although the sheath dissipation closure was concluded to be more successful in replicating the 3D dynamics, it is emphasized that the vorticity advection closure may still be relevant for situations where the parallel current is inhibited from closing through the sheath due to effects such as strong magnetic shear around X points or increased resistivity near the targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗