一种求解粘性绕流的三棱柱%2f非结构%2f直角坐标混合网格技术.pdfVIP

第十届全国高超声速气动力(热)学术交流会 一种求解粘性绕流的三棱柱／ 非结构／直角坐标混合网格技术+ 张来平呙超杨永健张涵信高树椿 中国空气动力研究与发展中心，四川绵阳621000 捅要 本文发展j’一种求解粘性绕流的三棱柱，非结构／多层次直角坐标混台网格技术． 在物面三角形生成以后。用层推进方法生成有一定压缩比的三棱柱网格．然后由四面体潮 格过渡到外场的多层次直角坐标网格。该混台网格方案充分发挥了各种网格的优越性．既 能模拟复杂外形，又提高丁整体的网格教率。在网格技术的基础j：，我们采用基于守恒变 量的隐式NND有限体积格式，求解，取椭球的高超声速绕流问题，给出r较好的结果． 表明本文的混合嗍格技术在三维复杂外形的绕流计算方面将有广阔的应用前景． 关键词 粘性流动，非结构阀恪，混合网格 0引言 为了求解高Re数的粘性绕流，人们一般采用有一定压缩比的结构网格。但是对于非 常复杂的三维外形，生成统一的贴体网格往往比较困难，大多需要采用分块或重叠技术。 非结构网格在离散复杂外形方面具有无比的优越性。但是要在牯性层内生成拉伸比很大 的非结构网格也并不容易。因此，文献“剖中提出将结构网格和菲结构网格结合起来， 构成局部或整体的半结构网格(三棱柱网格)。 本文则在作者以往发展的多层次直角坐标网格和非结构网格的无粘混合网格技术“4 的基础上，引入上述半结构网格的思想．发展了一种物面附近为三棱柱网格，然后由非 结构网格过渡到外场的多层次直角坐标网格的粘性混合网格技术。这种混合网格方案一 方面发挥了三角形网格描述复杂外形的优越性，同时有一定压缩比的三棱柱网格较之全 非结构网格．可以更有效地模拟粘性边界层；另一方面外场用易于生成的直角坐标网格， 可以大幅减少全场的网格数量(一个六面体网格可以划分为5-6个四面体网格)，从而 降低对计算机资源的要求。本文的三棱柱网格是在物面三角形网格生成之后，用层推进 方法生成的．非结构网格则是由阵面推进法生成的。而多层次的直角坐标网格是由八叉 树方法生成的。网格生成和数值计算的实例表明该方法在二维复杂外形的粘性绕流计算 方面将有广阔的应用前景。 +固家自然科学基金资助项目 239 1 粘性混合网格生成方法 1 l 阐格生成的总体框架 本文的混合网格由层推进方法、阵面推进法和八分树方法生成。八分树方法用于生 成层移：结构的直角坐标网格：阵面推进法用于生成物面三角形网格和空间四面体网格： 层推进方法用于生成三棱柱网格。网格生成的总体框架如下： (1、I复杂外形的定义与描述。 ； (2、}建立矩形结构背景网格。 (3·表面三角形网格自动生成。 ’ (4)层推进生成有一定压缩比的三棱柱网格。 (5_l按全场最大的网格步长，用直角坐标网格覆盖全流场，并建立树状层次数据结构。 (61按背景网格步长控制参数的要求不断细分各直角坐标网格，直到网格分布达到所 需要的步长分布。 (7J删除物面附近三棱柱网格外的数层直角坐标网格(一般两三层即可)，构成包古 物体的“洞”，并将洞边界(由矩形组成)划分为三角形，由洞边界和物面网椿 组成封闭的初始“阵面”。 (8)阵面推进自动生成“洞”内的非结构网格．随即进行非结构网格的优化。 关于外形描述、表面三角形网格生成、矩形背景网格技术、八叉树方法生成多层次 直角坐标网格、空间阵面推进方法等，作者已在以往的：[作中做过详细的论述”．，j．或 者可以参阅文献(8一lO)。这里仅就层推进技术作较为详细的介绍。 1．2层推进技术 对于任意复杂外形，要得到台适的三棱柱网格，“层推进”必须满足以下要求：(1) 确定适当的推进步长以保证粘性网格层的光滑度，并能做到与外部的四面体网格光滑连 结；(2)确定适当的推进方向以保证“层推进”过程中网格不相交：(3)各网格