不可压缩翼型绕流数值研究

使用标准k-ω模型及与色散模型相耦合的k-ω模型分别计算了NACA0012翼型和NACA 4412翼型的低速绕流问题.NACA 0012翼型计算了其来流雷诺数为2.88×106,攻角从0°到15°范围内的流动结构、翼型表面压力分布和升力、阻力特性;NACA 4412翼型计算了临界雷诺数为1.52×106,攻角为13.87°时的流动分离和翼型表面压力系数 ,并与实验数据进行对比.结果表明:在同等条件下,使用与色散模型相耦合的k-ω模型计算得到的NACA 0012翼型的升力和阻力系数比标准k-ω模型提高精度约5%,NACA 4412翼型的表面压力系数精度提高了约3%,进一步验证了其可信性,可将其进一步应用到低速飞行器的气动计算中.      

外边界位置对数值计算翼型绕流的影响

用数值方法计算翼型绕流时,外边界位置一般只取在距翼型有限的距离。本文用不可压位流理论找出该问题的精确解,从而得到外边界位置对计算的升力、阻力和表面压力分布的影响。     

M_∞=1及其附近翼型绕流的计算结果

本文全速势方程的有限差分数值计算结果,给出了0°和2°攻角下NACA0012翼型当地马赫数分布冻结时的自由流马赫数范围,以及M_∞稍大于1变至1.30时前方脱体激波的变化位置。     

NACA4412翼型低速绕流数值计算中湍流模型对比

使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。     

低雷诺数下二维翼型绕流的流场特性分析

采用高精度有限差分格式,对低雷诺数下二维翼型绕流进行了直接数值模拟,计算了雷诺数为1.0×104,NACA0012翼型0°,4°以及10°攻角下的流场。计算结果表明:在0°和4°攻角条件下,翼型绕流尾迹区的统计特性相似,0°攻角下的统计量值具有很好的对称性;在距翼型尾缘0.3弦长以后的尾迹区,旋涡排列成类似涡街的结构,涡量的极值、压力的极小值都位于旋涡中心,沿着流向,涡量极值的绝对值逐渐减小,压力的极小值逐渐增大。10°攻角下,翼型上表面从前缘开始分离,尾迹区统计分析结果所得图象与0°和4°完全不同,数值上较后者结果大;在翼型尾缘处,涡量的卷吸,高压力区域的形成,是旋涡脱落的条件,正向和反向旋涡的交替脱落,形成了类似涡街的结构。      

电磁力控制翼型体绕流的数值研究

为考察包覆电磁激活板的翼型体绕流特性,采用数值模拟的方法研究了低雷诺数下电磁力对翼型体流场结构和绕流过程的影响,探讨了不同迎角和不同电磁力大小对翼型体尾流涡街形态和流体边界层结构的影响。计算结果表明,作用于翼型体表面的电磁力可以明显改变流体边界层的结构,抑制边界层分离,并有效提高翼型体的升力。     

绕翼型跨音速大扰动对称势流的插值混合差分法

 本文将跨音速定常小扰动势流混合差分法推广到跨音速大扰动定常势流,提出了在局部速度坐标系中求解跨音速精确势流方程的插值混合差分法。作为算例,计算了双圆弧翼型和NACA0015翼型对称问题压力分布,并与已知实验值和双圆弧翼型小扰动混合差分法计算值进行比较,结果接近。试算表明,本文提出的插值混合差分格式是稳定和收敛的。本文解决了M_∞趋近于1的计算难点。     

不可压缩实际粘性流体绕流孤立翼型的数值研究

本文提出了基于微分-积分方程组求解N-S方程的有限差分法求解不可压缩实际粘性流体绕孤立翼型流动。计算中采用了K-ε湍流模型结合考虑逆压梯度的壁面函数法,并导出任意非正交曲线坐标系下的压力修正方法和适用于大流通量的混合格式。分别计算了层流和湍流情况下带分离翼型绕流,与流动显示及实验数据比较表明,在预测分离点位置,翼型表面压力、平均速度分布、雷诺切应力uv方面,两者吻合,结果是满意的。     

在正交流线坐标系中数值计算二维升力翼型跨音速绕流

本文将文献[1]提出的流线迭代法进一步应用于计算正交流线坐标系中有攻角时的二维跨音速绕流。用该算法计算了α=1°情况下的6％双圆弧翼型和NACA0012翼型,计算结果与实验数据符合较好。     

充气式翼型绕流大涡模拟分析

针对多梁式充气机翼阿翼型,采用大涡模拟方法对其表面的复杂流动进行了数值仿真研究,得到了其迎角在0°-20°时的流动情况和气动力特性,并与其他仿真结果和风洞试验结果进行了比较分析。结果表明：采用大涡模拟方法得到的计算结果与风洞试验结果较为吻合,可以更好地模拟充气式翼型的表面流动情况。     

