稀薄空气、等离子体环境外系统电磁脉冲的数值模拟

为研究大气环境对系统电磁脉冲(SGEMP)的影响,针对海拔50～100 km的X射线能量沉积区,分别应用3维PIC程序及3维PIC-MCC程序各自开展预电离等离子体和稀薄空气条件下外SGEMP的建模与模拟研究,针对3种不同的X射线注量(4×10^-3 J/cm²、 4×10^-2 J/cm²、 0.4 J/cm²),分别取对应两种不同海拔高度(70 km和80～90 km)的本底等离子体及海拔56 km的稀薄空气条件进行模拟计算,并和真空中的计算结果进行对比,得出预电离等离子体及稀薄空气对外SGEMP的影响规律：当X射线注量较低时,等离子体使得磁场增大,电场减小,而稀薄空气对外SGEMP效应影响不明显;随着X射线注量增大,空间电荷非线性效应越来越明显,等离子体及稀薄空气都使得电场、磁场同时增大,且稀薄空气的增大效应更显著。     

快点火中质子的能量沉积和神光Ⅱ升级装置上的质子束的产生

基于神光Ⅱ升级装置激光条件,利用流体程序、粒子模拟程序和Fokker-Placnck程序,模拟研究质子快点火中所需质子束的品质以及产生所需质子束的激光条件.首先根据快点火靶的条件,利用Fokker-Planck方程模拟快点火所需的质子束的能量范围,模拟表明当背景等离子密度为300 g/cm³时,能量为7—12 MeV的质子束适合点火;当背景等离子体密度为400 g/cm³时,能量为8—18 MeV的质子束适合点火.再根据神光Ⅱ升级装置实验条件研究质子束所需的激光参数,通过利用粒子模拟程序,结合流体程序给出的预等离子体,分别模拟研究了加预等离子体和不加预等离子体两种情况下的质子加速,在有预等离子体时得到的质子束最大能量约为22 MeV,没有预等离子体时得到的质子束最大能量为17.5 MeV,具体分析了两种情况下质子加速的物理机制,其结果跟等离子体自由膨胀模型结果符合得很好.     

类铝离子钟跃迁能级的超精细结构常数和朗德g因子

本文利用多组态Dirac-Hartree-Fock方法计算了类铝等电子序列从Si⁺到Kr²³⁺离子基组态3s²3p ²P_{1/2, 3/2}能级的超精细结构常数和朗德g因子.通过系统评估电子关联效应对Si⁺和Co¹⁴⁺离子中所关心原子参数的影响,尤其是与内壳层电子相关的关联效应,构建了可靠精确的计算模型,除Si⁺离子外,超精细结构常数和g因子的计算误差分别控制在1%左右和10^–5的量级.此外,进一步分析了超精细结构常数中电子部分矩阵元和g因子随原子序数Z的变化规律,并拟合了这些物理量与Z的定量依赖关系,利用拟合公式可以快速计算类铝离子在14≤Z≤54区间内任意同位素的超精细结构常数和g因子.     

欧拉流体力学数值方法中的界面计算

本文所提算法适用于二维和三维多介质流体力学两步欧拉数值方法中输运计算的混合网格(包括自由面网格)界面处理。在一个混合网格中，界面被近似地看作直线(二维)或平面(三维)。整个方法分为三步：(1)第一步，用混合网格周围的八个网格的介质面积份额(二维)或二十六个网格的介质体积份额(三维)确定界面的法线方向；第二步，用混合网格的本身的介质面积份额(二维)或体积份额(三维)确定界面的方程(位置)；第三步，用此直线方程求出通过网格边界的流以及下一时刻网格的面积份额(二维)或体积份额(三维)。最后给出了用此方法所做的一些数值计算及与SLIC算法的比较。     

长方体的热爆炸临界参数估计

本文给出了各侧面散热条件不同情况下的长方体反应系统临界F-K点火参数估计公式,并用已有的研究结果做了比较,表明该公式给出了较好的临界参数估计值。     

计算发散爆轰波传播的DSD方法

首先对爆轰冲击动力学（ＤＳＤ）方法的理论基础进行了探讨，完整地给出了它的数学表达，并研究了它的具体计算问题，指出了它的难点在于动边界问题的处理、然后针对两种特殊的几何条件，提出了用变量替换方法变动边界为定边界的思路。最后，根据计算结果，讨论了各种参数对波形的影响，特别是钝感炸药和敏感炸药所体现出的差异。     

非对称射流形成的几何理论(二)

将平板非对称正碰撞形成射流满足的定常方程组封闭条件推广到非正碰情况，提出了平面二维条件下非对称射流形成的一般理论，给出了射流形成的临界条件和出流射流的理论公式。对斜碰撞情形下非对称射流的出流方向、厚度随碰撞角和来流初始位形的变化规律进行了理论预测和数值模拟，理论解与二维欧拉程序的仿真结果符合较好。     

改进的 Morse 势和冲击压缩

用分子动力学方法模拟凝聚介质中一维分子链上的激波传播，以得到与冲击压缩实验相一致的结果。通过比较分子链上的孤立波解与激波的实验结果，提出对分子间Ｍｏｒｓｅ作用势进行与考虑量子效应相一致的修正，从而使分子动力学模拟得到的激波Ｕｓ－Ｕｐ关系与实验结果在波后粒子速度达到声速的范围内一致。     

爆轰波与金属的斜相互作用

应用二种截然不同的炸药PBX 95 0 1和T2 ,计算了爆轰波在不同入射角下与金属平板的斜相互作用。在正规反射区 ,计算结果与激波极线理论基本一致 ;应用燃烧模型 ,分别计算了Fe球壳装药JB90 0 3(HE)及JB90 14 (IHE)散心爆轰波的绕射传播 ,计算结果与实验很好地符合。在非正规反射区 ,二维拉氏程序计算结果明显地不同于经典理论结果 ,计算中没有出现Mach反射。计算结果显示 ,毗邻介质影响节点附近的爆轰波阵面形状及爆轰波速度 ;不同的反应率函数计算的节点图像不同。     

一种岩石损伤本构模型在地下强爆炸中的应用

为了描述地下爆炸波的传播,建立了一种考虑塑性硬化、剪胀和损伤软化效应的岩石本构模型,并应用于地下强爆炸自由场的数值计算,获得的速度和位移波形以及峰值速度衰减曲线等与国外地下试验测量数据和同类计算相比十分接近,从而验证了本构模型的有效性。