pith. sign in

arxiv: 2501.13216 · v2 · pith:VJE5O4FPnew · submitted 2025-01-22 · 🧮 math.NA · cs.NA

On a linear DG approximation of chemotaxis models with damping gradient nonlinearities

Daniel Acosta-Soba , Alessandro Columbu , J. Rafael Rodr\'iguez-Galv\'an This is my paper

Pith reviewed 2026-05-23 05:22 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords chemotaxis modelsdiscontinuous Galerkinpositivity preservingdamping gradient nonlinearitynonlinear diffusionchemotactic collapsenumerical approximation
0
0 comments X

The pith

A linear upwind DG approximation of chemotaxis models with damping gradient nonlinearities preserves positivity and prevents blow-up.

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

This paper develops a linear upwind discontinuous Galerkin method for chemotaxis models that incorporate damping gradient nonlinearities. The models include nonlinear diffusion along with chemoattraction, chemorepulsion, and logistic growth, considered in both local and nonlocal variants. The method is positivity preserving, and experiments on chemotactic collapse confirm that the damping term avoids blow-up, matching the analysis. Readers care because it offers a practical numerical tool for simulating cell movement and aggregation in biology without numerical instabilities from blow-up.

Core claim

The central claim is that a novel linear positivity-preserving upwind discontinuous Galerkin approximation can be constructed for chemotaxis models with damping gradient nonlinearities. For both local and nonlocal models, the scheme maintains positivity, and the damping gradient term is shown to prevent blow-up through numerical experiments that align with prior analysis.

What carries the argument

The upwind discontinuous Galerkin discretization of the damping gradient nonlinearity, which provides the regularization and positivity preservation in the discrete setting.

If this is right

  • The scheme applies to local and nonlocal models while preserving positivity.
  • The damping gradient term prevents finite-time blow-up in simulations of chemotactic collapse.
  • Numerical results are consistent with the theoretical analysis of the approximation.
  • The method handles nonlinear diffusion, chemoattraction, chemorepulsion, and logistic growth.

Where Pith is reading between the lines

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

  • This linear approach could be adapted to other types of nonlinear PDEs with damping mechanisms to ensure stability.
  • It may facilitate long-term simulations of pattern formation in chemotaxis systems.
  • Similar positivity properties might be investigated in finite element or finite volume alternatives.
  • The framework could be tested on three-dimensional or more complex biological scenarios.

Load-bearing premise

The damping gradient nonlinearity serves as the regularizing mechanism in both the continuous model and the discrete scheme to control blow-up and ensure positivity.

What would settle it

A simulation of the scheme applied to a standard chemotactic collapse problem without the damping gradient term, which would be expected to exhibit either positivity violation or finite-time blow-up if the assumption holds.

Figures

Figures reproduced from arXiv: 2501.13216 by Alessandro Columbu, Daniel Acosta-Soba, J. Rafael Rodr\'iguez-Galv\'an.

