对流扩散问题的一种有限差分区域分解算法.pdfVIP

对流扩散问题的一种有限差分对流扩散问题的一种有限差分 对流扩散问题的对流扩散问题的一种有限差分一种有限差分 区域分解算法区域分解算法 区域分解算法区域分解算法 王传丽 (吉林大学数学研究所应用数学专业, 吉林长春) 摘摘 要要 摘摘 要要 本文提出了对流扩散问题的一种区域分解差分算法 即在区域分解的基础之上 对各个子域采用精度 较高的交替方向隐格式差分方法来求解 将所求得子域的解拚接即得整个区域的解 此方法有很强的并行 性 文中给出了该算法的收敛性分析, 数值实例的结果证实该方法是收敛的 关键词关键词 区域分解算法 交替方向隐格式 对流扩散问题 关键词关键词 1 引引 言言 引引 言言 考虑二维非稳定对流扩散问题(1) ∂u ∂2 u ∂2 u ∂u ∂u −a −a +c +c 0 ∈Ω ∈ x 2 y 2 x y (x , y ) t (0,T ) ∂t ∂x ∂y ∂x ∂y ∈∂Ω ∈ u(x , y , t) g (x , y , t) (x , y ) t (0,T ) _ u(x , y ,0) u0 (x ,y ) (x , y ) ∈Ω 其中 Ω [0,1] ×[0,1] g u0 是给定的充分光滑的函数 而a , a ,c ,c 为已知常数且 x y x y a 0, a 0, a , a 分别是x , y 方向的扩散系数 c ,c 为x , y 的对流扩散系数 对流扩散问 x y x y x y 题是水环境的数学模型及地下水资源的研究探索工作中经常遇到的模型 因此研究对流扩散 问题数值解法具有重要的实际意义 此类方程定解问题的解常常出现局部剧烈 大梯度 的 变化 如含有边界层 瞬变层等 这给问题的数值求解带来一定的困难 2 高阶交替方向隐格式方法高阶交替方向隐格式方法 高阶交替方向隐格式方法高阶交替方向隐格式方法 由于问题 1 具有小扰动的特点 使得用差分方法求解此类问题时必须采用相当精细 1 的网格剖分 这样待求节点大量增加引起求解矩阵的阶增高 导致计算量的急剧增加 另一 方面若采用一般隐式差分格式的求解矩阵在二维情形是五对角的 在多维情形 通常的显隐 格式都不方便 因此选择精度更高的 计算量较小的无条件稳定格式显得尤为重要 本文选 用文献[1]中所提出的高阶交替方向隐格式法 ADI 格式 它对时间变量具有二阶精度 而 对空间变量具有四阶精度 (下面所用记号同文献[1])