Standard versus Asymptotic Preserving Time Discretizations for the Poisson-Nernst-Planck System in the Quasi-Neutral Limit

Clarissa Astuto

arxiv: 2511.07964 · v2 · submitted 2025-11-11 · 🧮 math.NA · cs.NA

Standard versus Asymptotic Preserving Time Discretizations for the Poisson-Nernst-Planck System in the Quasi-Neutral Limit

Clarissa Astuto This is my paper

Pith reviewed 2026-05-18 00:07 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords Poisson-Nernst-Planck systemasymptotic preserving schemesquasi-neutral limitDebye lengthIMEX time discretizationnumerical stabilityion diffusion

0 comments

The pith

IMEX schemes for the Poisson-Nernst-Planck system remain asymptotically stable for every Debye length and need no special initial conditions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies time discretizations for the Poisson-Nernst-Planck system that models correlated diffusion of two ion species. It validates an asymptotic-preserving scheme that correctly handles the quasi-neutral limit reached when the Debye length becomes very small, at which point the system collapses to a single diffusion equation with an effective coefficient. Standard schemes encounter severe stability limits and initial-condition difficulties for small but nonzero Debye lengths. IMEX schemes, by contrast, stay stable for all Debye lengths and require no assumptions on the starting data. The work compares the different discretizations to exhibit these asymptotic properties.

Core claim

The IMEX time discretization schemes are proved to be asymptotically stable for all Debye lengths and do not require any assumption on the initial conditions, while standard schemes suffer stability restrictions when the Debye length is small but not negligible.

What carries the argument

The Asymptotic-Preserving (AP) IMEX time discretization for the Poisson-Nernst-Planck system with respect to the Debye length, which preserves the quasi-neutral limit reduction to a diffusion equation.

If this is right

For vanishing Debye lengths the scheme reduces exactly to the single diffusion equation with the effective coefficient.
Stability holds uniformly for every positive Debye length.
The method works without any special preparation of initial data.
Direct comparison demonstrates that only the AP schemes capture the correct asymptotic regime.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same stability property could simplify long-time simulations of ion transport whenever the Debye length is orders of magnitude smaller than other macroscopic scales.
Comparable AP constructions may prove useful for other singularly perturbed electrodiffusion models that reduce to diffusion equations in a limit.
The absence of initial-condition restrictions removes a practical barrier that often appears in quasi-neutral plasma or electrolyte calculations.

Load-bearing premise

The numerical scheme proposed in the cited prior work is asymptotically preserving with respect to the Debye length, and the quasi-neutral limit reduces the system to a single diffusion equation with an effective diffusion coefficient.

What would settle it

A numerical experiment with successively smaller Debye lengths that shows instability or wrong limiting behavior in a standard scheme but stable convergence to the correct diffusion equation in the IMEX scheme.

read the original abstract

In this paper, we investigate the correlated diffusion of two ion species governed by a Poisson-Nernst-Planck (PNP) system. Here we further validate the numerical scheme recently proposed in \cite{astuto2025asymptotic}, where a time discretization method was shown to be Asymptotic-Preserving (AP) with respect to the Debye length. For vanishingly Debye lengths, the so called Quasi-Neutral limit can be adopted, reducing the system to a single diffusion equation with an effective diffusion coefficient \cite{CiCP-31-707}. Choosing small, but not negligible, Debye lengths, standard numerical methods suffer from severe stability restrictions and difficulties in handling initial conditions. IMEX schemes, on the other hand, are proved to be asymptotically stable for all Debye lengths, and do not require any assumption on the initial conditions. In this work, we compare different time discretizations to show their asymptotic behaviors.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper validates the authors' prior IMEX scheme for the PNP system via direct comparisons in the quasi-neutral limit rather than introducing new theory.

read the letter

