Pith. sign in

REVIEW 3 major objections 5 minor 1 cited by

Analysis of a Poisson-Nernst-Planck cross-diffusion system with steric effects

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read The paper proves that the Poisson–Nernst–Planck system with Lennard-Jones steric repulsion has global weak solutions, a weak–strong uniqueness property, and exponential decay to thermal equilibrium under pure Neumann boundary conditions.

desk verdict Solid entropy-method paper closing a real gap in PNP steric models; Theorem 3's proof has a misstated elliptic estimate that is fixable. read the letter →

arxiv 2411.17399 v1 pith:KUCEKQVP submitted 2024-11-26 math.AP

classification math.AP MSC 35K5135A0235B4035Q92
keywords Poisson–Nernst–Planckequationsstericeffectscross-diffusionentropymethodweak–stronguniquenessexponentialdecayLennard-Jonesforcemixedboundaryconditions
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

This paper studies a model of ion transport in crowded electrolytes, where finite ion sizes add a Lennard-Jones repulsion and turn the Nernst–Planck fluxes into a cross-diffusion system. Its central claim is that this system is well posed in a weak sense for any number of species and for mixed Dirichlet–Neumann boundary conditions on the electric potential: global weak solutions exist, and any weak solution coincides with a strong solution with the same initial data as long as the strong solution exists. In the pure Neumann case the paper additionally proves exponential convergence to the thermal equilibrium state $u_\infty$ with $\Phi_\infty=0$. A sympathetic reader would care because these results extend earlier existence theory that was restricted to two species, to Dirichlet-only boundary conditions, or to models without a self-consistent electric potential. The proof is carried by an entropy structure: a free energy whose dissipation controls the $H^1(\Omega)$ norms of the concentrations.

What carries the argument

The load-bearing object is the Boltzmann–Rao entropy $H_{BR}(u)$ (equation (7) in the paper), which combines the Boltzmann entropy of the concentrations, the electric energy of the potential, and the quadratic steric energy $\tfrac12\sum a_{ij}u_iu_j$. The system is a formal gradient flow for this entropy, with entropy variable $w_i = \sigma\log u_i + z_i\Phi + p_i(u)$. The key estimate (Lemma 7) shows that the entropy production $\sum_i\int u_i|\nabla w_i|^2$ dominates $4\sigma^2|\nabla\sqrt{u_i}|^2 + \alpha\sigma|\nabla u_i|^2$ up to lower-order terms, and the positivity of $\alpha$ is exactly the positive definiteness of $(a_{ij})$. This single estimate yields the uniform $H^1$ bounds, the compactness needed for existence, and the relative-entropy inequalities used for weak–strong uniqueness and exponential decay.

What would settle it

A direct calculation shows the key estimate fails without positive definiteness: for $a_{ij}=1$ for all $i,j$ and two species, taking $\nabla u_2=-\nabla u_1$ gives $\sum_{i,j}a_{ij}\nabla u_i\cdot\nabla u_j=0$ while $\sum_i|\nabla u_i|^2>0$, so no $\alpha>0$ can satisfy Lemma 7. The assumption is therefore sharp for the entropy method as written; a concrete way to test whether positive definiteness is truly necessary for well-posedness is to run the same fixed-point construction with a semidefinite rank-one matrix and see whether uniform $H^1$ bounds or global weak solutions persist.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that the entropy structure inherited from the free energy $H_{BR}(u) = \int_\Omega ( \sigma \sum_i u_i(\log u_i -1) + \tfrac12 |\nabla(\Phi-\Phi_D)|^2 + \tfrac12 \sum_{i,j} a_{ij}u_iu_j + \sum_i z_i u_i \Phi_D )\,dx$ is strong enough to control the cross-diffusion terms and the mixed boundary data simultaneously. Theorem 1 states that, under assumptions (A1)–(A4), there exists a global weak solution with nonnegative concentrations, $\sqrt{u_i}, u_i, \Phi \in L^2(0,T;H^1(\Omega))$, and the entropy inequality $H_{BR}(u(t)) + \int_0^t\int_\Omega u_i|\nabla w_i|^2 \le H_{BR}(u_0)$, where $w_i = \sigma \log u_i + z_i\Phi + p_i(u)$ is the entropy variable, with $p_i(u)=\sum_j a_{ij}u_j$. Theorem 2 gives weak–strong uniqueness via the relative entropy $H_R(u|\bar u)$. Theorem 3 shows that under pure Neumann conditions, $\|u(t)-u_\infty\|_{L^2} + \|\nabla\Phi(t)\|_{L^2} \le H_{BR}(u_0|u_\infty) e^{-\lambda\sigma t}$. The exponential decay rate is proportional to the bare diffusion coefficient $\sigma$, and it vanishes when $\sigma=0$, consistent with the loss of $H^1$ control noted in the analysis.

