高阶并行差分格式解抛物型方程的算法研究

需积分: 9 0 下载量 134 浏览量 更新于2024-08-12 收藏 248KB PDF 举报
"抛物型方程的一种高阶并行差分格式 (2009年)" 抛物型方程在物理学、工程学以及许多其他领域都有广泛应用,它们常常用来描述热传导、流体流动等现象。这篇2009年的论文主要探讨了一种新的、高阶精度的并行差分格式,用于高效地求解这类方程。 在传统的数值方法中,有限差分法是一种常用的方法来近似抛物型方程的解。然而,随着问题规模的增大,计算量会急剧增加,因此并行计算成为了解决这一问题的有效手段。该论文作者孙凯和王文洽提出了一种基于四阶紧致格式的并行差分算法,该算法能够在并行计算环境下提高计算效率,同时保持较高的数值精度。 算法的核心在于利用前三个时间层内界点的值,结合四阶紧致格式来并行计算子区域内的数值解。四阶紧致格式意味着在空间离散上使用了更高阶的插值和差分操作,从而可以提供比二阶格式更高的精度。这种方法减少了对相邻网格点数据的依赖,有利于并行化处理。 在并行计算中,区域分解是一种常见的策略,即将大的计算区域划分为多个子区域,每个子区域可以在独立的处理器上并行计算。论文中的算法首先在子区域内并行计算,然后利用区域边界点的信息显式更新内界点的值。这一过程确保了算法的连贯性和稳定性。 论文进一步证明了所提算法的稳定性条件至少为√2/3+1/6,这意味着在满足这个条件时,算法能够保持数值解的稳定,不会因为迭代过程中的误差积累而发散。此外,作者还证明了该算法的收敛精度为四阶,这意味着随着网格分辨率的提高,解与精确解之间的误差将以四阶的速度减小,这通常优于二阶格式。 为了验证算法的有效性,作者进行了数值算例的分析,结果显示,提出的并行差分格式在稳定性、收敛性和精度方面均优于其他常见的数值方法。这表明,这种新算法在解决大规模抛物型方程问题时,既能保证数值解的准确性,又能充分利用现代并行计算资源,从而提高了计算效率。 这篇论文为抛物型方程的高效并行求解提供了一个创新的数值方法,对于数值模拟和科学计算领域具有重要的理论价值和实际应用前景。其提出的高阶并行差分格式在处理大型复杂问题时,有望成为一种有力的工具。