覆冰导线气动力特性风洞试验研究

针对新月形和D形两种典型断面的特(超)高压大截面厚覆冰导线进行了气动力系数风洞测试,获取了0°～180°以5°攻角为间隔各攻角下的气动三分力系数,用均匀流和均匀湍流两种风场研究湍流对覆冰导线气动力特性的影响,根据Den Hartog理论分别计算了两种断面覆冰导线可能发生舞动的风攻角范围。结合以往试验结果,归纳了两类覆冰冰型气动力特征。试验结果表明:新月形覆冰导线可能发生Den Hartog舞动的攻角范围为15°～30°和175°～180°,D形覆冰导线则为65°～85°及135°～150°,湍流对导线舞动气动稳定性存在影响。与以往试验结果比较表明:本试验结果能够反映大覆冰新月形和D形覆冰导线气动力特征,小覆冰新月形导线气动力特性则与本试验存在较大差异。     

冲击射流流动换热超大涡模拟研究

为了准确预测发动机热端部件中广泛采用的冲击射流冷却复杂的流动和换热特性,发展了基于BSL k-ω模型的超大涡模拟(VLES)高精度模拟方法,并对高雷诺数Re=4×104,两种不同射流距离2和6的单孔冲击射流及三孔冲击射流这一经典的流动传热问题进行三维非稳态高精度数值计算。同时,将分离涡方法 (DDES)和k-ωSST,RNG,Transition SST三种RANS方法的数值模拟和开发的超大涡模拟(VLES)方法进行对比。研究表明,VLES方法均能够准确捕捉冲击射流流场的复杂非稳态流动及传热特征,包括自由射流区、壁面射流区小尺度涡系和大尺度湍流结构的演化和破碎,同时冲击壁面的换热系数计算结果与实验值吻合较好。DDES方法未能准确捕捉流场复杂的小尺度湍流结构,壁面换热计算结果与实验值差异较大。RANS方法计算的换热结果与实验数据差异最大,基本未能预测到壁面换热特性。在相同的计算网格和计算方法下,VLES方法计算结果优于DDES方法,DDES方法一般好于RANS方法。这表明新开发的VLES方法能够准确地计算冲击射流相关的流动及换热问题。     

覆冰四分裂导线节段模型驰振风洞模拟试验

设计制作覆冰四分裂导线节段模型驰振风洞模拟试验系统,利用位移测量装置和数据采集系统记录导线节段模型驰振过程中垂直和水平位移的时间历程以及运动轨迹。试验得到典型厚度新月形和扇形覆冰四分裂导线节段模型在典型风速下的驰振响应。进一步建立覆冰导线节段模型的有限元模型,利用已获得的舞动数值模拟方法模拟其驰振过程,得到与试验一致的结果,从而验证了舞动数值模拟方法的正确性。试验结果为覆冰导线舞动研究提供了重要的参考数据。     

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

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

转捩对风力机翼型和叶片失速特性影响的数值模拟

采用基于k-ωSST湍流模型的Gamma-Theta转捩模型对S809翼型和NREL Phase VI叶片进行了考虑转捩的气动力数值模拟,研究了转捩对其气动特性特别是失速特性的影响。首先对S809翼型在迎角0°～30°范围内开展了数值模拟,比较了转捩模拟和全湍流模拟获得的翼型升力系数和流场特征,发现前缘层流分离泡的存在显著影响了翼型的失速特性;然后对NREL Phase VI叶片开展了类似的数值模拟,结果显示转捩对叶片失速特性和翼型失速特性产生影响的作用方式是相似的。经过对以上数值模拟结果加以分析认为,转捩对翼型和叶片失速特性的影响主要通过前缘层流分离泡的作用体现出来,前缘层流分离泡的存在使翼型和叶片更早地进入深失速。     

覆冰导线群的动态气动力特性

针对我国电力线路严重驰振区的实际情况进行了覆冰导线群的动态气动力特性测试。试验模拟驰振发生时的扭转振荡,比较两种典型的覆冰导线模型的静态和动态的气动力特性。讨论了准静态假设在覆冰导线群的驰振分析中的适用性。覆冰导线模型的静态气动力特性测试可以给出动态气动力的主要特征。但是在覆冰导线的驰振分析中确定驰振发生的最低风速等条件时,需要考虑动态气动力负阻尼与静态试验结果的差别。     

NASA TrapWing高升力标模数值模拟研究

