BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Eric Sonnendrücker

Eric Sonnendrücker appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysismath-phmath.APmath.MPphysics.comp-phphysics.plasm-phComputational Engineering, Finance, and Sciencephysics.class-ph

Catalog footprint

What is connected

14works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

14 published item(s)

preprint2022arXiv

A Hamiltonian model for the macroscopic Maxwell equations using exterior calculus

A Hamiltonian field theory for the macroscopic Maxwell equations with fully general polarization and magnetization is stated in the language of differential forms. The precise procedure for translating the vector calculus formulation into differential forms is discussed in detail. We choose to distinguish between straight and twisted differential forms so that all integrals be taken over densities (i.e. twisted top forms). This ensures that the duality pairings, which are stated as integrals over densities, are orientation independent. The relationship between functional differentiation with respect to vector fields and with respect to differential forms is established using the chain rule. The theory is developed such that the Poisson bracket is metric and orientation independent with all metric dependence contained in the Hamiltonian. As is typically seen in the exterior calculus formulation of Maxwell's equations, the Hodge star operator plays a key role in modeling the constitutive relations. As a demonstration of the kind of constitutive models this theory accommodates, the paper concludes with several examples.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A mimetic discretization of the macroscopic Maxwell equations in Hamiltonian form

A mimetic spectral element discretization, utilizing a novel Galerkin projection Hodge star operator, of the macroscopic Maxwell equations in Hamiltonian form is presented. The idea of splitting purely topological and metric dependent quantities is natural in the Hamiltonian modeling framework as the Poisson bracket is metric free with the Hamiltonian containing all metric information. This idea may be incorporated into the mimetic spectral element method by directly discretizing the Poincaré duality structure. This "split exterior calculus mimetic spectral element method" yields spatially discretized Maxwell's equations which are Hamiltonian and exactly and strongly conserve Gauss's laws. Moreover, the new discrete Hodge star operator is itself of interest as a partition of the purely topological and metric dependent portions of the Hodge star operator. As a simple test case, the numerical results of applying this method to a one-dimensional version of Maxwell's equations are given.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On Geometric Fourier Particle In Cell Methods

In this article we describe a unifying framework for variational electromagnetic particle schemes of spectral type, and we propose a novel spectral Particle-In-Cell (PIC) scheme that preserves a discrete Hamiltonian structure. Our work is based on a new abstract variational derivation of particle schemes which builds on a de Rham complex where Low's Lagrangian is discretized using a particle approximation of the distribution function. In this framework, which extends the recent Finite Element based Geometric Electromagnetic PIC (GEMPIC) method to a variety of field solvers, the discretization of the electromagnetic potentials and fields is represented by a de Rham sequence of compatible spaces, and the particle-field coupling procedure is described by approximation operators that commute with the differential operators in the sequence. In particular, for spectral Maxwell solvers the choice of truncated $L^2$ projections using continuous Fourier transform coefficients for the commuting approximation operators yields the gridless Particle-in-Fourier method, whereas spectral Particle-in-Cell methods are obtained by using discrete Fourier transform coefficients computed from a grid. By introducing a new sequence of spectral pseudo-differential approximation operators, we then obtain a novel variational spectral PIC method with discrete Hamiltonian structure that we call Fourier-GEMPIC. Fully discrete schemes are then derived using a Hamiltonian splitting procedure, leading to explicit time steps that preserve the Gauss laws and the discrete Poisson bracket associated with the Hamiltonian structure. These explicit steps share many similarities with standard spectral PIC methods. As arbitrary filters are allowed in our framework, we also discuss aliasing errors and study a natural back-filtering procedure to mitigate the damping caused by anti-aliasing smoothing particle shapes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Variational Framework for Structure-Preserving Electromagnetic Particle-In-Cell Methods

