高精度差分格式及多尺度流场特性的数值模拟

简要介绍了作者近年发展的高精度差分格式,给出了两类格式精度的分析方法,分析了粘性面差分逼近所引入的耗散量低亏效应。文中按群速度将格式分为快型,慢型和混合型三类。对每类格式分析了高频波的传播特性。     

基于离散曲率计算的三角网格模型优化调整

采用面积夹角加权的三角网格模型顶点法矢及三角片质心权值对Taubin的三角网格模型离散曲率计算方法进行了改进,在此基础上提出了一种新的三角网格模型优化调整方法。用该方法调整三角网格模型,在模型上曲率变化较平缓的平坦区域及曲率变化较剧烈的特征区域都能取得较好的调整效果。     

二维任意域内基于节点的局部网格生成算法

凸域内基于节点的局部网格生成算法,克服了基于节点的有限元方法的网格生成可能产生的不一致性。将该基于节点的局部网格生成算法的适用范围拓展到二维任意域。另外,提出了通过使用约束Delaunay路径来划分任意域的区域划分算法,该算法使得在并行实现网格生成的过程中各处理器之间无需通信,从而大大提高了节点给定情形下有限元方法网格生成的并行效率。     

高速磁浮列车会车压力波数值计算研究

根据可压缩粘性流体的N-S方程和k-ε双方程湍流模型,运用移动网格技术,采用有限容积法对磁浮列车以500 km/h速度运行时的会车压力波进行了数值计算研究,计算结果与现车试验结果相差10%以内,表明了本文计算模型和方法的正确性;文中还研究了交会过程中列车周围的流场分布,结果表明:当列车交会时,两车之间的流场相互叠加,而两车之外的流场分布在整个交会过程中变化不大。     

高超声速二元变几何进气道气动方案设计与调节规律研究

 针对二元高超声速进气道宽马赫数大攻角工作要求,研究了转动唇口变几何进气道调节方案,给出了一种高性能变几何原型进气道设计方法,并通过数值仿真获得了变几何进气道各工况下的调节规律及性能变化规律。研究表明:采用特殊曲面唇口设计的变几何进气道宽马赫数范围内流场结构较好,总体性能优越;以马赫数Ma=6.0为设计点的原型进气道采用转动唇口方案无需附面层抽吸即可在唇口开启过程中实现接力点自起动,且最低自起动马赫数降至Ma=3.5;低马赫数大攻角状态下,通过转动唇口合理控制喉部平均马赫数范围可保证进气道正常工作。     

两车速度对超车过程轿车气动特性的影响

利用动网格技术对超车过程中两种相对速度下车辆周围非定常流场进行了数值模拟,得到超车车辆和被超车辆受到的侧向力、侧倾力矩和横摆力矩在超车过程中的变化情况,并根据两种速度情况下的流场特性,分析其对气动力特性的影响.计算结果表明:超车速度越大,被超车辆最大侧向力、侧倾力矩和横摆力矩越大;超车车辆最大侧向力、侧倾力矩越大,横摆力矩越小.      

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

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

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

前后缘型线同时可控的乘波体设计

提出了一种前后缘型线同时可控的乘波体设计方法,在马赫数可控的外锥形曲面基准流场中,结合流线追踪技术和混合函数,实现了椭圆前缘转椭圆后缘的乘波体设计,并在设计点(Ma=6.0)和接力点(Ma=4.0)对其进行数值仿真研究。前后缘同时可控的乘波体在型面剧烈过渡处产生了较弱的激波,出口两侧存在高温高压区,后部对称面附近的激波形状由圆弧变为平直线且出口处流场基本均匀,非常有利于与进气道匹配设计。另外,该乘波体具有较高的容积率和预压缩效率,附面层修正后的容积率为0.24,设计点时乘波特性较好,接力点时前部完全乘波,具有较高的升阻比,有黏条件下设计点和接力点的升阻比分别为2.54和2.41。此外,与给定前缘的乘波体相比,其升力、阻力、俯仰力矩和出口增压比都有明显增加,但是升阻比和出口总压恢复系数有所降低,在设计点无黏升阻比由3.56降为3.00。以上研究表明,本文的设计方法可行且更加灵活,拓宽了乘波体的选择范围。     

化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。