Identifying early tumour states in a Cahn-Hilliard-reaction-diffusion model

TL;DR

This paper develops a mathematical and computational framework to reconstruct early tumour states from a single late-time measurement using a coupled Cahn-Hilliard-reaction-diffusion model, addressing the ill-posed inverse problem.

Contribution

It introduces a novel approach combining uniqueness proofs, Tikhonov regularisation, and finite element methods for reconstructing tumour evolution in a complex non-linear model.

Findings

The method accurately reconstructs tumour states in test cases.

The proposed regularisation approach stabilizes the inverse problem.

Numerical results demonstrate the effectiveness of the algorithm.

Abstract

In this paper, we tackle the problem of reconstructing earlier tumour configurations starting from a single spatial measurement at a later time. We describe the tumour evolution through a diffuse interface model coupling a Cahn-Hilliard-type equation for the tumour phase field to a reaction-diffusion equation for a key nutrient proportion, also accounting for chemotaxis effects. We stress that the ability to reconstruct earlier tumour states is crucial for calibrating the model used to predict the tumour dynamics and also to identify the areas where the tumour initially began to develop. However, backward-in-time inverse problems are well-known to be severely ill-posed, even for linear parabolic equations. Moreover, we also face additional challenges due to the complexity of a non-linear fourth-order parabolic system. Nonetheless, we can establish uniqueness by using logarithmic convexity methods under suitable a priori assumptions. To further address the ill-posedness of the inverse problem, we propose a Tikhonov regularisation approach that approximates the solution through a family of constrained minimisation problems. For such problems, we analytically derive the first-order necessary optimality conditions. Finally, we develop a computationally efficient numerical approximation of the optimisation problems by employing standard $C^0$-conforming first-order finite elements. We conduct numerical experiments on several pertinent test cases and observe that the proposed algorithm consistently meets expectations, delivering accurate reconstructions of the original ground truth.