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

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

高升力翼型复杂流动模拟中湍流模型性能评估

采用软件Fluent中工程常用的7个涡黏湍流模型对某种高升力3段翼的降落阶段绕流进行了数值模拟.通过和试验及经过试验校验过的延迟脱体模拟（delayed detached eddy simulation,DDES）结果进行详细的对比分析,包括翼型压力系数、马赫数、涡量场和湍动能分布等,系统地研究了常用涡黏湍流模型对该高升力翼型的模拟性能.结果表明：对于平均流场,standard k-ω模型的模拟性能最好,能较好预测翼型压力系数、襟翼分离区位置和大小等;SST （shear stress transport）k-ω模型性能也较好,一方程SA （Spalart-Allmaras）模型和四方程v²-f模型具有相近的性能,而k-ε系列模型性能相对较差.对于湍流场的模拟性能,和上述平均流场对比所得结论相一致,但所有湍流模型均未预测出襟翼分离区附近湍动能“最大”的分布特征.     

湍流模型系数不确定度对翼型绕流模拟的影响

使用非嵌入式多项式混沌方法研究了湍流模型系数的不确定度对RAE2822跨声速翼型绕流模拟的影响。计算中关注了数值模拟的积分量（升力系数、阻力系数）和局部量（壁面压力、摩擦系数和空间马赫数分布）的不确定度量化结果。首先,从单输入变量入手,研究卡门常数的不确定度对数值模拟的影响。然后,同时考虑Spalart-Allmaras模型中9个参数的不确定度带来的影响。通过多项式混沌展开,得到系统输出对不确定输入变量的响应,由此可以得到输出的统计特性,包括平均值、方差和极值等信息。最后,在多变量不确定度量化过程中,通过Sobol指标来量化每个输入变量的不确定度对输出不确定度的贡献程度。本文计算只考虑了RAE2822跨声速翼型模拟的单一计算状态,影响规律是否可以推及其他工况和算例需要进一步检验。     

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

使用标准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%,进一步验证了其可信性,可将其进一步应用到低速飞行器的气动计算中.      

Gao-Yong模型用于吹气环量控制翼型的研究

应用Gao-Yong湍流模型模拟了吹气对翼型表面分离流控制的影响.结果表明,该模型不仅能够对翼型绕流的分离点、表面压力分布、升阻特性等做出较好预测,而且还能模拟出控制分离的吹气效果:①有效消除翼面分离涡;②"裹携"翼面来流进一步提高升力系数.提高吹气动量系数,升力曲线上移的幅度也相应增加;而吹气角度对吹气射流"裹携"作用的强弱也有一定的影响,当吹气方向相对于翼弦偏上时,吹气射流"裹携"来流与消除分离涡的作用增强,升力提升更为明显.      

基于浸入式边界方法的二维结冰机翼的数值模拟

应用浸入式边界方法模拟干净/结冰NACA0012翼形绕流问题.采用SIMPLEC方法离散求解二维不可压雷诺平均Navier-Stokes方程, 湍流模型采用一方程Spalart-Allmaras模型.计算结果与以往的实验结果吻合得很好, 证明了浸入式边界方法和所发展的程序的可靠性.      

多段翼型缝翼前缘结冰大迎角分离流动数值模拟

应用基于SST(Shear-Stress-Transport)湍流模型的IDDES(Improved Delayed Detached Eddy Simulation)方法,对大迎角状态下多段翼型缝翼前缘典型角状冰引起的复杂分离流动进行了数值模拟研究。采用后台阶流动标准算例和干净无冰多段翼型分离流动算例对数值方法的可靠性和适用性进行了验证。缝翼结冰状态下的数值模拟结果表明:来流迎角较大时,前缘角状冰将会导致结构相对稳定的流动分离泡产生,分离泡的非定常尾迹会对主翼前缘附近流场产生较为强烈的干扰,抑制了缝道流动的加速效应,使得缝翼增升效率降低。在失速点附近,由于分离泡回流强度随来流迎角而增长,同时脱落旋涡的输运方向逐渐向远离壁面方向偏移,使得尾迹影响区域范围和强度均有所增加。     

运用GAO-YONG湍流模型对扩压器内跨声速流动的数值模拟

