电磁内爆X射线数值计算

 叙述一维辐射磁流体力学方程组（三温）的数值模拟,用于电磁内爆X射线计算。动量方程、三温方程组采用分裂法求解,以克服强耦合非线性可能引起的数值不稳定性。采用隐式差分格式避免小时间步长限制。隐式格式采用Newton-Raphson法迭代求解。使用了分层网格、加权平均,迭代初值公式,自动调整时间步长,奇点处理等特殊处理方法。最后给出Pegasus l装置一个实验模型的计算结果,获得15 J的X射线能量,峰值功率为16 GW。     

电磁内爆套筒屈曲的数值模拟

 采用瞬态非线性有限元方法对电磁内爆套筒屈曲进行了数值模拟研究。介绍了瞬态非线性有限元方法的基本原理、电磁内爆套筒模拟的基本方法及对实验的模拟结果。通过分析比较认为,建立的数值模拟方法较好地反映了实验结果,可以作为研究分析电磁内爆套筒屈曲的基本方法之一。     

电磁内爆的准一维数值模拟

 提出一个用于电磁内爆过程数值模拟的准一维模型，对1 MJ电容器组电磁内爆优化方案进行了计算，并将计算结果与零维模型和一维模型的结果进行了比较。     

电流脉冲前沿对电磁驱动固体套筒内爆的影响

固体套筒内爆是采用实验方法研究高能量密度状态下的材料力学性能的重要加载手段之一,国内已经建立起若干开展电磁内爆研究的驱动器。从电流脉冲前沿对固体套筒内爆性能影响的角度进行分析,为如何选择现有的实验装置开展固体套筒内爆实验研究提供依据。采用不可压缩零维模型进行计算,获得了套筒内爆速度受套筒尺寸、电流幅值以及电流脉冲前沿的影响情况。计算结果表明,开展固体套筒内爆的实验研究应选择电流脉冲前沿大于2 s的装置,这也为未来设计驱动能力更强的固体套筒内爆实验装置奠定了基础。     

电磁内爆的一维显式数值模拟

 提出了一个用于电磁内爆过程数值模拟的一维磁流体力学显式模型，对中物院流体物理研究所1 MJ电容器组电磁内爆优化方案进行了计算，并将计算结果与一维德式模型的结果进行了比较。     

Z箍缩内爆的MARED程序数值模拟分析

MARED程序是模拟Z箍缩内爆过程的二维三温辐射磁流体力学程序,它适用于不同装置条件和不同负载参数.利用MARED程序对Z箍缩内爆进行模拟,结合丝阵Z箍缩实验分析表明:相同负载质量条件下,钨丝阵内爆产生的X射线辐射功率远大于铝丝阵产生的X射线功率;相同负载电流条件下,负载质量越大,计算得到的X射线功率越低;X射线功率随着负载电流增加而增加.MARED程序能够较好地反映Z箍缩内爆动力学过程,特别是不稳定性发展的二维图像,能够给出与不稳定性简化模型的理论分析及实验结果定性一致的演化规律.MARED程序模拟丝阵填充泡沫形成辐射场的初步计算得到了与Sandia实验室模拟Z装置上丝阵填充泡沫定性一致的结果.     

电磁内爆模拟技术实验研究进展

开展了中俄Z—pinch联合实验和强光一号实验，重点对双层丝阵内爆动力学模式、单层丝阵填充泡沫构建动力学黑腔及能谱相对强度的测量等进行了研究，取得了一批重要的物理研究成果。实现了双层丝阵内外层等离子体的同时内爆聚心，产生5．6TW的峰值辐射功率，     

点火靶高熵内爆的数值模拟

使用一维多群输运程序RDMG与二维少群扩散程序LARED-S对点火靶高脚与低脚内爆进行数值模拟.相对于低熵内爆,高脚高熵内爆通过提高预脉冲的辐射温度使得烧蚀面与物质界面的流体稳定性得到明显的改善,能够抑制流体不稳定的增长与热斑混合的发展.同时,高熵设计导致燃料的压缩变差,阻滞时刻燃料的压缩密度与面密度相应降低,中子产额降低.因此,高脚高熵内爆是通过牺牲燃料的高压缩,来换取靶丸内爆流体稳定性能的改善.     

FP-1装置铝套筒内爆动力学过程的一维磁流体力学模拟

作为一种重要的柱面会聚冲击和准等熵压缩加载源,磁驱动固体套筒内爆技术已广泛应用于高能量密度物理实验研究.针对FP-1装置驱动的固体套筒内爆动力学过程,建立了含强度的一维磁流体力学模型,并对典型实验进行了模拟.计算获得的套筒内爆速度同实验结果较为相符.模拟结果显示,该装置在40 kV充压条件下,可以将直径3 cm,厚0.5 mm的铝套筒加速至1.1 km/s,内壁速度超过1.5 km/s,同时保持大部分材料为固体状态.内爆套筒与相同材料靶筒碰撞产生的冲击压力约9 GPa.改变靶筒内部填充气体的压力,可以获得不同的靶筒运动速度、轨迹以及反弹半径,以满足不同类型实验的研究需要.     

二维弹塑性磁流体力学数值模拟

为了研究物质弹塑性对磁驱动实验运动过程、不稳定性发展等的影响,在MDSC2程序的基础上,增加了弹塑性模块,研制了包括弹塑性的磁流体力学程序,并进行了弹塑性项影响的数值模拟和分析。数值模拟表明：没有初始扰动时,弹塑性项几乎不影响套筒内外半径的运动轨迹;有初始扰动时,弹塑性项对磁驱动固体套筒的Rayleigh-Tayor不稳定性有明显的抑制作用。     

