Source author record

David Coster

David Coster appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.plasm-ph physics.comp-ph

Catalog footprint

What is connected

4works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

Gyrokinetic full-f particle-in-cell simulations on open field lines with PICLS

While in recent years, gyrokinetic simulations have become the workhorse for theoretical turbulence and transport studies in the plasma core, their application to the edge and scrape-off layer (SOL) region presents significant challenges. In particular, steep density and temperature gradients as well as large fluctuation amplitudes call for a "full-f" treatment. To specifically study problems in the SOL region, the gyrokinetic particle-in-cell (PIC) code PICLS has been developed. The code is based on an electrostatic full-f model with linearised field equations and uses kinetic electrons. Here, the well-studied parallel transport problem during an edge-localized mode (ELM) in the SOL shall be investigated for one spatial dimension. The results are compared to previous gyrokinetic continuum and fully kinetic PIC simulations and show good agreement.

preprint2019arXiv

Collisional gyrokinetic full-f particle-in-cell simulations on open field lines with PICLS

Applying gyrokinetic simulations for theoretical turbulence and transport studies to the plasma edge and scrape-off layer (SOL) presents significant challenges. To in particular account for steep density and temperature gradients in the SOL, the "full-f" code PICLS was developed. PICLS is a gyrokinetic particle-in-cell (PIC) code and is based on an electrostatic model with a linearized field equation and uses kinetic electrons. In previously published results we were applying PICLS to the well-studied 1D parallel transport problem during an edge-localized mode (ELM) in the SOL without collisions. As an extension to this collision-less case and in preparation for 3D simulations, in this work a collisional model will be introduced. The implemented Lenard-Bernstein collision operator and its Langevin discretization will be shown. Conservation properties of the collision operator as well as a comparison of the collisional and non-collisional case will be discussed.

preprint2019arXiv

Treatment of Advanced Divertor Configurations in the Flux-Coordinate Independent turbulence code GRILLIX

Advanced divertor configurations modify the magnetic geometry of the diverter to achieve a combination of strong magnetic flux expansion, increased connection length and higher divertor volume - to improve detachment stability, neutral/impurity confinement and heat-channel broadening. In this paper, we discuss the modification of the Flux-Coordinate Independent (FCI) turbulence code GRILLIX to treat generalised magnetic geometry, to allow for the investigation of the effect of magnetic geometry on turbulent structures in the edge and SOL. The development of grids and parallel operators from numerically-defined magnetic equilibria is discussed, as is the application of boundary conditions via penalisation, with the finite-width method generalised to treat complex non-conformal boundaries. Initial testing of hyperbolic (advection) and parabolic (diffusion) test cases is presented for the Snowflake scenario.

preprint2015arXiv

The field line map approach for simulations of magnetically confined plasmas

In the presented field line map approach the simulation domain of a tokamak is covered with a cylindrical grid, which is Cartesian within poloidal planes. Standard finite-difference methods can be used for the discretisation of perpendicular (w.r.t.~magnetic field lines) operators. The characteristic flute mode property $\left(k_{\parallel}\ll k_{\perp}\right)$ of structures is exploited computationally by a grid sparsification in the toroidal direction. A field line following discretisation of parallel operators is then required, which is achieved via a finite difference along magnetic field lines. This includes field line tracing and interpolation or integration. The main emphasis of this paper is on the discretisation of the parallel diffusion operator. Based on the support operator method a scheme is constructed which exhibits only very low numerical perpendicular diffusion. The schemes are implemented in the new code GRILLIX, and extensive benchmarks are presented which show the validity of the approach in general and GRILLIX in particular. The main advantage of the approach is that it does not rely on field/flux-aligned, which become singular on the separatrix/X-point. Most tokamaks are based on the divertor concept, and the numerical treatment of the separatrix is therefore of importance for simulations of the edge and scrape-off layer.