网格自适应与并行计算在气动力计算中的应用
计算的准确性与效率一直是计算流体力学研究的重点之一。本文从网格自适应与并行计算两方面对计算的准确性与效率进行了研究。针对自适应方法中的关键问题，提出了新的捕捉流场中关键区域的方法;针对现存串行程序与硬件条件，发展了一种非结构并行计算算法，并结合自适应技术发展了一套并行自适应方法。 基于非结构有限体积法，发展了一套适用于不可压/可压、定常/非定常、无粘/有粘、二绀/三维流场的数值模拟求解器，数值计算了二维翼型、三维机翼与翼身组合体的气动力。空间离散采用二阶精度的Jameson格心格式。时间离散，对定常流动采用五步Runge-Kutta与LU-SGS方法;对非定常流动采用双时间步长法。对流通量采用JST格式。湍流模型选取一方程SA模型。合理运用了一些加速收敛措施，并针对不可压流引入低速预处理加速收敛技术。最后，通过算例对准确性、收敛性进行了验证。 为了提高计算准确性、减少计算资源，采用了h型网格自适应技术。首先给出二维四边形、三维六面体网格为初始稀疏网进行流场的计算。在初始网格流场解的基础上，本文提出速度紊乱度构造指示函数对粘性流动中的附面层、湍流区与尾迹区进行捕捉并局部加密网格。为了简洁高效地捕捉激波与涡结构区域，文中提出以单位体积总能为指示变量构造统一的指示函数捕捉上述区域。对于局部网格的加密，二维情况下，一个四边形网格单元将被切分为四个小的四边形网格单元;三维情况下，一个六面体网格单元将被切分为八个小的六面体网格单元。但是对二维仨维物面第一层网格加密时，将不对法向方向进行切分。自适应后的网格不同层级网格间会出现悬挂点。文中采用按边循环的非结构求解器将解决悬挂点问题，同时保证了边/面上的通量守恒。新产生的网格单元处的流场值从初始流场解插值获得。这样自适应后的流场计算可在初始流场解的基础上进行，节约了计算时间。 针对三维数值模拟中网格数量巨大的特点与目前的计算机硬件条件，流场迭代计算中使用了并行计算方法。本文提出只在流场迭代部分采取并行措施，其余部分全部由0进程完成;初始网格与自适应后网格的划分采用在程序中调用METIS的方法;并行计算时，进程间信息的交换采用MPICH2函数接口。重点研究了分区边界处理策略，数据高效交换及快速的输入输出问题。METIS网格划分软件开源且易于获取，针对不同的网格数提供不同的调用接口，具有分区质量高、速度快、最小化子域边界数量的优点。MPICH2是目前主流的并行函数接口标准。因此上述方法对现有串行程序的改造简单快速，适用性广。由于自适应前后都保证了负载的平衡，此方法具有较高的并行效率。通过算例对并行自适应程序进行了验证。
非结构有限体积法;网格自适应;并行计算;气动力
南京航空航天大学
博士
流体力学
陈迎春;陆志良
2013
中文
O354
145
2014-06-05（万方平台首次上网日期，不代表论文的发表时间）