计算流体动力学中网格生成技术的发展.doc

计算流体动力学中网格生成技术的发展 计算流体动力学作为计算机科学、流体力学、偏微分方程数学理论、计算几何、数 值分析等学科的交叉融合，它的发展除依赖于这些学科的发展外，更直接表现于对网格生 成技术、数值计算方法发展的依赖。 在计算流体动力学中，按照一定规律分布于流场中的离散点的集合叫网格（Grid） ，分布这些网格节点的过程叫网格生成（Grid Generation）。网格生成对CFD至关重要， 直接关系到CFD计算问题的成败。1974年Thompson等提出采用求解椭圆型方程方法生成贴体 网格，在网格生成技术的发展中起到了开创作用。随后Steger等又提出采用求解双曲型方 程方法生成贴体网格。但直到二十世纪八十年代中期，相比于计算格式和方法的飞跃发展 ，网格生成技术未能与之保持同步发展。因而从二十世纪八十年代开始，各国计算流体和 工业界都十分重视网格生成技术的研究。二十世纪九十年代以来迅速发展的非结构网格和 自适应笛卡尔网格等方法，使复杂外形的网格生成技术呈现出了更加繁荣发展的局面。现 在网格生成技术已经发展成为CFD的一个重要分支，它也是计算流体动力学近二十年来一个 取得较大进展的领域。也正是网格生成技术的迅速发展，才实现了流场解的高质量，使工 业界能够将CFD的研究成果——求解Euler/NS方程方法应用于型号设计中。 随着CFD在实际工程设计中的深入应用，所面临的几何外形和流场变得越来越复杂， 网格生成作为整个计算分析过程中的首要部分，也变得越来越困难，它所需的人力时间已 达到一个计算任务全部人力时间的60%左右。在网格生成这一“瓶颈”没有消除之前，快速 地对新外形进行流体力学分析，和对新模型的实验结果进行比较分析还无法实现。尽管现 在已有一些比较先进的网格生成软件，如ICEM、Gridgen、Gambit等等，但是对一个复杂的 新外形要生成一套比较合适的网格，其需要的时间还是比较长，而对于设计新外形的工程 人员来说，一两天是他们可以接受的对新外形进行一次分析的最大周期。CFD已经成功地缩 短了新外形设计中所需要的风洞实验时间，但在CFD对任意外形成为一种适时的分析工具以 前，新外形设计中所需要的风洞实验时间依然非常巨大。要将CFD从专业的研究团体中脱离 出来，并且能让工程设计人员应用到实际的设计中去，就必须首先解决网格生成的自动化 、即时性问题，R.Consner等人在他们的一篇文章中，详细地讨论了这些方面的问题，并提 出：CFD研究人员的关键问题是“你能把整个设计周期缩短多少天？”。而缩短设计周期的 主要途径就是缩短网格生成时间和流场计算时间。因此，生成复杂外形网格的自动化和及 时性已成为应用空气动力学、计算流体力学最具挑战性的任务之一。 当今，有众多研究人员对复杂外形的网格生成技术从分区结构网格、非结构网格和 笛卡尔网格三个不同的方向展开研究。 1．2．1分区结构网格方法 分区结构网格方法将原始的物理区域按不同的空间拓扑结构分成若干区域块，每个 单块网格的拓扑结构简单，易于生成贴体网格，然后合并这些单域贴体网格来形成复杂外 形的空间网格。常用的传统单域贴体网格生成方法可分为代数网格生成方法、求解椭圆微 分方程生成方法和求解双曲微分方程生成方法，以及从求解椭圆微分方程生成方法发展而 成的求解抛物微分方程生成方法。近十多年来发展了不少新的分区结构网格。比较成熟的 分区结构网格方法有： 1）分区对接网格方法。 2）分区重叠网格方法。 它们都属于分区网格生成方法，且与流场的分区计算方法以及并行计算方法密切相 关，并构成这些算法的基础。对于分区对接网格方法，虽然网格块之间的关系和数据交换 比较简单，但相连网格块之间必须共用网格线（二维）或网格面（三维），为了提高计算 的精度，要求块与块之间的连接应尽量光滑，这给网格生成带来诸多限制和不便。随着分 区重叠网格方法的引入，网格线（面）的公用限制被打破，网格生成的难度得以降低，但 随之而来的是网格块之间数据交换变得复杂，对于三维复杂外形绕流，划分的子块可能上 百块，合理处理块与块之间的边界就变得非常繁琐，往往需要付出大量的手工劳动。 采用结构网格方法的优势在于它易于生成物面附近的边界层网格、有许多成熟的计 算方法和比较好的湍流计算模型，因此它仍然是目前复杂外形飞行器气动力数值模拟的主 要方法，计算技术最成熟。但是比较长的物面离散时间、单块网格边界条件的确定以及网 格块之间各种相关信息的传递，又增加了快速计算分析的难度，而且对于不同的复杂外形 ，它得构造不同的网格拓扑结构，因而无法实现网格生成的“自动”，生成网格费时费力 。其发展方向是朝着减少工作量，实现网格的自动生成和自适应加密，具有良好的人机