自由面被烧蚀磁驱动飞片的数值模拟

为了理解磁驱动飞片自由面烧蚀的机理、自由面被烧蚀后飞片的物质状态、自由面被烧蚀后激光速度干涉仪(VISAR)测量的机理等,采用磁驱动数值模拟程序MDSC2对聚龙一号装置上PTS-151发次磁驱动飞片实验中370 m厚飞片进行了模拟和分析。数值模拟表明,334 ns之前飞片自由面部分保持固体状态,334 ns之后飞片自由面部分已经熔化,到340 ns后整个飞片都被熔化。飞片自由面烧蚀的主要机制是电流焦耳加热,热扩散和压缩做功的贡献很小。数值模拟的固体密度反射面速度历史和VISAR测量的速度历史一致。飞片自由面熔化后,VISAR测量的速度是距离自由面最近的固体密度反射面的速度。     

磁驱动飞片二维磁流体力学数值模拟

针对磁驱动高速飞片发射技术,建立了磁驱动飞片的二维磁流体力学数学模型,并考虑了焦耳热对飞片的影响。在四边形网格的基础上,采用算子分裂法,把磁流体力学方程依次分成热扩散、磁扩散、理想流体力学三个物理过程进行求解,研制了磁驱动飞片二维磁流体力学数值模拟程序。对美国Sandia国家实验室Z装置上的一个磁驱动飞片发射实验进行了数值模拟,并分析了不同时刻焦耳加热对飞片的烧蚀情况,计算得到的飞片自由面速度曲线与实验激光速度干涉仪测量的结果相吻合。     

磁驱动飞片二维磁流体力学数值模拟

针对磁驱动高速飞片发射技术,建立了磁驱动飞片的二维磁流体力学数学模型,并考虑了焦耳热对飞片的影响。在四边形网格的基础上,采用算子分裂法,把磁流体力学方程依次分成热扩散、磁扩散、理想流体力学三个物理过程进行求解,研制了磁驱动飞片二维磁流体力学数值模拟程序。对美国Sandia国家实验室Z装置上的一个磁驱动飞片发射实验进行了数值模拟,并分析了不同时刻焦耳加热对飞片的烧蚀情况,计算得到的飞片自由面速度曲线与实验激光速度干涉仪测量的结果相吻合。     

“聚龙一号”装置上磁驱动铝飞片实验的数值模拟

采用二维磁驱动数值模拟程序MDSC2,对大电流脉冲装置聚龙一号上的151发次磁驱动370 m厚铝飞片实验、164发次磁驱动330 m厚铝飞片实验进行了数值模拟和分析。数值模拟表明：370 m厚铝飞片和330 m厚铝飞片的磁驱动过程中, VISAR测量的速度不是飞片自由面的速度,而是飞片中邻近自由面最近的固体反射面的速度。这是由于磁驱动飞片发射过程中,飞片自由面部分被烧蚀,密度低于固体密度状态,而飞片自由面和加载面中间的飞片还保持固体密度状态。VISAR测量的激光将穿过自由面的低于固体密度状态的飞片部分,到达飞片自由面最近的固体密度位置再反射回去,获得这一位置的飞片速度。数值模拟的飞片固体反射面速度历史和VISAR测量的速度历史相吻合。     

磁驱动飞片发射实验结构系数初步研究

为了确定磁驱动飞片发射实验结构系数的范围、影响因素、结构系数与影响因素的关系,对聚龙一号装置上的磁驱动飞片发射实验进行了数值模拟和分析。数值模拟表明,磁流体力学程序能正确模拟聚龙一号装置上各个磁驱动飞片发射实验;磁驱动双侧飞片发射实验的结构系数为0.7～0.8;磁驱动单侧飞片发射实验的结构系数为0.80～0.85。磁驱动飞片发射实验的结构系数与实验加载电流无关,仅由磁驱动飞片发射实验的负载结构决定。磁驱动飞片发射实验的结构系数取决于阴阳电极极板的初始宽度、阴阳电极之间的初始间隙以及阴阳电极上飞片厚度之和等三个因素。在磁驱动飞片发射实验中,电极初始宽度、阴阳电极之间的初始间隙不变的情况下,结构系数由阴阳电极上飞片厚度之和确定,阴阳电极上飞片厚度之和越大,结构系数越大。     

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

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

“阳”加速器磁驱动平面飞片实验和速度计算校验

