Abstract:In this paper, the design and analysis of high order accurate IMEX finite volume schemes for the compressible Euler–Poisson (EP) equations in the quasineutral limit is presented. As the quasineutral limit is singular for the governing equations, the time discretisation is tantamount to achieving an accurate numerical method. To this end, the EP system is viewed as a differential algebraic equation system (DAEs) via static condensation. As a consequence of this vantage point, high order linearly semi-implicit (SI) time discretisation are realised by employing a novel combination of the direct approach used for implicit discretisation of DAEs and, two different classes of IMEX-RK schemes: the additive and the multiplicative. For both the time discretisation strategies, in order to account for rapid plasma oscillations in quasineutral regimes, the nonlinear Euler fluxes are split into two different combinations of stiff and non-stiff components. The high order scheme resulting from the additive approach is designated as a classical scheme while the one generated by the multiplicative approach possesses the asymptotic preserving (AP) property. Time discretisations for the classical and the AP schemes are performed by standard IMEX-RK and SI-IMEX-RK methods, respectively so that the stiff terms are treated implicitly and the non-stiff ones explicitly. In order to discretise in space a Rusanov-type central flux is used for the non-stiff part, and simple central differencing for the stiff part. Results of numerical experiments are presented, which confirm that the high order schemes based on the SI-IMEX-RK time discretisation achieve uniform second order convergence with respect to the Debye length and are AP in the quasineutral limit.

A Decoupled, Linearly Implicit and High-Order Structure-Preserving Scheme for Euler-Poincaré Equations

A New High-Order Compact Scheme Of Unstructured Finite Volume Method

High-Order Linearly Implicit Structure-Preserving Exponential Integrators for the Nonlinear Schrödinger Equation

An Efficient High-Order Explicit Scheme for Solving Hamiltonian Nonlinear Wave Equations

Structure-preserving numerical schemes for Lindblad equations

High Resolution Optimized High-Order Schemes for Discretization of Non-Linear Straight and Mixed Second Derivative Terms

The Onsager principle and structure preserving numerical schemes

High-order conservative and accurately dissipative numerical integrators via auxiliary variables

A Linearly-Fitted Conservative (dissipative) Scheme for Efficiently Solving Conservative (dissipative) Nonlinear Wave PDEs

Decoupled structure-preserving discretization of incompressible MHD equations with general boundary conditions

A structure preserving difference scheme with fast algorithms for high dimensional nonlinear space-fractional Schrödinger equations

High Order Asymptotic Preserving and Classical Semi-implicit RK Schemes for the Euler–Poisson System in the Quasineutral Limit

A Linearly Implicit Structure-Preserving Scheme for the Camassa–Holm Equation Based on Multiple Scalar Auxiliary Variables Approach

Order-lifted data inversion/retrieval method of neighbor cells to implement general high-order schemes in unstructured-mesh-based finite-volume solution framework

A structure-preserving discontinuous Galerkin scheme for the Cahn-Hilliard equation including time adaptivity

Two linear energy-preserving compact finite difference schemes for coupled nonlinear wave equations

High-Order Schemes of Exponential Time Differencing for Stiff Systems with Nondiagonal Linear Part

A positivity-preserving implicit-explicit scheme with high order polynomial basis for compressible Navier-Stokes equations

Global error estimates of high-order fully decoupled schemes for the Cahn-Hilliard-Navier-Stokes model of Two-Phase Incompressible Flows

Constraint-preserved numerical schemes with decoupling structure for the Ericksen–Leslie model with variable density