“聚龙一号”装置上磁驱动铝飞片实验的数值模拟
|
摘要
采用二维磁驱动数值模拟程序MDSC2,对大电流脉冲装置"聚龙一号"上的151发次磁驱动370μm厚铝飞片实验、164发次磁驱动330μm厚铝飞片实验进行了数值模拟和分析。数值模拟表明:370μm厚铝飞片和330μm厚铝飞片的磁驱动过程中,VISAR测量的速度不是飞片自由面的速度,而是飞片中邻近自由面最近的固体反射面的速度。这是由于磁驱动飞片发射过程中,飞片自由面部分被烧蚀,密度低于固体密度状态,而飞片自由面和加载面中间的飞片还保持固体密度状态。VISAR测量的激光将穿过自由面的低于固体密度状态的飞片部分,到达飞片自由面最近的固体密度位置再反射回去,获得这一位置的飞片速度。数值模拟的飞片固体反射面速度历史和VISAR测量的速度历史相吻合。
Magnetically driven flyer plate experiments,shot 151 with a 370μm thick aluminum flyer plate and shot 164 with a 330μm thick aluminum flyer plate,carried out in the large current pulse device PTS are simulated and analyzed with the two-dimensional magneto-hydrodynamics code MDSC2.Numerical simulation shows that the material near the free surface of the flyer plate melts or evaporates due to ablation,with the remaining part between the free surface and the loading surface staying solid.This finding tells us that the velocity measured by VISAR in the experiments may not be that of the free surface as expectation but be the one of the liquid-solid interface beyond the free surface,since the laser of the VISAR will penetrate through the melt or evaporated part and reflect back at the liquid-solid interface.This idea is confirmed by the coincidence of the simulated velocity of the liquid-solid interface and the measured one by the VISAR.
Magnetically driven flyer plate experiments,shot 151 with a 370μm thick aluminum flyer plate and shot 164 with a 330μm thick aluminum flyer plate,carried out in the large current pulse device PTS are simulated and analyzed with the two-dimensional magneto-hydrodynamics code MDSC2.Numerical simulation shows that the material near the free surface of the flyer plate melts or evaporates due to ablation,with the remaining part between the free surface and the loading surface staying solid.This finding tells us that the velocity measured by VISAR in the experiments may not be that of the free surface as expectation but be the one of the liquid-solid interface beyond the free surface,since the laser of the VISAR will penetrate through the melt or evaporated part and reflect back at the liquid-solid interface.This idea is confirmed by the coincidence of the simulated velocity of the liquid-solid interface and the measured one by the VISAR.
引文
[1]Lemke R W,Knudson M D,Robinson A C,et al.Self-consistent,two-dimensional,magneto-hydrodynamic simulations of magnetically driven flyer plates[J].Phys Plasmas,2003,10(5):1867-1874.
[2]Lemke R W,Knudson M D,Davis J P.Magnetically driven hyper-velocity launch capability at the Sandia Z accelerator[J].International Journal of Impact Engineering,2011,38(6):480-485.
[3]Matzen M K,Sweeney M A,Adams R G,et al.Plused-power-driven high energy density physics and inertial confinement fusion research[J].Phys Plasmas,2005,12:055503.
[4]Knudson M D,Hanson D L,Bailey J E,et al.Equation of state measurements in liquid deuterium to 70GPa[J].Phys Rev Lett,2001,87:225501.
[5]Lemke R W,Knudson M D,Bliss D E,et al.Magnetically accelerated,ultrahigh velocity flyer plates for shock wave experiments[J].J Appl Phys,2005,98:073530.
[6]Knudson M D,Lemke R W,Hayes D B,et al.Near-absolute Hugoniot measurements in aluminum to 500GPa using a magnetically accelerated flyer plate technique[J].J Appl Phys,2003,94(7):4420-4431.
[7]Knudson M D,Hanson D L,Bailey J E,et al.Use of a wave reverberation technique to infer the density compression of shocked liquid deuterium to 75GPa[J].Phys Rev Lett,2003,90:035505.
[8]Knudson M D,Hanson D L,Bailey J E,et al.Principal Hugoniot,reverberating wave,and mechanical reshock measurements of liquid deuterium to 400GPa using plate impact techniques[J].Phys Rev B,2004,69:144209.
[9]Lemke R W,Knudson M D,Hall C A,et al.Characterization of magnetically accelerated flyer plates[J].Phys Plasmas,2003,10(4):1092-1099.
[10]夏明鹤,计策,王玉娟,等.PTS装置工作模式及波形调节[J].物理学报,2012,57(5):3027-3037.(Xia Minghe,Ji Ce,Wang Yujuan,et al.Operation models and waveform shaping of primary test stand.Acta Physica Sinica,2012,57(5):3027-3037)
[11]阚明先,蒋吉昊,王刚华,等.套筒内爆ALE方法二维MHD数值模拟[J].四川大学学报,2007,44(1):91-96.(Kan Mingxian,Jiang Jihao,Wang Ganghua,et al.ALE simulation 2D MHD for liner.Journal of Sichuan University,2007,44(1):91-96)
[12]阚明先,王刚华,赵海龙,等.磁驱动飞片二维磁流体力学数值模拟[J].强激光与粒子束,2013,25(8):2137-2141.(Kan Mingxian,Wang Ganghua,Zhao Hailong,et al.Two dimensional magneto-hydrodynamic simulations of magnetically accelerated flyer plates.High Power Laser and Particle Beams,2013,25(8):2137-2141)
[13]阚明先,王刚华,张红平,等.磁驱动高速飞片模拟中滑移界面处理[J].强激光与粒子束,2015,27:015002.(Kan Mingxian,Wang Ganghua,Zhang Hongping,et al.Sliding interface processing in simulation on magnetically driving high speed flyer.High Power Laser and Particle Beams,2015,27:015002)
[14]阚明先,王刚华,赵海龙,等.金属电阻率模型[J].爆炸与冲击,2013,33(3):282-286.(Kan Mingxian,Wang Ganghua,Zhao Hailong,et al.Electrical resistivity model for metals.Explosion and Shock Waves,2013,33(3):282-286.
[15]Davis J P,Brown J L,Knudson M D,et al.Analysis of shockless dynamic compression data on solids to multi-megabar pressures:Application to tantalum[J].J Appl Phys,2014,116:204903.
[2]Lemke R W,Knudson M D,Davis J P.Magnetically driven hyper-velocity launch capability at the Sandia Z accelerator[J].International Journal of Impact Engineering,2011,38(6):480-485.
[3]Matzen M K,Sweeney M A,Adams R G,et al.Plused-power-driven high energy density physics and inertial confinement fusion research[J].Phys Plasmas,2005,12:055503.
[4]Knudson M D,Hanson D L,Bailey J E,et al.Equation of state measurements in liquid deuterium to 70GPa[J].Phys Rev Lett,2001,87:225501.
[5]Lemke R W,Knudson M D,Bliss D E,et al.Magnetically accelerated,ultrahigh velocity flyer plates for shock wave experiments[J].J Appl Phys,2005,98:073530.
[6]Knudson M D,Lemke R W,Hayes D B,et al.Near-absolute Hugoniot measurements in aluminum to 500GPa using a magnetically accelerated flyer plate technique[J].J Appl Phys,2003,94(7):4420-4431.
[7]Knudson M D,Hanson D L,Bailey J E,et al.Use of a wave reverberation technique to infer the density compression of shocked liquid deuterium to 75GPa[J].Phys Rev Lett,2003,90:035505.
[8]Knudson M D,Hanson D L,Bailey J E,et al.Principal Hugoniot,reverberating wave,and mechanical reshock measurements of liquid deuterium to 400GPa using plate impact techniques[J].Phys Rev B,2004,69:144209.
[9]Lemke R W,Knudson M D,Hall C A,et al.Characterization of magnetically accelerated flyer plates[J].Phys Plasmas,2003,10(4):1092-1099.
[10]夏明鹤,计策,王玉娟,等.PTS装置工作模式及波形调节[J].物理学报,2012,57(5):3027-3037.(Xia Minghe,Ji Ce,Wang Yujuan,et al.Operation models and waveform shaping of primary test stand.Acta Physica Sinica,2012,57(5):3027-3037)
[11]阚明先,蒋吉昊,王刚华,等.套筒内爆ALE方法二维MHD数值模拟[J].四川大学学报,2007,44(1):91-96.(Kan Mingxian,Jiang Jihao,Wang Ganghua,et al.ALE simulation 2D MHD for liner.Journal of Sichuan University,2007,44(1):91-96)
[12]阚明先,王刚华,赵海龙,等.磁驱动飞片二维磁流体力学数值模拟[J].强激光与粒子束,2013,25(8):2137-2141.(Kan Mingxian,Wang Ganghua,Zhao Hailong,et al.Two dimensional magneto-hydrodynamic simulations of magnetically accelerated flyer plates.High Power Laser and Particle Beams,2013,25(8):2137-2141)
[13]阚明先,王刚华,张红平,等.磁驱动高速飞片模拟中滑移界面处理[J].强激光与粒子束,2015,27:015002.(Kan Mingxian,Wang Ganghua,Zhang Hongping,et al.Sliding interface processing in simulation on magnetically driving high speed flyer.High Power Laser and Particle Beams,2015,27:015002)
[14]阚明先,王刚华,赵海龙,等.金属电阻率模型[J].爆炸与冲击,2013,33(3):282-286.(Kan Mingxian,Wang Ganghua,Zhao Hailong,et al.Electrical resistivity model for metals.Explosion and Shock Waves,2013,33(3):282-286.
[15]Davis J P,Brown J L,Knudson M D,et al.Analysis of shockless dynamic compression data on solids to multi-megabar pressures:Application to tantalum[J].J Appl Phys,2014,116:204903.