面搭接嵌套网格对变前掠翼的计算范文

《空气动力学学报》2014年第三期

1面搭接技术

面搭接要求不同的计算域之间有公共的边界面，即搭接面。在CATIAV5三维建模软件中，利用CATIA中图形的“隐藏／显示”功能，使同一界面在两侧分别显示，可保证交界面的准确无缝搭接。更重要的是，面搭接在物理意义上站得住脚的关键在于搭接面两侧流场在计算过程中通量的守恒传递。以二维情况为例，其原理参考图1［10］。图1中Block1和Block2分别表示搭接面两侧的网格单元集合，（i，j）和（m，n）分别标记两侧的单元格，（i＋1，j）为（i，j）单元格在Block2中的镜像的虚拟单元格。

2控制方程与离散

对飞行器的操纵面进行数值模拟时采用三维N－S方程。在笛卡儿坐标系（x1，x2，x3）中，定义速度分量（u1，u2，u3），采用求和约定惯例，无热源的三维N－S方程守恒形式为：湍流模型选择SSTk－ω模型，利用有限体积法将控制方程（4）离散，对流项选用二阶迎风差分格式离散［11］。物面为无滑移条件，远场为自由流条件，计算残差收敛精度为1×10－5。

3算例验证

1938年，LangleyMemorialAeronautical实验室的Wenzinger和Bamber［12］在试验段为7ft×10ft的低速风洞中，对带有开缝式襟翼的NACA23012翼型进行吹风试验，得到了该翼型不同襟翼开度下的气动参数。该试验段吹风速度约80mile／h，动压为16．37lbf／ft2，Re数约为2．19×106。为简明，仅以文献［12］中襟翼下偏30°的NACA23012翼型为例加以验证。计算模型如图2所示，其展长为弦长的1％。分别采用常规分块网格和嵌套网格方法剖分，其展长方向视图的网格如图3和图4所示。边界层内第一层网格高度为0．001mm，伸展比为1．1，y＋值控制在2附近［13］，整体网格数量均约为280万，两侧为对称边界条件。嵌套网格的优点在于不同舵偏时主体网格可多次使用，只需更改舵面偏转部分的网格即可，简单便捷。采用图3的常规分块网格和图4的嵌套网格，襟翼下偏30°时NACA23012翼型升阻特性的数值计算和试验结果如图5所示。如图5（a），采用嵌套网格和常规网格时的升力系数CL数值解基本相当；与试验值相比，二者能够较好地反应升力作用规律，线性良好，仅在失速迎角附近曲线趋势不明显，这主要在于数值算法对分离流模拟能力的差异，与是否采用嵌套网格无关。如图5（b），两种网格对阻力系数CD的计算精度基本无影响。二者在正负大迎角时的模拟值与试验值的差异较大，但基本趋势一致。因阻力系数与升力系数相比小一个数量级，其数值模拟也一直是一个难题［14］。图6和图7为两种网格状态迎角α＝0°时计算模型的马赫数云图和襟翼缝隙处的流线图。经比较，嵌套网格的马赫数云图与常规网格基本一致，嵌套网格的结果在搭接面处速度场过度均匀、平缓，效果较好。如图7（b），搭接面两侧的流线连续性较好，气流秩序流动；嵌套网格同常规网格一样，均能捕捉到主翼型下端有形成涡的趋势，襟翼上表面气流分离位置一致。综上可见，采用面搭接嵌套网格计算方法与常规网格计算结果几乎一致，而嵌套方法的网格处理却要简便得多。

4嵌套网格在变前掠翼开裂式方向舵性能计算中的应用

变前掠翼主要操纵面设计如图8所示。该布局为无尾式，在起降与巡航状态采用平直翼构型，能充分利用直机翼在亚音速阶段的高升阻比特性；在空中格斗时机翼前掠转换为前掠翼构型，能充分利用前掠翼在亚音速阶段缓失速、大迎角机动的优良飞行品质。变前掠翼开裂式方向舵的性能计算，需在不同的舵偏量下计算升阻特性，计算量较大；同时，其特征符合面搭接嵌套网格应用的要求。因而，采用上述计算方法处理。

4．1计算模型开裂式方向舵单侧舵面弦线与主机翼截面弦线的夹角为舵偏量δ，计算模型与舵偏量如图9所示，设计舵偏量δmax＝30°，操纵时可与升降副翼等联动。

4．2计算网格变前掠翼开裂式方向舵的计算网格如图10所示，半模网格数量约320万，边界层内第一层网格高度控制在0．001mm，以满足机体表面粘性边界层的计算要求［15］。操纵面偏转部分采用嵌套网格处理，并设置面搭接条件以保证搭接面通量守恒。

4．3计算结果图11和图12分别为平直翼在Ma＝0．8巡航时和前掠翼在Ma＝0．6格斗时变前掠翼随开裂式方向舵偏转的升阻特性曲线。由图11可见，平直翼巡航时升力系数随迎角增加而增大，但随阻力舵的舵偏量基本无变化；其最大升阻比约在4°迎角附近，此时升阻比随舵偏量的增加而降低，但在大迎角时升阻比随舵偏量变化不大。由图12，前掠翼的失速迎角接近40°，可见前掠翼的高机动优势明显，常规迎角下其升力系数变化类似平直翼，但大迎角下其值随舵偏量的增加而略有降低；前掠翼升阻比变化规律同平直翼，大迎角时升阻比更趋于稳定。平直翼在迎角α＝6°和前掠翼在α＝36°时的三维流线图如图13所示。由图13可见，开裂式方向舵打开，平直翼时气流经过操纵面后，在后缘形成稳定的涡流，对周围流场影响相对较小；前掠翼大迎角时开裂式方向舵表面气流产生分离，其流动由翼尖指向翼根，避免了翼尖提前失速引起的操纵效率急剧下降问题，这也可由大迎角时升阻比基本恒定反映出。

在变前掠翼开裂式方向舵数值计算中，机体外形结构较为复杂，若采用常规的分块网格剖分，不仅需对阻力舵偏转10°，20°和30°的三个状态进行全机网格剖分，而且操纵面处结构更加复杂，剖分更困难。而本文采用嵌套网格的处理方法，仅需对一个状态进行全机网格剖分，而在内空间的开裂式方向舵因拓扑结构简单，网格剖分的工作量甚小，其工作量仅为常规网格的1／3，计算的偏转状态越多，其优势越加明显。因此，本文所采用的嵌套网格剖分方法，可显著提高效率，并起到节约研究费用的目的。5结论主要结论如下：（1）面积加权法能够保证搭接面两侧的通量守恒传递。（2）嵌套网格剖分简单、便捷，且主体网格可多次使用，效率高。（3）基于面搭接的嵌套式计算方法可行，计算精度与常规网格相当，二者均能较好模拟升阻特性。（4）嵌套网格方法在变前掠翼开裂式方向舵计算中应用良好，前掠翼机翼气流流动由翼尖指向翼根，避免了翼尖提前失速引起的操纵效率急剧下降的问题，表明了开裂式方向舵在大迎角飞行中使用的有效性。基于面搭接原理嵌套网格计算方法的应用，不仅仅在于翼型等简单构型偏转的计算，更在于对复杂飞行器的各种操纵面偏转时性能的计算，显见其简单便捷，不失为一种高效的方法。

作者：王旭于冲苏新兵陈鹏单位：空军工程大学71003893286部队