-
强爆炸光辐射是一种重要的热毁伤效应源[1-4]。大气层强爆炸过程中光辐射所占能量接近35%,是目标毁伤过程不可忽视的因素。每平方厘米百焦强度的光辐射会导致有机材料等易燃物发生燃烧等显著的热破坏效应[5-6]。对于金属、陶瓷、复合材料等热稳定性较好的材料,虽然光辐射不一定产生宏观的热破坏效应,但由于材料温升导致的力学性能下降,会对目标后续力学效应带来显著影响。
从物理过程来说,材料在光辐射环境下首先发生能量沉积,随之产生热效应,因而材料的光辐射能量沉积是决定光辐射效应的主要因素。通常用能量耦合特性来表征光与物质相互作用中能量沉积的大小[7-8]。在强激光辐照效应研究中,材料的耦合系数主要决定于单色光波长[9-12]。而对于光辐射来说,由于其光谱涵盖紫外、可见和红外部分(约0.3~4.5 µm),光谱耦合特性与实际光辐射谱分布密切相关。现阶段有关光辐射理论研究中[1-3],以太阳光谱或6000 K等效黑体谱作为光辐射的近似光谱,这种光谱近似对于材料光辐射耦合特性的影响暂无定量分析;文献[13-14]中采用0.7和0.2分别作为钢和有机材料光辐射吸收系数,计算了目标光辐射热效应,但并未说明吸收系数的来源及其可靠性。
本文主要针对强爆炸光辐射作用下材料的能量耦合特性开展研究,构建了完整的辐射源参数计算、大气传输计算物理模型,获取了目标位置处光辐射能谱特征。在此基础上,通过材料光谱反射率测量,结合光辐射平均吸收系数计算方法给出了几类材料的能量耦合系数。定量分析了采用实际光辐射和等效黑体辐射所得到的耦合系数差异,并进一步给出了材料耦合系数随强爆炸当量及目标距离的变化关系,为光辐射效应研究提供有力支撑。
-
从能量沉积过程来看,光辐射在材料中主要发生反射、吸收和透射3个过程。对于波长为λ的单色光,当材料为不透明介质时,材料吸收率Aλ、反射率Rλ满足以下关系:
$${A_\lambda } + {R_\lambda } = 1\;\;\;\;.$$ (1) 对于实际光辐射而言,其光谱范围较宽,涵盖紫外、可见和红外部分,因此光辐射的耦合系数与谱分布密切相关。定义平均吸收系数为:
$$A = \frac{{\displaystyle\sum {{A_\lambda }{D_\lambda }{\rm{d}}\lambda } }}{{\displaystyle\sum {{D_\lambda }{\rm{d}}\lambda } }} = \frac{{\displaystyle\sum {(1 - {R_\lambda }){D_\lambda }{\rm{d}}\lambda } }}{{\displaystyle\sum {{D_\lambda }{\rm{d}}\lambda } }}\;\;\;\;,$$ (2) 其中,Dλ为光谱中波长为λ的能量比例,当其为实际光辐射谱分布时,A即光辐射耦合系数。显然,要获取材料的光辐射耦合系数,必须同时给出光谱分布和材料的光谱反射率。
-
目标位置的光谱分布取决于两个因素:一是辐射源的特征,也即爆炸火球的辐射特征,其从根本上决定了光源的初始参数;另一个是大气传输。由于光辐射属于连续谱,涵盖波长范围较宽,其在大气中传输时不同波长的光辐射衰减特征有所不同,导致目标位置处的光辐射与初始辐射特征出现差异。为获得目标位置的光谱分布,需要针对上述两个过程分别开展。
-
火球是强爆炸光辐射的辐射源,其发展演化包含了辐射输运与流体强耦合的过程,因而,针对火球发展的数值模拟,必须建立在辐射流体力学[15]的基础上。基于灰体近似[16-17]的处理方法给出了火球发展的半径、辐射总强度等参数规律,但对于能谱特征来说,灰体近似遇到了本质困难。本文构建了多群辐射流体力学模型[18-19],为获取辐射源能谱特征奠定基础。
在局域热动力平衡(LTE)假定下,描述强爆炸火球发展的多群辐射流体力学方程,如式(3)所示:
$$\left\{ \begin{aligned} & \dfrac{{\partial \rho }}{{\partial t}} + \nabla \cdot (\rho {{v}}) = 0 \\ & \dfrac{\partial }{{\partial t}}(\rho {{v}}) + \nabla p = - {{{M}}^1} \\ & \dfrac{\partial }{{\partial t}}\left(\dfrac{1}{2}\rho {{{v}}^2} + \rho {e^I}\right) + \nabla \cdot \left[\left(\dfrac{1}{2}\rho {{{v}}^2} + \rho {e^I} + p\right){{v}}\right] = - {M^0} \\ & {M^0} = \displaystyle\sum\limits_g {M_g^0} = \displaystyle\sum\limits_g {\dfrac{{\partial {E_g}}}{{\partial t}} + \nabla \cdot {{{F}}_g}} \\ & {{{M}}^1} = \displaystyle\sum\limits_g {{{M}}_g^1} = \sum\limits_g {\dfrac{1}{{{c^2}}}\dfrac{{\partial {{{F}}_g}}}{{\partial t}} + \nabla \cdot {{{P}}_g}} \end{aligned} \right.,$$ (3) 其中,ρ为空气密度,v为空气速度,p为空气压强,eI为空气内能。方程后两项为各群内光辐射输运的总能量(M0)和总动量(M1),下标g为1,2,3...即分群数目。Eg、Fg、Pg、M0g、M1g分别为第g群辐射能密度、辐射能流、辐射压强张量及光辐射输运能量和动量。具体计算如式(4)所示:
$$\left\{\!\! \begin{array}{l} \dfrac{{\partial {E_g}}}{{\partial t}}\;\; + \nabla \cdot {{{F}}_g} = c{\kappa _g}\displaystyle\int\limits_{{\gamma _{g - 1}}}^{{\gamma _g}} {\left[ {\dfrac{{4{\text{π }}}}{c}B\left( \gamma \right) - E(\gamma )} \right]} {\rm{d}}\gamma + \dfrac{{{v}}}{c}{\kappa _g}\displaystyle\int\limits_{{\gamma _{g - 1}}}^{{\gamma _g}} {{{F}}(\gamma ){\rm{d}}\gamma } \\ \quad \quad \quad \quad \quad \;\;\, = {\rm{c}}{\kappa _g}\left( {{B_g} - {E_g}} \right) + \dfrac{{{\kappa _g}}}{c}{{v}} \cdot {{{F}}_g} \\ \dfrac{1}{{{c^2}}}\dfrac{{\partial {{{F}}_g}}}{{\partial t}} + \nabla \cdot {{{P}}_g} = - \dfrac{1}{c}{\kappa _g}\displaystyle\int\limits_{{\gamma _{g - 1}}}^{{\gamma _g}} {{{F}}(\gamma ){\rm{d}}\gamma } + \dfrac{{4{\text{π }}}}{{{c^2}}}{\kappa _g}{{v}}\displaystyle\int\limits_{{\gamma _{g - 1}}}^{{\gamma _g}} {B\left( \gamma \right)} {\rm{d}}\gamma + \dfrac{{{v}}}{c}{\kappa _g}\displaystyle\int\limits_{{\gamma _{g - 1}}}^{{\gamma _g}} {{{P}}(\gamma )} {\rm{d}}\gamma \\ \quad \quad \quad \quad \quad \quad\, = \dfrac{{{\kappa _g}}}{c}\left( { - {{{F}}_g} + {{v}}{B_g} + {{v}} \cdot {{{P}}_g}} \right) \end{array} \right.\quad,$$ (4) 其中,c为光速,kg为第g群吸收系数,γ为光子频率。由于各群光子能量具有上限和下限,因此采用γg、γg-1分别标记第g群光子频率的上限和下限。Eγ、Fγ和Pγ分别为辐射能密度、辐射能流及辐射压强张量函数,B(γ)为黑体辐射谱分布,
${B_g} = \dfrac{{4{\text{π}} }}{c}\displaystyle\int\limits_{{\gamma _{{{g - 1}}}}}^{{\gamma _{{g}}}} {B(\gamma )} {\rm{d}}\gamma$ 为黑体辐射下第g群辐射能密度。空气的状态方程采用实际空气状态方程[17]。光辐射按照能量分为21群,不同群(group, g)所对应的光子能量(photon energy,pe)范围见表1。
表 1 不同群内光子能量范围(21群)
Table 1. Photon energy in different groups (21)
g 1 2 3 4 5 6 7 8 9 10 11 12 13 pe (eV) 0.01−0.5 0.5−1.0 1.0−1.8 1.8−2.1 2.1−2.5 2.5−3.1 3.1−4.0 4.0−7.0 7.0−10 10−20 20−40 40−70 70−100 g 14 15 16 17 18 19 20 21 pe (eV) 100−200 200−400 400−1000 1000−2 000 2 000−5000 5000−10000 10000−20000 20000−80000 针对上述模型,采用算子分裂的数值求解方法,直接获得火球半径rf(t)以及不同波长范围的辐射能流Fg。火球辐射过程中,光谱主要集中在0.3~4.0 µm波段,对应于2~8群的辐射能量。由于分群数目及其对应的辐射参数(kg)限制,计算出的光谱分布特征不够精细。考虑到火球在整个后期发展阶段[1-2]其光谱与黑体辐射接近,通过式(5)计算得到火球的有效温度,从而对其能谱特征进行细致表征。10 kt当量下的火球半径及有效温度计算结果见图1。
$$\sigma T_e^4 = \sum\limits_g {{{{F}}_g}({r_f})} \quad\quad.$$ (5) -
光辐射从火球表面向整个空间传输,在辐射源近似为球面光源的情况下,距离爆心距离R处、波长为λ的光辐射强度为[20-21]:
$${E_\lambda }{\rm{ = }}\frac{{P_\lambda ^{{\rm{total}}}}}{{4{\text{π}} {R^2}}}{\tau _\lambda },$$ (6) 其中,Eλ是波长为λ、距离爆点R处的光辐射通量,τλ为对应波长下的大气透过率。Pλtotal为火球表面辐射出的波长λ附近的光辐射总功率。当火球半径为rf(t)、有效温度为Te时,其计算方法为:
$$\begin{split} P_\lambda ^{{\rm{total}}} =&\, \displaystyle\int_{{t_0}}^{{t_\infty }} {4{\text{π}} r_f^2(t){P_\lambda }} (t){\rm{d}}t \\ =&\, \displaystyle\int_{{t_0}}^{{t_\infty }} {4{\text{π}} r_f^2(t)\Biggl[\displaystyle\int_{{\lambda _1}}^{{\lambda _2}} {\dfrac{{2{\hbar }{c^2}}}{{{\lambda ^5}{e^{{\hbar }c/k\lambda {T_e}}} - 1}}{\rm{d}}\lambda \Biggl]} } {\rm{d}}t, \end{split} $$ (7) 其中λ1和λ2为包含波长λ的某个积分区间;Pλ(t)为对应区间的光辐射强度,
$\hbar$ 为普朗克常数,k为波尔兹曼常数。t0和t∞分别为光辐射的初始和结束时间。不同波长下的大气透过率τλ对光辐射传输具有重要影响。严格计算大气透过率,需要依据吸收和散射理论[22]分别计算大气各粒子对光的衰减效应。早期工程计算采用与波长无关的大气透过率经验计算方法[1]。为提高模型计算可信度,采用MODTRAN大气传输模块实现不同大气条件下、具有一定光谱精度的大气透过率计算[23-24]。
-
利用上述方法可以完整计算出爆炸条件下的光辐射谱分布特征。大气条件对光传输具有较大影响,为了便于与实测数据做对比分析,本文大气模式统一为中纬度夏季,能见度取23 km,采用沙漠气溶胶类型,爆点位于1 km高度、测点高度0 km,其他大气参数按照MODTRAN标准大气模式进行设定。由于缺少全谱范围内的细致测量结果,依据光辐射在紫外(小于0.4 μm)、可见(0.4~0.76 μm)和红外(大于0.76 μm)部分所占比例的测量结果[1, 3],对光谱特征计算方法进行验证。计算结果与实测对比见图2。由图2可以看出:除个别距离(约10 km和25 km)处可见和红外比例数据与实测数据有一定偏离外,其他总体符合的都较好。
图 2 光辐射可见、红外和紫外波段能量比例的实测结果与计算结果的对比
Figure 2. Comparison of measured and calculated energy ratios in visible, infrared and ultraviolet bands of thermal radiation
进一步根据该方法分析不同当量下光辐射谱特征的变化规律。图3给出了当量在1 kt~2 000 kt下距离爆心投影点1 km处的归一化光辐射能谱。由图3可以看出:随着当量的增加,光辐射强度峰值所对应的中心波长(λM)有一定红移,相应地短波方向的光能量有所减少,而长波方向的光辐射能量则逐渐增加。
根据维恩位移定律:
${\lambda _M} = \dfrac{{0.29 \times {{10}^{ - 2}}}}{T}$ 知,中心波长越长、对应的等效黑体温度越低,也即意味着:随着爆炸当量的增加,相同距离处的光辐射能谱所对应的等效温度是降低的。表2列举出了几种当量下光辐射中心波长及其对应的等效黑体温度。对比来看:从20 kt小当量变化到2000 kt大当量时,等效温度变化约为11%。图 3 不同当量下距爆心投影点1 km的归一化光辐射能谱
Figure 3. Normalized spectral distribution under different explosion yields at a range of 1 km from the center of burst projection point
表 2 不同当量下光辐射中心波长及其对应的等效黑体温度
Table 2. Central wavelength of thermal radiation under different explosion yields and their equivalent blackbody temperatures
当量/kt 20 100 2000 中心波长/µm 0.46 0.48 0.52 等效温度/K 6.3×103 6.0×103 5.6×103 进一步考虑爆心距离对于光谱的影响。图4给出了20 kt下不同距离处光辐射归一化能谱分布图。由图4可以看出:随着距离的增加,光辐射能谱的变化趋势与当量增加时相似。从能量构成来看:距离越远,短波范围的光辐射能量越小,而长波方向的光辐射能量越多,其主要原因还在于大气对短波吸收较强,而长波衰减较小。
-
材料反射率的测量主要依靠积分球,积分球测量原理[9, 12]及测量光路见图5。
图 5 材料光谱反射率测量原理及示意图
Figure 5. Principle and schematic diagram of the spectral reflectivity measurement system for different materials
实验中波长范围为0.2~2.5 µm,涵盖了光辐射的主要光谱区域,能量约占光辐射总能量的97%以上。上述范围的光谱反射率测量能够反映材料在光辐射下的能量耦合特性。选取常用的金属(铝、铜和银)、陶瓷(TiO2)以及碳纤维环氧复合材料进行分析,反射率测量结果见图6。
-
由铝材料反射率计算得到的吸收系数和光辐射能谱分布见图7。光谱计算中波长间隔dλ为0.05 µm,因此按照0.05 µm的等波长间距,将光谱范围进行重新划分,相应的材料吸收率取对应波长范围内的平均值。进一步根据式(2),计算得到铝材料的光辐射耦合系数为0.19。
图 7 铝材料光谱吸收率及光辐射能谱分布
Figure 7. Spectral absorptivity of aluminum material and its thermal spectrum distribution
按照同样方法,计算了其他几类材料的光辐射耦合系数,见表3。作为光辐射耦合特性研究的近似方法,通常将光辐射能谱等效为6000 K黑体辐射(与太阳光谱相近)。为定量分析上述光谱近似所带来的差异,对比分析了分别采用实际光谱和等效黑体辐射谱计算的耦合系数,结果见表3。
表 3 采用实际光辐射和6000 K黑体辐射谱计算的耦合系数
Table 3. Coupling coefficients calculated from the actual thermal radiation spectrum and 6000 K black-body radiation spectrum
材料类型 光辐射耦合系数 实际光谱 6000 K黑体辐射 相差(%) 金属 Al 0.19 0.17 −10% Cu 0.28 0.25 −10% Ag 0.21 0.18 −14% 陶瓷 TiO2 0.25 0.23 −8% 复合材料 C/E 0.92 0.89 −3.3% 可以看出:对于所列举的几类材料,金属、陶瓷材料的耦合系数相对较小,这是由于这些材料在整个可见、红外范围内都具有较高的反射率所造成的,而碳纤维环氧复合材料的耦合系数可达0.92。利用6000 K黑体辐射计算的耦合系数要比实际光谱的结果低,最大相差可达14%(Ag材料)。其原因是由于实际光辐射比6000 K等效黑体谱在可见光以上波长(大于0.4 µm)部分的能量更高(见表4),长波方向的能量占据了光辐射和黑体谱的绝大部分,也决定了材料光谱耦合特性的主要特征。
表 4 实际光辐射以及6000 K黑体在不同光谱区间的能量份额
Table 4. Energy proportion in different spectral intervals for actual thermal radiation and 6000 K blackbody
所占能量比例(%) <0.4 µm 0.4~0.76 µm >0.76 µm 光辐射 11 40 49 6000 K黑体辐射 14 43 43 以光辐射辐照下铝材料的实测温升为例,校验光辐射耦合系数计算的准确性。材料温度计算模型主要是傅立叶热传导方程[25-26],铝板厚度为0.8 mm,在约140 J/cm2的辐照下,背表面温升实测结果与采用耦合系数分别为0.19和0.17的计算结果对比见图8。由图8可以看出:计算给出的材料温升与实测数据总体较为符合;相对而言,采用实际光辐射耦合系数(0.19)的计算结果比6000 K黑体谱耦合系数(0.17)更加接近实测数据。
图 8 铝材料背表面温升实测结果与采用耦合系数分别为0.19和0.17的计算结果对比
Figure 8. Comparison of the measured temperature rise on the back surface of aluminum material with the calculated results under coupling coefficients of 0.19 and 0.17
如前所述,光辐射谱特征是与爆炸当量及爆心距离有关的。下面将以铝材料为例,分析光辐射耦合系数随当量及距离的变化情况。图9分别给出了材料耦合系数随当量的变化情况(距离为1 km)及随爆心距离(当量20 kt)的变化情况。
图 9 铝材料光辐射耦合系数与当量(a)及爆心距离(b)的关系
Figure 9. Relationship between the coupling coefficients of aluminum material and the explosion yield (a) and burst center distance (b)
由图9可以看出:铝材料光辐射耦合系数随当量以及爆心距离增加均逐渐减小,但总体变化幅度不大。当量从1 kt增加到2 000 kt,平均吸收系数从0.201减小到0.180,减小约10%;距离从1 km增加到9 km,平均吸收系数从0.193减小到0.169,减小约13%。对于同一材料而言,平均吸收系数仅取决于光谱特征,而光谱特征随当量及距离的变化趋势是相似的(见图3和图4)。正是由于这一点,导致材料耦合系数随当量和距离增加呈现出相同的变化趋势。
-
本文提出了一套完整的光辐射耦合系数计算方法,获取了目标位置处的光谱特征,计算给出了几类典型材料的光辐射耦合吸收系数。
研究发现:金属、陶瓷材料在可见和红外部分受高反射率的影响,它们的耦合系数相对较小,而碳纤维环氧复合材料的耦合系数可达0.92。采用实际光辐射能谱计算的耦合系数比6000 K黑体谱的结果要高,最大约为14%。
由于光辐射谱特征与爆炸当量及爆心距离相关,铝材料光辐射耦合系数随当量及爆心距离增加表现出逐渐减小的趋势,但总体变化幅度不大。
Energy coupling characteristic of materials under thermal radiation produced by strong explosion
-
摘要: 为了获取强爆炸光辐射作用下材料的能量耦合特性,发展了强爆炸辐射源参数以及光辐射传输的物理模型和计算方法,计算给出了不同条件下目标位置处的光辐射谱特征。利用材料光谱反射率测量方法,结合光辐射耦合系数计算方法获取了几类材料的能量耦合系数。结果显示:金属、陶瓷材料的光辐射耦合系数相对较小,而碳纤维环氧复合材料的耦合系数可达0.92;采用实际光辐射能谱计算的耦合系数比近似6000 K黑体谱的结果要高,最大约14%。以铝材料为例,光辐射耦合系数随当量及爆心距离增加均表现出逐渐减小的趋势,但总体变化幅度不大。Abstract: To obtain the energy coupling characteristic of materials under strong explosive thermal radiation, a physical model for calculating radiation source parameters and atmospheric transmission is constructed, and the characteristics of the radiation spectrum at the target location under different conditions are obtained. The energy coupling coefficients of several kinds of materials are produced by spectral reflectance measurement and by calculating the average absorption coefficient of thermal radiation. The coupling coefficients of metal and ceramic materials are relatively small while it can be as high as 0.92 for carbon fiber epoxy composites. The coupling coefficient calculated from the actual thermal radiation spectrum is higher than that calculated from 6000 K blackbody radiation spectrum, and the maximum difference is about 14%. Taking aluminum material as an example, the coupling coefficient of thermal radiation decreases gradually with the increase of explosion yield and distance, but the overall variation is small.
-
表 1 不同群内光子能量范围(21群)
Table 1. Photon energy in different groups (21)
g 1 2 3 4 5 6 7 8 9 10 11 12 13 pe (eV) 0.01−0.5 0.5−1.0 1.0−1.8 1.8−2.1 2.1−2.5 2.5−3.1 3.1−4.0 4.0−7.0 7.0−10 10−20 20−40 40−70 70−100 g 14 15 16 17 18 19 20 21 pe (eV) 100−200 200−400 400−1000 1000−2 000 2 000−5000 5000−10000 10000−20000 20000−80000 表 2 不同当量下光辐射中心波长及其对应的等效黑体温度
Table 2. Central wavelength of thermal radiation under different explosion yields and their equivalent blackbody temperatures
当量/kt 20 100 2000 中心波长/µm 0.46 0.48 0.52 等效温度/K 6.3×103 6.0×103 5.6×103 表 3 采用实际光辐射和6000 K黑体辐射谱计算的耦合系数
Table 3. Coupling coefficients calculated from the actual thermal radiation spectrum and 6000 K black-body radiation spectrum
材料类型 光辐射耦合系数 实际光谱 6000 K黑体辐射 相差(%) 金属 Al 0.19 0.17 −10% Cu 0.28 0.25 −10% Ag 0.21 0.18 −14% 陶瓷 TiO2 0.25 0.23 −8% 复合材料 C/E 0.92 0.89 −3.3% 表 4 实际光辐射以及6000 K黑体在不同光谱区间的能量份额
Table 4. Energy proportion in different spectral intervals for actual thermal radiation and 6000 K blackbody
所占能量比例(%) <0.4 µm 0.4~0.76 µm >0.76 µm 光辐射 11 40 49 6000 K黑体辐射 14 43 43 -
[1] 乔登江. 核爆炸物理概论[M]. 北京: 国防工业出版社, 2003. QIAO D J. Conspectus on Physics of Nuclear Weapon[M]. Beijing: National Defense Industry Press, 2003. (in Chinese) [2] BRODE H L, HILLENDAHL R W, LANDSHOFF R K. Thermal radiation phenomena. Volume V. Radiation hydrodynamics of high temperature air[R]. Sunnyvale: Lockheed Missiles and Space Co., 1967. [3] GLASSTONE S, DOLAN P J. The Effects of Nuclear Weapons[M]. 3rd ed. Washington, D.C.: United State Department of Defense and the United States Department of Energy, 1977. [4] 肇文丽, 张恩山, 李润新. 星载核爆炸光辐射能量分析[C]. 第十届全国核电子学与核探测技术学术年会论文集, 中国核学会, 中国电子学会, 2000. ZHAO W L, ZHANG E SH, LI R X. The optical radiation energy analysis of nuclear explosion in the satellite[C]. Proceedings of the Tenth National Conference on Nuclear Electronics & Nuclear Exploration Technology, China Nuclear Society, China Electronics Society, 2000. (in Chinese) [5] MARRS R E, MOSS W C, WHITLOCK B. Thermal radiation from nuclear detonations in urban environments[R]. Livermore: Lawrence Livermore National Lab., 2007. [6] 庞磊, 张奇. 无约束气云爆炸热辐射伤害效应研究[J]. 北京理工大学学报,2010,30(10):1147-1150. PANG L, ZHANG Q. Study into injury effect of thermal radiation from unconfined vapor cloud explosion[J]. Transactions of Beijing Institute of Technology, 2010, 30(10): 1147-1150. (in Chinese [7] 孙承伟. 激光辐照效应[M]. 北京: 国防工业出版社, 2002. SUN CH W. Laser Irradiation Effect[M]. Beijing: National Defense Industry Press, 2002. (in Chinese) [8] 刘峰, 王立君, 王玉恒, 等. 金属材料表面激光耦合系数的反演[J]. 光学 精密工程,2011,19(2):421-428. doi: 10.3788/OPE.20111902.0421 LIU F, WANG L J, WANG Y H, et al. Inversion for laser coupling coefficient on metal material surfaces[J]. Optics and Precision Engineering, 2011, 19(2): 421-428. (in Chinese doi: 10.3788/OPE.20111902.0421 [9] 吴丽雄, 叶锡生, 王立君, 等. 激光作用下环氧/硅树脂复合涂层的1.3 μm反射特性研究[J]. 光学学报,2011,31(8):0814002. doi: 10.3788/AOS201131.0814002 WU L X, YE X SH, WANG L J, et al. Reflectance of epoxy/organsilicone composite coating at 1.3 μm under laser irradiation[J]. Acta Optica Sinica, 2011, 31(8): 0814002. (in Chinese doi: 10.3788/AOS201131.0814002 [10] 齐立涛. 真空条件下不同波长固体激光烧蚀单晶硅的实验研究[J]. 中国光学,2014,7(3):442-448. QI L T. Different wavelength solid-state laser ablation of silicon wafer in vacuum[J]. Chinese Optics, 2014, 7(3): 442-448. (in Chinese [11] 王贵兵, 刘仓理. 芳纶纤维复合材料对激光的吸收特性研究[J]. 强激光与粒子束,2003,15(11):1065-1067. WANG G B, LIU C L. Absorption characteristics of the composite reinforced by polyaryl amide fibers irradiated by high power laser[J]. High Power Laser and Particle Beams, 2003, 15(11): 1065-1067. (in Chinese [12] 张翔宇. 激光辐照下金属样品的能量耦合特性研究[D]. 长沙: 国防科技大学, 2016. ZHANG X Y. Research on the energy coupling characteristic of metal samples under laser irradiation[D]. Changsha: National University of Defense Technology, 2016. (in Chinese) [13] 绪梅, 朱继, 王良厚, 等. 核爆炸光辐射对物体热效应的数值模拟[C]. 第十七届全国核电子学与核探测技术学术年会论文集, 中国电子学会, 2014. XU M, ZHU J, WANG L H, et al.. Thermal effect simulation of nuclear explosion light radiation on objects[C]. Proceedings of the Seventeenth National Conference on Nuclear Electronics and Nuclear Exploration Technology, Gansu Lanzhou, 2014. (in Chinese) [14] 关志远, 绪梅, 徐鸿明. 橡胶材料的核爆炸热毁伤效应分析[J]. 工兵装备研究,2019,38(3):45-47. GUAN ZH Y, XU M, XU H M. Analysis on thermal damage effect of rubber materials by nuclear explosion[J]. Engineer Equipment Research, 2019, 38(3): 45-47. (in Chinese [15] SYMBALISTY E M D, ZINN J, WHITAKER R W. RADFLO physics and algorithms[R]. Washington: Los Alamos National Lab., 1995. [16] 陈健华, 王心正, 谢龙生, 等. 均匀大气中的强爆炸一维辐射流体力学数值解[J]. 爆炸与冲击,1981,1(2):37-49. CHEN J H, WANG X ZH, XIE L SH, et al. An one-dimensional radiation hydrodynamic numerical solution for a strong explosion in uniform atmosphere[J]. Explosion and Shock Waves, 1981, 1(2): 37-49. (in Chinese [17] 田宙, 乔登江, 郭永辉. 不同高度强爆炸早期火球数值研究[J]. 兵工学报,2009,30(8):1078-1083. TIAN ZH, QIAO D J, GUO Y H. Numerical investigation of early fireball of strong explosion for different altitudes[J]. Acta Armamentarii, 2009, 30(8): 1078-1083. (in Chinese [18] 高银军, 闫凯, 田宙, 等. 强爆炸早期火球光辐射能谱的数值计算[J]. 爆炸与冲击,2015,35(3):289-295. GAO Y J, YAN K, TIAN ZH, et al. Numerical calculation of early fireball radiation spectrum in strong explosion[J]. Explosion and Shock Waves, 2015, 35(3): 289-295. (in Chinese [19] SHESTAKOV A I, OFFNER S S R. A multigroup diffusion solver using pseudo transient continuation for a radiation-hydrodynamic code with patch-based AMR[J]. Journal of Computational Physics, 2008, 227(3): 2154-2186. doi: 10.1016/j.jcp.2007.09.019 [20] 吴健辉. 核爆炸光辐射特性及探测技术的理论与实验研究[D]. 武汉: 华中科技大学, 2009. WU J H. Study on Theory and experiment of the characteristics and detection technology of nuclear explosion radiation[D]. Wuhan: Huazhong University of Science & Technology, 2009. (in Chinese) [21] 高银军, 关奇, 陈志华, 等. 强爆炸热辐射在大气中的传输计算与分析[J]. 现代应用物理,2014,5(3):191-195. GAO Y J, GUAN Q, CHEN ZH H, et al. Calculation and analysis of atmospheric transmission of thermal radiation from strong explosion[J]. Modern Applied Physics, 2014, 5(3): 191-195. (in Chinese [22] LYU Y Y, GAO SH, XU Q J. Scattering characteristics of C@H2O composite particle based on Mie light scattering theory[J]. Chinese Journal of Luminescence, 2019, 40(3): 298-303. [23] 吴健辉, 杨坤涛, 张南洋生. 核爆炸光辐射探测中的大气传输性能研究[J]. 应用光学,2008,29(5):815-820. WU J H, YANG K T, ZHANG N Y SH. Atmosphere transmission in detection of nuclear explosion light radiation[J]. Journal of Applied Optics, 2008, 29(5): 815-820. (in Chinese [24] 王飞翔, 郭杰, 许方宇, 等. 不同海拔地区红外大气透过率的计算和测量[J]. 中国光学,2019,12(4):843-852. doi: 10.3788/co.20191204.0843 WANG F X, GUO J, XU F Y, et al. Calculation and measurement of infrared atmospheric transmittance at different altitudes[J]. Chinese Optics, 2019, 12(4): 843-852. (in Chinese doi: 10.3788/co.20191204.0843 [25] 刘志伟, 路远, 侯典心, 等. 激光辐照VO2薄膜温度场分布及透射特性研究[J]. 发光学报,2018,39(11):1604-1612. doi: 10.3788/fgxb20183911.1604 LIU ZH W, LU Y, HOU D X, et al. Temperature field distribution and transmittance characteristics of VO2 films irradiated by laser[J]. Chinese Journal of Luminescence, 2018, 39(11): 1604-1612. (in Chinese doi: 10.3788/fgxb20183911.1604 [26] 孙毅, 高云国, 邵帅. 高功率激光热效应对合束系统的影响[J]. 光学 精密工程,2015,23(11):3097-3106. doi: 10.3788/OPE.20152311.3097 SUN Y, GAO Y G, SHAO SH. Influence of high power laser thermal effect on beam combination system[J]. Optics and Precision Engineering, 2015, 23(11): 3097-3106. (in Chinese doi: 10.3788/OPE.20152311.3097 -