二维翼型俯仰运动非定常运动的数值模拟方法研究.pdf

第九届(2005年)全国空气弹性学术交流会论文 二维翼型俯仰运动非定常运动 的数值模拟方法研究 许和勇王刚叶正寅 (西北工业大学翼型叶栅空气动力学国防科技重点实验室，陕西西安710072) 摘要：’本文研究了非结构对接运动网格下的翼型俯仰非定常运动的数值计算方法。在网格生成方面， 本文根据对接网格思想生成了一种简单高效的运动网格。在算法方面，本文采用双时间方法进行实时间推 进，运用通量线化假设和最大特征值方法进行雅可比矩阵分裂来实现隐式高斯一赛德尔拟时间迭代。本文 利用NACA0012翼型的两个俯仰运动状态作为算例，验证了本文提出的运动网格方案和引入的拟时间二维 隐式迭代算法对于模拟翼型俯仰运动具有一定的应用价值。 关键词：运动网格；对接网格；非定常计算；隐式迭代算法 引言 运用CFD技术模拟翼型的俯仰运动过程，一般是在运动网格的基础上求解非定常控制方 程。目前，生成用于非定常流场计算的运动网格主要有两种方法：刚性运动网格和弹簧动网 格。对于刚性运动网格，在计算网格运动速度时比较方便，但是外边界条件需要特殊处理， 比较复杂；而弹簧动网格在用于流场计算时，外边界网格速度为零，外边界条件可以等同于 定常求解器的外边界条件，但是在计算流场区域的网格速度时，要求事先计算出每个网格边 的倔强系数…，再通过求解与流场区域总节点数目同价的线性方程组得出新的节点坐标。本 文针对翼型俯仰运动的特点，在综合了上述两种运动网格的优点后，利用分区生成对接网格 的思路，生成了一种使用更加方便的非结构运动网格。这种运动网格在计算网格边速度时比 较方便和准确，而且外边界条件等同于定常流场求解器的外边界条件，也比较方便。 利用非结构网格进行流场求解的技术已经日益成熟，但是目前大部分基于非结构有限体 积法在时间推进上仍然是显式方法，其中经典的龙格一库塔时间推进应用最多。Batina旧1运 用通量线化假设和雅可比矩阵的正负特征值分裂的方法，构造了一种隐式高斯一赛德尔迭代 裂，提高了计算效率。本文采用中心有限体积法H1进行空间离散，采用全隐式的双时间方法 利用当地时间步长和隐式残值光顺等措施加速收敛，运用本文提出的运动对接网格方法，对 NACA0012翼型的俯仰运动进行了非定常无粘数值模拟。 ．．127．． 1．网格方案 于圆形区域中心。翼型绕25％弦长处做简谐振动，其瞬间攻角a(t)运动方程为： 口O)=口o+口。sin(o)t) 式中口。为平均攻角’口m为最大振幅攻角’减缩角频靴2爰一表示翼型弦长醒 角频率，圪为无穷远来流速度。 本文把计算区域分成两个区：以翼型的25％弦长处为圆心，以2为半径的圆域(挖去 翼型部分)为第一区，其余的为第二区。然后分别生成网格，图l为第一区的网格生成图， 图2为第二区的网格生成图。在计算中将两个区的网格进行对接，便形成了整体的计算网格， 如图4所示。在进行非定常计算时，第一区的网格速度为网格转动时的当地周向速度，而第 二区的网格速度除了与第一区网格对接的一圈网格外都为零。 非定常计算时，按如下方法对两个区的网格进行对接：(1)计算出第一区网格绕25％ 弦长处转动小角度Aa后的新的网格节点坐标，而网格编号顺序保持不变；(2)对于第二区 网格的每个对接边网格节点，从第一区网格的对接边网格节点中寻找到它的一个最近点作为 对接点，由于两个区的对接边节点都是等距的，所以只要在找到一个对接点后，其他的可以 按顺序直接确定对接点；(3)在找出所有的对接点后，把第二区网格的对接边节点全部用找 出的对接点替换，即形成了非定常计算中下一时刻的计算网格。 从上面的网格对接方法中可以看出，第一区网格只是位置发生了变化，而形状和面积并 没有发生任何变化，不需要重复计算，第二区网格除了对接边处的网格发生了小变形以外， 其他的网格没有发生任何变化。图5给出了两个非定常状态的网格重叠显示图，从中可以清 楚的看出网格的变形情况，注意的是变形区域的最外一层网格是原第二区网格的对接边界上 的网格。 图1第一区网格生成图(外圈为对接边) 图2第二区网格生成图(内圈为对接边) ．128． 图3第二区网格对接边处网格放大图 图4整体网格图