{
  "abstract": "In this paper, we present a high-order, stable summation-by-parts (SBP) finite difference method for solving the Vlasov-Maxwell system in a 2D2V phase space. Central SBP operators for the advection terms are not stable when the solution becomes non-smooth and fine-scale filamentary structures develop, as is typical in high-dimensional Vlasov-Maxwell simulations. To address this issue, the method is stabilized using high-order upwind SBP operators. High-order explicit Runge-Kutta methods are employed for time integration. We prove that the fully discrete scheme exactly conserves mass and preserves momentum up to truncation error. Furthermore, we present a matrix-free implementation of the method on modern GPU architectures. A range of challenging benchmark problems is solved to demonstrate the accuracy, robustness, and performance of the proposed scheme.",
  "abstract_shown": true,
  "affiliations": [],
  "article_type": null,
  "arxiv": "2609.06452",
  "author_total": 3,
  "authors": [
    {
      "name": "Robin Dymér"
    },
    {
      "name": "Ken Mattsson"
    },
    {
      "name": "Murtazo Nazarov"
    }
  ],
  "bibtex": "https://research.draconitic.com/paper/fa0h06wd3p9f3.bib",
  "cited_by": 0,
  "doi": "10.48550/arxiv.2609.06452",
  "fulltext_url": null,
  "id": "fa0h06wd3p9f3",
  "issns": [],
  "landing_url": "https://arxiv.org/abs/2609.06452",
  "language": null,
  "licence": "CC-BY-4.0",
  "licence_label": "CC BY 4.0",
  "licence_source": "article",
  "licence_url": "http://creativecommons.org/licenses/by/4.0/",
  "metadata_licence": "CC0-1.0",
  "native_subjects": [
    [
      "arxiv",
      "cs.NA"
    ],
    [
      "arxiv",
      "math.NA"
    ],
    [
      "arxiv",
      "physics.comp-ph"
    ],
    [
      "arxiv",
      "physics.plasm-ph"
    ],
    [
      "openalex",
      "2206|Computational Fluid Dynamics and Aerodynamics"
    ]
  ],
  "openalex": "W7211999016",
  "pdf_url": "https://arxiv.org/pdf/2609.06452",
  "pmcid": null,
  "pmid": null,
  "published": "2026-09-06",
  "publisher": null,
  "ris": "https://research.draconitic.com/paper/fa0h06wd3p9f3.ris",
  "sources": [
    {
      "first_seen": "2026-10-30T12:19:48Z",
      "native_id": "2609.06452",
      "source": "arxiv",
      "url": "https://arxiv.org/abs/2609.06452"
    }
  ],
  "subjects": [
    {
      "field": "mathematics",
      "origin": "arxiv:math.NA",
      "subfield": "numerical-analysis"
    },
    {
      "field": "physics",
      "origin": "arxiv:physics.comp-ph",
      "subfield": "computational"
    },
    {
      "field": "physics",
      "origin": "arxiv:physics.plasm-ph",
      "subfield": "plasma"
    }
  ],
  "title": "High-Order Structure-Preserving SBP Finite Difference Methods for the Vlasov-Maxwell System on Matrix-Free GPUs",
  "title_text": "High-Order Structure-Preserving SBP Finite Difference Methods for the Vlasov-Maxwell System on Matrix-Free GPUs",
  "type": "preprint",
  "updated_at": "2026-10-30T12:19:58Z",
  "url": "https://research.draconitic.com/paper/fa0h06wd3p9f3",
  "venue": "arXiv",
  "versions": []
}