激光触发SF6-N2混合气体开关的数值模拟

 在气体放电物理的基础上，对SF₆和N₂采用双光子电离模型，对碳氢化合物的光电离采用3能级模型，并考虑了混合气体的热电离和激光对气体的欧姆加热作用，建立了激光触发SF₆-N₂混合气体开关的数值模型，模拟了激光触发SF6-N2混合气体多级间隙开关实验。激光触发延迟时间的计算值与实验结果符合较好。理论计算表明：激光触发SF₆-N₂混合气体间隙开关的延迟时间随SF₆含量的增加呈上升趋势，而随激光脉冲能量、充气压力等的上升呈下降趋势。     

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

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

等离子体断路开关的一维磁流体力学数值模拟

 采用磁流体力学的方法对等离子体端路开关进行了一维数值模拟，通过对给定的进入等离子体的电流波形模型的计算，给出了轴向上等离子体的运动及其不同时刻的密度分布，分析了磁场在等离子体内部的异常渗透过程，指出考虑Hall效应是产生这一异常现象的最为主要的因素。     

电爆炸丝1维磁流体模型数值模拟

 给出了电爆炸丝1维磁流体模型的具体描述和基本方程，并利用数值方法求解。通过对含电爆炸丝断路开关的电感储能脉冲功率电路模型的计算，给出了电爆炸丝电流、电压随时间的变化，以及电爆炸过程中丝的密度空间分布，膨胀半径随时间的变化，讨论了电爆炸丝长度和初始电压对电爆炸丝电流、电压的影响。     

准等熵压缩实验中金属铝的热力学耗散分析

为获取高压下材料的纯热力学压力-比容参考线和完全物态方程,减去应力-应变曲线中的其它信息,对准等熵压缩实验中由加载应变率引起的黏性耗散和热传导引起的热耗散做了分析讨论。基于反积分计算和流体动力学积分计算相结合的方法,根据激光加载(约108 s-1)和磁驱动准等熵压缩(约105 s-1)的实验数据,对材料声速、应力-应变曲线、温度和熵增等物理量进行计算,分析了不同应变率与该物理量的关系;还对热传导和SCG本构模型进行了计算,分析了热传导引起的温度变化对材料屈服强度、剪切模量和拉格朗日声速的影响。结果表明:激光加载实验中,应变率引起的温升差异约为180K,熵增差异约为250J/(kg·K),热传导导致温度下降40K;磁驱动准等熵压缩应变率较低,引起的熵增变化小于8J/(kg·K)。     

衬套内爆ALE方法二维MHD数值模拟

早期关于二维磁流体力学方程组的描述都是在Lagrangian或Eulerian坐标内进行的。在Lagrangian描述中，计算网格固定在物体上随物体一起运动，网格点与物质点在物体的变形过程中始终保持一致，因此能准确描述物体的移动界面。但对于大变形问题，物质的扭曲将导致计算网格的畸形而使计算无法进行下去。在Eulerian描述中，网格固定在空间中，因而计算网格在物体的变形过程中保持不变，因此很容易处理物质的扭曲。但对于运动界面，需要引入非常复杂的数学映射，将可能导致较大的误差。     

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

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

磁化套筒惯性聚变初步探索研究与实验可行性分析北大核心CSCD

磁化套筒惯性聚变(MagLIF)构型可充分利用现有大型脉冲功率驱动装置,如聚龙一号等。基于磁流体力学方程组和1∶1比例氘氚(DT)混合燃料聚变模型,开发了零维MagLIF数值模拟程序并进行了初步探索研究。计算结果表明初始负载参数(如轴向磁场强度,预加热温度、时刻,负载半径等)与聚变产额之间有着密切的联系,在给定条件下,可依据计算给出的定性关系进行负载优化设计。值得注意的是,根据计算结果,即使在理想条件下,氘氚燃料要实现能量收支平衡,则驱动器的电流必须大于21.2 MA。这意味着聚龙一号装置(10 MA)无法开展集成化的MagLIF实验,进一步的校验计算验证了上述观点,并在此基础上提出铝套筒分解实验的建议和负载设计参数。所取得的计算结果有利于加深对MagLIF套筒压缩阶段物理过程的认知和理解。     

双层钨衬套Z箍缩内爆数值模拟

 介绍了在强脉冲电流的作用下双层钨衬套Z箍缩动力学过程的数值模拟,计算了一个双层钨衬套模型,分析了双层衬套的碰撞过程,指出了双层衬套下电磁内爆的现象和特点,给出了碰撞后电流在内外衬筒上的分布,并对箍缩功率和箍缩产生的能量进行了预估。     

Load optimal design for a primary test stand facility based on a zero-dimensional load model

In order to couple the numerical simulation of a primary test stand driver with an optimal load design, a zero-dimensional wire array load model is designed based on the Saturn load model using PSPICE, which is an upgraded version of the Simulation Program with Integrated Circuit Emphasis (SPICE) designed by the ORCAD Corporation to perform circuit simulations. This paper calculates different load parameters and discusses factors influencing the driving current curve. With appropriate driving current curves chosen, further magneto-hydrodynamic calculations are carried out and discussed to provide the best results for experiments. The suggested optimal load parameters play an important role in experimental load design.