运用GAO-YONG可压缩湍流方程组,采用同位网格SIMPLE算法,对扩压器跨声速流动中的二维激波/湍流边界层干扰现象进行了数值模拟。将计算得到的流场的时均参数与实验值进行比较,数值模拟结果在激波强度、壁面压力分布以及分离点和再附点位置等方面,与实验值吻合较好,表明GAO-YONG可压缩湍流方程组能够比较准确的模拟较强激波/湍流边界层干扰流动,从而进一步为GAO-YONG湍流模型的正确性及其在可压缩流场模拟方面的适用性提供了佐证。     

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

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

基于S-A湍流模型和间歇因子输运方程的转捩流数值模拟

考虑到γ-Re_θt转捩模型中间隙因子输运方程并不依赖于具体的湍流模型,因而提出了耦合S-A湍流模型和间歇因子输运方程的转捩流模拟方法.其主要思想是构造与间歇因子相关的两个耦合函数,分别作用到S-A模型的生成项和耗散项用以控制湍流的产生与发展,从而实现层流到湍流的数值转捩.结果表明:该方法中模型方程具有完全当地性的特点,易于实施,可直接用以求解三维转捩流动.4个典型算例表明该方法模拟结果与γ-Re_θt模型准确度相当,但少求解了两个输运方程,计算耗时减少了15%.      

壁面湍流模型对湍流分离流动数值模拟的影响

本文用三种近壁湍流模型计算了二个二维的湍流分离流动：单包流动和多包流动。结果表明：二层模型和低Reynolds数模型具有类似的特性。它们基本上给出合理的分离流动结构。壁函数的方法由于对数律的假定基本上不适合于计算湍流的分离流动。     

常气压辉光放电等离子体控制翼型失速的数值模拟研究

基于Shyy提出的大气压下均匀辉光放电等离子体与空气干扰的物理模型,通过求解电位势方程得到电场分布,得到了作用于流体上的电场力.通过数值求解考虑等离子体作用的流体运动控制方程,以NACA0015翼型低速绕流为对象,研究了常压辉光放电等离子体位置和个数控制翼型绕流分离的影响.位于分离点上游的等离子体能够有效地抑止翼型分离,而在分离区的等离子体对流动影响很弱,同实验结果一致,并给出了等离子体对翼型气动力影响的规律.     

基于IDDES方法的翼型结冰失速分离流动数值模拟

应用基于k-ωSST湍流模型的IDDES(Improved Delayed Detached Eddy Simulation)方法,就失速点附近翼型前缘典型双角状积冰导致的复杂分离流动进行了数值模拟研究.通过与风洞试验结果进行对比,表明对于此类分离流动问题,IDDES方法能够在壁面附近取得良好的速度预测结果,有效解析分离区域内的中小尺度湍流结构,较为准确地描述大尺度时均分离泡的再附位置和形态特征,适用于翼型结冰后复杂流动的精细分析.同时计算结果显示当此带冰翼型位于失速点附近时,角状冰后方脱落剪切层内部的旋涡不稳定析出和输运过程促进了外部流动与回流区域流动间的掺混,将导致流动发生非定常再附现象.     

翼型低雷诺数层流分离泡数值研究

翼型在不可压和低雷诺数条件下工作时的性能越来越引起了人们的研究兴趣。本文采用Rogers发展的双时间步Roe格式,求解拟压缩性修正不可压N-S方程。数值模拟了低雷诺数条件下(Re=6.0×104,1.0×105,2.0×106),E387翼型在不同攻角条件下(α=0°、4°、7°)翼上表面后缘部分的层流分离现象,对应于长分离泡。时均化的结果同McGhee[12]的试验结果进行比较,验证了计算结果的正确性。研究了低雷诺数层流分离的非定常和时均化特性,并对漩涡脱落中主涡、二次涡以及出现的漩涡对并的周期性过程和对气动力脉动造成的影响进行了较为细致的研究。非定常计算结果表明,低雷诺数条件下的层流分离现象,是周期性的旋涡脱落过程。所谓长层流分离泡是其时均化积分的结果。     

细长体大迎角湍流流场的数值模拟

通过数值方法对大迎角细长体超声速流场的模拟，探求正确模拟导弹大迎角绕流湍流流场的简单而有效的方法。数值模拟求解一般曲线坐标系下的三维可压缩Navier-Stokes方程，时间离散采用Euler向后差分，无粘项的空间离散采用二阶TVD格式，分别研究了Bald-win-Lomax代数模型及其修正形式(BLDS)对大迎角分离流动的模拟能力。数值试验表明修正的Baldwin-Lomax模型更精确地预测了流场的旋涡与分离情况，给出合理的表面压力分布，因而更准确地模拟了存在横向分离的导弹流场。     

