弹道枪不同水深下全淹没式发射膛口流场的数值分析

张京辉 余永刚

引用本文:
Citation:

弹道枪不同水深下全淹没式发射膛口流场的数值分析

    作者简介: 张京辉(1995- ),男,博士研究生,zjhsg@vip.qq.com;
    通讯作者: 余永刚, yygnjust801@163.com
  • 中图分类号: O354.5; TJ6

Numerical investigation on the muzzle flow field of an underwater submerged launched ballistic gun at different water depths

    Corresponding author: YU Yonggang, yygnjust801@163.com
  • CLC number: O354.5; TJ6

  • 摘要: 为了解弹道枪水下全淹没发射时,水深对膛口流场演化特性的影响,建立了二维轴对称非稳态膛口流场模型。采用流体体积函数多相流模型、标准k-ε湍流模型和Schnerr-Sauer空化模型,结合动网格及用户自定义函数技术,对水下全淹没发射膛口流场演变全过程进行了数值模拟。搭建了弹道枪水下可视化射击实验平台,对12.7 mm口径弹道枪在水中全淹没式发射时膛口流场演化过程进行了观测,并验证了数值模型的合理性。在此基础上,对比了不同水深下(h=1~100 m)膛口流场的演化特性。通过对比发现:在不同水深条件下,在膛口流场影响范围内,弹丸膛外行程随时间的变化均满足指数函数规律;水越深,膛口流场典型波系结构形成所需时间越长,且燃气在膛口轴向马赫盘处的温度和压力峰值越低,压力振荡幅度也越小,更快趋于平稳,但在径向上,水越深,压力振荡持续时间越长。
  • 图 1  计算域及边界条件

    Figure 1.  Calculation domain and boundary conditions

    图 2  不同网格数下膛口中心速度沿轴向的变化(t=1 ms)

    Figure 2.  Velocity of the muzzle center varying along the axial direction for different grid quantities at t=1 ms

    图 3  实验装置示意图

    Figure 3.  Schematic diagram of the experimental device

    图 4  实验照片阴影图与模拟相图对比情况

    Figure 4.  Comparison of experimental shadow photos and simulated phase diagrams

    图 5  射流轴向、径向最大位移对比

    Figure 5.  Comparison of the maximum axial and radial displacements of the jet flow between experimental and simulated results

    图 6  不同水深下弹丸膛外行程

    Figure 6.  Displacement-time curves of the projectile at different water depths

    图 7  空气中发射时典型膛口流场流谱示意图

    Figure 7.  Schematic diagram of the typical muzzle flow field during air launch

    图 8  全淹没发射时膛口流场流谱示意图

    Figure 8.  Schematic diagram of the muzzle flow field during submerged launch

    图 9  h=1 m时膛口中心剖面马赫数云图与等值线图

    Figure 9.  Mach number cloud map and contour map at h=1 m

    图 10  h=50 m时膛口中心剖面马赫数云图和等值线图

    Figure 10.  Mach number cloud map and contour map at h=50 m

    图 11  h=100 m时膛口中心剖面马赫数云图和等值线图

    Figure 11.  Mach number cloud map and contour map at h=100 m

    图 12  h=50 m是膛口中心剖面压力云图

    Figure 12.  Pressure cloud map of the muzzle center section at h=50 m

    图 13  膛口压力分布曲线

    Figure 13.  Pressure distribution curves of muzzle

    表 1  膛口初始参数

    Table 1.  Initial parameters for the muzzle

    h/ml/mv0/(m·s−1)pk0/MPa
    1123014.5
    50122015.5
    100120720.5
    下载: 导出CSV

    表 2  拟合参数

    Table 2.  Fitting parameters

    h/mx0/mx1/mt1/ms
    11.091.094.55
    500.840.843.55
    1000.790.793.42
    下载: 导出CSV
  • [1] 李鸿志, 姜孝海, 王杨, 等. 中间弹道学[M]. 北京: 北京理工大学出版社, 2015: 10.
    [2] 姜孝海, 范宝春, 李鸿志. 膛口流场动力学过程数值研究 [J]. 应用数学和力学, 2008, 29(3): 316–324. DOI: 10.3879/j.issn.1000-0887.2008.03.006.
    JIANG X H, FAN B C, LI H Z. Numerical investigations on the dynamic process of the muzzle flow [J]. Applied Mathematics and Mechanics, 2008, 29(3): 316–324. DOI: 10.3879/j.issn.1000-0887.2008.03.006.
    [3] 吴伟, 许厚谦, 王亮, 等. 含化学反应膛口流场的无网格数值模拟 [J]. 爆炸与冲击, 2015, 35(5): 625–632. DOI: 10.11883/1001-1455(2015)05-0625-08.
    WU W, XU H Q, WANG L, et al. Numerical simulation of a muzzle flow field involving chemical reactions based on gridless method [J]. Explosion and Shock Waves, 2015, 35(5): 625–632. DOI: 10.11883/1001-1455(2015)05-0625-08.
    [4] 陈川琳, 黄陈磊, 许辉, 等. 小口径步枪弹头后效期运动特性试验与数值研究 [J]. 兵工学报, 2019, 40(2): 265–275. DOI: 10.3969/j.issn.1000-1093.2019.02.006.
    CHEN C L, HUANG C L, XU H, et al. Experimental and numerical research on motion characteristics of a small caliber bullet in muzzle flows [J]. Acta Armamentarii, 2019, 40(2): 265–275. DOI: 10.3969/j.issn.1000-1093.2019.02.006.
    [5] 张欣尉, 余永刚. 水下发射对机枪膛口温度场影响的数值分析 [J]. 含能材料, 2017, 25(11): 932–938. DOI: 10.11943/j.issn.1006-9941.2017.11.008.
    ZHANG X W, YU Y G. Numerical analysis for the effect of underwater launch on the temperature field of machine gun muzzle [J]. Chinese Journal of Energetic Materials, 2017, 25(11): 932–938. DOI: 10.11943/j.issn.1006-9941.2017.11.008.
    [6] 张欣尉, 余永刚, 莽珊珊. 装药参数对水下机枪密封式膛口流场影响的数值分析 [J]. 兵工学报, 2018, 39(1): 18–27. DOI: 10.3969/j.issn.1000-1093.2018.01.002.
    ZHANG X W, YU Y G, MANG S S. Numerical analysis of influence of charge parameters on flow field around sealed muzzle of underwater machine gun [J]. Acta Armamentarii, 2018, 39(1): 18–27. DOI: 10.3969/j.issn.1000-1093.2018.01.002.
    [7] HU Z T, YU Y G. Expansion characteristics of multiple wall jets in cylindrical observation chamber [J]. Applied Thermal Engineering, 2017, 113: 1396–1409. DOI: 10.1016/j.applthermaleng.2016.11.140.
    [8] ZHAO J J, YU Y G. Flow structure of conical distributed multiple gas jets injected into a water chamber [J]. Journal of Mechanical Science and Technology, 2017, 31(4): 1683–1691. DOI: 10.1007/s12206-017-0316-9.
    [9] ZHOU L L, YU Y G. Study on interaction characteristics between multi gas jets and water during the underwater launching process [J]. Experimental Thermal and Fluid Science, 2017, 83: 200–206. DOI: 10.1016/j.expthermflusci.2017.01.007.
    [10] 郝宗睿, 王乐勤, 吴大转. 水下喷气推进高速射流场数值研究 [J]. 浙江大学学报(工学版), 2010, 44(2): 408–412. DOI: 10.3785/j.issn.1008-973X.2010.02.036.
    HAO Z R, WANG L Q, WU D Z. Numerical simulation of high-speed jet flow field of underwater jet propulsion craft [J]. Journal of Zhejiang University (Engineering Science), 2010, 44(2): 408–412. DOI: 10.3785/j.issn.1008-973X.2010.02.036.
    [11] 唐云龙, 李世鹏, 谢侃, 等. 有相变的水下超音速燃气射流数值模拟 [J]. 哈尔滨工程大学学报, 2016, 37(9): 1237–1243. DOI: 10.11990/jheu.201506010.
    TANG Y L, LI S P, XIE K, et al. Numerical simulation of underwater supersonic gas jets with phase transitions [J]. Journal of Harbin Engineering University, 2016, 37(9): 1237–1243. DOI: 10.11990/jheu.201506010.
    [12] 张焕好, 郭则庆, 王瑞琦, 等. 水下超声速气体射流的初始流动特性研究 [J]. 振动与冲击, 2019, 38(6): 88–93, 131. DOI: 10.13465/j.cnki.jvs.2019.06.013.
    ZHANG H H, GUO Z Q, WANG R Q, et al. Initial flow characteristics of an underwater supersonic gas jet [J]. Journal of Vibration and Shock, 2019, 38(6): 88–93, 131. DOI: 10.13465/j.cnki.jvs.2019.06.013.
  • [1] 吴伟许厚谦王亮薛锐 . 含化学反应膛口流场的无网格数值模拟. 爆炸与冲击, 2015, 35(5): 625-632. doi: 10.11883/1001-1455(2015)05-0625-08
    [2] 赵丽俊焦志刚李晓婕黄晓杰赵东志朱小平唐辉王昌富田振喻薪宁 . 预制破片侵彻均质装甲钢的极限穿透速度. 爆炸与冲击, 2018, 38(1): 183-190. doi: 10.11883/bzycj-2016-0116
    [3] 赵成功王聪魏英杰张孝石 . 细长体水下运动空化流场及弹道特性实验. 爆炸与冲击, 2017, 37(3): 439-446. doi: 10.11883/1001-1455(2017)03-0439-08
    [4] 尤国钊魏琪 . 膛口流场的数值模拟. 爆炸与冲击, 1989, 9(3): 254-260.
    [5] 刘东尧周彦煌 . 液体发射药电热化学发射内弹道一维两相流模型及数值模拟. 爆炸与冲击, 2002, 22(2): 158-162.
    [6] 贠来峰芮筱亭侯日升何斌 . 考虑两相流内弹道的自行火炮发射动力学计算. 爆炸与冲击, 2007, 27(1): 12-17. doi: 10.11883/1001-1455(2007)01-0012-06
    [7] 贠来峰芮筱亭王国平陈涛 . DCD格式在破碎发射药床两相流内弹道计算中的应用. 爆炸与冲击, 2010, 30(3): 295-300. doi: 10.11883/1001-1455(2010)03-0295-06
    [8] 黄守龙王红马毓义 . 含内插管扩张型抑制装置的膛口流场数值研究. 爆炸与冲击, 1998, 18(2): 155-161.
    [9] 余永刚方谋鑫汪庆永 . 膛口流场一维简化非傅里叶传热模型. 爆炸与冲击, 1994, 14(4): 359-362.
    [10] 栗保明金志明 . 再生式液体发射药火炮贮液室流体响应特性的分析. 爆炸与冲击, 1995, 15(4): 356-362.
    [11] 戴荣栗保明卢顺祥 . 固体发射药等离子体点火过程及弹道效应分析. 爆炸与冲击, 1999, 19(1): 84-89.
    [12] 夏正友张河陈家安 . 一种非火药驱动气体炮内弹道模型及发射诸元协调. 爆炸与冲击, 1999, 19(2): 146-150.
    [13] 张磊阮文俊王浩王鹏新 . 轻型火箭发射噪声场的分布特性. 爆炸与冲击, 2015, 35(5): 735-740. doi: 10.11883/1001-1455(2015)05-0735-06
    [14] 张博孜王浩王珊珊 . 含扰流装置的中心燃气式抛撒内弹道过程的数值模拟. 爆炸与冲击, 2015, 35(2): 208-214. doi: 10.11883/1001-1455(2015)02-0208-07
    [15] 朱珠罗松卢丙举于勇 . 旋转射弹高速倾斜入水多相流场与弹道数值模拟. 爆炸与冲击, 2019, 39(11): 113901-1-113901-11. doi: 10.11883/bzycj-2018-0315
    [16] 周彦煌张兆钧官汉章栗保明 . 一种整装式液体发射药燃烧控制技术. 爆炸与冲击, 1995, 15(3): 275-281.
    [17] 夏倩倩张文毅祁兵杨风波乐贵高 . 湿式发射对同心筒瞬态热冲击的影响机理. 爆炸与冲击, 2018, 38(1): 37-49. doi: 10.11883/bzycj-2016-0141
    [18] 陈二云赵改平杨爱玲 . 环形激波聚焦流场特性的数值研究. 爆炸与冲击, 2012, 32(3): 291-296. doi: 10.11883/1001-1455(2012)03-0291-06
    [19] 杨艺何学秋刘建章王从银 . 瓦斯爆燃火焰内部流场分形特性研究. 爆炸与冲击, 2004, 24(1): 30-36.
    [20] 李润之司荣军 . 低温环境下甲烷爆炸流场特性模拟. 爆炸与冲击, 2015, 35(6): 901-906. doi: 10.11883/1001-1455(2015)06-0901-06
  • 加载中