This paper is mainly a numerical validation of an IMEX time discretization the authors proposed earlier for the Poisson-Nernst-Planck system. They compare it to standard schemes as the Debye length shrinks toward the quasi-neutral regime, where the system reduces to a single diffusion equation with an effective coefficient. The central point is that the IMEX approach stays stable for any Debye length and works with arbitrary initial data, while standard methods hit stability walls and need compatible starting conditions. The comparisons of asymptotic behavior are the concrete contribution here. They give a practical sense of how the schemes perform when the Debye length is small but not zero, which is the regime that matters for many simulations. The work does a reasonable job using the limit equation as a benchmark to check consistency. The soft spot is that the asymptotic-preserving property and the stability proof are taken from the earlier citation. This manuscript adds the side-by-side tests but does not appear to supply independent derivations or broad error estimates on its own. If the numerical examples are limited in scope or resolution, that would be the main place a referee might ask for more. Readers who run ion-transport codes in electrochemistry or biophysics would get the most out of it, especially if they already know the prior paper and want to see how the scheme holds up in practice. It is not a foundational new framework, but the targeted comparisons address a real computational headache. I would send it to peer review. The topic is relevant and the evidence is at least plausible from the abstract and structure, even if the full strength depends on the details in the numerics section.

Referee Report

2 major / 2 minor

Summary. The manuscript compares standard and IMEX time discretizations for the Poisson-Nernst-Planck system governing correlated diffusion of two ion species. It validates the asymptotic-preserving (AP) property of an IMEX scheme from prior work with respect to the Debye length, showing that the scheme remains stable for arbitrarily small Debye lengths in the quasi-neutral limit without requiring compatible initial conditions, while standard schemes encounter severe stability restrictions. Numerical comparisons illustrate that the IMEX approach reproduces the limiting single diffusion equation with effective coefficient as the Debye length vanishes.

Significance. If the stability and asymptotic behavior hold under the reported conditions, the work supplies concrete numerical evidence supporting the use of AP-IMEX schemes for stiff PNP problems in the quasi-neutral regime. This is relevant for applications in electrokinetics and semiconductor modeling where small Debye lengths induce numerical stiffness. The manuscript appropriately credits the prior AP proof and focuses on direct scheme comparisons, adding value through reproducible validation of uniform stability.

major comments (2)

[§3] §3 (or the section presenting the IMEX scheme): the statement that IMEX schemes are 'proved to be asymptotically stable for all Debye lengths, and do not require any assumption on the initial conditions' is load-bearing for the central claim but is asserted by reference to the prior work without a self-contained outline of the key stability argument or the precise role of the initial-data compatibility condition; a short recap would allow readers to verify the no-assumption claim directly.
[Numerical results] Numerical results section (e.g., around the tables or figures comparing schemes for λ = 10^{-4} to 10^{-6}): the reported time-step restrictions and qualitative asymptotic behaviors are shown, but quantitative error measures (such as L2 deviation from the quasi-neutral diffusion solution or observed convergence rates as λ → 0) are not provided; without these, it is difficult to confirm that the IMEX scheme accurately captures the effective diffusion coefficient from the cited limit reduction.

minor comments (2)

Figure captions should explicitly label which curves correspond to the standard scheme, the IMEX scheme, and the quasi-neutral reference solution for each Debye length tested.
[Introduction] The introduction could briefly recall the explicit form of the effective diffusion coefficient in the quasi-neutral limit (from CiCP-31-707) rather than only citing it, to make the comparison metric clearer.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment and constructive comments on our manuscript. We address each major comment below and will incorporate the suggested improvements in the revised version.

read point-by-point responses

Referee: [§3] §3 (or the section presenting the IMEX scheme): the statement that IMEX schemes are 'proved to be asymptotically stable for all Debye lengths, and do not require any assumption on the initial conditions' is load-bearing for the central claim but is asserted by reference to the prior work without a self-contained outline of the key stability argument or the precise role of the initial-data compatibility condition; a short recap would allow readers to verify the no-assumption claim directly.