针对NASA TrapWing高升力全展襟翼构型，采用计算流体力学（ CFD）方法进行三维复杂流场仿真模拟，考察网格尺度和湍流模型对高升力模型气动特性的影响。采用“超立方体”概念，生成绕TrapWing模型的不同网格密度高质量多块结构网格，通过求解雷诺平均Navier-Stokes方程，研究网格尺度对高升力模型气动特性的影响。在此基础上选取中等规模计算网格，考察Spalart Allmaras和Menter k-ωSST湍流模型对高升力全展构型流场模拟能力，分析湍流模型对气动特性的影响。研究结果表明：网格策略具有较好的网格收敛性；湍流模型对翼稍附近位置上翼面压力系数的预测稍有影响，SA湍流模型预测的压力系数较SST更接近实验值。确认研究工作为大型飞机增升装置数值模拟提供了一定的参考。     

基于非结构/混合网格的脱体涡模拟算法

为了提高二阶精度有限体积算法的湍流数值模拟能力,在原始Roe格式基础上建立了与脱体涡模拟(DES)方法相匹配的二阶混合耗散自适应格式,根据流场信息自动调节格式耗散,并分别基于Spalart-Allmaras一方程湍流模型和k-ω剪切应力输运(SST)两方程湍流模型,发展了基于非结构/混合网格的DES方法。采用该方法计算了雷诺数为3 900的圆柱绕流和NACA 0021翼型60°大迎角分离流两个典型算例,通过与试验数据以及其他数值结果的对比验证了该方法的可行性。同时开展了不同数值格式、湍流模型的对比分析,研究结果表明:采用混合格式的DES算法能够解析更小尺度的湍流涡结构、计算数据更接近试验值;本文的DES类算法受其基准湍流模型影响较小。     

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

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

Sajben扩压器复杂流动中湍流模型的性能评估

用8个常用的湍流模型对Sajben扩压器中跨声速流动进行了数值模拟，评估了Spalart-Allmaras， 标准k-ε， RNG (re-normalization group) k-ε，realizable k-ε，标准k-ω，SST(shear stress transport) k-ω，v2-f，Reynolds stress共8个湍流模型对激波/湍流边界层相互作用的模拟预测能力.通过与实验数据比较发现:SST k-ω模型和v2-f模型比其他模型模拟的更准确，其中SST k-ω模型比v2-f更能准确地预测壁面压力，然而对于分离点、再附点以及分离区长度v2-f比SST k-ω预测得更准确.      

覆冰四分裂导线静态气动力特性试验

覆冰四分裂导线的空气动力特性,由于覆冰子导线尾流的相互干扰可能不同于覆冰单导线。针对两种不同冰型制作覆冰四分裂导线模型,通过风洞试验测试获得不同冰厚和不同风速下覆冰四分裂导线静态空气动力系数随攻角的变化曲线。结果表明,在一定的攻角下,子导线尾流相互干扰对空气动力系数有明显影响。所得试验结果为覆冰四分裂导线的舞动及其防止技术的研究提供必要的数据。     

临近空间螺旋桨低雷诺数高效翼型数值分析

针对临近空间翼型边界层变厚更容易发生分离的特点,在修正了湍流模型后,通过分析比较13种低雷诺数翼型的升阻系数得到了比较适合用于临近空间螺旋桨叶素的高效翼型。结果显示,攻角在-6°～8°、8°～14°、14°～22°三种工况下,分别采用SST k-ω、RNG k-ε、Realizable k-ε湍流模型,可以得到与国外实验比较接近的合理结果;工况条件对翼型气动性能好坏的影响很大,翼型S-1223和FX63-137在所研究的工况内都具有较好的气动性能;因此可以选择这两种翼型作为临近空间螺旋桨用高效翼型。     

冲压翼伞流场与气动操纵特性的数值模拟

采用有限体积法求解shear stress transport(SST) k-ω二方程湍流模型下的Navier-Stokes(N-S)控制方程,对冲压翼伞的气动特性进行数值模拟,分析翼伞的流场机理和气动操纵特性.模拟得到的升阻特性与试验数据较吻合,在此基础上分析前缘切口、弧形下反以及稳定幅对升阻特性的影响.通过软件Fluent实现了非定常流动模拟,研究了翼伞的非定常升力特性,其升力系数的脉动受脱体涡的非定常过程影响,当迎角为16°时,翼伞升力变化周期为0.36s.最后分析了翼伞稳定滑翔阶段的纵向静稳定性,相比于单边后缘下拉方式,通过闭合翼伞一侧进气口实现航向操纵更稳定有效.      

降雨条件下翼型的气动特性

