A structure-preserving LDG discretization of the Fisher-Kolmogorov equation for modeling neurodegenerative diseases
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_ | 1866915510191390720 |
|---|---|
| author | Antonietti, Paola F. Corti, Mattia Gómez, Sergio Perugia, Ilaria |
| author_facet | Antonietti, Paola F. Corti, Mattia Gómez, Sergio Perugia, Ilaria |
| contents | This work presents a structure-preserving, high-order, unconditionally stable numerical method for approximating the solution to the Fisher-Kolmogorov equation on polytopic meshes, with a particular focus on its application in simulating misfolded protein spreading in neurodegenerative diseases. The model problem is reformulated using an entropy variable to guarantee solution positivity, boundedness, and satisfaction of a discrete entropy-stability inequality at the numerical level. The scheme combines a local discontinuous Galerkin method on polytopal meshes for the space discretization with a $ν$-step backward differentiation formula for the time integration. Implementation details are discussed, including a detailed derivation of the linear systems arising from Newton's iteration. The accuracy and robustness of the proposed method are demonstrated through extensive numerical tests. Finally, the method's practical performance is demonstrated through simulations of $α$-synuclein propagation in a two-dimensional brain geometry segmented from MRI data, providing a relevant computational framework for modeling synucleopathies (such as Parkinson's disease) and, more generally, neurodegenerative diseases. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2504_05784 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | A structure-preserving LDG discretization of the Fisher-Kolmogorov equation for modeling neurodegenerative diseases Antonietti, Paola F. Corti, Mattia Gómez, Sergio Perugia, Ilaria Numerical Analysis 65M60, 65N22, 35Q92 This work presents a structure-preserving, high-order, unconditionally stable numerical method for approximating the solution to the Fisher-Kolmogorov equation on polytopic meshes, with a particular focus on its application in simulating misfolded protein spreading in neurodegenerative diseases. The model problem is reformulated using an entropy variable to guarantee solution positivity, boundedness, and satisfaction of a discrete entropy-stability inequality at the numerical level. The scheme combines a local discontinuous Galerkin method on polytopal meshes for the space discretization with a $ν$-step backward differentiation formula for the time integration. Implementation details are discussed, including a detailed derivation of the linear systems arising from Newton's iteration. The accuracy and robustness of the proposed method are demonstrated through extensive numerical tests. Finally, the method's practical performance is demonstrated through simulations of $α$-synuclein propagation in a two-dimensional brain geometry segmented from MRI data, providing a relevant computational framework for modeling synucleopathies (such as Parkinson's disease) and, more generally, neurodegenerative diseases. |
| title | A structure-preserving LDG discretization of the Fisher-Kolmogorov equation for modeling neurodegenerative diseases |
| topic | Numerical Analysis 65M60, 65N22, 35Q92 |
| url | https://arxiv.org/abs/2504.05784 |