翼型—扰流片的分离气动特性计算

本文介绍了用涡面元法模拟带扰流片的翼型低速无粘分离绕流。在翼型和扰流片的面元上分布线性变化涡。在扰流片后的上下分离流线的面元上分布等强度的涡。上分离流线始自扰流片的梢部，下分离流线自翼型的后缘引出。分离所泡由两离散涡结尾。气泡内总压为常值，它与涡强大小一同求解。分离气泡的形状在迭代求解过程中确定。压强分布和升力系数的计算值与现存文献的数值结果和实验数据是一致的。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

正交流线坐标系中计算二维跨音速流的流线迭代法

本文提出了流线迭代法,用以实现文献[7]所提出的正交流线坐标系中的二维定常无粘气体动力学方程组的数值求解。 此方法的基本做法是:从物面起到流场远边界逐条流线地求出流动速度值,同时求出下一条流线的形状;按此方式对全流场进行迭代。迭代的收敛性很好。 用本文方法在零攻角和跨音速绕流的条件下计算了双圆弧翼型、NACA 0012翼型和R.A.E.101翼型。双圆弧翼型结果与实验数据符合较好。钝头翼型结果与其它方法结果相比是令人满意的。     

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

用粘流—无粘流相互作用方法计算翼型的分离和失速

本文用势流－边界层相互作用方法计算低速翼型的分离和失速。势流用对称面元法。边界层用改进的滞后掺混法，考虑了高阶项影响，适用于计算分离。文中对粘流－无粘流耦合方法作了改进，改进的半反－局部联立耦合方法。考虑了相邻点之间的作用，收敛性较好。计算了ＮＡＣＡ４４１２翼型在不同迎角下的压力分布和气动力。计算结果与实验符合良好，算例表明，翼型高升力状态计算必须包括尾流的作用，也应用对势流计算压力与实际压力之间     

用粘流－无粘流相互作用方法计算翼型的分离和失速

本文用势流－边界层相互作用方法计算低速翼型的分离和失速。势流用对称面元法。边界层用改进的滞后掺混法，考虑了高阶项影响，适用于计算分离。文中对粘流－无粘流耦合方法作了改进。改进的半反－局部联立耦合方法，考虑了相邻点之间的作用，收敛性较好。计算了ＮＡＣＡ４４１２翼型在不同迎角下的压力分布和气动力。计算结果与实验符合良好。算例表明，翼型高升力状态计算必须包括尾流的作用，也应当对势流计算压力与实际压力之间的差别进行修正。     

翼型动态失速的数值研究

用不可压缩流动的求解算法，结合WilcoxDC提出的k－ω模式和k－ωSST湍流模式，对翼型的动态失速进行了数值模拟。通过对典型的振荡翼型轻失速和深失速算例的计算结果分析可以看出：（1）绕动态失速翼型的流场结构十分复杂，轻失速和深失速在流动特性上有很大区别。计算结果显示：轻失速主要是由于后缘分离引起，分离涡的影响范围主要是在后缘附近。而深失速则首先形成很大的前缘分离涡，该分离涡在翼型表面上运动，并诱发出二次分离涡，引起翼型升、阻力系数的显著变化。（2）对于动态失速的翼型绕流，k－ωSST湍流模式是较为有效的，计算出的气动力系数迟带曲线变化趋势与实验结果符合得比较好。     

椭圆形控制环量翼型气动力的计算

本文用边界层理论与计及分离尾迹影响的位流理论相结合的方法计算控制环量翼型的绕流,用离散涡模型模拟分离尾迹的影响。由上、下翼面分离点处引入离散涡。只是翼型附近尾迹中的离散涡可以自由移动。离翼型较远的离散涡假设为按来流速度移动。忽略射流出口上游的边界层对射流的影响。这些简化使本文的方法十分有效。与实验比较表明,本文的方法能较准确地计算出翼型的分离点位置和分离点前的压强分布。计算出的升力系数与射流动量系数的关系与实验相符很好。     

低雷诺数翼型绕流的电磁控制

对低雷诺数不可压翼型绕流在Lorentz力(电磁力)作用下的流场变化情况进行了实验探讨.发现当施加的Lorentz力与翼型表面流向相同时,其表面流体的分离被抑制,从而使尾流涡消失,且翼型阻力减小而升力增加;当Lorentz力与流向相反时,则其受阻力变大而升力降低.     

二维不可压分离流计算的约束迭代法

本文提出约束迭代法以改进不可压分离流计算中的迭代收敛性。在每个时间步的迭代中,对涡量传输方程的边界涡量,预先进行近似估计。因此它相当于一种隐式的边值算法。 两个翼型的算例表明,约束迭代法具有高的收敛速度。由于无需寻找最优松弛参数,因而节省人力,便于使用。     

计及边界层汇流效应的多段翼型失速特性计算

 本文给出计及边界层汇流效应的多段翼型失速特性的解法。用高阶奇点分布面元法（Panel Method）求位流解,然后进行各翼段粘性尾迹形状迭代,并解出正常边界层和汇流边界层特性。当翼段上有后缘分离时,还要确定分离尾迹的形状,用位移厚度当量源（汇）模拟粘性效应。进行粘/位流迭代直至收敛。上述方法在超过多段翼型失速迎角时仍然有效。计算结果与实验数据比较,符合良好。     

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

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