Authors: We agree that a brief self-contained recap would strengthen the presentation. In the revised manuscript we will add a short paragraph in the section introducing the IMEX scheme that outlines the main steps of the asymptotic stability argument from the cited prior work and clarifies that the IMEX discretization does not impose compatibility requirements on the initial data, in contrast to standard schemes. revision: yes
Referee: [Numerical results] Numerical results section (e.g., around the tables or figures comparing schemes for λ = 10^{-4} to 10^{-6}): the reported time-step restrictions and qualitative asymptotic behaviors are shown, but quantitative error measures (such as L2 deviation from the quasi-neutral diffusion solution or observed convergence rates as λ → 0) are not provided; without these, it is difficult to confirm that the IMEX scheme accurately captures the effective diffusion coefficient from the cited limit reduction.

Authors: The referee is correct that quantitative error measures are not currently reported. Although the present numerical study emphasizes uniform stability and qualitative reproduction of the limiting behavior, we will augment the numerical results section with L2-error tables measuring deviation from the quasi-neutral diffusion solution together with observed convergence rates as λ → 0, thereby providing direct quantitative confirmation that the effective diffusion coefficient is recovered. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper cites prior work by the same author to establish the asymptotic-preserving property of the IMEX scheme but performs independent numerical comparisons of asymptotic behaviors across different time discretizations for small but nonzero Debye lengths. These comparisons directly demonstrate stability and lack of initial-condition restrictions without reducing any claim to the self-citation by construction. No self-definitional steps, fitted inputs renamed as predictions, or load-bearing uniqueness theorems appear in the argument structure. The derivation chain for the reported behaviors remains self-contained via the present manuscript's experiments.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the asymptotic-preserving property established in the referenced prior work and on the standard reduction of the PNP system to an effective diffusion equation in the quasi-neutral limit.

axioms (1)

domain assumption The quasi-neutral limit reduces the PNP system to a single diffusion equation with an effective diffusion coefficient.
Stated directly in the abstract when describing the limit behavior.

pith-pipeline@v0.9.0 · 5466 in / 1127 out tokens · 69065 ms · 2026-05-18T00:07:45.198952+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

44 extracted references · 44 canonical work pages · 1 internal anchor

[1]

U. M. Ascher, S. J. Ruuth, and R. J. Spiteri. Implicit-explicit Runge-Kutta methods for time- dependent partial differential equations.Applied Numerical Mathematics, 25(2-3):151–167, 1997

work page 1997
[2]

Self-regulated biological transportation structures with general entropy dissipation: 2D case and leaf-shaped domain

C. Astuto, P. Markowich, S. Portaro, and G. Russo. Self-regulated biological transportation structures with general entropy dissipation: 2D case and leaf-shaped domain.arXiv preprint arXiv:2408.15680, 2024. 20 Clarissa Astuto [0000−0002−9417−2959]

work page arXiv 2024
[3]

A comparison of the coco-russo scheme and ghost-fem for elliptic equations in arbitrary domains

C. Astuto, A. Coco, and U. Zerbinati. A comparison of the Coco-Russo scheme and ghost-FEM for elliptic equations in arbitrary domains.arXiv preprint arXiv:2405.16582, 2024

work page arXiv 2024
[4]

Asymptotic Preserving and Accurate scheme for Multiscale Poisson-Nernst-Planck (MPNP) system

C. Astuto and G. Russo. Asymptotic Preserving and Accurate scheme for Multiscale Poisson- Nernst-Planck (MPNP) system.arXiv preprint arXiv:2507.01402, 2025

work page internal anchor Pith review Pith/arXiv arXiv 2025
[5]

C. Astuto. High order multiscale methods for advection-diffusion equation in highly oscillatory regimes: application to surfactant diffusion and generalization to arbitrary domains.Communica- tions in Computational Physics, 2024

work page 2024
[6]

Astuto, A

C. Astuto, A. Raudino, and G. Russo. Multiscale modeling of sorption kinetics.Multiscale Modeling & Simulation, 21(1):374–399, 2023

