热流CFD计算中格式和网格效应若干问题研究

通过数值实验、理论分析、试验数据和经验公式的综合研究,探讨了CFD热流计算的格式效应及网格效应:研究了Roe的FDS格式、Van Leer的FVS格式、AUSM 格式及中心差分格式的热流计算分辨率,探讨了不同格式对热流计算精度的影响,认为AUSM 格式在热流计算精确性方面具有优势;在热流绝对值的模拟上,网格的影响很大,能否计算出正确的驻点热流值将是获得精确物面热流解的关键;通过大量的数值计算、理论分析及其同试验数据的对比研究,对可能影响热流计算的网格因素进行了细致的分析,揭示了各因素对热流计算结果的影响规律。     

NS方程激波计算的摄动有限差分方法

摄动有限差分(PFD)方法从一阶迎风差分格式出发,将差分系数展开为网格步长的幂级数,通过提高修正微分方程的逼近精度来获得更高精度的差分格式。由于格式基于一阶迎风格式,因此具有迎风效应、网格节点少等特点。本文首先通过对Burgers方程的摄动差分格式的推导,将摄动有限差分格式引入时间相关法的计算,并构造了守恒形式的摄动有限差分格式,然后推广到一维Navier-Stokes方程组的计算。数值比较研究表明:本文构造的NS方程摄动有限差分格式具有比一阶迎风较高的精度和分辨率,而且保持了一阶迎风格式的无振荡性质。     

涵道尾桨气动特性及翼型CFD分析

以Jameson/TVD混合格式对低速时的OAF和NACA63a312翼型的升阻特性进行了计算。该方法的特点是采用有限体积法进行离散,并用龙格-库塔时间推进法进行求解。通过将计算的结果与NASA公布的试验曲线进行对比,验证了该混合格式用于旋翼气动力计算的可行性。为了进一步研究旋翼的拉力和功率,本文建立了以微圆环动量理论和叶素理论为基础的计算模型,并用该模型计算了涵道尾桨在不同条件下的拉力和功率情况,计算的结果具有一定的代表性。      

高阶精度格式WCNS-E-5在亚跨声速流动中的应用研究

采用高阶精度非线性紧致加权格式WCNS-E-5和Baldwin-Lomax模型,求解雷诺平均Navier-Stokes方程,开展了典型翼型与机翼的湍流流动数值模拟研究.对方程中粘性项采用的四阶精度差分近似以及网格导数求解与边界格式的四阶精度,保证了高精度算法的实现.计算结果表明:本文算法能够准确地模拟这些翼型与机翼的亚跨声速流场,得到与实验测量十分吻合的壁面压力分布,计算结果对网格的依赖性小.     

超声速底部喷流干扰流场数值模拟

数值模拟了不同马赫数,不同喷流压比下的轴对称超声速底部喷流干扰流场,采用LU隐式算法进行数值求解并引入了Baldwin-Lomax代数湍流模型.采用分区网格将弹身与底部区域合为一个整体进行计算,得到了清晰的流场结构和弹体表面及底部的压力分布,并与试验结果进行了比较.数值模拟结果表明超声速底部喷流干扰流场结构复杂,有、无喷流时底部流场有很大不同, 喷流对底压分布有明显影响,进而对轴向力系数影响显著.     

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

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

高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

预处理AUSMDV格式在微型扑翼绕流数值模拟中的应用

由于微型扑翼的低速低雷诺特性,对其粘性绕流的数值模拟存在较大困难。传统的中心空间离散格式不具有格式耗散而需要加入带有自由度的人工耗散,给计算带来不确定因素;此外,在低速低雷诺条件下,RANS方程组的条件数太大,导致计算收敛缓慢甚至不收敛,需要引入预处理方法加速收敛。本文使用具有低耗散特性的AUSMDV空间离散格式结合预处理方法对扑翼的非定常粘性绕流进行了数值模拟,非定常推进采用了含双时间迭代的LU-SGS隐式时间推进,湍流模型采用了BL模型。计算结果表明,加入预处理能够克服AUSMDV格式难以收敛的缺点,并且与中心格式相比,基于预处理的AUSMDV格式能够更准确地模拟扑翼的非定常粘性绕流,通过将本文计算结果与文献结果对比验证了本文方法的有效性。在此基础上研究了扑翼相关参数对气动特性的影响,对进一步了解扑翼气动力产生机理,指导扑翼的气动设计具有参考价值。     

气冷涡轮级流场的数值模拟方法与试验验证

采用具有TVD(total variation diminishing)性质的三阶精度Godunov格式、自由型曲面网格生成技术以及分区网格算法,对某型涡轮级进行考虑冷气掺混的全三维Navier-Stokes(N-S)方程数值求解,并将所得的结果与试验数据进行对比分析.通过分析数值模拟结果研究了该涡轮级所具有的气动特点.结果表明:所开发的快速、高效冷气掺混网格生成方法以及任意分区的流场求解算法可以满足工程上对气冷涡轮级的总体性能的快速估算及流场结构的详细描述.