pith. sign in

arxiv: 2604.27512 · v1 · submitted 2026-04-30 · 🧮 math.NA · cs.NA

Discontinuous Galerkin IMEX Pressure Correction Scheme for the Poisson-Nernst-Planck-Navier-Stokes Equations

Bikram Bir , Amiya K. Pani This is my paper

Pith reviewed 2026-05-07 09:52 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords discontinuous GalerkinIMEX schemepressure correctionPoisson-Nernst-PlanckNavier-Stokeserror estimatesmass conservationelectrokinetic flows
0
0 comments X

The pith

A discontinuous Galerkin IMEX pressure-correction scheme for the Poisson-Nernst-Planck-Navier-Stokes equations achieves optimal error estimates in L2 and energy norms along with discrete mass conservation for the ions.

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

The paper introduces a fully discrete numerical scheme that combines discontinuous Galerkin finite element methods for the spatial discretization with an implicit-explicit pressure correction approach for time integration applied to the coupled system of Poisson-Nernst-Planck and Navier-Stokes equations. It establishes optimal error bounds in the L2 norm and energy norms for the ion concentrations, electrostatic potential, fluid velocity, and pressure, in addition to proving that the discrete scheme conserves the mass of positive and negative ions separately. These results are important because they offer theoretical justification for using the method in simulations of electrokinetic flows where accuracy and conservation properties directly impact the reliability of predictions for physical systems like microfluidic devices or biological ion channels. Numerical experiments are included to validate the error orders and conservation properties.

Core claim

Based on a discontinuous Galerkin method in the spatial directions and an improved implicit-explicit pressure-correction scheme in the temporal direction, the fully discrete scheme for the Poisson-Nernst-Planck-Navier-Stokes equations achieves optimal error estimates in L2 and energy norms for the concentrations of positive and negative ions, the electrostatic potential, the fluid velocity, and the L2 norm of the fluid pressure. The discrete mass conservation properties of both ions are established.

What carries the argument

The combination of discontinuous Galerkin spatial discretization and improved implicit-explicit pressure-correction temporal discretization, which separates stiff and non-stiff terms to enable both error analysis and proof of discrete conservation laws.

If this is right

  • Optimal error estimates in L2 and energy norms guarantee that the discrete solutions converge to the exact solution at the expected rates under mesh refinement and time-step reduction.
  • Discrete mass conservation for each ion species ensures that the scheme respects the physical invariance of total ion count in the absence of external sources.
  • The pressure-correction splitting decouples the velocity-pressure solve from the ion transport, supporting stable long-time integration of the coupled system.
  • Numerical simulations confirm the theoretical error orders and conservation properties, validating the scheme for practical electrokinetic modeling.

Where Pith is reading between the lines

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

  • Similar analysis techniques could be applied to derive error bounds for extensions of the scheme to three space dimensions or variable-coefficient problems.
  • The stability restrictions on the time step may motivate development of adaptive or implicit-explicit variants to handle stiff regimes without severe step-size limits.
  • The conservation properties might make the scheme particularly suitable for long-time simulations of closed systems where global invariants control qualitative behavior.

Load-bearing premise

The error analysis assumes that the exact solution possesses sufficient regularity and that the time step satisfies stability restrictions arising from the IMEX pressure-correction splitting.

What would settle it

Numerical experiments on a problem with a known smooth exact solution that fail to exhibit the predicted optimal convergence rates in the L2 or energy norms, or that show ion masses deviating from their initial values beyond roundoff error, would falsify the claims.

Figures

Figures reproduced from arXiv: 2604.27512 by Amiya K. Pani, Bikram Bir.

