A Conservative Finite Element ALE Scheme for Mass-Conserving Reaction-Diffusion Equations on Evolving Two-Dimensional Domains

Christopher F. Rowlatt; John A. Mackenzie; Robert H. Insall

arxiv: 1910.02282 · v1 · pith:433WMLSInew · submitted 2019-10-05 · 🧮 math.NA · q-bio.CB

A Conservative Finite Element ALE Scheme for Mass-Conserving Reaction-Diffusion Equations on Evolving Two-Dimensional Domains

John A. Mackenzie , Christopher F. Rowlatt , Robert H. Insall This is my paper

Pith reviewed 2026-05-18 09:08 UTC · model grok-4.3

classification 🧮 math.NA q-bio.CB

keywords ALE finite element methodmass conservationreaction-diffusion equationsevolving domainsmoving mesh PDEcell migrationbulk-surface coupling

0 comments

The pith

A finite element ALE scheme on moving meshes conserves total mass exactly for reaction-diffusion systems on deforming domains, independent of mesh velocity and time-step size.

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

The paper develops a conservative arbitrary Lagrangian-Eulerian finite element method for bulk-surface reaction-diffusion equations whose solutions must preserve global mass. Mass conservation is shown to hold exactly for the fully discrete solution on any valid sequence of meshes generated by the MMPDE moving-mesh procedure. The method is motivated by models of protein activation at the cell membrane that cycle between cytoplasmic and membrane-bound forms while the cell itself changes shape. Because the conservation identity does not depend on the particular ALE velocity or the size of the time step, the scheme can be applied directly to cell-migration problems driven by external chemotactic signals.

Core claim

Global conservation of the fully discrete finite-element solution is established independently of the ALE velocity field and the time-step size; numerical tests confirm that the scheme remains second-order accurate and exactly mass-conserving on deforming domains.

What carries the argument

Arbitrary Lagrangian-Eulerian finite-element discretization of bulk-surface reaction-diffusion equations on MMPDE-generated meshes that move with the evolving domain.

If this is right

The scheme can be used for any bulk-surface reaction-diffusion model that requires strict global mass conservation on a moving two-dimensional domain.
Second-order spatial and temporal accuracy is retained while mass is conserved exactly.
The method applies directly to models of single-cell migration under chemotactic gradients.

Where Pith is reading between the lines

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

Because conservation holds for arbitrary time steps, the integrator can safely take larger steps when the reaction terms become stiff.
The same conservation proof supplies a template for constructing conservative schemes on three-dimensional deforming domains once robust volumetric MMPDE generators become available.

Load-bearing premise

The MMPDE mesh generators continue to produce non-tangling bulk and surface meshes for the domain deformations that arise in the target cell-migration problems.

What would settle it

A single numerical run in which the computed total mass deviates from its initial value by more than the truncation error for some admissible ALE velocity field and time-step size.

read the original abstract

Mass-conservative reaction-diffusion systems have recently been proposed as a general framework to describe intracellular pattern formation. These systems have been used to model the conformational switching of proteins as they cycle from an inactive state in the cell cytoplasm, to an active state at the cell membrane. The active state then acts as input to downstream effectors. The paradigm of activation by recruitment to the membrane underpins a range of biological pathways - including G-protein signalling, growth control through Ras and PI 3-kinase, and cell polarity through Rac and Rho; all activate their targets by recruiting them from the cytoplasm to the membrane. Global mass conservation lies at the heart of these models reflecting the property that the total number of active and inactive forms, and targets, remains constant. Here we present a conservative arbitrary Lagrangian Eulerian (ALE) finite element method for the approximate solution of systems of bulk-surface reaction-diffusion equations on an evolving two-dimensional domain. Fundamental to the success of the method is the robust generation of bulk and surface meshes. For this purpose, we use a moving mesh partial differential equation (MMPDE) approach. Global conservation of the fully discrete finite element solution is established independently of the ALE velocity field and the time step size. The developed method is applied to model problems with known analytical solutions; these experiments indicate that the method is second-order accurate and globally conservative. The method is further applied to a model of a single cell migrating in the presence of an external chemotactic signal.

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.

Desk Editor's Note private letter to a colleague

The paper gives a practical ALE finite-element scheme for bulk-surface reaction-diffusion on deforming 2-D domains together with a clean global-conservation proof that holds independently of the mesh velocity.

read the letter