Load-bearing premise

The argument rests on the steric-repulsion matrix $(a_{ij})$ being positive definite; if that matrix is only semidefinite, the entropy method's central estimate fails and the paper's existence and decay conclusions are not established.

Editorial extensions

If this is right

  • Any number of ion species with local Lennard-Jones repulsion now falls under the same existence theory; the earlier restrictions to two species, to Dirichlet-only boundary conditions, or to fixed potentials are removed.
  • Weak solutions are unique whenever a strong solution exists, so any numerical scheme that converges to a weak solution of (1)–(4) converges to the unique solution selected by the data.
  • Under pure Neumann boundary conditions the approach to equilibrium is exponential with rate $\lambda\sigma$, so smaller bare diffusivity means slower relaxation; the rate degenerates as $\sigma\to0$.
  • The entropy inequality (11) is a ready-made a priori estimate for designing structure-preserving discretizations, which the paper states as future work.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: the positive-definiteness of $(a_{ij})$ is doing more than technical work—it selects the strictly parabolic regime. In the rank-one semidefinite case the system is hyperbolic-parabolic, and the paper's Remark 12 suggests a change of variables to a drift-diffusion equation for the total concentration; global weak solutions there would require a different notion of solution, pos
  • A testable extension: one could replace the constant matrix $(a_{ij})$ by concentration-dependent repulsion strengths $a_{ij}(u)$. The entropy method should survive as long as the quadratic form remains uniformly positive definite along the flow, but the gradient-flow identity and the invertibility of $u\mapsto w$ would need reworking.
  • Mixed boundary conditions are the physically realistic case for ion channels, and the numerical experiment shows convergence to a nonconstant steady state there; a Lyapunov argument targeting that non-equilibrium stationary state, rather than $\Phi_\infty=0$, would be the natural next step.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. The paper analyzes a transient n-species Poisson–Nernst–Planck system with a symmetric positive definite steric-interaction matrix, no-flux boundary conditions for the ion concentrations, and mixed Dirichlet–Neumann boundary conditions for the electric potential. The authors prove global existence of weak solutions for d ≤ 4 via the boundedness-by-entropy method with an implicit time discretization and higher-order regularization, prove a weak–strong uniqueness result by a relative-entropy argument, and claim exponential decay to the thermal equilibrium in the pure-Neumann case. A numerical experiment illustrates convergence to a nonconstant steady state under mixed boundary conditions.

Significance. If the results are correct, the paper fills a genuine gap: previous existence results covered either two-species models or models without potential terms, so the n-species mixed-boundary case is new. The entropy structure is elegant, and the entropy-production estimate in Lemma 7 is a clean key step. The paper is self-contained modulo standard tools (the boundedness-by-entropy framework and elliptic regularity), and it contains no fitted parameters. However, the proof of Theorem 3 contains two load-bearing gaps: the 'usual elliptic estimate' used in the decay proof is false as stated, and the final norm bound does not follow from the entropy decay that is actually proved. These issues appear repairable, but the theorem as stated and its proof need correction.

