任一鹏 刘 枫 林 其 胡华伟

∗(北京宇航系统工程研究所，北京100076)

†(中国空气动力研究与发展中心，四川绵阳621000)

∗∗(成都铭峰新源科技有限公司，成都610042)

真实飞行器外形复杂，所涉及的流动结构和现象对数值模拟中所需的计算网格提出了较高的要求。一方面，飞行器本身构型多样，各种部件无论从尺度还是位置都存在着多种变化和组合，将传统的参数化的网格生成方法应用其中存在重大障碍；另一方面，真实飞行条件多变，所产生的流动结构性质各异，对计算网格的要求也不尽相同。在流动的数值模拟中，一个熟知的原则就是网格分布应该尽可能符合当地的流动状态。因此，建立能够适应复杂外形的网格生成方法是必要的。

国内外一些学者在复杂外形的网格生成方法[1-2]上进行了诸多探索。在过去的二十年中，许多重要的网格生成技术得到建立和发展，比如多块对接/拼接结构网格、重叠/嵌套网格、非结构(四面体/金字塔)网格等。近年来混合网格(包含有多种形态的网格)越来越受到关注，得到了较快的发展。Zhang等[3-4]建立了基于各向异性四面体凝聚的混合网格生成方法，并成功运用于大雷诺数的复杂飞行器数值模拟中；Esquieu[5]和Kallinderis等[6]建立了三棱柱/四面体混合网格生成算法，并运用于喷气式客机绕流流场的计算；Timothy[7]建立了三棱柱/金字塔/四面体/六面体混合网格生成方法；Kannan等[8]和Wang[9]建立了自适应笛卡尔网格生成算法。众所周知，边界层网格的质量对计算结果有重大影响。……