超临界翼型跨声速激波振荡数值模拟

应用DES(detached eddy simulation)方法研究了SC(2)-0714超临界翼型跨声速自维持激波振荡现象.使用的DES方法预测了激波在翼型上表面的周期性运动,分析了雷诺数对激波运动的减缩频率、平均激波位置、脉动压力系数等非定常特性的影响.结果表明:雷诺数对激波振荡运动特性的影响呈非线性关系;而雷诺数为6×106时的激波振荡运动特性与雷诺数分别为15×106及30×106时的特性存在明显差异.      

壁面粗糙度湍流扩展模型及流动数值模拟

采用有限体积法数值求解雷诺平均Navier-Stokes流动控制方程,数值模拟表面粗糙度对气动特性和热力学特性的影响。通过修正Knopp的基于切应力输运k-ω湍流模型的粗糙度湍流扩展模型中边界条件ω的值,提高了壁面附近的涡黏度,增加了总壁面切应力,从而引入壁面粗糙度的影响。通过对不同粗糙度的平板流动和NACA 65_2215翼型绕流进行数值模拟,分析粗糙度对表面摩阻系数、斯坦顿数和升力系数的影响。计算结果表明:粗糙度使得平板表面摩阻系数和斯坦顿数增大,翼型流动分离提前,升力系数减小。流经粗糙壁面的流动处于过渡状态时,改进的粗糙度湍流扩展模型改善了原模型速度偏移量与经验公式相比偏小的缺点,同时,表面摩阻系数、斯坦顿数和升力系数与实验结果吻合得更好。该改进模型不需要对壁面附近网格进行加密,能够快速、准确地模拟流经粗糙壁面的流动。     

等离子体气动激励抑制翼型失速分离的仿真研究

通过求解表面放电的二维流体体力模型,建立了翼型等离子体流动控制的数学模型,得到等离子体气动激励诱导的体力和热量分布,与Navier-Stokes方程耦合求解.进行了低雷诺数条件下,等离子体气动激励抑制NACA0009翼型失速分离的数值仿真研究,研究了等离子体激励的强度、激励电极数目和激励位置对流动分离抑制和翼型升阻特性的影响.在雷诺数为58000、攻角为24°的情况下,施加等离子体激励后,升力系数由0.7449增大到1.2404;阻力系数由0.4012减小到0.3503.      

翼型和机翼阵风响应的数值模拟

采用有限体积法来数值求解非定常Euler和NS方程,并在非定常计算中通过引入"网格速度"来计入阵风的影响,实现了对刚性翼型和机翼阵风响应的有效数值模拟。首先采用该方法来对翼型攻角突然增大的阵风、陡边型阵风这两种典型阵风情况下的翼型气动力响应进行了计算,结果表明:计算结果与现有的一些理论结果都符合较好,与文献中计算结果也吻合,NS方程的计算结果与Euler方程计算结果在变化趋势上是一致的,但具体的数值有差别。最后针对攻角突然增大的阵风,对不同马赫数下不同展弦比平直机翼的阵风响应进行了数值模拟,得出的计算结果与文献中的计算结果一致。     

二维机翼混合相结冰数值模拟

近年来冰晶引发的飞机结冰问题逐渐引起人们的重视。针对冰晶和过冷水滴同时存在的混合相结冰问题,通过数值模拟的方法实现结冰冰形预测。空气流场和对流换热的计算中采用了转捩剪切应力输运（SST）湍流模型,基于欧拉法计算冰晶和水滴收集系数。在Messinger结冰热力学模型基础上进行扩展,分析了二维结冰部件表面在混合相气象条件下的传热传质过程,建立了适用于混合相的结冰热力学模型,同时考虑冰晶的黏附效应,添加黏附效率经验公式。利用FLUENT的用户自定义函数（UDF）编程求解混合相热力学模型,计算了霜冰和明冰条件下NACA0012翼型表面结冰情况,与国外风洞试验结果进行比较,验证了计算模型和方法的有效性。结果表明,冰晶黏附效应对混合相结冰量及冰形有很大影响,明冰条件结冰形状偏向楔形。