Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBORAY: An Axisymmetric Ray Tracing Code Supports Both Closed and Open Field Lines Plasmas
272
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Ray tracing codes are useful to study the electromagnetic wave propagation and absorption in the geometrical optics approximation. In magnetized fusion plasma community, most ray tracing codes assume the plasma density and temperature be functions of the magnetic flux and study waves only inside the last closed flux surface, which are sufficient for the present day tokamak. However, they are difficult to be used for configurations with open magnetic field line plasmas, such as mirror machine and field-reversed-configuration (FRC). We develop a ray tracing code in cylindrical coordinates $(r,phi,z)$ to support arbitrary axisymmetric configurations with both closed and open field lines plasmas. For wave propagation, the cold plasma dispersion relation is usually sufficient, and we require the magnetic field ${bf B}(r,z)$ and species densities $n_{s0}(r,z)$ profiles as input. For wave absorption, we require a further temperature $T_{s0}(r,z)$ profile to solve a hot kinetic plasma dispersion relation. In difference to other ray tracing codes which calculate the imaginary part of wave vector ${bf k}_{perp,i}$ for wave absorption, we calculate the imaginary part of wave frequency $omega_i$, which is shown to be equivalent with the former technique under weak damping approximation. The code can use either numerical or analytical equilibrium. Examples and benchmarks with electron cyclotron wave, lower hybrid wave and ion cyclotron wave for tokamak, spherical tokamak (ST), FRC and mirror machine are shown.
rate research
Read More
Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently has it 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. of 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 which 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, as well as 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.
Turbulence induced by the ion temperature gradient (ITG) is investigated in the helical and axisymmetric plasma states of a reversed field pinch device by means of gyrokinetic calculations. The two magnetic configurations are systematically compared, both linearly and nonlinearly, in order to evaluate the impact of the geometry on the instability and its ensuing transport, as well as on the production of zonal flows. Despite its enhanced confinement, the high-current helical state demonstrates a lower ITG stability threshold compared to the axisymmetric state, and ITG turbulence is expected to become an important contributor to the total heat transport.
The properties of the boundary plasma in a tokamak are now recognized to play a key role in determining the achievable fusion power and the lifetimes of plasma-facing components. Accurate quantitative modeling and improved qualitative understanding of the boundary plasma ultimately require five-dimensional gyrokinetic turbulence simulations, which have been successful in predicting turbulence and transport in the core. The additional challenges of boundary-plasma simulation necessitate the development of new gyrokinetic codes or major modifications to existing core gyrokinetic codes. In this thesis, we develop the first gyrokinetic continuum code capable of simulating plasma turbulence on open magnetic field lines, which is a key feature of a tokamak scrape-off layer. In contrast to prior attempts at this problem, we use an energy-conserving discontinuous Galerkin discretization in space. To model the interaction between the plasma and the wall, we design conducting-sheath boundary conditions that permit local currents into and out of the wall. We start by designing spatially one-dimensional kinetic models of parallel SOL dynamics and solve these systems using novel continuum algorithms. By generalizing these algorithms to higher dimensions and adding a model for collisions, we present results from the first gyrokinetic continuum simulations of turbulence on two types of open-field-line systems. The first simulation features uniform and straight field lines, such as found in some linear plasma devices. The second simulation is of a hypothetical model we developed of the NSTX scrape-off layer featuring helical field lines. These developments comprise a major step towards a gyrokinetic continuum code for quantitative predictions of turbulence and transport in the boundary plasma of magnetic fusion devices.
An axisymmetric fully ionized plasma rotates around its axis when a charge separation between magnetic surfaces is produced from DC fields or RF waves. On each magnetic surface both electrons and ions obey the isorotation law and perform an azimuthal E cross B rotation at the same angular velocity. When Coulomb collisions are taken into account such a flow displays no Ohmic current short circuiting of the charge separation and thus no linear dissipation. A nonlinear Ohmic response appears when inertial effects are considered, providing a dissipative relaxation of the charge separation between the magnetic surfaces. This nonlinear conductivity results from an interplay between Coriolis, centrifugal and electron-ion collisional friction forces. This phenomena is identified, described and analyzed. In addition, both the quality factor of angular momentum storage as well as the efficiency of wave driven angular momentum generation are calculated and shown to be independent of the details of the charge separation processes.
It is demonstrated that finite-pressure, approximately quasi-axisymmetric stellarator equilibria can be directly constructed (without numerical optimization) via perturbations of given axisymmetric equilibria. The size of such perturbations is measured in two ways, via the fractional external rotation and, alternatively, via the relative magnetic field strength, i.e. the average size of the perturbed magnetic field, divided by the unperturbed field strength. It is found that significant fractional external rotational transform can be generated by quasi-axisymmetric perturbations, with a similar value of the relative field strength, despite the fact that the former scales more weakly with the perturbation size. High mode number perturbations are identified as a candidate for generating such transform with local current distributions. Implications for the development of a general non-perturbative solver for optimal stellarator equilibria is discussed.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
