Multiphase modeling and qualitative analysis of the growth of tumor cords

In this paper a macroscopic model of tumor cord growth is developed, relying on the mathematical theory of deformable porous media. Tumor is modeled as a saturated mixture of proliferating cells, extracellular fluid and extracellular matrix, that occupies a spatial region close to a blood vessel whence cells get the nutrient needed for their vital functions. Growth of tumor cells takes place within a healthy host tissue, which is in turn modeled as a saturated mixture of non-proliferating cells. Interactions between these two regions are accounted for as an essential mechanism for the growth of the tumor mass. By weakening the role of the extracellular matrix, which is regarded as a rigid non-remodeling scaffold, a system of two partial differential equations is derived, describing the evolution of the cell volume ratio coupled to the dynamics of the nutrient, whose higher and lower concentration levels determine proliferation or death of tumor cells, respectively. Numerical simulations of a reference two-dimensional problem are shown and commented, and a qualitative mathematical analysis of some of its key issues is proposed.

💡 Research Summary

This paper presents a macroscopic, continuum‑scale model for the growth of tumor cords—thin, vascular‑adjacent sheets of proliferating cancer cells—by employing the theory of deformable porous media. The authors treat the tumor region as a saturated mixture of three constituents: proliferating cells, extracellular fluid, and extracellular matrix (ECM). The surrounding healthy tissue is modeled as a saturated mixture of non‑proliferating cells and fluid. Mass and momentum balance equations are written for each constituent, and interaction terms (drag between cells and fluid, pressure transmission through the ECM, etc.) are introduced to capture the mechanical coupling between the two regions.

A key simplifying assumption is that the ECM behaves as a rigid, non‑remodeling scaffold. By “weakening” the role of the ECM (i.e., neglecting its deformation and remodeling), the full multiphase system collapses to two coupled nonlinear partial differential equations (PDEs) governing (i) the cell volume fraction φ(x,t) and (ii) the nutrient concentration c(x,t). The nutrient field represents oxygen, glucose, or any essential substrate supplied by a nearby blood vessel. The governing equations can be written in a compact form:

∂tφ = ∇·(Dφ∇φ) + G(c) φ,
∂tc = ∇·(Dc∇c) – λ φ c,

where Dφ and Dc are effective diffusion coefficients for cells and nutrient, respectively; λ is the rate at which cells consume the nutrient; and G(c) is a growth‑death function that is positive when c exceeds a critical threshold c* (promoting proliferation) and negative when c falls below c* (inducing apoptosis/necrosis).

Boundary conditions reflect the physiological setting: the vessel wall (one side of the domain) imposes a fixed nutrient concentration c = c0 and zero cell volume fraction (φ = 0) because tumor cells have not yet invaded that surface; all other boundaries are taken as no‑flux (Neumann) to represent a closed tissue block. Nondimensionalisation introduces a small set of dimensionless groups—nutrient supply ratio, growth‑inhibition threshold, and cell‑nutrient coupling strength—that control the qualitative behaviour of the system.

The authors perform a qualitative mathematical analysis. They prove the existence of steady‑state solutions and examine their linear stability. When the nutrient supply is sufficiently high, a stable steady state exists in which φ decays exponentially away from the vessel, reproducing the observed thin proliferative layer. If the supply drops below a critical value, the steady state loses stability, leading to the formation of a necrotic core where φ collapses to zero. Bifurcation analysis shows how varying the dimensionless parameters can trigger transitions between sustained cord growth, limited growth with a central dead zone, and complete tumor regression.

Numerical simulations are carried out on a two‑dimensional rectangular domain (≈200 µm × 400 µm) with the vessel placed along the lower edge. An initial small perturbation of φ is introduced near the vessel. Over time, the simulations reproduce the hallmark features of tumor cords: a narrow, high‑φ proliferative band hugging the vessel, and a downstream region where nutrient concentration falls below c* and φ sharply declines, forming a necrotic core. Parameter sweeps illustrate that increasing the nondimensional nutrient supply shrinks the necrotic region, whereas decreasing it expands the dead zone and can eventually eradicate the cord.

Overall, the paper demonstrates that a reduced two‑equation porous‑media model can capture the essential physics of tumor‑cord growth—mechanical coupling, nutrient‑driven proliferation, and spatial heterogeneity—while remaining analytically tractable. By bridging detailed multiphase mixture theory with a parsimonious PDE system, the work offers a solid foundation for future extensions such as ECM remodeling, angiogenesis, and drug delivery modeling, thereby providing a valuable tool for both theoretical oncology and the design of therapeutic strategies.