major comments (3)
  1. [Section 4, after Eq. (35)] The assertion after Eq. (35) that -2σ∫Ω(Σ_i z_i u_i)^2 dx ≤ -2σ∫Ω|∇Φ|^2 dx is not a valid elliptic estimate. For the pure-Neumann problem -ΔΦ = ρ with zero mean, the correct estimate is ‖∇Φ‖_{L²}^2 ≤ C_P^2 ‖ρ‖_{L²}^2, i.e. ∫ρ^2 dx ≥ C_P^{-2}∫|∇Φ|^2 dx, where C_P is the Poincaré constant. The constant cannot be omitted: on Ω = (0,L) with L > π, Φ = cos(πx/L) gives ρ = (π/L)^2 cos(πx/L) and ∫ρ^2/∫|∇Φ|^2 = (π/L)^2 < 1. Replacing the coefficient 2 by 2 C_P^{-2} repairs the Gronwall argument and keeps a positive decay rate, so the result is likely true, but the proof as written is incomplete.
  2. [Theorem 3 statement and final paragraph of Section 4] The proof establishes a differential inequality for the relative entropy H_BR(u(t)|u∞), yielding exponential decay of the entropy. But the theorem claims a linear bound for ‖u(t)-u∞‖_{L²} + ‖∇Φ(t)‖_{L²}. The entropy controls the square of the L² distance through the positive definite quadratic term and (1/2)‖∇Φ‖_{L²}^2, so Gronwall gives at best a square-root estimate with half the decay rate. For a small perturbation u0 = u∞ + εφ with ∫φ dx = 0 and ε << 1, the left-hand side is of order ε, while H_BR(u0|u∞) is of order ε², so the stated inequality cannot hold generically. The theorem statement and the final step of Section 4 need to be corrected, for example by replacing the right-hand side by a constant times √H_BR(u0|u∞) e^{-λσt/2}, unless an additional argument proves a stronger decay.
  3. [Section 3, Eq. (34)] The central inequality for weak-strong uniqueness is introduced by 'After some computations' and 'similar arguments.' Since Theorem 2 rests entirely on this display, the authors should either present the full computation or put it in an appendix; at minimum the cancellation of the mixed terms in identities (32)–(33) should be shown step by step so that the sign conventions can be verified.
minor comments (5)
  1. [Introduction, after Eq. (4)] The line 'ΓD ∪ ΓN = ∂Ω and ΓD ∩ ΓD = ∅' contains a typo: it should read 'ΓD ∩ ΓN = ∅'.
  2. [Theorem 2] Theorem 2 says both solutions satisfy the entropy inequality (11), but the proof uses the Rao-type inequality (30); the statement should refer to the correct entropy inequality.
  3. [Remark 13] Remark 13 explicitly leaves the limit δ → 0 'to the reader'; since no proof is given, this remark should be phrased as a heuristic discussion or an open problem rather than as a claimed extension.
  4. [Section 4, derivation of (35)] The parenthetical note that (35) does not follow directly from (11) is important; please give the limiting argument (or a citation) that justifies replacing H_BR(u0) by H_BR(u(s)) and passing to the limit.
  5. [Section 5] The numerical section describes convergence 'close to the equilibrium state' without a quantitative measure; a plot of the decay of the error or of the entropy would make the illustration more informative, though this is not required for the analytical claims.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: the existence, weak-strong uniqueness, and exponential-decay theorems are derived from the stated assumptions using standard entropy, compactness, and fixed-point arguments; self-citations are methodological references, not load-bearing inputs.

full rationale

The paper's derivation chain is self-contained. The existence proof (Theorem 1) starts from the approximate problem (14)-(15) and derives the approximate entropy inequality (Lemma 6), the entropy production estimate (Lemma 7), and uniform bounds (Lemma 8) directly from the equations; the entropy inequality (11) is a consequence, not an assumption. The invertibility of the entropy-variable mapping is proved in Lemma 4. Theorem 2 uses the relative entropy (8) and derives inequality (34) from the equations and the assumed strong solution; it does not invoke the desired uniqueness as an input. Theorem 3 differentiates the relative entropy (10), uses the already-proved entropy inequality, the positive definiteness of (aij) (an explicit Assumption (A2)), the Poisson equation, and standard logarithmic Sobolev and Poincare inequalities; the decay rate is existential, not fitted. The many self-citations ([8], [24], [25], etc.) concern established methods or standard inequalities and are not used to assert this paper's theorems. The only flagged mathematical concern is the 'usual elliptic estimate' -2 sigma integral (sum z_i u_i)^2 <= -2 sigma integral |grad Phi|^2 in Section 4, which appears to be a genuine inequality error; however, this is a correctness issue, not circularity, since the step is not equivalent to the target result by construction and the proof is repairable. Remarks 12 and 13 explicitly acknowledge deferred semidefinite and sigma = 0 cases; these are stated limitations, not circular reasoning.

Assumptions & free parameters 0 free parameters · 7 assumptions · 0 invented entities

The paper's central claims rest on model assumptions (cross-diffusion form, positive definite interaction matrix, boundary data conditions A4/A5) and standard PDE tools. No free parameters are fitted: sigma, zi, aij, Phi_D, u0 are physical inputs, and the decay rate lambda is an existential constant. No new entities are introduced.

