Pith. sign in

REVIEW 2 minor 47 references

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

T0 review · 0 major / 2 minor · reviewed 2026-05-24 · grok-4.3

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Four finite volume schemes for the unipolar degenerated drift-diffusion model provide stability and existence, with convergence for two.

desk verdict Four flux formulations tailored to this degeneracy, with stability and existence for all four schemes plus convergence for two. read the letter →

arxiv 1907.11126 v2 pith:ALHCY3RS submitted 2019-07-25 math.NA cs.NAmath.APphysics.comp-ph

Clément Cancès , Claire Chainais-Hillairet , Jürgen Fuhrmann , Benoît Gaudeul This is my paper · ORCID
classification math.NAcs.NAmath.APphysics.comp-ph
keywords finitevolumemethoddrift-diffusionequationsdegenerateddiffusionstabilityanalysisconvergenceproofnumericalschemeschemicalpotentialunipolarmodel
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

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.

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.

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.

Watch

Extended reading notes

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.

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.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

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 · score 0.0 of 10

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.

Assumptions & free parameters 0 free parameters · 2 assumptions · 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.

assumptions (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.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A numerical analysis focused comparison of several Finite Volume schemes for an Unipolar Degenerated Drift-Diffusion Model." pith.science (2026). https://pith.science/paper/ALHCY3RS

@misc{pith2026190711126,
  author       = {Pith},
  title        = {Pith review of: A numerical analysis focused comparison of several Finite Volume schemes for an Unipolar Degenerated Drift-Diffusion Model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ALHCY3RS}},
  note         = {Machine review of arXiv:1907.11126}
}
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.

Figures

Figures reproduced from arXiv: 1907.11126 by the authors.

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 N = {σ ⊂ Γ… view at source ↗
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) = +∞. 3.2. Uni… view at source ↗
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 10 3 10 2 … view at source ↗
Figures from the paper (10 more)
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 …
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 …
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]
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_FUL…
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 …
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…
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.00…
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 …
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:…
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 …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

47 extracted references · 47 canonical work pages

  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

Show all 47 references
  1. [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

    2018

  2. [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,

  3. [11]

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

  4. [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

    2016

  5. [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

    2017

  6. [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

    2018

  7. [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

    2003

  8. [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

    2005

  9. [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

    1999

  10. [18]

    Deimling.Nonlinear functional analysis

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

    1985

  11. [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

    2014

  12. [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

    2013

  13. [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

    2014

  14. [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

    2017

  15. [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

    2018

  16. [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

    2003

  17. [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

    2014

  18. [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

    2000

  19. [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

    2017

  20. [28]

    Fuhrmann

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

    2015

  21. [29]

    Fuhrmann

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

    2016

  22. [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

    2019 doi

  23. [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

    2019 doi

  24. [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 ...

  25. [33]

    Springer International Publishing

  26. [34]

    Gallouët

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

    2007

  27. [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

    2016

  28. [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

    2010

  29. [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

    2017

  30. [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

    1995

  31. [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

    2013

  32. [40]

    Leray and J

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

    1934

  33. [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

    2010

  34. [42]

    Revels, M

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

    2016 arXiv

  35. [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

    1969

  36. [44]

    Selberherr.Analysis and simulation of semiconductor devices

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

    2012

  37. [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,

  38. [46]

    Solid State Electrochem

    to appear in J. Solid State Electrochem

  39. [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. Lill...

    1998

Pith tools

Reviewed May 24, 2026 · model on record in the stance chip above.