The main contribution is a fully discrete ALE finite-element method for mass-conserving bulk-surface reaction-diffusion systems on evolving domains, paired with an MMPDE moving-mesh generator. The authors prove that the scheme conserves the total mass exactly, regardless of the chosen ALE velocity or time-step size. That property is useful for the cell-polarity models they target, where total protein number must stay fixed. Numerical tests on problems with exact solutions indicate second-order accuracy and confirm the conservation statement to machine precision. The method is also run on a simple migrating-cell example with an external chemotactic cue. These are incremental but concrete advances over earlier separate treatments of conservation or moving meshes for bulk-surface problems. The chief limitation is that only the abstract is available, so the actual discrete equations, the details of the conservation proof, and the error tables cannot be checked. The robustness claim for the MMPDE mesh generator also rests on the deformations that arise in the target applications; whether it remains non-tangling for more extreme cell shapes is not yet shown. The work is aimed at numerical analysts and modelers who already work with reaction-diffusion systems on deforming domains, especially those in mathematical biology. It is solid enough to merit a serious referee report rather than a desk reject, provided the full proofs and tables hold up.

Referee Report

1 major / 1 minor

Summary. The manuscript presents an arbitrary Lagrangian-Eulerian (ALE) finite-element discretization of mass-conserving bulk-surface reaction-diffusion systems on evolving two-dimensional domains. Mesh motion is generated by a moving-mesh PDE (MMPDE) approach. The central theoretical claim is that the fully discrete solution satisfies a global conservation identity that holds independently of the ALE mesh velocity and the time-step size. Numerical experiments on problems with analytic solutions are stated to confirm second-order accuracy and exact conservation; the scheme is then applied to a model of a single migrating cell under chemotactic signaling.

Significance. If the conservation identity and second-order accuracy are rigorously established, the method supplies a structure-preserving spatial-temporal discretization that is directly applicable to the class of intracellular pattern-formation models whose global mass balance is biologically essential. The independence from mesh velocity is a non-trivial property that would distinguish the scheme from many existing ALE formulations.

major comments (1)

The abstract asserts that global conservation of the fully discrete solution is proved independently of the ALE velocity and time-step size, yet supplies neither the discrete weak form nor the telescoping argument. Without these steps the load-bearing claim cannot be verified.

minor comments (1)

The MMPDE mesh-regularity assumptions that guarantee non-tangling for the target cell-migration deformations are stated only qualitatively; a precise statement of the admissible deformation class would strengthen the applicability section.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and the constructive comment. We address the single major point below and will revise the manuscript to improve clarity of the conservation argument.

read point-by-point responses

Referee: The abstract asserts that global conservation of the fully discrete solution is proved independently of the ALE velocity and time-step size, yet supplies neither the discrete weak form nor the telescoping argument. Without these steps the load-bearing claim cannot be verified.

Authors: We agree that the abstract, being concise, does not contain the discrete weak formulation or the telescoping summation. The full proof appears in Section 3 of the manuscript, where the discrete weak form is stated (Eq. (3.7)) and global conservation is obtained by choosing the constant test function, yielding exact cancellation independent of the ALE velocity and time-step size. In the revised version we will add a short paragraph immediately after the abstract (or in the introduction) that sketches the discrete weak form and the key telescoping step so that the central claim can be verified without first reading the entire analysis section. revision: yes

Circularity Check

0 steps flagged

Derivation self-contained; conservation proved independently

full rationale

The abstract states that global conservation of the fully discrete ALE-FEM solution is established independently of mesh velocity and time-step size, presented explicitly as a theorem rather than a fitted or self-referential quantity. No equations, self-citations, or ansatzes are supplied that would allow any claimed result to reduce to its own inputs by construction. The conservation identity is the expected discrete telescoping property of a properly formulated conservative weak form and does not rely on the MMPDE mesh generator or any prior author work for its validity.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The method rests on standard finite-element theory and the well-posedness of the MMPDE mesh generator; no free parameters or invented entities are mentioned.

axioms (2)

standard math Standard Sobolev-space setting for linear finite elements on polygonal domains
Implicit in any finite-element convergence claim for second-order accuracy.
domain assumption The MMPDE produces non-degenerate meshes for the cell-shape evolutions considered
Required for the ALE mapping to remain valid throughout the simulation.

pith-pipeline@v0.9.0 · 5788 in / 1240 out tokens · 20551 ms · 2026-05-18T09:08:14.369627+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

LedgerCanonicality charge_conserved echoes

?

echoes
ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

Global conservation of the fully discrete finite element solution is established independently of the ALE velocity field and the time step size.
Atomicity sequential_preserves_conservation echoes

?

echoes
ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

The developed method is applied to model problems with known analytical solutions; these experiments indicate that the method is second-order accurate and globally conservative.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.