In this article we apply a discrete action principle for the Vlasov--Maxwell equations in a structure-preserving particle-field discretization framework. In this framework the finite-dimensional electromagnetic potentials and fields are represented in a discrete de Rham sequence involving general finite element spaces, and the particle-field coupling is represented by a set of projection operators that commute with the differential operators. With a minimal number of assumptions which allow for a variety of finite elements and shape functions for the particles, we show that the resulting variational scheme has a general discrete Poisson structure and thus leads to a semi-discrete Hamiltonian system. By introducing discrete interior products we derive a second type of space discretization which is momentum preserving, based on the same finite elements and shape functions. We illustrate our method by applying it to spline finite elements, and to a new spectral discretization where the particle-field coupling relies on discrete Fourier transforms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Geometric Particle-in-Cell Simulations of the Vlasov--Maxwell System in Curvilinear Coordinates

Numerical schemes that preserve the structure of the kinetic equations can provide stable simulation results over a long time. An electromagnetic particle-in-cell solver for the Vlasov-Maxwell equations that preserves at the discrete level the non-canonical Hamiltonian structure of the Vlasov-Maxwell equations has been presented in [Kraus et al. 2017]. Whereas the original formulation has been obtained for Cartesian coordinates, we extend the formulation to curvilinear coordinates in this paper. For the discretisation in time, we discuss several (semi-)implicit methods either based on a Hamiltonian splitting or a discrete gradient method combined with an antisymmetric splitting of the Poisson matrix and discuss their conservation properties and computational efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Energy conservation and numerical stability for the reduced MHD models of the non-linear JOREK code

In this paper we present a rigorous derivation of the reduced MHD models with and without parallel velocity that are implemented in the non-linear MHD code JOREK. The model we obtain contains some terms that have been neglected in the implementation but might be relevant in the non-linear phase. These are necessary to guarantee exact conservation with respect to the full MHD energy. For the second part of this work, we have replaced the linearized time stepping of JOREK by a non-linear solver based on the Inexact Newton method including adaptive time stepping. We demonstrate that this approach is more robust especially with respect to numerical errors in the saturation phase of an instability and allows to use larger time steps in the non-linear phase.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulations of Kinetic Electrostatic Electron Nonlinear (KEEN) Waves with Variable Velocity Resolution Grids and High-Order Time-Splitting

KEEN waves are nonlinear, non-stationary, self-organized asymptotic states in Vlasov plasmas outside the scope or purview of linear theory constructs such as electron plasma waves or ion acoustic waves. Nonlinear stationary mode theories such as those leading to BGK modes also do not apply. The range in velocity that is strongly perturbed by KEEN waves depends on the amplitude and duration of the ponderomotive force used to drive them. Smaller amplitude drives create highly localized structures attempting to coalesce into KEEN waves. These cases have much more chaotic and intricate time histories than strongly driven ones. The narrow range in which one must maintain adequate velocity resolution in the weakly driven cases challenges xed grid numerical schemes. What is missing there is the capability of resolving locally in velocity while maintaining a coarse grid outside the highly perturbed region of phase space. We here report on a new Semi-Lagrangian Vlasov-Poisson solver based on conservative non-uniform cubic splines in velocity that tackles this problem head on. An additional feature of our approach is the use of a new high-order time-splitting scheme which allows much longer simulations per computational e ort. This is needed for low amplitude runs which take a long time to set up KEEN waves, if they are able to do so at all. The new code's performance is compared to uniform grid simulations and the advantages quanti ed. The birth pains associated with KEEN waves which are weakly driven is captured in these simulations. These techniques allow the e cient simulation of KEEN waves in multiple dimensions which will be tackled next as well as generalizations to Vlasov-Maxwell codes which are essential to understanding the impact of KEEN waves in practice.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An exponential integrator for a highly oscillatory Vlasov equation

In the framework of a Particle-In-Cell scheme for some 1D Vlasov-Poisson system depending on a small parameter, we propose a time-stepping method which is numerically uniformly accurate when the parameter goes to zero. Based on an exponential time differencing approach, the scheme is able to use large time steps with respect to the typical size of the fast oscillations of the solution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Some numerical aspects of the conservative PSM scheme in a 4D drift-kinetic code