Figure 1
Figure 1. Figure 1: Πh 1u m at different time steps in Test 5.1 (c = 0, χ = 5, ξ = 0) 11 [PITH_FULL_IMAGE:figures/full_fig_p011_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: v m at different time steps in Test 5.1 (c = 0, χ = 5, ξ = 0) To the best of our knowledge, there are no theoretical results providing conditions for blow-up in the model (2) with τ = 1. Indeed, the existing literature on blow-up phenomena for fully parabolic chemotaxis models is very limited. While progress has been made in the simplified parabolic-elliptic cases, the techniques used there do not easily e… view at source ↗
Figure 3
Figure 3. Figure 3: ∥u m∥L∞(Ω) over time for different values of c and γ in Test 5.1 (χ = 5, ξ = 0) 13 [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: ∥u m∥L∞(Ω) over time for different values of c and γ in Test 5.2 (χ = 5, ξ = 1) 5.3 Attraction-repulsion nonlocal model Now, we consider the nonlocal model (3) and compute a numerical experiment in the two-dimensional domain Ω = {(x, y) ∈ R 2 : x 2 + y 2 < 1}. We take the initial condition for the cell density u0(x, y, z) := 100e −35(x 2+y 2 ) , 14 [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Πh 1u m at different time steps in Test 5.3 (c = 0, α = 1.5) With this choice, we are under the assumptions in [11, Theorem 3], where for χ, ξ, λ, µ, n2 = n3 > 0, ρ = 1, k > 1 and α > β and ( n2 + α > max{n1 + d 2 k, k} if n1 ≥ 0, n2 + α > max{ d 2 k, k} if n1 < 0, that guarantee the existence of a blowing-up solution in the case c = 0 for a certain initial condition that satisfies R Ω u0 ≥  λ µ |Ω| k−1 … view at source ↗
Figure 6
Figure 6. Figure 6: Πh 1u m (first row), v m (second row) and w m (third row) at t = 3 · 10−3 for different values of γ and c = 10−3 in Test 5.3 (α = 1.5) 16 [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Πh 1u m (first row), v m (second row) and w m (third row) at t = 3 · 10−3 for different values of γ and c = 10−1 in Test 5.3 (α = 1.5) 17 [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
read the original abstract

In this work we present a novel linear and positivity preserving upwind discontinuous Galerkin (DG) approximation of a class of chemotaxis models with damping gradient nonlinearities. In particular, both a local and a nonlocal model including nonlinear diffusion, chemoattraction, chemorepulsion and logistic growth are considered. Some numerical experiments in the context of chemotactic collapse are presented, whose results are in accordance with the previous analysis of the approximation and show how the blow-up can be prevented by means of the damping gradient term.

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 / 3 minor

Summary. The manuscript presents a novel linear upwind discontinuous Galerkin (DG) approximation for a class of chemotaxis models with damping gradient nonlinearities. Both local and nonlocal variants are considered, incorporating nonlinear diffusion, chemoattraction, chemorepulsion, and logistic growth. The central claims are that the scheme is positivity-preserving (via upwind fluxes and the damping structure) and that numerical experiments on chemotactic collapse scenarios align with prior analysis while demonstrating that the damping gradient term prevents blow-up.

Significance. If the positivity preservation holds as claimed, the work supplies a practical, linear DG scheme for simulating chemotaxis models without artificial blow-up, which is valuable in mathematical biology. The explicit use of upwind fluxes combined with the damping nonlinearity to establish positivity is a clear strength, as is the reported consistency between the discrete experiments and the continuous-model analysis. These elements provide a falsifiable numerical check on the regularization role of the damping term.

minor comments (3)
  1. [Abstract] Abstract: the reference to 'previous analysis of the approximation' is ambiguous as to whether the analysis appears in the present manuscript or a cited prior work; a brief clarification or citation would improve precision.
  2. The numerical section would be strengthened by the inclusion of at least one table or set of quantitative metrics (e.g., maximum cell density over time with/without damping) to make the blow-up prevention claim more directly verifiable.
  3. Notation for the local versus nonlocal models could be introduced with a short comparative table early in the model section to aid readers.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading and positive assessment of the manuscript. The recommendation for minor revision is noted. No specific major comments appear in the report, so we provide no point-by-point responses below. We remain ready to address any minor points or clarifications requested by the editor.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained

full rationale

The manuscript introduces a new linear upwind DG scheme for the given chemotaxis models, derives positivity and stability directly from the discrete flux structure and the damping gradient nonlinearity in the continuous equations, and confirms consistency via numerical experiments. No load-bearing step reduces to a fitted parameter renamed as prediction, self-citation chain, or self-definitional equivalence; the damping term's regularizing effect is analyzed from the model PDEs rather than presupposed by the discretization itself.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The abstract invokes standard DG theory for hyperbolic/parabolic problems and the existence of a damping term that regularizes the continuous model; no free parameters, new entities, or ad-hoc axioms are named.

axioms (1)
  • domain assumption Standard assumptions of DG theory for convection-diffusion equations guarantee positivity and stability when upwinding and damping are present.
    Invoked implicitly when the authors state that the approximation is positivity preserving and that results accord with analysis.

pith-pipeline@v0.9.0 · 5617 in / 1191 out tokens · 62783 ms · 2026-05-23T05:22:16.626894+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

40 extracted references · 40 canonical work pages

  1. [1]

    Acosta-Soba, F

    D. Acosta-Soba, F. Guill´ en-Gonz´ alez, and J. R. Rodr´ ıguez-Galv´ an. An Unconditionally En- ergy Stable and Positive Upwind DG Scheme for the Keller–Segel Model. Journal of Scientific Computing, 97(1):18, 2023. doi: 10.1007/s10915-023-02320-4

    work page doi:10.1007/s10915-023-02320-4 2023

  2. [2]

    Acosta-Soba, F

    D. Acosta-Soba, F. Guill´ en-Gonz´ alez, and J. R. Rodr´ ıguez-Galv´ an. An upwind DG scheme preserving the maximum principle for the convective Cahn–Hilliard model. Numerical Algo- rithms, 92(3):1589–1619, 2022. doi: 10.1007/s11075-022-01355-2

    work page doi:10.1007/s11075-022-01355-2 2022

  3. [3]

    M. Aida, T. Tsujikawa, M. Efendiev, A. Yagi, and M. Mimura. Lower estimate of the attractor dimension for a chemotaxis growth system. Journal of the London Mathematical Society , 74(2):453–474, 2006. doi: 10.1112/S0024610706023015

    work page doi:10.1112/s0024610706023015 2006

  4. [4]

    Almeda, H

    R. Almeda, H. van Someren Gr´ eve, and T. Kiørboe. Behavior is a major determinant of predation risk in zooplankton. Ecosphere, 8(2):e01668, 2017. doi: 10.1002/ecs2.1668

    work page doi:10.1002/ecs2.1668 2017

  5. [5]

    M. S. Alnaes, A. Logg, K. B. Ølgaard, M. E. Rognes, and G. N. Wells. Unified Form Lan- guage: A domain-specific language for weak formulations of partial differential equations.ACM Transactions on Mathematical Software, 40(2):1–37, 2014. doi: 10.1145/2566630

    work page doi:10.1145/2566630 2014

  6. [6]

    Andreianov, M

    B. Andreianov, M. Bendahmane, and M. Saad. Finite volume methods for degenerate chemo- taxis model. Journal of Computational and Applied Mathematics , 235(14):4015–4031, 2011. doi: 10.1016/j.cam.2011.02.023

    work page doi:10.1016/j.cam.2011.02.023 2011

  7. [7]

    Badia, J

    S. Badia, J. Bonilla, and J. V. Guti´ errez-Santacreu. Bound-preserving finite element approxi- mations of the Keller–Segel equations. Mathematical Models and Methods in Applied Sciences, 33(03):609–642, 2023. doi: 10.1142/S0218202523500148

    work page doi:10.1142/s0218202523500148 2023

  8. [8]

    Bellomo, A

    N. Bellomo, A. Bellouquid, Y. Tao, and M. Winkler. Toward a mathematical theory of Keller– Segel models of pattern formation in biological tissues. Mathematical Models and Methods in Applied Sciences, 25(09):1663–1763, 2015. doi: 10.1142/S021820251550044X

    work page doi:10.1142/s021820251550044x 2015

  9. [9]

    Chertock and A

    A. Chertock and A. Kurganov. A second-order positivity preserving central-upwind scheme for chemotaxis and haptotaxis models. Numerische Mathematik , 111(2):169–205, 2008. doi: 10.1007/s00211-008-0188-0

    work page doi:10.1007/s00211-008-0188-0 2008

  10. [10]

    P. G. Ciarlet. The Finite Element Method for Elliptic Problems . Classics in Applied Mathe- matics. SIAM, 2002. doi: 10.1137/1.9780898719208

    work page doi:10.1137/1.9780898719208 2002

  11. [11]

    Columbu, S

    A. Columbu, S. Frassu, and G. Viglialoro. Properties of given and detected unbounded so- lutions to a class of chemotaxis models. Studies in Applied Mathematics , 151(4):1349–1379,

    work page

  12. [12]

    doi: 10.1111/sapm.12627. 18

    work page doi:10.1111/sapm.12627

  13. [13]

    Columbu, R

    A. Columbu, R. D. Fuentes, and S. Frassu. Uniform-in-time boundedness in a class of local and nonlocal nonlinear attraction–repulsion chemotaxis models with logistics. Nonlinear Analysis: Real World Applications, 79:104135, 2024. doi: 10.1016/j.nonrwa.2024.104135

    work page doi:10.1016/j.nonrwa.2024.104135 2024

  14. [14]

    D. A. Di Pietro and A. Ern. Mathematical Aspects of Discontinuous Galerkin Methods , num- ber 69 in Math´ ematiques et Applications. Springer, 2012.doi: 10.1007/978-3-642-22980-0

    work page doi:10.1007/978-3-642-22980-0 2012

  15. [15]

    Ern and J.-L

    A. Ern and J.-L. Guermond. Theory and Practice of Finite Elements , number 159 in Applied mathematical sciences. Springer, 2004. doi: 10.1007/978-1-4757-4355-5

    work page doi:10.1007/978-1-4757-4355-5 2004

  16. [16]

    L. C. Evans. Partial differential equations , number 19 in Graduate Studies in Mathematics. American Mathematical Society, 2010. isbn: 978-1-4704-6942-9. doi: 10.1090/gsm/019

    work page doi:10.1090/gsm/019 2010

  17. [17]

    Geipel, C

    I. Geipel, C. E. Kernan, A. S. Litterer, G. G. Carter, R. A. Page, and H. M. Ter Hofstede. Predation risks of signalling and searching: bats prefer moving katydids. Biology Letters , 16(4):20190837, 2020. doi: 10.1098/rsbl.2019.0837

    work page doi:10.1098/rsbl.2019.0837 2020

  18. [18]

    Guill´ en-Gonz´ alez and G

    F. Guill´ en-Gonz´ alez and G. Tierra. Finite element numerical schemes for a chemo-attraction and consumption model. Journal of Computational and Applied Mathematics , 441:115676,

    work page

  19. [19]

    doi: 10.1016/j.cam.2023.115676

    work page doi:10.1016/j.cam.2023.115676 2023

  20. [20]

    Huang and J

    F. Huang and J. Shen. Bound/Positivity Preserving and Energy Stable Scalar auxiliary Variable Schemes for Dissipative Systems: Applications to Keller–Segel and Poisson–Nernst– Planck Equations. SIAM Journal on Scientific Computing , 43(3):A1832–A1857, 2021. doi: 10.1137/20M1365417

    work page doi:10.1137/20m1365417 2021

  21. [21]

    J. D. Hunter. Matplotlib: A 2D graphics environment. Computing in Science & Engineering , 9(3):90–95, 2007. doi: 10.1109/MCSE.2007.55

    work page doi:10.1109/mcse.2007.55 2007

  22. [22]

    Ishida, J

    S. Ishida, J. Lankeit, and G. Viglialoro. A Keller-Segel type taxis model with ecological inter- pretation and boundedness due to gradient nonlinearities.Discrete and Continuous Dynamical Systems-B, 29(9):3955–3969, 2024. doi: 10.3934/dcdsb.2024029

    work page doi:10.3934/dcdsb.2024029 2024

  23. [23]

    Preytaxis

    P. Kareiva and G. Odell. Swarms of Predators Exhibit “Preytaxis” if Individual Predators Use Area-Restricted Search. The American Naturalist, 130(2):233–270, 1987. doi: 10.1086/ 284707

    work page 1987

  24. [24]

    E. F. Keller and L. A. Segel. Initiation of slime mold aggregation viewed as an instability. Journal of Theoretical Biology, 26(3):399–415, 1970. doi: 10.1016/0022-5193(70)90092-5

    work page doi:10.1016/0022-5193(70)90092-5 1970

  25. [25]

    E. F. Keller and L. A. Segel. Model for chemotaxis. Journal of Theoretical Biology, 30(2):225– 234, 1971. doi: 10.1016/0022-5193(71)90050-6

    work page doi:10.1016/0022-5193(71)90050-6 1971

  26. [26]

    T. Li, D. Acosta-Soba, A. Columbu, and G. Viglialoro. Dissipative gradient nonlinearities prevent δ-formations in local and nonlocal attraction-repulsion chemotaxis models. arXiv preprint, 2024. doi: 10.48550/arXiv.2405.03586

    work page doi:10.48550/arxiv.2405.03586 2024

  27. [27]

    T. Li, S. Frassu, and G. Viglialoro. Combining effects ensuring boundedness in an attraction– repulsion chemotaxis model with production and consumption. Zeitschrift f¨ ur angewandte Mathematik und Physik , 74(3):109, 2023. doi: 10.1007/s00033-023-01976-0

    work page doi:10.1007/s00033-023-01976-0 2023

  28. [28]

    Mimura and T

    M. Mimura and T. Tsujikawa. Aggregating pattern dynamics in a chemotaxis model including growth. Physica A: Statistical Mechanics and its Applications , 230(3-4):499–543, 1996. doi: 10.1016/0378-4371(96)00051-9

    work page doi:10.1016/0378-4371(96)00051-9 1996

  29. [29]

    Quarteroni and A

    A. Quarteroni and A. Valli. Numerical approximation of partial differential equations , vol- ume 23. Springer Science & Business Media, 2008. doi: 10.1007/978-3-540-85268-1

    work page doi:10.1007/978-3-540-85268-1 2008

  30. [30]

    N. Saito. Conservative upwind finite-element method for a simplified Keller–Segel system modelling chemotaxis. IMA Journal of Numerical Analysis , 27(2):332–365, 2007. doi: 10 . 1093/imanum/drl018. 19

    work page 2007

  31. [31]

    M. W. Scroggs, I. A. Baratta, C. N. Richardson, and G. N. Wells. Basix: a runtime finite element basis evaluation library. Journal of Open Source Software , 7(73):3982, 2022. doi: 10.21105/joss.03982

    work page doi:10.21105/joss.03982 2022

  32. [32]

    M. W. Scroggs, J. S. Dokken, C. N. Richardson, and G. N. Wells. Construction of Arbitrary Order Finite Element Degree-of-Freedom Maps on Polygonal and Polyhedral Cell Meshes. ACM Transactions on Mathematical Software , 48(2):18, 2022. doi: 10.1145/3524456

    work page doi:10.1145/3524456 2022

  33. [33]

    Shen and J

    J. Shen and J. Xu. Unconditionally Bound Preserving and Energy Dissipative Schemes for a Class of Keller–Segel Equations. SIAM Journal on Numerical Analysis , 58(3):1674–1695,

    work page

  34. [34]

    doi: 10.1137/19M1246705

    work page doi:10.1137/19m1246705

  35. [35]

    P. Souplet. Finite Time Blow-up for a Non-linear Parabolic Equation with a Gradient Term and Applications. Mathematical methods in the applied sciences, 19(16):1317–1333, 1996. doi: 10.1002/(SICI)1099-1476(19961110)19:16<1317::AID-MMA835>3.0.CO;2-M

    work page doi:10.1002/(sici)1099-1476(19961110)19:16 1996

  36. [36]

    C. B. Sullivan and A. Kaszynski. PyVista: 3D plotting and mesh analysis through a stream- lined interface for the Visualization Toolkit (VTK).Journal of Open Source Software, 4(37):1450,

    work page

  37. [37]

    doi: 10.21105/joss.01450

    work page doi:10.21105/joss.01450

  38. [38]

    A. W. Visser. Motility of zooplankton: fitness, foraging and predation. Journal of Plankton research, 29(5):447–461, 2007. doi: 10.1093/plankt/fbm029

    work page doi:10.1093/plankt/fbm029 2007

  39. [39]

    Zhou and N

    G. Zhou and N. Saito. Finite volume methods for a Keller–Segel system: discrete energy, error estimates and numerical blow-up analysis. Numerische Mathematik, 135(1):265–311, January

    work page

  40. [40]

    doi: 10.1007/s00211-016-0793-2. 20

    work page doi:10.1007/s00211-016-0793-2