Electromagnetic pulse driven liner implosion and compression of magnetized target
-
摘要: 采用磁流体方程和有限差分法，对内爆过程中高能量密度状态下的磁场对带电粒子和压缩过程的作用机制进行研究。结果显示：内爆过程中的各项参数为电子离子温度（50 keV）、压强（1 TPa）、粒子数密度（1024 cm-3）。套筒材料对约束时间、点火条件有重要影响；同时当磁感应强度大于5 T时，电子热传导系数比无磁场时减小2个数量级，离子热传导系数也出现了明显下降，在压缩峰值处，磁感应强度超过5 T时α粒子能量沉积密度比磁感应强度为0和1 T时相对增加约200倍。磁化在一定程度上也会阻碍内爆压缩过程。Abstract: In this paper, we studied the program implosion process and the influence of the magnetic field under a high-energy density state on charged particles and the mechanism of the compression process, using the MHD equations and finite difference method. The results show that the implosion process parameters meet the following conditions:the electron temperature plasma (50 keV), the pressure (1 TPa), the particle number density (1024 cm-3). The liner material had a major influence on the confined time and ignition conditions; at the same time, the electronic thermal conductivity was reduced by more than 100 times in the absence of a magnetic field when the magnetic field was greater than 5 T, and the ion thermal conductivity also experienced a significant decrease in the compression peak. The alpha particle energy deposition density in the 5 T embedded magnetic field was relatively raised by more than about 200 times than in the 1 T embedded magnetic field. Also, it was found that magnetization hindered implosion compression process to a certain extent.
-
Key words:
- electromagnetically driven implosion /
- magnetized target /
- Z-pinch /
- MHD simulation
-
飞机座舱风挡、座舱盖和弦舱材料为高聚物材料, 在飞机飞行时需要承受一定的温度和载荷作用。在起飞和着陆时, 还有可能会遇到飞鸟的撞击, 军用飞机可能在空战中被非爆弹丸或高爆碎片击中, 这相当于对这些位置的材料施加一定的高速率载荷。同时, 由于飞机在低、高空飞行时, 环境温度相差很大, 高聚物材料的力学性能受外界温度的影响很大。因此, 开展对高聚物材料在高温和高应变率下的力学性能研究具有重要意义。
C.G’Sell等[1]对聚合物在常应变率下的变形行为进行了研究, 发现应变率影响材料的应变硬化现象。陈建桥等[2]、吴衡毅等[3]通过恒应变率的实验对有机玻璃进行了研究, 结果表明, 在较低应变率下, 有机玻璃的力学性能具有明显的应变率相关性。E.M.Arruda等[4]、O.S.Lee等[5]分别研究了准静态及中、高应变率下温度和应变率对PMMA材料力学性能的影响, 发现随着温度的升高, 材料的应力应变行为发生了显著变化, 同时伴随有热软化现象发生。G.W.Adams等[6]、Machida等[7]通过研究温度和应变率对PC材料力学性能的影响后发现:温度在293~338K之间时, 材料变形过程中储存的能量对实验温度的依赖性并不明显。S.M.Walley等[8]在一定的应变率范围内对一些高聚物材料进行了力学性能测试, 解释了在高应变率下由于绝热效应的原因, 使得材料出现应变软化行为。索涛等[9-10]对几种航空有机玻璃进行了研究, 认为在某一应变率和某一温度下, 由于试样内部的微裂纹使得有机玻璃表现出明显的应变软化行为; 在相同温度下, 随着应变率的升高, 有机玻璃的流动应力增大, 而破坏应变减小, 并建立了考虑温度效应的本构模型。
本文中对YB-2航空有机玻璃在不同温度下的准静态和动态力学行为进行系统的研究。分别采用Hopkinson分离式压杆和CSS44100型电子万能试验机对材料进行初始温度分别为218、253、293、343、373 K, 应变率分别为0.001、0.1、700、3 000 s-1下的压缩实验。同时, 结合实验数据, 基于朱-王-唐模型, 拟合考虑应变率效应以及温度效应的本构模型, 并给出模型预测和实验所得的应力应变曲线。
1. 实验过程与方法
1.1 导热时间的测定实验
有机玻璃是一种高分子材料, 对温度等环境比较敏感。由于聚合体材料的导热系数很低, 约是金属的1%, 即聚合体的导热系数为0.2 mm/min, 试样在加热或冷却的过程中, 如果加温或冷却的时间不充分, 试样表面温度和心部温度就会不一致, 导致实验结果不准确。因此, 在涉及温度的实验时, 为了确保实验中试样受热均匀和热平衡, 从而保证实验数据的准确性, 首先测试了YB-2有机玻璃在不同温度下试样表面和试样内部达到相同温度时所需要的时间。实验中, 采用正方体试样, 边长为10 mm。在试样中心沿着厚度方向钻孔, 孔的深度为5 mm。通过孔洞, 在试样心部引入一对热电偶。将另外一对热电偶缠绕在试样表面。用密闭的电热炉, 对试样进行加热或在密闭的容器中喷洒低温氮气对试样进行冷却。当经过标定的两组热电偶显示的温度相同时, 认为试样内部和外部的温度达到了平衡。
1.2 低应变率实验
低应变率压缩实验在CSS44100型电子万能试验机上进行, 压缩试样为直径6 mm、长4 mm的圆柱体, 如图 1所示。两种低应变率分别为0.1、0.001 s-1, 温度分别为218、253、293、343、373 K。实验前将压头降温或加热至相对应的实验温度后对试验机进行空压, 得到位移随着压力变化关系, 以此标定试验机刚度, 用于消除试样位移测量中试验机位移的影响。低温通过液氮循环获得, 高温实验采用电阻丝加热法对试样进行加热, 并用一个与试样相连的热电偶测量和显示温度, 通过电磁继电器的通断控制温度, 使之保持在实验值的±5 ℃误差范围内。
1.3 高应变率实验
700、3 000 s-1时的高应变率压缩实验是在一套直径12.7 mm的分离式Hopkinson压杆上进行的, 压杆所用材料为马氏体时效钢(18Ni(C250)), 实验装置见图 2。采取同低应变率下相同的试样。利用撞击杆撞击入射杆产生的应力波对放置于入射杆和透射杆之间的试样进行加载; 通过粘贴在两根杆子上与惠斯通电桥相连接的应变片采集输出的电压信号; 利用超动态应变仪和高频数据采集系统记录到实测波形, 并根据一维应力波原理可以得到材料的应力应变关系。高低应变下低温和高温的获得采用相同的方法。
2. 实验结果及分析
2.1 导热时间的测定
图 3(a)和(b)为加热或者冷却温度与保温时间的曲线以及不同的实验温度所需最短时间值。可以看出图 3(a)中, 当经过一定时间后, 无论是高温加热还是低温冷却, 试样内部和外部的温度达到一致, 记录加热或冷却的时间, 就得到了该温度下试样内外温度达到平衡所需的时间。从图 3(b)可以看出, 低温冷却试样内外温度平衡所需时间比高温时所需时间要短, 随着温度的逐渐升高, 平衡所需时间呈现逐渐变短的趋势。当温度高于有机玻璃的玻璃化转变温度后, 平衡时间变长。但如果温度继续升高, 平衡所需的时间几乎相同。
图 3 不同温度下保温时间曲线和保温所需最小时间Figure 3. Time-temperature curves and minimum times at different temperatures2.2 低应变率和高应变率压缩实验
YB-2航空有机玻璃有别于定向有机玻璃, 是未经过取向的高聚物, 分子链段的排列是随机的, 因此理论上材料体现出各向同性。在开始实验之前, 为了测试该材料的各向性能, 首先对其进行了考虑方向性的准静态压缩实验。试样分别取自于垂直于板材和平行于板材两个方向, 见图 1。实验结果(见图 4)表明:沿2个方向取样所测出的应力应变曲线基本吻合, 这说明材料力学各向同性。在后续的实验中采用的试样沿V方向取样。
图 5给出了在不同温度下试样的真实应力应变曲线。可以看出, 高分子材料的力学性能强烈依赖于温度和应变率。图 5(a)中, 当应变率为0.001 s-1时, 随着温度的升高, 材料的弹性模量呈现降低的趋势, 屈服强度减小, 塑性变形能力增加。在其余的3种应变率下(见图 5(b)~(d)), 材料呈现相似的趋势。在室温及以下的高应变率时, 材料的破坏应变随着温度的升高而增大, 此时, 材料为玻璃体呈现出典型的脆性破坏形式, 这是因为在低温高应变率下, 材料处于玻璃态, 且加载时间非常短暂, 材料的分子结构没有足够的时间发生弛豫, 因而断裂特性表现为脆性; 而当实验温度在材料的玻璃化转变温度(90 ℃)附近时, 材料处于高弹态, 其承载能力急剧下降, 变形能力增强, 不再呈现脆性破坏, 这是由于分子结构的活化能降低、自由体积增大, 因此受力后分子结构中的链段运动更加容易[11]。
图 6给出了同一温度下不同应变率时的应力应变曲线。在218 K时(见图 6(a)), 应变率对材料的弹性模量和屈服应力的影响并不明显, 但屈服应力仍呈现随应变率升高而升高的趋势, 破坏应变则反之。当温度升高时(见图 6(b)~(d)), 应变率对材料的弹性模量的影响逐渐变大, 屈服应力以及破坏应变呈现同218K时相似的趋势。4种应变率下, 材料均出现不同程度的应变软化。低应变率时, 由于加载脉冲宽度有限, 材料在实验中并未发生破坏。已有研究结果[12]表明:低温时, 聚合物处于玻璃态, 分子链的链段被冻结发生松弛的时间是无穷大, 此时, 加载速率(即应变率)的改变对聚合物材料变形的影响作用非常微弱, 表现为弹性模量对应变率的不敏感; 而在高温时, 聚合物处于高弹态, 其变形实际上是一种强迫高弹形变的过程, 分子链的链段已经被激活, 并且分子链的链段运动受到松弛时间的影响非常大。由于应变率的变化对松弛时间的影响很大, 因此在高温时弹性模量对应变率较为敏感。
图 7为不同应变率时温度-应力峰值曲线。当实验温度为218 K时, 不同应变率下试样的应力峰值差值较小。随着温度的升高, 应力峰值的差值变大, 在室温时应力峰值的差值达到最大。随后, 该差值又逐渐变小。这表明:在室温时材料的应变率敏感性达到最强。在高应变率加载下, 应变率对材料应力峰值的影响逐渐弱化。
3. YB-2有机玻璃的本构模型
朱兆祥等[13]对典型工程塑料进行的一系列研究发现:在变形小于8%, 应变率在10-4~103 s-1范围内, 朱-王-唐模型(ZWT)可以较好地描述高聚物材料的本构模型。该模型描述如下:
$$ \sigma(\varepsilon, \dot{\varepsilon})=E_{0} \varepsilon+\alpha \varepsilon^{2}+\beta \varepsilon^{2}+E_{1} \int_{0}^{t} \dot{\varepsilon} \exp \left(-\frac{t-\tau}{\theta_{1}}\right) \mathrm{d} t+E_{2} \int_{0}^{t} \dot{\varepsilon} \exp \left(-\frac{t-\tau}{\theta_{2}}\right) \mathrm{d} t $$ (1) 式中:E0、α和β为材料的弹性常数。2个积分项分别用来描述低应变率和高应变率下的黏弹性响应, E1、θ1和E2、θ2分别是所对应的Maxwell单元的弹性常数和松弛时间。
为描述温度对材料的影响, 索涛等[10]对朱-王-唐模型进行修正, 并分别对4种航空有机玻璃实验所得的真实应力应变曲线进行拟合, 实验与模型预测的曲线吻合较好。修正的模型中加入的温度效应项, 其形式为
/f0(ε, T0), 它随温度的升高而减小, 随应变的增大而增大。改进后的模型描述为:
$$ \sigma(\varepsilon, \dot{\varepsilon}, T)=\sigma(\varepsilon, \dot{\varepsilon}) f(\varepsilon, T) / f\left(\varepsilon, T_{0}\right) $$ (2) $$ f(\varepsilon, T)=a(\varepsilon) \exp [b(\varepsilon) T] $$ (3) 式中:
)为式(1)对应的温度和应变率下得到的应力值, T0为选定的参考温度。
根据公式(1)~(3)可得到温度影响项的表达式:
, 由实验数据拟合以多项式
给出。
基于该模型, 以293 K作为参考温度, 针对YB-2试样, 经过拟合实验曲线得到常温时的6个参数:E0=3.26 GPa、α=-20.12 GPa、E1=2.4 GPa、θ1=0.205 s、E2=2.622 GPa、θ2=6 μs。
在此基础上, 采用最小二乘法拟合可以得到YB-2航空有机玻璃温度效应影响项中b(ε)的估计表达式为:
$$ b(\varepsilon ) = 0.6429{\varepsilon ^2} - 0.1336\varepsilon - 0.004$$ (4) 对于拟合参数优度的评估, 我们将对回归参数的方差和置信区间进行分析。根据实验结果, 对公式(4)中各系数的方差和置信区间的分析如下:假定置信系数为0.95, 计算得出拟合函数的方差和置信区间见表 1。
表 1 温度影响项中各系数的方差和置信区间Table 1. Variances and confidence intervals of the coefficients for thermo-item参数 方差 置信区间 a 0.381 131 (0.527 425, 1.273 1) b 0.009 826 (-0.235 72, -0.103 88) c 0.000 295 (-0.004 56, -0.003 54) 图 8给出了拟合结果与实验结果的比较。不难看出, 图 8(a)~(c)中拟合曲线同实验曲线吻合较好; 对于图 8(d), 拟合曲线和实验曲线有较大差异, 这是因为在PMMA的热变形温度(70 ℃)左右时, 材料将开始发生玻璃态与高弹态之间的相态转变。而ZWT模型主要用来描述材料的黏弹性行为, 因此在拟合343 K下材料的变形行为时表现出较大的误差。总体来说, 改进后的ZWT模型可以较好的描述一定温度范围内材料的非线性黏弹性力学行为。
图 8 修正的ZWT模型拟合结果与实验结果的比较Figure 8. Comparison of prediction using modified ZWT model with experimental results4. 结论
通过对YB-2航空有机玻璃的研究得出以下结论:
(1) YB-2有机玻璃的力学性能强烈依赖于应变率和温度:在低应变率下, 随着温度的升高, 材料的弹性模量和流动应力逐渐减小, 弹性变形能力降低, 塑性变形能力增加; 在高应变率下, 材料的流动应力随着温度的升高逐渐减小, 破坏应变随着温度的升高而增大。在室温及以下时, 材料为玻璃体, 呈现出典型的脆性破坏, 弹性模量对应变率不敏感; 当实验温度在材料的玻璃化转变温度附近时, 高应变率下材料破坏不再呈现脆性破坏, 弹性模量对应变率较为敏感。
(2) 随着应变率的升高, 材料的流动应力显著升高, 但破坏应变减小。在室温时, 材料的应变率敏感性达到最强。在高应变率下, 材料的屈服强度受应变率的影响逐渐弱化。
(3) 加入温度效应的ZWT模型能够较好地拟合应变范围在8%以内的实验曲线, 但由于模型中未考虑损伤以及高分子高温下自由体积和复杂的分子链运动的影响, 因而不能较准确地描述应变软化以及大应变时的应变硬化现象。
-
图 6 磁场对电子热传导系数的影响
Figure 6. Effect of magnetic field on the electron heat conduction coefficient
图 7 磁场对离子热传导系数的影响
Figure 7. Effect of magnetic field on the ion heat conduction coefficient
图 8 磁场对α粒子能量沉积的影响
Figure 8. Effect of magnetic field on the alpha particle energy deposition
图 9 等离子体加热能量随磁场的变化
Figure 9. Plasma heating energy varying with the change of magnetic field
-
[1] TACCETTI J M, INTRATOR T P, WURDEN G A, et al. FRX-L:A field-reversed configuration plasma injector for magnetized target fusion[J]. Review of Scientific Instruments, 2003, 74(10):4314-4323. doi: 10.1063/1.1606534 [2] DEGNAN J H, AMDAHL D J, Domonkos M, et al. Recent Magneto-inertial fusion experiments on FRCHX[C]//24nd IAEA Fusion Energy Conference. San Diego: 2012. [3] WUEDEN G A, GRABOWSKI T C, DEGNAN J H, et al. Increased FRC lifetimes using a longer trap[C]//Bulletin of the American Physical Society (APS Meeting), 2013: 58. [4] SLUTZ S A, HERRMANN M C, VESEY R A, et al. Pulsed-power-driven cylindrical liner implosions of laser preheated fuel magnetized with an axial field[J]. Physics of Plasmas, 2010, 17(5):263-52. [5] CUNEO M E, HERRMANN M C, SINARS D B, et al. Magnetically driven implosions for inertial confinement fusion at sandia national laboratories[J]. IEEE Transactions on Plasma Science, 2012, 40(12):3222-3245. doi: 10.1109/TPS.2012.2223488 [6] GOTCHEV O V, KNAUER J P, CHANG P Y, et al. Seeding magnetic fields for laser-driven flux compression in high-energy-density plasmas[J]. Review of Scientific Instruments. 2009, 80(4):495. http://cat.inist.fr/?aModele=afficheN&cpsidt=21489875 [7] HOHENBERGER M, CHANG P Y, FIKSEL G, et al. Inertial confinement fusion implosions with imposed magnetic field compression using the OMEGA Lasera[J]. Physics of Plasmas, 2012, 19(5):139. https://www1.psfc.mit.edu/research/hedp/Home%20Page/Papers/Hohenberger_PoP-2012.pdf [8] LABERGE M. Experimental results for an acoustic driver for MTF[J]. Journal of Fusion Energy, 2009, 28(2):179-182. doi: 10.1007/s10894-008-9181-y [9] HSU S C, WITHERSPOON F D, CASSIBRY J T, et al. Overview of the plasma liner experiment (PLX)[C]//Bulletin of the American Physical Society (APS Meeting), 2009: 56. [10] GARANIN S F, MAMYSHEC V I, YAKUBOV V B. The MAGO system:Current status[J]. Plasma Science, IEEE Transactions, 2006, 34(5):2273-2278. doi: 10.1109/TPS.2006.878368 [11] 孙奇志, 方东凡, 刘伟, 等."荧光-1"实验装置物理设计[J].物理学报, 2013, 62(7):000507. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=wlxb201307074SUN Qizhi, FANG Dongfan, LIU Wei, et al. Physical design of the "Ying-Guang 1" device[J]. Acta Physica Sinica, 2013, 62(7):000507. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=wlxb201307074 [12] 李璐璐, 张华, 杨显俊.反场构型的二维磁流体力学描述[J].物理学报, 2014, 63(16):165202. doi: 10.7498/aps.63.165202LI Lulu, ZHANG Hua, YANG Xianjun.Two-dimensional magneto-hydrodynamic description of field reversed configuration[J]. Acta Physica Sinica, 2014, 63(16):165202. doi: 10.7498/aps.63.165202 [13] GIBBS W W. Triple-threat method sparks hope for fusion[J]. Natrue, 2014, 505(7481):9-10. https://www.scientificamerican.com/article/triple-threat-method-sparks-hope-for-nuclear-fusion-energy/ [14] GOMEZ M R, SLUTZ S A, SEFKOW A B, et al. Experimental verification of the magnetized liner inertial fusion (MagLIF) concept[C]//IEEE International Conference on Plasma Sciences, 2014: 1. [15] 邓建军, 王勐, 谢卫平, 等.面向Z箍缩驱动聚变能源需求的超高功率重复频率驱动器技术[J].强激光与粒子束, 2014, 26(10):100201. http://industry.wanfangdata.com.cn/dl/Detail/Periodical?id=Periodical_qjgylzs201410001DENG Jianjun, WANG Meng, XIE Weiping, et al. Super-power repetitive Z-pinch driver for fusion-fission reactor[J].High Power Laser and Particle Beams, 2014, 26(10):100201. http://industry.wanfangdata.com.cn/dl/Detail/Periodical?id=Periodical_qjgylzs201410001 [16] BASKO M. A1-D 3-T hydrodynamic code for simulating ICF targets driven by fast ion beams[R]. Version 4. Institute for Theoretical and Experimental Physics, Moscow, 2001. [17] 邓爱东, 张华, 杨显俊.电磁驱动产生超强磁场的参数优化设计[J].高压物理学报, 2015, 29(2):123-128. doi: 10.11858/gywlxb.2015.02.006DENG Aidong, ZHANG Hua, YANG Xianjun. Parameters optimization of the strong magnetic gield generation driven by electromagnetic force[J]. Chinese Journal of High Pressure Physics, 2015, 29(2):123-128. doi: 10.11858/gywlxb.2015.02.006 [18] 宁成, 丰志兴, 薛创.Z箍缩驱动动态黑腔中的基本能量转移特征[J].物理学报, 2014, 63(12):125208. doi: 10.7498/aps.63.125208NING Cheng, FENG Zhixing, XUE Chuang. Basic characteristics of kinetic energy transfer in the dynamic hohlraums of Z-pinch[J]. Acta Physica Sinica, 2014, 63(12):125208. doi: 10.7498/aps.63.125208 [19] ANDREAS J K. Magnetized cylindrical implosions driven by heavy ion beams[R]. Max Planck Institute of Quantum Optics, 2001. [20] 谷同祥, 安恒斌, 刘兴平, 等.迭代法和预处理技术:上册[M].北京:科学出版社, 2015:78-90. [21] MEYER-TER-VEHN J. 惯性聚变物理[M]. 沈柏飞, 译. 北京: 科学出版社, 2008: 28-31. [22] BASKO M M, KEMP A J, MEYER-TER-VEHN J. Ignition conditions for magnetized target fusion in cylindrical geometry[J]. Nuclear Fusion, 2000, 40(1):41-45. [23] 刘斌, 李成, 邓爱东, 等.脉冲驱动磁化等离子体内爆升温点火的数值模拟[J].强激光与粒子束, 2016, 28(7):075010. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=qjgylzs201607011LIU Bin, LI Cheng, DENG Aidong, et al. The numerical modeling about ignition and implosion heating process of magnetized plasma driven by pulsed-power[J]. High Power Laser and Particle Beams, 2016, 28(7):075010. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=qjgylzs201607011 期刊类型引用(1)
1. 张世文，龙建华，贾宏志，刘仓理. 柱壳约束对散心冲击波空间分布的影响. 爆炸与冲击. 2018(02): 345-352 . 
本站查看其他类型引用(3)
-
下载:
下载:
计量
- 文章访问数: 5353
- HTML全文浏览量: 1926
- PDF下载量: 169
- 被引次数: 4