Numerical Simulation of 2-D Radiation-Drive Ignition Implosion Process

A Lagrangian compatible radiation hydrodynamic algorithm and the nuclear dynamics computing module are developed and implemented in the LARED Integration code, which is a radiation hydrodynamic code based on the 2-D cylindrical coordinates for the numerical simulation of the indirect-drive Inertial Confined Fusion. A number of 1-D and 2-D ignition implosion numerical simulations by using the improved LARED Integration code (ILARED) are presented which show that the 1-D numerical results are consistent with those computed by the 1-D radiation hydrodynamic code RDMG, while the simulation results of the 2-D low-mode radiative asymmetry and hydrodynamic instability growth,according to the physical analysis and anticipation, are satisfactory. The capsules driven by the sources from SGII experiments are also simulated by ILARED, and the fuel shapes agree well with the experimental results. The numerical simulations demonstrate that ILARED can be used in the simulation of the 1-D and 2-D ignition capsule implosion using the multi-group diffusion model for radiation.     

聚心火焰在共振腔作用下引发爆轰的数值研究

 基于带化学反应的二维轴对称Euler方程，利用带有Superbee限制函数的波传播算法，对共振腔中的氢气-空气预混气的聚心燃烧进行了数值模拟，讨论了共振腔不同抛物面对起爆的影响。数值结果表明，在开始阶段，燃烧诱导的激波在轴心、火焰和固壁的反射，使火焰失稳，随后共振腔中的抛物壁面上产生一定频率和强度的反射激波，不断穿越火焰，使火焰进一步失稳，加剧了燃烧速度，最终导致爆轰的形成。同时，火焰在与激波的作用过程中，形状扭曲变形，呈封闭端小敞口端大的扁平头部蘑菇云。共振腔抛物面的不同形状引起激波聚焦位置的变化，会影响激波和火焰的相互作用，使起爆提前或推迟，甚至不起爆。     

增强后效复合药型罩结构的数值模拟

针对增强聚能射流的破甲后效问题,设计了等壁平顶锥形铜铝双层复合药型罩装药结构,采用冲击波物理显示欧拉动力学软件SPEED开展复合射流成型及对钢-铝间隔靶侵彻过程的数值模拟,分析内外双层药型罩高度比ε、药型罩锥角α等参数对复合射流成型和间隔靶侵彻性能的影响规律。研究结果表明:复合射流的头部速度随ε增大呈先减小后增加的趋势,在ε约为1/2时,可形成具有相近速度的铜铝同轴复合射流微元,利于铝射流微元与目标相互作用实现后效增强毁伤;且当α在50°~60°范围内时,复合射流中段为集中的铝射流微元,更利于侵彻后的爆炸或爆燃反应。对优化参数的复合药型罩结构数值模拟结果与文献公布的实验结果吻合较好。研究结果对增强后效聚能装药设计具有参考价值。     

辐射驱动内爆温度源脉冲整形研究

辐射驱动温度源脉冲整形是获得高压缩比的重要途径。本文从脉冲整形的基本原理出发，采用数值模拟方法，给出适合于神光Ⅲ装置内爆实验的辐射驱动温度源脉冲整形曲线。     

靶丸支撑膜对ICF内爆性能影响的模拟研究

使用二维多群辐射扩散流体力学程序LARED-S,模拟研究靶丸支撑膜在惯性约束聚变氘氘（DD）气体靶内爆过程中的扰动演化过程及其对内爆性能的影响.二维模拟表明：靶丸支撑膜显著降低DD气体靶内爆的中子产额,二维模拟产额为一维结果的55.2%.内爆性能下降的主要物理机制是支撑膜使靶壳生成大幅度的尖钉深入DD气体区,使烧蚀层与DD气体之间物质界面处的电子热传导漏失功率显著增大,导致DD核反应速率显著降低,中子反应速率峰值时刻（bang-time时刻）提前.相比一维理想内爆的模拟结果,支撑膜引入的扰动显著降低bang-time时刻DD气体压强与内爆动能转化为DD气体内能的效率,壳层剩余动能相应大幅增加.     

不同罩材聚能装药对多层介质侵彻的实验与数值模拟

 为实现聚能装药对多层介质的大破孔侵彻,提出了钛合金药型罩聚能装药设计方案。采用实验与数值模拟相结合的方法,对钛合金、低碳钢及紫铜罩聚能装药侵彻多层介质进行了研究,分析了钛合金聚能侵彻体相对于紫铜和低碳钢侵彻体在成型过程中,其动能、头部速度及射流长度等的差异,并对侵彻过程中应力波的传播特性进行了分析。结果表明：相对于紫铜和低碳钢,钛合金罩聚能侵彻体的能量转换率高,所获得的动能大,头尾速度梯度小,外形更为短粗;虽对多层介质侵彻时侵彻深度有所减小,但漏斗坑尺寸明显增大,且平均破孔孔径提高了约20%。     

惯性约束聚变中内爆混合的模型构建

从热斑质量方程和能量守恒方程入手,重新计算考虑混合后聚变燃料的比内能和比热容等热力学参数,分析混合效应在轫致辐射损失等能量输运方面的作用,构建有杂质混合情况的热斑燃烧动力学模型.根据静态模型中的热斑燃烧的功率平衡条件,研究烧蚀层杂质混合比例与点火阈值和热斑自持燃烧的关系.理论分析和数值计算表明,混合效应导致热斑中的轫致辐射增强是点火失败的重要因素之一.通过调整不同掺杂材料、混合浓度及混合方式,得到壳层混合与热斑面密度、热斑离子温度的演化之间的关系.最后,基于模拟结果给出两种降低混合影响的方法.