Analysis of characteristic control parameters of long-rod penetration
-
摘要: 针对理想长杆侵彻,通过对长杆侵彻Alekseevskii-Tate模型近似解进行分析,指出单一的无量纲速度衰减系数α(deceleration index)不足以完全表征长杆高速侵彻的准定常阶段。在此基础上,重新定义了2个无量纲特征参量:Johnson破坏数ΦJp和特征时间系数β,2个参量之间的关系为α=β/ΦJp。分析表明,ΦJp和β(或α和β)可实现对长杆高速侵彻准定常阶段的弹尾速度的完全表征;若再引入长杆弹相对临界速度vc*,则可完全表征长杆侵彻的准定常阶段。此外,还证明了α能够判定侵彻过程偏离定常状态的程度,并指出通过确定ΦJp和β(或α和β),可针对攻防需求对长杆弹侵彻设计进行指导。
-
关键词:
- 长杆侵彻 /
- Alekseevskii-Tate模型 /
- 特征参量 /
- Johnson破坏数 /
- 特征时间系数 /
- 无量纲速度衰减系数
Abstract: For ideal long-rod penetration, by analyzing the approximate solutions of the Alekseevskii-Tate model for long-rod penetration, it is pointed out that the single deceleration index α is not sufficient to fully describle the quasi-steady process of long-rod penetration. This paper redefines two dimensionless parameters, namely Johnson demage parameter ΦJp and characteristic time parameter β, and α=β/ΦJp. The analysis shows that two characteristic parameters ΦJp and β (or α and β) can completely characterize the impact velocity of the projectile tail in the quasi-steady process of long-rod penetration. If the dimensionless critical impact velocity vc* is introduced, the quasi-steady process of long-rod penetration can be fully characterized. In addition, this paper strictly proves that the degree of deviation from the steady state in the penetration process can be determined by α, and confirms that by determining ΦJp and β (or α and β), the design of long-rod penetration can be guided for offensive and defensive needs. -
高压气体在运输、储存、加工和使用过程中常发生爆炸事故, 造成巨大的人员伤亡和财产损失, 同时, 由于爆炸过程还蕴含了许多诸如界面不稳定性以及激波与界面相互作用等复杂流体物理现象, 因此对气云物理爆炸特性的研究不仅具有重要的实际意义, 还具有重要的理论价值。
自从G.I.Taylor[1]提出爆炸波冲击理论以来, 不少学者做过类似研究[2-3]。事实上, 气相爆炸中包含了极其丰富的物理现象, 特别是高压、高密度的气体云爆炸, 其包含了在激波作用下自由剪切层的发展、Richtmyer-Meshkov (RM)不稳定性以及湍流的转捩、气动噪声等基本的流体物理现象。其中RM不稳定性是指不同密度的流体分界面在运动激波的作用下, 在界面附近形成不断扩展的湍流混合层, 这种复杂的非线性现象在超燃冲压发动机的混合燃烧[4]、爆轰[5-6]、惯性约束聚变[7]以及天体物理的超新星爆发等问题中起重要作用, 因而RM失稳现象得到了广泛的关注, 研究人员做了许多激波与密度界面相互作用的研究[8-9]。C.A.Zoldi等[10]在激波管上开展了激波加载SF6气柱RM界面不稳定性实验, 并采用RAGE程序模拟了SF6气柱在空气冲击波作用下界面的演化、发展过程以及后期空气和SF6气体的混合。G.Layes等[11]利用高速摄影方法研究了激波对不同密度气泡作用后的变形过程。邹立勇等[12]利用高速摄影测试技术, 对弱激波冲击下空气中的SF6气柱和气帘界面的演变过程进行了实验研究, 对演变过程中涡的产生和发展进行了初步解释。而事实上, 对于气相爆炸, 人们关注更多的是由球形激波引起的RM失稳。J.G.Zheng等[13]采用数值模拟的方法研究了由内聚球形激波引起的RM失稳现象, 他们的研究显示了平面激波和球形激波作用的不同特征。
柱形或球形气云爆炸所产生的RM不稳定性由于气云内部反射激波的作用变得更复杂, 而且其不稳定性后期非线性演化可加速可压气体的湍流转捩。本文中, 采用大涡模拟(large eddy simulation, LES)方法以及高精度混合格式对高压、高密度的SF6气体云冷爆炸进行数值模拟, 以研究其爆炸过程中RM失稳以及二次激波与其作用后期非线性演化的湍流过程。
1. 控制方程和数值方法
三维可压Navier-Stokes(N-S)方程经过Favre滤波后, 其标准的可压三维LES方程作为控制方程。连续性方程、动量方程和能量方程分别为:
∂−ρ∂t+∂∂xi(ˉρ˜ui)=0 (1) ∂ˉρ˜ui∂t+∂∂xj(ˉρ˜ui˜uj)=−∂ˉpδij∂xi+∂σij∂xj−∂τsgsij∂xj (2) ∂−ρE∂t+∂∂xj(ˉρ˜E+ˉp)˜uj=∂∂xj(ˉk∂˜T∂xj−qT−ggsj+˜uiσij) (3) 式中:上标“-”表示空间平均, “~”表示Favre滤波; ρ为密度, p为压力, T为温度, ui为i方向速度, δij为Kronecker函数; 滤波黏性应力张量
; 滤波压力
为理想气体常数,
为摩尔质量
为亚格子应力张量, 总能
为热通量。
在大涡模拟中, 亚格子湍流应力项无法直接从可解尺度中计算求得, 因此亚格子应力需要引入模型进行封闭。这里采用拉伸涡亚格子模型[14]对小尺度输运项建模。拉伸涡亚格子模型由D.I.Pullin[14]提出, 他假设亚网格内的流动是由亚网格内的涡导致的, 并通过对称涡结构来构建小尺度涡量。拉伸涡亚格子模型已被成功用于亚格子能量谱的估算以及湍流混合模拟。
以上LES方程可用有限体积法进行数值求解。由于爆炸过程包括激波与气体相互作用加速诱导湍流转捩, 因此对方程中的对流项离散要求很高, 需数值格式同时具备对激波与湍流的高分辨率。一般采取高阶的中心差分与迎风型格式的混合形式, 本文中采用TCD(tuned centered-difference)/WENO(weighted essentially nonoscillatory)混合数值方法[15]。时间推进采用三阶精度的Runge-Kutta法。
2. 计算模型
在爆炸初期, 流场结构呈现对称情形, 而由于RM失稳, 后期的流场渐渐转为非对称结构。为描述RM失稳以及后期二次激波的作用过程, 选择初始状态为球形气云的SF6气体作为研究对象。SF6气体和空气的密度分别为5.97和1.18 kg/m3, 两者的比热比和气体常量分别为γ(SF6)=1.09, R(SF6)=56.92; γ(air)=1.40, R(air)=287。选择标准状态的空气参量, p0=101.3 kPa, T0= 287 K。SF6气体的初始压力为50p0, SF6气云初始半径为R0=0.1 m, 爆炸空间径向大小为15R0。计算采用笛卡尔网格模拟气云球体表面, 以该网格扰动作为初始随机扰动, 球体表面网格总量约为800万, 整个计算域网格总量约为4 500万。
3. 结果与讨论
当激波从重质气体传向轻质气体时, 会产生向内传播的稀疏波和向外传播的透射激波, 从而引起失稳, 加速界面混合。激波作用不同密度气体界面后, 界面演化一般经历3个阶段[16-17]:(1)振幅线性增长阶段; (2)非线性增长阶段, 轻质流体发展成“气泡”, 重质流体发展成“尖钉”; (3)湍流混合阶段, “尖钉”开始破碎, 不同尺度的涡相互吞并, 不同流体达到湍流混合状态。从图 1中可以清晰地看到气体界面的3种状态:在膨胀初期, 随机扰动在气体界面发生, 扰动振幅在界面线性增长, 见图 1(a); 随着膨胀的加剧, 振幅经历非线性增长, “气泡”与“尖钉”结构出现, 见图 1(b); 随着膨胀的结束以及球心处低压区的形成, 圆球界面开始向球心回流, “尖钉”结构破碎, 不同流体开始形成湍流混合状态, 见图 1(c)。
为了能更清晰地展现爆炸过程中流场内部激波与球体界面相互作用的过程, 图 2给出了高压球体爆炸过程中对称面不同时刻的密度纹影。由图 2可看出, 开始时透射激波向外传播、稀疏波向内传播以及界面失稳过程中的详细变化。开始时, 内部的SF6气体由于高压往外膨胀, 其形成的入射激波和气体分界面作用分成透射激波和反射的稀疏波, 透射激波一直往外传播, 而稀疏波则向内传播, 见图 2(a)~(b)。此时为振幅线性增长阶段, 气体分界面在透射激波的作用下加速往外膨胀, 且气体分界面上产生的随机增长呈线性变化。随后为非线性阶段, 重质的SF6气体进入空气中形成尖钉结构, 而轻质的空气演化为气泡结构, 见图 2(c), 且发展过程中界面RM不稳定性呈现非线性增长, 振幅越来越大。随着稀疏波的向内传播以及气体分界面向外的膨胀, 气体分界面内部的气体压力迅速降低, 导致气体分界面开始向内急剧收缩, 见图 2(d)~(e), 从而加速界面部分流场的层流向湍流转捩。至t=10 ms时, 反射的稀疏波在圆心处汇聚形成反射的二次激波, 见图 2(f)。
图 3给出了不同时刻对称面上的径向压力p曲线, 其中p1为透射激波前缘压力, p2为反射稀疏波压力。初始时刻, SF6气体在内部高压下急剧膨胀, 经过气体界面产生透射激波和反射稀疏波, 透射激波在p1的作用下向外传播。从t=1.4 ms到t=3.0 ms, 随着透射激波向外传播, p1衰减较快, 但是透射激波波速很高。而在t=3.0 ms之后, 图中已无p1轨迹, 由此说明透射激波此时已传出计算域。反观稀疏波, 其压力远低于透射激波压力:在形成初期, 稀疏波随气体界面向外传播, p2逐渐减小(t=0~3.0 ms); 而随着SF6气体急剧膨胀, 在球心形成极低压区; 随后, 稀疏波反向向内回传, 最终在t=10 ms时在原球心位置会聚形成二次激波, 此时p2达至其极大值。
图 4显示了气体分界面、透射激波, 稀疏波前缘和后缘以及后期形成的二次激波的变化轨迹。从图中可看出, 透射激波以约600 m/s的速度在空气中往外传播, 在该激波的作用下气体分界面加速膨胀, 在t=40 ms之前, 气体界面一直处于线性膨胀阶段, 之后, 由于球心处的低压导致其速度逐渐降低, 最终在t=65 ms时开始向内收缩。同时从图中可以看到稀疏波的传播:稀疏波前缘在初始时向内传播, 在球心汇聚后加速向外传播, 而其后缘则一直低速向外传播; 稀疏波的前缘约在t=2.8 ms时追上其后缘, 从而汇合成同一道波, 与气体界面类似, 两者汇聚后一直向外传播, 直到t=5 ms时开始向内传播, 最终在t=10 ms时在球心再次会聚碰撞, 形成二次激波; 该激波线性向外膨胀, 直到t=124 ms时与向内传播的气体界面相互作用, 导致后期流场湍流发展。
图 5则刻画了在圆心处会聚反射的二次激波往外传播及其与流场的过程。此阶段反射激波的强度很高, 见图 3(f), 且气体界面已不甚明显, 轻质“气泡”和重质“尖钉”处的流场基本形成湍流。此时反射激波往外传播, 而湍流区界面则向内传播, 因此两者相互作用强烈, 见图 5(c), 最终导致整个流场呈现完全湍流状态, 见图 5(d)。
4. 结论
采用大涡模拟方法, 结合高阶混合格式, 对SF6重质气云在空气中的爆炸过程进行了数值模拟。模拟结果清晰地再现了界面演化的整个过程。爆炸产生的激波经过气体分界面时分为透射激波和反射稀疏波。透射激波向外传播, 导致气体分界面处RM失稳增强, 加速了2种气体的混合, 而反射稀疏波后缘低速向外传播, 其前缘在初始时向内传播, 在球心汇聚后加速向外传播, 直至与稀疏波前缘汇合, 两者随着气体界面一直向外膨胀。之后, 由于球心处的低压导致回流, 最终在球心处形成二次激波。在后期的发展中, 该强激波与气体分界面相互作用, 使整个流场区域呈现完全湍流状态。
-
表 1 长杆侵彻设计工况中的相关参数
Table 1. Related parameters of long-rod penetration design
v0/(km·s−1) L/mm ρp/(g·cm−3) ρt/(g·cm−3) Yp/GPa Rt/GPa 1.5 81.7 17.4 7.8 2.0 4.94 表 2 设计工况中相关参数
Table 2. Summary of parameters in the designed cases
工况 v0/
(km·s−1)L/
mmρp/
(g·cm−3)ρt/
(g·cm−3)Yp/GPa α/% ΦJp β 1 1.5 100 19.00 9.00 2.0 11.48 21.375 2.453 2 1.5 50 19.00 9.00 2.0 11.48 21.375 2.453 3 3.0 100 19.00 9.00 2.0 2.87 85.500 2.453 4 3.0 100 13.00 9.00 2.0 3.76 58.500 2.202 5 3.0 100 13.00 5.00 2.0 4.47 58.500 2.612 6 1.5 100 19.00 9.00 1.0 5.74 42.750 2.453 7 1.5 100 23.75 5.56 2.0 11.48 26.719 3.067 8 1.5 100 19.00 12.56 2.2 11.48 19.432 2.230 9 1.5 100 19.00 9.00 0 0 ∞ 2.453 -
[1] 焦文俊, 陈小伟. 长杆高速侵彻问题研究进展 [J]. 力学进展, 2019, 49(1): 312–391. DOI: 10.6052/1000-0992-17-021.JIAO W J, CHEN X W. Review on long-rod penetration at hypervelocity [J]. Advances in Mechanics, 2019, 49(1): 312–391. DOI: 10.6052/1000-0992-17-021. [2] ALLEN W A, ROGERS J W. Penetration of a rod into a semi-infinite target [J]. Journal of the Franklin Institute, 1961, 272(4): 275–284. DOI: 10.1016/0016-0032(61)90559-2. [3] ALEKSEEVSKII V P. Penetration of a rod into a target at high velocity [J]. Combustion, Explosion and Shock Waves, 1966, 2(2): 63–66. DOI: 10.1007/BF00749237. [4] TATE A. A theory for the deceleration of long rods after impact [J]. Journal of the Mechanics and Physics of Solids, 1967, 15(6): 387–399. DOI: 10.1016/0022-5096(67)90010-5. [5] TATE A. Further results in the theory of long rod penetration [J]. Journal of the Mechanics and Physics of Solids, 1969, 17(3): 141–150. DOI: 10.1016/0022-5096(69)90028-3. [6] 焦文俊. 长杆高速侵彻的1D和2D理论模型 [D]. 合肥: 中国科学技术大学, 2019. [7] JIAO W J, CHEN X W. Analysis of the velocity relationship and deceleration of long-rod penetration [J]. Acta Mechanica Sinica, 2019, 35(4): 852–865. DOI: 10.1007/s10409-019-00862-1. [8] JIAO W J, CHEN X W. Approximate solutions of the Alekseevskii-Tate model of long-rod penetration [J]. Acta Mechanica Sinica, 2018, 34(2): 334–348. DOI: 10.1007/s10409-017-0672-9. [9] WALTERS W, WILLIAMS C, NORMANDIA M. An explicit solution of the Alekseevski-Tate penetration equations [J]. International Journal of Impact Engineering, 2006, 33(1−12): 837–846. DOI: 10.1016/j.ijimpeng.2006.09.057. [10] ANDERSON Jr C E, WALKER J D. An examination of long-rod penetration [J]. International Journal of Impact Engineering, 1991, 11(4): 481–501. DOI: 10.1016/0734-743X(91)90015-8. -