The purpose of this work is simulation of magnetised plasmas in the ITER project framework. In this context, kinetic Vlasov-Poisson like models are used to simulate core turbulence in the tokamak in a toroidal geometry. This leads to heavy simulations because a 6D dimensional problem has to be solved, even if reduced to a 5D in so called gyrokinetic models. Accurate schemes, parallel algorithms need to be designed to bear these simulations. This paper describes the numerical studies to improve robustness of the conservative PSM scheme in the context of its development in the GYSELA code. In this paper, we only consider the 4D drift-kinetic model which is the backbone of the 5D gyrokinetic models and relevant to build a robust and accurate numerical method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accuracy of unperturbed motion of particles in a gyrokinetic semi-Lagrangian code

Inaccurate description of the equilibrium can yield to spurious effects in gyrokinetic turbulence simulations. Also, the Vlasov solver and time integration schemes impact the conservation of physical quantities, especially in long-term simulations. Equilibrium and Vlasov solver have to be tuned in order to preserve constant states (equilibrium) and to provide good conservation property along time (mass to begin with). Several illustrative simple test cases are given to show typical spurious effects that one can observes for poor settings. We explain why Forward Semi-Lagrangian scheme bring us some benefits. Some toroidal and cylindrical GYSELA runs are shown that use FSL.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Analysis of a new class of Forward Semi-Lagrangian schemes for the 1D Vlasov-Poisson Equations

The Vlasov equation is a kinetic model describing the evolution of charged particles, and is coupled with Poisson's equation, which rules the evolution of the self-consistent electric field. In this paper, we introduce a new class of forward Semi-Lagrangian schemes for the Vlasov-Poisson system based on a Cauchy Kovalevsky (CK) procedure for the numerical solution of the characteristic curves. Exact conservation properties of the first moments of the distribution function for the schemes are derived and a convergence study is performed that applies as well for the CK scheme as for a more classical Verlet scheme. The convergence in L1 norm of the schemes is proved and error estimates are obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Homogenization of the Vlasov Equation and of the Vlasov - Poisson System with a Strong External Magnetic Field

Motivated by the difficulty arising in the numerical simulation of the movement of charged particles in presence of a large external magnetic field, which adds an additional time scale and thus imposes to use a much smaller time step, we perform in this paper a homogenization of the Vlasov equation and of the Vlasov-Poisson system which yield approximate equations describing the mean behavior of the particles. The convergence proof is based on the two scale convergence tools introduced by N'Guetseng and Allaire. We also consider the case where, in addition to the magnetic field, a large external electric field orthogonal to the magnetic field and of the same magnitude is applied.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Test of some numerical limiters for the conservative PSM scheme for 4D Drift-Kinetic simulations

The purpose of this work is simulation of magnetised plasmas in the ITER project framework. In this context, Vlasov-Poisson like models are used to simulate core turbulence in the tokamak in a toroidal geometry. This leads to heavy simulation because a 6D dimensional problem has to be solved, 3D in space and 3D in velocity. The model is reduced to a 5D gyrokinetic model, taking advantage of the particular motion of particles due to the presence of a strong magnetic field. However, accurate schemes, parallel algorithms need to be designed to bear these simulations. This paper describes a Hermite formulation of the conservative PSM scheme which is very generic and allows to implement different semi-Lagrangian schemes. We also test and propose numerical limiters which should improve the robustness of the simulations by diminishing spurious oscillations. We only consider here the 4D drift-kinetic model which is the backbone of the 5D gyrokinetic models and relevant to build a robust and accurate numerical method.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

A Forward semi-Lagrangian Method for the Numerical Solution of the Vlasov Equation

This work deals with the numerical solution of the Vlasov equation. This equation gives a kinetic description of the evolution of a plasma, and is coupled with Poisson's equation for the computation of the self-consistent electric field. The coupled model is non linear. A new semi-Lagrangian method, based on forward integration of the characteristics, is developed. The distribution function is updated on an eulerian grid, and the pseudo-particles located on the mesh's nodes follow the characteristics of the equation forward for one time step, and are deposited on the 16 nearest nodes. This is an explicit way of solving the Vlasov equation on a grid of the phase space, which makes it easier to develop high order time schemes than the backward method.

0 citations0 reviewsNo code linkedOpen access