Many neurodegenerative diseases are connected to the spreading of misfolded prionic proteins. In this paper, we analyse the process of misfolding and spreading of both $\alpha$-synuclein and Amyloid-$\beta$, related to Parkinson's and Alzheimer's diseases, respectively. We introduce and analyze a positivity-preserving numerical method for the discretization of the Fisher-Kolmogorov equation, modelling accumulation and spreading of prionic proteins. The proposed approximation method is based on the discontinuous Galerkin method on polygonal and polyhedral grids for space discretization and on $\vartheta-$method time integration scheme. We prove the existence of the discrete solution and a convergence result where the Implicit Euler scheme is employed for time integration. We show that the proposed approach is structure-preserving, in the sense that it guaranteed that the discrete solution is non-negative, a feature that is of paramount importance in practical application. The numerical verification of our numerical model is performed both using a manufactured solution and considering wavefront propagation in two-dimensional polygonal grids. Next, we present a simulation of $\alpha$-synuclein spreading in a two-dimensional brain slice in the sagittal plane. The polygonal mesh for this simulation is agglomerated maintaining the distinction of white and grey matter, taking advantage of the flexibility of PolyDG methods in the mesh construction. Finally, we simulate the spreading of Amyloid-$\beta$ in a patient-specific setting by using a three-dimensional geometry reconstructed from magnetic resonance images and an initial condition reconstructed from positron emission tomography. Our numerical simulations confirm that the proposed method is able to capture the evolution of Parkinson's and Alzheimer's diseases.

翻译：许多神经退行性疾病与错误折叠的朊蛋白传播有关。本文分析了分别与帕金森病和阿尔茨海默病相关的α-突触核蛋白和β-淀粉样蛋白的错误折叠与传播过程，提出并分析了一种保持正性的数值方法，用于离散模拟朊蛋白积累与传播的Fisher-Kolmogorov方程。该近似方法基于多边形/多面体网格上的间断伽辽金方法进行空间离散，并采用ϑ-方法时间积分格式。我们证明了离散解的存在性，以及采用隐式欧拉格式进行时间积分时的收敛性结果。研究表明该方法具有结构保持特性，即能确保离散解的非负性，这一特征在实际应用中至关重要。数值验证通过制造解方法及二维多边形网格中的波前传播测试完成。进一步，我们模拟了矢状面二维脑切片中α-突触核蛋白的传播过程：利用多面体间断伽辽金方法在网格构建中的灵活性，通过聚合网格保持白质与灰质区域的区分。最后，基于磁共振图像重建的三维几何结构与正电子发射断层扫描重建的初始条件，我们模拟了患者特异性模型中β-淀粉样蛋白的传播过程。数值模拟结果证实，该方法能够有效捕捉帕金森病与阿尔茨海默病的演变过程。