work page 2023
[7]

Raudino, A

A. Raudino, A. Grassi, G. Lombardo, G. Russo, C. Astuto, and M. Corti. Anomalous Sorption Kinetics of Self-Interacting Particles by a Spherical Trap.Communications in Computational Physics, 31(3):707–738, 2022

work page 2022
[8]

Astuto, A

C. Astuto, A. Coco, and G. Russo. A finite-difference ghost-point multigrid method for multi-scale modelling of sorption kinetics of a surfactant past an oscillating bubble.Journal of Computational Physics, 476:111880, 2023

work page 2023
[9]

Astuto, M

C. Astuto, M. Lemou, and G. Russo. Time multiscale modeling of sorption kinetics I: uniformly accurate schemes for highly oscillatory advection-diffusion equation.Multiscale Modeling & Simulation, 23(2):753—776, 2025

work page 2025
[10]

Raudino, D

A. Raudino, D. Raciti, A. Grassi, M. Pannuzzo, and M. Corti. Oscillations of Bubble Shape Cause Anomalous Surfactant Diffusion: Experiments, Theory, and Simulations.Langmuir, 2016

work page 2016
[11]

Astuto, J

C. Astuto, J. Haskovec, P. Markowich, and S. Portaro. Self-regulated biological transportation structures with general entropy dissipations, part I: The 1D case.Journal of Dynamics and Games, 2023

work page 2023
[12]

Astuto, D

C. Astuto, D. Boffi, G. Russo, and U. Zerbinati. A nodal ghost method based on variational formulation and regular square grid for elliptic problems on arbitrary domains in two space dimensions.Computer Methods in Applied Mechanics and Engineering, 443:118041, 2025, Elsevier

work page 2025
[13]

Boscarino, F

S. Boscarino, F. Filbet, and G. Russo. High order semi-implicit schemes for time dependent partial differential equations.Journal of Scientific Computing, 68:975–1001, 2016

work page 2016
[14]

Boscarino, L

S. Boscarino, L. Pareschi, and G. Russo. Implicit-explicit methods for evolutionary partial differ- ential equations.SIAM, 2024

work page 2024
[15]

Brezzi, M.-O

F. Brezzi, M.-O. Bristeau, L. P. Franca, M. Mallet, and G. Rog´e. A relationship between stabilized finite element methods and the Galerkin method with bubble functions.Computer Methods in Applied Mechanics and Engineering, 96(1):117–129, 1992

work page 1992
[16]

F. F. Chen.Introduction to Plasma Physics. Springer, 2012

work page 2012
[17]

Christie, D

I. Christie, D. F. Griffiths, A. R. Mitchell, and O. C. Zienkiewicz. Finite element methods for second order differential equations with significant first derivatives.International Journal for Numerical Methods in Engineering, 10(6):1389–1396, 1976

work page 1976
[18]

High order asymptotic preserving penalized numerical schemes for the euler-poisson system in the quasi-neutral limit

N. Crouseilles, G. Dimarco, and S. Samantaray. High order Asymptotic Preserving penal- ized numerical schemes for the Euler-Poisson system in the quasi-neutral limit.arXiv preprint arXiv:2409.04807, 2024

work page arXiv 2024
[19]

