scholarly journals Development of a Gyrokinetic Particle-in-Cell Code for Whole-Volume Modeling of Stellarators

Plasma ◽  
2019 ◽  
Vol 2 (2) ◽  
pp. 179-200 ◽  
Author(s):  
Toseo Moritaka ◽  
Robert Hager ◽  
Michael Cole ◽  
Samuel Lazerson ◽  
Choong-Seock Chang ◽  
...  

We present initial results in the development of a gyrokinetic particle-in-cell code for the whole-volume modeling of stellarators. This is achieved through two modifications to the X-point Gyrokinetic Code (XGC), originally developed for tokamaks. One is an extension to three-dimensional geometries with an interface to Variational Moments Equilibrium Code (VMEC) data. The other is a connection between core and edge regions that have quite different field-line structures. The VMEC equilibrium is smoothly extended to the edge region by using a virtual casing method. Non-axisymmetric triangular meshes in which triangle nodes follow magnetic field lines in the toroidal direction are generated for field calculation using a finite-element method in the entire region of the extended VMEC equilibrium. These schemes are validated by basic benchmark tests relevant to each part of the calculation cycle, that is, particle push, particle-mesh interpolation, and field solver in a magnetic field equilibrium of Large Helical Device including the edge region. The developed code also demonstrates collisionless damping of geodesic acoustic modes and steady states with residual zonal flow in the core region.

2020 ◽  
Vol 86 (3) ◽  
Author(s):  
A. Biancalani ◽  
A. Bottino ◽  
P. Lauber ◽  
A. Mishchenko ◽  
F. Vannini

Numerical simulations of Alfvén modes driven by energetic particles are performed with the gyrokinetic (GK) global particle-in-cell code ORB5. A reversed shear equilibrium magnetic field is adopted. A simplified configuration with circular flux surfaces and large aspect ratio is considered. The nonlinear saturation of beta-induced Alfvén eigenmodes (BAE) is investigated. The roles of the wave–particle nonlinearity of the different species, i.e. thermal ions, electrons and energetic ions are described, in particular for their role in the saturation of the BAE and the generation of zonal flows. The nonlinear redistribution of the electron population is found to be important in increasing the BAE saturation level and the zonal flow amplitude.


2009 ◽  
Vol 80 (12) ◽  
pp. 123501 ◽  
Author(s):  
Peter Drewelow ◽  
Torsten Bräuer ◽  
Matthias Otte ◽  
Friedrich Wagner ◽  
Andreas Werner

2009 ◽  
Vol 5 (H15) ◽  
pp. 434-435
Author(s):  
A. Lazarian ◽  
G. Kowal ◽  
E. Vishniac ◽  
K. Kulpa-Dubel ◽  
K. Otmianowska-Mazur

AbstractA magnetic field embedded in a perfectly conducting fluid preserves its topology for all times. Although ionized astrophysical objects, like stars and galactic disks, are almost perfectly conducting, they show indications of changes in topology, magnetic reconnection, on dynamical time scales. Reconnection can be observed directly in the solar corona, but can also be inferred from the existence of large scale dynamo activity inside stellar interiors. Solar flares and gamma ray busts are usually associated with magnetic reconnection. Previous work has concentrated on showing how reconnection can be rapid in plasmas with very small collision rates. Here we present numerical evidence, based on three dimensional simulations, that reconnection in a turbulent fluid occurs at a speed comparable to the rms velocity of the turbulence, regardless of the value of the resistivity. In particular, this is true for turbulent pressures much weaker than the magnetic field pressure so that the magnetic field lines are only slightly bent by the turbulence. These results are consistent with the proposal by Lazarian & Vishniac (1999) that reconnection is controlled by the stochastic diffusion of magnetic field lines, which produces a broad outflow of plasma from the reconnection zone. This work implies that reconnection in a turbulent fluid typically takes place in approximately a single eddy turnover time, with broad implications for dynamo activity and particle acceleration throughout the universe. In contrast, the reconnection in 2D configurations in the presence of turbulence depends on resistivity, i.e. is slow.


