三维流体动力学并行软件包NaSt3DGP及其应用.pdf

　　中　国　科　学　院　上　海　天　文　台　年　刊 ２００７年　总第２８期　　ＡＮＮＡＬＳＯＦＳＨＡＮＧＨＡＩＡＳＴＲＯＮＯＭＩＣＡＬＯＢＳＥＲＶＡＴＯＲＹ，ＣＡＳ 　 Ｎｏ．２８，　２００７ 三维流体动力学并行软件包 ＮａＳｔ３ＤＧＰ及其应用 常英 立 （上海水产大学，上海２０００９０） 提　　　要 该文介绍了解算直角坐标系下不可压缩 ＮａｖｉｅｒＳｔｏｋｅｓ方程的大型有限差分并行软件包 ＮａＳｔ３ＤＧＰ。介绍了其采用的区域分解算法、差分格式及计算流程；对并行效率进行了测试，作出了 处理器数目和计算效率的关系曲线；把此程序应用到旋转的行星流体动力学中，分别用自检验和互 检验的方法证明了修改后程序的可靠性，并且得到了满意的非线性数值模拟结果。 主题词：：流体动力学 — 软件包 — 区域分解 — 有限差分 — 数值模拟 分类号：ＴＰ３１７，Ｏ３５１．２ １　引　　言 大型有限差分并行软件包 ＮａＳｔ３ＤＧＰ由德国波恩大学 ＭｉｃｈａｅｌＧｒｉｅｂｅｌ、ＲｏｂｅｒｔｏＣｒｏｃｅ、 ＦｒａｎｋＫｏｓｔｅｒ和ＭｉｃｈａｅｌＭｅｙｅｒ用Ｃ＋＋语言编写而成，主要用于解算三维直角坐标系下不可压缩 的ＮａｖｉｅｒＳｔｏｋｅｓ方程。这是一个开源软件，很容易进行修改，以便应用于求解更复杂的流体动力 ［１］ 学问题。例如，在行星流体动力学研究中，应用此软件可以得到行星流场的非线性数值解 。 ２　程序的实现及使用方法 该程序求解的方程为： ｕ ｇ １  ２ ＋ｕ· ｕ＝ － ｐ＋ ｕ （１）    ｔ Ｆｒ Ｒｅ  ·ｕ＝０ （２） Θ １ ２ ＋ｕ·Θ ＝ Θ （３） ( ｔ ) ＲｅＰｒ  式（１）是速度场方程。式中ｕ表示流体的运动速度；ｇ表示重力加速度；Ｆｒ称为弗劳德 （Ｆｒｏｕｄｅ）数，其物理意义是表示惯性力与重力之比，流体动力相似的流动其 Ｆｒ值相同；Ｒｅ称 为雷诺（Ｒｅｙｎｏｌｄｓ）数，其物理意义为惯性力与粘性力的比值。两种流动的粘性力作用相似，它 收稿日期：２００７０６０８；　修回日期：２００７０６２９ 总第２８期　　　　　　　　三维流体动力学并行软件包ＮａＳｔ３ＤＧＰ及其应用 １７３ 们的雷诺数必定相等。式（２）表示不可压缩流的连续性方程。式（３）表示温度（或能量）方 程。式中 表示热力学温度，无量钢Ｐｒ称为普朗特（Ｐｒａｎｄｔｌ）数，它表征了流体中动量扩散与 Θ 热量扩散能力的比值。此程序可以求解出任意时刻的速度场和温度场。 ＮａＳｔ３ＤＧＰ采用的是并行算法。并行计算的原理是：整个计算区域分割成多个子区域，然 后分配给多个处理器同时计算。子区域的初始信息和几何信息分别装载入各子区域对应处理 器的内存中，在每一个处理器中启动计算进程