pith. sign in

arxiv: 1907.11126 · v2 · pith:ALHCY3RSnew · submitted 2019-07-25 · 🧮 math.NA · cs.NA· math.AP· physics.comp-ph

A numerical analysis focused comparison of several Finite Volume schemes for an Unipolar Degenerated Drift-Diffusion Model

Cl\'ement Canc\`es , Claire Chainais-Hillairet , J\"urgen Fuhrmann , Beno\^it Gaudeul This is my paper

Pith reviewed 2026-05-24 16:04 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath.APphysics.comp-ph
keywords finite volume methoddrift-diffusion equationsdegenerated diffusionstability analysisconvergence proofnumerical schemeschemical potentialunipolar model
0
0 comments X

The pith

Four finite volume schemes for the unipolar degenerated drift-diffusion model provide stability and existence, with convergence for two.

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

The paper examines an unipolar degenerated drift-diffusion system defined by the relation between concentration c and chemical potential h(c) = log(c/(1-c)). Four finite volume schemes are developed, each using a different formulation of the numerical fluxes. Stability analysis and proofs of existence of solutions are provided for every scheme. Convergence with respect to the discretization parameters is established specifically for two of the schemes. The behavior of all schemes is illustrated through numerical experiments.

Core claim

The authors design four finite volume schemes based on four different formulations of the fluxes for the unipolar degenerated drift-diffusion system. They prove stability and existence results for all four schemes and establish convergence proofs for two of them with respect to the discretization parameters.

What carries the argument

Four different formulations of the fluxes used to define the finite volume schemes, chosen to preserve key properties like positivity from the continuous model.

Load-bearing premise

The continuous unipolar degenerated drift-diffusion system with the given relation for h(c) is well-posed, allowing the discrete schemes to inherit its structural properties on admissible meshes.

What would settle it

Finding a sequence of meshes and time steps where one of the schemes fails to produce a bounded positive solution or where the two convergent schemes do not approach the continuous solution.

Figures

Figures reproduced from arXiv: 1907.11126 by Beno\^it Gaudeul, Claire Chainais-Hillairet, Cl\'ement Canc\`es, J\"urgen Fuhrmann.

