A Numerical Simulation Analysis of Mono-Temperature and Tri-Temperature Models by MDSC Program in Z-Pinch Implosion
-
摘要: 用于模拟Z箍缩物理过程的辐射磁流体力学程序(MDSC)存在单温和三温两种模型, 了解其差别有助于更好地使用该工具。介绍了程序中单温模型和三温模型在物理建模上的差异, 设计了3个算例, 采用两种模型分别予以模拟计算。结果表明:相对于单温模型, 三温模型获得更高的碰轴速度, 该差异在低温大质量负载上表现得更为明显。结合Z箍缩的实际物理过程, 分析认为:碰轴速度的差异主要由状态方程引起, 低温下单温模型中的原子状态方程与三温模型中的离子状态方程有显著差异, 导致加载面的膨胀不一致, 最终使获得的动能不一致, 而高温下两种状态方程的差异减小, 因此碰轴速度的差异不明显。状态方程的差异还导致三温模型中更容易形成冲击波和先驱等离子体, 从而缩短了碰轴时间。动能的不同是引起计算差异的主要因素。Abstract: Given the mono-temperature model and tri-temperature model as the two radiation-magneto-hydrodynamic models in the physical process of the Z-pinch, it is useful to understand the difference between the two models for a better use of these tools.In this paper the difference between the two models in the process of modeling is explained and three calculating models are designed to study the calculation difference resulting from the two models.The simulation results achieved show that the tri-temperature model has a faster spindle-touching speed, which is more obvious in low-temperature and heavy loads.The research by taking into account the physical reality shows that the difference between the two models is mainly caused by the equation of state (EOS):at low temperatures, there is a significant difference between atomic EOS and ionic EOS, hence the loaded surface has an uneven inflation, which eventually leads to the difference of the received kinetic energy; at high temperatures, this difference is reduced.Moreover, the EOS makes it more apt for impulse waves to be formed in the tri-temperature program than in the mono-temperature program and for the pioneering plasma to be formed as well.However, the calculation difference caused by the kinetic energy plays a major role.
-
表 1 3个计算模型的初始参数
Table 1. Initial parameters for three calculation models
Model No. Material mtot/(mg) H/(mm) D/(mm) δ/(mm) ρ0/(g/cm3) T0/(K) 1 W 1.746 20 21.0 1.0 1.389 10 000 2[16] Al 0.656 20 18.2 1.0 0.607 10 000 3 Al 85.536 12 5.0 0.6 2.700 300 表 2 3个计算模型的MDSC程序模拟结果
Table 2. MDSC calculation results of the three models
Model No. Material timp/(ns) vmax/(km/s) Exp. 1T model 3T model 1T model 3T model 1 W 111 114 111 426 612 2 Al 84[16] 83 82 593 713 3 Al 448 418 12 18 -
[1] PETERKIN R E Jr, DEGNAN J H, HUSSEY T W, et al.A long conduction time compact torus plasma opening switch[J].IEEE Trans Plasma Sci, 1993, 21(5):522-528. doi: 10.1109/27.249638 [2] HELLINGER P, MANGENEY A, MATTHEWS A.Whistler waves in 3D hybrid simulations of quasiperpendicular shocks[J].Geophys Res Lett, 1996, 23(6):621-624. doi: 10.1029/96GL00453 [3] CHITTENDEN J P, LEBEDEV S V, JENNINGS C A, et al.X-ray generation mechanisms in three-dimensional simulations of wire array Z-pinches[J].Plasma Phys Control Fusion, 2004, 46(12B):457-476. doi: 10.1088/0741-3335/46/12B/039 [4] 阚明先, 胡熙静, 王刚华.二维磁流体力学ALE数值模拟[J].高能量密度物理, 2009(2):56-60. http://www.cqvip.com/Main/Detail.aspx?id=30778709KAN M X, HU X J, WANG G H.Two dimensional MHD ALE numerical simulation[J].High Energy Density Physics, 2009(2):56-60. http://www.cqvip.com/Main/Detail.aspx?id=30778709 [5] 丁宁, 邬吉明, 戴自换, 等.Z箍缩内爆的MARED程序数值模拟分析[J].物理学报, 2010, 59(12):8707-8716. doi: 10.7498/aps.59.8707DING N, WU J M, DAI Z H, et al.Numerical simulation analysis of Z-pinch implosion using MARED code[J].Acta Physica Sinica, 2010, 59(12):8707-8716. doi: 10.7498/aps.59.8707 [6] ROBINSON A C, GARASI C J.Three-dimensional z-pinch wire array modeling with ALEGRA-HEDP[J].Comput Phys Commun, 2004, 164(1/2/3):408-413. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=9452d46b8e23e77426d775907430be86 [7] 阚明先, 蒋吉昊, 王刚华, 等.衬套内爆ALE方法二维MHD数值模拟[J].四川大学学报(自然科学版), 2007, 44(1):91-96. doi: 10.3969/j.issn.0490-6756.2007.01.020KAN M X, JIANG J H, WANG G H, et al.ALE simulation of 2D MHD for liner[J].Journal of Sichuan University (Natural Science Edition), 2007, 44(1):91-96. doi: 10.3969/j.issn.0490-6756.2007.01.020 [8] MARTIN M R, LEMKE R W, MCBRIDE R D, et al.Solid liner implosions on Z for producing multi-megabar, shockless compressions[J].Phys Plasmas, 2012, 19(5):056310. doi: 10.1063/1.3694519 [9] MCBRIDE R D, MARTIN M R, LEMKE R W, et al.Beryllium liner implosion experiments on the Z accelerator in preparation for magnetized liner inertial fusion[J].Phys Plasmas, 2013, 20(5):056309. doi: 10.1063/1.4803079 [10] 阚明先, 王刚华, 赵海龙, 等.金属电阻率模型[J].爆炸与冲击, 2013, 33(3):282-286. doi: 10.3969/j.issn.1001-1455.2013.03.010KAN M X, WANG G H, ZHAO H L, et al.Electrical resistivity model for metals[J].Explosion and Shock Waves, 2013, 33(3):282-286. doi: 10.3969/j.issn.1001-1455.2013.03.010 [11] 王刚华, 胡熙静, 孙承纬.双层钨衬套Z箍缩内爆数值模拟[J].高压物理学报, 2004, 18(4):364-367. doi: 10.3969/j.issn.1000-5773.2004.04.013WANG G H, HU X J, SUN C W.Simulation of magnetohydrodynamics for plasma jetting on wire pinch[J].Chinese Journal of High Pressure Physics, 2004, 18(4):364-367. doi: 10.3969/j.issn.1000-5773.2004.04.013 [12] 阚明先, 王刚华, 赵海龙, 等.磁驱动飞片二维磁流体力学数值模拟[J].强激光与粒子束, 2013, 25(8):2137-2140. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs201308052KAN M X, WANG G H, ZHAO H L, et al.Two-dimensional magneto-hydrodynamic simulations of magnetically accelerated flyer plates[J].High Power Laser and Particle Beams, 2013, 25(8):2137-2140. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs201308052 [13] 王刚华, 孙承纬, 赵剑衡, 等.磁驱动平面飞片的一维磁流体力学计算[J].爆炸与冲击, 2008, 28(3):261-264. doi: 10.3321/j.issn:1001-1455.2008.03.011WANG G H, SUN C W, ZHAO J H, et al.One-dimensional, magnetohydrodynamic simulations of magnetically driven flyer plates[J].Explosion and Shock Waves, 2008, 28(3):261-264. doi: 10.3321/j.issn:1001-1455.2008.03.011 [14] HAINS M G.A heuristic model of the wire array Z-pinch[J].IEEE Trans Plasma Sci, 1998, 26(4):1275-1281. doi: 10.1109/27.725160 [15] 徐家鸾, 金尚宪.等离子体物理学[M].北京:原子能出版社, 1981:37-40.XU J L, JIN S X.Plasma physics[M].Beijing:Atomic Energy Press, 1981:37-40. [16] SANDFORD T W L, ALLSHOUSE G O, MARDER B M, et al.Improved symmetry greatly increases X-ray power from wire-array Z-pinches[J].Phys Rev Lett, 1996, 77(25):5063-5066. doi: 10.1103/PhysRevLett.77.5063 [17] COOPER N G.An invitation to participate in the LASL equation of state library: LASL-79-62[R].Los Alamos, USA: Los Alamos Scientific Laboratory, 1979. [18] LYON S P, JOHNSON J D.T-1 handbook of the SESAME equation of state library: LA-CP-98-100[R].Los Alamos, USA: Los Alamos Scientific Laboratory, 1998.