pith. sign in

arxiv: 1907.06590 · v1 · pith:3KUCBRQ3new · submitted 2019-07-11 · 🧬 q-bio.TO

A tractable mathematical model for tissue growth

Joe Eyles , John F. King , Vanessa Styles This is my paper

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

classification 🧬 q-bio.TO
keywords tissue growthfree boundary problemmean curvature flowasymptotic analysistumor modelingkinetic under-coolingdiffuse interface approximation
0
0 comments X

The pith

Formal asymptotic methods derive a free boundary problem for tissue growth as forced mean curvature flow driven by an interior PDE.

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

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.

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.

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.

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.

Where Pith is reading between the lines

These are 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.

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.

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.

Figures

Figures reproduced from arXiv: 1907.06590 by Joe Eyles, John F. King, Vanessa Styles.

Figure 1
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 ζ(ϕ)… view at source ↗
Figure 2
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. Figure 3: Comparison of the radius (upper plots) and the pressure [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗
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 with α = 0.65 such that 3αβ < 2 (linearly unstable). The instability manifest in the right plot is associated with the higher perimeter–length… view at source ↗
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] view at source ↗
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 scheme (bottom row) [PITH_FULL_IMAGE:figures/full_fig_p020_6.png] view at source ↗
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–interface solutions are at t = 0, 4, 5 and the parametric solutions are at t = 0, 3.5, 4.4 [PITH_FULL_IMAGE:figures/full_fig_p022_7.png] view at source ↗
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] view at source ↗
Figure 9
Figure 9. Figure 9: Comparison of the radius (upper plots) and the pressure [PITH_FULL_IMAGE:figures/full_fig_p023_9.png] view at source ↗
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] view at source ↗
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–growth model. 6.5.5 R 2 : β = 0, α = 0.1, γ = 0.1 [PITH_FULL_IMAGE:figures/full_fig_p025_11.png] view at source ↗
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 row) [PITH_FULL_IMAGE:figures/full_fig_p026_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Simulation using the parametric scheme with [PITH_FULL_IMAGE:figures/full_fig_p026_13.png] view at source ↗
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 the oblate spheroid with equation x 2 1.0 2 + y 2 0.5 2 + z 2 1.0 2 = 1 . In the parametric examples the mesh size was taken to be h ≈ 0.15 a… view at source ↗
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_p028_15.png] view at source ↗
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/full_fig_p029_16.png] view at source ↗
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.

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

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.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 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.

pith-pipeline@v0.9.0 · 5616 in / 1205 out tokens · 47244 ms · 2026-05-24T22:58:25.150340+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

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

  9. [9]

    An algorithm for evolutionary surfaces

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

    work page 1991

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

    work page 2016

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

    work page 2019

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

    work page 2016

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

    work page 2009

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

    work page 1972

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

    work page 1976

  16. [16]

    Visualization Handbook

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

    work page 2005

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

    work page 2017

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

    work page 2004

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

    work page 2006

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

    work page 2007

  21. [21]

    Design of adaptive finite element software

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

    work page 2006

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

    work page 2003

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

    work page 1992

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

    work page 2006