不可压缩势流绕多段翼型流动的快速迭代解法

利用保角变换计算不可压缩势流绕翼型的流动,具有计算速度快和要求计算机内存小的优点,但主要用于单段翼型问题,而对于多段翼型,尤其是三段和三段以上的翼型则很难应用。本文在单段翼型的保角变换的基础上,提出了一种利用迭代算法计算任意多段翼型的不可压缩势流的绕流问题的方法。本文方法和目前常用的边界元法相比,具有计算速度快的优点。本方法还可进一步推广到考虑无粘流/边界层的迭代计算中。     

翼型跨音速绕流的数值计算

本文选用精确速势方程作为翼型跨音速无粘绕流的数学模型。在变换过的直角座标中,用有限差分(亚音速区用中心差分,超音速区用旋转差分)对精确速势方程离散化。差分方程形成的代数方程组用列松弛迭代法求解。 为了计入粘性效应,利用精确速势方程和附面层动量积分方程进行联合迭代求解。这一点对超临界翼型的计算显得特别重要。 算例与其它数值解及实验作了比较,对于一般翼型和超临界翼型,本文结果是良好的。     

NLR-7301高升力翼型绕流评估计算

探讨了在二维高升力多段翼型绕流计算中附面层网格采用各向异性三角形网格的技术。通过对经典的NLR-7301两段翼型绕流的计算,从数值计算角度评估了自研混合网格雷诺平均N-S方程求解软件对该类型附面层网格的解算能力。通过与确认试验数据对比,在附面层采用各向异性三角形网格进行二维高升力多段翼型绕流计算,其结果具有较高的可信度。在此基础上,采用同样策略,评估了多段翼型绕流计算中翼型后缘的几何处理、湍流模型选取等关键CFD要素对计算结果的影响。     

跨音速非定常二维翼型绕流计算

引言 为了计算跨音速非定常气动载荷和了解其流场的性质,需要对翼型绕流进行数值模拟。在构造一个跨音速非定常流动的计算方法时,重要的是要考虑能够处理流场中不断运动的激波;其次要能够处理各种类型的翼型运动,并具有高的计算效率。     

可压缩性粘流绕翼型的摩阻计算

本文研究了一种可压缩性粘流绕翼型的摩擦阻力计算方法。这种方法包括用修正了的Thwaites~[1]法计算头部层流边界层以及按Green~[4]法计算湍流边界层。按照上述方法计算了在不同雷诺数和不同马赫数下两个翼型剖面的摩阻系数。经与文献[5]比较,吻合较好。     

欧拉方程多段翼型绕流计算的嵌套网格方法

本文采用嵌合体埋入技术，将主段翼型和襟翼分成两子域，对多段翼型的复杂构型进行了分区迭代求解，在有限体体积离散的基础上，求解欧拉方程计算各子域流场的解，通过各流场的耦合迭代到多段翼型绕流的解，典型的带３０％襟翼ＧＡ（Ｗ）－１翼型算例表明，该计算方法稳定，能较地模拟多段翼型中主段翼型与襟翼交接处的复杂流动现象计算结果与实验值有较好的一致性。     

多段翼型粘性绕流计算研究

针对多段翼型外形结构复杂,结构网格生成难度大的缺点,运用割补法对重叠网格技术进行了研究,生成贴体及与边界正交的高质量的结构化网格,在识别贡献单元方面,介绍了一种基于kd树的找重方法.流场求解应用雷诺平均Navier-Stokes方程结合Spalart-Allmaras一方程湍流模型在所生成的重叠网格上对多段翼型的粘性绕流进行了模拟,为了兼顾流场计算效率问题,在重叠网格的各个子块上采用多重网格方法.计算结果表明,用该方法可以很好地预计压力分布,并且计算效率得到很大的提高.     

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

机翼跨声速绕流的Galerkin有限元解法

本文给出了翼型和机翼跨声速绕流的Galerkin带权余数有限元解法。利用一系列适当形状的元素将计算区域离散,在元素内采用线性插值函数。对于超临界流动,文中采用了适合于有限元法的超声速上风技术,具有捕捉激波的能力。将这种技术应用于线松弛迭代解法中,在计算钝头翼型和机翼绕流时,都获得了成功。     

多段翼型绕流的数值模拟研究

采用C—H型混合的结构化网格,选用Wilcox提出的k-ω模型求解可压的N-S方程,模拟了多段翼型的流场。结果表明：在4°和6°迎角下的非定常计算结果与试验结果的对比,计算结果与试验值吻合的较好。本文应用了亚迭代时间法,正确地描述了30P30N多段翼型前缘结构区域的流动特征和流动机理。