https://doi.org/10.1051/epjconf/202634601002
A structure-preserving spline finite element solver for the cold-plasma model
1 Max Planck Institute for Plasma Physics, Garching, Germany.
2 Department of Mathematics, Technical University of Munich
* Corresponding author: This email address is being protected from spambots. You need JavaScript enabled to view it.
** Corresponding author: This email address is being protected from spambots. You need JavaScript enabled to view it.
*** Corresponding author: This email address is being protected from spambots. You need JavaScript enabled to view it.
**** Corresponding author: This email address is being protected from spambots. You need JavaScript enabled to view it.
Published online: 7 January 2026
We present a spline finite element solver, which preserves the Hamiltonian structure of the coldplasma model as well as several physical invariants, such as the energy, the total charge and the zero divergence of the magnetic field. The scheme is naturally adapted to Cartesian and curvilinear geometries. A key feature of the scheme is that in the presence of a time-harmonic source, it is consistent with a high-order approximation of the associated time-harmonic solution: this makes the solver intrinsically stable and long simulation runs cannot develop unphysical effects. Our implementation relies on PSYDAC, an isogeometric B-splines finite-element library that can be used to build efficient solvers based on modern numerical methods. In this paper we include an overview of the library and present an example of implementation.
© The Authors, published by EDP Sciences, 2026
This is an Open Access article distributed under the terms of the Creative Commons Attribution License 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
