SAMR网格上扩散方程有限体格式的逼近性和两层网格算法.pdfVIP

第4期 舒 适等：SAMR网格上扩散方程有限体格式的逼近性与两层网格算法 391 地减少网格规模 ，节省计算量和计算资源．目前，SAMR自适应计算 已经广泛应用于辐射流体力学 、天体 物理H “等实际问题 的数值模拟 中．并且 ，基于 SAMR 自适应计算的应用软件支撑框架也得到快速发 展 ． 对于流体力学方程组的显式离散格式 ，Berger等人提出的基于细化模式的时间积分相对成熟 I4．在该 模式 中，时间积分从最粗到最细 网格层逐层进行 ，相对于粗 网格层 的一个时间步，细网格层基于粗 网格层提 供 的边界条件 ，积分若干个时间步．对于扩散方程的隐式离散格式 ，通常采用组合模式进行时间积分Ⅲ ． 在组合模式中，所有 网格层组合在一起形成组合 网格 ¨ ，空间离散和时间积分均在组合 网格上进行 ，此时， 所有网格层具有相 同的时间步长．多物理耦合问题的数值模拟和大规模并行计算表明，组合模式具有很强的 局 限性 ，因此 ，近年来也发展 了求解扩散方程的细化模式隐式时间积分算法 ’ ． 给定 SAMR网格 ，无论采用细化模式还是组合模式 ，扩散方程的求解主要涉及两个关键问题 ：离散格式 的构造和离散代数系统的求解 ．对于细化模式 ，离散格式只需针对单层均匀网格进行设计 ，但需要基于粗细 界面处 的细化插值 (为细网格层提供边界条件)和加密区域 的粗化插值 (为粗 网格层提供校正量)在网格层 之间进行多层迭代 ，直到收敛 ]．对于组合模式，离散格式于组合网格需要处理粗细界面的 “一对多”非协调 单元 ．对非协调单元的离散格式 ，本质上也可视为在粗细界面处通过引入细化插值和粗化插值进行构造． 因此 ，在多层迭代收敛的前提下，细化模式与组合模式在离散格式层面是等价的．简单起见 ，本文仅针对组合 模式． 由于本文的讨论仅涉及到SAMR网格上的空间离散格式，为简单起见，考虑如下二维定常扩散问题 r一7·(，c7 “)+AM=f in = (0，0)×(0，b)， {v【--：g 。。n00, 其 中，l为a 的单位外法向量 ，o和b均为正常数 ，，c=diag(，c。，，c)，这里K、K：、A、f和g均为已知的连续函 1 SAMR网格上的有限体格式 1．1 混合五点格式 设Q 为对 作n ×ny等距剖分得到的矩形网格(如图2所示)， 其沿 和 Y方向的剖分步长分别为 hx：a／nx和 hy=b／ny，h= (f， 1) ● max(hx，hy)为Q 的尺寸．记剖分单元 ，，= ( 一hx／2， +hx／2)× Q。 (Yj—hy／2，Yj+hy／2)，其 中 =(i一1／2)hx，Yj= (一1／2)hy，i= ● (i-1，) (f√) (-1 1，2，…，n ， =1，2，…，ny．若d na =咖，则称力 为 内部单元 ， ， ． ● 否则称为边界单元．