High order Asymptotic Preserving penalized numerical schemes for the Euler-Poisson system in the quasi-neutral limit

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Crouseilles, Nicolas, Dimarco, Giacomo, Samantaray, Saurav
Format: Preprint
Veröffentlicht: 2024
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866912018159632384
author Crouseilles, Nicolas
Dimarco, Giacomo
Samantaray, Saurav
author_facet Crouseilles, Nicolas
Dimarco, Giacomo
Samantaray, Saurav
contents In this work, we focus on the development of high-order Implicit-Explicit (IMEX) finite volume numerical methods for plasmas in quasineutral regimes. At large temporal and spatial scales, plasmas tend to be quasineutral, meaning that the local net charge density is nearly zero. However, at small time and spatial scales, measured by the the Debye length, quasineutrality breaks down. In such regimes, standard numerical methods face severe stability constraints, rendering them practically unusable. To address this issue, we introduce and analyze a class of penalized IMEX Runge-Kutta methods for the Euler-Poisson (EP) system, specifically designed to handle the quasineutral limit. These schemes are uniformly stable with respect to the Debye length and degenerate into high-order methods as the quasineutral limit is approached. Several numerical tests confirm that the proposed methods exhibit the desired properties.
format Preprint
id arxiv_https___arxiv_org_abs_2409_04807
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle High order Asymptotic Preserving penalized numerical schemes for the Euler-Poisson system in the quasi-neutral limit
Crouseilles, Nicolas
Dimarco, Giacomo
Samantaray, Saurav
Numerical Analysis
In this work, we focus on the development of high-order Implicit-Explicit (IMEX) finite volume numerical methods for plasmas in quasineutral regimes. At large temporal and spatial scales, plasmas tend to be quasineutral, meaning that the local net charge density is nearly zero. However, at small time and spatial scales, measured by the the Debye length, quasineutrality breaks down. In such regimes, standard numerical methods face severe stability constraints, rendering them practically unusable. To address this issue, we introduce and analyze a class of penalized IMEX Runge-Kutta methods for the Euler-Poisson (EP) system, specifically designed to handle the quasineutral limit. These schemes are uniformly stable with respect to the Debye length and degenerate into high-order methods as the quasineutral limit is approached. Several numerical tests confirm that the proposed methods exhibit the desired properties.
title High order Asymptotic Preserving penalized numerical schemes for the Euler-Poisson system in the quasi-neutral limit
topic Numerical Analysis
url https://arxiv.org/abs/2409.04807