assumptions (7)
  • domain assumption The model equations (1)-(4) with localized Lennard-Jones force adequately describe ion transport with steric effects.
    The paper takes the cross-diffusion PNP system as given from prior physics literature [21,30]; no derivation is provided.
  • domain assumption (aij) is symmetric and positive definite (Assumption A2).
    Essential for parabolicity and for the H1 estimates in Lemma 7; the semidefinite case is deferred to Remark 12.
  • domain assumption Assumption (A4): Phi_D extends to H1 cap Linfty with Delta Phi_D = 0 and nabla Phi_D . nu = 0 on Gamma_N.
    Needed to define the Boltzmann-Rao entropy and entropy variables; a technical condition on boundary data.
  • domain assumption Assumption (A5): H2 regularity for the mixed boundary Poisson problem.
    Used for the Rao entropy inequality (20) and weak-strong uniqueness; standard but not proved in the paper.
  • standard math Sobolev embedding H1 subset L4 for d <= 4.
    Used to bound ui grad wi in L4/3 and to justify test functions.
  • standard math Logarithmic Sobolev and Poincare-Wirtinger inequalities.
    Used in Theorem 3 for exponential decay; constants absorbed into lambda.
  • standard math Leray-Schauder fixed-point theorem, Lax-Milgram, Aubin-Lions lemma, discrete Gronwall, Stampacchia truncation.
    Standard tools for the existence proof.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Analysis of a Poisson-Nernst-Planck cross-diffusion system with steric effects." pith.science (2026). https://pith.science/paper/KUCEKQVP

@misc{pith2026241117399,
  author       = {Pith},
  title        = {Pith review of: Analysis of a Poisson-Nernst-Planck cross-diffusion system with steric effects},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KUCEKQVP}},
  note         = {Machine review of arXiv:2411.17399}
}
read the original abstract

A transient Poisson-Nernst-Planck system with steric effects is analyzed in a bounded domain with no-flux boundary conditions for the ion concentrations and mixed Dirichlet-Neumann boundary conditions for the electric potential. The steric repulsion of ions is modeled by a localized Lennard-Jones force, leading to cross-diffusion terms. The existence of global weak solutions, a weak--strong uniqueness property, and, in case of pure Neumann conditions, the exponential decay towards the thermal equilibrium state is proved. The main difficulties are the cross-diffusion terms and the different boundary conditions satisfied by the unknowns. These issues are overcome by exploiting the entropy structure of the equations and carefully taking into account the electric potential term. A numerical experiment illustrates the long-time behavior of the solutions when the potential satisfies mixed boundary conditions.

Figures

Figures reproduced from arXiv: 2411.17399 by the authors.

Figure 1
Figure 1. Ion concentrations u1 (left column), u2 (middle column), and u3 (right column) at time steps N = 0 (top row), N = 30 (middle row), and N = 380 (bottom row). [4] L. Chen, E. Daus, and A. J¨ungel. Rigorous mean-field limit and cross diffusion. Z. Angew. Math. Phys. 70 (2019), article no. 122, 21 pages. [5] X. Chen and A. J¨ungel. Weak–strong uniqueness of renormalized solutions to reaction–cross-diffusion systems. Mat… view at source ↗
Figure 2
Figure 2. Electric potential Φ at time steps N = 0 (left), N = 30 (middle), and N = 380 (right). [11] G. Galiano and V. Selgas. On a cross-diffusion segregation problem arising from a model of interacting particles. Nonlin. Anal. Real World Appl. 18 (2014), 34–49. [12] G. Galiano, S. Shmarev, and J. Velasco. Existence and multiplicity of solutions to a cell-growth contact inhibition problem. Discrete Cont. Dyn. Sys. 35 (2015)… view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Global Finite-Energy Weak Solutions and Sharp Entropy Decay for a Poisson-Nernst-Planck System with Interspecies Drag and Steric Effects

    math.AP 2026-07 accept novelty 7.0 of 10

    Global finite-energy weak solutions exist, and near equilibrium the optimal entropy-decay constant is an explicit spectral quantity, for a two-species Poisson–Nernst–Planck system with steric interactions and interspe...

Reference graph

Works this paper leans on