二维弹性扑翼沉浮运动流动特性的数值研究

采用非定常势流板块法和Euler-Bernoul1i梁振动微分方程,对后缘带弹性薄板的翼型在作沉浮运动时翼型绕流和薄板弹性振动之间的流-固耦合过程进行数值模拟,并计算流动特性、非定常气动载荷,尤其是沉浮运动的推力效应.用有限差分法求解薄板的弹性变形运动,气动弹性分析采用弱耦合迭代方式.通过对翼型非定常气动力(升力和推力)、推进效率、尾迹脱落涡计算结果的分析发现,通过调整弹性薄板的刚度,可以使得翼型产生最佳的非定常推力.     

低速两段翼型的数值模拟

在多段翼型风洞实验中,很难观察翼型不同迎角下主翼、襟翼上的绕流及缝道流动。使用了O-H型混合的结构化网格,采用k-ε二方程湍流模型求解可压的N—S方程来模拟二元风洞中两段翼型的流场。计算结果与实验进行了比较,结果表明：计算与实验结果吻合良好,说明本方法可以较好地模拟两段翼型的绕流,并能很好地显示流场的变化情况。另外,从结果中可以看出：随着迎角的变化,主翼尾流和缝道间的流动对襟翼附面层分离有很大的影响作用。     

跨音速翼型分离绕流的粘性/无粘干扰迭代计算

本文用有粘/无粘干扰迭代的概念计算了跨音速任意翼型的绕流问题。位流的速位方程用AF2格式求解,而边界层微分方程用C-S盒式法求解,逆算法的引用可以克服边界层方程在分离点处的奇性问题,对分离区湍流代数模型的修正可以得到与实验更吻合的结果。计算结果表明有粘/无粘干扰迭代概念在小分离泡的情况下也是适用的。     

多段翼型失速分离特性的计算

本文用高阶奇点分布板块法(Panel Method》及边界层理论,对多段翼型的失速性能进行计算。采用分离尾迹的模型,按尾迹是流线用迭代法求出分离尾迹边界(此为内迭代)。在尾迹上配置奇点,然后把尾迹考虑在内进行粘/位流迭代(此为外迭代)。收敛后可得多段翼型空气动力特性,且可计算到失速以后的情况。算例有单段翼型GA(W)-1,NACA 4412,和F111飞机的三段翼剖面,所得结果与实验值能很好地符合。     

旋翼翼型-桨尖涡干扰的数值模拟

采用Kirchhoff方法计算NACA0012翼型的远场气动噪声,Kirchhoff面为包围整个翼型的一层网格。首先采用欧拉方程计算程序求解得到翼型的Kirchhoff面上在桨尖涡干扰下的近场气动流场解,再根据Kirchhoff方法求解远场声压。从而体现出在计算桨涡干扰噪声时,旋翼CFD技术与Kirchhoff方法相结合的良好效果。      

适用于整体求解吊舱推进器的定常面元法

为了方便开展吊舱推进器设计和混合式CRP(Contra-Rotating propeller)的面元法分析,提出了一种将吊舱和螺旋桨作为整体进行定常面元法求解的可行办法,对相关的面元法基本积分公式进行了分析,说明了定常面元法整体求解的结果与实际结果之间的关系。从吊舱桨尾流区的实际流场状况和计算效果出发,分析了整体求解吊舱推进器的面元法尾涡模型,并对此种尾涡模型的修正予以详细说明。在进行相关的理论分析后,编制了整体求解吊舱推进器的定常面元法程序,对吊舱推进器进行计算分析,并与试验结果和RANS(Reynolds-averaged Navier-Stokes)方法计算结果进行了对比,最后对本文所提出的定常整体面元法求解方法进行了误差分析。计算结果表明,整体面元法求解吊舱推进器的误差与试验值相比在5%以内;在计算耗时方面,将吊舱和螺旋桨分开的迭代求解面元法是整体求解面元法的1.7倍。     

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

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

基于响应面法的低速翼型气动优化设计

响应面方法相较于其它直接优化方法有其高效、实用的优势,此前的研究更多地将响应面方法用于超音速和跨音速翼型的减阻优化设计中。本文将此方法应用于低速翼型优化设计中,进行了基于RANS(Reynolds-Aver-aged Navier-Stokes)方程和自由转捩预测耦合求解的低速翼型气动优化设计。通过计算附面层方程得到附面层参数并用en方法计算转捩位置,并考虑了T-S波和层流分离造成的转捩。RANS方程计算中,使用了转捩过渡区模型,以保证附面层外边界压力分布的精度。RANS方程和转捩预测迭代进行至转捩位置收敛。在响应面模型计算中,使用不含二阶交叉项的二阶多项式模型,减少了构造模型所需的计算量;合理的选择设计空间保证了构造的响应面模型具有较高的精度。使用三个设计点的多目标优化设计,保证了设计的合理性。通过对NACA64(1)-112翼型优化计算结果表明,本文的方法可以有效地进行低速翼型的气动优化,各设计点上转捩位置也得到了改善,有较好的工程实用前景。