2. 西北工业大学 航海学院, 陕西 西安 710072
超空泡减阻技术可大大降低航行体水下运动所受阻力, 提升水下兵器的打击速度和打击距离。在实际作战中, 为提高打击目标的有效性, 往往需要向水中进行密集射击, 而这一过程中就涉及超空泡射弹异步并联倾斜入水问题, 异步并联入水的超空泡射弹会相互影响, 进而影响射弹运动稳定性, 降低射击精度, 因此对超空泡射弹异步并联倾斜入水问题开展研究对提升入水武器打击效果具有重要意义。
近年来, 众多学者针对倾斜入水问题的研究大多集中在单发射弹入水, 其中李佳川等[1]对弹体模型以不同扰动角速度入水过程进行了弹道仿真; 陈晨等[2]对小型运动体高速倾斜入水过程中流体动力特性、流场结构特性与空泡发展进行了分析; 宋武超等[3]通过数值模拟对不同入水速度和不同入水角度回转体低速倾斜入水过程中的流场特性进行了分析; 路丽睿等[4]开展了不同头型射弹低速倾斜入水的对比试验, 得到了射弹头型对弹道特性的影响规律; 秦杨等[5]通过数值模拟分析了不同入水角下弹体空泡形态发展规律、弹道特性及流体动力特性的变化规律; 施瑶等[6]研究了弹头形状、入水速度和入水角度对射弹入水时流体动力学的影响; 高英杰等[7]对回转体倾斜高速入水的载荷分布、空泡形态、冲击载荷等开展了数值研究; 宋志杰等[8]通过实验和数值方法研究了圆柱形射弹高速倾斜入水下空泡的溃灭及重新产生过程。上述研究针对单发射弹倾斜入水过程中空泡演化、流场特性和弹道特性已进行了较为完善的分析, 然而, 目前对于射弹并联入水问题的研究仍较少, 鹿麟等[9]对超空泡射弹异步并联垂直入水的空泡演化和运动特性展开了研究; 王聪等[10]对运动体并联入水过程的运动特性进行了研究, 并提出运动体相对运动趋势的预测公式; 魏英杰等[11]研究了不同入水速度对圆柱体并联入水运动特性的影响; 余德磊等[12]开展了不同初始参数的回转体高速并联入水模拟计算, 得到了不同参数下的流场及运动特征。但以上研究中运动体大多为简单回转体, 且均为垂直入水, 缺乏对更贴合实际的并联弹丸倾斜入水问题的研究。射弹倾斜入水初期射弹空泡不对称, 且倾斜运动过程中所受水动力与垂直入水也有所不同, 导致超空泡射弹倾斜并联入水过程与垂直并联入水区别较大。
本文基于VOF多相流模型与重叠网格技术建立了超空泡射弹异步并联倾斜入水数值方法, 通过开展单发射弹倾斜入水实验, 对数值方法的有效性进行了验证, 利用此数值方法对超空泡射弹异步并联倾斜入水进行数值模拟, 分析研究了前后两射弹横向间距一定时, 不同入水顺序及不同前后间距下的射弹入水空泡演化及运动特性。
1 数值计算方法 1.1 控制方程超空泡射弹倾斜入水过程涉及水、水蒸气和空气三相均质流体, 本文采用VOF多相流模型对射弹倾斜入水过程进行描述。
混合介质的连续性方程为
|
(1) |
动量守恒方程为
|
(2) |
式中: ρm为混合密度, ρm=αlρl+αvρv+αgρg, ρl, ρv和ρg分别为水、水蒸气和空气的密度, αl, αv和αg分别为三相的体积分数, 三者关系为αl+αv+αg=1。xi, xj为笛卡尔坐标下位移分量; ui, uj为速度分量; p为压力; μm和μt分别为混合介质的动力黏度和湍流动力黏度。
1.2 湍流模型和空化模型本文采用Realizable k-ε湍流模型, 其对于射流、旋转、边界处有较大压力梯度及分离流的计算更加准确, 其表达式为
|
(3) |
|
(4) |
式中：Gk为速度梯度引起的湍动能；Gb为浮力引起的湍动能；YM为脉动扩张在湍流耗散率所占当量；C1ε, C2ε和C3ε为经验系数；σk, σε为k和ε的普朗特常数；Sk和Sε为源项。
本文空化模型采用Schnerr-Sauer空化模型, 水和水蒸气两相间质量传递方程为
|
(5) |
式中：αnuc=5×10-4为不可凝结气体体积分数；RB为气核半径，RB=1×10-6 m；p为远场压力；pv为饱和蒸汽压；Fvap=50, Fcond=0.001。
1.3 数值方法验证通过开展单发射弹倾斜入水实验对本文所用数值方法的准确性和合理性进行验证, 实验装置示意图及实验射弹模型如图 1所示。
|
| 图 1 实验装置及实验射弹示意图 |
实验装置包括水箱、高压发射装置、LED照明灯板、高速摄像机、击发装置及控制电脑, 实验射弹采用弹头为圆锥空化器的超空泡射弹, 射弹直径为D=7.62 mm, 弹长为L=59.4 mm, 空化器锥角为100°, 材质为钢, 质量为12.6 g。数值计算模型中射弹初始位于空气域, 入水初始速度及入水初始角度分别为200 m/s和60°, 与实验射弹初始入水参数相同。图 2和图 3分别为射弹入水空泡形态与速度衰减的数值计算和实验结果对比, 可看出数值计算结果与实验结果基本吻合, 说明本文所选数值方法可有效模拟超空泡射弹倾斜入水。
|
| 图 2 单发射弹倾斜入水实验与仿真空泡形态对比图(单位：mm) |
|
| 图 3 射弹速度变化曲线 |
本文采用射弹模型与前文所述实验射弹完全一致, 图 4为前后两射弹相对位置示意图, 两射弹速度方向与自由液面初始夹角β为60°, 两射弹内侧之间距离始终为D, 前后两弹尖之间距离为S。定义位于虚线左侧的区域为上侧, 位于虚线右侧的区域为下侧; 前发射弹弹尖接触自由液面时t1=0;无量纲量Δ=S/L; 前后两射弹轴向左侧面为背水面, 轴向右侧面为迎水面。
|
| 图 4 射弹相对位置示意图 |
本文所用边界条件及计算域设置如图 5所示, 本文所采用的超空泡射弹为对称回转体, 因此计算域采用1/2计算域, 在对称面设置对称边界条件, 侧面设置固壁边界条件, 顶部及底部设置压力出口条件, 使用FLUENT仿真软件进行数值计算。
|
| 图 5 边界条件及计算域设置示意图 |
本文采用重叠网格技术, 如图 6所示, 整个计算域由背景域和子域构成, 在背景域对射弹运动区域和自由液面进行网格加密, 对距离射弹较远的外层区域使用较疏的网格, 同时在射弹子域的射弹表面添加边界层网格。
|
| 图 6 计算域网格 |
为保证仿真结果独立性, 按照上文网格疏密分布规律, 调整网格尺寸, 开展网格无关性验证, 分别使用网格总数为95万, 125万, 155万的3套网格对Δ=3.0、上侧射弹先入水工况进行数值计算, 图 7为不同计算域网格下的前发射弹速度变化曲线。可以看出, 125万网格和155万网格的速度变化规律基本相同, 故本文采用125万网格的计算域。
|
| 图 7 不同网格数下上侧射弹速度变化曲线 |
图 8~10依次为Δ=0.5, 1.0, 3.0时不同入水顺序下的空泡演化规律, 不同工况下前后两发射弹的初速均为200 m/s, 入射角均为60°。从图中可以看出, 不同前后间距和入水顺序下前后两射弹的空泡演化具有一定差异。对于前发射弹, 当Δ=0.5时, 如图 8a)~8b)所示, 随着后发射弹的入水, 前发射弹附近空泡内侧被挤压收缩, 导致空泡无法完全包裹前发射弹, 导致射弹发生沾湿。但不同入水顺序下, 前发射弹的沾湿情况不同, 当上侧射弹先入水时, 前发射弹在0.23 ms左右内侧尾部即被流体冲击并发生沾湿, 且随着入水深度的增加, 前发射弹迎水面持续沾湿, 射弹最终失稳; 而当下侧射弹先入水时, 前发射弹在0.63 ms左右内侧尾部才与空泡壁接触发生部分沾湿, 随后在空泡内发生一次尾拍, 在1.43 ms左右射弹外侧与空泡壁接触, 沾湿面积持续增大, 最终失稳并对后发射弹运动造成影响。可以发现2种入水顺序下前发射弹运动差异性较大, 将在后文重点分析。观察图 9~10可发现, 随着两射弹前后距离的不断增加, 前发射弹周围空泡受后发射弹的影响不断减弱, 不同入水顺序下前发射弹周围空泡均扩张良好, 射弹运动平稳。
|
| 图 8 Δ=0.5时空泡演化过程(单位：ms) |
|
| 图 9 Δ=1.0时空泡演化过程(单位：ms) |
|
| 图 10 Δ=3.0时空泡演化过程(单位：ms) |
而对于后发射弹, 由图 8~10可观察到, 不同工况下后发射弹空泡沿射弹轴线两侧始终不对称, 自身运动稳定性持续受到影响, 在Δ=0.5时向射弹内侧偏转; Δ=1.0时偏转幅度很小, 射弹运动平稳; Δ=3.0时向射弹外侧偏转。
首先对Δ=0.5时前发射弹运动的差异进行分析。图 11为射弹倾斜入水时空泡壁面附近流体Y方向速度矢量图, 可以发现当射弹倾斜入水时, 空泡壁面上指向Y轴正方向的速度矢量箭头数目更多且数值更大。这是因为当射弹倾斜入水时, 在Y轴正方向上有速度分量Vy, 导致射弹弹尖与流体接触时, 射弹背水面附近的流体会获得更多动能, 此时背水面空泡壁面附近流体速度值更大, 且速度方向指向Y轴正方向, 这也导致射弹倾斜入水时形成的空泡口向Y轴正方向扩张距离更远, 有助于入水前期射弹背水面空气随射弹进入空泡促进背水面空泡扩张。最终导致倾斜入水射弹的迎水面附近流体动能和空泡直径均较小, 而背水面附近流体动能更大, 空泡直径也更大。
|
| 图 11 倾斜入水速度分量示意图 |
当上侧射弹先入水时, 后发射弹的背水面空泡与前发射弹迎水面空泡相互挤压; 而当下侧射弹先入水时, 后发射弹的迎水面空泡与前发射弹的背水面空泡相互挤压, 由于射弹倾斜入水时迎水面和背水面附近流体动能和空泡直径的差异, 当上侧射弹先入水时, 后发射弹对前发射弹空泡的挤压更加剧烈。从图 12也可以观察到, 当上侧射弹先入水时, 两射弹内侧流体Y方向速度矢量均由后发射弹指向Y轴正方向, 前发射弹内侧空泡被严重挤压, 前发射弹内侧尾部此时已发生沾湿; 而当下侧射弹先入水时, 两射弹内侧流体Y方向速度矢量除了有从后发射弹指向Y轴负方向的, 还有从前发射弹指向Y轴正方向的, 前发射弹内侧空泡被挤压程度较小, 空泡此时仍能包裹射弹。同时发现上侧射弹先入水后0.23 ms发生沾湿时, 其还未完全入水形成稳定空泡, 流体的冲击和射弹沾湿会导致其运动稳定性被严重破坏, 无法依靠自身恢复运动, 故其内侧持续沾湿快速失稳; 结合前文知, 当下侧射弹在先入水时, 其在入水后0.63 ms才发生沾湿, 此时其已完全入水生成稳定空泡, 自身稳定性较强, 所以在沾湿后先发生一次尾拍, 之后才持续沾湿并最终失稳。上述两方面原因最终导致了Δ=0.5时2种入水顺序下前发射弹的运动差异。
|
| 图 12 0.23 ms时Δ=0.5工况流体Y方向速度矢量示意图 |
结合图 13~15对上述后发射弹的偏转原因进行分析。图 13a)~13b)为不同入水顺序下后发射弹弹尖附近流线及压力图, 可以发现随着Δ的增大, 不同入水顺序下后发射弹弹尖附近流线及压力图的变化规律一致, 均为弹尖内侧流体由运动受阻变为流速加快, 弹尖高压区也由弹尖内侧逐步向弹尖外侧移动, 且由前文可知不同入水顺序下后发射弹偏转情况同样一致, 故下文仅对上侧射弹先入水时的后发射弹展开分析。
|
| 图 13 后发射弹弹尖附近流线及弹尖压力 |
|
| 图 14 上侧射弹先入水时空泡壁不同特征位置(Δ=0.5, 1.0, 2.0, 3.0)处截面速度矢量图 |
|
| 图 15 上侧射弹先入水时发射弹受力偏转示意图 |
图 14为上侧射弹入水运动4L时空泡壁不同特征位置(Δ=0.5, 1.0, 2.0, 3.0)处截面速度矢量图, 其中A0.5, A1.0, A2.0, A3.0分别表示上侧射弹先入水Δ=0.5, 1.0, 2.0和3.0时的后发射弹。从图 14可以看出, 随着Δ增大, 前发射弹空泡壁面上的速度值逐渐减小。这是由于射弹入水后会将动能不断地传递给周围流体并促进空泡扩张, 随着入水距离的增加, 在内外压差力及流体阻力的作用下, 空泡扩张速度不断衰减, 同时距离弹头越远的位置空泡扩张速度越小。观察图 13a)流线图可看出, 不同特征位置处前发射弹空泡扩张速度不同, 对不同特征位置处的后发射弹产生的影响不同, Δ=0.5处前发射弹空泡扩张速度较大, 后发射弹A0.5弹尖内侧流体受到剧烈挤压形成回流, 随着射弹前后间距逐渐增加, Δ=1.0处前发空泡扩张速度减小, 后发弹尖内侧回流现象随之减弱, 当Δ增大到2.0和3.0时, 由于前发空泡扩张速度已衰减至很小, 后发弹尖内侧排开流体运动未受到前发射弹的阻碍, 同时在前发射弹空泡内部低压区的促进作用下, 流速相对较快。由图 13a)弹尖压力图可发现, 不同工况下后发射弹弹尖压力由于流体运动影响均呈不对称分布, 弹尖内外两侧的压力差形成了促使弹丸失稳旋转的偏转力Fp, Δ=0.5时内侧流速较小导致A0.5弹尖内侧压力较大, 偏转力Fp指向射弹外侧; Δ=1.0时由于阻碍作用减小, A1.0弹尖两侧压力基本对称, 偏转力Fp的大小约等于零; 而Δ=2.0和3.0时内侧流体流速加快, 导致A2.0, A3.0弹尖内侧所受压力减小, 偏转Fp指向射弹内侧。如图 15所示, 偏转力Fp会产生偏转力矩, 最终导致后发射弹在Δ=0.5时向射弹内侧偏转; Δ=1.0时偏转幅度很小, 射弹运动基本平稳; Δ=2.0和3.0时向射弹外侧偏转。
2.2 异步并联斜入水运动特性分析为分析射弹运动过程中的稳定性, 定义无量纲偏转角Kdeg, 表达式为
|
(6) |
式中：θ为射弹运动过程中的偏转角；β为射弹初始入水角度。
横向阻力系数Cy定义为
|
(7) |
式中: Fy为射弹在Y方向所受流体阻力; Fp为射弹所受偏转力; ρw为水的密度; V0为射弹入水速度; A1为射弹侧截面面积;
图 16~17分别为不同入水顺序下前发射弹的无量纲偏转角Kdeg变化曲线、横向阻力系数Cy变化曲线和无量纲偏转角Kdeg随Δ的变化曲线。可以发现对于前发射弹而言, 只有在Δ=0.5时, 无量纲参数Kdeg和横向阻力系数曲线呈现较大的非线性变化趋势, 其他工况下曲线变化与单发斜入水基本一致, 射弹运动平稳。通过对比不同入水顺序下Δ=0.5时前发射弹运动曲线变化规律还可以发现, 当上侧射弹入水发生沾湿后, 其偏转角持续增大, 相应的横向阻力系数也不断增大; 而下侧射弹入水沾湿后, 其偏转角同样持续增大, 但其横向阻力系数变化曲线不是持续增大的, 而是在0.63 ms和1.63 ms左右有2次峰值, 且第二次峰值时横向阻力系数变化的绝对值大于第一次峰值的绝对值。结合前文分析其原因, 当Δ=0.5且上侧射弹先入水时, 前发射弹在入水后0.23 ms左右即发生沾湿, 而此时前发射弹自身稳定性较差, 导致其发生持续沾湿, 横向阻力系数不断增大, 最终失稳; 而下侧射弹先入水时, 其在入水后0.63 ms左右发生沾湿, 但此时射弹自身稳定性较强, 故未直接失稳, 而是在1.43 ms左右第二次沾湿后稳定性被破坏, 最终持续沾湿失稳。图 17b)中下侧射弹先入水时的2次峰值分别对应射弹的前后2次沾湿, 第二次峰值时横向阻力系数变化的绝对值更大的原因是射弹在第二次沾湿前的偏转过程中积攒了更多的惯性势能, 更大的惯性势能最终也导致射弹稳定性被破坏, 射弹持续沾湿失稳, 横向阻力系数持续增大。
|
| 图 16 上侧射弹先入水时前发射弹运动变化规律 |
|
| 图 17 下侧射弹先入水时前发射弹运动变化规律 |
由于不同入水顺序下前发射弹偏转角变化规律Δ=0.5与Δ=1.0, 2.0和3.0存在显著区别, 为深入分析两射弹前后间距对前发射弹运动稳定性的影响规律, 增加了Δ=0.6, 0.8, 1.5和2.5的工况, 并提取了各工况下前发射弹的偏转角, 如图 16c)和图 17c)所示。图 16c)和图 17c)为不同入水顺序下前发射弹偏转角随Δ的变化曲线, 包括5个不同的入水时刻(t1=0.5, 1.0, 1.5, 2.0, 2.5 ms)。由于射弹大幅度偏转后会发生失稳, 数值仿真停止, 因此图中部分入水时刻数据点缺失。观察16c)和图 17c)可发现, 随着Δ的增大, 不同入水顺序下前发射弹偏转角逐渐减小至曲线平稳。其中Δ=1.0为前发射弹偏转角变化曲线的转折点, 即当Δ≥1.0时, 前发射弹偏转角基本一致且较小, 同时随着入水时间的增加, 偏转角增长缓慢; 而当Δ < 1.0时, 前发射弹偏转角随着Δ的减小而逐渐增大, 且随着入水时间的增加, 偏转角增长速度较快。值得注意的是, 当下侧射弹先入水时, 前发射弹偏转角增长速度相较于上侧射弹先入水时较慢。
综合前后文分析, 可得出以下结论: 当Δ≥1.0时, 前发射弹偏转角均较小, 运动稳定性良好; Δ < 1.0时, 前发射弹偏转角均较大, 其中下侧射弹先入水时, 前发射弹发生失稳的时间较上侧射弹先入水时更晚, 其稳定性稍好。
图 18和图 19为不同入水顺序下后发射弹的无量纲偏转角Kdeg变化曲线、横向阻力系数Cy变化曲线和无量纲偏转角Kdeg随Δ的变化曲线, 定义后发射弹弹尖接触水面时t2=0, 下面对不同工况下后发射弹的偏转原因进行分析。
|
| 图 18 上侧射弹先入水时后发射弹运动变化规律 |
|
| 图 19 上侧射弹先入水时后发射弹运动变化规律 |
对比图 17和图 18中不同入水顺序下后发射弹的无量纲偏转角可以看出, 当Δ=1.0时, 不同入水顺序下后发射弹与单发射弹斜入水基本一致, 弹丸偏转角均较小且运动平稳; 当Δ=0.5时, 相较于下侧射弹先入水, 上侧射弹先入水时后发射弹偏转角变化幅度更大; 而当Δ=2.0和3.0时, 下侧射弹先入水时后发射弹偏转角变化幅度更大。综合前文图 13和图 15进行分析, 由前文可知上侧射弹先入水时, Δ=0.5工况下后发射弹所受Fp的方向与射弹所受流体阻力Y方向分量FY方向相同, Δ=2.0和3.0工况下Fp与FY方向相反; 而下侧射弹先入水时, 后发射弹在不同工况下所受Fp的方向与上侧射弹先入水时相反, 即Δ=0.5时后发射弹所受Fp与FY方向相反, Δ=2.0和3.0时Fp与FY方向相同。方向相反时力会相互抵消, 反之则不会。两力的相互抵消作用会导致此时后发射弹所受偏转力矩减小, 后发射弹偏转角减小, 进而导致Δ=0.5时上侧射弹先入水时的后发射弹偏转角更大, 而当Δ=2.0和3.0时下侧射弹先入水的后发射弹偏转角更大。
通过对比图 18b)和图 19b), 可观察到不同入水顺序下后发射弹横向阻力系数变化规律与上文分析基本一致, 即随着Δ从0.5增大到3.0, 上侧射弹先入水时, 后发射弹所受Fp与FY的方向由相同变为相反, 导致后发射弹横向阻力系数逐渐减小; 而下侧射弹先入水时, 后发射弹所受Fp与FY的方向由相反变为相同, 后发射弹横向阻力系数逐渐增大。
图 18c)和图 19c)给出了后发射弹在5个不同入水时刻(t2=0.5, 1.0, 1.5, 2.0, 2.5 ms)偏转角随Δ的变化曲线, 可以明显观察到随着Δ的增大, 2种入水顺序后发射弹偏转角均呈先减小后增大的趋势, 其中Δ=1.0为后发射弹偏转角的极小值点; 对于其他工况, 当Δ < 1.0时, 上侧射弹先入水时后发射弹偏转角随时间的增长速度较下侧射弹先入水时更快, 偏转角更大; 当Δ>1.0时, 下侧射弹先入水时后发射弹偏转角增长速度更快, 偏转角更大。
综上所述, 对于后发射弹来说, Δ=1.0为后发射弹偏转角的极小值点, 不同入水顺序下后发射弹均运动平稳, 其他工况下后发射弹均发生偏转; 当Δ < 1.0时, 下侧射弹先入水时后发射弹偏转角较上侧射弹先入水时更小, 运动稳定性较好; 而当Δ>1.0时, 上侧射弹先入水时后发射弹偏转角较小。
3 结论本文以超空泡射弹异步并联倾斜入水问题为探究对象, 运用数值模拟的方法分别对不同工况下前后两弹丸的空泡演化及运动特性进行了对比分析, 研究了异步并联射弹初始前后间距和入水顺序的不同对弹丸入水特性的影响规律, 主要得到了以下结论:
1) 异步并联倾斜入水超空泡射弹的初始前后间距对射弹空泡演化影响明显。超空泡射弹异步并联倾斜入水时, 后发射弹会对前发射弹内侧空泡形成挤压, 导致前发射弹空泡不对称, 同时后发射弹空泡在前发射弹空泡内部低气压的正激励作用下, 空泡向射弹内侧扩张。随着Δ从0.5增加至3.0, 前发射弹附近空泡受后发射弹空泡的影响逐渐减弱, 而后发射弹附近空泡受前发射弹空泡的影响持续存在。
2) 异步并联倾斜入水超空泡射弹的入水顺序对射弹空泡演化影响明显。Δ=0.5时, 由于射弹倾斜入水时存在水平方向速度分量, 射弹入水前期背水面附近流体动能和空泡直径相较于迎水面更大, 与自下侧入水相比, 自上侧入水的前发射弹空泡所受后发射弹的挤压更加剧烈; 随着Δ的增大, 入水顺序对射弹空泡演化的影响减弱。
3) 异步并联倾斜入水超空泡射弹的初始前后间距和入水顺序对射弹运动稳定性影响明显。对于前发射弹, 当Δ=0.5时不同入水顺序下射弹均较快发生失稳, 其中前发射弹自上侧入水时, 沾湿及失稳较下侧入水发生更快; 随着初始前后间距的增大, 2种入水顺序下前发射弹均逐渐运动平稳。对于后发射弹, 2种入水顺序下后发射弹运动规律基本一致, 在Δ=1.0时偏转角最小, 运动平稳; 其他工况偏转角均较大, 在Δ=0.5时向射弹内侧偏转失稳; 在Δ=2.0和3.0时向射弹外侧偏转失稳。
| [1] |
李佳川, 魏英杰, 王聪, 等. 不同扰动角速度高速射弹入水弹道特性[J]. 哈尔滨工业大学学报, 2017, 49(4): 131-136.
LI Jiachuan, WEI Yingjie, WANG Cong, et al. Water entry trajectory characteristics of high-speed projectiles with various turbulent angular velocity[J]. Journal of Harbin Institute of Technology, 2017, 49(4): 131-136. (in Chinese) |
| [2] |
陈晨, 魏英杰, 王聪, 等. 小型运动体高速倾斜入水空泡流动数值研究[J]. 兵工学报, 2019, 40(2): 334-340.
CHEN Chen, WEI Yingjie, WANG Cong, et al. Computational analysis of cavity flow induced by high-speed oblique water-entry of axisymmetric body[J]. Acta Armamentarii, 2019, 40(2): 334-340. (in Chinese) |
| [3] |
宋武超, 王聪, 魏英杰. 运动参数对回转体倾斜入水空泡演化影响[J]. 哈尔滨工业大学学报, 2018, 50(10): 130-136.
SONG Wuchao, WANG Cong, WEI Yingjie. Influence of motion parameters on the flow field characteristics of revolution body with oblique water-entry[J]. Journal of Harbin Institute of Technology, 2018, 50(10): 130-136. (in Chinese) |
| [4] |
路丽睿, 魏英杰, 王聪, 等. 不同头型射弹低速倾斜入水空泡及弹道特性试验研究[J]. 兵工学报, 2018, 39(7): 1364-1371.
LU Lirui, WEI Yingjie, WANG Cong, et al. Experimental investigation into the cavity and ballistic characteristics of low-speed oblique water entry of revolution body[J]. Acta Armamentarii, 2018, 39(7): 1364-1371. (in Chinese) |
| [5] |
秦杨, 易文俊, 管军. 超空泡射弹高速倾斜入水的空化流动数值模拟[J]. 兵器装备工程学报, 2019, 40(7): 99-104.
QIN Yang, YI Wenjun, GUAN Jun. Numerical simulation for cavitation flow of high-speed rotating body oblique entrying into water[J]. Journal of Ordnance Equipment Engineering, 2019, 40(7): 99-104. (in Chinese) |
| [6] | SHI Y, WANG G, PAN G. Experimental study on cavity dynamics of projectile water entry with different physical parameters[J]. Physics of Fluids, 2019, 31(6): 067103. DOI:10.1063/1.5096588 |
| [7] |
高英杰, 孙铁志, 张桂勇, 等. 回转体高速倾斜入水的空泡演化及结构响应[J]. 爆炸与冲击, 2020, 40(12): 123301.
GAO Yingjie, SUN Tiezhi, ZHANG Guiyong, et al. Flow characteristics and structure response of high-speed oblique water-entry for a revolution body[J]. Explosion and Shock Waves, 2020, 40(12): 123301. (in Chinese) |
| [8] | SONG Z J, DUAN W Y, XU G D, et al. Experimental and numerical study of the water entry of projectiles at high oblique entry speed[J]. Ocean Engineering, 2020, 211: 107574. DOI:10.1016/j.oceaneng.2020.107574 |
| [9] | LU L, YAN X, LI Q, et al. Numerical study on the water-entry of asynchronous parallel projectiles at a high vertical entry speed[J]. Ocean Engineering, 2022, 250: 111026. DOI:10.1016/j.oceaneng.2022.111026 |
| [10] |
王聪, 何超杰, 余德磊. 回转体并联入水运动状态预测[J]. 哈尔滨工业大学学报, 2021, 53(6): 21-26.
WANG Cong, HE Chaojie, YU Delei. Prediction of the state of revolution bodies entering water in parallel[J]. Journal of Harbin Institute of Technology, 2021, 53(6): 21-26. (in Chinese) |
| [11] |
路丽睿, 魏英杰, 王聪, 等. 双圆柱体低速并联入水过程空泡及运动特性试验研究[J]. 振动与冲击, 2019, 38(7): 42-49.
LU Lirui, WEI Yingjie, WANG Cong, et al. Tests for cavities and motion characteristics in process of two-cylinder in parallel water entry at low speed[J]. Journal of Vibration and Shock, 2019, 38(7): 42-49. (in Chinese) |
| [12] | 余德磊, 王聪, 何超杰. 回转体并联入水过程空泡及运动特性数值模拟[J]. 哈尔滨工业大学学报, 2021, 53(12): 23-32. |
2. School of Marine Science and Technology, Northwestern Polytechnical University, Xi'an 710072, China
