Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStepped pressure equilibrium with relaxed flow and applications in reversed-field pinch plasmas
61
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The Multi-region Relaxed MHD (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The Stepped Pressure Equilibrium Code (SPEC) [S.R. Hudson et al., Phys. Plasmas 19, 112502 (2012)] was developed to calculate MRxMHD equilibria numerically. In this work, we have extended the SPEC code to compute MRxMHD equilibria with field-aligned flow and rotation, following the theoretical development to incorporate cross-helicity and angular momentum constraints. The code has been verified for convergence and compared to a Grad-Shafranov solver in 2D. We apply our new tool to study the flow profile change before and after the sawtooth crash of a reversed-field pinch discharge, in which data of the parallel flow is available. We find the promising result that under the constraints of cross-helicity and angular momentum, the parallel flow profile in post-crash SPEC equilibrium is flat in the plasma core and the amplitude of the flow matches experimental observations. Finally, we provide an example equilibrium with a 3D helical field structure as the favoured lower energy state. This will be the first 3D numerical equilibrium in which the flow effects are self-consistently calculated.
rate research
Read More
In the reversed field pinch RFX-mod strong electron temperature gradients develop when the Single-Helical-Axis regime is achieved. Gyrokinetic calculations show that in the region of the strong temperature gradients microtearing instabilities are the dominant turbulent mechanism acting on the ion Larmor radius scale. The quasi-linear evaluation of the electron thermal conductivity is in good agreement with the experimental estimates.
A new equilibrium pressure profile extending the Rigid-Rotor (RR) model with a simple unified expression $P=P(psi;beta_{s},alpha, sigma)$ for both inside and outside the separatrix is proposed, in which the radial normalized field-reversed configuration (FRC) equilibrium profiles for pressure, magnetic field, and current can be determined by only two dimensionless parameters $beta_sequiv P_s/2mu_0B_e^2$ and $delta_sequiv L_{ps}/R_s$, where $P_s$ is the thermal pressure at the separatrix, $B_e$ is the external magnetic field strength, $L_{ps}$ is the pressure profile scale length at the separatrix, and $R_s$ is the separatrix radius. This modified rigid rotor (MRR) model has sufficient flexibility to accommodate the narrow scrape of layer (SOL) width and hollow current density profiles, and can be used to fit experimental measurements satisfactorily. Detailed one-dimensional (1D) characteristics of the new MRR model are investigated analytically and numerically, and the results are also confirmed in two-dimensional (2D) numerical equilibrium solutions.
Till now the magnetohydrodynamic (MHD) simulation of the reversed field pinch (RFP) has been performed by assuming axis-symmetric radial time independent dissipation profiles. In helical states this assumption is not correct since these dissipations should be flux functions, and should exhibit a helical symmetry as well. Therefore more correct simulations should incorporate self-consistent dissipation profiles. As a first step in this direction, the case of uniform dissipation profiles was considered by using the 3D nonlinear visco-resistive MHD code SpeCyl. It is found that a flattening of the resistivity profile results in the reduction of the dynamo action, which brings to marginally-reversed or even non-reversed equilibrium solutions. The physical origin of this result is discussed in relation to the electrostatic drift explanation of the RFP dynamo. This sets constraints on the functional choice of dissipations in future self-consistent simulations.
A new tool (GSEQ-FRC) for solving two-dimensional (2D) equilibrium of field-reversed configuration (FRC) based on fixed boundary and free boundary conditions with external coils included is developed. Benefiting from the two-parameter modified rigid rotor (MRR) radial equilibrium model and the numerical approaches presented by [Ma et al, Nucl. Fusion, 61, 036046, 2021], GSEQ-FRC are used to study the equilibrium properties of FRC quantitatively and will be used for fast FRC equilibrium reconstruction. In GSEQ-FRC, the FRC equilibrium can be conveniently determined by two parameters, i.e., the ratio between thermal pressure and magnetic pressure at the seperatrix $beta_s$, and the normalized scrape of layer (SOL) width $delta_s$. Examples with fixed and free boundary conditions are given to demonstrate the capability of GSEQ-FRC in the equilibrium calculations. This new tool is used to quantitatively study the factors affecting the shape of the FRC separatrix, revealing how the FRC changes from racetrack-like to ellipse-like.
The Stepped Pressure Equilibrium Code (SPEC) [Hudson et al., Phys. Plasmas 19, 112502 (2012)] has been successful in the construction of equilibria in 3D configurations that contain a mixture of flux surfaces, islands and chaotic magnetic field lines. In this model, the plasma is sliced into sub-volumes separated by ideal interfaces, and in each volume the magnetic field is a Beltrami field. In the cases where the system is far from possessing a continuous symmetry, such as in stellarators, the existence of solutions to a stepped-pressure equilibrium with given constraints, such as a multi-region relaxed MHD minimum energy state, is not guaranteed but is often taken for granted. Using SPEC, we have studied two different scenarios in which a solution fails to exist in a slab with analytic boundary perturbations. We found that with a large boundary perturbation, a certain interface becomes fractal, corresponding to the break up of a Kolmogorov-Arnold-Moser (KAM) surface. Moreover, an interface can only support a maximum pressure jump while a solution of the magnetic field consistent with the force balance condition can be found. An interface closer to break-up can support a smaller pressure jump. We discovered that the pressure jump can push the interface closer to being non-smooth through force balance, thus significantly decreasing the maximum pressure it can support. Our work shows that a convergence study must be performed on a SPEC equilibrium with interfaces close to break-up. These results may also provide insights into the choice of interfaces and have applications in finding out the maximum pressure a machine can support.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
