Fully explicit stabilized multirate (mRKC) methods are well-suited for the numerical solution of large multiscale systems of stiff ordinary differential equations thanks to their improved stability properties. To demonstrate their efficiency for the numerical solution of stiff, multiscale, nonlinear parabolic PDE's, we apply mRKC methods to the monodomain equation from cardiac electrophysiology. In doing so, we propose an improved version, specifically tailored to the monodomain model, which leads to the explicit exponential multirate stabilized (emRKC) method. Several numerical experiments are conducted to evaluate the efficiency of both mRKC and emRKC, while taking into account different finite element meshes (structured and unstructured) and realistic ionic models. The new emRKC method typically outperforms a standard implicit-explicit baseline method for cardiac electrophysiology. Code profiling and strong scalability results further demonstrate that emRKC is faster and inherently parallel without sacrificing accuracy.
翻译:完全显式稳定化多速率(mRKC)方法凭借其改进的稳定性特性,特别适用于求解大型多尺度刚性常微分方程系统。为证明该方法在求解刚性、多尺度、非线性抛物型偏微分方程中的数值效率,我们将mRKC方法应用于心脏电生理学中的单域方程。在此基础上,我们提出了一种针对单域模型专门优化的改进版本,即显式指数多速率稳定化(emRKC)方法。通过考虑不同有限元网格(结构化和非结构化)及真实离子模型,开展了多项数值实验以评估mRKC与emRKC方法的效率。新型emRKC方法在心脏电生理计算中通常优于标准的隐式-显式基线方法。代码性能分析和强可扩展性结果进一步表明,emRKC方法在保持精度的同时具有更快的计算速度和固有的并行性。