Error analysis of a fully discrete structure-preserving finite element scheme for a diffuse-interface model of tumour growth
Fuente:
arXiv
Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866908605549117440 |
|---|---|
| author | Soenjaya, Agus L. Lin, Ping Tran, Thanh |
| author_facet | Soenjaya, Agus L. Lin, Ping Tran, Thanh |
| contents | We develop a linear fully discrete structure-preserving finite element method for a diffuse-interface model of tumour growth. The system couples a Cahn--Hilliard type equation with a nonlinear reaction-diffusion equation for nutrient concentration and admits a dissipative energy law at the continuous level. For the discretisation, we employ a scalar auxiliary variable (SAV) formulation together with a mixed finite element method for the Cahn--Hilliard part and standard conforming finite elements for the reaction-diffusion equation in space, combined with a first-order Euler time-stepping scheme. The resulting method is unconditionally energy-stable, mass-preserving, and inherits a discrete energy dissipation law associated with the SAV-based approximate energy functional, while requiring the solution of only linear systems at each time step. Under suitable regularity assumptions on the exact solution, we derive rigorous error estimates in $L^2$, $H^1$, and $L^\infty$ norms, establishing first-order accuracy in time and optimal-order accuracy in space. A key step in this analysis is the proof of boundedness of the numerical solutions in $L^\infty$. Numerical experiments validate the theoretical convergence rates and demonstrate the robustness of the method in capturing characteristic phenomena such as aggregation and chemotactic tumour growth. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2509_14486 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Error analysis of a fully discrete structure-preserving finite element scheme for a diffuse-interface model of tumour growth Soenjaya, Agus L. Lin, Ping Tran, Thanh Numerical Analysis 65M12, 65M15, 65M60, 35Q92 We develop a linear fully discrete structure-preserving finite element method for a diffuse-interface model of tumour growth. The system couples a Cahn--Hilliard type equation with a nonlinear reaction-diffusion equation for nutrient concentration and admits a dissipative energy law at the continuous level. For the discretisation, we employ a scalar auxiliary variable (SAV) formulation together with a mixed finite element method for the Cahn--Hilliard part and standard conforming finite elements for the reaction-diffusion equation in space, combined with a first-order Euler time-stepping scheme. The resulting method is unconditionally energy-stable, mass-preserving, and inherits a discrete energy dissipation law associated with the SAV-based approximate energy functional, while requiring the solution of only linear systems at each time step. Under suitable regularity assumptions on the exact solution, we derive rigorous error estimates in $L^2$, $H^1$, and $L^\infty$ norms, establishing first-order accuracy in time and optimal-order accuracy in space. A key step in this analysis is the proof of boundedness of the numerical solutions in $L^\infty$. Numerical experiments validate the theoretical convergence rates and demonstrate the robustness of the method in capturing characteristic phenomena such as aggregation and chemotactic tumour growth. |
| title | Error analysis of a fully discrete structure-preserving finite element scheme for a diffuse-interface model of tumour growth |
| topic | Numerical Analysis 65M12, 65M15, 65M60, 35Q92 |
| url | https://arxiv.org/abs/2509.14486 |