数值涡流发生器技术在平板绕流数值模拟中的应用研究

本文将数值源项模型应用于叶片式涡流发生器的数值模拟中,以叶片在当地流场中产生的升力作为源项,该源项由叶片的坐标位置、安装方向、高度等参数根据Plandd升力理论确定。数值模拟中,叶片位于平板上,其高度约为当地边界层厚度的1／6。在相同计算条件下,对比分析采用该源项模型数值模拟结果和相关资料的对应数据,证明采用本文源项模型可以得到与资料采用全网格以及采用其涡流发生器模型基本一致的诱导速度场。     

CRM-WBPN风洞试验模型数值模拟

基于雷诺平均Navier-Stokes（RANS）方程和多块对接结构网格技术，采用亚跨超CFD软件平台（TRIP 3.0），研究了共同研究模型（CRM）风洞模型支撑装置对CRM翼/身/架/舱（CRM-WBPN）模型压力分布及气动特性的影响。构造了一个包含小、粗、中、细网格的4套网格，并进行了网格收敛性研究，在升力系数为0.50条件下，开展了固定升力系数下模型支撑对CRM-WBPN模型和CRM-WB模型压力分布及挂架短舱阻力增量的影响，以及不同来流迎角下模型支撑对CRM-WBPN模型气动特性的影响。通过与NASA国家跨声速风洞设备（NTF）试验结果和CRM-WB模型数值模拟结果的对比分析，表明模型支撑对固定升力系数下挂架短舱阻力增量影响较小；在计算构型中考虑模型支撑的影响，引起升力系数和阻力系数下降、俯仰力矩系数增加，数值模拟结果更加接近相应的试验结果。     

基于动量源方法的直升机旋翼/机身流场数值模拟

建立了一套基于动量源方法的直升机旋翼/机身流场数值模拟方法.将旋翼对流场的作用以动量源项的形式代表,建立了一个包含动量源项的CFD方法.针对直升机旋翼/机身流场的特点以及动量源方法对网格系统的要求,采用了一种混合网格生成方法.通过结合旋翼桨叶的运动方式、几何特征及气动特性,建立一个包含动量源项的N-S方程的旋翼流场计算方法和迭代流程,并采用该方法进行直升机旋翼、机身和旋翼/机身下洗流场的数值模拟,得到了关于旋翼/机身干扰流场的一些有意义的结论.     

结冰模拟软件AERO-ICE中的关键数值方法

自然积冰对于航空飞行安全造成重大隐患,飞机在穿越富含过冷水汽的云层时冰形将按一定的物理规律积聚生长。介绍了一款三维结冰数值模拟软件AERO-ICE,该软件由网格自动生成、空气流场RANS计算、水滴场欧拉方法计算、结冰热力学分析四个模块组成。在空气流场计算方面,采用SPF k-v²-ω湍流模型,该模型引入湍流非平衡特性修正,预测的带冰翼型最大升力系数和失速攻角相对SA和SST模型有显著的提高。水滴场欧拉方程由于源项较大,迭代求解时容易发散,AERO-ICE软件采用流场光顺、二阶MUSCL空间离散以及LU-SGS隐式时间推进方法改善了数值稳定性。在结冰热力学分析模块,AERO-ICE软件同时具有Messinger与Myers模型,并将Messinger模型预估的壁面温度作为Myers模型的边界条件,从而解决了Myers模型温度设置的经验性问题。AERO-ICE软件支持多块网格、多重网格加速技术与大规模并行计算,其冰形计算结果得到了初步的验证。     

微型扑翼飞行器非定常运动对平尾的影响

以西北工业大学自行研制的微型扑翼飞行器ASN211为研究对象,利用其简化的二维扑翼及平尾串列翼模型进行了非定常数值模拟,分析了扑翼俯仰运动及沉浮运动对平尾气动性能的影响。在数值模拟模块中,模型的俯仰运动及沉浮运动由动网格技术实现。通过计算流体力学(CFD)软件Fluent对此非定常流场进行数值计算,重点研究了扑翼非定常运动尾流对平尾气动效率的影响。定常状态与非定常时均条件下平尾升力曲线的对比分析表明,扑翼的非定常运动能够增大平尾的失速迎角及最大升力系数,因而使平尾的失速特性得到改善。     

CRM翼身组合体模型高阶精度数值模拟

基于五阶空间离散精度的WCNS格式,开展了CRM翼身组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细、极细四套网格,网格规模从"粗网格"的2 578 687个网格点逐渐扩展到"极细网格"的65 464 511个网格点。研究了设计升力系数下,网格规模对气动特性、压力分布和翼根后缘局部分离区的影响,采用"中等网格"开展了抖振特性的数值模拟研究。通过与二阶精度的计算结果、DPW V统计结果和部分试验结果的对比分析,高阶精度数值模拟结果表明,阻力系数计算结果与DPW V统计平均结果吻合较好;网格密度对机翼上表面的激波位置和翼身结合部后缘局部分离区略有影响;迎角为4°时,升力系数下降的主要原因是机翼上表面激波诱导分离区和翼身结合部后缘局部分离区的增加。     

