数值模拟辐射驱动2维内爆压缩过程

采用2维三温非平衡辐射流体力学程序LARED-H数值模拟辐射驱动内爆过程。针对2维三温能量方程九点差分格式离散后的线性方程组,采用了高效的Krysolv子空间迭代解法,改进了代数解法器。将1维间接驱动内爆总体程序CFJ与LARED-H程序的计算结果进行比对,验证了LARED-H程序数值模拟1维内爆问题的正确性。并数值模拟了不同腔长辐射温度源驱动下的2维靶球运动,数值结果显示:随着腔长的增加,高压缩内爆燃料区分别被压缩成香肠形、球形和铁饼形,数值模拟结果与神光Ⅱ的实验结果定性上相同。     

基于子网格边界近似Riemann解的人为粘性方法

在交错网格型Lagrange(拉格朗日)流体力学算法中,通常采用人为粘性捕捉激波,人为粘性的好坏对于计算结果至关重要.研究了一种基于子网格边界处近似Riemann解的新型人为粘性.新人为粘性能够满足动量守恒和熵不等式.利用子网格边界速度差中引入的限制器,新人为粘性能够区分激波和等熵压缩,并能满足球对称问题中的波面不变性.新人为粘性在典型模型数值模拟及惯性约束聚变黑腔整体数值模拟中取得了较好的结果.     

点火内爆靶辐射驱动不对称性全过程数值模拟

为满足分层掺杂点火内爆靶辐射驱动不对称性全过程物理分析的需求,在激光聚变二维总体程序LARED集成上发展了辐射输运建模下的多介质ALE方法-RTALE(Radiative Transfer Arbitrary-Lagrangian-Eulerian)。为提高多介质ALE方法的健壮性,发展了驰豫网格重构算法,该重构算法生成的新网格能自适应流场的变化。数值模拟了激波与气柱相互作用的RM不稳定性实验,模拟的气泡变形程度与试验结果基本一致,其中驰豫网格重构算法中的驰豫因子能够很好地反映流场密度梯度。基于辐射多群输运建模的LARED集成程序能够完整模拟辐射驱动不对称性条件下掺杂点火靶二维内爆过程,克服了传统ALE方法计算不下去和算不好的困难,界面变形程度也符合物理分析。     

岩石、混凝土和土抗侵彻能力数值计算与分析

利用LS-DYNA3D软件数值计算了弹体侵彻岩石、混凝土和土问题,分析在不同碰撞速度条件下的弹体响应和靶体抗侵彻能力。碰撞速度小于900 m/s时,弹体侵彻岩石的减加速度峰值约是侵彻混凝土的2倍,而侵彻混凝土的减加速度峰值约是侵彻土的6倍。减加速度峰值高则稳态侵彻过程短,弹体能量消耗很快。碰撞速度超过1.5 km/s时,随靶体材料的强度、密度逐渐减小,侵彻深度和孔径逐渐缓慢增加,岩石、混凝土和土3种靶体材料相比,最大侵彻深度增加41%~62%,最大扩孔口径增加16%~25%。     

斜入射激光平面靶耦合的数值模拟

LARED-H是一个柱对称二维三温非平衡辐射流体力学程序，主要用于黑腔物理研究，也可模拟正入射激光平面靶耦合过程，但无法模拟斜入射平面靶问题。实际上斜入射现象会经常发生，因此需要扩充LARED-H程序在直角坐标系下的功能，用于模拟斜入射激光平面靶耦合过程。引进拟直角坐标系，把直角坐标系与柱坐标系下的总体方程的差分格式统一，实现了直角坐标系下LARED-H程序的模拟能力。     

靶丸变形实验的多群扩散整体数值模拟

神光Ⅱ装置上靶丸压缩实验中,黑腔长度的变化会影响靶丸赤道和两极受到的辐照强度,从而导致燃料区被压缩成不同的形状。介绍了辐射多群扩散建模下的辐射流体力学方程组及其数值方法。采用最近研制的二维总体LARED集成程序,对神光Ⅱ不同长度黑腔靶丸压缩变形实验进行了整体数值模拟。结果表明辐射多群扩散建模可以反映腔内辐照均匀性变化情况,靶丸压缩形状与神光Ⅱ实验定性一致。     

保持对称性的新积分梯度格式

针对柱坐标系下拉氏流体力学的动量方程,提出一种积分梯度格式IGTSP(Integral Gradient Total Symmetry-Preserving),它具备现有积分梯度格式IGA(Integral Gradient Average)和IGT(Integral Gradient Total)的优点,不仅克服了IGT格式不能保持柱坐标系下的一维球对称性的缺点,而且系统的总动量守恒误差为O(h),比IGA格式更好地保持系统的总动量守恒.数值试验进一步显示了该格式理论分析的优点.     

激光腔靶耦合二维辐射输运数值模拟

间接驱动激光聚变过程中,黑腔内物质处于非局域热动平衡(non-LTE)状态,而且辐射传输具有非平衡、各向异性等特点。为了精确描述黑腔辐射场的演化及其与物质的相互作用,最新研制的激光聚变二维总体LARED集成程序,基于non-LTE的多群辐射输运建模,首次实现了激光黑腔靶实验的全过程数值模拟。数值结果表明,辐射输运计算较好地反映了黑腔辐射场均匀性变化,腔壁光斑区与非光斑区X光发射强度比与实验测量值接近,靶丸压缩形状与实验图像定性一致。     

激光间接驱动惯性约束聚变二维总体程序——LARED集成程序

介绍了辐射流体力学程序LARED集成程序的物理背景、模型方程、数值方法和数值算例。该程序主要应用于激光间接驱动惯性约束聚变的二维整体模拟,兼顾激光直接驱动、辐射驱动靶丸内爆过程和流体不稳定性等物理过程的数值模拟。通过与实验数据、一维辐射流体力学程序进行比对,验证了程序的可靠性。该程序实现了多群输运建模下NIF点火靶的全过程数值模拟,并已应用于惯性约束聚变的物理研究。     

自适应于流场变化的网格重构方法

提出九点网格重构方法,不仅生成的网格正交性和光滑性好,而且重构后的网格与物理问题解的空间分布相适应,实现了重构后的网格品质好.理论上证明了九点网格重构方法逼近椭圆型方程,并具有球对称性,数值试验表明,该方法具有很好的适应复杂区域的稳健性,在爆轰驱动或高速冲击ALE(Arbitrary Lagrangian Eulerian)数值模拟中具有较高的应用价值.