Conservative discontinuous Galerkin interpolation: sheared boundary conditions

Manaure Francisquez,Noah R. Mandell,Ammar Hakim,Gregory W. Hammett
DOI: https://doi.org/10.1016/j.cpc.2024.109109
IF: 4.717
2024-01-31
Computer Physics Communications
Abstract:Local studies of accretion disks and laboratory magnetized plasmas employ analytical coordinate mappings that introduce sheared boundary conditions (BCs). We present a discontinuous Galerkin (DG) algorithm to apply such BCs based on projections and quadrature-free integration. The procedure is high-order accurate, preserves moments exactly and works in multiple dimensions. Tests of increasing complexity are provided, beginning with translations of one and two dimensional fields, followed by 3D and 5D simulations with sheared (twist-shift) BCs. Results show that the algorithm is (p+1) -order accurate in the DG representation and (p+2) -order accurate in the cell averages, with p being the order of the polynomial basis. Quantification of the algorithm's hyperdiffusion and discussion of aliasing errors are given. This technique enables conservative local simulations of plasma turbulence with DG, not possible until now.
physics, mathematical,computer science, interdisciplinary applications
What problem does this paper attempt to address?