38 extracted references · 36 canonical work pages · cited by 1 Pith paper

  1. [22]

    C.-Y. Hsieh. Global existence of solutions for the Poisson–Nernst–Planck system with steric effects. Nonlin. Anal. Real World Appl. 50 (2019), 34–54

  2. [27]

    J¨ ungel, S

    A. J¨ ungel, S. Portisch, and A. Zurek. Nonlocal cross-diffusion systems for multi-species populations and networks. Nonlin. Anal. 219 (2022), no. 112800, 26 pages. ANALYSIS OF A POISSON–NERNST–PLANCK SYSTEM 29

  3. [28]

    J¨ ungel and M

    A. J¨ ungel and M. Vetter. A convergent entropy-dissipating BDF2 finite-volume scheme for a popula- tion cross-diffusion system. Comput. Meth. Appl. Math. 24 (2024), 725–746

  4. [1]

    Berendsen, M

    J. Berendsen, M. Burger, V. Ehrlacher, and J.-F. Pietschmann. Uniqueness of strong solutions and weak–strong stability in a system of cross-diffusion equations. J. Evol. Eqs. 20 (2020), 459–483

  5. [2]

    Bertsch, M

    M. Bertsch, M. Gurtin, D. Hilhorst, and L. Peletier. On interacting populations that disperse to avoid crowding: preservation of segregation. J. Math. Biol. 23 (1985), 1–13

  6. [3]

    Burger, M

    M. Burger, M. Di Francesco, J.-F. Pietschmann, and B. Schlake. Nonlinear cross-diffusion with size exclusion. SIAM J. Math. Anal. 42 (2010), 2842–2871. ANALYSIS OF A POISSON–NERNST–PLANCK SYSTEM 27 Figure 1. Ion concentrations u1 (left column), u2 (middle column), and u3 (right column) at time steps N = 0 (top row), N = 30 (middle row), and N = 380 (bottom row)

  7. [4]

    L. Chen, E. Daus, and A. J¨ ungel. Rigorous mean-field limit and cross diffusion. Z. Angew. Math. Phys. 70 (2019), article no. 122, 21 pages

  8. [5]

    Chen and A

    X. Chen and A. J¨ ungel. Weak–strong uniqueness of renormalized solutions to reaction–cross-diffusion systems. Math. Models Meth. Appl. Sci. 29 (2019), 237–270

