DLR-F6翼身组合体的高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6翼身组合体的高阶精度数值模拟研究。主要目的是确认WCNS模拟跨声速典型运输机构型的能力。采用粗、中、细3套网格开展了网格收敛性研究,从气动特性、压力分布、表面流态等方面研究了网格密度对DLR-F6翼身组合体气动特性的影响;采用中等网格开展了来流迎角对气动特性的影响研究。通过与试验数据、CFL3D软件和TRIP软件计算结果的对比,表明网格密度主要影响激波位置和压差阻力系数,同时对翼身结合部分离区大小有一定影响;采用高阶精度计算方法显著提高了气动力系数的模拟精度,力矩系数数值模拟结果与试验的差异有待进一步分析。     

风沙两相流数值模拟研究进展

大气表面层中的风沙运动是典型的高雷诺数湍流气固两相流动,同时具有可蚀地表上的冲击溅起过程、风-沙-电多场耦合、多尺度与跨尺度等复杂特性,使其定量预测极具挑战性。回顾了近30年风沙运动数值模拟的研究进展以及取得的成果,主要包括基于湍流雷诺平均Navier-Stokes （RANS）方法的风沙运动多场耦合模拟与风成地貌耦合尺度模型建立,基于湍流大涡模拟（LES）与高雷诺数湍流理论的风沙运动模拟方法。最后明确了风沙运动现有数值模拟方法的不足和值得进一步研究的方向。     

战斗机垂尾脉动压力数值模拟

在亚跨超计算流体力学（CFD）软件平台（TRIP）上开发了基于RANS/LES混合思路的IDDES流动模拟技术,并通过NACA0021翼型60°大迎角分离流动与串列圆柱绕流模拟对RANS/LES混合方法的精确度进行了验证,针对某战斗机外形的垂尾脉动压力开展了数值模拟研究。战斗机来流马赫数为0.1,基于全机长度的雷诺数为2×10⁶,模型迎角为20°、30°和40°。分别通过脉动压力系数、脉动压力功率谱密度、空间流动结构以及侧向力响应曲线等对战斗机的垂尾脉动压力进行了分析。脉动压力模拟结果表明：当垂尾完全沉浸在边条翼脱体涡破碎后的宽频湍流脉动气流中时,垂尾翼梢位置的脉动压力会发生明显的增大。     

6种典型飞行器的RANS计算及大分离区域的DES分析

提出了一种将RANS(Reynolds averaged Navier-Stokes)与DES(detached eddy simulation)相结合计算流场的工程新算法,并用于第一代载人飞船Mercury、第二代载人飞船Gemini、人类第一枚成功到达火星上空的Fire-Ⅱ探测器、具有丰富风洞实验数据(来流Mach数从0.50变到2.86)的巡航导弹、高升阻比的Waverider(乘波体)以及具有大容积效率与高升阻比的CAV(common aero vehicle)等6种国际上著名飞行器的绕流计算.在流场计算中,采用分区技术,即首先对全流场采用RANS计算,然后对分离较大或者分离较严重的那些区域采用DES分析技术.文中所完成的6个典型算例总共63个工况的数值计算表明:仅在分离较大的区域采用DES分析技术,虽然从严格意义上讲这样的处理并非真正意义上的DES,但这样近似处理后所得到的流场数值结果(其中包括气动力和气动热)与相关实验数据较为贴近,并且流场的计算效率较全场进行DES计算时要高得多.      

一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

内管射流占优时轴对称共轴射流结构分析

以轴对称Reynolds-averaged Navier-Stokes(RANS)方程为控制方程模拟了内流占优条件下轴对称共轴喷嘴射流流场,并基于湍动能、速度和压强等流场参数的变化规律,将这类流场划分为交汇段、初始段和主流段等三个区段,给出三段流场的结构特征和流动机理,发现初始段内湍流动能大大低于单管射流,将导致共轴射流噪声低于单管射流;主流段内湍流动能则高于单管射流,将加强射流的掺混作用.将流量计算结果与试验结果进行了对比,证明使用的算法是有效的.      

基于Delaunay图映射的动弹网格技术及其应用

在多块混合网格的基础上发展了一套基于Delaunay图映射方法的动弹网格技术,实现了针对气动弹性计算的网格变形运动方法,并使用基于多块混合网格的RANS方程的有限体积法耦合结构静平衡方程,针对大展弦比无人机在实际使用中可能出现的显著的静气动弹性问题进行了计算、分析和研究。研究结果表明,所发展的动弹网格方法为粘性混合网格的变形提供了较好的解决方案;RANS方程解算器具有较高的精度和较强的适用性,具备求解粘性流体与结构耦合等复杂航空气动弹性问题的能力。     

用强耦合RANS方法模拟旋翼悬停流场

在旋转坐标系下,将Spalart-Allmaras(S-A)一方程湍流模型和Reynolds-averaged NavierStokes(RANS)方程耦合成一个新的RANS方程,并发展了基于多块重叠网格的强耦合RANS求解方法,用于直升机旋翼悬停流场的数值模拟.为了提高计算效率,针对多重网格方法在多块重叠网格上实施的困难,提出了一种基于重叠网格的多重网格实施方法.通过对Caradonna-Tung(C-T)和ONERA 7A旋翼悬停算例验证了发展的强耦合RANS方法和基于重叠网格的多重网格实施方法的有效性.研究结果表明:发展的基于重叠网格的多重网格方法有较高的计算效率,3层网格的加速比约为7.7;强耦合RANS法的计算精度明显高于传统的松耦合RANS方法,特别是在与阻力相关性能参数的预测中,强耦合RANS方法的预测结果更加精确.     

基于混合笛卡尔网格的翼型绕流数值模拟

在对重叠网格技术研究的基础上,一套以笛卡尔网格为主体、贴体的结构网格为特殊固壁边界的混合笛卡尔网格(Hybrid Cartesian Grid,HCG)生成方法得到了发展,同时基于流场特征的笛卡尔网格自适应技术和有限体积方法的雷诺平均Navier-Stokes(RANS)的数值求解方法也得到了开发。通过对经典跨音速RAE2822翼型、双NACA0012翼型和高升力两段NLR-7301翼型进行绕流数值模拟,获得了较好的结果,验证了方法的可行性与精确性,为研究二维复杂外形流动问题提供了一种解决方法。     

k-kL两方程湍流模型的改进及验证研究

提高湍流数值模拟的准确性,从而明确湍流模型对数值模拟结果的影响具有重要的意义,应用K.S.Abdol-Hamid给出的尺度自适应k-kL两方程模型封闭RANS方程,并修改von Karman长度尺度的限制方法,通过平板、翼型、后台阶等流动的模拟,考察k-kL模型在湍流模拟中的准确性,及其反映主要流动特征的能力和网格收敛性,并对影响流动模拟准确性的因素进行讨论.结果表明:改进长度尺度限制之后的k-kL两方程模型无论是对附着流动还是分离流动都可以给出比较准确的结果.