LJKDeterministic Models and Algorithms: EDPMOISE Seminar

On Thursday January 23 2014 at 11h00 in waiting for a room

Seminary of Mme Aurore BACK (Centre de Physique ThÃ©orique, Marseille)

A generalization of the magnetic geometry in GYSELA code : A 4D VlasovPoisson solver an arbitrary curvilinear grid

Summary

The large magnetic field in a Tokamak generates a huge anisotropy in the physics along and across magnetic field lines. For this reason aligning the grid on the magnetic surfaces and possibly on the magnetic field lines can considerably increase the accuracy for a given resolution. For Tokamaks with a circular poloidal cross section, magnetic surfaces are circular and thus standard toroidal coordinates are naturally aligned on magnetic surfaces. This is not the case for more general equilibria with an Xpoint, where a numerical definition of the mesh is required. For this reason, we need a Vlasov solver that can handle such a mesh.
The method of CAODAO can help us describe the geometry of magnetic surfaces including the Xpoint. It will give us a mesh of the poloidal plane as a spline surface. This is defined as a small collection of patches where a logical grid is defined and mapped to the actual computational domain. Each patch defines a curvilinear grid. We solve all equations in the patch using the coordinate transformation defined by the mapping. So we must develop a numerical method for the Vlasov equation in four dimensional phasespace that can handle this change of coordinates. Note that the velocity grid remains cartesian.
The semiLagrangian method consists in two steps: 1) Following the characteristics, which are particle trajectories, originating from grid points, 2) Interpolating back on the grid. The second step needs to be performed on the patch. For the first part, we have three options.
The first is to move the particles in the physical domain. Then we need to know the inverse of the change of coordinates in order to transform back on the patch for interpolation. That could become complicated and costly because we don't know in general an analytical form for the inverse of the change of coordinates.
The second is to move directly all the particles in the patch. In this case, we need to solve the equations of motion in the patch coordinates, where they become very involved and more costly to solve.
The third option proposed, and that we have chosen, is to move the positions of the particles in the patch and the velocities of the particles in the physical domain. With this strategy, we do not need to know the inverse of the change of coordinates and the equations of motion are simpler.