脉冲功率驱动源作为磁驱动加载的重要手段,通过调整其电路参数可调节负载电流波形,实现对样品无冲击准等熵加载。在阳加速器上,开展了一系列轴对称结构和带状结构构型的磁驱动平面飞片发射实验,电极材料采用不锈钢和LY-12铝。实验中测量了进入负载的电流历史和电极后自由面速度历史,并通过时序控制将两者时间关联起来。本文以测量到的电流历史数据为基础,引入负载电流分布系数,并结合已知的LY-12铝的状态方程数据,计算电极后自由面速度历史和飞片速度历史。通过实验测量自由面速度历史校验负载各个位置的电流分布系数。另外,基于装置参数和实验数据确定了考虑负载电感变化的装置等效电路模型,形成了计算样品压力加载历史和电极后自由面速度历史估算程序。此外,初步分析不同厚度电极的自由面速度历史,获取了电极材料的准等熵加载波剖面信息,观察到一系列准等熵加载下材料动力学性能引起的物理现象。     

“阳”加速器磁驱动平面飞片实验和速度计算校验

脉冲功率驱动源作为磁驱动加载的重要手段,通过调整其电路参数可调节负载电流波形,实现对样品无冲击准等熵加载。在"阳"加速器上,开展了一系列轴对称结构和带状结构构型的磁驱动平面飞片发射实验,电极材料采用不锈钢和LY-12铝。实验中测量了进入负载的电流历史和电极后自由面速度历史,并通过时序控制将两者时间关联起来。本文以测量到的电流历史数据为基础,引入负载电流分布系数,并结合已知的LY-12铝的状态方程数据,计算电极后自由面速度历史和飞片速度历史。通过实验测量自由面速度历史校验负载各个位置的电流分布系数。另外,基于装置参数和实验数据确定了考虑负载电感变化的装置等效电路模型,形成了计算样品压力加载历史和电极后自由面速度历史估算程序。此外,初步分析不同厚度电极的自由面速度历史,获取了电极材料的准等熵加载波剖面信息,观察到一系列准等熵加载下材料动力学性能引起的物理现象。     

磁驱动高速飞片模拟中滑移界面处理

在二维单介质磁驱动数值模拟程序(MDSC2)的基础上,采用局部结构整体非结构的网格拼接方法,研制了处理滑移界面的二维多介质程序。二维多介质程序解决了单介质程序无法处理的电极与飞片之间相互作用的模拟问题,能更准确地模拟磁驱动发射实验的飞片状态。计算结果表明:发射飞片的自由面部分始终保持固体密度状态;自由面速度历史和VISAR测量的速度曲线相吻合;飞片中间部分在电流加载过程中始终具有良好的平面性,飞片两端出现拖尾质量,产生不稳定性,这是由于飞片两端与飞片中间部分的磁场强度分布不同造成的。     

磁驱动高速飞片模拟中滑移界面处理

在二维单介质磁驱动数值模拟程序（MDSC2）的基础上,采用局部结构整体非结构的网格拼接方法,研制了处理滑移界面的二维多介质程序。二维多介质程序解决了单介质程序无法处理的电极与飞片之间相互作用的模拟问题,能更准确地模拟磁驱动发射实验的飞片状态。计算结果表明：发射飞片的自由面部分始终保持固体密度状态;自由面速度历史和VISAR测量的速度曲线相吻合;飞片中间部分在电流加载过程中始终具有良好的平面性,飞片两端出现拖尾质量,产生不稳定性,这是由于飞片两端与飞片中间部分的磁场强度分布不同造成的。     

脉冲激光烧蚀Ge产生等离子体特性的数值模拟

针对激光烧蚀半导体材料Ge初期的特点,建立了1维的热传导和流体动力学模型。对波长为248 nm、脉宽为17 ns、峰值功率密度为4×108 W/cm2的KrF脉冲激光在133.32 Pa氦气环境下烧蚀Ge产生等离子体的特性进行了数值模拟。结果表明：单个激光脉冲对靶的烧蚀深度达到55 nm,蒸气膨胀前端由于压缩背景气体产生压缩冲击波, 波前的速度最大,温度很高。从不同时刻的电离率分布图中得出,在靶面附近区域,Ge的1阶电离始终占优势;在中心区域,脉冲作用时间内,Ge的2阶电离率比1阶电离率大,脉冲结束后,Ge的2阶电离率下降,1阶电离率逐渐变大。     

脉冲激光烧蚀Ge产生等离子体特性的数值模拟

 针对激光烧蚀半导体材料Ge初期的特点,建立了1维的热传导和流体动力学模型。对波长为248 nm、脉宽为17 ns、峰值功率密度为4×10⁸ W/cm²的KrF脉冲激光在133.32 Pa氦气环境下烧蚀Ge产生等离子体的特性进行了数值模拟。结果表明：单个激光脉冲对靶的烧蚀深度达到55 nm,蒸气膨胀前端由于压缩背景气体产生压缩冲击波, 波前的速度最大,温度很高。从不同时刻的电离率分布图中得出,在靶面附近区域,Ge的1阶电离始终占优势;在中心区域,脉冲作用时间内,Ge的2阶电离率比1阶电离率大,脉冲结束后,Ge的2阶电离率下降,1阶电离率逐渐变大。