Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTemperature gradient driven heat flux closure in fluid simulations of collisionless reconnection
68
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Recent efforts to include kinetic effects in fluid simulations of plasmas have been very promising. Concerning collisionless magnetic reconnection, it has been found before that damping of the pressure tensor to isotropy leads to good agreement with kinetic runs in certain scenarios. An accurate representation of kinetic effects in reconnection was achieved in a study by Wang et al. (Phys. Plasmas, volume 22, 2015, 012108) with a closure derived from earlier work by Hammett and Perkins (PRL, volume 64, 1990, 3019). Here, their approach is analyzed on the basis of heat flux data from a Vlasov simulation. As a result, we propose a new local closure in which heat flux is driven by temperature gradients. That way, a more realistic approximation of Landau damping in the collisionless regime is achieved. Previous issues are addressed and the agreement with kinetic simulations in different reconnection setups is improved significantly. To the authors knowledge, the new fluid model is the first to perform well in simulations of the coalescence of large magnetic islands.
rate research
Read More
In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment. In this paper, 2D fluid simulations are presented of this merging process in order to understand the underlying physics, and better interpret the experimental data. These simulations examine the individual and combined effects of tight-aspect ratio geometry and two-fluid physics on the merging. The ideal self-driven flux-rope dynamics are coupled to the diffusion layer physics, resulting in a large range of phenomena. For resistive MHD simulations, the flux-ropes enter the sloshing regime for normalised resistivity eta < 1E-5. In Hall-MHD three regimes are found for the qualitative behaviour of the current sheet, depending on the ratio of the current sheet width to the ion-sound radius. These are a stable collisional regime, an open X-point regime, and an intermediate regime that is highly unstable to tearing-type instabilities. In toroidal axisymmetric geometry, the final state after merging is a MAST-like spherical tokamak with nested flux-surfaces. It is also shown that the evolution of simulated 1D radial density profiles closely resembles the Thomson scattering electron density measurements in MAST. An intuitive explanation for the origin of the measured density structures is proposed, based upon the results of the toroidal Hall-MHD simulations.
The kinetic features of secondary magnetic reconnection in a single flux rope undergoing internal kink instability are studied by means of three-dimensional Particle-in-Cell simulations. Several signatures of secondary magnetic reconnection are identified in the plane perpendicular to the flux rope: a quadrupolar electron and ion density structure and a bipolar Hall magnetic field develop in proximity of the reconnection region. The most intense electric fields form perpendicularly to the local magnetic field, and a reconnection electric field is identified in the plane perpendicular to the flux rope. An electron current develops along the reconnection line in the opposite direction of the electron current supporting the flux rope magnetic field structure. Along the reconnection line, several bipolar structures of the electric field parallel to the magnetic field occur making the magnetic reconnection region turbulent. The reported signatures of secondary magnetic reconnection can help to localize magnetic reconnection events in space, astrophysical and fusion plasmas.
We investigate the linear theory of the ion-temperature-gradient (ITG) mode, with the goal of developing a general understanding that may be applied to stellarators. We highlight the Wendelstein 7X (W7-X) device. Simple fluid and kinetic models that follow closely from existing literature are reviewed and two new first-principle models are presented and compared with results from direct numerical simulation. One model investigates the effect of regions of strong localized shear, which are generic to stellarator equilibria. These shear spikes are found to have a potentially significant stabilizing affect on the mode; however, the effect is strongest at short wavelengths perpendicular to the magnetic field, and it is found to be significant only for the fastest growing modes in W7-X. A second model investigates the long-wavelength limit for the case of negligible global magnetic shear. The analytic calculation reveals that the effect of the curvature drive enters at second order in the drift frequency, confirming conventional wisdom that the ITG mode is slab-like at long wavelengths. Using flux tube simulations of a zero-shear W7-X configuration, we observe a close relationship to an axisymmetric configuration at a similar parameter point. It is concluded that scale lengths of the equilibrium gradients constitute a good parameter space to characterize the ITG mode. Thus, to optimize the magnetic geometry for ITG mode stability, it may be fruitful to focus on local parameters, such as the magnitude of bad curvature, connection length, and local shear at locations of bad curvature (where the ITG mode amplitude peaks).
Experimental data from the Mega Amp Spherical Tokamak (MAST) is used to show that the inverse gradient scale length of the ion temperature R/LTi (normalized to the major radius R) has its strongest local correlation with the rotational shear and the pitch angle of the magnetic field (or, equivalently, an inverse correlation with q/{epsilon}, the safety factor/the inverse aspect ratio). Furthermore, R/LTi is found to be inversely correlated with the gyro-Bohm-normalized local turbulent heat flux estimated from the density fluctuation level measured using a 2D Beam Emission Spectroscopy (BES) diagnostic. These results can be explained in terms of the conjecture that the turbulent system adjusts to keep R/LTi close to a certain critical value (marginal for the excitation of turbulence) determined by local equilibrium parameters (although not necessarily by linear stability).
Results of the first validation of large guide field, $B_g / delta B_0 gg 1$, gyrokinetic simulations of magnetic reconnection at a fusion and solar corona relevant $beta_i = 0.01$ and solar wind relevant $beta_i = 1$ are presented, where $delta B_0$ is the reconnecting field. Particle-in-cell (PIC) simulations scan a wide range of guide magnetic field strength to test for convergence to the gyrokinetic limit. The gyrokinetic simulations display a high degree of morphological symmetry, to which the PIC simulations converge when $beta_i B_g / delta B_0 gtrsim 1$ and $B_g / delta B_0 gg 1$. In the regime of convergence, the reconnection rate, relative energy conversion, and overall magnitudes are found to match well between the PIC and gyrokinetic simulations, implying that gyrokinetics is capable of making accurate predictions well outside its regime of formal applicability. These results imply that in the large guide field limit many quantities resulting from the nonlinear evolution of reconnection scale linearly with the guide field.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
