Identifying early tumour states in a Cahn-Hilliard-reaction-diffusion model
Fuente:
arXiv
Saved in:
| Main Authors: | , , , , |
|---|---|
| Format: | Preprint |
| Published: |
2024
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866913514617044992 |
|---|---|
| author | Agosti, Abramo Beretta, Elena Cavaterra, Cecilia Fornoni, Matteo Rocca, Elisabetta |
| author_facet | Agosti, Abramo Beretta, Elena Cavaterra, Cecilia Fornoni, Matteo Rocca, Elisabetta |
| contents | 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. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2409_15925 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Identifying early tumour states in a Cahn-Hilliard-reaction-diffusion model Agosti, Abramo Beretta, Elena Cavaterra, Cecilia Fornoni, Matteo Rocca, Elisabetta Analysis of PDEs Numerical Analysis Optimization and Control 35G31, 35Q92, 35R30, 49K20, 65M32, 92C50 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. |
| title | Identifying early tumour states in a Cahn-Hilliard-reaction-diffusion model |
| topic | Analysis of PDEs Numerical Analysis Optimization and Control 35G31, 35Q92, 35R30, 49K20, 65M32, 92C50 |
| url | https://arxiv.org/abs/2409.15925 |