Pith. sign in

REVIEW 4 major objections 4 minor 56 references

A posteriori existence for the Keller-Segel model via a finite volume - finite element scheme

T0 review · 4 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash

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

Pith's one-line read This paper shows that conditional a posteriori error estimates for a finite volume scheme can certify existence of weak solutions to the Keller-Segel system up to the simulated time.

desk verdict A serious conditional a posteriori existence framework for Keller-Segel, with real novelty in the residual analysis, but the advertised fully computable condition depends on an unproved vector maximum-norm bound and an unstated constant. read the letter →

arxiv 2509.17710 v2 pith:7BQAL5Q7 submitted 2025-09-22 math.NA cs.NA

Marc Hoffmann , Jan Giesselmann This is my paper · ORCID
classification math.NAcs.NA MSC 65M1565M0835K40
keywords conditionalaposteriorierrorestimatesKeller-Segelsystemfinitevolumeschemeexistenceweaksolutionsblow-upMorley-typereconstructionelementmethod
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 establishes conditional a posteriori error estimates for a finite volume-finite element discretization of the parabolic-elliptic Keller-Segel system, a model of chemotaxis. The central claim is that when a fully computable smallness condition holds for the reconstructed numerical solution, the estimator controls the error to the true weak solution in $L^\infty(0,T;L^2)$ and $L^2(0,T;H^1)$, and as a consequence the weak solution actually exists on $[0,T]$. This converts a good simulation into a rigorous existence proof for a model whose solutions can blow up in finite time. The authors provide two stability frameworks, one based on an Adjusted Generalized Gronwall Lemma and one based on a local-in-time continuation argument, and they show numerically that the second condition is satisfiable on longer time horizons in two dimensions. The practical reach is limited by the sharpness of explicit constants, especially in three dimensions.

What carries the argument

The load-bearing object is the reconstructed pair $(\tilde\rho,\tilde c)$: a Morley-type interpolant (a nonstandard $C^0$ finite element whose degrees of freedom are vertex values and face-normal flux integrals) of the cell-centered finite-volume density, linearly interpolated in time, together with $\tilde c$ defined as the exact elliptic solution with right-hand side $\tilde\rho$, so that no residual enters the elliptic equation. Stability is supplied by two frameworks: an Adjusted Generalized Gronwall Lemma (Proposition 3.1) and a local-in-time continuation argument, each converting a small residual into an error bound under a smallness condition. The non-computable term $\|\nabla\tilde c(t)\|_{L^\infty}$ is replaced by a computable upper bound using a $P_1$ finite element approximation $q_h$ plus a maximum-norm a posteriori estimator with explicit constant $C_\infty$ (Remark 3.4). The residual is decomposed into diffusive, temporal, and convective parts, each estimated elementwise using explicit constants from Poincaré, trace, Sobolev-embedding, and elliptic-regularity inequalities.

What would settle it

Run the two-dimensional manufactured example from Section 7.1 on successively finer meshes until condition (7) or (9) holds, and compare the certified bound ($\delta A E$ or $\Psi_{N_t-1}$) with the exactly known error $\sup_t\|\rho-\tilde\rho\|_{L^2}^2 + \int_0^T |\rho-\tilde\rho|_{H^1}^2\,dt$; a single case in which the condition holds and the true error exceeds the bound would contradict Theorem 6.7.

Watch

Extended reading notes

Core claim

The paper's central claim is Theorem 6.7: for the Morley-type reconstruction $(\tilde\rho,\tilde c)$ of the finite-volume/finite-element approximation, the full residual in $H^{-1}$ is bounded by a computable estimator that is a sum over mesh elements and time steps (18). If either the Gronwall-type condition (7) or the local-continuation condition (9) holds, then $\sup_{t\in[0,T]}\|\rho(t)-\tilde\rho(t)\|_{L^2}^2 + \int_0^T |\rho(t)-\tilde\rho(t)|_{H^1}^2\,dt$ is bounded by $\delta A E$ or by $\Psi_{N_t-1}$, respectively. By Proposition 2.5, that boundedness forces $T<T_{\max}$, so in the paper's phrasing '$(\rho,c)$ is a weak solution up to time $T$'. Existence is thereby inferred a posteriori from the numerical approximation rather than assumed a priori.

Load-bearing premise

The whole method depends on explicit, sufficiently sharp upper bounds for the Sobolev embedding constants, the elliptic regularity constant, and the maximum-norm constant $C_\infty$; in three dimensions the known bounds are too rough for the conditions to hold on feasible meshes, and without an explicit $C_\infty$ the term $\|\nabla\tilde c(t)\|_{L^\infty}$ in the stability coefficient is not computable.

Editorial extensions