Figure 1
Figure 1. Figure 1: Illustration of an admissible mesh as in Definition 2. We denote by mK the d-dimensional Lebesgue measure of the control volume K. The set of the faces is partitioned into two subsets: the set Eint of the interior faces defined by Eint = {σ ∈ E | σ = K|L for some K, L ∈ T } , and the set Eext of the exterior faces defined by Eext = {σ ∈ E | σ ⊂ ∂Ω} , which can also be partitioned into E D = {σ ⊂ ΓD} and E … view at source ↗
Figure 2
Figure 2. Figure 2: Evolution of the face concentration C(cK, cL, ΦK, ΦL) as a function of the jump of the potential ΦL − ΦK for the choice cK = 0.3 and cL = 0.7. Lemma 3.2. The face dissipation functional defined by (3.4) and either (2.8), (2.9), (2.10) or (2.11) satisfies the following dissipation property: given δ ∈ (0, 1) and M ∈ R, there holds, for Ψ and Υ as defined in (3.6): limcL→0 Ψδ,M(cL) = +∞, limcL→1 Υδ,M(cL) = +∞… view at source ↗
Figure 3
Figure 3. Figure 3: Left: time evolution of solution on domain Ω = (0, 50) from constant initial value c = 1 2 with Dirichlet boundary conditions Φ(0) = 10, Φ(50) = 0, c dp = − 1 2 and homogeneous Neumann boundary conditions for c. Right: Evolution of the relative free energy according to (1.11). 0 20 40 x 10 4 10 3 10 2 10 1 10 0 10 1 10 2 10 3 10 4 t : min=-62.5000 max=-0.0000 65 52 39 26 13 0 13 26 39 52 65 0 20 40 x 10 4 … view at source ↗
Figure 4
Figure 4. Figure 4: Left: time evolution of solution on domain Ω = (0, 50) from constant initial value c = 0.3 with Dirichlet boundary conditions Φ(0) = 0, Φ(50) = 0, c dp = − 1 2 and homogeneous Neumann boundary conditions for c. Right: Evolution of the relative free energy according to (1.11). decreases during time evolution for all four schemes discussed in this paper. We also remark that for zero applied potential, the co… view at source ↗
Figure 5
Figure 5. Figure 5: Left: time evolution of solution on domain Ω = (0, 50) from constant initial value c = 0.7 with Dirichlet boundary conditions Φ(0) = 0, Φ(50) = 0, c dp = − 1 2 and homogeneous Neumann boundary conditions for c. Right: Evolution of the relative free energy according to (1.11). 10 4 10 3 10 2 10 1 h 10 5 10 4 10 3 10 2 10 1 10 0 L 2 Error Centered Activity Bess-Ch Sedan O(h 2 ) 10 4 10 3 10 2 10 1 h 10 3 10 … view at source ↗
Figure 6
Figure 6. Figure 6: Convergence behavior of the different schemes for the case depicted in [PITH_FULL_IMAGE:figures/full_fig_p020_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Stationary solution with Dirichlet boundary conditions for c and Φ. 10 4 10 3 10 2 10 1 h 10 5 10 4 10 3 10 2 10 1 10 0 10 1 L 2 Error Centered Activity Bess-Ch Sedan O(h 2 ) 10 4 10 3 10 2 10 1 h 10 2 10 1 10 0 H 1 Error Centered Activity Bess-Ch Sedan O(h) [PITH_FULL_IMAGE:figures/full_fig_p021_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Convergence behavior of the different schemes. Left: L 2 -error, right: H1 error. Correspondence to the equation in the paper: “Centered”: (2.8), “Sedan”: (2.9), “Activity”: (2.10), “Bess-Ch”: (2.11). 0.00 0.02 0.04 0.06 0.08 0.10 0.000 0.005 0.010 40 20 0 20 40 0 1 2 3 4 5 Centered Activity Bess-Ch Sedan [PITH_FULL_IMAGE:figures/full_fig_p021_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Discretization grid of refinement level nref = 1 (left) and corresponding I-V curves for different discretization schemes (right). following boundary conditions at the contacts:  Φ c  =  −5 1 2  at Γsource = (0, 0.2 · L) × H  Φ c  =  5 1 2  at Γdrain = (0.8 · L, L) × H  ∇Φ · n J · n  =  − 1 d (Φ − Ugate) 0  at Γgate = (0.3 · L, 0.7 · L) × H  ∇Φ · n J · n  =  0 0  at ∂Ω \ (Γgate ∪ Γsource ∪ … view at source ↗
Figure 10
Figure 10. Figure 10: Electrostatic potential (left) and concentration (right) for closed gate (Ugate = 50). 0.00 0.02 0.04 0.06 0.08 0.10 x 0.000 0.002 0.004 0.006 0.008 0.010 y : min=-10.0000 max=10.0000 25 20 15 10 5 0 5 10 15 20 25 0.00 0.02 0.04 0.06 0.08 0.10 x 0.000 0.002 0.004 0.006 0.008 0.010 y c: min=8.221e-02 max=1-2.827e-03 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Electrostatic potential (left) and concentration (right) for Ugate = 0). 0.00 0.02 0.04 0.06 0.08 0.10 x 0.000 0.002 0.004 0.006 0.008 0.010 y : min=-24.9376 max=10.0000 25 20 15 10 5 0 5 10 15 20 25 0.00 0.02 0.04 0.06 0.08 0.10 x 0.000 0.002 0.004 0.006 0.008 0.010 y c: min=5.000e-01 max=1-2.326e-13 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Electrostatic potential (left) and concentration (right) for open gate (Ugate = −50), with concentration in the channel reaching the saturation value 1. 10 4 10 3 10 2 hx 10 2 10 1 10 0 10 1 ||IV IVref||2 Centered Activity Bess-Ch Sedan O(h) O(h 2 ) [PITH_FULL_IMAGE:figures/full_fig_p022_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Convergence of the I-V curves calculated using the different discretization schemes. Here, Φgate ∈ (−50, 50) is the gate voltage, and d = 0.1 · H is the gate thickness. We introduce a slightly anisotropic rectangular grid nx × ny grid with nx = 10 × 2 nref and ny = 5 × 2 nref , where nref is the refinement level. Each cell in the rectangular grid is subdivided into two triangles, see [PITH_FULL_IMAGE:fig… view at source ↗
read the original abstract

In this paper, we consider an unipolar degenerated drift-diffusion system where the relation between the concentration of the charged species $c$ and the chemical potential $h$ is $h(c)=\log \frac{c}{1-c}$. We design four different finite volume schemes based on four different formulations of the fluxes. We provide a stability analysis and existence results for the four schemes. The convergence proof with respect to the discretization parameters is established for two of them. Numerical experiments illustrate the behaviour of the different schemes.

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

0 major / 2 minor

Summary. The paper considers the unipolar degenerated drift-diffusion system with the relation h(c) = log(c/(1-c)). It constructs four finite-volume schemes from distinct flux formulations, establishes stability and existence results for all four schemes, proves convergence with respect to discretization parameters for two of them, and presents numerical experiments comparing their behavior.

Significance. If the stability, existence, and convergence claims hold, the manuscript supplies a systematic comparison of flux formulations for a degenerate parabolic system arising in semiconductor or biological modeling. The convergence proofs for two schemes constitute a concrete contribution that can guide scheme selection; the numerical experiments provide practical illustration of the theoretical distinctions.

minor comments (2)
  1. The abstract states that convergence is proved for two schemes but does not identify which two or the precise sense of convergence (e.g., strong L^1 or weak-*); adding this information would improve readability.
  2. Notation for the admissible meshes and the discrete gradient operators should be introduced once in a dedicated preliminary section rather than repeated across scheme definitions.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of the manuscript, the clear summary of its contributions, and the recommendation for minor revision. No specific major comments were listed in the report, so there are no individual points requiring point-by-point response or manuscript changes.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper constructs four finite volume schemes from explicit flux formulations for the given degenerate drift-diffusion model with h(c) = log(c/(1-c)). Stability, existence, and convergence results follow from direct a-priori estimates, positivity preservation, and compactness arguments on admissible meshes. No load-bearing step reduces by definition to its inputs, no fitted parameters are relabeled as predictions, and no uniqueness or ansatz is imported via self-citation chains. The modeling premise (well-posedness of the continuous system) is an external assumption standard to the field and does not create internal circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard mathematical assumptions from parabolic PDE theory and finite-volume discretization; no free parameters are fitted and no new physical entities are postulated.

axioms (2)
  • domain assumption The continuous unipolar degenerated drift-diffusion system admits solutions satisfying 0 < c < 1
    Implicit when the model is introduced via the relation h(c)=log(c/(1-c)) and when discrete schemes are required to preserve these bounds.
  • standard math Meshes are admissible in the sense required by standard finite-volume convergence theory
    Invoked for the convergence statements of the two schemes.

pith-pipeline@v0.9.0 · 5640 in / 1458 out tokens · 45300 ms · 2026-05-24T16:04:21.537390+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

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

  1. [1]

    Ait Hammou Oulhaj, C

    A. Ait Hammou Oulhaj, C. Cancès, and C. Chainais-Hillairet. Numerical analysis of a nonlinearly stable and positive Control Volume Finite Element scheme for Richards equation with anisotropy. ESAIM: Mathematical Modelling and Numerical Analysis, 52(4):1532–1567, 2018

    work page 2018

  2. [2]

    Andreianov, C

    B. Andreianov, C. Cancès, and A. Moussa. A nonlinear time compactness result and applications to discretization of degenerate parabolic–elliptic PDEs.J. Funct. Anal., 273(12):3633–3670, 2017

    work page 2017

  3. [3]

    Bessemoulin-Chatard

    M. Bessemoulin-Chatard. A finite volume scheme for convection-diffusion equations with nonlinear diffusion derived from the Scharfetter-Gummel scheme.Numer. Math., 121(4):637–670, 2012

    work page 2012

  4. [4]

    Bessemoulin-Chatard, C

    M. Bessemoulin-Chatard, C. Chainais-Hillairet, and F. Filbet. On discrete functional inequalities for some finite volume schemes. IMA J. Numer. Anal., 35:1125–1149, 2015

    work page 2015

  5. [5]

    Bezanson, Al

    J. Bezanson, Al. Edelman, S. Karpinski, and Viral B Shah. Julia: A fresh approach to numerical computing.SIAM review, 59(1):65–98, 2017

    work page 2017

  6. [6]

    Blakemore

    J.S. Blakemore. Approximations for fermi-dirac integrals, especially the function f12 (η) used to describe electron density in a semiconductor.Solid-State Electronics, 25(11):1067–1076, 1982

    work page 1982

  7. [7]

    Brenner, C

    K. Brenner, C. Cancès, and D. Hilhorst. Finite volume approximation for an immiscible two-phase flow in porous media with discontinuous capillary pressure.Comput. Geosci., 17(3):573–597, 2013

    work page 2013

  8. [8]

    Brochard, J

    F. Brochard, J. Jouffroy, and P. Levinson. Polymer-polymer diffusion in melts.Macromolecules, 16(10):1638–1641, 1983

    work page 1983

  9. [9]

    C. Cancès. Energy stable numerical methods for porous media flow type problems. Oil & Gas Science and Technology-Rev.IFPEN, 73:1–18, 2018

    work page 2018

  10. [10]

    Cancès, C

    C. Cancès, C. Chainais-Hillairet, and S. Krell. Numerical analysis of a nonlinear free-energy diminishing Discrete Duality Finite Volume scheme for convection diffusion equations.Computational Methods in Applied Mathematics,

    work page

  11. [11]

    Special issue on ”Advanced numerical methods: recent developments, analysis and application”

    work page

  12. [12]

    Cancès and C

    C. Cancès and C. Guichard. Convergence of a nonlinear entropy diminishing Control Volume Finite Element scheme for solving anisotropic degenerate parabolic equations.Math. Comp., 85(298):549–580, 2016

    work page 2016

  13. [13]

    Cancès and C

    C. Cancès and C. Guichard. Numerical analysis of a robust free energy diminishing finite volume scheme for parabolic equations with gradient structure.Found. Comput. Math., 17(6):1525–1584, 2017

    work page 2017

  14. [14]

    Cancès, F

    C. Cancès, F. Nabet, and M. Vohralík. Convergence and a posteriori error analysis for energy-stable finite element approximations of degenerate parabolic equations. HAL: hal-01894884, 2018

    work page 2018

  15. [15]

    Chainais-Hillairet, J.-G

    C. Chainais-Hillairet, J.-G. Liu, and Y.-J. Peng. Finite volume scheme for multi-dimensional drift-diffusion equations and convergence analysis.ESAIM: M2AN, 37(2):319–338, 2003

    work page 2003

  16. [16]

    R.Coehoorn, W.F.Pasveer, P.A.Bobbert, andM.A.J.Michels.Charge-carrierconcentrationdependenceofthehopping mobility in organic materials with gaussian disorder.Physical Review B, 72(15):155206, 2005

    work page 2005

  17. [17]

    Coudière, J.-P

    Y. Coudière, J.-P. Vila, and P. Villedieu. Convergence rate of a finite volume scheme for a two dimensional convection- diffusion problem.ESAIM Mathematical Modelling and Numerical Analysis, 33:493–516, 1999

    work page 1999

  18. [18]

    Deimling.Nonlinear functional analysis

    K. Deimling.Nonlinear functional analysis. Springer-Verlag, Berlin, 1985

    work page 1985

  19. [19]

    Dreyer, C

    W. Dreyer, C. Guhlke, and M. Landstorfer. A mixture theory of electrolytes containing solvation effects. Electrochemistry Communications, 43:75–78, 2014. FINITE VOLUME SCHEMES FOR IONIC LIQUID 29

    work page 2014

  20. [20]

    Dreyer, C

    W. Dreyer, C. Guhlke, and R. Müller. Overcoming the shortcomings of the nernst–planck model.Physical Chemistry Chemical Physics, 15(19):7075–7086, 2013

    work page 2013

  21. [21]

    J. Droniou. Finite volume schemes for diffusion equations: introduction to and review of modern methods.Math. Models Methods Appl. Sci., 24(8):1575–1620, 2014

    work page 2014

  22. [22]

    Droniou and R

    J. Droniou and R. Eymard. The asymmetric gradient discretisation method. In C. Cancès and P. Omnes, editors,Finite volumes for complex applications VIII - methods and theoretical aspects, volume 199 ofSpringer Proc. Math. Stat., pages 311–319, Cham, 2017. Springer

    work page 2017

  23. [23]

    Droniou, R

    J. Droniou, R. Eymard, T. Gallouët, C. Guichard, and R. Herbin.The Gradient Discretisation Method, volume 42 of Mathématiques et Applications. Springer International Publishing, 2018

    work page 2018

  24. [24]

    Eymard and T

    R. Eymard and T. Gallouët. H-convergence and numerical schemes for elliptic problems.SIAM J. Numer. Anal., 41(2):539–562, 2003

    work page 2003

  25. [25]

    Eymard, T

    R. Eymard, T. Gallouët, C. Guichard, R. Herbin, and R. Masson. TP or not TP, that is the question.Comput. Geosci., 18:285–296, 2014

    work page 2014

  26. [26]

    Eymard, T

    R. Eymard, T. Gallouët, and R. Herbin. Finite volume methods. Ciarlet, P. G. (ed.) et al., in Handbook of numerical analysis. North-Holland, Amsterdam, pp. 713–1020, 2000

    work page 2000

  27. [27]

    Farrell, T

    P. Farrell, T. Koprucki, and J. Fuhrmann. Computational and analytical comparison of flux discretizations for the semiconductor device equations beyond Boltzmann statistics.J. Comput. Phys., 346:497–513, 2017

    work page 2017

  28. [28]

    Fuhrmann

    J. Fuhrmann. Comparison and numerical treatment of generalised Nernst–Planck models.Comput. Phys. Commun., 196:166–178, 2015

    work page 2015

  29. [29]

    Fuhrmann

    J. Fuhrmann. A numerical strategy for Nernst-Planck systems with solvation effect.Fuel Cells, 16, 12 2016

    work page 2016

  30. [30]

    Fuhrmann

    J. Fuhrmann. UnipolarDriftDiffusion.jl - Numerical examples for finite volume schemes for unipolar drift-diffusion problems, 2019. URL:https://github.com/j-fu/UnipolarDriftDiffusion.jl, DOI:10.5281/zenodo.3351467

    work page doi:10.5281/zenodo.3351467 2019

  31. [31]

    Fuhrmann

    J. Fuhrmann. VoronoiFVM.jl - Solver for coupled nonlinear partial differential equations based on the Voronoi finite volume method, 2019. URL:https://github.com/j-fu/VoronoiFVM.jl, DOI:10.5281/zenodo.3351449

    work page doi:10.5281/zenodo.3351449 2019

  32. [32]

    Fuhrmann and C

    J. Fuhrmann and C. Guhlke. A finite volume scheme for Nernst-Planck-Poisson systems with ion size and solvation effects. In C. Cancès & P. Omnes, editor,Finite Volumes for Complex Applications VIII - Hyperbolic, Elliptic and Parabolic Problems, volume 200 ofSpringer Proceedings in Mathematics & Statistics, pages 497–505, Lille, France,

    work page

  33. [33]

    Springer International Publishing

    work page

  34. [34]

    Gallouët

    T. Gallouët. Nonlinear methods for linear equations. InProceedings of Tamtam’07conference held in Tipaza, 2007

    work page 2007

  35. [35]

    Gavish and A

    N. Gavish and A. Yochelis. Theory of phase separation and polarization for pure ionic liquids.The journal of physical chemistry letters, 7(7):1121–1126, 2016

    work page 2016

  36. [36]

    Glitzky and J

    A. Glitzky and J. A. Griepentrog. Discrete Sobolev-Poincaré inequalities for Voronoi finite volume approximations. SIAM J. Numer. Anal., 48:372–391, 2010

    work page 2010

  37. [37]

    Hackbusch.Elliptic Differential Equations: Theory and Numerical Treatment

    W. Hackbusch.Elliptic Differential Equations: Theory and Numerical Treatment. Springer Series in Computational Mathematics 18. Springer-Verlag Berlin Heidelberg, 2nd edition, 2017

    work page 2017

  38. [38]

    R. Herbin. An error estimate for a finite volume scheme for a diffusion–convection problem on a triangular mesh. Numer. Methods Partial Differential Equations, 11(2):165–173, 1995

    work page 1995

  39. [39]

    Koprucki and K

    T. Koprucki and K. Gärtner. Discretization scheme for drift-diffusion equations with strong diffusion enhancement. Optical and Quantum Electronics, 45(7):791–796, 2013

    work page 2013

  40. [40]

    Leray and J

    J. Leray and J. Schauder. Topologie et équations fonctionnelles.Ann. Sci. École Norm. Sup. (3), 51:45–78, 1934

    work page 1934

  41. [41]

    Paasch and S

    G. Paasch and S. Scheinert. Charge carrier density of organics with gaussian density of states: analytical approximation for the Gauss–Fermi integral.Journal of Applied Physics, 107(10):104501, 2010

    work page 2010

  42. [42]

    Forward-Mode Automatic Differentiation in Julia

    J. Revels, M. Lubin, and T. Papamarkou. Forward-mode automatic differentiation in julia.arXiv:1607.07892 [cs.MS], 2016

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  43. [43]

    D. L. Scharfetter and H. K. Gummel. Large-signal analysis of a silicon read diode oscillator.Electron Devices, IEEE Transactions on, 16(1):64–77, 1969

    work page 1969

  44. [44]

    Selberherr.Analysis and simulation of semiconductor devices

    S. Selberherr.Analysis and simulation of semiconductor devices. Springer, 2012

    work page 2012

  45. [45]

    Vágner, C

    P. Vágner, C. Guhlke, V. Miloš, R. Müller, and J. Fuhrmann. A continuum model for yttria-stabilised zirconia incor- porating triple phase boundary, lattice structure and immobile oxide ions. WIAS preprint 2583, Weierstrass Institute,

    work page

  46. [46]

    Solid State Electrochem

    to appear in J. Solid State Electrochem

    work page

  47. [47]

    Yu and R

    Z. Yu and R. Dutton. SEDAN III. www-tcad.stanford.edu/tcad/programs/sedan3.html, 1998. Clément Cancès ( clement.cances@inria.fr): Inria, Univ. Lille, CNRS, UMR 8524 - Laboratoire Paul Painlevé, F-59000 Lille Claire Chainais-Hillairet ( claire.chanais@univ-lille.fr): Univ. Lille, CNRS, UMR 8524, Inria - Labo- ratoire Paul Painlevé, F-59000 Lille Jürgen Fuh...

    work page 1998