Author(s):  
Andrew L Haynes ◽  
Clare E Parnell ◽  
Klaus Galsgaard ◽  
Eric R Priest

The heating of the solar corona is probably due to reconnection of the highly complex magnetic field that threads throughout its volume. We have run a numerical experiment of an elementary interaction between the magnetic field of two photospheric sources in an overlying field that represents a fundamental building block of the coronal heating process. The key to explaining where, how and how much energy is released during such an interaction is to calculate the resulting evolution of the magnetic skeleton. A skeleton is essentially the web of magnetic flux surfaces (called separatrix surfaces) that separate the coronal volume into topologically distinct parts. For the first time, the skeleton of the magnetic field in a three-dimensional numerical magnetohydrodynamic experiment is calculated and carefully analysed, as are the ways in which it bifurcates into different topologies. A change in topology normally changes the number of magnetic reconnection sites. In our experiment, the magnetic field evolves through a total of six distinct topologies. Initially, no magnetic flux joins the two sources. Then, a new type of bifurcation, called a global double-separator bifurcation , takes place. This bifurcation is probably one of the main ways in which new separators are created in the corona (separators are field lines at which three-dimensional reconnection takes place). This is the first of five bifurcations in which the skeleton becomes progressively more complex before simplifying. Surprisingly, for such a simple initial state, at the peak of complexity there are five separators and eight flux domains present.


Author(s):  
Kyung Sun Park

We performed high-resolution three-dimensional global MHD simulations to determine the impact of weak southward interplanetary magnetic field (IMF) (Bz = −2 nT) and slow solar wind to the Earth’s magnetosphere and ionosphere. We considered two cases of differing, uniform time resolution with the same grid spacing simulation to find any possible differences in the simulation results. The simulation results show that dayside magnetic reconnection and tail reconnection continuously occur even during the weak and steady southward IMF conditions. A plasmoid is generated on closed plasma sheet field lines. Vortices are formed in the inner side of the magnetopause due to the viscous-like interaction, which is strengthened by dayside magnetic reconnection. We estimated the dayside magnetic reconnection which occurred in relation to the electric field at the magnetopause and confirmed that the enhanced electric field is caused by the reconnection and the twisted structure of the electric field is due to the vortex. The simulation results of the magnetic field and the plasma properties show quasi-periodic variations with a period of 9–11 min between the appearances of vortices. Also the peak values of the cross-polar cap potential are both approximately 50 kV, the occurrence time of dayside reconnections are the same, and the polar cap potential patterns are the same in both cases. Thus, there are no significant differences in outcome between the two cases.


2016 ◽  
Vol 82 (5) ◽  
Author(s):  
Z. Akbari ◽  
M. Hosseinpour ◽  
M. A. Mohammadi

In a three-dimensional non-null magnetic reconnection, the process of magnetic reconnection takes place in the absence of a null point where the magnetic field vanishes. By randomly injecting a population of 10 000 protons, the trajectory and energy distribution of accelerated protons are investigated in the presence of magnetic and electric fields of a particular model of non-null magnetic reconnection with the typical parameters for the solar corona. The results show that protons are accelerated along the magnetic field lines away from the non-null point only at azimuthal angles where the magnitude of the electric field is strongest and therefore particles obtain kinetic energies of the order of thousands of MeV and even higher. Moreover, the energy distribution of the population depends strongly on the amplitude of the electric and magnetic fields. Comparison shows that a non-null magnetic reconnection is more efficient in accelerating protons to very high GeV energies than a null-point reconnection.


2008 ◽  
Vol 26 (8) ◽  
pp. 2471-2483 ◽  
Author(s):  
K. G. Tanaka ◽  
A. Retinò ◽  
Y. Asano ◽  
M. Fujimoto ◽  
I. Shinohara ◽  
...  

