使用高阶逆风通量差分裂格式的悬停旋翼流场数值模拟

为减少Jameson二阶中心差分有限体积法导致的旋翼尾迹的数值耗散,将三阶逆风格式(MUSCL)与通量差分裂方法相结合,建立了一个三维雷诺平均N-S方程数值模拟悬停状态旋翼流场的方法。为了充分考虑旋翼尾迹对流场的影响,采用周期性边界条件和由动量理论导出的远场边界条件。为进一步减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。然后,进行了算例计算,给出了桨叶表面的压强分布,与可得到的试验数据及二阶中心差分方法的计算结果进行了对比,并针对几种不同桨尖形状的旋翼悬停流场进行计算,数值结果显示:后掠桨尖可减弱超临界流动。此外,还计算和分析了旋翼下方不同位置上的涡量分布。计算结果表明,本文方法能够有效地减少旋翼尾迹的数值耗散。      

径向双旋流燃烧室流场结构大涡模拟研究

为深入了解真实航空发动机燃烧室内流场结构,在自有CFD平台上采用动态亚网格湍流模型对一种径向双旋流环形燃烧室的单个头部构型冷态流场进行了大涡模拟。为保证模拟精度,没有对模型进行常规简化处理,对包括全部气膜冷却小孔在内的所有精细结构均进行了完全仿真。计算验证了程序对高度复杂流场的模拟能力,结果表明,大涡模拟能较为全面地反映燃烧室内复杂流场从静止启动到统计定常的非定常发展过程,并成功捕捉到流场中心回流区等各种大尺度结构及涡旋破碎泡等旋流特征;大涡模拟所获得的时间平均流场结构与已有PIV试验结果定性一致,与RANS计算相比更接近试验测量值。     

旋转机翼飞机旋翼模式前飞状态干扰气动特性

和传统直升机相比,旋转机翼（CRW）飞机在旋翼模式前飞时各部件之间存在更为严重的气动干扰。为了获得旋转机翼/机身/鸭翼/平尾之间的非定常气动干扰规律,基于运动嵌套网格技术,通过求解三维非定常雷诺平均Navier-Stokes（URANS）方程,建立了旋翼前飞流场数值模拟方法。首先对传统直升机旋翼/机身干扰模型进行了计算,验证了方法的可靠性,然后对某旋转机翼飞机全机在旋翼模式前飞状态下的非定常流场进行了数值模拟,并对各个气动部件上的非定常气动力和力矩的变化进行了分析。结果表明：飞机在旋翼模式前飞时,机身部件对旋转机翼的干扰较弱,在经过机身上方时拉力峰值仅略有增加;旋转机翼对鸭翼和垂尾干扰较弱,对机身和平尾干扰较强,随着前飞速度增大,旋转机翼对平尾的干扰会产生较大的升力损失和抬头力矩,需要引起重视。计算结果为该类飞行器的总体综合设计提供了参考。     

基于高阶WENO格式的旋翼非定常涡流场数值模拟

为提高直升机旋翼黏性涡流场计算流体力学（CFD）的计算精度及降低数值耗散,发展了一套基于五阶加权本质无振荡（WENO）格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统,以Navier-Stokes方程作为主控方程,湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度,采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值,选取适当的加权因子,通过多个重构模板的凸组合来构造五阶WENO格式,在交界面附近获得具有五阶精度的状态变量。为了提高计算效率,采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）推进方法和并行计算策略。最后,运用该方法分别对悬停状态下C-T（Caradonna-Tung）旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟,并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明：在计算时间增长10%～20%的前提下,五阶WENO格式能够捕捉到更精确的涡流场特性,反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。     

基于伴随方程的网格自适应及误差修正

基于流场方程的离散伴随优化理论和三维非结构网格,建立了网格自适应技术和目标函数误差修正方法。详细研究了用流动伴随变量进行目标函数的误差估计和修正技术,构造了适用于格心格式有限体积法的流场变量插值技术和网格单元剖分判据,初步实现了网格物面投影和空间单元优化,发展了适用于有限体积法的整套网格自适应方法。对NACA0012翼型和ONERA-M6机翼绕流进行了自适应数值模拟,并对升、阻力等目标函数进行了误差修正。数值结果表明,本文自适应方法能正确地捕捉到影响目标函数计算精度的敏感区域,网格自适应和误差修正两项技术显著提高了升、阻力等目标函数的计算精度。     