P. Degond. Asymptotic-Preserving Schemes for Fluid Models of Plasmas.Numerical Models for Fusion, Panoramas et Synth`eses, 39-40:1–90, 2013. Soci´et´e Math´ematique de France

work page 2013
[20]

Multigrid methods for the ghost finite element approximation of elliptic problems

H. Dilip and A. Coco. Multigrid methods for the ghost finite element approximation of elliptic problems.arXiv preprint arXiv:2505.051052025

work page arXiv
[21]

Eisenberg and W

B. Eisenberg and W. Liu. Poisson–Nernst–Planck systems for ion channels with permanent charges. SIAM Journal on Mathematical Analysis, 2007

work page 2007
[22]

J. C. Heinrich, P. S. Huyakorn, O. C. Zienkiewicz, and A. R. Mitchell. An “upwind” finite element scheme for two-dimensional convective transport equation.International Journal for Numerical Methods in Engineering, 11(1):131–143, 1977

work page 1977
[23]

J. C. Heinrich and O. C. Zienkiewicz, Quadratic finite element schemes for two-dimensional convective-transport problems,International Journal for Numerical Methods in Engineering, 11(12):1831–1844, 1977

work page 1977
[24]

L. J. Henderson. The fitness of the environment, an inquiry into the biological significance of the properties of matter.The American Naturalist, 47(554):105–115, 1913. Title Suppressed Due to Excessive Length 21

work page 1913
[25]

H. K. Henisch. Semiconductor contacts: an approach to ideas and models. 1984

work page 1984
[26]

S. Jin. Efficient asymptotic-preserving (AP) schemes for some multiscale kinetic equations.SIAM Journal on Scientific Computing, 21(2):441–454, 1999

work page 1999
[27]

J¨ ungel and Y.-J

A. J¨ ungel and Y.-J. Peng. A hierarchy of hydrodynamic models for plasmas. Quasi-neutral limits in the drift-diffusion equations.Asymptotic Analysis, 28, 2000

work page 2000
[28]

Lu and Y

B. Lu and Y. C. Zhou. Poisson–Nernst–Planck equations for simulating biomolecular diffusion- reaction processes II: Size effects on ionic distributions and diffusion-reaction rates.Biophysical Journal, 2011

work page 2011
[29]

Lundstrom

M. Lundstrom. Fundamentals of carrier transport, 2nd edn.Measurement Science and Technology, 13(2):230, 2002

work page 2002
[30]

Gasser, C

I. Gasser, C. D. Levermore, P. A. Markowich, and C. Schmeiser. The initial time layer problem and the quasineutral limit in the semiconductor drift-diffusion model.European Journal of Applied Mathematics, 12(4):497–512, 2001

work page 2001
[31]

P. A. Markowich, C. A. Ringhofer, and C. Schmeiser,Semiconductor Equations, Springer Science & Business Media, 2012

work page 2012
[32]

E. A. Mason and E. W. McDaniel.Transport Properties of Ions in Gases. Wiley, 1988

work page 1988
[33]

Nastasi and V

G. Nastasi and V. Romano. An efficient GFET structure.IEEE Transactions on Electron Devices, 68(9):4729–4734, 2021

work page 2021
[34]

Nastasi and V

G. Nastasi and V. Romano. A full coupled drift-diffusion-Poisson simulation of a GFET.Com- munications in Nonlinear Science and Numerical Simulation, 87:105300, 2020

work page 2020
[35]

Osher and R

S. Osher and R. P. Fedkiw.Level Set Methods and Dynamic Implicit Surfaces. Springer, Applied Mathematical Sciences, vol. 153, 2002

work page 2002
[36]

Pareschi and G

L. Pareschi and G. Russo. Implicit-explicit Runge-Kutta schemes for stiff systems of differential equations.Recent Trends in Numerical Analysis, 3:269–289, 2000

work page 2000
[37]

Pareschi and G

L. Pareschi and G. Russo, High order asymptotically strong-stability-preserving methods for hyperbolic systems with stiff relaxation, inHyperbolic Problems: Theory, Numerics, Applications, Proc. 9th Int. Conf. on Hyperbolic Problems, CalTech, Pasadena, March 25–29, 2002, pp. 241–251, Springer, 2003

work page 2002
[38]

Pareschi and G

L. Pareschi and G. Russo, Implicit–explicit Runge–Kutta schemes and applications to hyperbolic systems with relaxation,Journal of Scientific Computing, 25:129–155, 2005, Springer

work page 2005
[39]

Russo and P

G. Russo and P. Smereka. A remark on computing distance functions.Journal of Computational Physics, 163(1):51–67, 2000

work page 2000
[40]

J. A. Sethian.Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science. Cambridge Monographs on Applied and Computational Mathematics 3, 2nd edition, Cambridge University Press, 1999

work page 1999
[41]

R. M. Smith. Finite element solutions of the energy equation at high Peclet number.Computers & Fluids, 8(3):335–350, 1980

work page 1980
[42]

Sussman, P

M. Sussman, P. Smereka, and S. Osher. A level set approach for computing solutions to incom- pressible two-phase flow.Journal of Computational Physics, 114(1):146–159, 1994

work page 1994
[43]

Wesseling.Principles of Computational Fluid Dynamics

P. Wesseling.Principles of Computational Fluid Dynamics. Springer, 2001

work page 2001
[44]

The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam

Richardson, Lewis Fry.IX. The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam. Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character, 1911

work page 1911

[1] [1]

U. M. Ascher, S. J. Ruuth, and R. J. Spiteri. Implicit-explicit Runge-Kutta methods for time- dependent partial differential equations.Applied Numerical Mathematics, 25(2-3):151–167, 1997

work page 1997

[2] [2]

Self-regulated biological transportation structures with general entropy dissipation: 2D case and leaf-shaped domain

C. Astuto, P. Markowich, S. Portaro, and G. Russo. Self-regulated biological transportation structures with general entropy dissipation: 2D case and leaf-shaped domain.arXiv preprint arXiv:2408.15680, 2024. 20 Clarissa Astuto [0000−0002−9417−2959]

work page arXiv 2024

[3] [3]

A comparison of the coco-russo scheme and ghost-fem for elliptic equations in arbitrary domains

C. Astuto, A. Coco, and U. Zerbinati. A comparison of the Coco-Russo scheme and ghost-FEM for elliptic equations in arbitrary domains.arXiv preprint arXiv:2405.16582, 2024

work page arXiv 2024

[4] [4]

Asymptotic Preserving and Accurate scheme for Multiscale Poisson-Nernst-Planck (MPNP) system

C. Astuto and G. Russo. Asymptotic Preserving and Accurate scheme for Multiscale Poisson- Nernst-Planck (MPNP) system.arXiv preprint arXiv:2507.01402, 2025

work page internal anchor Pith review Pith/arXiv arXiv 2025

[5] [5]

C. Astuto. High order multiscale methods for advection-diffusion equation in highly oscillatory regimes: application to surfactant diffusion and generalization to arbitrary domains.Communica- tions in Computational Physics, 2024

work page 2024

[6] [6]

Astuto, A

C. Astuto, A. Raudino, and G. Russo. Multiscale modeling of sorption kinetics.Multiscale Modeling & Simulation, 21(1):374–399, 2023

work page 2023

[7] [7]

Raudino, A

A. Raudino, A. Grassi, G. Lombardo, G. Russo, C. Astuto, and M. Corti. Anomalous Sorption Kinetics of Self-Interacting Particles by a Spherical Trap.Communications in Computational Physics, 31(3):707–738, 2022

work page 2022

[8] [8]

Astuto, A

C. Astuto, A. Coco, and G. Russo. A finite-difference ghost-point multigrid method for multi-scale modelling of sorption kinetics of a surfactant past an oscillating bubble.Journal of Computational Physics, 476:111880, 2023

work page 2023

[9] [9]

Astuto, M

C. Astuto, M. Lemou, and G. Russo. Time multiscale modeling of sorption kinetics I: uniformly accurate schemes for highly oscillatory advection-diffusion equation.Multiscale Modeling & Simulation, 23(2):753—776, 2025

work page 2025

[10] [10]

Raudino, D

A. Raudino, D. Raciti, A. Grassi, M. Pannuzzo, and M. Corti. Oscillations of Bubble Shape Cause Anomalous Surfactant Diffusion: Experiments, Theory, and Simulations.Langmuir, 2016

work page 2016

[11] [11]

Astuto, J

C. Astuto, J. Haskovec, P. Markowich, and S. Portaro. Self-regulated biological transportation structures with general entropy dissipations, part I: The 1D case.Journal of Dynamics and Games, 2023

work page 2023

[12] [12]

Astuto, D

C. Astuto, D. Boffi, G. Russo, and U. Zerbinati. A nodal ghost method based on variational formulation and regular square grid for elliptic problems on arbitrary domains in two space dimensions.Computer Methods in Applied Mechanics and Engineering, 443:118041, 2025, Elsevier

work page 2025

[13] [13]

Boscarino, F

S. Boscarino, F. Filbet, and G. Russo. High order semi-implicit schemes for time dependent partial differential equations.Journal of Scientific Computing, 68:975–1001, 2016

work page 2016

[14] [14]

Boscarino, L

S. Boscarino, L. Pareschi, and G. Russo. Implicit-explicit methods for evolutionary partial differ- ential equations.SIAM, 2024

work page 2024

[15] [15]

Brezzi, M.-O

F. Brezzi, M.-O. Bristeau, L. P. Franca, M. Mallet, and G. Rog´e. A relationship between stabilized finite element methods and the Galerkin method with bubble functions.Computer Methods in Applied Mechanics and Engineering, 96(1):117–129, 1992

work page 1992

[16] [16]

F. F. Chen.Introduction to Plasma Physics. Springer, 2012

work page 2012

[17] [17]

Christie, D

I. Christie, D. F. Griffiths, A. R. Mitchell, and O. C. Zienkiewicz. Finite element methods for second order differential equations with significant first derivatives.International Journal for Numerical Methods in Engineering, 10(6):1389–1396, 1976

work page 1976

[18] [18]

High order asymptotic preserving penalized numerical schemes for the euler-poisson system in the quasi-neutral limit

N. Crouseilles, G. Dimarco, and S. Samantaray. High order Asymptotic Preserving penal- ized numerical schemes for the Euler-Poisson system in the quasi-neutral limit.arXiv preprint arXiv:2409.04807, 2024

work page arXiv 2024

[19] [19]

P. Degond. Asymptotic-Preserving Schemes for Fluid Models of Plasmas.Numerical Models for Fusion, Panoramas et Synth`eses, 39-40:1–90, 2013. Soci´et´e Math´ematique de France

