磁驱动飞片的一维磁流体动力学数值研究

磁驱动高速飞片技术是近年来发展的一种新型实验技术,在冲击波物理领域得到应用。该过程伴随着磁场扩散,并由此引起焦耳加热,使得飞片加载面的相状态发生变化,这决定了飞片厚度的范围。基于拉格朗日坐标系,利用磁流体动力学方程组、电阻率方程和状态方程数据库,对磁驱动铝飞片进行了一维磁流体动力学数值计算,获得了不同时刻铝飞片密度、温度的剖面分布,得到了磁场扩散速率随加载电流密度的变化关系。文章所选取的电导率方程只考虑到汽化点为止,对于等离子体形成的过程无法描述,如果要精确描述更高电流密度下的驱动过程,需考虑更为普适的电导率方程。磁场扩散速率随加载电流密度的变化存在转折点,在转折点前后可分别用两个线性关系表达式加以刻画。利用这些关系和冲击波物理相关知识,对磁压加载等熵驱动飞片实验样品厚度的选择进行了研究。     

斜波压缩下HMX晶体的弹塑性行为

开展了（010）、（011）晶向HMX晶体的斜波压缩实验,获得了约15 GPa压力下的速度响应剖面。实验结果表明,HMX单晶存在明显弹塑性转变行为,且速度波形有下降趋势,这是材料的黏性效应导致,材料的弹性极限随着样品厚度增加而变化,不同晶向的材料动力学特性存在差异。结合Hobenemser-Prager黏弹塑性本构关系和三阶Birch-Murnaghan物态方程开展了HMX晶体斜波压缩物理过程的数值模拟,计算结果可以很好地描述HMX晶体的弹塑性转变这一物理过程。     

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

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

准等熵压缩流场反演技术研究

 与一般流体力学计算中将流体力学方程组进行空间离散、再进行时间推进计算不同,采用Hayes提出的流场反演方法,将流体力学方程组进行时间离散,然后将VISAR测量得到的样品后自由面速度历史作为输入数据,再进行空间的反演计算。与一般流体力学计算进行了比较,并与柏劲松设计的Pillow飞片撞击铜样品实验的计算结果作了比较,该方法反演出的样品加载面处的压力历史与其计算结果基本一致,通过参数优化方法多次进行反积分计算调试出的等熵压缩线具有较高的可信度。     

磁驱动准等熵平面压缩和超高速飞片发射实验技术原理、装置及应用

利用脉冲大电流产生随时间平滑上升的磁压,实现对样品的准等熵平面压缩和超高速飞片发射,是近十年来发展和完善起来的一种新型的强动态斜波加载技术(ramp wave loading).本文简述了其原理、加载装置及数据处理方法等方面的研究进展,同时着重评述利用该技术和方法开展高压物态方程、材料动力学响应方面的研究进展,并对该技术在冲击动力学、天体物理和高能量密度物理等方面的应用前景进行了展望.      

窗口声阻抗对锆相变动力学的影响

基于磁驱动加载装置CQ-4开展了锆的斜波压缩相变实验,研究了锆样品后表面窗口声阻抗对相变波形的影响.实验结果显示,锆后表面为较低声阻抗窗口(自由面和LiF窗口)时,相变起始对应的特征粒子速度约331.0 m/s,而高阻抗蓝宝石窗口时,特征粒子速度约301.9 m/s,特征速度对应的压力从约9.14 GPa下降到8.27 GPa.相变对应的速度特征拐点是与多种因素相关的实验信息,因此它对应的压力并不是材料属性参数相变压力.结合基于热力学Helmholtz自由能的多相状态方程和非平衡相变动力学方程开展了锆的相变动力学数值模拟研究,相变弛豫时间为30 ns,计算结果与三种情况的实验结果符合良好,可以较好地模拟斜波压缩下锆的弹塑性转变、相变等物理过程.在压力-比容和温度-压力热力学平面,相变前锆的准等熵线与冲击绝热线差异很小,相变后准等熵线都位于冲击绝热线下方,随着压力的增加准等熵线和冲击绝线偏差越来越大,温度-压力平面中在20 GPa时相差约100 K.相变开始后,由于相变引起比容的间断,导致锆的拉氏声速迅速下降约7%,相变完成后拉氏声速恢复到体波声速.     

周期性铜线密排结构的冲击压缩特性研究

对铜线密排结构材料开展了一维平面应变冲击压缩实验与数值模拟研究。利用激光位移干涉仪测量了样品/窗口界面速度剖面,获取了有效的结构冲击压缩实验数据。从样品界面速度曲线可以推断,密排结构冲击压缩时没有形成致密结构,在卸载过程中发生了分层。应用光滑粒子流体动力学方法(SPH)建立了平面铜线密排结构的三维计算模型,计算得到样品冲击加载下的压缩特性,实验和数值模拟得到的界面速度和压力结果吻合较好,为后续开展柱面铜线密排结构的冲击压缩过程研究奠定了基础。     

D201×4树脂对金(Ⅲ)的吸附

大量的金是从金矿中获得,但需要很多工序。从含金的废液、废料中再生回收金也是获得金的重要来源,如:化学电镀法、电解法、硫脲法等都比较简单易行,但有各自的局限性和缺陷。20世纪70年代,     

较高场畸变系数的大电流三电极气体开关

气体开关作为脉冲功率装置的关键部件,其自击穿概率以及触发放电延时抖动对整个脉冲功率系统具有至关重要的影响。降低开关工作系数有利于提高开关稳定性,但延时抖动会随之增大。针对用于磁驱动实验的10 MA级大电流装置应用需求,设计了一种具有较高场畸变系数、能在较低工作系数条件下稳定工作的三电极气体开关,并开展了该开关的性能研究。模拟与实验结果表明：在触发电压与充电电压相当的条件下,开关的场畸变系数接近4,开关工作系数高于60%时,开关具有较低的延时抖动,抖动均方根小于3 ns。结合该开关设计了一个两级Marx储能模块,充电电压±50 kV条件下短路放电,模块回路放电电流峰值达到150 kA、周期2μs。上千次放电实验后,开关电极表面未发生明显烧蚀,工作正常。工作系数68.5%时,共计4 000发实验中未出现自放电现象,自击穿概率低于2.5×10-4。上述结果表明该开关可满足300～400只开关同时工作的大电流装置需求。     

斜波压缩下锡的相变动力学特性

利用磁驱动加载装置(CQ-4)和高精度速度测试装置(DPV),开展了斜波加载下锡的动态压缩实验。实验结果表明:锡在加载阶段经历了弹塑性转变和相变等物理过程,相变压力约为7.5 GPa。β–γ相变对应的特征速度随着锡厚度的增加,从676.3 m/s减小到636.8 m/s,对应的压力从7.62 GPa降低到7.11 GPa。结合Hayes多相状态方程和非平衡相变动力学模型,对锡的斜波压缩实验过程进行了模拟,数值计算结果可以较好地描述锡在加载阶段的弹塑性转变和相变等物理过程。讨论了体模量在不同热力学过程中的物理形式,计算结果显示,斜波压缩过程需考虑压力对体模量的修正。分析了相变弛豫时间、体模量等典型物理参数对速度波形的影响,结果表明,相变弛豫时间和各相初始自由能主要影响混合区部分速度波形,γ相的体模量参数只影响相变后的速度波形,β相的体模量参数会影响整体速度波形。