覆盖分区算法的理论基础和应用.pdfVIP

覆盖分区算法的理论基础与应用 吴子牛，邹辉。 (北京航空航天大学一国家计算流体力学实验室，北京100083) 摘要本文对覆盖分区算法所涉及的各种理论问题和应用问题进行总结和稚广 1引言 覆盖分区方法在计算流体力学中被用来解决绕复杂物体的流动问题。覆盖分区方法又 可以方便地用来进行并行计算。尤其在使用隐式格式时，各区之间的联立迭代通过覆 盖可以加快收敛速度。 在覆盖分区计算中，子区之间的数据交换是通过内边界上的插值即耦合条件实现 的，对精度、稳定性、守恒性、收敛性以及解的唯—性等均有影响。另外，在覆盖分区 网格上对隐式格式实行并行计算时，因必须对耦合条件进行时阿滞后处理，从而影响 定常问题的收敛速度和非定常问题的时间精度。 这里讨论和推广该方法的某些重要理论研究结论，同时结合理论研究将覆盖分区方 法应用到并行计算和航空工程问题的数值模拟中。主要涉及如下几个问题： 11线性研究的综合讨论。线性研究包含稳定性、收敛性、解的唯一性，等。这些问 题已经在文献【2，3，4，5]中单独研究，研究结果已经在文献[7】中详细介绍。这里主要 揭示稳定性、收敛性和解的唯—性之间的统一关系，并且阐明这些特性与守恒性之间 是相互矛盾的。 2)非线性研究的成果的讨论与推广。线性研究指明：守恒的覆盖分区算法稳定性、 收敛性和解的唯一性等很差。因此似乎不存在能同时保证守恒、稳定、收敛和解唯一 的覆盖分区算法。这与大量实际工程问题的成功计算相矛盾。为了解决这一既有理论 意义又有实用价值的问题，文献[6]6对与守恒相关联的激波与分区界面的相互干扰进行 了非线性研究，得出了一条重要结论：如果在激波内格式的数值粘性不小于Roe格式 的数值粘性，则在分区边界上可使用非守恒插值而不影响守恒性。这一结论为构造守 恒稳定的覆盖分区算法指明了方法。文献[61的理论结论是以一维标量方程为依据的， 这里将最主要的结论推广到双曲系统。 3)分区与并行计算。当对隐式格式进行分区并行处理时，在分区边界上的简化处理 (显式化)会严重影响收敛到定常解的收敛速度和非定常问题的时间精度。这里以文献 f4】的最优覆盖长度理论为起点，实现具有最佳并行效率的覆盖分区并行算法。 2覆盖分区分析模型与耦合条件 一般可压缩流动问题由三维Ns方程tifi述。而影响本文感兴趣问题的关键来自方程的非 ■ 线性部分。这里只介绍一维欧拉方程的分析模型。它属于一般双曲系统： tm+，(b=0，t0，一1zl (1) ·由自然科学基金，国家教委优秀青年教师基金与航空基金资助。 239 其中”为未知量，，=，(u)为通量函数且其雅可比矩阵可以对角化。因上式可以对角 化，所以在大部分情况下可以用标量方程代替。为了唯一确定方程(1)的解，还需给出 适当的初始条件和外边界条件。对于非定常问题，任意时刻的解都依赖初始条件。对 于定常问题，通常把该问题祝为非定常问题，在时闻足够大且解不再随时间变化时， 便认为获得了定常解。解的唯—性要求定常解与初始条件无关。 原始的单区算法将上述初边值问题构造统一网格进行求解。在覆盖分区算法中，将 整体区域分解成两个(或更多)的相互重叠的子区，在每个子区上独立构造网格和求 解。 ％ -筹越截； 讲矿—弋r≮芦i 图1：覆盖分区网格示意图。 覆盖分区算法具有两个明显特征：首先是流场的分块独立求解和解的耦合，其次是 在覆盖区流场重复计算。 以两个