work page 2013

[20] [20]

Multigrid methods for the ghost finite element approximation of elliptic problems

H. Dilip and A. Coco. Multigrid methods for the ghost finite element approximation of elliptic problems.arXiv preprint arXiv:2505.051052025

work page arXiv

[21] [21]

Eisenberg and W

B. Eisenberg and W. Liu. Poisson–Nernst–Planck systems for ion channels with permanent charges. SIAM Journal on Mathematical Analysis, 2007

work page 2007

[22] [22]

J. C. Heinrich, P. S. Huyakorn, O. C. Zienkiewicz, and A. R. Mitchell. An “upwind” finite element scheme for two-dimensional convective transport equation.International Journal for Numerical Methods in Engineering, 11(1):131–143, 1977

work page 1977

[23] [23]

J. C. Heinrich and O. C. Zienkiewicz, Quadratic finite element schemes for two-dimensional convective-transport problems,International Journal for Numerical Methods in Engineering, 11(12):1831–1844, 1977

work page 1977

[24] [24]

L. J. Henderson. The fitness of the environment, an inquiry into the biological significance of the properties of matter.The American Naturalist, 47(554):105–115, 1913. Title Suppressed Due to Excessive Length 21

work page 1913

[25] [25]

H. K. Henisch. Semiconductor contacts: an approach to ideas and models. 1984