图(13)表(2)
计量
  • 文章访问数:  1564
  • HTML全文浏览量:  653
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-25
  • 录用日期:  2020-03-02
  • 网络出版日期:  2020-09-21
  • 刊出日期:  2020-10-05

弹道枪不同水深下全淹没式发射膛口流场的数值分析

    作者简介:张京辉(1995- ),男,博士研究生,zjhsg@vip.qq.com
    通讯作者: 余永刚, yygnjust801@163.com
  • 南京理工大学能源与动力工程学院,江苏 南京 210094

摘要: 为了解弹道枪水下全淹没发射时,水深对膛口流场演化特性的影响,建立了二维轴对称非稳态膛口流场模型。采用流体体积函数多相流模型、标准k-ε湍流模型和Schnerr-Sauer空化模型,结合动网格及用户自定义函数技术,对水下全淹没发射膛口流场演变全过程进行了数值模拟。搭建了弹道枪水下可视化射击实验平台,对12.7 mm口径弹道枪在水中全淹没式发射时膛口流场演化过程进行了观测,并验证了数值模型的合理性。在此基础上,对比了不同水深下(h=1~100 m)膛口流场的演化特性。通过对比发现:在不同水深条件下,在膛口流场影响范围内,弹丸膛外行程随时间的变化均满足指数函数规律;水越深,膛口流场典型波系结构形成所需时间越长,且燃气在膛口轴向马赫盘处的温度和压力峰值越低,压力振荡幅度也越小,更快趋于平稳,但在径向上,水越深,压力振荡持续时间越长。