Figure 5.1
Figure 5.1. Figure 5.1: Numerical errors with k = 1 for Example 5.1 with DG method. 0.0313 0.0625 0.125 0.25 0.5 h 10-7 10-6 10-5 10-4 10-3 10-2 10-1 100 Error e u -L2 e c 1 -L2 e c 2 -L2 e -L2 h 3 0.0313 0.0625 0.125 0.25 0.5 h 10-5 10-4 10-3 10-2 10-1 100 101 Error e u -H1 e c 1 -H1 e c 2 -H1 e -H1 e p -L2 h 2 view at source ↗
Figure 5.2
Figure 5.2. Figure 5.2: Numerical errors with k = 2 for Example 5.1 with DG method. In order to verify the convergence rates in temporal direction, we chose a uniform partition of the time interval [0, 0.1] with ∆t = 0.1 × 2 −i , i = 1, 2, . . . , 6. The mesh parameter h has to be fixed as O(∆t 1 k+1 ) for L 2 -norm and O(∆t 1 k ) for H1 -norm. The view at source ↗
Figure 5.3
Figure 5.3. Figure 5.3: Numerical errors with respect to time variable for Examp view at source ↗
Figure 5.4
Figure 5.4. Figure 5.4: Mass (left) and deviation of mass (right) for both the po view at source ↗
Figure 5.5
Figure 5.5. Figure 5.5: Minimum and maximum values at time t for both the positive and negative charge ions concen￾tration (left) and electric potential and total energy (right) of the system with DG method. positive profile on the lower left side of the channel and an equal negative profile on the upper right portion of the channel. The initial concentration profile is given by ci(x, 0) = r0 2πR2 exp ( − ((x − Lx 2 + qi 8 ) 2 … view at source ↗
Figure 5.6
Figure 5.6. Figure 5.6: Domain for ion spreading. Set up of the geometry, bound view at source ↗
Figure 5.7
Figure 5.7. Figure 5.7: The evolution of net charge of ion concentrations ( view at source ↗
Figure 5.8
Figure 5.8. Figure 5.8: The evolution of electrostatic potential ( view at source ↗
Figure 5.9
Figure 5.9. Figure 5.9: The evolution of streamlines of fluid velocity in different time view at source ↗
Figure 5.10
Figure 5.10. Figure 5.10: Mass (left) and deviation of mass (right) for both the p view at source ↗
Figure 5.11
Figure 5.11. Figure 5.11: Minimum and maximum values at time t for both the positive and negative charge ions con￾centration (left) and electric potential and total energy (right) of the system. [16] Y. He and H. Chen. Efficiently high-order time-stepping R-GSAV schemes for the Navier–Stokes– Poisson–Nernst–Planck equations. Physica D: Nonlinear Phenomena, 466:134233, 2024. [17] Y. He and H. Chen. Stability and error analysis of… view at source ↗
read the original abstract

Based on a discontinuous Galerkin method in the spatial directions and an improved implicit-explicit pressure-correction scheme in the temporal direction, this paper discusses a fully discrete scheme for the Poisson-Nernst-Planck-Navier-Stokes equations. Optimal error estimates are derived in $L^2$ and in the energy norms for the concentrations of positive and negative ions, the electrostatic potential, the fluid velocity, and the $L^2$ norm of the fluid pressure. The discrete mass conservation properties of both ions are established. Finally, numerical simulations are performed, whose results confirm our theoretical findings.

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.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a fully discrete scheme for the Poisson-Nernst-Planck-Navier-Stokes equations that combines a discontinuous Galerkin spatial discretization with an improved IMEX pressure-correction time integrator. It derives optimal L² and energy-norm error estimates for the positive and negative ion concentrations, electrostatic potential, fluid velocity, and the L² norm of the pressure; establishes discrete mass conservation for both ion species; and presents numerical simulations that are said to confirm the theoretical results.

Significance. A rigorous error analysis together with proven mass conservation for a high-order DG-IMEX scheme on this coupled electrohydrodynamic system would be a valuable contribution to the numerical analysis of multiphysics problems. Such results could guide the design of stable, accurate, and structure-preserving methods for applications in microfluidics and electrokinetics, provided the regularity assumptions can be justified or relaxed.

major comments (2)
  1. [error analysis section] The error analysis (presumably §4) establishes optimal rates only under the assumption that the exact solution lies in sufficiently high Sobolev spaces (typically H^{k+2} or equivalent) uniformly in time. For the PNP-NS system this regularity is not guaranteed near boundaries or at moderate-to-high Péclet numbers; the manuscript does not discuss how the constants in the error bounds behave when this assumption is violated, which directly affects the claimed optimality.
  2. [stability and error analysis] The IMEX pressure-correction splitting treats convective and nonlinear terms explicitly, inducing a time-step restriction for stability. No explicit statement or derivation of this restriction appears in the stability or error sections, nor is it verified numerically for the fully coupled system; without it the error estimates cannot be guaranteed to hold for arbitrary Δt.
minor comments (2)
  1. [abstract and introduction] The abstract and introduction would benefit from a brief statement of the precise DG polynomial degree k and the precise form of the pressure-correction splitting (e.g., which terms are treated implicitly versus explicitly).
  2. [preliminaries] Notation for the broken Sobolev spaces, the DG bilinear forms, and the projection operators should be collected in a single preliminary subsection for easier reference during the error analysis.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and indicate the changes we will make.

read point-by-point responses

  1. Referee: [error analysis section] The error analysis (presumably §4) establishes optimal rates only under the assumption that the exact solution lies in sufficiently high Sobolev spaces (typically H^{k+2} or equivalent) uniformly in time. For the PNP-NS system this regularity is not guaranteed near boundaries or at moderate-to-high Péclet numbers; the manuscript does not discuss how the constants in the error bounds behave when this assumption is violated, which directly affects the claimed optimality.

    Authors: We agree that the optimal error estimates are obtained under the standard assumption of sufficient Sobolev regularity of the exact solution, which is required to close the estimates at the stated rates. The constants in the bounds depend on these higher norms, so reduced regularity would generally yield suboptimal rates or larger prefactors. This is a common limitation in the analysis of nonlinear multiphysics systems. We will add a clarifying remark in the introduction and at the conclusion of the error analysis section that explicitly states the regularity hypotheses and notes their implications for the PNP-NS system, including possible degradation near boundaries or at elevated Péclet numbers. revision: partial

  2. Referee: [stability and error analysis] The IMEX pressure-correction splitting treats convective and nonlinear terms explicitly, inducing a time-step restriction for stability. No explicit statement or derivation of this restriction appears in the stability or error sections, nor is it verified numerically for the fully coupled system; without it the error estimates cannot be guaranteed to hold for arbitrary Δt.

    Authors: The referee is correct that the explicit treatment of convection induces a CFL-type restriction on Δt for stability. While the stability analysis in Section 3 implicitly requires a sufficiently small time step to control the explicit terms, we did not state the restriction explicitly. In the revised version we will derive the precise time-step condition from the stability estimates and insert it into the stability section. We will also add a short numerical check in Section 5 confirming that the time steps used in the reported experiments satisfy the derived restriction for the coupled system. revision: yes

Circularity Check

0 steps flagged

No circularity in derivation chain

full rationale

The paper defines a DG-IMEX pressure-correction scheme for the PNP-NS system and performs standard a priori error analysis to obtain optimal L2 and energy-norm bounds under explicit regularity assumptions on the exact solution. Discrete mass conservation follows directly from the weak formulation of the scheme by summation of test functions. Neither the error bounds nor the conservation properties reduce to self-definitions, fitted parameters renamed as predictions, or load-bearing self-citations; the derivation relies on external approximation theory for DG spaces and stability restrictions induced by the IMEX splitting, which are independent of the target results.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on standard assumptions from finite-element analysis for time-dependent PDEs; no free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption The exact solution of the continuous PNP-NS system possesses sufficient regularity to attain optimal convergence rates.
    Standard assumption invoked for deriving optimal error estimates in DG methods for evolutionary PDEs.

pith-pipeline@v0.9.0 · 5399 in / 1259 out tokens · 76593 ms · 2026-05-07T09:52:04.169709+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    D. N. Arnold. An interior penalty finite element method with discont inuous elements. SIAM J. Numer. Anal., 19(4):742–760, 1982

    work page 1982

  2. [2]

    Bajpai, D

    S. Bajpai, D. Goswami, and K. Ray. A priori error estimates of a d iscontinuous Galerkin method for the Navier-Stokes equations. Numer. Algorithms , 94(2):937–1002, 2023

    work page 2023

  3. [3]

    M. Z. Bazant, K. Thornton, and A. Ajdari. Diffuse-charge dyna mics in electrochemical systems. Phys. Rev. E. , 70(2):021506, 2004

    work page 2004

  4. [4]

    B. Bir, H. Hutridurga, and A. K. Pani. On a completely discrete disc ontinuous Galerkin method for incompressible Chemotaxis-Navier-Stokes Equations. J. Sci. Comput. , 105(2):1–41, 2025

    work page 2025

  5. [5]

    Bothe, A

    D. Bothe, A. Fischer, and J. Saal. Global well-posedness and sta bility of electrokinetic flows. SIAM J. Math. Anal. , 46(2):1263–1316, 2014

    work page 2014

  6. [6]

    P. G. Ciarlet. The Finite Element Method for Elliptic Problem . North-Holland Publ. Comp., Amsterdam, 1978

    work page 1978

  7. [7]

    Constantin, M

    P. Constantin, M. Ignatova, and F. N. Lee. Nernst–Planck–Na vier–Stokes system near equilibrium. Pure Appl. Func. Anal. , 7(1):175–196, 2022

    work page 2022

  8. [8]

    C. I. Correa, G. N. Gatica, E. Henr ´ ıquez, R. Ruiz-Baier, and M. Solano. Banach spaces-based mixed finite element methods for the coupled Navier–Stokes and Poisson– Nernst–Planck equations. Calcolo, 61(2):31, 2024

    work page 2024

  9. [9]

    D. A. Di Pietro and A. Ern. Mathematical aspects of discontinuous Galerkin methods , volume 69. Springer Science & Business Media, 2011. 26 0 0.2 0.4 0.6 0.8 10 0.5 1 1.5 2 c1-c2 7 6 5 4 3 2 1 0 -1 -2 -3 -4 -5 -6 -7 0 0.2 0.4 0.6 0.8 10 0.5 1 1.5 2 c1-c2 7 6 5 4 3 2 1 0 -1 -2 -3 -4 -5 -6 -7 0 0.2 0.4 0.6 0.8 10 0.5 1 1.5 2 c1-c2 6 5 4 3 2 1 0 -1 -2 -3 -4 ...

    work page 2011

  10. [10]

    Douglas and T

    J. Douglas and T. Dupont. Interior penalty procedures for ellip tic and parabolic Galerkin methods. In Computing methods in applied sciences (Second Internat. Sy mpos., Versailles, 1975) , Lecture Notes in Phys., Vol. 58, pages 207–216. Springer, Berlin-New York, 1976

    work page 1975

  11. [11]

    Girault and P

    V. Girault and P. A. Raviart. Finite element approximation of the Navier-Stokes equatio ns, volume 749. Springer Berlin, 1979. 27 0 0.2 0.4 0.6 0.8 10 0.5 1 1.5 2 4.5 4 3.5 3 2.5 2 1.5 1 0.5 0 -0.5 -1 -1.5 -2 -2.5 -3 -3.5 -4 0 0.2 0.4 0.6 0.8 10 0.5 1 1.5 2 4.5 4 3.5 3 2.5 2 1.5 1 0.5 0 -0.5 -1 -1.5 -2 -2.5 -3 -3.5 -4 0 0.2 0.4 0.6 0.8 10 0.5 1 1.5 2 4 3....

    work page 1979

  12. [12]

    Girault, B

    V. Girault, B. Rivi` ere, and M. F. Wheeler. A discontinuous Galer kin method with nonoverlapping domain decomposition for the Stokes and Navier-Stokes problems. Math. Comput. , 74(249):53–84, 2005

    work page 2005

  13. [13]

    Girault, B

    V. Girault, B. Rivi` ere, and M. F. Wheeler. A splitting method usin g discontinuous Galerkin for the transient incompressible Navier-Stokes equations. M2AN Math. Model. Numer. Anal. , 39(6):1115–1147, 28 0 0.2 0.4 0.6 0.8 10 0.5 1 1.5 2 Streamline 2.5E-02 2.3E-02 2.0E-02 1.8E-02 1.6E-02 1.3E-02 1.1E-02 8.8E-03 6.5E-03 4.2E-03 1.9E-03 -3.8E-04 -2.7E-03 -5...

    work page 2005

  14. [14]

    He and P

    M. He and P. Sun. Mixed finite element analysis for the Poisson–Ne rnst–Planck/Stokes coupling. J. Comput. Appl. Math. , 341:61–79, 2018

    work page 2018

  15. [15]

    He and P

    M. He and P. Sun. Mixed finite element method for modified Poisson –Nernst–Planck/Navier–Stokes equations. J. Sci. Comput. , 87(3):80, 2021. 29 0 1 2 3 4 5 6 7 8 9 10 Time 1.3586 1.35865 1.3587 1.35875 1.3588 1.35885 1.3589 1.35895 1.359 Mass c1(x,t) dx c10(x) dx c2(x,t) dx c20(x) dx 0.26 0.28 0.3 1.358841214042 1.358841214044 1.358841214046 Zoomed-In View ...

    work page 2021

  16. [16]

    He and H

    Y. He and H. Chen. Efficiently high-order time-stepping R-GSA V s chemes for the Navier–Stokes– Poisson–Nernst–Planck equations. Physica D: Nonlinear Phenomena , 466:134233, 2024

    work page 2024

  17. [17]

    He and H

    Y. He and H. Chen. Stability and error analysis of pressure-cor rection scheme for the Navier–Stokes– Planck–Nernst–Poisson equations. ESAIM: Math. Model. Numer. Anal. , 59(3):1471–1504, 2025

    work page 2025

  18. [18]

    F. Hecht. New development in freefem++. J. Numer. Math. , 20(3-4):251–265, 2012

    work page 2012

  19. [19]

    A. T. Hill and E. S¨ uli. Approximation of the global attractor for the incompressible Navier-Stokes equations. IMA J. Numer. Anal. , 20(4):633–667, 2000

    work page 2000

  20. [20]

    J. W. Jerome. Analytical approaches to charge transport in a moving medium. Transp. Theory Stat. Phys., 31(4-6):333–366, 2002

    work page 2002

  21. [21]

    F. N. Lee. Global regularity for Nernst–Planck–Navier–Stoke s systems with mixed boundary conditions. Nonlinearity, 36(1):255, 2022

    work page 2022

  22. [22]

    Li and Z

    M. Li and Z. Li. Error estimates for the finite element method of the Navier-Stokes-Poisson-Nernst- Planck equations. Appl. Numer. Math. , 197:186–209, 2024

    work page 2024

  23. [23]

    M. Li, G. Zou, and B. Wang. A new fully-decoupled energy-stable bdf2-fem scheme for the electro- hydrodynamic equations. Math. and Comput. Simul. , 239:172–191, 2026

    work page 2026

  24. [24]

    X. Li, J. Shen, and Z. Liu. New sav-pressure correction metho ds for the Navier-Stokes equations: stability and error analysis. Math. Comput. , 91(333):141–167, 2022

    work page 2022

  25. [25]

    Linga, A

    G. Linga, A. Bolet, and J. Mathiesen. Transient electrohydrod ynamic flow with concentration-dependent fluid properties: modelling and energy-stable numerical schemes. J. Comput. Phys. , 412:109430, 2020. 30

    work page 2020

  26. [26]

    J. G. Liu and J. Wang. Global existence for Nernst–Planck–Nav ier–Stokes system in Rn. Commun. Math. Sci. , 18(6), 2020

    work page 2020

  27. [27]

    Liu and C

    X. Liu and C. Xu. Efficient time-stepping/spectral methods for the Navier-Stokes-Nernst-Planck-Poisson equations. Commun. Comput. Phys. , 21(5):1408–1428, 2017

    work page 2017

  28. [28]

    Masri, C

    R. Masri, C. Liu, and B. Riviere. Improved a priori error estimat es for a discontinuous Galerkin pressure correction scheme for the Navier-Stokes equations. Numer. Methods Partial Differential Equations , 39(4):3108–3144, 2023

    work page 2023

  29. [29]

    M. Pan, S. Liu, W. Zhu, F. Jiao, and D. He. A linear, second-orde r accurate, positivity-preserving and unconditionally energy stable scheme for the Navier–Stokes–Poiss on–Nernst–Planck system. Commun. Nonlinear Sci. Numer. , 131:107873, 2024

    work page 2024

  30. [30]

    Prohl and M

    A. Prohl and M. Schmuck. Convergent finite element discretiza tionsof the Navier-Stokes-Nernst-Planck- Poisson system. ESAIM: Math. Model. Numer. Anal. , 44(3):531–571, 2010

    work page 2010

  31. [31]

    Qin and C

    Y. Qin and C. Wang. A second-order accurate, positivity-pres erving numerical scheme for the Poisson– Nernst–Planck–Navier–Stokes system. IMA J. Numer. Anal. , page draf027, 2025

    work page 2025

  32. [32]

    W. H. Reed and T. R. Hill. Triangular mesh methods for the neutro n transport equation. Technical report, Los Alamos Scientific Lab., N. Mex.(USA), 1973

    work page 1973

  33. [33]

    Rivi` ere

    B. Rivi` ere. Discontinuous Galerkin methods for solving elliptic and pa rabolic equations: theory and implementation. SIAM, 2008

    work page 2008

  34. [34]

    R. J. Ryham. Existence, uniqueness, regularity and long-term behavior for dissipative systems modeling electrohydrodynamics. arXiv preprint arXiv:0910.4973 , 2009

    work page arXiv 2009

  35. [35]

    M. Schmuck. Analysis of the Navier–Stokes–Nernst–Planck–P oisson system. Math. Models Methods Appl. Sci. , 19(06):993–1014, 2009

    work page 2009

  36. [36]

    M. Wang, G. Zou, B. Wang, and W. Zhao. Unconditionally energy- stable discontinuous Galerkin method for the chemo-repulsion-Navier-Stokes system. Comput. Math. Appl. , 150:132–155, 2023

    work page 2023

  37. [37]

    Zhang and Z

    Z. Zhang and Z. Yin. Global well-posedness for the Navier–Stok es–Nernst–Planck–Poisson system in dimension two. Appl. Math. Lett. , 40:102–106, 2015

    work page 2015

  38. [38]

    Zhao and Y

    J. Zhao and Y. Li. Global existence of large solutions for the thr ee-dimensional incompressible Navier– Stokes–Poisson–Nernst–Planck equations. Math. Method. Appl. Sci. , 47(15):11933–11952, 2024

    work page 2024

  39. [39]

    Zhou and C

    X. Zhou and C. Xu. Efficient time-stepping schemes for the Navie r-Stokes-Nernst-Planck-Poisson equa- tions. Comput. Phys. Commun. , 289:108763, 2023. 31

    work page 2023