work page 1984

[26] [26]

S. Jin. Efficient asymptotic-preserving (AP) schemes for some multiscale kinetic equations.SIAM Journal on Scientific Computing, 21(2):441–454, 1999

work page 1999

[27] [27]

J¨ ungel and Y.-J

A. J¨ ungel and Y.-J. Peng. A hierarchy of hydrodynamic models for plasmas. Quasi-neutral limits in the drift-diffusion equations.Asymptotic Analysis, 28, 2000

work page 2000

[28] [28]

Lu and Y

B. Lu and Y. C. Zhou. Poisson–Nernst–Planck equations for simulating biomolecular diffusion- reaction processes II: Size effects on ionic distributions and diffusion-reaction rates.Biophysical Journal, 2011

work page 2011

[29] [29]

Lundstrom

M. Lundstrom. Fundamentals of carrier transport, 2nd edn.Measurement Science and Technology, 13(2):230, 2002

work page 2002

[30] [30]

Gasser, C

I. Gasser, C. D. Levermore, P. A. Markowich, and C. Schmeiser. The initial time layer problem and the quasineutral limit in the semiconductor drift-diffusion model.European Journal of Applied Mathematics, 12(4):497–512, 2001

work page 2001

[31] [31]

P. A. Markowich, C. A. Ringhofer, and C. Schmeiser,Semiconductor Equations, Springer Science & Business Media, 2012