采用计算流体力学方法对降雨条件下翼型的气动特性进行了研究.通过求解定常可压的Navier-Stokes方程来计算流场,采用拉格朗日法对流场中的雨滴轨迹进行了模拟跟踪,得到了翼型表面的雨滴收集率;建立了水膜层数学模型,并假定水膜层沿着翼型表面的法向方向增长,预测了降雨引起的翼型外形的变化,并得到了水膜层表面的粗糙度高度;采用k-ωSST(shear-stress transport)两方程湍流模型结合增强的壁面函数,研究了降雨对翼型气动特性的影响.结果表明,在低迎角范围内,降雨对翼型的升力系数和阻力系数影响很小;当达到失速迎角后,降雨会引起上翼面边界层气流的提前分离,造成翼型气动性能的严重损失.      

一种分布式电动飞机螺旋桨滑流影响机理

提出了一种分布式电推进螺旋桨飞机，采用二阶精度求解RANS方程的k-ω SST （Shear-Stress-Transport）湍流模型，基于多参考系（MRF）方法，针对低速特性进行数值模拟，得到了分布式螺旋桨滑流效应对全机气动特性的影响规律，重点对螺旋桨后方速度场及机翼表面压力分布进行分析。结果表明有滑流状态增加了全机升力和阻力，升力系数最大增量超过65%，且升力增量随迎角的增加而增大，改善了失速性能，增加了低头力矩；螺旋桨旋转增加了周向速度，改变了径向速度分布，增加轴向速度，高能量螺旋桨滑流改变了机翼当地升阻力特性；螺旋桨桨叶向上旋转一侧气流受上洗影响而局部迎角增加，另一侧局部迎角降低，越靠近桨盘位置，受螺旋桨洗流带来的影响越大；螺旋桨的旋转方向对螺旋桨两侧的机翼表面压力分布有较大影响，尤其是翼尖螺旋桨对全机气动性能影响较大。     

平流层飞艇大迎角气动特性分析

为了获得大迎角下的平流层飞艇的气动特性,基于N-S方程及k-ωSST湍流模型,对某平流层飞艇的大迎角气动特性进行数值模拟分析.结果表明：平流层飞艇的大迎角失速特性及力矩特性与飞机具有较大的区别,飞艇的失速迎角较大（45°左右）;俯仰力矩特性在失速前呈现完全不同的两种趋势,小迎角范围内其静不稳定,迎角增加到一定程度后变为静稳定且维持该状态直至失速.研究结果对平流层飞艇的大迎角操稳设计、飞控设计具有重要的参考意义.     

Pω增强型k-ω湍流模型在三角翼旋涡流动的应用

在三角翼旋涡绕流数值模拟中,标准 Wilcox k-ω湍流模型生成项未考虑旋度的影响而导致预测的旋涡强度较弱。通过引入探测因子区分剪切层和涡核,在旋涡流动的高旋度区域增加ω方程生成项的方法,基于结构化网格上的 RANS 求解器,加入了 Pω增强型 k-ω湍流模型,对绕尖前缘三角翼亚声速和跨声速旋涡流场进行了数值模拟。计算结果与 NASA 的 NTF 风洞和 DLR 的 DNW-TWG 风洞试验数据进行了对比分析,结果表明：不论在亚声速还是跨声速自由来流条件下,Pω增强型 k-ω湍流模型计算的压力分布、涡破裂位置均与试验数据吻合良好,准确地预测出了三角翼上翼面的主涡、二次涡结构,特别是跨声速条件下激波干扰导致的涡破裂的临界迎角及涡破裂位置,表明 Pω增强型 k-ω湍流模型在绕三角翼旋涡流动数值模拟中具有良好的适用性。     

临界雷诺数区准椭圆形覆冰导线风压特性研究

采用刚性模型进行测压试验,得到了不同雷诺数下准椭圆形覆冰导线的风压分布规律,通过对比平均风荷载、脉动风荷载及风荷载谱等参数,分析了雷诺数对风荷载以及横风向驰振稳定性的影响。当雷诺数达到临界区,与亚临界区的对应值相比,平均阻力系数下降、平均升力系数随风向角变化幅度大且在某些对称工况产生横风向平均升力系数;平均风压系数分布对风向角等参数更为敏感。旋涡脱落由亚临界区的规则脱落变为不规则脱落,周向风压相关性减弱,特征频率消失。临界区内平均升力系数急剧的下降段使得结构更易发生横风向驰振。     

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

使用标准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湍流方程组计算翼型分离流

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