旋翼非定常气动特性CFD模拟的通用运动嵌套网格方法

针对直升机旋翼非定常气动特性CFD模拟中的网格生成难题,提出了一套高效、通用的运动嵌套网格生成方法.首先,基于Poisson方程求解和翻折法生成旋翼桨叶的正交贴体网格.其次,针对旋翼桨叶的扭转分布及变距、挥舞等复杂运动,建立了一套通用的洞单元识别的扰动衍射法;为保证洞包络面的封闭性,完善了挖洞过程中网格加密策略;在洞边界确立基础上,提出了一种高效、鲁棒的最小距离法贡献单元搜索的改进方法.在此基础上,建立了基于RANS(Reynolds-averaged Navier-Stokes)方程的旋翼非定常流场CFD模拟方法.最后,采用所建立的方法分别对悬停和前飞状态下的C-T(Caradonna-Tung)旋翼和7A(Helishape 7A)旋翼的气动特性、桨尖涡的位置进行了计算,计算结果与试验值误差小于5%,验证了该运动嵌套网格生成方法在旋翼非定常气动特性CFD模拟中的有效性.      

小展弦比飞翼标模雷诺数影响数值模拟研究

采用数值模拟方法开展小展弦比飞翼标模的雷诺数影响研究。使用自行研制的多块对接结构网格亚跨超声速流场解算器程序Mbflow,计算了试验雷诺数、二倍试验雷诺数和飞行雷诺数等三种雷诺数情况下小展弦比飞翼标模的流场。通过对计算结果的分析研究,得到了不同马赫数(Ma=0.2、0.8和1.5)和攻角情况下雷诺数对小展弦比飞翼标模的气动特性和流场特征的影响规律。小攻角情况下,雷诺数主要影响摩阻的大小;而大攻角时,雷诺数对压阻也有明显影响。在亚声速时,雷诺数主要影响分离涡的起始位置和强度;跨声速和超声速时,雷诺数还会影响到激波的位置和强度。进一步研究了小展弦比飞翼标模的自准区雷诺数问题,发现试验雷诺数接近于小展弦比飞翼标模的自准区雷诺数。     

轴流压气机引气后流场的数值模拟

选取了NASA Stage 35作为算例进行压气机引气数值模拟研究.通过对比验证了网格无关性和计算的准确性,应用两种数值实现方法分别进行圆孔和槽引气计算,两种数值方法的结果基本吻合.进行了两种引气方式（即圆孔和槽）的对比.结果表明:槽引气转子出口总压损失小,效果好于圆孔引气.应用NUMECA软件自带引气模块操作简单,容易实现.绘制引气管路网格增加了流道网格绘制的难度和计算边界条件设置的复杂度,但更加符合压气机实际的引气情况.引气后压气机效率 流量曲线整体上移,总压比 流量曲线下降,失速裕度有较大提升.      

适用于倾转旋翼机气动分析的新型嵌套网格方法

