Researcher profile

Stephan Brunner

Stephan Brunner contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

physics.plasm-ph physics.comp-ph physics.flu-dyn

Trust snapshot

Quick read

Trust 17 - Baseline

4works

0followers

3topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

preprint2020arXiv

Eliminating turbulent self-interaction through the parallel boundary condition in local gyrokinetic simulations

In this work, we highlight an issue that may reduce the accuracy of many local nonlinear gyrokinetic simulations - turbulent self-interaction through the parallel boundary condition. Given a sufficiently long parallel correlation length, individual turbulent eddies can span the full domain and "bite their own tails," thereby altering their statistical properties. Such self-interaction is only modeled accurately when the simulation domain corresponds to a full flux surface, otherwise it is artificially strong. For Cyclone Base Case parameters and typical domain sizes, we find that this mechanism modifies the heat flux by roughly 40% and it can be even more important. The effect is largest when using kinetic electrons, low magnetic shear, and strong turbulence drive (i.e. steep background gradients). It is found that parallel self-interaction can be eliminated by increasing the parallel length and/or the binormal width of the simulation domain until convergence is achieved.

preprint2020arXiv

Gyrokinetic Simulations on Many- and Multi-core Architectures with the Global Electromagnetic Particle-In-Cell Code ORB5

Gyrokinetic codes in plasma physics need outstanding computational resources to solve increasingly complex problems, requiring the effective exploitation of cutting-edge HPC architectures. This paper focuses on the enabling of ORB5, a state-of-the-art, first-principles-based gyrokinetic code, on modern parallel hybrid multi-core, multi-GPU systems. ORB5 is a Lagrangian, Particle-In-Cell (PIC), finite element, global, electromagnetic code, originally implementing distributed parallelism through MPI-based on domain decomposition and domain cloning. In order to support multi/many cores devices, the code has been completely refactored. Data structures have been re-designed to ensure efficient memory access, enhancing data locality. Multi-threading has been introduced through OpenMP on the CPU and adopting OpenACC to support GPU acceleration. MPI can still be used in combination with the two approaches. The performance results obtained using the full production ORB5 code on the Summit system at ORNL, on Piz Daint at CSCS and on the Marconi system at CINECA are presented, showing the effectiveness and performance portability of the adopted solutions: the same source code version was used to produce all results on all architectures.

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.

preprint2012arXiv

Locality properties of the free energy fluxes in gyrokinetic turbulence

The nature of the nonlinear interactions in gyrokinetic (GK) turbulence, driven by an ion-temperature gradient instability, is investigated using numerical simulations of single ion species plasma in three-dimensional flux tube geometry. To account for the level of separation existing between scales involved in an energetic interaction, the degree of locality of the free energy scale flux is analyzed employing Kraichnan's infrared (IR) and ultraviolet locality functions. Due to the nontrivial dissipative nature of GK turbulence, an asymptotic level for the locality exponents, indicative of a universal dynamical regime for GK's, is not recovered and an accentuated non-local behavior of the IR interactions is found instead, in spite of the local energy cascade observed.