Show all 38 references
  1. [6]

    D. Clark. Short proof of a discrete Gronwall inequality. Discrete Appl. Math. 16 (1987), 279–281

  2. [7]

    De Nitti and N

    N. De Nitti and N. Zamponi. Fractional cross-diffusion in a bounded domain: existence, weak–strong uniqueness, and long-time asymptotics. Preprint, 2024. arXiv:2407.19824

  3. [8]

    Dreher and A

    M. Dreher and A. J¨ ungel. Compact families of piecewise constant functions in Lp(0, T; B). Nonlin. Anal. 75 (2012), 3072–3077

  4. [9]

    Druet and A

    P.-E. Druet and A. J¨ ungel. Analysis of cross-diffusion systems for fluid mixtures driven by a pressure gradient. SIAM J. Math. Anal. 52 (2020), 2179–2197

  5. [10]

    Gajewski

    H. Gajewski. On existence, uniqueness and asymptotic behavior of solutions of the basic equations for carrier transport in semiconductors. Z. Angew. Math. Mech. 65 (1985), 101–108. 28 P. HIR VONEN AND A. J ¨UNGEL Figure 2. Electric potential Φ at time steps N = 0 (left), N = 3...

  6. [11]

    Galiano and V

    G. Galiano and V. Selgas. On a cross-diffusion segregation problem arising from a model of interacting particles. Nonlin. Anal. Real World Appl. 18 (2014), 34–49

  7. [12]

    Galiano, S

    G. Galiano, S. Shmarev, and J. Velasco. Existence and multiplicity of solutions to a cell-growth contact inhibition problem. Discrete Cont. Dyn. Sys. 35 (2015), 1479–1501

  8. [13]

    N. Gavish. Poisson–Nernst–Planck equations with steric effects – non-convexity and multiple station- ary solutions. Physica D 368 (2018), 50–65

  9. [14]

    N. Gavish. Poisson–Nernst–Planck equations with high-order steric effects. Physica D 411 (2020), article no. 132536, 10 pages

  10. [15]

    Gerstenmayer and A

    A. Gerstenmayer and A. J¨ ungel, Analysis of a degenerate parabolic cross-diffusion system for ion transport. J. Math. Anal. Appl. 461 (2018), 523–543

  11. [16]

    Gillespie

    D. Gillespie. A review of steric interactions of ions: Why some theories succeed and others fail to account for ion size. Microfluid Nanofluid 18 (2015), 717–738

  12. [17]

    Glitzky and R

    A. Glitzky and R. H¨ unlich. Global estimates and asymptotics for electro–reaction–diffusion systems in heterostructures. Appl. Anal. 66 (1997), 205–226

  13. [18]

    Giunta, T

    V. Giunta, T. Hillen, M. Lewis, and J. Potts. Positivity and global existence for nonlocal advection- diffusion models of interacting populations. Submitted for publication, 2023. arXiv:2312.09692

  14. [19]

    Hopf and M

    K. Hopf and M. Burger. On multi-species diffusion with size exclusion. Nonlin. Anal. 224 (2022), article no. 113092, 27 pages

  15. [20]

    Hopf and A

    K. Hopf and A. J¨ ungel. Convergence of a finite volume scheme and dissipative measure-valued– strong stability for a hyperbolic–parabolic cross-diffusion system. Submitted for publication, 2023. arXiv:2304.00787

  16. [21]

    Horng, T.-C

    T.-L. Horng, T.-C. Lin, C. Liu, and B. Eisenberg. PNP equations with steric effects: a model of ion flow through channels. J. Phys. Chem. B 116 (2012), 11422–11441

  17. [23]

    X. Huo, A. J¨ ungel, and A. Tzavaras. Weak–strong uniqueness for Maxwell–Stefan systems. SIAM J. Math. Anal. 54 (2022), 3215–3252

  18. [24]

    J¨ ungel

    A. J¨ ungel. The boundedness-by-entropy method for cross-diffusion systems. Nonlinearity 28 (2015), 1963–2001

  19. [25]

    J¨ ungel

    A. J¨ ungel. Entropy Methods for Diffusive Partial Differential Equations . BCAM Springer Briefs, Springer, Switzerland, 2016

  20. [26]

    J¨ ungel and A

    A. J¨ ungel and A. Massimini. Analysis of a Poisson–Nernst–Planck–Fermi system for charge transport in ion channels. J. Differ. Eqs. 395 (2024), 38–68

  21. [29]

    Kilic, M

    M. Kilic, M. Bazant, and A. Ajdari. Steric effects in the dynamics of electrolytes at large applied voltages. II. Modified Poisson–Nernst–Planck equations. Phys. Rev. E 75 (2007), article no. 021503, 11 pages

  22. [30]

    Lin and B

    T.-C. Lin and B. Eisenberg. A new approach to the Lennard–Jones potential and a new model: PNP-steric equations. Commun. Math. Sci. 12 (2014), 149–173

  23. [31]

    W. Nernst. Die elektromotorische Wirksamkeit der Jonen. Z. Phys. Chem. 4 (1889), 129–181

  24. [32]

    M. Planck. Ueber die Erregung von Electricit¨ at und W¨ arme in Electrolyten.Ann. Phys. Chem. 275 (1890), 161–186

  25. [33]

    Sch¨ oberl

    J. Sch¨ oberl. C++11 implementation of finite elements in NGSolve. ASC Report 30/2014, TU Wien, 2014, 23 pages. http://hdl.handle.net/20.500.12708/28346

  26. [34]

    Schuss, B

    Z. Schuss, B. Nadler, and B. Eisenberg. Derivation of Poisson and Nernst–Planck equations in a bath and channel from a molecular model. Phys. Rev. E 64 (2001), article no. 036116, 14 pages

  27. [35]

    Siddiqua, Z

    F. Siddiqua, Z. Wang, and S. Zhou. A modified Poisson–Nernst–Planck model with excluded volume effect: theory and numerical implementation. Commun. Math. Sci. 16 (2018), 251–271

  28. [36]

    Simpson, K

    M. Simpson, K. Landman, and B. Hughes. Multi-species simple exclusion processes. Physica A 388 (2009), 399–406

  29. [37]

    R. Temam. Infinite-Dimensional Dynamical Systems in Mechanics and Physics . 2nd edition. Springer, New York, 1997

  30. [38]

    Troianiello

    G. Troianiello. Elliptic Differential Equations and Obstacle Problems . Plenum Press, New York, 1987. Institute of Analysis and Scientific Computing, TU Wien, Wiedner Hauptstraße 8–10, 1040 Wien, Austria Email address : peter.hirvonen@tuwien.ac.at Institute of Analysis and Sci...

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.