Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOverview of HiFi -- implicit spectral element code framework for multi-fluid plasma applications
77
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
An overview of the algorithm and a sampling of plasma applications of the implicit, adaptive high order finite (spectral) element modeling framework, HiFi, is presented. The distinguishing capabilities of the HiFi code include adaptive spectral element spatial representation with flexible geometry, highly parallelizable implicit time advance, and general flux-source form of the partial differential equations and boundary conditions that can be implemented in its framework. Early algorithm development and extensive verification studies of the two-dimensional version of the code, known as SEL, have been previously described [A.H. Glasser & X.Z. Tang, Comp. Phys. Comm., 164 (2004); V.S. Lukin, Ph.D. thesis, Princeton University (2008)]. Here, substantial algorithmic improvements and extensions are presented together with examples of two- and three- dimensional applications of the HiFi framework. These include a Cartesian two-dimensional incompressible magnetohydrodynamic simulation of low dissipation magnetic reconnection in a large system, a two-dimensional axisymmetric simulation of self-similar compression of a magnetic plasma confinement configuration, and a three-dimensional Hall MHD simulation of spheromak tilting and relaxation. Some planned efforts to further improve and expand the capabilities of the HiFi modeling framework are discussed.
rate research
Read More
A new modular code called BOUT++ is presented, which simulates 3D fluid equations in curvilinear coordinates. Although aimed at simulating Edge Localised Modes (ELMs) in tokamak X-point geometry, the code is able to simulate a wide range of fluid models (magnetised and unmagnetised) involving an arbitrary number of scalar and vector fields, in a wide range of geometries. Time evolution is fully implicit, and 3rd-order WENO schemes are implemented. Benchmarks are presented for linear and non-linear problems (the Orszag-Tang vortex) showing good agreement. Performance of the code is tested by scaling with problem size and processor number, showing efficient scaling to thousands of processors. Linear initial-value simulations of ELMs using reduced ideal MHD are presented, and the results compared to the ELITE linear MHD eigenvalue code. The resulting mode-structures and growth-rate are found to be in good agreement (BOUT++ = 0.245, ELITE = 0.239). To our knowledge, this is the first time dissipationless, initial-value simulations of ELMs have been successfully demonstrated.
We report on the algorithms and numerical methods used in Viriato, a novel fluid-kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations [Schekochihin et al., Astrophys. J. Suppl. 182:310 (2009)]. Two main applications of these equations are magnetised (Alfvenic) plasma turbulence and magnetic reconnection. Viriato uses operator splitting (Strang or Godunov) to separate the dynamics parallel and perpendicular to the ambient magnetic field (assumed strong). Along the magnetic field, Viriato allows for either a second-order accurate MacCormack method or, for higher accuracy, a spectral-like scheme composed of the combination of a total variation diminishing (TVD) third order Runge-Kutta method for the time derivative with a 7th order upwind scheme for the fluxes. Perpendicular to the field Viriato is pseudo-spectral, and the time integration is performed by means of an iterative predictor-corrector scheme. In addition, a distinctive feature of Viriato is its spectral representation of the parallel velocity-space dependence, achieved by means of a Hermite representation of the perturbed distribution function. A series of linear and nonlinear benchmarks and tests are presented, including a detailed analysis of 2D and 3D Orszag-Tang-type decaying turbulence, both in fluid and kinetic regimes.
The document describes a numerical algorithm to simulate plasmas and fluids in the 3 dimensional space by the Euler method, in which the spatial meshes are fixed to the space. The plasmas and fluids move through the spacial Euler mesh boundary. The Euler method can represent a large deformation of the plasmas and fluids. On the other hand, when the plasmas or fluids are compressed to a high density, the spatial resolution should be ensured to describe the density change precisely. The present 3D Euler code is developed to simulate a nuclear fusion fuel ignition and burning. Therefore, the 3D Euler code includes the DT fuel reactions, the alpha particle diffusion, the alpha particle deposition to heat the DT fuel and the DT fuel depletion by the DT reactions, as well as the thermal energy diffusion based on the three-temperature compressible fluid model.
Here we present details of a mixed implicit-explicit numerical scheme for the solution of the gyrokinetic-Poisson system of equations in the local limit. This scheme has been implemented in a new code called $texttt{stella}$, which is capable of evolving electrostatic fluctuations with full kinetic electron effects and an arbitrary number of ion species in general magnetic geometry. We demonstrate the advantages of this mixed approach over a fully explicit treatment and provide linear and nonlinear benchmark comparisons for both axisymmetric and non-axisymmetric magnetic equilibria.
We have developed an efficient algorithm for steady axisymmetrical 2D fluid equations. The algorithm employs multigrid method as well as standard implicit discretization schemes for systems of partial differential equations. Linearity of the multigrid method with respect to the number of grid points allowed us to use $256times 256$ grid, where we could achieve solutions in several minutes. Time limitations due to nonlinearity of the system are partially avoided by using multi level grids(the initial solution on $256times 256$ grid was extrapolated steady solution from $128times 128$ grid which allowed using long integration time steps). The fluid solver may be used as the basis for hybrid codes for DC discharges.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
