Pith. sign in

REVIEW 3 minor 24 references

A tractable mathematical model for tissue growth

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

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

Pith's one-line read Formal asymptotic methods derive a free boundary problem for tissue growth as forced mean curvature flow driven by an interior PDE.

desk verdict The paper derives a closed free-boundary model for tissue growth via asymptotics, then checks it with stability analysis and finite-element runs on a diffuse version. read the letter →

arxiv 1907.06590 v1 pith:3KUCBRQ3 submitted 2019-07-11 q-bio.TO

Joe Eyles , John F. King , Vanessa Styles This is my paper · ORCID
classification q-bio.TO
keywords tissuegrowthfreeboundaryproblemmeancurvatureflowasymptoticanalysistumormodelingkineticunder-coolingdiffuseinterfaceapproximation
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 uses formal asymptotic methods to derive a free boundary problem for the growth and death of tumors or biological tissue. This model describes a closed interface that evolves by forced mean curvature flow, regularized by kinetic under-cooling, with the forcing given by the solution of a PDE in the enclosed domain. Linear stability analysis is performed on the model, and a diffuse-interface approximation is derived. Finite-element methods discretize two related models, yielding computational results for comparison. Such a reduction provides a simple mathematical framework for studying tissue dynamics.

What carries the argument

Forced mean curvature flow of a closed interface with kinetic under-cooling regularization, forced by the solution of an interior PDE.

What would settle it

A direct numerical comparison of interface evolution and stability thresholds between this reduced model and simulations of more detailed underlying equations.

Watch

Extended reading notes

Core claim