If this is right

  • When condition (7) holds for a computed run, the error to the exact weak solution is bounded by $\delta A E$ and, by Corollary 3.6, the weak solution exists beyond the simulated time $T$.
  • When condition (9) has a root on every time step, the same conclusion holds with the local estimator $\Psi_{N_t-1}$; in the paper's two-dimensional tests this condition is satisfied on substantially longer time horizons than the Gronwall condition.
  • Because the residual estimator and the factor $A$ shrink under mesh refinement for a convergent scheme, the conditional estimates become available on sufficiently fine meshes.
  • In the three-dimensional manufactured test, the residual estimator scales approximately linearly in the mesh size, matching the order of the $H^1$ error that the estimator bounds.

Reading between the lines

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

  • Beyond the paper, condition (9) could serve as a run-time certificate in three dimensions: if an explicit $C_\infty$ and sharper Sobolev constants are supplied, verifying the condition on a feasible mesh would settle existence for concrete initial data whose long-time behavior is analytically open.
  • The same 'norm-blow-up criterion plus computable conditional estimator' template seems transferable to other nonlinear evolution equations with finite-time blow-up, such as semilinear heat equations, whenever a stability framework and a flux-preserving reconstruction are available.
  • A natural testable extension not implemented here is adaptive mesh refinement driven by the elementwise estimators $\Theta_{K,n}$; refining where the estimator is large may extend the certified time horizon more efficiently than uniform refinement.
  • Using higher-order reconstructions or sharper stability constants would likely enlarge the attainable horizons more than raw mesh refinement, since the conditions in (7) and (9) depend exponentially on the stability coefficient through $E$.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 4 minor

Summary. The manuscript derives two conditional a posteriori error estimates for a cell-centered finite volume / P1 finite element approximation of the parabolic-elliptic Keller-Segel system in two and three dimensions. The numerical density is reconstructed with a Morley-type interpolant, the chemical concentration with an elliptic reconstruction, and the residual is estimated in L2(0,T;H^-1). Two stability frameworks are combined with these residual bounds: one based on a Generalized Gronwall Lemma and one based on a local-in-time continuation argument. The central result, Theorem 6.7, states that if a fully computable condition involving the computable residual estimator and a maximum-norm bound on the chemical gradient is satisfied, then the a posteriori error estimator bounds the error to the exact solution, and consequently a weak solution exists up to time T. Numerical experiments in 2D compare the availability times of the two estimators, and experiments in 3D report convergence rates of the residual estimator and the reconstruction error under mesh refinement.

Significance. If the fully computable conditions are rigorously established, this is a valuable contribution to computer-assisted existence for a PDE with finite-time blow-up: the conditional inference from numerical behavior to existence is conceptually sound, and the paper correctly avoids circularity by computing the residual and the stability condition from the reconstruction alone. The paper also provides reproducible code and explicit comparisons of two stability frameworks. The residual estimates and the Morley-type reconstruction are worked out in considerable detail. However, the central claim of 'fully computable' conditions is not currently supported: the maximum-norm a posteriori bound for the vector-valued elliptic problem is asserted but not proved, the constant C∞ is never made explicit, and the unisolvence of the Morley-type element is left unproved. In addition, the numerical validations solve non-homogeneous manufactured problems whose source terms are not accounted for in the theoretical residual estimates. These are load-bearing gaps in the manuscript as it stands.

