N-S方程隐式分区并行计算.pdfVIP

cHINEsE oFcoMPuTATIoNAL 2001年5月荛。毒品一 cH。NE。E，。u篡AL嚣c。M：k焉NA。，HY。。c。JouRNAL PHYslcs％：?j揣i：：。{ioi 019907 [文章编号]1001．246x(2001)03 N．S方程隐式分区并行计算 邹 辉1， 李劲菁1， 吴子牛2 z [摘要】将wu N和zouH提出的一套适合三点块三对角隐式格式的定常和非定常可压无粘流的覆盖分区并 行计算方法推广到Ns方程的计算在二维定常亚音速、跨音速和非定常情况下都得到了较好的绝对并行效率 [关键词]覆盖耦合条件；肝行计算；Nav·盱s佃kes方程 1；V211l3 [中图分类号]、恐11 [文献标识码】A 区处理的收敛一样快；而对于上游差分格式这一类 0 引言 强耗散格式，覆盖长度几乎不影响收敛速度，考虑到 目前，并行计算在计算流体力学(cFD)占有非汁算时间，应取小的覆盖长度 常重要的地位在CFD问题里进行并行最自然的方 对于非定常问题，时间滞后处理一定会削弱交 法就是分区【3，…．这对下显式格式来说相对容易，因 界面附近的时间精度．对于类似于三对角并行解法 为显式格式被定义在点方向t是自然并行的但对 于隐式格式，离散方程是空间相关的，显然情况更复 schur补充法，或其它类似的方法，可能会增加计算 杂．本文只考虑在每一空曲】方向上牵涉到三点的三 时间或增加算法的复杂性．对于离散问题，波峰以有 对角隐式格式．在高维问题中，隐式系统会通过近似 限的速度移动．既然cFL数是基于最大波峰速度 因子法分裂为一维系统 (特征值)，那么格式的当地扰动(因为时间滞后)移 对隐式系统来说，并行计算的主要困难是如何 动的距离，以网格点数目的形式表示，不会比CFL 用并行的方法来对隐式系统求逆传统上，我们使用 数大根据这一点，有可能选择与cFL数成比例的 覆盖长度，并建立起容易实现和保证时间精确的耦 并行三对角解法或schwam和schur补充法_41在 并行三对角解法里，用于对三对角系统求逆的传统 合条件． 的托马斯算法被修改了schwar【z法是基于有长度 wuzN和zouH最近提出了一套适合三点块 的分区，并且在每一时间步上需要一些叠代以对隐 t对角隐式格式的定常和非定常可压无粘流的覆盖 分区并行计算方法，并对方法的性能进行了分析证 式系统求逆schur补充法不需要子区在交界面上 有覆盖，但需要与交界面相关的子系统的解．本文使 明_2J．这套方法如果推广到N—S方程，对定常问题 用覆盖多块网格来研究并行计算，考虑的覆盖法不 的收敛性和非定常问题的时间精度有什么影响，是 象schur算法，因为后者在每一时间步上需要叠代．本文要讨论的这种方法对于定常和非定常问题都 对于定常和非定常问题，在交界面上使用覆盖的原 不需要在每一时间步上叠代，也不需要对隐式部分 因是不同的 的托马斯算法进行修正．在本文中，我们对数值效率 对于定常状态的问题，最简单的方法是使用时 比对并行的其它特性更感兴趣 间滞后耦合条件．这尼然很简单，但似乎会减慢收敛 到定常的速度对于特定的格式，大量的研究1。2得 1适合并行计算的覆盖分区算法 出了一个使人吃惊的结果：根据正规模态分析，对于