One-Dimensional Magneto-Hydrodynamics Calculation and Analysis of Implosion Magnetic Flux Compression Process
-
摘要: 在一维内爆磁流体力学编码MC11D的基础上,引入了解析形式的状态方程组。利用改造后的编码对样品为铜材料的柱面内爆磁通量压缩实验进行了模拟,得到了样品管内壁的速度历史,材料的p-V关系曲线以及温升曲线。与原编码利用列表式状态方程数据库计算的结果比较,由两种状态方程模拟得到的速度曲线符合程度很好,所得p-V曲线的最大相对偏差为6.6%,温升曲线基本吻合,且模拟的速度曲线与实验结果符合较好。通过对MC11D编码的改造,拓展了其应用范围,使其可对未包含于现有列表式状态方程数据库中的材料进行模拟,为实验结果分析及数据处理提供了一种方法。Abstract: In this work we modified the 1D implosion and the MHD code—MC11D by introducing the equation of state (EOS) in analytical form into the program.Then, using the modified code, we simulated a cylindrical implosion magnetic flux compression experiment on the copper sample and obtained the velocity histories of the sample tube's inner wall, the p-V lines of the copper sample, and the temperature histories of the material.The results from our simulation are closely consistent with those simulated with the original code using the tabular EOS database.In addition, the velocity curves obtained from the simulation are also found to accord with the experimental result.Our work expands the application scope of the MC11D code.Using the modified code, the materials which are not included in the tabular EOS database can also be simulated.It will serve as a useful method for the analysis of experimental results and data processing.
-
表 1 常见金属的相关状态方程参数
Table 1. Equation of state parameters for common metal
Material Q/(GPa) Q′ l M/(g/mol) β0/(10-7J·g-1·K-2) ρ0/(g/cm3) ρ0K/(g/cm3) Al 37.890 44 8.418 61 6 26.982 415.30 2.71 2.746 Fe 39.964 32 11.280 27 8 55.845 193.90 7.85 7.960 Cu 59.716 60 9.888 54 9 63.546 174.40 8.93 9.050 W 129.588 00 8.583 59 8 183.840 83.85 19.20 19.310 Pb 14.884 78 11.270 50 30 207.200 104.50 11.34 11.560 -
[1] SAKHAROV A D, LYUGAEV R Z, SMIRNOV E N, et al.Magnetic cumulation[J]. Uspekhi Fizicheskikh Nauk, 1991, 161(5):47-51. doi: 10.3367/UFNr.0161.199105f.0047 [2] ALTGILBERS L L, BROWM M D J, GRISHNAEV I, et al.磁通量压缩发生器[M].孙承纬, 周之奎, 译.北京: 国防工业出版社, 2008. [3] GU Z W, LUO H, ZHANG H D, et al.Cylindrical isentropic compression by ultrahigh magnetic field[J]. J Phys Confer Ser, 2014, 500(14):142018-142022. doi: 10.1088/1742-6596/500/14/142018 [4] ROBINSON A C, BRUNNER T A, CARROLL S, et al.ALEGRA: an arbitrary Lagrangian-Eulerian multimaterial, multiphysics code[C]//46th AIAA Aerospace Sciences Meeting and Exhibit.Reno, Nevada, 2008: 1235. https://www.researchgate.net/publication/269044194_ALEGRA_An_Arbitrary_Lagrangian-Eulerian_Multimaterial_Multiphysics_Code [5] RIEBEN R N, WHITE D A, WALLIN B K, et al.An arbitrary Lagrangian-Eulerian discretization of MHD on 3D unstructured grids[J]. J Comp Phys, 2007, 226(1):534-570. http://www.sciencedirect.com/science/article/pii/S0021999107001714 [6] OLIPHANT T A.Raven physics manual: LA-8802-M[R]. Los Alamos: Los Alamos National Laboratory, 1987. Raven physics manual: LA-8802-M [7] KOZLOV M B, ASEEVA V V, BORRISKOV G V, et al.Computational investigation of operating conditions of the cascade generator MC-1 with large HE charge[C]//Digest of Technical Papers, Pulsed Power Plasma Science 2001.Las Vegas, Nevada, 2001, 2: 1182-1184. https://www.researchgate.net/publication/3947101_Computational_investigation_of_operating_conditions_of_the_cascade_generator_MC-1_with_large_HE_charge [8] DOLOTENKO M I, ASEEVA V V, BORISKOV G V, et al.Numerical models of magnetocumulative generator MC-1 testing[C]//Digest of Technical Papers, Pulsed Power Plasma Science 2001.Las Vegas, Nevada, 2001, 2: 1185-1188. https://www.researchgate.net/publication/3947103_Numerical_models_of_magnetocumulative_generator_MC-1_testing [9] FRESE S D, FRESE M H.Recent improvements to MACH2 and MACH3 for fast Z-pinch modeling[J]. IEEE International Conference on Plasma Science, 2002, 651(1):380-383. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=CC026868665 [10] MACGILLIVRAY J, PETERKIN R E, BURKE D A, et al.MACH3: a CHSSI code for computational magnetohydrodynamics of general materials in generalized coordinates[C]//User Group Conference.Piscataway, New Jersy: IEEE Publishing, 2003: 228-235. https://www.researchgate.net/publication/4037342_MACH3_A_CHSSI_code_for_computational_magnetohydrodynamics_of_general_materials_in_generalized_coordinates [11] 王玉芝, 王贵荣, 董志伟, 等.二维磁流体力学程序MFCG Ⅱ与MFCG Ⅰ计算结果的对比分析[J].强激光与粒子束, 2003, 15(9):899-904. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200410020WANG Y Z, WANG G R, DONG Z W, et al.Contrast and analysis of the results computed by 2-D cylindrical coordinate system magnetohydrodynamics code MFCG Ⅱ and MFCG Ⅰ[J]. High Power Laser and Particle Beams, 2003, 15(9):899-904. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200410020 [12] 王玉芝, 王贵荣, 董志伟, 等.二维磁流体力学程序MFCG Ⅲ计算结果与测试结果的对比分析[J].强激光与粒子束, 2004, 16(10):1307-1312. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200410020WANG Y Z, WANG G R, DONG Z W, et al.Contrast and analysis of the results computed by 2-D magnetohydrodynamics code MFCG Ⅲ and experiment results[J]. High Power Laser and Particle Beams, 2004, 16(10):1307-1312. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200410020 [13] 丁宁, 邬吉明, 杨震华, 等.Z箍缩内爆的MARED程序1维模拟分析[J].强激光与粒子束, 2008, 20(2):212-218. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200802008DING N, WU J M, YANG Z H, et al.Simulation of Z-pinch implosion using MARED code[J]. High Power Laser and Particle Beams, 2008, 20(2):212-218. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200802008 [14] 孙奇志, 孙承纬, 刘伟.同轴型爆磁压缩发生器1维磁扩散模型[J].强激光与粒子束, 2007, 19(10):1755-1760. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200710039SUN Q Z, SUN C W, LIU W.1-D diffusion model of coaxial explosive magnetic generator[J]. High Power Laser and Particle Beams, 2007, 19(10):1755-1760. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200710039 [15] 张恒第.MC-1型爆磁压缩发生器一维内爆——磁流体力学模拟[D].绵阳: 中国工程物理研究院, 2012. http://cdmd.cnki.com.cn/Article/CDMD-82818-1012499394.htmZHANG H D.1D implosion & MHD simulation of MC-1 flux compression generator[D]. Mianyang: China Academy of Engineering Physics, 2012. http://cdmd.cnki.com.cn/Article/CDMD-82818-1012499394.htm [16] 赵继波, 孙承纬, 谷卓伟, 等.平面和圆柱面构型的磁流体力学计算[J].爆炸与冲击, 2016, 36(1):9-16. http://www.cnki.com.cn/Article/CJFDTotal-BZCJ201601002.htmZHAO J B, SUN C W, GU Z W, et al.Magnetic hydrodynamic calculation on planar and cylindrical configurations[J]. Explosion and Shock Waves, 2016, 36(1):9-16. http://www.cnki.com.cn/Article/CJFDTotal-BZCJ201601002.htm [17] BURGESS T J.Electrical resistivity model of metals: SAND 86-1093C[R]. Albuquerque, California: Sandia National Laboratory, 1986. https://www.researchgate.net/publication/252412772_Electrical_resistivity_model_of_metals [18] 徐锡申, 张万箱.实用物态方程理论引导[M].北京:科学出版社, 1986:443-444.XU X S, ZHANG W X.Introduction to the theory of equation of state[M]. Beijing:Science Press, 1986:443-444. [19] 汤文辉.固体零温零压密度的估算[J].国防科技大学学报, 1996(1):135-140. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199600107447TANG W H.Estimation of solid density under zero temperature and zero pressure[J]. Journal of National University of Defense Technology, 1996(1):135-140. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199600107447