English Abstract

  • 膛口流场一般指枪炮发射时,膛内气体在膛口外膨胀而形成的随时间变化的气流区域[1],是非定常、带有强激波间断的复杂流场。众多学者针对枪炮在空气中发射的膛口流场进行了大量实验和数值研究。姜孝海等[2]运用有限体积法,对空气中膛口流场形成的动力学过程进行了数值模拟。吴伟等[3]基于无网格方法,对包含大位移运动边界和非平衡化学反应的膛口流场进行了数值模拟。陈川琳等[4]结合实验与数值模拟,研究了弹丸在膛口流场中的受力情况和运动规律。

    而在水下枪炮发射过程中,由于水的物性特点与空气大不相同,且不同水深下发射的环境压力也会不同,这不仅使得水下航行体的运动阻力极大,而且会对水下枪炮发射所形成的膛口流场的演化特性产生复杂的影响,这对枪炮等身管武器在水下进行发射提出了极大的挑战。在水下膛口流场方面,张欣尉等[5-6]运用Fluent软件,针对机枪在空气中和水下密封式发射所形成的膛口流场进行了数值对比分析。水下膛口流场涉及气液两相流相互作用问题,围绕燃气射流与水的作用,许多学者开展了实验与数值模拟。Hu等[7]建立了燃气射流在受限空间内扩展的三维非稳态数理模型,对贴壁燃气射流在充液圆管内的扩展特性进行了数值模拟。Zhao等[8]建立了多股燃气射流数值求解模型,对锥形分布燃气射流的干涉、汇聚特性进行了分析。在前人的基础上,Zhou等[9]对静态和准动态条件下气幕式发射的内弹道及气液相互作用过程进行了数值模拟。而在水下超声速射流方面,郝宗睿等[10]采用VOF (volume of fluid)两相流模型结合动网格方法对系留状态和运动状态下的燃气喷射流场进行数值模拟,分析了两种状态下燃气喷射形成的复杂流场结构。唐云龙等[11]设计了一种解决水下超音速燃气射流复杂相变过程的计算模型,并对典型工况的水下高温、高速燃气射流问题中相变过程进行了数值分析。张焕好等[12]建立了二维轴对称水下超声速气体射流的数值计算模型,并进行了相关的数值模拟,得到了水下超声速气体射流的初始流动结构。

    已有研究主要针对空气中膛口流场和水下燃气射流问题,但对水下全淹没发射(即身管内充满水),尤其是不同发射水深条件下的膛口流场演化特性的研究较少。本文中,建立二维轴对称非稳态膛口流场模型,对水下全淹没发射膛口流场演变全过程进行数值模拟,并搭建水下可视化射击实验平台,利用高速录像观测12.7 mm口径弹道枪在水中全淹没式发射时膛口流场演变的过程。在此基础上,针对12.7 mm口径弹道枪在3种水深下全淹没发射膛口流场演化过程进行数值模拟,以期研究结果对深入了解水下膛口流场特性以及降低膛口气流危害有所帮助,同时对新型水下枪炮设计有一定的参考价值。

    • 根据弹道枪水下全淹没式发射特点,采用下列物理模型:

      (1)火药燃烧遵循几何燃烧定律。药粒均在平均压力下燃烧,且遵循指数燃速定律。单位质量火药燃烧所放出的热量及生成的燃气温度均为定值,在以后的膨胀做功过程中,不考虑燃气组分变化,火药力fp、余容αp及比热比kp等均视为常数。

      (2)枪管内膛表面的热散失用减小火药力fp或增大比热比kp的方法间接修正,用系数φ来考虑其他的次要功。

      (3)考虑到在膛口流场中气相密度远小于液相密度,重力场影响很小,以二维效应为主,因此本文中采用二维轴对称模型对膛口流场进行研究,湍流模型采用标准k-ε模型。

      (4)将火药燃气看作理想气体,不考虑燃气在膛口的二次燃烧,忽略体积力的影响。

      (5)假设水为不可压缩相,密度取998.2 kg/m3。因膛口燃气射流与水的作用时间短暂,忽略燃气射流对水的汽化作用。

      (6)对于膛口以及弹丸表面附近水的空化现象,采用Schnerr-Sauer空化模型进行计算。

    • 针对上述物理模型建立下列数学模型:

      (1)连续性方程:

      $\frac{{\partial \left( {{\alpha _q}{\rho _q}} \right)}}{{\partial t}} + \nabla \times \left( {{\alpha _q}{\rho _q}{{ U}_q}} \right) = {S_{\rm{m}}}$

      式中:q=1,2,3分别代表燃气相、液相和水蒸气相,ρq为对应项的密度;α1α2α3为各相的体积分数,且α1+α2+α3=1;Uq为流体速度;${S_{\rm{m}}}$为单位体积内气液两相的质量流量,这一项在水与蒸汽相发生质量输运时才有,当各相之间不发生质量输运时,此项为零。

      (2)动量守恒方程:

      $\frac{{\partial \left( {{\alpha _q}{\rho _q}{{ U}_q}} \right)}}{{\partial t}} + \nabla \times \left( {{\alpha _q}\left( {{\rho _q}{{ U}_q} \times {{ U}_q}} \right)} \right) = - {\alpha _q}\nabla {p_q} + {\mu _q}{\nabla ^2}{{ U}_q} + {{ S}_{\rm{p}}}$

      式中:混合物密度$\rho = \displaystyle\sum\limits_i^{\rm{3}} {{\alpha _q}{\rho _q}} $pq为对应项的压力;${{ S}_{\rm{p}}}$为汽液两相间质量输运引起的动量变化,此项在水与蒸汽相发生质量输运时才有,当各相之间不发生质量输运时,此项为零。

      (3)能量守恒方程:

      $\frac{\partial }{{\partial t}}\left( {\rho E} \right) + \nabla \times \left[ {{{ U}_q}\left( {\rho E + p} \right)} \right] = \nabla \times \left( {{k_{\rm{eff}}}\nabla T} \right) + {S_{\rm{e}}}$

      $E = {{\sum\limits_q^{\rm{3}} {{\alpha _q}{\rho _q}{E_q}} } \Bigg/ {\left( {\sum\limits_q^{\rm{3}} {{\alpha _q}{\rho _q}} } \right)}}$

      $T = {{\sum\limits_q^{\rm{3}} {{\alpha _q}{\rho _q}{T_q}} } \Bigg/ {\left( {\sum\limits_q^{\rm{3}} {{\alpha _q}{\rho _q}} } \right)}}$

      式中:$E$为平均能量;p为混合相压力;$T$为平均温度;keff为有效热传导率;Se为能量方程源项,即空化引起的流场能量的改变,其大小等于气液之间质量交换量与汽化潜热的乘积。

      (4)气体状态方程:

      $p = \rho RT$

      式中:R=319.6 J/(kg·K),为火药气体常数。

      (5)$k{\text{-}}\varepsilon $湍流方程:

      $\frac{\partial }{{\partial t}}\left( {\rho k} \right) + \frac{\partial }{{\partial {x_i}}}\left( {\rho k{u_i}} \right) = \frac{\partial }{{\partial {x_j}}}\left[ {\left( {\mu + \frac{{{\mu _t}}}{{{\sigma _{\rm{k}}}}}} \right)\frac{{\partial k}}{{\partial {x_j}}}} \right] - \rho \overline {u_i^{'}u_j^{'}} \frac{{\partial {u_j}}}{{\partial {x_i}}} - \rho \varepsilon $

      $\frac{\partial }{{\partial t}}\left( {\rho \varepsilon } \right) + \frac{\partial }{{\partial {x_i}}}\left( {\rho \varepsilon {u_i}} \right) = \frac{\partial }{{\partial {x_j}}}\left[ {\left( {\mu + \frac{{{\mu _t}}}{{{\sigma _\varepsilon }}}} \right)\frac{{\partial \varepsilon }}{{\partial {x_j}}}} \right] - {C_{\varepsilon 2}}\rho \frac{{{\varepsilon ^2}}}{k} + {C_{\varepsilon 1}}\frac{\varepsilon }{k}\rho \overline {u_i^{'}u_j^{'}} \frac{{\partial {u_j}}}{{\partial {x_i}}}$

      式中:$ k $为湍动能,$\varepsilon $为耗散率;ij为自由指标,$u_i^{}$$u_j^{}$为速度矢量,μ为黏度系数,$\overline {u_i^{'}u_j^{'}} $为雷诺压力项;常数σk=1.0为湍动能对应的普朗特数;σε=1.3为耗散率对应的普朗特数;μt=Cμk2/ε表示湍流黏性系数,常数Cε1=1.44,Cε2=1.92和Cμ=0.08为经验系数。

      (6)Schnerr-Sauer空化模型:

      ${\varGamma _{\rm{e}}} = \frac{{{\rho _{\rm{v}}}{\rho _{\rm{l}}}}}{\rho }{\alpha _{\rm{v}}}(1 - {\alpha _{\rm{v}}})\frac{3}{{{R_{\rm{b}}}}}\sqrt {\frac{{2\left( {{p_{\rm{v}}} - p} \right)}}{{{\rho _{\rm{l}}}}}} \quad\quad\quad {p_{\rm{v}}} {\text{>}} p$

      ${\varGamma _{\rm{c}}} = \frac{{{\rho _{\rm{v}}}{\rho _{\rm{l}}}}}{\rho }{\alpha _{\rm{v}}}(1 - {\alpha _{\rm{v}}})\frac{3}{{{R_{\rm{b}}}}}\sqrt {\frac{2}{3}\frac{{\left( {p - {p_{\rm{v}}}} \right)}}{{{\rho _{\rm{l}}}}}} \quad\quad\quad {p_{\rm{v}}} {\text{<}} p$

      式中:${\varGamma _{\rm{e}}}$${\varGamma _{\rm{c}}}$分别为气泡生长和溃灭时质量传递源项;${\alpha _{\rm{v}}}$为空泡的体积分数,c为单位体积内的气泡数,${\rho _{\rm{v}}}$${\rho _{\rm{l}}}$分别为汽相和液相的密度,$\rho $为混合相的密度,${p_{\rm{v}}}$为蒸汽分压力,${R_{\rm{b}}}$为气泡的半径。气泡半径的表达式为:

      ${R_{\rm{b}}} = {\left( {\frac{{{\alpha _{\rm{v}}}}}{{1 - {\alpha _{\rm{v}}}}}\cdot\frac{3}{{4{\text{π}} }}\cdot\frac{1}{c}} \right)^{\frac{1}{3}}}$

      膛口流场的计算中需要耦合内弹道方程组,内弹道方程如下:

      (1)形状函数:

      $\psi = \left\{ {\begin{array}{*{20}{l}} {\chi Z\left( {1 + \lambda Z + {\mu _{\rm{c}}}{Z^2}} \right)}&{ Z {\text{<}} 1} \\ {{\chi _{\rm{s}}}\dfrac{Z}{{{Z_{\rm{k}}}}}\left( {1 + {\lambda _{\rm{s}}}\dfrac{Z}{{{Z_{\rm{k}}}}}} \right)}&{{\rm{ 1}} {\text{≤}} Z {\text{<}} {Z_{\rm{k}}}} \\ 1&{{\rm{ }}Z {\text{≥}} {Z_{\rm{k}}}} \end{array}} \right.$

      式中:ψ为火药燃烧质量百分比,χλμc均为火药形状函数,Z为火药燃烧相对厚度,Zk为当颗粒燃尽时燃烧相对厚度,${\chi _{\rm{s}}}$${\lambda _{\rm{s}}}$为在火药表面缩小阶段的形状函数。

      (2)燃速方程:

      $\frac{{{\rm{d}}Z}}{{{\rm{d}}t}} = \frac{{{u_1}{p^n}}}{{{e_1}}}$

      式中:${u_1}$为火药燃速系数,n为火药燃速指数,${e_1}$为火药半弧厚,p为膛内火药燃气平均压力,t为时间。

      (3)弹丸运动方程:

      $\int_0^A {({p_{\rm{d}}} - {p_{\rm{f}}}} ){\rm{d}}A = \varphi m\frac{{{\rm{d}}v}}{{{\rm{d}}t}}$

      式中:${p_{\rm{d}}}$${p_{\rm{f}}}$分别为弹底和弹前压力,A为弹丸横截面积,φ为次要功系数,m为弹丸质量,v为弹丸运动速度。

      (4)内弹道基本方程:

      $\frac{{Ap({l_\psi } + x)}}{\theta } = \frac{{fw\psi }}{\theta } - \frac{\varphi }{2}m{v^2} - \int_0^x {\int_0^A {{p_{\rm{f}}}} {\rm{d}}A} {\rm{d}}x$

      ${l_\psi } = {l_0}\left[ {1 - \frac{\varDelta }{{{\rho _{\rm{p}}}}}(1 - \psi ) - {\alpha _{\rm{p}}}\varDelta \psi } \right]$

      式中:${l_\psi }$为药室自由容积缩径长,绝热指数$\theta = {k_{\rm{p}}} - 1$w为装药质量,x为弹丸行程,Δ为装填密度,αp为单位质量气体分子体积有关的修正量,称为余容。

      (5)弹丸速度与行程关系式:

      $\frac{{{\rm{d}}x}}{{{\rm{d}}t}} = v$

      方程(12)~(17)组成内弹道方程组。

    • 本文中借助FLUENT软件开展数值模拟,结合用户自定义函数(user defined function,UDF)技术,耦合内弹道方程组求解得到弹丸速度、位移和膛压;采用流体体积函数(volume of fluid,VOF)多相流模型来描述气液相互作用;利用PRESTO!方法对压力项离散;动量和能量的离散采用一阶迎风格式;采用PISO算法对压力与速度进行耦合;计算采用的时间步长控制在0.1 µs内。

    • 图1所示,计算域分为身管区域和膛口流场区域两个部分,其中身管区域为长1 m、半径为6.35 mm的圆柱形区域,流场区域为长1 m、半径0.3 m的圆柱形区域。以结构化网格为主,对膛口附近的网格进行局部加密,最小网格尺寸为0.3 mm×0.3 mm。采用了动网格技术中的层铺法来实现弹丸的运动,随着弹丸向前运动,当靠近弹底的网格长度被拉长至0.6 mm时,网格会分裂成两个网格,同理,当靠近弹头的网格长度被压缩至0.2 mm时,该网格会与相邻的网格合并。

      图  1  计算域及边界条件

      Figure 1.  Calculation domain and boundary conditions

      本文计算的边界条件为:燃烧室为压力入口,其初始参数由内弹道方程组求解得出;身管壁面及弹丸为固壁边界;计算域外边界为压力出口,压力为环境压力。将弹前身管区域和流场区域初始化为充满水,温度和压力初始化为环境条件,即温度为300 K,压力根据不同水深分别设置为111.5、608.0和1 114.6 kPa。

    • 为尽量避免网格粗细程度对计算结果产生较大影响,分别采用20万、15万和10万的网格数(N)进行试算,选择弹丸出膛后1 ms时刻膛口中心轴向速度分布作为网格无关性验证的特征参数,无关性验证结果如图2所示。相对于20万网格数的计算结果,15万网格数的最大相对误差仅为6%,而10万网格数的最大相对误差达到30.6%。综合考虑计算效率和计算精度,现选择网格数为15万的网格进行计算。

      图  2  不同网格数下膛口中心速度沿轴向的变化(t=1 ms)

      Figure 2.  Velocity of the muzzle center varying along the axial direction for different grid quantities at t=1 ms

    • 为观察水下全淹没式膛口流场中气、汽、液三相的掺混过程以及验证数值模型的合理性,以12.7 mm口径的弹道枪为平台,搭建了全淹没式发射的可视化实验测试系统,如图3所示。实验系统由水箱、12.7 mm口径弹道枪、水泥平台以及脉冲电点火器组成,其中水箱上开有射击孔,供身管伸入,水箱侧面开有观察窗。采用FASTCAM-Ultima APX高速摄像机和计算机进行流场演化过程的观察和记录,其中高速摄像机最大分辨率为1 024 ×512像素,频率为4 000 s−1

      图  3  实验装置示意图

      Figure 3.  Schematic diagram of the experimental device

      在实验准备阶段,将弹丸填充在枪膛底部,通过弹带进行密封以阻止水进入弹后燃烧室;为保证全淹没发射环境,将身管以及水箱内注满水;调整高速摄像拍摄的焦距与角度,使其能透过观察窗清楚的拍摄到膛口流场区域;当仪器调整好后,将称量好的火药装入燃烧室,装好电底火,并将其连接到脉冲电点火器上。所有装置就绪后,用脉冲电点火器放电进行点火。实验后,用计算机将高速摄像拍摄的结果进行截取和后处理。

      本次模型验证所选取的实验工况水深为0.3 m,对应的数值模拟工况水深也为0.3 m。采用的实心弹丸质量为45 g,装药采用6 g的4/7单基火药,实验测得弹丸初速为264.2 m/s。实验阴影照片与数值模拟相图的对比情况如图4所示,将弹丸出膛口时刻定义为t=0时刻。

      图  4  实验照片阴影图与模拟相图对比情况

      Figure 4.  Comparison of experimental shadow photos and simulated phase diagrams

      图4可见,在全淹没发射条件下,身管内与外界环境都充满了水。图4(a)(b)对应点火后,弹丸在膛内如同活塞般不断推动弹前水柱的过程,其中图4(a)对应弹丸出膛前0.8 ms时刻。水柱在弹丸推动下不断加速,使得膛口附近的局部压力低于水的饱和蒸汽压力,膛口附近的水不断汽化为水蒸汽,在膛口逐渐堆积成一个蒸汽团。图4(b)对应弹丸即将出膛口瞬间,弹丸头部被侧面的蒸汽团完全包裹。图4(c)(e)对应弹丸逐渐飞离膛口的过程。在图4(c)中,高压火药燃气由膛内喷出,与蒸汽团迅速掺混,初步形成了泰勒空腔,同时提高了膛口附近的压力,使得蒸汽团迅速溃灭。随着弹丸运动到射流头部,对射流在轴向上的约束减弱,使射流在轴向上迅速扩展,而在径向上扩展仍受周围水介质的限制,因此头部呈锥形,同时气液之间的强烈湍流掺混引起的Kelvin-Helmholtz不稳定效应使得射流边界变得不规则。弹丸穿过泰勒空腔后相对水介质做高速运动,使得弹丸头部表面的水的压强降低,初步形成自然空泡。在图4(d)(e)中弹丸已完全超越射流,弹丸出膛1.2 ms后,随着弹丸继续向前飞行,弹丸表面的空泡长度逐渐变长,形成超空泡,包裹着弹丸前进,同时在弹丸尾部会留下细长的气柱。

      图4中可以看出实验阴影照片与数值模拟相图中的膛口流场的外部轮廓基本吻合,为进一步验证数值模型的合理性,做出了射流轴向、径向最大位移对比图,如图5所示。高速录像在进行拍摄时,由于光路传播原因,照片上水箱前后透明观察窗固定螺钉的间距为等比例缩小。在实验中,枪管位于水箱对称面位置,根据梯形中位线特征,在进行膛口数据处理时,前后观察窗的螺钉间距可以直接从照片中量取,两者的平均值与实际螺钉间距的比值即为量取长度与实际长度的比值,再结合照片中量取的射流轴向、径向最大位移从而得到实际射流的最大位移。为了减小误差,对每种工况实验数据进行了多次测量求平均值。由图5可知,模拟结果中不同时刻燃气射流的轴向、径向最大位移与实验测量值吻合较好,最大相对误差为6.4%,说明本文中采用的数值模型基本合理。

      图  5  射流轴向、径向最大位移对比

      Figure 5.  Comparison of the maximum axial and radial displacements of the jet flow between experimental and simulated results

    • 为研究发射水深对全淹没式发射膛口流场的影响,针对12.7 mm口径弹道枪在3种水深h下全淹没发射膛口流场演化过程进行了数值模拟和对比分析。3种工况的初始参数如表1所示,表中l为身管长度,v0为弹丸初速,pk0为膛口压力。图6给出了不同水深条件下弹丸的膛外行程对比曲线。

      h/ml/mv0/(m·s−1)pk0/MPa
      1123014.5
      50122015.5
      100120720.5

      表 1  膛口初始参数

      Table 1.  Initial parameters for the muzzle

      图  6  不同水深下弹丸膛外行程

      Figure 6.  Displacement-time curves of the projectile at different water depths

      表1可见,随着水深增加,弹丸初速降低明显,膛口压力也随之升高。由图6可见,水深越深,弹丸速度越慢。对膛口流场影响范围内,弹丸行程随时间变化特性进行拟合,发现其在不同水深条件下均满足指数函数规律,即:

      ${x_t} = {x_0} - {x_1}{\rm{ }}{{\rm{e}}^{ - t/{t_1}}}$

      式中:xt为弹丸膛外行程,m;t为时间,ms;x0x1t1为弹丸膛外行程随时间变化的拟合参数,如表2所示。由表2可见,x0x1t1都会随着水深增大而减小。

      h/mx0/mx1/mt1/ms
      11.091.094.55
      500.840.843.55
      1000.790.793.42

      表 2  拟合参数

      Table 2.  Fitting parameters

      全淹没发射条件下,由于膛内的火药燃气处于高温高压状态,当燃气从膛口喷出后,会迅速膨胀与高密度水介质以及弹丸相互作用,从而形成复杂的波系结构。为了便于对不同发射环境条件下的膛口流场流谱结构进行对比分析,图7给出了空气中枪炮发射时的膛口流场流谱示意图[1], 图8给出了全淹没发射时膛口流场流谱示意图。为研究水深对膛口流场演化特性的影响,图911分别给出了3种水深条件下的膛口中心剖面的马赫数云图和等值线图。

      图  7  空气中发射时典型膛口流场流谱示意图

      Figure 7.  Schematic diagram of the typical muzzle flow field during air launch

      图  8  全淹没发射时膛口流场流谱示意图

      Figure 8.  Schematic diagram of the muzzle flow field during submerged launch

      图  9  h=1 m时膛口中心剖面马赫数云图与等值线图

      Figure 9.  Mach number cloud map and contour map at h=1 m

      图  10  h=50 m时膛口中心剖面马赫数云图和等值线图

      Figure 10.  Mach number cloud map and contour map at h=50 m

      图  11  h=100 m时膛口中心剖面马赫数云图和等值线图

      Figure 11.  Mach number cloud map and contour map at h=100 m

      图7可见,空气中膛口流场将射流理论边界包围的整个流场区域划分为了5个区域,分别为:(1)区,即瓶状激波内的自由膨胀区,气流主要在此瓶状区内膨胀,压力剧降、速度激增;(2)区,即相交激波与射流边界之间的超声速区;(3)区,即马赫盘下游的亚声速区,气流经过正激波后速度剧降、压力和温度陡增;(4)区,经过两次斜激波后,流动情况变得复杂,压力与(3)区相同,但为超声速区,两区之间形成类似于拉瓦尔喷管形状的切向间断面;(5)区,即湍流混合过渡区。

      图911可见,水下全淹没式发射时,由于燃气射喷出膛口后迅速膨胀,膨胀波系在气液边界以及弹底发生反射形成激波,约在t=0.4 ms时刻,演化成了如图8所示的波系结构。类似于空气中膛口流场典型结构,入射激波、马赫盘以及反射激波交汇形成了三波点,通过正激波后的亚声速气流与经过斜激波后速度仍为超声速的气流之间有速度差,因此在三波点沿流向形成切向间断。瓶状激波虽然已经形成,但受高密度水介质的影响,流场波系结构的分布特点存在一些不同:(1)区(自由膨胀区)受气液边界的影响,形状尚不饱满,入射斜激波在径向上向内收缩,马赫盘在轴向上向外突出且直径较小;受气液边界和(1)区影响,(2)区被挤压呈狭长状被未能充分发展至将(1)区包裹起来;(3)区也受高密度水介质的压缩影响,在径向上大大被压缩。在后续的发展中(t>0.4 ms),随着弹丸完全脱离泰勒空腔,对膛口近场的波系结构不再有影响,同时膛内后续气体不断喷出,使得泰勒空腔在径向上扩展明显,对瓶状激波在径向上的限制减小。整个瓶状激波系逐渐饱满,最终形成如图7所示的膛口流场典型波系结构。

      图911的对比中可以看出,在3种水深下,弹丸出膛后0.4 ms,包括入射激波、反射激波、马赫盘和三波点等结构的瓶状激波都已经形成。当水深为1 m时,在t=0.4~1.4 ms期间,瓶状激波在径向上趋于饱满,马赫盘直径也相应增大,最终在t=1.4 ms形成与如图所示的膛口流场典型波系结构。而在水深为50 m时,在t=0.4~1.4 ms期间,瓶状激波在径向上有一定扩展,马赫盘沿轴向逐渐向下游移动,使得瓶状激波的“瓶口”仍向外突出。在t=1.4 ms时,马赫盘已经到达距膛口最大位移并逐渐向膛口方向移动,在t=1.8 ms时,瓶状激波才趋于饱满,形成膛口流场典型波系结构。而在水深为100 m的条件下,在t=2.4 ms时,流场内才形成膛口流场典型波系结构。

      这说明相同装药条件下,水深越大,对应的弹丸初速越低,脱离膛口流场区域所需的时间越长,对流场内的波系结构演化的限制时间更长,同时也由于水深处的环境压力较大,使得水深越大,瓶状激波在径向上需要更长时间趋于饱满,形成膛口流场典型波系结构。

      为研究流场内压力变化特性,图12给出了水深为50 m的条件下膛口中心剖面的压力云图,图13则给出了不同水深下,压力沿膛口中心轴线的分布曲线和瓶状激波内(x=1.02 m处)压力沿径向的分布曲线。

      图  12  h=50 m是膛口中心剖面压力云图

      Figure 12.  Pressure cloud map of the muzzle center section at h=50 m

      图  13  膛口压力分布曲线

      Figure 13.  Pressure distribution curves of muzzle

      图12可见,弹丸出膛后,膛内的火药燃气喷出膛口,迅速膨胀,在t=0.4 ms已经形成了包含瓶状激波在内的复杂激波结构,从而影响了流场内的压力分布。瓶状激波内是燃气射流的主要膨胀区,射流在瓶区内迅速膨胀,使得瓶区内的压力剧降。而在马赫盘的下游,由于射流受到周围高密度水介质与弹底的压缩作用,使得压力骤升,形成高压区域。随后由于燃气自身被压缩后的膨胀作用,以及受到的压缩作用,使得马赫盘后高压区域与低压区域交替出现。在t=1.4 ms后,由于弹丸远离膛口区域,以及后续燃气喷出,使得马赫盘下游的压力趋于平均,不再有突变,而自由膨胀区内的压力仍然较低,瓶状激波结构仍清晰可见。

      由图13(a)(b)可见,在水深为1 m时,在弹丸出膛0.4 ms时刻,此时弹丸仍处在膛口近场区域内,燃气从膛口喷出后迅速膨胀,在轴向上压力降至约0.2 MPa。随后由于水介质与弹底的压缩作用,压力在马赫盘处分别骤升至1.6 MPa。由于射流自身被压缩后的膨胀作用,以及水介质与弹底的压缩作用,流场内的压力在马赫盘至弹底区间内多次震荡,最终在弹底骤升至约60 MPa。而在径向上,燃气压力在小幅度下降后,于瓶状激波的斜激波处压力骤升至1.3 MPa,随后在斜激波至气液边界区间内震荡,最终缓慢下降。

      随着弹丸逐渐远离膛口(t=1.0 ms),膨胀区内的最低温度与最低压力有所降低,马赫盘处的温度与压力峰值也相应减小。在马赫盘下游,压力在轴向上的震荡幅度有所减弱,但仍十分明显,而在径向上,压力在斜激波处到达压力峰值后呈缓慢下降趋势,不再震荡。弹丸出膛1.4 ms后,马赫盘后压力的震荡幅度进一步减弱,几乎消失,压力呈缓慢下降趋势。

      对比不同发射水深下的压力曲线,可以发现,在轴向上,在弹丸出膛后相同时刻,相对于1 m水深时,50和100 m水深发射时,在轴向上,射流出膛口后在马赫盘处的压力峰值偏小,且马赫盘后压力震荡程度大大减弱,震荡次数也相对减少,在弹丸出膛1.4 ms后,马赫盘后的压力几乎不再震荡,沿程缓慢下降。在径向上,由于环境压力的增大,与1 m水深时不同,50和100 m水深发射时,弹丸出膛1.0 ms后,燃气压力于瓶状激波的斜激波处压力骤升后,在斜激波至气液边界区间内仍有较大幅度的震荡,直到t=1.4 ms,压力才趋于平稳,几乎不再震荡。

      整体而言,在不同发射水深条件下,燃气射流出膛后迅速膨胀,压力在轴向与径向上迅速降低,随后分别在马赫盘与斜激波处骤升。随后轴向与径向压力都多次震荡,且震荡幅度随时间逐渐减小,最终呈缓慢下降趋势。当发射水深越深时,在轴向上,燃气在马赫盘处的温度与压力峰值相应越低,压力震荡程度也越小,且更快趋于平稳,而在径向上,压力的震荡持续的时间更长,需要更长时间趋于平稳。

    • 针对弹道枪水下全淹没发射膛口流场问题建立了二维轴对称非稳态模型,借助Fluent软件,结合动网格和用户自定义函数技术,对3种水深下的膛口流场演变过程进行数值模拟,得到如下结论:

      (1)全淹没发射条件下,在弹丸出膛前,身管内水柱在弹丸推动下不断加速,在膛口空化形成蒸汽团;高压火药燃气由膛内喷出后,与蒸汽团迅速掺混,初步形成了泰勒空腔;穿过泰勒空腔后,弹丸头部初步形成自然空泡;随着弹丸继续向前飞行,弹丸表面形成超空泡,包裹着弹丸前进;实验结果与数值模拟结果吻合较好,验证了数值模型的合理性。

      (2)同空气中膛口流场典型波系结构相比,全淹没发射条件下,弹丸出膛后0.4 ms瓶状激波虽已形成,但在径向上尚不饱满;在后续发展中,瓶状激波在径向上逐渐扩展,最终形成膛口流场典型波系结构,且水深越深,形成膛口流场典型波系结构所需时间越长。

      (3)在膛口流场影响范围内,弹丸膛外行程随时间在不同水深条件下均满足指数函数规律,即:$x(t) = {x_0} + {x_1}{{\rm{e}}^{ - t/{t_1}}}$,且水深越深燃气在膛口轴向马赫盘处的温度与压力峰值越低,压力振荡幅度也越小,更快趋于平稳,但在径向上,水深越深,压力振荡的持续时间越长。

参考文献 (12)

目录

    /

    返回文章
    返回