结构网格方法对高升力构型的应用研究

采用"亚跨超CFD软件平台"(TRIP)数值模拟了梯形翼全展长与半展长高升力构型的复杂流场,并与试验做了对比分析。对应风洞试验是在NASA Langley 14×22英尺亚声速风洞(FST)和NASA Ames 12英尺增压风洞(PWT)中完成的。计算中采用一方程SA湍流模型和MUSCL-ROE格式,并综合运用对接/重叠/拼接网格技术,数值模拟了两种高升力构型的气动特性,给出了典型站位的压力分布,并对比研究了不同结构网格技术对此类高升力构型的计算差异。研究表明,与修正后的试验数据相比,数值模拟得到的气动力系数和典型剖面的压力分布均与试验结果吻合良好。     

倾转旋翼机旋翼对机翼气动干扰的建模及分析

建立了一个耦合旋翼自由尾迹模型和机翼面元模型的悬停状态旋翼/机翼气动干扰迭代计算方法,用于分析旋翼对机翼的气动干扰影响.在该方法中,为较好地模拟大桨盘载荷及大扭转桨叶的气动特征,旋翼桨叶采用Weissinger-L升力面模型;为考虑厚度效应及机翼的升力影响,建立了包含源面元和偶极子面元的厚度机翼模型;为正确模拟旋翼桨尖涡与机翼表面间的贴近干扰,采用了一个"分析数值匹配法"的"贴近涡/面干扰模型".应用上述方法,对单独倾转旋翼下洗流分布以及旋翼对机翼的气动干扰影响进行了计算.结果表明,在旋翼下洗流场的干扰下,机翼各剖面都产生向下载荷,但并非简单地随拉力系数的增大而增大;机翼受到的旋翼干扰影响与旋翼下洗流沿桨叶展向变化密切相关.     

旋翼近地飞行时诱导速度的N—S方程计算

采用有限体积法用三维Ｎ－Ｓ方程对旋翼近地飞行时的速度场进行了数值计算，对旋翼桨盘处的诱导速度进行了分析。旋翼对流场的作用由分布在特定区域内的动量源项模拟，动量源项的强度由旋翼的几何尺寸，叶片剖面的空气动力特性及当地流场的性质迭代定出。本文讨论了旋翼近地飞行时，在不同前进比及不同离地高度情况下的诱导速度随径向距离及方位角的变化，并将数值计算结果与实验结果进行了比较。     

应用等离子体进行翼型升力控制的数值模拟研究

将等离子体对中性气体的作用模型化为彻体力矢量,求解带源项的Navier-Stokes方程,数值模拟了在NLF(1)-0213翼型上表面60%弦长处安装等离子体激励器对升力的控制效果.彻体力为净电荷在外加电场作用下的电场力.解拉普拉斯方程得到外加电场分布,等离子体中的净电荷分布由泊松方程给出.升力线计算结果与实验值吻合,激励器工作时,升力线向上平移,控制效果与襟翼类似.     

基于DOE方法超紧凑S型进气道被动流动控制优化设计

采用实验设计（DOE）方法对影响超紧凑进气道流动控制效果的五个因素进行优化设计。采用数值源项模型代替真实物理叶片进行流场求解，分析了27组涡流发生器应用于进气道流动控制的效果，从中找出各影响因素变量对进气道总压恢复、流场畸变不同影响，应用响应曲面法给出最佳影响因素组合，为超紧凑进气道被动流动控制设计提供技术参考。     

涡流发生器数值计算方法研究

利用RANS方程和SA湍流模型解算器,采用多重网格法和预处理技术,对叶片式涡流发生器进行了数值模拟研究,初步探讨了涡流发生器的安装方式、剖面形状、几何尺度等因素在机翼分离流动控制中的影响规律和设计原则,并初步研究了安装涡流发生器对超临界机翼气动性能的影响。     

一组近耦合鸭式布局的低速气动力数值模拟

本文使用位流中前缘离体涡模拟的数值计算方法,对于不可压流动,大迎角情况下的气流流经一组近耦合鸭式布局的流动,进行了数值模拟。分析表明,在大迎角下,在一定的主翼-鸭翼的参数选择和位置配置下,鸭式布局的升力较之单独主翼为高的主要原因是因为鸭翼有推迟主翼离体涡破碎的作用,鸭翼离体涡在主翼翼面上形成的负压以及鸭翼离体涡流动造成的主翼流场的变化,也是提高主翼升力的因素。     

旋成体导弹小展弦比舵面大偏度对称状态下非对称流动机理

