平面和圆柱面构形的磁流体力学计算
|
摘要
在一维弹塑性反应流体力学Lagrange编码SSS的基础上,研制了多介质、多空腔、与集中参数外电路耦合的实验构形弹塑性反应磁流体力学一维计算编码SSS-MHD。在编码中,通过对构形边界与外电路的耦合以及不同构形磁场边界条件的统一处理等措施,解决了其中的关键问题。本文中编码的优势体现在两个方面:能够统一处理不同的磁流体力学实验构形,包括平面及圆柱面几何构形;可以同时进行磁流体力学和含能材料反应流动的计算。对平面准等熵压缩炸药和炸药爆轰驱动套筒压缩磁通量等两类典型实验进行了数值计算,结果与实验结果一致。
In this work we developed a one-dimensional elastic-plastic reactive magnetic hydrodynamic code SSS/MHD with multi-medium,multi-cavities and coupling with external circuits on experimental configurations,based on the Lagrange code SSS.The key problems in coding were resolved through defining the configuration boundary coupling with circuits and unifying the treatment of the magnetic field boundary conditions of various configurations.The advantages of our coding are mainly shown in that,on the one hand,the various magnetic hydrodynamic experiment configurations-including both planar and cylindrical ones-can be handled in a unified fashion and,on the other,that the magnetic and reactive hydrodynamic calculation can be simultaneously conducted.Thus,following our way of coding,two typical experiments were calculated,i.e.the magnetic driven isentropic compression planar explosive and the explosive detonation driven liner for flux compression,the results of which accord with the experimental results.
In this work we developed a one-dimensional elastic-plastic reactive magnetic hydrodynamic code SSS/MHD with multi-medium,multi-cavities and coupling with external circuits on experimental configurations,based on the Lagrange code SSS.The key problems in coding were resolved through defining the configuration boundary coupling with circuits and unifying the treatment of the magnetic field boundary conditions of various configurations.The advantages of our coding are mainly shown in that,on the one hand,the various magnetic hydrodynamic experiment configurations-including both planar and cylindrical ones-can be handled in a unified fashion and,on the other,that the magnetic and reactive hydrodynamic calculation can be simultaneously conducted.Thus,following our way of coding,two typical experiments were calculated,i.e.the magnetic driven isentropic compression planar explosive and the explosive detonation driven liner for flux compression,the results of which accord with the experimental results.
引文
[1]Robinson A C,Brunner T A,Carrol S,et al.ALEGRA:An arbitrary Lagrangian-Eulerian multimaterial,multiphysics code[R].AIAA 2008-1235,2008.
[2]Lemke R W,Knudson M D,Robinson A C,et al.Self-consistent,two-dimensional,magnetohydrodynamic simulations of magnetically driven flyer plates[J].Physics of Plasmas,2003,10(5):1867-1874.
[3]Knudson M D,Lemke R W,Hayes D B,et al.Near-absolute Hugoniot measurements in aluminum to 500GPa using a magnetically accerlerated flyer plate technique[J].Journal of Applied Physics,2003,94(7):4420-4431.
[4]Knudson M D,Hanson D L,Bailey J E,et al.Principal Hugoniot,reverberating wave,and mechanical reshock measurement of liquid deuterium to 400GPa using plate impact techniques[J].Physical Review B,2004,69(14):144209.
[5]Oliphant T A,Witte K H.RAVEN[R].NM:Los Alamos National Laboratory,1987.
[6]Sheppard M G,Fowler G V,Freeman B L.Generation of ultra-high magnetic fields for AGEX[R].NM:Los Alamos National Laboratory,1994.
[7]Kozlov M B,Aseeva V V,Boriskov G V,et al.Computational investigational of operating conditions of the cascade generator MC-1with large HE charge[C]∥Proceedings of 28th IEEE International Conference on Plasma Science and 13th IEEE International Pulsed Power Conference.Las Vegas,NV,2001:1182-1184.
[8]Frese M H.MACH2:A two-dimensional magnetohydrodynamic simulation code for complex experimental configurations[R].AMRC-R-874,1987.
[9]MacGillivray J,Peterkin R E,Burke D A,et al.MACH3:A CHSSI code for computational magnetohydrodynamics of general materials in generalized coordinates[C]∥IEEE Computer Society.Proceedings of the users group2003conference.Bellevue,WA,2003:9-13.
[10]Frese S D,Frese M H.Recent improvements to MACH2and MACH3for fast Z-pinch modeling[C]∥Proceedings of 5th International Conference on dense Z-pinch.Albuquerque,NM,2002:380-383.
[11]L’Eplattenier P,Cook G,Ashcraft C,et al.Introduction of an electromagnetism module in LS-DYNA for coupled mechanical-thermal-electromagnetic simulations[J].Steel Research International,2009,80(5):351-358.
[12]王刚华.磁驱动准等熵压缩和高速飞片实验、计算和反积分数据处理技术[D].绵阳:中国工程物理研究院,2008.
[13]张恒第.MC-1型爆磁压缩发生器一维内爆:磁流体力学模拟[D].绵阳:中国工程物理研究院,2012.
[14]张旭平.磁驱动超高速飞片发射技术研究[D].绵阳:中国工程物理研究院,2012.
[15]孙承纬.一维冲击波和爆轰波计算程序SSS[J].计算物理,1986,3(2):142-154.Sun Chengwei.SSS:A code computing one dimensional shock and detonation wave propagation[J].Chinese Journal of Computational Physics,1986,3(2):143-154.
[2]Lemke R W,Knudson M D,Robinson A C,et al.Self-consistent,two-dimensional,magnetohydrodynamic simulations of magnetically driven flyer plates[J].Physics of Plasmas,2003,10(5):1867-1874.
[3]Knudson M D,Lemke R W,Hayes D B,et al.Near-absolute Hugoniot measurements in aluminum to 500GPa using a magnetically accerlerated flyer plate technique[J].Journal of Applied Physics,2003,94(7):4420-4431.
[4]Knudson M D,Hanson D L,Bailey J E,et al.Principal Hugoniot,reverberating wave,and mechanical reshock measurement of liquid deuterium to 400GPa using plate impact techniques[J].Physical Review B,2004,69(14):144209.
[5]Oliphant T A,Witte K H.RAVEN[R].NM:Los Alamos National Laboratory,1987.
[6]Sheppard M G,Fowler G V,Freeman B L.Generation of ultra-high magnetic fields for AGEX[R].NM:Los Alamos National Laboratory,1994.
[7]Kozlov M B,Aseeva V V,Boriskov G V,et al.Computational investigational of operating conditions of the cascade generator MC-1with large HE charge[C]∥Proceedings of 28th IEEE International Conference on Plasma Science and 13th IEEE International Pulsed Power Conference.Las Vegas,NV,2001:1182-1184.
[8]Frese M H.MACH2:A two-dimensional magnetohydrodynamic simulation code for complex experimental configurations[R].AMRC-R-874,1987.
[9]MacGillivray J,Peterkin R E,Burke D A,et al.MACH3:A CHSSI code for computational magnetohydrodynamics of general materials in generalized coordinates[C]∥IEEE Computer Society.Proceedings of the users group2003conference.Bellevue,WA,2003:9-13.
[10]Frese S D,Frese M H.Recent improvements to MACH2and MACH3for fast Z-pinch modeling[C]∥Proceedings of 5th International Conference on dense Z-pinch.Albuquerque,NM,2002:380-383.
[11]L’Eplattenier P,Cook G,Ashcraft C,et al.Introduction of an electromagnetism module in LS-DYNA for coupled mechanical-thermal-electromagnetic simulations[J].Steel Research International,2009,80(5):351-358.
[12]王刚华.磁驱动准等熵压缩和高速飞片实验、计算和反积分数据处理技术[D].绵阳:中国工程物理研究院,2008.
[13]张恒第.MC-1型爆磁压缩发生器一维内爆:磁流体力学模拟[D].绵阳:中国工程物理研究院,2012.
[14]张旭平.磁驱动超高速飞片发射技术研究[D].绵阳:中国工程物理研究院,2012.
[15]孙承纬.一维冲击波和爆轰波计算程序SSS[J].计算物理,1986,3(2):142-154.Sun Chengwei.SSS:A code computing one dimensional shock and detonation wave propagation[J].Chinese Journal of Computational Physics,1986,3(2):143-154.