The electromagnetic radiation produced by hypervelocity impact
-
摘要: 超高速碰撞产生的电磁辐射是固体物质在强冲击作用下的重要物理响应，在深空探测、航天器对空间碎片的防护设计、武器毁伤评估应用广泛。本文中概述了超高速碰撞产生的电磁辐射现象，总结了不同碰撞条件下，超高速碰撞产生微波和闪光的时频特性；从超高速碰撞产生材料破碎和产生等离子体两个方面，分析了超高速碰撞产生微波的辐射模型；归纳了超高速碰撞下的发光机理，并阐述了超高速碰撞产生连续光谱和线谱的辐射模型，指出了超高速碰撞产生电磁辐射研究存在的不足与发展趋势。Abstract: The electromagnetic radiation from hypervelocity impact is an important physical phenomenon of solid matter under strong impact loadings, and the research results have important application value in fields of deep space exploration, protection design of spacecraft against space debris, assessment of weapon damage, etc. In this paper, the electromagnetic radiation caused by hypervelocity impact was briefly summarized. The time-frequency characteristics of microwave and flash by hypervelocity impact under various collision conditions were respectively provided. The radiation models of microwave generated by hypervelocity impact were analyzed from two aspects of material fragmentation and plasma phase transition. The luminescence mechanisms during hypervelocity impact were described as a whole, and the radiation models of continuous spectrum and line spectrum were performed. The shortcomings and development trend of electromagnetic radiation by hypervelocity impact were pointed out. The results show that the intensity and frequency of microwave and flash are closely dependent on the target thickness, impact material, environmental pressure as well as collisional velocity and angle. The microwave is in the form of pulses lasting from a few to hundreds of microseconds. The intensity of the flash, however, accumulates rapidly and then decays slowly. When the plasma is not accomplished during the hypervelocity impact, the microwave radiation is mainly formed because of the movement of ionized debris. Once the plasma is formed, the effects of the collision radiation in the plasma, namely the bremsstrahlung and recombination radiation, and the expansion of the plasma should also be taken into account, and the specific spectrum lies on the characteristic parameters of plasma. In vacuum, only the gasification and plasma due to the impact are necessary to be considered in the flash spectrum, while the gasification and even plasma phase transition resulted from the ablation between the gas and debris should be involved at high environmental pressure.
-
Key words:
- hypervelocity impact /
- electromagnetic radiation /
- plasma /
- flash /
- microwave
-
激波与两种不同密度气体界面相互作用时，会在界面附近发生反射、折射及绕射等现象，同时在界面处还会产生Richtmyer-Meshkov(RM)和Kelvin-Helmholtz(KH)不稳定性，导致界面扭曲变形。该现象广泛存在于超燃冲压发动机中的燃料混合、惯性约束核聚变、水中炸药爆炸及天体物理中的超新星爆炸等领域，因而具有重要的研究价值。
目前，对激波与气体界面相互作用的研究主要集中以下两方面：一是气体界面在RM不稳定性与KH不稳定性作用下的演变过程及其相关复杂流动机理；二是界面演化后期的失稳、破碎及湍流混合问题。K.A.Meyer等[1]对激波与气体界面作用时的RM不稳定性进行了数值模拟，得到了与Meshkov实验定性一致的结果；Q.Zhang等[2]研究了二维和三维初始扰动的发展过程，得到了气泡、尖钉速度及界面发展的总体增长率的计算公式；N.N.Anuchina等[3]对初始扰动幅度下相同入射强度激波与正弦扰动气体界面的作用过程进行了数值模拟，结果表明三维数值计算得到的扰动增长率高于二维结果；G.A.Ruev等[4]对马赫数为3.5的平面入射激波与正弦扰动气体界面的作用过程进行了数值模拟，给出了气体界面RM不稳定性的发展过程及两种气体混合区的宽度变化；J.H.J.Niederhaus等[5]对不同强度平面入射激波与多种不同介质的气柱、气泡相互作用过程进行了大量的数值研究，重点分析了作用过程中流场环量的变化和两种气体的混合机理；同时，B.Thornber等[6]对低马赫数的平面入射激波与R22方柱作用过程的数值模拟，结果与Holder等的实验结果完全相吻；B.Hejazialhosseini等[7]对不同马赫数下平面入射激波与氦气气泡作用过程进行了数值研究，得到了气泡的变形过程及气泡表面在KH不稳定作用下次级涡的产生过程；O.Schilling等[8]采用数值方法研究了平面入射激波与Air/R22气体界面的相互作用过程，并分析了反射激波与分界面作用前后流场结构和混合层演变；B.Tian等[9]利用ALE(arbitrary Lagrangian-Eulerian)和GALE(global arbitrary Lagrangian-Eulerian)方法对多种激波与气体界面作用问题进行了数值模拟；S.K.Shankar等[10]采用数值方法模拟了平面入射激波与重气柱作用过程，结果表明流场初始压力和密度梯度是影响两个主涡产生的关键因素；C.Bailie等[11]对不同强度平面入射激波与Air/R22正弦小扰动界面的作用过程进行了数值模拟，并分析入射激波强度对作用过程的影响。
目前，对界面RM不稳定性控制方法的研究表明，合适的磁场对气体界面RM不稳定性的产生与发展具有一定的抑制作用，但仍处在理论研究阶段。S.Chandrasekhar等[12]和V.Wheatley等[13-14]采用线性理论方法推导了垂直或平行于界面的磁场对不可压缩流体界面不稳定性的影响，结果表明：平行于界面的磁场能够抑制RT和RM不稳定性；但是在线性理论下，平面流场中的RM不稳定性与垂直于平面的磁场无关，与实际情况不符。为进一步研究磁场效应，M.Khan等[15]在势流理论的基础上建立了新模型，发现考虑非线性效应后，磁场能够影响界面不稳定性的发展。J.Cao等[16]同时考虑了初始平行于界面的磁场和剪切流对界面不稳定性的影响，发现磁场能够抑制界面不稳定性，而剪切流将会促进界面不稳定性的发展，两者之间存在竞争机制。M.S.Shin等[17]对不同方向磁场作用下，平面激波与球形气泡相互作用过程进行了数值模拟，发现在发展初期，磁场效应很小，而在发展后期，即使是弱磁场，其对界面形状的影响也很明显。此外，磁场越强，界面内外物质混合率越低，并越能抑制后期湍流的发展。李源等[18]基于理论分析的方法，综合考虑黏性和表面张力的影响，考察了磁场中非理想流体中RT不稳定性气泡的演化过程，并推导出描述二维非理想磁流体RT不稳定性气泡运动的控制方程组，得到了流体黏性、表面张力和磁场对气泡界面不稳定性发展的影响。
综上所述，目前对磁场控制不稳定性的研究主要集中在基于不可压理想磁流场的解析推导及磁场影响下所产生现象的数值模拟两个方面。本文中，基于磁流体动力学(magneto-hydro-dynamics, MHD)方程，采用CTU+CT算法对激波冲击磁化R22气柱进行数值模拟，并以磁压力为切入点研究磁场对不稳定性及射流的影响。
1. 数值方法及计算模型
1.1 数值方法
在满足流体力学方程的条件下，带电荷的流体团更有利于形成团块，主要原因如下：(1)带电粒子间的碰撞频率远大于中性粒子之间的碰撞频率，因此其碰撞平均自由程较短；(2)有磁场存在时，磁场对带电粒子有侧向约束。洛伦兹力使带电粒子在垂直于磁场方向只能保持在回转半径范围内运动。等离子体通常比中性粒子有更低的密度和更高的温度，由于整个过程在短时间内完成，本文中考虑高磁雷诺数的情况，使用MHD方程来研究激波冲击R22气柱中出现的不稳定性：
∂ρ∂t+∇⋅[ρv]=0 (1) ∂(ρv)∂t+∇⋅[ρvv−BBμ+p∗]=−∇⋅Π (2) ∂E∂t+∇⋅[(E+p∗)v−B(B⋅v)μ]=−∇⋅Q−∇⋅(Π⋅v) (3) ∂B∂t−∇×(v×B−ηJ−ηH(J×B)−ηAD(J×B)×B)=0 (4) p∗=p+B⋅B2μE=p/(γ−1)+ρ(v⋅v)2+B⋅B2μ (5) J=1μ∇×BQ=−κ0∇T−κ∥ˆbˆb⋅∇T (6) Π=−ν0∇v−3ν∥(ˆbˆb−13I)(ˆbˆb−13I):∇v (7) 式中：ρ为质量密度，v为速度矢量，p为压力，E为能量， μ为磁导率，B为磁感应强度，ˆb为B的单位矢量，Π为包含各向同性与各向异性部分的黏性应力张量，分别由系数ν0和ν‖控制，γ为比热比，T为温度，I为单位张量，符号“:”表示双点积运算。在各向异性的情况下，黏性通量限制在平行于磁场线方向。Q为热通量，同样包含各向同性与各向异性的情况，分别由系数κ0和κ‖控制。在各向异性的情况下，热通量同样限制在平行于磁场线方向。同时，感应方程包含多种非理想MHD效应，包括欧姆损耗(由电阻率η控制)，霍尔效应(由系数ηH控制)和双极性扩散(由系数ηAD控制)。
本文中采用CTU+CT算法对MHD方程组进行求解，其中CTU算法是一种用于求解守恒双曲系统的不分离的二维有限体积算法，其三维变形则由J.Saltzman[19]提出。另外，由于传统的12-solve CTU算法忽略了MHD方程中始终存在的约束▽·B=0，不能作为求解MHD方程的有效算法，因此需要构建一个Godunov通量中网格单元边缘的平均电场算法，通常称为CT算法，该算法可以描述为以Godunov电场为预测值，以CT电场为校正值的预测过程。鉴于在文献[20]中ε c CT已显示的最佳性能，本文中采用ε c CT算法。该算法包含逆风偏差(根据接触模式)，可以将网格对齐的平行平面流动降低到正确的Godunov电动势。
1.2 计算模型
图 1为激波与R22气柱相互作用的计算模型，其中R22气泡直径为0.05m，气柱中心与左边界的距离为0.05m，与右边界的距离为0.25m，计算域高度为0.089m。初始时，气柱内采用R22作为重质气体，外部介质为空气，气体参数如表 1所示。气柱内外压力均为101325Pa，入射激波从左向右传播，来流马赫数Ma=1.22，上、下边界设为固体反射边界，而右边界设为出口条件。计算域网格采用均匀分布的笛卡尔网格，经网格收敛性测试后，区域网格总数设置为3000×890。本文中，设初始磁场强度为0.01T(β≈2500)，令平行于来流的方向为x方向，垂直于来流的方向为y方向。
表 1 气体参数Table 1. Gas parameters used in this paper气体 比热比 密度/(kg·m-3) 摩尔质量/(g·mol-1) 当地声速/(m·s-1) R22 1.185 3.69 86.468 183 Air 1.400 1.18 29.000 345 为了使气体受磁场影响，首先通过高温气体电离，温度为12000K。由于含等离子体，电导率、热导率、黏性系数、霍尔系数以及双极扩散系数由以下方法确定。电导率采用基于玻尔兹曼方程得到的Spitzer公式进行计算，即由玻尔兹曼方程的碰撞积分出发，忽略短程碰撞的影响，考虑剩余的碰撞项，得到弱电场情况下由离子与电子所组成的等离子体电导率(σ)为：
σ=2mγ03π3/2e2j3lnΛ (8) 式中：j= m/(2kB T)，其中m为质量，kB为Boltzmann常数；γ0=0.1 Z+0.49；e为电子电荷量；Z为离子电荷数。
如果碰撞积分到德拜长度为止，则库仑对数项为：
Λ=3/π1/2n1/2ee3(kBT/2)3/2 (9) 计算可得，初始电导率约为107S/m。
当气体中相当部分的原子被电离后，静电力将对粒子间碰撞起主要作用，若考虑远程作用，则离子碰撞将不能视为两体方式。通过求解玻尔兹曼方程得到一阶近似电离气体热导率为：
{kc=5νcV/2ν=1.75×10−20ˉm12[kB(T−273.15)]52z2e4Λ (10) 式中：ν为一阶近似黏度，m为等离子体平均质量。采用动力论有关比热关系式，cV=N0kB/(2m)，其中N0为常数。计算可得热导率的表达式为：
kc=2.97×10−20(2π)12kB[kB(T−273.15)]ˉm12e4z2Λ52 (11) 根据S.Chapman等[22]的理论可以推导出kc、ν的二阶近似值。尽管二阶近似值与一阶近似值有显著的差异，但其后各阶近似所增加的修正量都是很微小的。假定黏度均由重离子引起，电子黏度很小，由于电子的热运动速度远大于离子的热运动速度，因此热扩散效应主要是由电子引起。经过估算：本文中黏性系数约为3.72×10-5Pa·s，热导率约为1.4W/(m·K)。
由于准电中性条件，电子和离子在等离子体内不可能独立运动，当电子快速离开某等离子体单元时，单元内将形成由正电荷产生的电场，该电场阻碍电子继续离开所在体积单元，并促使离子更快地离去，这种扩散状态称为双极型状态。由于电子的扩散系数远大于离子，因此：
ηAD≈ηi(1+TeTi) (12) 式中：ηi为正离子的热扩散系数，Te、Ti分别为电子、离子温度。由此可见，双极扩散系数恒小于电子自由扩散系数，而大于离子扩散系数，因此双极性电场极大影响了电子的定向速度。本文中取双极扩散系数为(10-5/ ρ)m2·Pa·s-1。
Hall效应可在磁流体通道中产生横向洛伦兹力，而且体力的不均匀分布可导致二次流的产生，对流体的运动产生重要影响，因此需要考虑霍尔效应。此处，霍尔系数ηH取为10-7m3/C。
2. 结果与讨论
为了验证数值方法与计算模型的可靠性，图 2给出了不加磁场时入射激波与R22气柱相互作用过程计算结果与文献[21]中实验结果的对比，两者具有相同的初始条件。由图 2可知：计算结果中激波的反射、绕射及发展过程均与实验结果完全相符，两者气柱的变形过程也相似；但是，由于初始时两种气体接触界面采用间断面形式，因此计算结果中接触界面上的KH不稳定性更明显，从而导致界面上、下两侧的湍流转捩过程(见图 2(e))与实验结果稍有不同。
图 3为后期激波与R22气柱作用过程的计算纹影。从图 3可以看出：由于激波对不同密度气体界面的作用，界面产生了RM不稳定性；在作用初期，随机扰动在气体界面发生，扰动振幅在界面线性增长；随着变形的加剧，振幅经历非线性增长，“气泡”与“尖钉”结构出现；随着运动的继续进行，在KH不稳定性的影响下，“尖钉”破碎，“气泡”被拉长，最终形成湍流混合状态。
图 3 后期激波与R22气柱作用过程计算纹影Figure 3. Numerical results of the interaction of shock wave with the R22 air column in later period.图 4为加入均匀磁场后，激波与R22气柱相互作用过程的密度纹影，其中上方子图的磁场方向垂直于来流方向(即磁场方向沿y轴正向)，而下方子图中磁场方向平行于来流方向。由图 4可知，与无磁场情况(见图 2)相比，磁作用力能明显抑制界面处的KH不稳定性，且垂直于来流方向的磁作用力对不稳定性的控制较平行于来流方向的磁作用力更明显。为了更清晰地描述此现象，将动量方程中与磁场相关的项拆分为磁压力项▽(B·B)/(2 μ)与磁张力项▽•(BB)/ μ两部分，这两项均为二阶张量。然而，由于计算结果中磁压力值远大于磁张力值，下文只对磁压力的变化情况进行讨论。
图 4 B=0.01T时不同磁场方向对激波与R22气柱作用过程的影响Figure 4. Numerical results of the interaction process of shock wave and R22 air column when B=0.01T由式(4)可知，垂直于速度方向的磁场分量梯度较大，磁场变化较快，该效应会引起磁场在某些位置(如垂直于磁场方向的界面处)的集中。因此，随着时间的推移，气泡界面处的磁场强度不断增大，妨碍轻质气体和重质气体的相互混合。由图 4可知：当初始磁场垂直于来流方向时，磁压力作用在x轴方向，可以有效阻止气柱两侧的卷起；在t=450 μs之后，由于斜压作用加剧，气柱尾部界面开处始形成涡结构；然而，沿x轴方向的磁压力抑制了不稳定性的进一步发展，最终在气柱尾部界面附近形成了上、下两个突起，如图 4(f)所示。当磁场平行于来流方向时，磁压力作用在y轴方向，因上侧界面处x轴向速度的差异，在界面处卷起涡串，如图 4(c)和图 4(d)所示；同时，由于斜压效应，气柱上、下两侧卷起两个涡环，如图 4(e)所示；随后，沿y轴方向的磁压力阻止了主涡的进一步发展，最终在主涡位置形成“尖钉”状结构，见图 4(f)。由于磁化气柱的热导率远大于空气的，界面处气体混合速度较无磁场时更快。
图 5和图 6分别给出了磁场方向垂直或平行于来流方向时气柱对称轴上压力(p)与x方向速度(u)的分布。由图 5和图 6可知，与无磁场的情况[23]相比，初始压力与速度变化不大。当t=200 μs时，透射激波在x轴上发生聚焦，并且在汇聚处产生局部高压区：初始磁场沿y轴正方向时，压力曲线峰值pm=700kPa，初始磁场沿x轴正方向时，pm=760kPa，均小于无磁场情况。此时，x方向速度同样处于峰值，但施加沿x轴方向磁场时的速度(u=153m/s)远大于无磁场(u=100m/s)及施加沿y轴方向磁场时的速度(u=125m/s)。随着相交点向下游移动，聚焦点的压力峰值不断下降，速度逐渐增大，并在t=250 μs时达到最大值(u≈170m/s)。随后，在流体阻力作用下，速度不断衰减。综上所述，磁场的加入对压力影响不大，但却能提升透射激波聚焦时速度峰值，而磁场沿x轴方向时影响更剧烈。此外，在激波衰减过程中，沿y方向磁场较沿x方向磁场更能减缓射流的衰减。
图 5 初始磁场沿y轴正方向时入射激波诱导射流产生过程Figure 5. Generation process of jet induced by incident shock wave when initial magnetic field is along y axis
图 6 初始磁场沿x轴正方向时入射激波诱导射流产生过程Figure 6. Generation process of jet induced by incident shock wave when initial magnetic field is along x axis图 7给出了磁场沿y轴正方向或x轴正方向时激波与R22气柱作用过程的计算纹影。由图 7可知，磁场能够较好抑制RM不稳定性。当初始磁场垂直于来流方向时，界面处的横向磁压力对气体混合的抑制效果较好，同时气柱截面被压扁拉长，两个主涡在磁场作用下受到抑制，最终被拉成长条状；当初始磁场平行于来流方向时，后期气柱宽度小于磁场垂直于来流方向的情形，虽然形成两个主涡，但磁场抑制了气柱的破碎及次级涡的生成。由此可见，加入磁场后，气柱界面处没有在KH不稳定性的作用下持续失稳并卷起“珠状”次级涡，气柱顶端也未出现RM不稳定性，因此，两种方向的初始磁场对不稳定性都有较好的抑制作用。
图 7 后期激波与R22气柱作用过程密度纹影Figure 7. Numerical results of interaction process of shock wave with R22 air column in later period为更清楚地看到磁场对气柱界面的影响，图 8给出了初始磁场沿y轴正方向时不同时刻磁压力和磁能量的分布。由图 8可知，初始时刻，透射激波在气柱内部向下游传播，圆弧反射激波向上游传播，磁压力较大处集中在气柱界面处，磁能量较大处集中在界面处及透射激波尾端。当t=165 μs时，x方向磁压力沿x轴对称分布，而y方向磁压力沿x轴反对称分布；同时x方向磁压力对气体界面内侧产生向右拉力，对气体界面外侧则产生向左拉力，而y方向磁压力对气体界面内侧产生向气柱中心的拉力，而对气体界面外侧产生向外的拉力。另外，当入射激波与气柱作用时，因气柱内外气体声阻抗的不同，气柱内外会产生速度差。气柱外部界面处，空气x方向速度大于界面速度，而气柱内部界面处，空气x方向速度小于当前界面速度。由于磁场变化率与速度差密切相关，因此这两处的速度梯度会带来磁场的变化，由此产生的磁压力进一步抑制气柱界面受内外剪切力的影响(同理，对空气y方向速度，也有类似结论)。由图 7(c)和图 7(d)可知，发展后期，磁能量与磁压力分布类似，包裹在界面处。综上所述，当初始磁场垂直于来流方向时，磁压力始终包裹着界面，阻止界面受周围环境的进一步影响，从而防止界面上涡量的生成和扩散，抑制界面的不稳定性。
图 8 初始磁场沿y轴正方向时磁压力及磁能量的分布(每幅子图中，左上为磁压力，左下为磁能量，中间和右侧分别为磁压力在x、y方向的分量)Figure 8. Magnetic pressure and magnetic energy when initial magnetic field is along y axis (upper left: magnetic pressure; lower left: magnetic energy; middle: magnetic pressure along x axis; right: magnetic pressure along y axis)图 9为初始磁场沿y轴正方向时磁压力(pm)和磁能量(Em)沿横向(x方向)及纵向(y方向)的分布情况。由图 9可知：当透射激波在气柱内传播时，磁能量在两个位置取值较大，一处是气柱左侧界面处，另一处位于透射激波尾部，与之前分析相符；界面处的磁压力远大于其他位置。由图 9可知，气柱上、下界面处的磁压力及磁能量远大于左、右界面处。图 9(c)显示：当透射及入射激波先后聚焦后，会在气柱尾端射流位置形成局部高磁压力区域，同时，该区域磁能量很低，接近于零。由于图 9(c)和图 9(d)中射流位置的磁压力及磁能量很低，因此磁场对射流并无明显影响；但是射流前、后的磁能量可以一定程度上抑制射流的衰减。此外，图 9中不同时刻磁压力和磁能量分布曲线的形状无明显变化，但上、下界面处磁压力和磁能量的峰值随时间不断增加，这意味着界面处磁场强度随时间不断增大。
图 9 初始磁场沿y轴正方向时的磁压力及磁能量分布Figure 9. Distribution of magnetic pressure and magnetic energy when initial magnetic field is along y axis图 10为初始磁场沿x轴正方向时磁压力及磁能量的分布。由图 10可知，磁压力与磁能量的分布类似，其峰值均集中在界面附近。与初始磁场沿y轴正方向的情况相似，x方向磁压力沿x轴对称分布，y方向磁压力沿x轴反对称分布，同时磁压力与磁能量的作用区域更紧贴界面。由于磁压力较弱且分布不均匀，图 11给出了边界处磁压力矢量分布。由图 11(a)可知，当磁场平行于来流方向时，初始磁压力分布于界面处，且与剪切力及外界对其压力的合力方向相反。随后，在剪切力的作用下，气柱表面开始卷起涡，磁压力的方向随着界面变形而改变，如图 11(b)所示。在主涡生成时，磁压力较大的位置依在界面附近，此对涡中心部分影响不大，如图 11(c)和图 11(d)所示。由此可见：当涡量大于零时(涡逆时针方向)，磁压力指向涡动中心；而当涡量小于零时，磁压力指向外部。
图 10 初始磁场沿x轴正方向时磁压力及磁能量分布(每幅子图中，左上为磁压力，左下为磁能量，中间和右侧分别为磁压力在x、y方向的分量)Figure 10. Magnetic pressure and magnetic energy when initial magnetic field is along x axis (upper left: magnetic pressure; lower left: magnetic energy; middle: magnetic pressure along x axis; right: magnetic pressure along y axis)
图 11 初始磁场沿x轴正方向时边界处磁压力矢量图Figure 11. Vector diagram of magnetic pressure on the interface when initial magnetic field is along x axis图 12给出了初始磁场沿x轴正方向时磁压力及磁能量沿横向及纵向的分布。由图 12(a)和图 12(b)可知，初始时刻，横轴和纵轴上磁压力及磁能量均小于初始磁场沿y轴正方向的情况。当t=165 μs时，横轴上的磁压力集中在透射激波尾部，磁能量在顶部界面后迅速减小，接近透射激波时快速增大并达到峰值；而纵轴上磁能量和磁压力的峰值均在上、下界面处。当t=200 μs时，横轴上的磁压力集中在气柱尾端，磁能量仍然在透射激波尾部，纵轴上的磁压力不断增大；此时由于涡在界面处开始生成，纵轴上磁压力的分布不对称。当t=370 μs时，透射激波聚集，在尾部形成射流，圆弧形二次激波向外传播，并最终碰到上、下界面发生反射；此时，横轴上磁压力与磁能量最大值位置相同，均在气柱尾部界面处达到峰值，而纵轴上的磁压力峰值则继续增大，如图 12(c)所示。当两个主涡形成后，纵轴处主涡位置的磁能量较为明显，但此时磁能量及磁压力最大位置仍在界面附近；横轴上磁压力峰值减小，磁能量继续增加，并集中在气柱尾部，如图 12(d)所示。
图 12 初始磁场为x轴正方向时磁压力及磁能量分布Figure 12. Distribution of magnetic pressure and magnetic energy when initial magnetic field is along x axis3. 结论
基于CTU+CT算法，对磁场环境下激波冲击重质气体进行数值模拟，并与不加磁场的情况进行对比，所得结论如下。
(1) R22气柱与入射激波作用后，由于KH不稳定性，气柱上、下两侧不断卷起最终形成两个主涡。同时，由于RM不稳定性，后期气柱左端界面处产生了“尖钉”与“气泡”状结构。磁场能够抑制界面不稳定性，平行磁场通过减少界面上涡量的生成减缓不稳定性的增长；垂直磁场则通过包裹界面的磁力线压缩界面，抑制剪切力对界面的作用，使界面处不再卷起，从而抑制界面不稳定性。
(2) 垂直磁场和平行磁场都可以很好地抑制RM不稳定性，对于KH不稳定性，垂直磁场的控制效果更好，不仅可以抑制界面上涡串的卷起，还可以阻止主涡的发展，而平行磁场做不到后者。磁场对射流影响不大，射流处的磁能量可以在一定程度上抑制射流的衰减，同时，垂直磁场可以减小射流聚焦时压力及速度峰值。
(3) 当磁场强度较小(β≈2500)时，初始磁场对激波与R22气柱作用过程没有显著影响，随后磁压力及磁能量不断增加，从而抑制不稳定性的发展。磁压力及磁能量对界面影响很大，在上、下界面处达到峰值。同时，磁压力对界面没有压缩作用，而是产生一个反向作用力减小外界对界面的影响。当涡形成后，磁场对涡量影响不大，但可以通过抑制外界的剪切力阻止其进一步的发展。
-
图 1 各波段电磁波对应的波长和频率
Figure 1. The wavelength and frequency of electromagnetic wave in each band
图 15 材料破碎时产生微波辐射模型
Figure 15. The model of microwave radiation when materials are damaged
图 18 超高速碰撞产生可见光的机理
Figure 18. The mechanism of visible light caused by hypervelocity impact
表 1 铝原子和铝离子的共振线
Table 1. Resonance lines of aluminum atom and aluminum ion
元素 跃迁能级 波长/nm 对基态的能量/eV Al Ⅰ 3s23p2P1/2～3s23d2D3/2 308.215 4.02 Al Ⅰ 3s23p2P1/2～3s24s2S1/2 394.401 3.15 Al Ⅱ 3s21S0～3s3p3P2 265.007 4.67 Al Ⅱ 3s21S0～3s3p3P1 266.916 4.65 
下载: 导出CSV
-
[1] 经福谦. 超高速碰撞现象 [J]. 爆炸与冲击, 1990, 10(3): 279–288.JING F Q. Hypervelocity impact phenomena [J]. Explosion and Shock Waves, 1990, 10(3): 279–288. [2] PERKINS M A, SIMPSON J A, TUZZOLINO A J. A cometary and interplanetary dust experiment on the Vega spacecraft missions to Halley’s Comet [J]. Nuclear Instruments and Methods in Physics Research A, 1985, 239(2): 310–323. DOI: 10.1016/0168-9002(85)90731-4. [3] LEDERER S M, JENSEN E, FANE M, et al. Unveiling clues from spacecraft missions to comets and asteroids through impact experiments [C] // Division for Planetary Sciences Meeting Abstracts, 2016: 331.01 https://ui.adsabs.harvard.edu/abs/ 2016DPS....4833101L. [4] A’HEARN M F, BELTON M J S, DELAMERE A, et al. Deep impact: a large-scale active experiment on a cometary nucleus [J]. Space Science Reviews, 2005, 117(1−2): 1–21. DOI: 10.1007/s11214-005-3387-3. [5] ERNST C M, SCHULTZ P H. Evolution of the deep impact flash: implications for the nucleus surface based on laboratory experiments [J]. Icarus, 2007, 190(2): 334–344. DOI: 10.1016/j.icarus.2007.03.030. [6] CRAWFORD D A, SCHULTZ P H. Laboratory observations of impact−generated magnetic fields [J]. Nature, 1988, 336(6194): 50–52. DOI: 10.1029/91JE02012. [7] ANZ-MEADOR P. Orbital debris quarterly news [R]. Texas: NASA Oribital Debris Program Office, 2020. [8] WU Q, ZHANG Q M, LONG R R, et al. Potential space debris shield structure using impact-initiated energetic materials composed of polytetrafluoroethylene and aluminum [J]. Applied Physics Letters, 2016, 108(10): 101903. DOI: 10.1063/1.4943584. [9] ZHANG P L, GONG Z Z, TIAN D B, et al. Comparison of shielding performance of Al/Mg impedance-graded-material-enhanced and aluminum whipple shields [J]. International Journal of Impact Engineering, 2019, 126: 101–108. DOI: 10.1016/j.ijimpeng.2018.12.007. [10] FLETCHER A, MATHIAS D L, CLOSE S. Susceptibility of spacecraft to impact-induced electromagnetic pulses [C] // 2015 Annual Reliability and Maintainability Symposium. 2015: 1−6. DOI: 10.1109/RAMS.2015.7105118. [11] THOMSON G M, MCNEIR M R. Impact flash: a tool for rapid battle damage assessment [C] // Infrared Technology and Applications XXX. 2004, 5406: 690−700. DOI: 10.1117/12.538654. [12] LAWRENCE R J, REINHART W D, CHHABILDAS L C, et al. Spectral measurements of hypervelocity impact flash [J]. International Journal of Impact Engineering, 2006, 33(1-12): 353–363. DOI: 10.1016/j.ijimpeng.2006.09.010. [13] RUDOLPH M. Review of radio frequency emission from hypervelocity impact plasmas [J]. Procedia Engineering, 2013, 58(1): 409–417. DOI: 10.1016/j.proeng.2013.05.047. [14] SUGITA S, SCHULTZ P H, ADAMS M A. Spectroscopic measurements of vapor clouds due to oblique impacts [J]. Journal of Geophysical Research Planets: E, 1998, 103(8): 19427–19441. DOI: 10.1029/98je02026. [15] TAKANO T, MUROTANI Y, MAKI K, et al. Microwave emission due to hypervelocity impacts and its correlation with mechanical destruction [J]. Journal of Applied Physics, 2002, 92(9): 5550–5554. DOI: 10.1063/1.1513885. [16] MA Z X, SHI A H, LI J L, et al. Radiation mechanism analysis of hypervelocity impact ejecta cloud [J]. International Journal of Impact Engineering, 2020, 141: 103560. DOI: 10.1016/j.ijimpeng.2020.103560. [17] SUGITA S, SCHULTZ P H. Interactions between impact-induced vapor clouds and the ambient atmosphere: 1. spectroscopic observations using diatomic molecular emission [J]. Journal of Geophysical Research: E, 2003, 108(6): 5051. DOI: 10.1029/2002je001959. [18] CLOSE S, LINSCOTT I, LEE N, et al. Detection of electromagnetic pulses produced by hypervelocity micro particle impact plasmas [J]. Physics of Plasmas, 2013, 20(9): 2102. DOI: 10.1063/1.4819777. [19] HEW Y M, GOEL A, CLOSE S, et al. Hypervelocity impact flash and plasma on electrically biased spacecraft surfaces [J]. International Journal of Impact Engineering, 2018, 121: 1–11. DOI: 10.1016/j.ijimpeng.2018.05.008. [20] CLOSE S, KELLEY M, FLETCHER A, et al. RF signatures of hypervelocity impacts on spacecraft [C] // 3rd AIAA Atmospheric Space Environments Conference. 2011: 3150. DOI: 10.2514/6.2011-3150. [21] SUGITA S, SCHULTZ P H. Spectroscopic characterization of hypervelocity jetting: comparison with a standard theory [J]. Journal of Geophysical Research Planets: E, 1999, 104(12): 30825–30845. DOI: 10.1029/1999JE001061. [22] TAKANO T, IKEDA H, MAEDA T. Consideration of the mechanism of microwave emission due to material destruction [J]. Journal of Applied Physics, 2010, 108(8): 083722. DOI: 10.1063/1.3499291. [23] JEAN B, ROLLINS T L. Radiation from hypervelocity impact generated plasma [J]. AIAA Journal, 2012, 8(10): 1742–1748. DOI: 10.2514/3.5984. [24] MA Z X, SHI A H, LI J L, et al. Radiation evolution characteristics of the ejecta cloud produced by aluminum projectiles hypervelocity impacting aluminum plates [J]. International Journal of Impact Engineering, 2020, 138: 103480. DOI: 10.1016/j.ijimpeng.2019.103480. [25] MA Z X, HUANG J, SHI A H, et al. Analysis technique for ejecta cloud temperature using atomic spectrum [J]. International Journal of Impact Engineering, 2016, 91: 25–33. DOI: 10.1016/j.ijimpeng.2015.12.008. [26] MA Z X, HUANG J, SHI A H, et al. The analysis technique for ejecta cloud temperature based on atomic spectrum [J]. Procedia Engineering, 2015, 103: 357–364. DOI: 10.1016/j.proeng.2015.04.033. [27] ZHANG K, ZHANG Q M, LONG R R. The experiment study on flash spectrum produced by hypervelocity impact [J]. Applied Mechanics & Materials, 2015, 782: 197–203. DOI: 10.4028/www.scientific.net/AMM.782.197. [28] XUE Y J, ZHANG Q M, LIU D Y, et al. Hypersonic impact flash characteristics of a long-rod projectile collision with a thin plate target [J]. Defence Technology, 2021, 17(2): 368–376. DOI: 10.1016/j.dt.2020.02.011. [29] ZHANG Q M, GONG L F, MA Y F, et al. The electromagnetic properties of plasma produced by hypervelocity impact [J]. Physics of Plasmas, 2018, 25(2): 022906. DOI: 10.1063/1.5009067. [30] LI H W, HAN J W, CAI M H, et al. The analogue experiment of small space debris impact inducing solar array discharge [J]. International Journal of Impact Engineering, 2020, 143: 103582. DOI: 10.1016/j.ijimpeng.2020.103582. [31] 张凯. 超高速碰撞LY12铝靶产生电磁辐射实验研究[D]. 北京: 北京理工大学, 2016. DOI: CNKI:CDMD:1.1018.811986. [32] ZHANG K, LONG R R, ZHANG Q M, et al. Flash characteristics of plasma induced by hypervelocity impact [J]. Physics of Plasmas, 2016, 23(8): 083519. DOI: 10.1063/1.4960297. [33] LAWRENCE R J, REINHART W D, CHHABILDAS L C, et al. Hypervelocity impact flash at 6, 11, and 25 KM/S [C] // AIP Conference Proceedings. 2006, 845(1): 1349−1352. DOI: 10.1063/1.2263574. [34] DUGGER P, HENDRIX R. Measurements of transient hypervelocity impact phenomena at the AEDC [C] // 32nd Aerospace Sciences Meeting and Exhibit. 1994: 87. DOI: 10.2514/6.1994-87. [35] REINHART W D, THORNHILL T F, CHHABILDAS L C, et al. Temperature measurements of expansion products from shock compressed materials using high-speed spectroscopy [J]. International Journal of Impact Engineering, 2008, 35(12): 1745–1755. DOI: 10.1016/j.ijimpeng.2008.07.062. [36] 石安华, 柳森, 黄洁, 等. 铝弹丸超高速撞击铝靶光谱辐射特性实验研究 [J]. 宇航学报, 2008, 29(2): 715–717. DOI: 10.3873/j.issn.1000-1328.2008.02.061.SHI A H, LIU S, HUANG J, et al. Spectra measurement of radiation produced by aluminum projectiles impacting aluminum targets at hypervelocity speeds [J]. Journal of Astronautics, 2008, 29(2): 715–717. DOI: 10.3873/j.issn.1000-1328.2008.02.061. [37] 石安华, 柳森, 黄洁, 等. 超高速碰撞可见光谱辐射强度测量技术 [J]. 实验流体力学, 2007, 21(4): 83–85. DOI: 10.3969/j.issn.1672-9897.2007.04.017.SHI A H, LIU S, HUANG J, et al. Visible spectral radiant intensity measuring technology for hypervelocity impact phenomena [J]. Journal of Experiments in Fluid Mechanics, 2007, 21(4): 83–85. DOI: 10.3969/j.issn.1672-9897.2007.04.017. [38] MEYER-VERNET N, MAKSIMOVIC M, CZECHOWSKI A, et al. Dust detection by the wave instrument on STEREO: nanoparticles picked up by the solar wind? [J]. Solar Physics, 2009, 256(1−2): 463–474. DOI: 10.1007/s11207-009-9349-2. [39] MEYER-VERNET N, LECACHEUX A, KAISER M L, et al. Detecting nanoparticles at radio frequencies: Jovian dust stream impacts on Cassini/RPWS [J]. Geophysical Research Letters, 2009, 36(3). DOI: 10.1029/2008GL036752. [40] GURNETT D A, KURTH W S, KIRCHNER D L, et al. The Cassini radio and plasma wave investigation [J]. Space Science Reviews, 2004, 114(1−4): 395–463. DOI: 10.1007/s11214-004-1434-0. [41] WARWICK J W, PEARCE J B, EVANS D R, et al. Planetary radio astronomy observations from Voyager 1 near Saturn [J]. Science, 1981, 212(4491): 239–243. DOI: 10.1126/science.212.4491.239. [42] WARWICK J W, EVANS D R, ROMIG J H, et al. Planetary radio astronomy observations from Voyager 2 near Saturn [J]. Science, 1982, 215(4532): 582–587. DOI: 10.1126/science.215.4532.582. [43] MAKI K, TAKANO T, FUJIWARA A, et al. Radio-wave emission due to hypervelocity impacts in relation to optical observation and projectile speed [J]. Advances in Space Research, 2004, 34(5): 1085–1089. DOI: 10.1016/j.asr.2003.02.032. [44] MAKI K, SOMA E, TAKANO T, et al. Dependence of microwave emissions from hypervelocity impacts on the target material [J]. Journal of Applied Physics, 2005, 97(10): 104911. DOI: 10.1063/1.1896092. [45] THOMA K, SCHÄFER F, HIERMAIER S, et al. An approach to achieve progress in spacecraft shielding [J]. Advances in Space Research, 2004, 34(5): 1063–1075. DOI: 10.1016/j.asr.2003.03.034. [46] HERMALYN B, SCHULTZ P H, SHIRLEY M, et al. Scouring the surface: ejecta dynamics and the LCROSS impact event [J]. Icarus, 2012, 218(1): 654–665. DOI: 10.1016/j.icarus.2011.12.025. [47] ERNST C M, SCHULTZ P H. Effect of initial conditions on impact flash decay [C] // Lunar and Planetary Science Conference. 2003: 2020. https://ui.adsabs.harvard.edu/abs/2003LPI....34.2020E. [48] ERNST C M, SCHULTZ P H. Effect of velocity and angle on light intensity generated by hypervelocity [C] // Proceedings, 33rd Lunar and Planetary Science Conference. 2002: 1782. https://ui.adsabs.harvard.edu/abs/2002LPI....33.1782E. [49] BAIRD J K, HOUGH G R, KING T R. Velocity dependence of impact fluorescence [J]. International Journal Of Impact Engineering, 1997, 19(3): 273–276. DOI: 10.1016/S0734-743X(96)00039-5. [50] BURCHELL M J, KAY L, RATCLIFF P R. Use of combined light flash and plasma measurements to study hypervelocity impact processes [J]. Advances in Space Research, 1996, 17(12): 141–145. DOI: 10.1016/0273-1177(95)00772-7. [51] MACCORMACK R W. Investigation of impact flash at low ambient pressures [C] // 6th Symposium on Hypervelocity Impact. 1963: 613−625. [52] EICHHORN G. Analysis of the hypervelocity impact process from impact flash measurements [J]. Planetary and Space Science, 1976, 24(8): 771–781. DOI: 10.1016/0032-0633(76)90114-8. [53] MISRA A. Electromagnetic effects at metallic fracture [J]. Nature, 1975, 254(5496): 133–134. DOI: 10.1038/254133a0. [54] FLETCHER A, CLOSE S, MATHIAS D. Simulating plasma production from hypervelocity impacts [J]. Physics of Plasmas, 2015, 22(9): 093504. DOI: 10.1063/1.4930281. [55] JU Y Y, ZHANG Q M, ZHANG D J, et al. Theoretical model for plasma expansion generated by hypervelocity impact [J]. Physics of Plasmas, 2014, 21(9): 092112. DOI: 10.1063/1.4895592. [56] ZEL’DOVICH Y B, RAIZER Y P. Physics of shock waves and high-temperature hydrodynamic phenomena [M]. Chicago: Courier Corporation, 2002. [57] GHISELLINI G. Radiative processes in high energy astrophysics [M]. New York: Springer, 2012: 873. DOI: 10.1007/978-3-319-00612-3. 期刊类型引用(7)
1. 赵浩霖，刘吉，武锦辉，李绍栋，郭大千，苏海鑫. 高速碰撞金属靶闪光辐射特征研究. 激光杂志. 2025(01): 46-49 . 
百度学术2. 司宇，刘吉，武锦辉，赵磊，闫晓燕. 基于闪光光谱的侵彻过程光学观察窗口分析. 光谱学与光谱分析. 2023(03): 718-723 . 百度学术3. 巨圆圆，张磊，张庆明. 超高速碰撞产生等离子体的机理分析. 兵器材料科学与工程. 2023(02): 19-24 . 百度学术4. 巨圆圆，张磊，张庆明. 超高速碰撞产生等离子体及电磁辐射效应的研究进展. 科学技术与工程. 2023(25): 10605-10616 . 百度学术5. 巨圆圆，张磊，张庆明. 超高速碰撞力电耦合的模型分析. 兵器材料科学与工程. 2023(06): 25-30 . 百度学术6. 刘芹芹，申旭辉，周炳红，李明涛，张庆明，徐文杰，龚自正. 近地小行星撞遇地球灾害链效应. 空间碎片研究. 2022(01): 1-17 . 百度学术7. 刘立恒，马兆侠，杜雪飞，柳森，范晓樯，石安华. 超高速撞击RDX光辐射特征测量试验. 兵器装备工程学报. 2022(09): 291-298 . 百度学术其他类型引用(0)
-
下载:
百度学术
计量
- 文章访问数: 1602
- HTML全文浏览量: 503
- PDF下载量: 139
- 被引次数: 7