major comments (4)
  1. [Remark 3.4 / Theorem 6.7] The claim that ||∇c̄(t)||L∞ ≤ ||qh(t)||L∞ + C∞η∞(t) is a computable upper bound is not established. The vector-valued problem (I−Δ)∇c̄ = ∇ρ̄ with ∇c̄·n = 0 is not the scalar reaction-diffusion problem covered by [19, Lemma 6], and the manuscript does not specify the finite element space for qh, the boundary conditions imposed on qh, or the definition of the jumps of ∇qh on a polygonal domain. Componentwise application of a scalar maximum-norm estimator is delicate because the normal and tangential components of ∇c̄ satisfy different boundary conditions on flat boundary parts. Moreover, C∞ is never given; the statement that it 'can be estimated explicitly by investigating the proof of [19, Lemma 6]' is not an explicit constant. Since a(t) in conditions (7) and (9) contains ||qh(t)||L∞ + C∞η∞(t), these conditions are not fully computable as written, and the a posteriori existence conclusion in Theorem 6.7 lacks a verified premise.
  2. [Section 5, Lemma 5.1] The unisolvence of the Morley-type element (K, P_K, Σ_K) is asserted with only 'Proof ... omitted for brevity', and no reference is given that covers this exact triple. The reconstruction ρ̃ is defined through (12)-(13) as the unique element of U_h, so a failure of unisolvence would invalidate the definition of ρ̃ and hence all subsequent residual estimates. A complete proof or a precise reference for this specific element is required.
  3. [Sections 7.1-7.2, Eqs. (19)-(22)] The numerical experiments use manufactured solutions that solve non-homogeneous systems with nonzero source terms f and g, whereas Algorithm 4.3, the residual estimate (18), and Theorem 6.7 concern the homogeneous Keller-Segel system (KS.1)-(KS.3). The manuscript does not state how f and g are incorporated into the numerical scheme or into the residual estimator. Unless the residual and the stability framework are extended to include these source terms, the convergence results in Table 2 and the availability horizons in Table 1 do not provide evidence for Theorem 6.7 as stated.
  4. [Corollary 3.6 / Proposition 2.5] The existence inference in Corollary 3.6 is stated as a direct consequence of Proposition 2.5, but the proof is not given. Since Theorem 3.3 is a stability estimate for weak solutions that are already assumed to exist, the argument should explicitly apply the estimate on [0,T'] for all T'<Tmax, use the fact that the condition (7) and the constants remain controlled as T'→Tmax, and then pass to the limit to contradict the blow-up criterion. This is likely fixable, but it is not written out and is needed for the conditional existence claim.
minor comments (4)
  1. [Abstract] The abstract contains a typo: 'a weak solution exits' should read 'a weak solution exists'.
  2. [Section 6.2, Theorem 6.7(ii)] In the final display of part (ii), the error bound is written with ∥ρ(t,·)−¯ρ(t,·)∥, but the reconstruction is denoted ˜ρ throughout; this appears to be a typo and should be corrected to ∥ρ(t,·)−˜ρ(t,·)∥.
  3. [Section 6.1, Lemma 6.3] The constant C_Ω is described only by 'depends only on the shape parameter c_usr of S_T_h, see [20, Section 5.6.2.2]'. Since the theorem claims fully computable bounds, an explicit expression or a precise derivation for C_Ω should be given.
  4. [Section 5, Remark 5.2] The claim about optimal orders of convergence of the Morley interpolation for d=2 is stated without proof or reference; if it is not needed for the main result, it could be removed or justified.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the conditional existence check is computed from the reconstruction and residual estimators, with an external blow-up criterion serving as the bridge to existence.

full rationale

The paper's derivation chain is genuinely conditional rather than circular. The residual bound (18) in Theorem 6.7 is assembled entirely from the numerical reconstruction (ρ̃,c̃) and computable quantities such as jumps, temporal differences, and the finite-element error estimator η_Ω; the exact solution (ρ,c) never enters the definition of the residual R_ρ̃. The stability frameworks in Theorems 3.3 and 3.7 assume a weak solution only as a working hypothesis and then show that if condition (7) or (9) holds, the error is bounded by δAE or Ψ_{N_t-1}. Proposition 2.5 then converts such a bound into an existence statement by contradiction with the blow-up criterion of Lemma 2.4: if the maximal existence time T_max were ≤ T, the verified error bound would force a uniform L^2 bound contradicting the asserted blow-up of ‖ρ‖_{L^2} at T_max. This is a standard continuation argument, not an instance of assuming the conclusion. The cited blow-up criterion [28, Lemma 2.5] comes from prior work by one of the present authors, but it is a parameter-free lemma with stated hypotheses that do not include the a posteriori condition, so it functions as independent external support rather than as a self-validating premise. The paper itself flags the non-explicit constant C_∞ in Remark 3.4 and the roughness of Sobolev constants for d=3 in Remark 3.8 and Section 7.1; these are practical verifiability and correctness concerns about whether condition (7)/(9) can actually be checked on feasible meshes, not circular reductions. No fitted parameter is later relabeled as a prediction, and no uniqueness theorem is imported from the authors' own work to force a particular choice. Accordingly, no load-bearing step reduces, by construction, to its own inputs.

Assumptions & free parameters 2 free parameters · 7 assumptions · 4 invented entities

The derivation relies on standard PDE theory (local weak existence, elliptic regularity, Sobolev embeddings), on the authors' prior blow-up criterion, and on new technical elements introduced here without full proof: the Morley-type element unisolvence and the maximum-norm a posteriori bound with explicit constant C∞. The practical bottleneck is the sharpness and computability of the constants, which the paper itself flags for d=3.

free parameters (2)
  • δ = 3/2 in Section 7 experiments
    User-selected parameter in conditions (7) and (9); the paper says δ=3/2 'leads to the best results in numerical experiments', so it is tuned rather than determined by the mathematics.
  • C∞
    Constant in Remark 3.4 for the maximum-norm a posteriori bound on ∇c̃; the paper asserts it can be estimated explicitly but does not provide a value, so the computable estimator depends on an unspecified quantity.
assumptions (7)
  • domain assumption Local existence of weak solutions for ρ0 ∈ L2 (Theorem 2.2, from [9])
    Proposition 2.5 and Corollaries 3.6/3.9 assume a maximal existence time Tmax > 0 for the exact solution; this is guaranteed by the cited existence theory for sufficiently smooth domains.
  • domain assumption Blow-up criterion: if Tmax < ∞ then ||ρ(t)||L∞ → ∞ as t ↗ Tmax (Lemma 2.4, from [28])
    The contradiction argument in Proposition 2.5 converts the uniform L2 bound into existence beyond T; this lemma is self-cited from [28].
  • domain assumption Elliptic regularity (1) with explicit constant C_ell for domains satisfying the H2 estimate
    Used throughout the stability estimates (equation 6) and in Lemmas 6.3/6.6; for convex polygons/polyhedra explicit C_ell exists by [6,29].
  • standard math Adjusted Generalized Gronwall Lemma (Proposition 3.1) holds as stated
    Adapted from [3, Proposition 6.2]; the proof is sketched in the paper and is standard.
  • ad hoc to paper Unisolvence and approximation properties of the Morley-type element (Lemma 5.1)
    The proof is omitted ('We omit it for brevity'); the element's flux-preservation property (13) is essential for the residual estimates in Lemmas 6.1-6.6.
  • ad hoc to paper Maximum-norm a posteriori bound ||∇c̃ - qh||L∞ ≤ C∞ η∞ with explicitly computable C∞ (Remark 3.4)
    Needed to turn the non-computable term ||∇c̃||L∞ in a(t) into a computable quantity; C∞ is not computed or bounded in the paper.
  • domain assumption Well-centeredness of the tetrahedral mesh and the associated dual mesh construction
    The FV scheme and the dual mesh for the FE concentration require circumcenters inside elements (Section 4); this restricts the class of admissible meshes.
invented entities (4)
  • Morley-type reconstruction ρ̃
    purpose: Maps piecewise constant FV solution into H1(Ω) preserving discrete diffusive fluxes (13), enabling residual estimation.
    New technical construct in Section 5; no falsifiable prediction outside the paper; validity rests on the unproved Lemma 5.1.
  • Dual mesh S_T_h
    purpose: Makes normal derivatives of P1-FE concentration well-defined on primal faces for the FV convective flux.
    Construct introduced in Section 4.2; its geometric validity depends on well-centeredness of the primal mesh.
  • Elliptic reconstruction c̃
    purpose: Reconstructs the concentration so that the perturbed system (3) has zero residual in the elliptic equation.
    Non-computable theoretical object (Remark 5.3); it is bounded via computable estimators.
  • Intersected mesh C_{T_h}^{S_T_h}
    purpose: Provides elements on which both div(ρ̃∇c_h) and Δρ̃ are well-defined for the convective residual estimates.
    Technical mesh construction in Section 6.1 used in Lemma 6.3.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A posteriori existence for the Keller-Segel model via a finite volume - finite element scheme." pith.science (2026). https://pith.science/paper/7BQAL5Q7

@misc{pith2026250917710,
  author       = {Pith},
  title        = {Pith review of: A posteriori existence for the Keller-Segel model via a finite volume - finite element scheme},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7BQAL5Q7}},
  note         = {Machine review of arXiv:2509.17710}
}
abstract