work page 2012

[32] [32]

E. A. Mason and E. W. McDaniel.Transport Properties of Ions in Gases. Wiley, 1988

work page 1988

[33] [33]

Nastasi and V

G. Nastasi and V. Romano. An efficient GFET structure.IEEE Transactions on Electron Devices, 68(9):4729–4734, 2021

work page 2021

[34] [34]

Nastasi and V

G. Nastasi and V. Romano. A full coupled drift-diffusion-Poisson simulation of a GFET.Com- munications in Nonlinear Science and Numerical Simulation, 87:105300, 2020

work page 2020

[35] [35]

Osher and R

S. Osher and R. P. Fedkiw.Level Set Methods and Dynamic Implicit Surfaces. Springer, Applied Mathematical Sciences, vol. 153, 2002

work page 2002

[36] [36]

Pareschi and G

L. Pareschi and G. Russo. Implicit-explicit Runge-Kutta schemes for stiff systems of differential equations.Recent Trends in Numerical Analysis, 3:269–289, 2000

work page 2000

[37] [37]

Pareschi and G

L. Pareschi and G. Russo, High order asymptotically strong-stability-preserving methods for hyperbolic systems with stiff relaxation, inHyperbolic Problems: Theory, Numerics, Applications, Proc. 9th Int. Conf. on Hyperbolic Problems, CalTech, Pasadena, March 25–29, 2002, pp. 241–251, Springer, 2003

work page 2002

[38] [38]

Pareschi and G

L. Pareschi and G. Russo, Implicit–explicit Runge–Kutta schemes and applications to hyperbolic systems with relaxation,Journal of Scientific Computing, 25:129–155, 2005, Springer

work page 2005

[39] [39]

Russo and P

G. Russo and P. Smereka. A remark on computing distance functions.Journal of Computational Physics, 163(1):51–67, 2000

work page 2000

[40] [40]

J. A. Sethian.Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science. Cambridge Monographs on Applied and Computational Mathematics 3, 2nd edition, Cambridge University Press, 1999

work page 1999

[41] [41]

R. M. Smith. Finite element solutions of the energy equation at high Peclet number.Computers & Fluids, 8(3):335–350, 1980

work page 1980

[42] [42]

Sussman, P

M. Sussman, P. Smereka, and S. Osher. A level set approach for computing solutions to incom- pressible two-phase flow.Journal of Computational Physics, 114(1):146–159, 1994

work page 1994

[43] [43]

Wesseling.Principles of Computational Fluid Dynamics

P. Wesseling.Principles of Computational Fluid Dynamics. Springer, 2001

work page 2001

[44] [44]

The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam

Richardson, Lewis Fry.IX. The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam. Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character, 1911

work page 1911