Pillow飞片气炮加载实验的数值模拟研究

 采用多介质弹塑性流体动力学计算方法,研制了适用于复杂加载的Lagrange计算程序MLEP,对劳伦斯利弗莫尔国家实验室（Lawrence Livermore National Laboratory,LLNL）19层和冲击波物理与爆轰物理国防科技重点实验室（Laboratory of Shock Wave and Detonation Physics,LSD）设计的29层Mg-Cu体系Pillow密度梯度飞片气炮加载实验过程进行了数值模拟和比较,获得的速度剖面计算结果与实验测试结果吻合一致,验证了流体动力学计算方法、不同材料体系混合模型以及计算程序的有效性和实用性,为进一步开展可控路径的复杂加载实验研究奠定了基础。     

冲击加载-准等熵加载过程的密度梯度飞片计算设计

 采用一维弹塑性流体动力学计算方法,通过对LLNL（Lawrence Livermore National Laboratory）Mg-Cu体系密度梯度飞片冲击加载-准等熵加载实验过程数值计算和比较,验证了流体动力学计算方法、不同材料体系混合模型以及计算程序的正确性和实用性。考虑到在飞片材料制备中,单层厚度最小值为0.2 mm和Mg-W最大阻抗混合质量分数的致密条件限制,开展了对新的Mg-W体系密度梯度飞片实现冲击加载-准等熵加载过程计算设计,给出了满足加载要求的飞片结构特征;从计算给出的粒子速度波剖面可见,密度梯度飞片波阻抗分布对加载过程和加载强度非常敏感,通过精心设计准连续型变阻抗的梯度飞片,可以进行不同复杂加/卸载过程的物理模型设计和实验研究。     

渐变阻抗飞片超高速发射数值模拟

 为在轻气炮上实现二级飞片的超高速发射,设计了波阻抗连续变化的渐变阻抗飞片。变阻抗飞片撞击面为镁,底面为钨,中间为镁、铜和钨的掺合材料。建立了变阻抗飞片超高速发射的计算模型,采用分成单元层描述阻抗变化。对变阻抗飞片超高速发射二级飞片过程进行了数值模拟计算,分析了阻抗分布、二级飞片直径和缓冲层厚度对二级飞片速度的影响。结果表明：在优化阻抗分布、二级飞片直径和缓冲层厚度的条件下,二级飞片的最大速度能够达到16 km/s以上。     

阻抗梯度飞片加载下的超高速发射二维数值模拟方法

阻抗梯度飞片准等熵加载和超高速发射的二维数值模拟,在计算方法上集中体现了多介质、多界面、大变形、高密度比等特点.采用多介质流体高精度PPM计算方法,以VOF为基础研制MFPPM2二维计算程序,数值模拟Sandia实验室的超高速发射实验模型,获得了与Sandia实验室数值计算一致的结果.     

激光驱动飞片过程中激光烧蚀深度的计算

利用激光体烧蚀模型,数值模拟了激光驱动飞片的加速过程,包括激光的吸收和飞片的速度历史等。在光强为GW/cm2量级的激光作用下,激光烧蚀产生的等离子体的流体力学运动可用改编的1维Lagrange流体力学计算程序SSS来描述。通过计算得到不同激光能量下的飞片密度剖面,由此给出金属薄膜的烧蚀深度与实验测量值进行比对,二者符合得较好。     

磁驱动飞片的三维数值模拟及分析

磁驱动准等熵加载和超高速飞片发射是一种全新的冲击动力学和高能量密度物理实验加载技术。利用三维磁流体动力学软件,模拟了磁驱动飞片的物理过程,计算得到的飞片自由面速度与实验结果符合较好。通过计算飞片横断面的温度、密度和磁场分布,得到了加载过程中磁扩散速度和飞片的剩余厚度。飞片加载过程中飞片边缘的卷曲变形严重,分析认为是由电流和磁场分布的不均匀导致飞片边侧受斜上方较大的加载力所致,并且电流分布的不均匀是主要因素。实验设计时,可利用极板构型的变化调节加载面的电流分布,从而提高飞片的平面性,减小边侧的卷曲变形。