We derive two forms of conditional a posteriori error estimates for a finite volume scheme approximating the parabolic-elliptic Keller-Segel system. The estimates control the error in the $L^\infty(0,T, L^2(\Omega))$- and $L^2(0,T;H^1(\Omega))$-norm and exhibit linear convergence in the mesh size, as observed in numerical experiments. Crucially, we show that, as long as the condition of the error estimate is satisfied, a weak solution exists. This means, as long as the numerical solution has good properties, we can rigorously infer existence of an exact solution.

Figures

Figures reproduced from arXiv: 2509.17710 by the authors.

Figure 1
Figure 1. Illustration of the construction of a single dual element (red) with respect to two [PITH_FULL_IMAGE:figures/full_fig_p013_1.png] view at source ↗
Figure 2
Figure 2. Illustration of the primal mesh (blue), its circumcenters (green) and the correspond [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

56 extracted references · 40 canonical work pages

  1. [19]

    Maximum-norm a posteriori error estimates for singu- larly perturbed elliptic reaction-diffusion problems

    A. Demlow and N. Kopteva. “Maximum-norm a posteriori error estimates for singu- larly perturbed elliptic reaction-diffusion problems”. In:Numerische Mathematik133.4 (2016), pp. 707–742.doi:https://doi.org/10.1007/s00211-015-0763-0

    work page doi:10.1007/s00211-015-0763-0 2016

  2. [1]

    An unconditionally energy stable and positive upwind DG scheme for the Keller–Segel model

    D. Acosta-Soba, F. Guill´ en-Gonz´ alez, and J.R. Rodr´ ıguez-Galv´ an. “An unconditionally energy stable and positive upwind DG scheme for the Keller–Segel model”. In:Journal of Scientific Computing97.1 (2023), p. 18.doi:https://doi.org/10.1007/s10915- 023-02320-4

    doi:10.1007/s10915- 2023

  3. [2]

    Two and three dimensional H2-conforming finite element approximations without C1-elements

    M. Ainsworth and C. Parker. “Two and three dimensional H2-conforming finite element approximations without C1-elements”. In:Computer Methods in Applied Mechanics and Engineering431 (2024), p. 117267.doi:https://doi.org/10.1016/j.cma.2024. 117267

    doi:10.1016/j.cma.2024 2024

  4. [3]

    Bartels.Numerical methods for nonlinear partial differential equations

    S. Bartels.Numerical methods for nonlinear partial differential equations. Vol. 47. Springer Series in Computational Mathematics. Springer, 2015.doi:https : / / doi . org / 10 . 1007/978-3-319-13797-1

    work page 2015

  5. [4]

    Bartels.Numerical methods for partial differential equations

    S. Bartels.Numerical methods for partial differential equations. Vol. 64. Texts in Applied Mathematics. Springer, 2015.doi:https://doi.org/10.1007/978-3-319-32354-1

    work page doi:10.1007/978-3-319-32354-1 2015

  6. [5]

    Robust a priori and a posteriori error analysis for the approximation of Allen–Cahn and Ginzburg–Landau equations past topologi- cal changes

    S. Bartels, R. M¨ uller, and C. Ortner. “Robust a priori and a posteriori error analysis for the approximation of Allen–Cahn and Ginzburg–Landau equations past topologi- cal changes”. In:SIAM journal on numerical analysis49.1 (2011), pp. 110–134.url: https://doi.org/10.1137/090751530

    work page doi:10.1137/090751530 2011

  7. [6]

    A note on the Poincar´ e inequality for convex domains

    M. Bebendorf. “A note on the Poincar´ e inequality for convex domains”. In:Zeitschrift f¨ ur Analysis und ihre Anwendungen22.4 (2003), pp. 751–756.doi:https://doi.org/ 10.4171/zaa/1170

    doi:10.4171/zaa/1170 2003

  8. [7]

    A finite volume scheme for convection–diffusion equations with nonlinear diffusion derived from the Scharfetter–Gummel scheme

    M. Bessemoulin-Chatard. “A finite volume scheme for convection–diffusion equations with nonlinear diffusion derived from the Scharfetter–Gummel scheme”. In:Numerische Mathematik121 (2012), pp. 637–670.doi:https://doi.org/10.1007/s00211-012- 0448-x

    doi:10.1007/s00211-012- 2012

Show all 56 references
  1. [8]

    Local and global solvability of some parabolic systems modelling chemotaxis

    P. Biler. “Local and global solvability of some parabolic systems modelling chemotaxis”. In:Advances in Mathematical Sciences and Applications8 (1998), pp. 715–743.url: https://api.semanticscholar.org/CorpusID:115296668

    1998

  2. [9]

    Existence and nonexistence of solutions for a model of grav- itational interaction of particles, I

    P. Biler and T. Nadzieja. “Existence and nonexistence of solutions for a model of grav- itational interaction of particles, I”. In:Colloquium Mathematicae. Vol. 66. 2. 1993, pp. 319–334.doi:https://doi.org/10.4064/cm-66-2-319-334

    1993 doi

  3. [10]

    S. C. Brenner and L. R. Scott.The mathematical theory of finite element methods. Springer, 2008.doi:https://doi.org/10.1007/978-0-387-75934-0

    2008 doi

  4. [11]

    A Posteriori Error Control for the Allen–Cahn Equation with Variable Mobility

    A. Brunk, J. Giesselmann, and M. Luk´ aˇ cov´ a-Medviˇdov´ a. “A Posteriori Error Control for the Allen–Cahn Equation with Variable Mobility”. In:SIAM Journal on Numerical Analysis63.4 (2025), pp. 1540–1560.doi:https://doi.org/10.1137/24M1646406

    2025 doi

  5. [12]

    Brunk, J

    A. Brunk, J. Giesselmann, and T. Tscherpel.A Posteriori Existence of strong solutions to the Navier-Stokes equations in 3D. (in preparation)

  6. [13]

    A numerical verification method for a system of FitzHugh-Nagumo type

    S. Cai, K. Nagatou, and Y. Watanabe. “A numerical verification method for a system of FitzHugh-Nagumo type”. In:Numerical functional analysis and optimization33.10 (2012), pp. 1195–1220.doi:https://doi.org/10.1080/01630563.2012.677918

    2012

  7. [14]

    Adaptivity and blow-up detection for nonlinear evolution problems

    A. Cangiani, E. H. Georgoulis, I. Kyza, and S. Metcalfe. “Adaptivity and blow-up detection for nonlinear evolution problems”. In:SIAM Journal on Scientific Computing 38.6 (2016), A3833–A3856.doi:https://doi.org/10.1137/16M106073X. 30

    2016 doi

  8. [15]

    A posteriori error estimates for un- steady convection–diffusion–reaction problems and the finite volume method

    N. Chalhoub, A. Ern, T. Sayah, and M. Vohral´ ık. “A posteriori error estimates for un- steady convection–diffusion–reaction problems and the finite volume method”. In:Finite Volumes for Complex Applications VI Problems & Perspectives. Springer Proceedings in Mathematics. Spri...

    2011 doi

  9. [16]

    Error estimates and blow-up analysis of a finite-element approximation for the parabolic-elliptic Keller-Segel system

    W. Chen, Q. Liu, and J. Shen. “Error estimates and blow-up analysis of a finite-element approximation for the parabolic-elliptic Keller-Segel system”. In:International Journal of Numerical Analysis and Modeling19.2–3 (2022), pp. 275–298.url:https://doi. org/10.48550/arXiv.2212.07655

    2022 doi arXiv

  10. [17]

    A second-order positivity preserving central-upwind scheme for chemotaxis and haptotaxis models

    A. Chertock and A. Kurganov. “A second-order positivity preserving central-upwind scheme for chemotaxis and haptotaxis models”. In:Numerische Mathematik111 (2008), pp. 169–205.doi:https://doi.org/10.1007/s00211-008-0188-0

    2008 doi

  11. [18]

    Convergence rate of a finite volume schemefor a two dimensional convection-diffusion problem

    Y. Coudi` ere, J.-P. Vila, and P. Villedieu. “Convergence rate of a finite volume schemefor a two dimensional convection-diffusion problem”. In:ESAIM: Mathematical Modelling and Numerical Analysis33.3 (1999), pp. 493–516.doi:https://doi.org/10.1051/ m2an:1999149

    1999

  12. [20]

    D. A. Di Pietro and A. Ern.Mathematical aspects of discontinuous Galerkin methods. Vol. 69. Math´ ematiques et Applications. Springer, 2011.doi:https://doi.org/10. 1007/978-3-642-22980-0

    2011

  13. [21]

    Green’s matrices of second order elliptic systems with measurable coefficients in two dimensional domains

    H. Dong and S. Kim. “Green’s matrices of second order elliptic systems with measurable coefficients in two dimensional domains”. In:Transactions of the American Mathemat- ical Society361.6 (2009), pp. 3303–3323.doi:https://doi.org/10.1090/S0002- 9947-09-04805-3

    2009 doi

  14. [22]

    Fully discrete analysis of a discontinuous finite element method for the Keller-Segel chemotaxis model

    Y. Epshteyn and A. Izmirlioglu. “Fully discrete analysis of a discontinuous finite element method for the Keller-Segel chemotaxis model”. In:Journal of Scientific Computing40 (2009), pp. 211–256.doi:https://doi.org/10.1007/s10915-009-9281-5

    2009 doi

  15. [23]

    New interior penalty discontinuous Galerkin methods for the Keller–Segel chemotaxis model

    Y. Epshteyn and A. Kurganov. “New interior penalty discontinuous Galerkin methods for the Keller–Segel chemotaxis model”. In:SIAM Journal on Numerical Analysis47.1 (2009), pp. 386–408.doi:https://doi.org/10.1137/07070423X

    2009 doi

  16. [24]

    A posteriori error estimation based on potential and flux reconstruction for the heat equation

    A. Ern and M. Vohral´ ık. “A posteriori error estimation based on potential and flux reconstruction for the heat equation”. In:SIAM Journal on Numerical Analysis48.1 (2010), pp. 198–223.doi:https://doi.org/10.1137/090759008

    2010 doi

  17. [25]

    Estimates for the Sobolev trace constant with critical exponent and applications

    J. Fern´ andez Bonder and N. Saintier. “Estimates for the Sobolev trace constant with critical exponent and applications”. In:Annali di Matematica Pura ed Applicata187.4 (2008), pp. 683–704.doi:https://doi.org/10.1007/s10231-007-0062-1

    2008 doi

  18. [26]

    A finite volume scheme for the Patlak–Keller–Segel chemotaxis model

    F. Filbet. “A finite volume scheme for the Patlak–Keller–Segel chemotaxis model”. In: Numerische Mathematik104 (2006), pp. 457–488.doi:https://doi.org/10.1007/ s00211-006-0024-3

    2006

  19. [27]

    A posteriori error analysis of a positivity preserving scheme for the power-law diffusion Keller–Segel model

    J. Giesselmann and N. Kolbe. “A posteriori error analysis of a positivity preserving scheme for the power-law diffusion Keller–Segel model”. In:IMA Journal of Numerical Analysis(2024), drae073.doi:https://doi.org/10.1093/imanum/drae073. 31

    2024 doi

  20. [28]

    A posteriori error control for a discontinuous Galerkin approximation of a Keller-Segel model

    J. Giesselmann and K. Kwon. “A posteriori error control for a discontinuous Galerkin approximation of a Keller-Segel model”. In:Advances in Computational Mathematics 50.6 (2024), pp. 1–32.doi:https://doi.org/10.1007/s10444-024-10212-w

    2024 doi

  21. [29]

    Grisvard.Elliptic Problems in Nonsmooth Domains

    P. Grisvard.Elliptic Problems in Nonsmooth Domains. Classics in Applied Mathematics. Society for Industrial and Applied Mathematics, 2011.doi:https://doi.org/10. 1137/1.9781611972030

    2011

  22. [30]

    Energy dissipative local discontinuous Galerkin methods for Keller–Segel chemotaxis model

    L. Guo, X. H. Li, and Y. Yang. “Energy dissipative local discontinuous Galerkin methods for Keller–Segel chemotaxis model”. In:Journal of Scientific Computing78 (2019), pp. 1387–1404.doi:https://doi.org/10.1007/s10915-018-0813-8

    2019 doi

  23. [31]

    Exact sequences on Worsey–Farin splits

    J. Guzm´ an, A. Lischke, and M. Neilan. “Exact sequences on Worsey–Farin splits”. In: Mathematics of Computation91.338 (2022), pp. 2571–2608.doi:https://doi.org/ 10.1090/mcom/3746

    2022 doi

  24. [32]

    TetGen, a Delaunay-based quality tetrahedral mesh generator

    S. Hang. “TetGen, a Delaunay-based quality tetrahedral mesh generator”. In:ACM Trans. Math. Softw.41.2 (2015), p. 11.doi:https://doi.org/10.1145/2629697

    2015 doi

  25. [33]

    Efficient Numerical Schemes for a Two-Species Keller-Segel Model and Investigation of Its Blowup Phenomena in 3D

    X. Huang and J. Shen. “Efficient Numerical Schemes for a Two-Species Keller-Segel Model and Investigation of Its Blowup Phenomena in 3D”. In:Acta Applicandae Math- ematicae190.1 (2024), p. 10.doi:https://doi.org/10.1007/s10440-024-00647-0

    2024 doi

  26. [34]

    A discrete boundedness-by-entropy method for finite-volume approximations of cross-diffusion systems

    A. J¨ ungel and A. Zurek. “A discrete boundedness-by-entropy method for finite-volume approximations of cross-diffusion systems”. In:IMA Journal of Numerical Analysis43.1 (2023), pp. 560–589.doi:https://doi.org/10.1093/imanum/drab101

    2023 doi

  27. [35]

    Initiation of slime mold aggregation viewed as an in- stability

    E. F. Keller and L. A. Segel. “Initiation of slime mold aggregation viewed as an in- stability”. In:Journal of Theoretical Biology26.3 (1970), pp. 399–415.doi:https : //doi.org/10.1016/0022-5193(70)90092-5

    1970 doi

  28. [36]

    A general approach to transforming finite elements

    R. C. Kirby. “A general approach to transforming finite elements”. In:The SMAI journal of computational mathematics4 (2018), pp. 197–224.doi:https://doi.org/ 10.5802/smai-jcm.33

    2018 doi

  29. [37]

    Analysis for time discrete approximations of blow-up solu- tions of semilinear parabolic equations

    I. Kyza and C. Makridakis. “Analysis for time discrete approximations of blow-up solu- tions of semilinear parabolic equations”. In:SIAM journal on numerical analysis49.1 (2011), pp. 405–426.doi:https://doi.org/10.1137/100796819

    2011 doi

  30. [38]

    Pointwise a posteriori error bounds for blow-up in the semi- linear heat equation

    I. Kyza and S. Metcalfe. “Pointwise a posteriori error bounds for blow-up in the semi- linear heat equation”. In:SIAM Journal on Numerical Analysis58.5 (2020), pp. 2609– 2631.doi:https://doi.org/10.1137/19M1264758

    2020 doi

  31. [39]

    Local discontinuous Galerkin method for the Keller- Segel chemotaxis model

    X. H. Li, C.-W. Shu, and Y. Yang. “Local discontinuous Galerkin method for the Keller- Segel chemotaxis model”. In:Journal of Scientific Computing73.2 (2017), pp. 943–967. doi:https://doi.org/10.1007/s10915-016-0354-y

    2017 doi

  32. [40]

    Space and time reconstructions in a posteriori analysis of evolution problems

    C. Makridakis. “Space and time reconstructions in a posteriori analysis of evolution problems”. In:ESAIM: Proceedings. Vol. 21. EDP Sciences. 2007, pp. 31–44.doi: https://doi.org/10.1051/proc:072104

    2007 doi

  33. [41]

    Elliptic reconstruction and a posteriori error esti- mates for parabolic problems

    C. Makridakis and R. H. Nochetto. “Elliptic reconstruction and a posteriori error esti- mates for parabolic problems”. In:SIAM Journal on Numerical Analysis41.4 (2003), pp. 1585–1594.doi:https://doi.org/10.1137/S0036142902406314. 32

    2003 doi

  34. [42]

    Solutions of the 3D Navier–Stokes equations for initial data in ˙H1/2: robustness of regularity and numerical verification of regularity for bounded sets of initial data in ˙H1

    P. Mar´ ın-Rubio, J.C. Robinson, and W. Sadowski. “Solutions of the 3D Navier–Stokes equations for initial data in ˙H1/2: robustness of regularity and numerical verification of regularity for bounded sets of initial data in ˙H1”. In:Journal of Mathematical Analysis and Applica...

    2013 doi

  35. [43]

    A method of verified computa- tions for solutions to semilinear parabolic equations using semigroup theory

    M. Mizuguchi, A. Takayasu, T. Kubo, and S. Oishi. “A method of verified computa- tions for solutions to semilinear parabolic equations using semigroup theory”. In:SIAM journal on numerical analysis55.2 (2017), pp. 980–1001.doi:https://doi.org/10. 1137/141001664

    2017

  36. [44]

    Estimation of Sobolev embedding constant on a domain dividable into bounded convex domains

    M. Mizuguchi, K. Tanaka, K. Sekine, and S. Oishi. “Estimation of Sobolev embedding constant on a domain dividable into bounded convex domains”. In:Journal of Inequali- ties and Applications(2017), pp. 1–18.doi:https://doi.org/10.1186/s13660-017- 1571-0

    2017 doi

  37. [45]

    On approximate solutions of the incompressible Euler and Navier–Stokes equations

    C. Morosi and L. Pizzocchero. “On approximate solutions of the incompressible Euler and Navier–Stokes equations”. In:Nonlinear Analysis: Theory, Methods & Applications 75.4 (2012), pp. 2209–2235.doi:https://doi.org/10.1016/j.na.2011.10.020

    2012 doi

  38. [46]

    Behavior of radially symmetric solutions of a system related to chemotaxis

    T. Nagai and T. Senba. “Behavior of radially symmetric solutions of a system related to chemotaxis”. In:Nonlinear Analysis: Theory, Methods & Applications30.6 (1997), pp. 3837–3842.doi:https://doi.org/10.1016/S0362-546X(96)00256-8

    1997 doi

  39. [47]

    A Posteriori Error Estimations of Some Cell Centered Finite Volume Meth- ods for Diffusion-Convection-Reaction Problems

    S. Nicaise. “A Posteriori Error Estimations of Some Cell Centered Finite Volume Meth- ods for Diffusion-Convection-Reaction Problems”. In:SIAM Journal on Numerical Anal- ysis44.3 (2006), pp. 949–978.doi:https://doi.org/10.1137/040611483

    2006 doi

  40. [48]

    A posteriori existence in numerical computations

    C. Ortner. “A posteriori existence in numerical computations”. In:SIAM Journal on Numerical Analysis47.4 (2009), pp. 2550–2577.doi:https : / / doi . org / 10 . 1137 / 060668183

    2009

  41. [49]

    A mixed finite element method for 2-nd order elliptic problems

    P.-A. Raviart and J.-M. Thomas. “A mixed finite element method for 2-nd order elliptic problems”. In:Mathematical Aspects of Finite Element Methods: Proceedings of the Conference Held in Rome, December 10–12, 1975. Springer. 2006, pp. 292–315.doi: https://doi.org/10.1007/BFb0064470

    1975 doi

  42. [50]

    Conservative upwind finite-element method for a simplified Keller–Segel sys- tem modelling chemotaxis

    N. Saito. “Conservative upwind finite-element method for a simplified Keller–Segel sys- tem modelling chemotaxis”. In:IMA journal of numerical analysis27.2 (2007), pp. 332– 365.doi:https://doi.org/10.1093/imanum/drl018

    2007 doi

  43. [51]

    A positivity-preserving finite element method for chemotaxis problems in 3D

    R. Strehl, A. Sokolov, D. Kuzmin, D. Horstmann, and S. Turek. “A positivity-preserving finite element method for chemotaxis problems in 3D”. In:Journal of Computational and Applied Mathematics239 (2013), pp. 290–303.doi:https://doi.org/10.1016/ j.cam.2012.09.041

    2013

  44. [52]

    A positivity preserving moving mesh finite element method for the Keller–Segel chemotaxis model

    M. Sulman and T. Nguyen. “A positivity preserving moving mesh finite element method for the Keller–Segel chemotaxis model”. In:Journal of Scientific Computing80 (2019), pp. 649–666.doi:https://doi.org/10.1007/s10915-019-00951-0

    2019 doi

  45. [53]

    Triangulation of simple 3D shapes with well-centered tetrahedra

    E. VanderZee, A. N. Hirani, and D. Guoy. “Triangulation of simple 3D shapes with well-centered tetrahedra”. In:Proceedings of the 17th international meshing roundtable. Springer, 2008, pp. 19–35.doi:https://doi.org/10.1007/978-3-540-87921-3_2

    2008 doi

  46. [54]

    Verf¨ urth.A posteriori error estimation techniques for finite element methods

    R. Verf¨ urth.A posteriori error estimation techniques for finite element methods. Oxford University Press, 2013.doi:https://doi.org/10.1093/acprof:oso/9780199679423. 001.0001. 33

    2013

  47. [55]

    Error estimates for somequasi-interpolation operators

    R. Verf¨ urth. “Error estimates for somequasi-interpolation operators”. In:ESAIM: Math- ematical Modelling and Numerical Analysis33.4 (1999), pp. 695–713.doi:https:// doi.org/10.1051/m2an:1999158

    1999 doi

  48. [56]

    Modeling chemotactic waves in saturated porous media using adaptive mesh refinement

    N.D. Ward, S. Falle, and M. S. Olson. “Modeling chemotactic waves in saturated porous media using adaptive mesh refinement”. In:Transport in Porous Media89 (2011), pp. 487–504.doi:https://doi.org/10.1007/s11242-011-9782-1. 34

    2011 doi

Pith tools

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