Abstract. The magnetopause (MP) reconnection is characterized by a density asymmetry across the current sheet. The asymmetry is expected to produce characteristic features in the reconnection layer. Here we present a comparison between the Cluster MP crossing reported by Retinò et al. (2006) and virtual observations in two-dimensional particle-in-cell simulation results. The simulation, which includes the density asymmetry but has zero guide field in the initial condition, has reproduced well the observed features as follows: (1) The prominent density dip region is detected at the separatrix region (SR) on the magnetospheric (MSP) side of the MP. (2) The intense electric field normal to the MP is pointing to the center of the MP at the location where the density dip is detected. (3) The ion bulk outflow due to the magnetic reconnection is seen to be biased towards the MSP side. (4) The out-of-plane magnetic field (the Hall magnetic field) has bipolar rather than quadrupolar structure, the latter of which is seen for a density symmetric case. The simulation also showed rich electron dynamics (formation of field-aligned beams) in the proximity of the separatrices, which was not fully resolved in the observations. Stepping beyond the simulation-observation comparison, we have also analyzed the electron acceleration and the field line structure in the simulation results. It is found that the bipolar Hall magnetic field structure is produced by the substantial drift of the reconnected field lines at the MSP SR due to the enhanced normal electric field. The field-aligned electrons at the same MSP SR are identified as the gun smokes of the electron acceleration in the close proximity of the X-line. We have also analyzed the X-line structure obtained in the simulation to find that the density asymmetry leads to a steep density gradient in the in-flow region, which may lead to a non-stationary behavior of the X-line when three-dimensional freedom is taken into account.


1995 ◽  
Vol 299 ◽  
pp. 153-186 ◽  
Author(s):  
P. A. Davidson

It is well known that the imposition of a static magnetic field tends to suppress motion in an electrically conducting liquid. Here we look at the magnetic damping of liquid-mental flows where the Reynolds number is large and the magnetic Reynolds number is small. The magnetic field is taken as uniform and the fluid is either infinite in extent or else bounded by an electrically insulating surface S. Under these conditions, we find that three general principles govern the flow. First, the Lorentz force destroys kinetic energy but does not alter the net linear momentum of the fluid, nor does it change the component of angular momentum parallel to B. In certain flows, this implies that momentum, linear or angular, is conserved. Second, the Lorentz force guides the flow in such a way that the global Joule dissipation, D, decreases, and this decline in D is even more rapid than the corresponding fall in global kinetic energy, E. (Note that both D and E are quadratic in u). Third, this decline in relative dissipation, D / E, is essential to conserving momentum, and is achieved by propagating linear or angular momentum out along the magnetic field lines. In fact, this spreading of momentum along the B-lines is a diffusive process, familiar in the context of MHD turbulence. We illustrate these three principles with the aid of a number of specific examples. In increasing order of complexity we look at a spatially uniform jet evolving in time, a three-dimensional jet evolving in space, and an axisymmetric vortex evolving in both space and time. We start with a spatially uniform jet which is dissipated by the sudden application of a transverse magnetic field. This simple (perhaps even trivial) example provides a clear illustration of our three general principles. It also provides a useful stepping-stone to our second example of a steady three-dimensional jet evolving in space. Unlike the two-dimensional jets studied by previous investigators, a three-dimensional jet cannot be annihilated by magnetic braking. Rather, its cross-section deforms in such a way that the momentum flux of the jet is conserved, despite a continual decline in its energy flux. We conclude with a discussion of magnetic damping of axisymmetric vortices. As with the jet flows, the Lorentz force cannot destroy the motion, but rather rearranges the angular momentum of the flow so as to reduce the global kinetic energy. This process ceases, and the flow reaches a steady state, only when the angular momentum is uniform in the direction of the field lines. This is closely related to the tendency of magnetic fields to promote two-dimensional turbulence.


Sign in / Sign up

Export Citation Format

Share Document