基于自适应混合网格的脱体涡模拟

基于混合网格和CGNS（CFD General Notation System）数据结构,建立了一种各向同性加密/稀疏的网格自适应方法。在悬空点的后处理中,让含有悬空点的单元转化为任意多面体,从而简化了自适应单元剖分模版,同时自适应网格单元之间可完全相容,自适应生成的网格能够直接用于可处理任意多面体的流场求解器。将该自适应方法与脱体涡模拟（DES）算法相结合,开展了65°后掠三角翼大迎角流动的数值模拟应用,并与初始网格的模拟结果进行了详细比较。对比表明：采用网格自适应方法适当增加局部网格量,能够以较小的成本迅速提高三角翼背风区的空间分辨率,增强数值模拟对小尺度涡系结构的解析能力,从而弥补了基于混合网格的脱体涡模拟中常用二阶格式计算的空间分辨率相对偏低、不利于湍流多尺度结构精细模拟的不足。     

空中加油系统的多体动力学和CFD仿真

介绍了作者基于Fortran程序开发的一个多体动力学求解器.依托多刚体系统理论、基于绝对节点坐标(ANCF)法描述的常质量/变质量索梁单元和计算多体动力学方程组求解技术,该求解器可用于分析飞行器拖曳结构中的索-刚体约束系统的动力学问题.本文基于并行CFD求解器PMB3D的动态重叠网格功能,构建了多体动力学和计算流体力学的非定常耦合计算程序.以软式空中加油系统为例,对真实外形飞行器拖曳系统的气动-动力学特性进行了仿真研究.探讨了战斗机构型进行软式加油对接时产生的气动干扰问题,研究了受油机的接近速度和初始偏移量对对接效果的影响.     

某型倾转旋翼机的旋翼桨叶气动优化设计

倾转旋翼机的旋翼气动外形设计需要对其直升机模式和固定翼模式下的不同要求进行综合考虑。通过基于Kriging模型的多目标遗传算法建立一套适用于倾转旋翼桨叶气动外形优化设计方法,采用拉丁超立方抽样试验设计方法得到样本点,并建立Kriging模型替代费时的流动数值模拟。以最大化旋翼地面悬停拉力和高空巡航效率为目标,以地面悬停功率不增和巡航拉力不减为约束条件,进行倾转旋翼桨叶平面形状优化设计,并经过非定常数值模拟和风洞试验验证。结果表明:数值模拟结果和风洞试验结果吻合良好,优化结果满足设计指标。     

物理质团法关于可压流计算的程序验证

为了评估无网格物理质团法(PECM)对可压流计算的可信性,结合物理认识对该方法的设计思想和数值格式进行了先验分析,基于9个一维模型数值解与参照解的对比信息,通过后验分析方法对PECM进行了程序验证,参照解包括无网格再生核粒子法(RKPM)的数值解以及精确解。验证结果表明:各种物理参数以及状态方程的强烈间断会导致RKPM严重的失稳或失真,而PECM在包括极大密度比在内的各种强间断条件下仍保持稳定和收敛并具有较高的准确度。      

离散伴随方法在气动优化设计中的应用

针对工程问题中大规模设计变量难以进行快速优化的难题,在自研流场解算器(WiseCFD-UG)基础之上,发展了基于离散伴随方法的梯度计算模块,实现了目标函数对大规模设计变量的快速精确梯度求解.此外,结合FFD参数化、线弹性体网格变形,优化算法等模块建立了气动优化设计系统.运用系统对M6机翼和CRM翼身组合体构型进行定升减阻优化设计,结果表明:机翼表面的压力分布得到了显著改善,表面激波得到不同程度的消除或减弱.气动优化设计系统有效可靠,具有一定的工程应用价值.