Using formal asymptotic methods we derive a free boundary problem representing one of the simplest mathematical descriptions of the growth and death of a tumour or other biological tissue. The mathematical model takes the form of a closed interface evolving via forced mean curvature flow (together with a `kinetic under-cooling' regularisation) where the forcing depends on the solution of a PDE that holds in the domain enclosed by the interface. We perform linear stability analysis and derive a diffuse-interface approximation of the model. Finite-element discretisations of two closely related models are presented, together with computational results comparing the approximate solutions.

Load-bearing premise

Formal asymptotic reduction of unspecified underlying biological or mechanical equations yields a faithful free-boundary description whose forcing term is given by the solution of a single interior PDE.

Editorial extensions

If this is right

  • Linear stability analysis can be carried out on the interface evolution.
  • A diffuse-interface approximation can be derived from the free boundary model.
  • Finite-element discretizations enable numerical computation of solutions.
  • Computational results allow comparison between approximate solutions of related models.

Reading between the lines

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

  • The interior PDE could be chosen to represent different biological processes such as nutrient diffusion.
  • Numerical experiments with the model might reveal conditions for stable tumor shapes.
  • The approach could be applied to other free-boundary problems in biology by varying the asymptotic assumptions.
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 / 3 minor

Summary. The paper uses formal asymptotic methods to derive a free-boundary problem for the growth and death of biological tissue (e.g., tumours). The resulting model is a closed interface evolving by forced mean-curvature flow with kinetic under-cooling regularisation, where the forcing term is supplied by the solution of a PDE posed inside the domain. The authors then carry out linear stability analysis, derive a diffuse-interface approximation, and present finite-element discretisations of two related models together with computational results.

Significance. If the asymptotic reduction is correct, the work supplies one of the simplest closed mathematical descriptions of tissue growth that still retains an interior PDE and a free boundary. The combination of a formal derivation, explicit linear stability calculation, diffuse-interface reformulation, and finite-element computations provides a self-contained framework that could serve as a baseline for further analysis in mathematical biology. The manuscript supplies the starting equations, the asymptotic steps, the stability dispersion relation, and the numerical scheme, which are positive features.

minor comments (3)
  1. [Abstract] The abstract states that the model is 'one of the simplest mathematical descriptions'; this phrasing is subjective and could be replaced by a more neutral statement such as 'a simple mathematical description'.
  2. [Numerical results] In the numerical section, the captions of the computational figures should explicitly list the values of all non-dimensional parameters used in each run so that the results can be reproduced without consulting the main text.
  3. [Diffuse-interface approximation] The diffuse-interface approximation is introduced after the sharp-interface model; a brief remark on the expected convergence rate as the interface thickness tends to zero would help readers assess the approximation quality.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of the manuscript, including recognition of the asymptotic derivation, explicit stability calculation, diffuse-interface reformulation, and numerical results as a self-contained framework. The recommendation for minor revision is noted. As the report lists no specific major comments, we have no points requiring detailed rebuttal or revision at this stage.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity; derivation self-contained

full rationale

The paper performs an explicit formal asymptotic reduction from supplied starting equations for tissue growth to a forced mean-curvature free-boundary problem, followed by stability analysis, diffuse-interface approximation, and finite-element numerics. All load-bearing steps (asymptotic matching, dispersion relation, discretization) are derived internally and presented in full; no fitted parameters are relabeled as predictions, no self-citations carry the central claim, and the reduction is not equivalent to its inputs by definition. The modeling premise that the asymptotics are faithful is an explicit choice, not a hidden circularity.

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

Only the abstract is available, so the ledger records only the modelling choices explicitly named; no numerical values or new entities are stated.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A tractable mathematical model for tissue growth." pith.science (2026). https://pith.science/paper/3KUCBRQ3

@misc{pith2026190706590,
  author       = {Pith},
  title        = {Pith review of: A tractable mathematical model for tissue growth},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3KUCBRQ3}},
  note         = {Machine review of arXiv:1907.06590}
}
read the original abstract

Using formal asymptotic methods we derive a free boundary problem representing one of the simplest mathematical descriptions of the growth and death of a tumour or other biological tissue. The mathematical model takes the form of a closed interface evolving via forced mean curvature flow (together with a `kinetic under-cooling' regularisation) where the forcing depends on the solution of a PDE that holds in the domain enclosed by the interface. We perform linear stability analysis and derive a diffuse-interface approximation of the model. Finite-element discretisations of two closely related models are presented, together with computational results comparing the approximate solutions.

Figures

Figures reproduced from arXiv: 1907.06590 by the authors.

Figure 1
Figure 1. Diffuse–interface configuration. To derive a diffuse–interface approximation of (1.1), (1.2) we first introduce the variational form: Z Ω ∇u · ∇η dx + 1 α Z Γ uη dx = Z Γ (Q − βκ)η dx − Z Ω η dx ∀η ∈ H 1 (Ω). (4.1) We now follow the techniques described in [1] to obtain a diffuse–interface approximation of (4.1). Setting u˜ ∈ D to be a diffuse–interface approximation to u ∈ Ω we have Z Ω u(x) dx ≈ Z D ζ(ϕ)˜u(x) dx a… view at source ↗
Figure 2
Figure 2. An enlarged section of the diffuse–interface mesh (left plot) and parametric mesh (right [PITH_FULL_IMAGE:figures/full_fig_p016_2.png] view at source ↗
Figure 3
Figure 3. Comparison of the radius (upper plots) and the pressure [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗
Figures from the paper (13 more)
Figure 4
Figure 4. Figure 4: Linear stability, (3.5), with Q = 1, γ = 0 and β = 0.1. The left plot displays the initial geometry R(θ, 0) = 2Q + 0.25 sin(θ) cos(9θ), the centre plot displays uh at t = 700 with α = 0.68 such that 3αβ > 2 (linearly stable) and the right plot displays uh at t = 700 wi…
Figure 5
Figure 5. Figure 5: Results with γ = 0, α = 0.1 and β = 1.0: uh given by the parametric scheme at t = 0, 7, 10, 12.6. This simulation exhibits repeated tip splitting akin to that arising in a variety of moving boundary problems [PITH_FULL_IMAGE:figures/full_fig_p019_5.png]
Figure 6
Figure 6. Figure 6: Results with γ = 0, α = 1.0 and β = 0.1 at t = 0, 11, 23 (columns 1 to 3): uh given by the parametric scheme (top row), u˜h given by the diffuse–interface scheme (middle row), ϕh (in red and blue) from the diffuse–interface scheme and Xh (in white) from the parametric …
Figure 7
Figure 7. Figure 7: Results with γ = 0, α = 0.1 and β = 0.1: uh given by the parametric scheme (top row), u˜h given by the diffuse–interface scheme (middle row), ϕh (in red and blue) from the diffuse–interface scheme and Xh (in white) from the parametric scheme (bottom row). The diffuse–i…
Figure 8
Figure 8. Figure 8: Simulation relating to the thin–film limit using the parametric scheme (5.2), (5.1), with [PITH_FULL_IMAGE:figures/full_fig_p022_8.png]
Figure 9
Figure 9. Figure 9: Comparison of the radius (upper plots) and the pressure [PITH_FULL_IMAGE:figures/full_fig_p023_9.png]
Figure 10
Figure 10. Figure 10: presents results with β = 0, α = 1 and γ = 0.1, displayed in the same format as the results in [PITH_FULL_IMAGE:figures/full_fig_p024_10.png]
Figure 11
Figure 11. Figure 11: Simulations for β = 0, α = 1.0 and γ = 0.1 at t = 0, 1, 20, 30 using the diffuse–interface scheme, with u˜h displayed in the top row and ϕh in the bottom row. Such simulations illustrate the complicated pathologies that can result even for the current simple tissue–gr…
Figure 12
Figure 12. Figure 12: Simulations with β = 0, α = 0.1 and γ = 0.1 at t = 0, 3, 7: uh given by the parametric scheme (top row), u˜h given by the diffuse–interface scheme (middle row), ϕh (in red and blue) from the diffuse–interface scheme and Xh (in white) from the parametric scheme (bottom…
Figure 13
Figure 13. Figure 13: Simulation using the parametric scheme with [PITH_FULL_IMAGE:figures/full_fig_p026_13.png]
Figure 14
Figure 14. Figure 14: Simulations for β = 0 and α = 1.0, with γ = 0.05 at t = 20 (left plot), γ = 0.025 at t = 14 (centre plot) and γ = 0.01 at t = 11.6 (right plot). 6.5.7 R 3 : β = 0, α = 1.0, γ = 0.1 We set β = 0, α = 1.0, γ = 0.1 and Q = 1.25 and the initial geometry Γ(0) is given by t…
Figure 15
Figure 15. Figure 15: Parametric scheme with β = 0, α = 1.0, γ = 0.1: looking down the x axis (first row), looking down the y axis (second row), cross section in the z = 0 plane (third row), and uh on the plane z = 0 (fourth row). Taken at t = 0, 7, 13 [PITH_FULL_IMAGE:figures/full_fig_p0…
Figure 16
Figure 16. Figure 16: Diffuse–interface scheme with β = 0, α = 1.0, γ = 0.1: looking down the x axis (first row), looking down the y axis (second row), cross section in the z = 0 plane (third row), and u˜h on the plane z = 0 (fourth row). Taken at t = 0, 7, 13 [PITH_FULL_IMAGE:figures/ful…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

  1. [1]

    L AM, K.F., & S TINNER , B

    A BELS , H. L AM, K.F., & S TINNER , B. Analysis of the diffuse domain approach for a bulk– surface coupled PDE system. SIAM J. Numer. Anal., 47(5):3687–3725, 2015

    work page 2015

  2. [2]

    B ARRETT , J.W., N ¨URNBERG , R., & S TYLES , V Finite element approximation of a phase– field model for void electromigration. SIAM J. Numer. Anal., 42(2):738–772, 2005

    work page 2005

  3. [3]

    A phase–field model with a double obstacle potential

    B LOWEY J.F., & E LLIOTT , C.M. A phase–field model with a double obstacle potential. 1994

    work page 1994

  4. [4]

    Analysis of the diffuse domain method for second order elliptic boundary value problems

    B URGER , M., S CHLOTTBOM , M., & E LVETUN , O.L. Analysis of the diffuse domain method for second order elliptic boundary value problems. Found. Comput. Math., 17(3):627–674, 2017

    work page 2017

  5. [5]

    Free boundary value problems associated with the growth and development of multicellular spheroids

    C HAPLAIN , M.A.J., & B YRNE , H.M. Free boundary value problems associated with the growth and development of multicellular spheroids. Eur. J. Applied Math. , 8(6):639–658, 1997

    work page 1997

  6. [6]

    Nonlinear simulations of solid tumor growth using a mixture model: invasion and branching

    C RISTINI , V., L I, X., L OWENGRUB , J.S., & W ISE , S.M. Nonlinear simulations of solid tumor growth using a mixture model: invasion and branching. J. Math. Biol., 58(4–5):723– 763, 2009

    work page 2009

  7. [7]

    Computation of geometric partial differ- ential equations and mean curvature flow

    D ECKELNICK , K., E LLIOTT , C.M., & D ZIUK , G. Computation of geometric partial differ- ential equations and mean curvature flow. Acta Numerica, 14:139–232, 2005

    work page 2005

  8. [8]

    An existence and uniqueness result for a phase-field model of diffusion-induced grain-boundary motion

    D ECKELNICK , K., & E LLIOTT , C.M. An existence and uniqueness result for a phase-field model of diffusion-induced grain-boundary motion. Proceedings of the Royal Society of Ed- inburgh, pages 1323–1344, 2001

    work page 2001

Show all 24 references
  1. [9]

    An algorithm for evolutionary surfaces

    D ZIUK , G. An algorithm for evolutionary surfaces. Numerisch Mathematik, 58:603–611, 1991

    1991

  2. [10]

    On algorithms with good mesh properties for problems with moving boundaries based on the Harmonic Map Heat Flow and the DeTurck trick

    E LLIOTT , C.M., & F RITZ , H. On algorithms with good mesh properties for problems with moving boundaries based on the Harmonic Map Heat Flow and the DeTurck trick. SMAI J. Comput. Math., 2:141–176, 2016

    2016

  3. [11]

    Numerical analysis and simulations of a tractable model for tumour growth

    E YLES , J. Numerical analysis and simulations of a tractable model for tumour growth. PhD thesis, University of Sussex, 2019

    2019

  4. [12]

    A Cahn–Hilliard–Darcy model for tumour growth with chemotaxis and active transport

    G ARCKE , H., L AM, K.F, S ITKA , E., & S TYLES , V. A Cahn–Hilliard–Darcy model for tumour growth with chemotaxis and active transport. Math. Models Meth. Appl. Sci. , 26(6), 2016

    2016

  5. [13]

    Gmsh: a three–dimensional finite element mesh generator with built–in pre– and post–processing facilities.Int

    G EUZAINE , C., & R EMACLE , J.F. Gmsh: a three–dimensional finite element mesh generator with built–in pre– and post–processing facilities.Int. J. Numer. Meth. Eng., 79(11):1309–1331, 2009. A tractable mathematical model for tissue growth 35

    2009

  6. [14]

    Models for the growth of a solid tumour by diffusion

    G REENSPAN , H.P. Models for the growth of a solid tumour by diffusion. Stud. Appl. Math., 52:317–340, 1972

    1972

  7. [15]

    On the growth and stability of cell cultures and solid tumours

    G REENSPAN , H.P. On the growth and stability of cell cultures and solid tumours. J. Theor. Biol., 56:229–242, 1976

    1976

  8. [16]

    Visualization Handbook

    H ANSEN , C.D., & J OHNSON , C.R. Visualization Handbook. Elsevier–Butterworth Heine- mann, 2005

    2005

  9. [17]

    Interface dynamics in biological tissue growth

    K ING , J.R. Interface dynamics in biological tissue growth. C.M. E LLIOTT , Y. G IGA , M. HINZE AND V. S TYLES (EDS.), (2017), Emerging Developments in Interfaces and Free Boundaries, Oberwolfach Reports, 14, 267–338

    2017

  10. [18]

    Mathematical analysis of some multi–dimensional tissue growth models

    K ING , J.R., & F RANKS , S.J. Mathematical analysis of some multi–dimensional tissue growth models. Eur. J. Appl. Math., 15, 273–295, 2004

    2004

  11. [19]

    Mathematical modelling of nutrient–limited tissue growth

    K ING , J.R., & F RANKS , S.J. Mathematical modelling of nutrient–limited tissue growth. Free Boundary Problems. International Series of Numerical Mathematics, 154, 273–282, 2006

    2006

  12. [20]

    Stability properties of some tissue growth models

    K ING , J.R., & F RANKS , S.J. Stability properties of some tissue growth models. 175–182, in Mathematical Modelling of Biological Systems, V olume I, Eds. A Deutsch, L Brusch, H Byrne, G Vries, H Herzel, Springer 2007

    2007

  13. [21]

    Design of adaptive finite element software

    S CHMIDT , A., & S IEBERT , K.G. Design of adaptive finite element software. Springer Berlin Heidelberg, 2006

    2006

  14. [22]

    Mushy regions in negative squeeze films

    O CKENDON , J.R., H OWISON , S.D., & L ACEY, A.A. Mushy regions in negative squeeze films. Quart. J. Mech. Appl. Math., 56(3), 361–379, 2003

    2003

  15. [23]

    Complex variable methods in Hele–Shaw moving boundary problems

    H OWISON , S.D. Complex variable methods in Hele–Shaw moving boundary problems.. Eur. J. Appl. Math., 3, 209–224, 1992

    1992

  16. [24]

    Modelling the response of spatially structured tumours to chemotherapy: Drug kinetics

    N ORRIS , E.S., K ING , J.R., & B YRNE , H.M. Modelling the response of spatially structured tumours to chemotherapy: Drug kinetics. Math. Comp. Model., 43, 820–837, 2006

    2006

Pith tools

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