基于非结构嵌套网格和逆风格式的旋翼悬停流场数值模拟

为更好地计及桨叶间相互影响,同时捕捉桨尖附近的流场细节,发展了一套非结构运动嵌套网格方法,并针对关键的贡献单元搜寻问题提出了一种高效的搜索方法。为了可方便地扩展到前飞状态的计算,同时为提高旋翼流场涡量的捕捉精度和效率,采用双时间法求解惯性坐标系下的非定常N-S主控方程,空间方向上使用二阶迎风格式,在此基础上,将不同精度格式的计算结果与试验值进行对比,并讨论了桨尖形状变化对旋翼气动特性的影响,结果表明本文建立的CFD方法对旋翼悬停流场的模拟是很有效的。

非结构网格迎风格式中的一种重构方法研究

研究了一种基于网格节点重构的二阶迎风格式在非结构网格无粘流场计算中的应用 ,提出了一种新的网格节点流场参数的重构方法。并通过 ONERA M6机翼、带副翼的直机翼、翼身组合体以及翼身组合体带增升装置无粘绕流流动计算对该重构方法进行了验证。算例表明新方法可以在确保格式精度和不显著增加流场计算时间的前提下 ,明显改善迎风格式的稳定性、收敛性和鲁棒性 ,减小了网格质量对流场求解精度的影响 。

非结构网格上的TV-van Leer混合迎风格式

为了改进迎风型对流通量格式的激波捕捉计算稳定性,对迎风格式产生数值激波不稳定现象的原因进行分析,并且在分析结果基础上构造新型迎风型格式。为此,首先设计了数值算例,分析非结构网格上迎风格式出现数值激波不稳定现象的主要影响因素。结果表明,在一定条件下,网格和流动状态共同作用导致产生不合理的数值结果,并影响计算的稳定。其中,网格间的动量通量对数值激波不稳定现象影响显著。因此,基于TV格式构造了一种自适应调节动量通量耗散项的混合格式,从而在激波附近引入耗散,抑制扰动的产生和传播。数值计算结果显示新型混合格式在非结构网格离散条件下能够稳定地捕捉激波,同时混合格式保持了TV格式原有的黏性流动计算分辨率。

非结构网格高超声速化学非平衡MHD流场数值模拟

在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学(MHD)方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。

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

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

半正定两相驱动问题的多步有限体积方法及其理论分析

考虑多维半正定两相驱动方程的初边值问题,在非结构网格上构造多步的迎风有限体积格式,利用微分方程先验估计理论证明了格式的离散模形式的误差估计为D(△t^2+h),其中△t和h分别表示时空步长．数值算例进一步验证了格式的有效性．