炮口二次燃烧的正激波点火模型

本文提出了考虑炮口激波结构中的正激波强度、化学动力学和湍流射流混合的预测炮口二次燃烧点火的模型,得到了与实验相一致的结果。     

含有大位移动边界的复杂流场的数值模拟

含有运动弹丸的膛口流场是典型的大位移动边界问题,同时弹丸在出膛过程中,流场的结构也会发生变化,增加了流场的复杂程度.在对该流场的数值模拟中,将它分为两个区:弹丸运动区和普通流场区,它们之间用一个特殊的分区边界联系,同时运用网格局部重构技术处理弹丸运动造成的网格变形问题,并将对称轴定义为网格变形边界,使得弹丸在对称轴上运动过程中不会导致对称轴上的网格体积为负.从计算结果可以看出整个膛口波系结构变化过程和弹丸先加速后减速的过程,从而表明该动网格处理方法是成功的.     

AUFS 格式在无网格方法中的应用简

将计算量小,激波分辨率高的AUFS (artificially upstream flux vector splitting) 格式应用于无网格方法. 所发展算法基于多项式基函数最小二乘无网格方法,采用线性基函数曲面拟合及AUFS 格式计算各离散点的空间导数,应用四阶Runge-Kutta 法进行时间显式推进. 为验证算法健壮性、精度以及计算效率,对Riemann 问题、超音速平面流动,以及不同攻角NACA0012 翼型跨音速流场进行了数值模拟,其结果同采用HLLC (Harten-Lax-van Leer-contact) 格式的无网格方法以及文献报道结果吻合较好,并且计算量较形式简单HLLC 格式减少约15%.     

AUFS格式在无网格方法中的应用

将计算量小,激波分辨率高的AUFS(artificially upstream flux vector splitting)格式应用于无网格方法.所发展算法基于多项式基函数最小二乘无网格方法,采用线性基函数曲面拟合及AUFS格式计算各离散点的空间导数,应用四阶Runge--Kutta法进行时间显式推进.为验证算法健壮性、精度以及计算效率,对Riemann问题、超音速平面流动,以及不同攻角NACA0012翼型跨音速流场进行了数值模拟,其结果同采用HLLC(Harten-Lax-van Leer-contact)格式的无网格方法以及文献报道结果吻合较好,并且计算量较形式简单HLLC格式减少约15%.     

含化学反应膛口流场的无网格数值模拟

基于无网格方法，对包含大位移运动边界和非平衡化学反应的膛口流场进行了数值模拟。所发展算法是基于线性基函数最小二乘显式无网格方法，忽略黏性及湍流的影响，对流场采用ALE(arbitrary Lagrangian-Eulerian)形式的Euler方程描述，对流通量和化学反应源项采用多组分HLLC(Harten-Lax-van Leer-Contact)格式和有限速率反应模型计算，对于运动边界造成的点云畸形采用局部点云重构方法处理，重构过程中采用虚拟边阵面推进。对圆柱绕流和激波诱导燃烧流场进行了数值模拟，验证了重构方法和化学反应计算的有效性。最后对12.7 mm口径机枪膛口流场进行了模拟，结果同实验照片、非结构网格方法结果吻合较好，数值结果清晰地再现了膛口初始冲击波、膛口冲击波、欠膨胀射流波系结构的动力学发展过程，以及膛口焰的时间、空间分布特征。     

计算含动边界非定常流动的无网格算法

在无网格算法中考虑了含动边界的流动问题,研究了可以计算处理包含一定位移及扭转动边界非定常流动的算法.创建了无网格算法的动点法则,并引入抗扭方法对弹簧方法进行改进来处理离散点运动,提高了方法的可用度及精度.发展了求解基于无网格的ALE方程组的算法,在点云离散的基础上采用曲面逼近计算空间导数及HLLC格式计算数值通量,运用四步龙格-库塔法进行时间推进.在跨、超音速条件下,计算模拟了典型翼型简谐振动流场,计算结果与实验结果及文献对比吻合,验证了该算法的正确性.     

激波点燃粉尘的数值模拟研究

建立了用于模拟入射激波后可燃粉尘颗粒点火的一维非定常两相化学反应流模型，该模型考虑了气固两相间的相互作用、粉尘颗粒的加速、加热和化学反应。粉尘颗粒着火前的化学反应用发生在颗粒外表面和内孔表面的非均相反应描述，颗粒内部的温度变化用一含有化学反应源项的非稳态热传导方程来描述，以颗粒外表面温度的突跃上升作为可燃粉尘颗粒点燃的着火条件。我们用该模型和ＰＳＩＣ方法，对由中等强度激波从纯气相传入煤粉－氧气混合物而引起的非定常两相流动现象，包括气固两相间的相互作用、粉尘颗粒的加速、加热以及点火过程进行了数值研究，计算了对应于不同载荷比、马赫数为４～５的入射激波后煤尘颗粒的点火延迟时间，分析了由于可燃粉尘颗粒的存在，入射激波及波后气固两相流动参数的变化规律。数值计算结果与实验数据符合较好。文中建立的模型和所用的基于ＰＳＩＣ算法的数值方法，用最自然的方式描述气固两相流动，即用连续流模型（欧拉方程）描述输运相（气相）的流动，用轨道颗粒模型（拉格朗日方程）描述分散相（颗粒相）的运动。用这种方法模拟含尘介质中激波后颗粒的点火是很有效的，它可以清楚地确定哪一个颗粒群最先着火，它的初始位置以及在整个点火延迟时间内     

驻定斜爆轰波并行数值模拟

采用多组分化学反应Euler方程组对驻定在高速飞行弹丸上的斜爆轰波流场进行了数值模拟。计算中分别采用TVD格式和基元反应模型,并基于并行编程模型MPI（message passing interface）实现了非结构网格上的并行计算,对流项和化学反应项用时间分裂法进行处理。计算结果表明并行计算能有效地提高计算速度,扩展计算规模,为进一步研究超驱爆轰推进技术奠定基础。     

高速动能弹温度场数值模拟

随着兵器发射技术和空气动力学技术的发展,动能弹的发射初速和飞行状态正从超声速向高超声速发展,由此产生了气动热问题.准确预测动能弹温度场是其气动力和热防护设计的关键技术.采用CFD预测温度场的方法,包括平衡流流动控制方程及差分格式,构造平衡流通量Jacob矩阵,在差分格式矢通量分裂过程中嵌入平衡流真实气体模型模拟温度场,获得平衡流气体状态方程.对典型高速动能弹热环境进行验证,考察方法的合理性.对设计的一种新型高超声速动能弹温度场进行数值模拟,为其气动设计及热防护提供了较可靠的数据.     

运动网格上的水下爆炸数值模拟方法

 在对水下爆炸过程的数值模拟中，把水-气介质界面当作一种特殊的内部边界，并且在该边界上的网格节点都有双重定义，分别对应边界两侧的不同介质状态，通过利用HLLC（Harten，Lax，van Leer，Contact）方法求解任意ALE（Arbitrary Largrangian Eulerian Method，拉格朗日-欧拉方程）方程组，以及运用界面跟踪和动网格方法，对一维球对称水下爆炸模型进行了数值模拟。结果表明，所得到的数值计算结果与实验结果吻合得较好，并且在扩展到二维柱对称水下爆炸模型时，同样得到了满意的结果。