求解二维半线性抛物方程的校正型显隐区域分解算法.pdfVIP

中国科学 A 辑 : 数学 2009 年 第 39 卷 第 6 期 : 749 ∼ 774 求解二维半线性抛物方程的校正型显隐区域 分解算法 ①②∗ ② ① 廖洪林 , 史汉生 , 孙志忠 ① 东南大学数学系, 南京 210096 ② 中国人 民解放军理工大学理学院应用数学与物理系, 南京 211101 E-mail: liaohl2003@, zzsun@ 收稿 日期: 2007-09-03; 接受 日期: 2008-11-21; * 通信作者 国家 自然科学基金 (批准号: 资助项 目 摘要 本文研究了利用分布式并行计算系统求解二维半线性抛物方程的内边界校正型显 隐区域分解 (CEIDD) 算法. 在实际问题 中通常利用简洁的直线 内边界 (SI) 将空间区域分 解成若干个相互不重叠的条状或块状子 区域. 利用 Leray-Schauder 不动点定理和离散能量 方法证 明了基于不交叉直线 内边界的 CEIDD-SI 算法的唯一可解性, 无条件稳定性和收敛 性, 并得到 了一个改进 的误差估计. 当直线 内边界在 区域 内部相互交叉时, 这种在 内边界上 追加 了隐式校正步的算法需要在每一个 时间层进行全局通信, 从而使算法的并行可扩展性 大为降低. 为克服这一缺点, 设计 了一种 由直线和锯齿形接点组合而成的复合 内边界 (CI). 分析表 明, 基于复合 内边界 的 CEIDD-CI 算法无条件稳定、通信效率高、可 以直接利用现 有的串行算法计算子 区域 的隐式解, 是一类可扩展 的并行算法. 为验证算法的稳定性和收 敛性, 文 中给 出了两个具体算例. 关键词 半线性抛物方程 显隐区域分解算法 Leray-Schauder 不动点定理 离散能量 法 收敛性与稳定性 MSC(2000) 主题分类 65M06, 65M12, 65M55, 68Y05 1 引言 区域分解是设计偏微分方程并行解法的重要技术, 已有很多的文献 [1, 2] 讨论定常和非 定常问题 的区域分解算法. 对于抛物方程, 显式格式是本性并行 的, 且很容易在并行计算系 统 中实现, 但时间积分步长往往受苛刻 的稳定性条件所 限制. 通常, 隐式格式没有稳定性条 件限制, 常用于非线性发展方程 的数值计算. 但是, 由于在每个时间层上需要求解一个椭 圆 型问题, 隐格式不是本性并行 的, 因而不能直接在并行系统 中应用. 隐格式并行化的一种 自然而然的策略是充分利用显隐交替格式. 周等 [3−6] 提出的所谓 本性并行差分格式就是这样一类方法. 他们构造 了一系列带参数 的本性并行差分格式. 适 引用格式: 廖洪林, 史汉生, 孙志忠. 求解二维半线性抛物方程的校正型显隐区域分解算法. 中国科学 A, 2009, 39(6): 749–774 Liao H L, Shi H S, Sun Z Z. Corrected explicit-implicit domain decomposition algorithms for two- dimensional semilinear parabolic equations. Sci China Ser A, 2009, 52, DOI: 10.1007/s11425-009-0040-8 廖洪林等: 求解二维半线性抛物方程的校正型显隐区域分解算法 当选择这些参数可以得到诸如交替分组显格式 [7] 、交替分块显隐格式 [8] 以及分块 ADI 格 式 [9]等具体的差分算法. 由于具有 良好 的稳定性, 可以使用大步长进行时间积分, 这些差分 格式是有效 的并行算法. 但是, 正如文献 [10, 11] 所指出, 这些方法往往不能直接利用现有 的串行计算程序. 隐格式并行化的另一种策略是将计算区域分解成若干个互相不重叠的小区域, 从前面几 个时间层的数值解 中获取子区域之间的界面信息, 而后在每个