针对跨声速条件下,小展弦比截尖三角翼尾舵的旋成体导弹在小迎角、零侧滑、大舵偏对称状态下呈现出的非对称流动现象,本文首次对其进行了分析研究。首先,通过一系列测力试验、表面油流试验及粒子图像测速（PIV）试验对该非对称流动现象进行了精准捕捉,并对其产生的原因进行了分析。然后,基于已获得的试验数据及流场观测结果,借助数值模拟方法对所述非对称流动的细节、拓扑结构、空间形态及舵面压力分布等问题做了深入研究,并进行了详细讨论。结果表明：旋成体导弹小展弦比舵面大偏度对称偏转时,舵面前缘产生的翼尖涡会因舵面相距较近而相互干扰,促使翼尖涡沿流向非对称发展,使得舵面压力分布不均,最终导致非对称流动和较大横向量的产生,影响导弹的气动性能。     

基于黏性涡模型的旋翼流场数值方法

 建立了一种适用于旋翼非定常流场特性分析的黏性涡数值方法。在该方法中:流场中的大尺度涡被离散为若干微小的涡元,通过求解涡量-速度形式的Navier-Stokes方程模拟涡元的输运等过程;黏性扩散效应采用高精度的粒子强度交换法进行计算,而桨叶附着涡以及新生涡环量采用了Weissinger-L升力面理论进行求解;为显著提高计算效率,在诱导速度及其梯度的计算中还引入了快速多极子算法(FMM)。应用上述方法,对悬停和前飞状态下的多个旋翼流场算例进行了计算,通过对比旋翼尾迹涡量特征和诱导速度分布等,验证了该方法的有效性。此外,还将本方法与旋翼计算流体力学(CFD)方法及传统的自由尾迹方法进行了比较,结果表明黏性涡方法在兼顾效率的同时,还能够更好地捕捉旋翼尾迹运动。     

DBD涡流发生器及其在角区流动控制中的数值研究

为明确介质阻挡放电(DBD)涡流发生器对马蹄涡的影响,在采用唯象模型的基础上,通过数值方法研究了DBD涡流发生器诱导产生流向涡的结构特性及其对马蹄涡的控制特性。研究结果表明,流动在激励器上游边缘处形成羊角涡,自由剪切层卷入流向涡的涡核并为其提供持续涡量;在柱体根部角区流动中,当对称面两侧激励器诱导流动指向对称面,诱导涡与马蹄涡环绕方向相反时,马蹄涡可以得到有效抑制,反之,则控制效果不佳。最后得出,诱导涡对下游马蹄涡的控制机制体现在其黏性扩散作用、掺混作用以及低压效应3个方面。     

基于多流路体积力模型的翼身融合布局飞机计算方法

针对翼身融合布局(BWB)飞机分布式系统一体化研究方法存在计算资源耗费巨大的问题,发展基于多流路体积力的飞机分布式系统一体化计算方法。从沿流线的体积力模型出发,将压气机对气流的作用力用沿流线的源项场代替。将单流路的计算模型加以扩展,综合考虑分布式推进系统与BWB的耦合影响,在分布式推进系统的压缩系统所在区域添加体积力源项,建立多流路体积力计算模型。将其应用于某大型BWB飞机进行一体化数值仿真。结果显示:相较于传统的全周定常数值模拟方法,采用多流路的体积力模型能够在保证一定计算精度的条件下提高至少10倍的计算效率。     

粒子分离器涡流叶片鸟撞击损伤试验

研究粒子分离器涡流叶片受鸟撞击损伤规律,提高粒子分离器抗外物损伤设计能力.通过外物冲击损伤试验方法,研究了模拟鸟撞击粒子分离器涡流叶片的过程.采用瞬态响应测试系统,获得了模拟鸟撞击涡流叶片的应变-时间历程.对采集到的曲线和数据进行了分析.利用三坐标测量仪,测得撞击前后涡流叶片叶型;经过比较计算,得到了涡流叶片的变形.模拟鸟以Ma=0.2速度撞击涡流叶片时,涡流叶片产生了塑性应变,变形很大.      

基于数值涡流发生器双S型进气道被动流动控制装置参数优化

采用数值涡流发生器代替真实涡流发生器叶片进行双S型进气道流动控制数值模拟,结合实验设计理论分析了27组涡流发生器应用于进气道流动控制计算结果,从中找出控制装置参数变化对进气道总压恢复、流场畸变的不同影响,应用响应面法给出最佳参数组合,为双S型进气道被动流动控制装置参数优化提供技术参考.     

旋涡发生器对叶片根部马蹄涡的影响

针对简化了的叶片-平板结构,在叶片上游平板上的某位置布置旋涡发生器影响马蹄涡.数值模拟了不同高度旋涡发生器作用下的角区湍流流动.结果表明:任何一种旋涡发生器都能不同程度地削弱马蹄涡系,高度为2mm的旋涡发生器能使这种削弱作用达到最大;另外,随着高度的增加,叶片两侧的涡腿有向叶片聚拢的趋势.提出的方法无需额外消耗能量,为一种被动的控制手段,有很大的实际应用前景.布置旋涡发生器的控制方法简单实用,工程中采用该法是不错的一种选择.