为克服围绕复杂倾转旋翼机生成高质量一体化结构网格的困难,将分块网格策略与嵌套网格方法融合,提出并建立了一套适用于倾转旋翼机气动干扰分析的新型嵌套网格生成方法。该方法中,将整体网格按照倾转旋翼机结构特性进行了分块生成,在保证插值精度的基础上降低了网格规模。为高效连接不同部件网格,首先,引入Inverse Map方法来辅助部件网格定位,采用新颖的多向投影方法保证了部件间洞边界的连续性;其次,基于独立挖洞嵌套策略完成机身与背景网格的组合,并使用了相应的辅助计算坐标进行背景边界网格标识。在贴体网格区域采用基于S A（Spalart Allmaras）湍流模型可压缩RANS(雷诺平均Navier Stokes）方程求解流场,背景网格区域采用欧拉方程,并运用了SPMD(single program, multiple data)模式的并行加速技术,建立了适合倾转旋翼机气动干扰特性分析的高效混合CFD方法。采用Robin直升机机身作为数值算例,验证了CFD方法的有效性。在此基础上着重对倾转旋翼机在悬停、过渡和前飞巡航模式下的干扰流场进行了数值分析研究,分别得出了悬停状态下倾转旋翼机典型的“喷泉效应”干扰现象、过渡状态下的非定常气动干扰特性以及巡航状态下的部件气动变化特性,计算结果表明建立的新型嵌套网格方法能够较好地表征倾转旋翼机外形特性,并能够有效地用于倾转旋翼机的气动特性分析,加速比能超过5.0。      

各向异性网格自适应计算在超燃模拟中的应用

采用各向异性网格自适应求解技术,将其应用于DLR超燃冲压发动机燃烧室中的超燃模拟.开展了3个算例,包括采用滑移壁面条件的燃烧室冷流场模拟、采用无滑移壁面条件的冷流场模拟及采用无滑移壁面条件的反应流模拟.模拟中,各向异性网格自适应计算捕捉到了如激波、射流、边界层、火焰面等具有各向异性特征的大梯度区域,并利用各向异性网格进行了很好的加密.对比利用各向同性网格的初始流场计算,各向异性网格自适应计算使基于滑移条件、无滑移条件的冷流计算及反应流计算的网格单元数量分别下降了36.2%,36.4%和36.8%,有效降低了计算规模,而且流场大梯度区域的计算结果更准确,辨析度更好.结果表明:对于像超燃这类具有各向异性特征的问题,各向异性网格系统比各向同性网格系统有更好的计算效率及准确性,同时也表明基于Mach数场构造的各向异性网格系统可以有效应用于超燃计算.      

垂直/短矩起降飞机机翼内埋式风扇布局气动特性分析

以垂直/短距起降飞机过渡飞行状态为背景,针对机翼内埋式风扇布局的自由来流/风扇喷流混合型流动,基于结构/非结构混合网格使用CFD方法进行了非定常数值模拟和分析.首先使用滑移网格技术对NASA涵道螺旋桨进行算例验证,其时均计算结果与实验值的误差为5.3%,证明了计算方法的可靠性和准确性,然后数值模拟了机翼内埋式风扇布局在不同迎角下的气动性能.结果表明:风扇喷流在机翼上产生了特有的“抽吸”和“堵塞”效应,引起了机翼总升阻力的显著增加,升力最大增量达到干净机翼升力的2.6倍,阻力最大增量为干净机翼阻力的3.2倍,混合流场在机翼后缘引起了升力损失并卷起对涡.      

基于冷态叶型生成叶片加工坐标的方法

建立了基于冷态叶型生成叶片加工坐标的方法,包括排序、插值加密、叠加位移量、去倒圆及端区叶型拟合等步骤.以某多级轴流压气机为例,应用该方法生成的加工坐标叶片除端区外主流通道叶高内误差在0.01mm以下,叶片表面波纹度和前后缘曲率良好.通过进口级的数值模拟表明:该方法生成的加工坐标叶片与冷态叶型的误差对气动性能影响小,流量-总压比-等熵效率特性变化微小,叶片进出口气流速度、方向及流场细节等变化较小,满足工程应用的需求.      

三种T700级碳纤维及其复合材料性能比较

对MT700、T700-A及T700-B三种碳纤维拉伸性能、表面形貌、单向板力学性能及网格加筋圆筒轴压稳定性进行逐级对比研究。结果表明:MT700碳纤维拉伸性能达到同级别进口碳纤维水平且具有高模量特征;MT700碳纤维表面均布沟槽的结构特点使得MT700/603复合材料体系表现出良好的界面性能和拉伸-压缩匹配性,单向板压缩强度、层剪强度及弯曲强度均明显高于T700-A/603和T700-B/603;MT700/603网格加筋圆筒轴压破坏强度及模量分别达到870 k N和108.2 GPa,相比于T700-B/603分别提高11.5%和33.1%。MT700碳纤维更适用于制备航天领域结构复杂承力构件。     

螺旋桨飞机滑流非定常数值模拟研究

采用非定常数值模拟方法,基于滑移网格技术,分析研究螺旋桨滑流对全机气动特性的干扰影响。用多块结构网格数值模拟软件,从搭接边界建立、搭接网格生成、流场插值传递3个方面对螺旋桨滑移搭接网格的数值模拟方法进行全面详细的介绍,并且阐述基于该技术的非定常数值模拟方法。在此基础上,开展螺旋桨飞机模型验证计算,分析研究螺旋桨滑流对全机流场和气动性能的影响,通过与风洞实验数据的对比分析,验证滑移搭接网格技术的可行性与准确性。研究结果表明,基于滑移搭接网格技术的螺旋桨滑流数值模拟,可以真实反映螺旋桨滑流流场特性,为螺旋桨类飞机的设计评估提供一种有效的数值模拟手段,在工程方面具有重要的实用价值。     

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。