蒙特卡罗方法初步—椭圆内均匀分布的抽样技巧

以简明通俗的语言，避开统计学中抽象的数学公式和符号，介绍了作为计算数学一个独立分支的蒙特卡罗方法的基本思想和概念、蒙特卡罗方法在粒子输运问题中的应用以及最基本的几种抽样方法，最后介绍一种椭圆内均匀分布的蒙特卡罗“压缩”抽样方法及其证明。     

强流相对论电子束的数值模拟

我们研制了CEBQ程序和ENVLOP程序，模拟稳态相对论电子束在二极管区的产生和加速，以及电子束在阳极通道内的传输过程。CEBQ程序模拟二极管区电子轨道，当阳极电压为1．5MV时，阴极总发射电流为17．9kA，进入阳极通道内的传输电流为4．3kA，阳极入口处D=4cm的束发射度为24．9cm．mrad。ENVLOP程序使用线性近似下的包络方程，描述电子束包络在螺旋管引导磁场作用下的变化，讨论了外加引导磁场对包络的影响。计算结果表明，在合适引导磁场作用下，束包络起伏较小，进入阳极通道的束流(4．3kA)能全部通过。     

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

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

Simulation for double shell pinch

Basic shock phenomena are presented in a composite pinch, a hybrid of the Z-pinch. The successive transfer of current within the plasma structure is demonstrated by our calculations. Properties of the shock wave are described.The current distribution between the two shells after the outer shell hitting the inner shell is also discussed.     

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

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

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

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

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

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

单丝Z箍缩的一维单温辐射磁流体力学计算

 介绍在正弦电流波形下铝丝的内爆动力学过程的数值模拟计算，计算采用一维单温辐射磁流体力学模型，能流中除了常规热传导外，也包括黑体辐射流项，以考虑辐射对内爆动力学过程的影响。计算了用于箍缩预脉冲的低电流（幅值10～50 kA）和较高电流（幅值400 kA）情况下单丝的内爆结果。     

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

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