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

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

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

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

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

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

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

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

不同辐射建模下球腔的整体数值模拟

不同辐射建模对于腔内辐射场描述的精确程度不同,需要分析不同建模对腔内辐射温度的影响。开展了三温建模与辐射多群输运建模下LARED集成程序数值模拟两孔球型黑腔模型,实现了球腔的完整数值模拟。数值模拟结果表明,三温与辐射多群输运模拟的等离子体状态接近,辐射温度存在差异。物理分析显示辐射温度差异的主要原因是使用的辐射不透明度,修改辐射不透明度参数后的三温计算结果与输运计算符合更好,从而可以用三温建模更快更准确地估计出所需的激光能量和功率。     

不同辐射建模下球腔的整体数值模拟

不同辐射建模对于腔内辐射场描述的精确程度不同,需要分析不同建模对腔内辐射温度的影响。开展了三温建模与辐射多群输运建模下LARED集成程序数值模拟两孔球型黑腔模型,实现了球腔的完整数值模拟。数值模拟结果表明,三温与辐射多群输运模拟的等离子体状态接近,辐射温度存在差异。物理分析显示辐射温度差异的主要原因是使用的辐射不透明度,修改辐射不透明度参数后的三温计算结果与输运计算符合更好,从而可以用三温建模更快更准确地估计出所需的激光能量和功率。     

任意四边形网格上解辐射输运方程的一种数值方法

解二维辐射输运问题时,对于矩形网格,用Sn菱形(diamond)格式处理十分简单有效,但该格式不能直接用于非矩形网格,而辐射流体力学计算往往要求在非矩形网格上进行.本文讨论任意四边形网格的辐射输运数值方法,给出一种全正的子网格平衡计算格式,并进行数值计算,与矩形网格Sn菱形格式结果比较,得出本方法可行性结论.     

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

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

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

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

辐射与物质强耦合问题的二维输运数值模拟

分析确定性输运数值方法模拟辐射与物质强耦合问题的难点,给出在任意四边形网格上求解该类问题的简单隅角平衡算法,并应用灰体输运综合加速算法提高输运方程的源迭代收敛速度,最后给出数值例子.     

多群辐射输运计算的输运综合加速方法

研究多群辐射输运计算的源迭代输运综合加速方法,通过优化灰体输运综合加速,构造了双层嵌套的输运综合加速方法.数值算例表明新方法较灰体输运综合加速进一步提高计算效率,且适用于吸收系数强间断的高维复杂几何问题.     

辐射输运与电子能量强耦合源项的一种高效计算方法

本文研究辐射输运和电子能量耦合方程组的数值方法. 在具有光性厚特征的应用问题中,这两类方程的耦合源项表现出强刚性,使得设计稳健高效的数值格式陷入困难.针对辐射输运多群模型和电子能量的耦合方程组的刚性源项,我们给出一种基于电子温度变化规律拟设（ansatz）的积分算法,其时间步长不受刚性源项限制,从而使得计算效率比传统显式方法或隐式非线性迭代获得本质的提高.在所基于的拟设有效时,算法确保给出具有物理意义的解,数值算例显示其给出的解具有较高的精确度.     

半透明梯度折射率介质内辐射熵传递方程及其数值模拟

在非相干辐射条件下,基于Planck光谱辐射熵强度定义,导出半透明梯度折射率介质内光谱辐射熵传递方程,以及局部辐射熵产率理论表达式.基于离散坐标法对辐射熵传递方程进行数值求解.以一维半透明梯度介质平板为例,对辐射熵方程及其算法进行验证.平板整体无因次辐射熵产的计算结果与宏观热力学定律的结果一致.     

High-Precision Direct Method for the Radiative Transfer Problems

It is the main aim of this paper to investigate the numerical methods of the radiative transfer equation.Using the five-point formula to approximate the differential part and the Simpson formula to substitute for integral part respectively, a new high-precision numerical scheme, which has 4-order local truncation error, is obtained. Subsequently,a numerical example for radiative transfer equation is carried out, and the calculation results show that the new numerical scheme is more accurate.     

辐射传输中的相函数展开方法研究

从辐射传输方程中相函数的展开方法出发,分别介绍了δ-M展开法和δ-fit展开法,并对这两种相函数展开方法进行了比较。结果表明,在相函数展开中,当所取勒让德多项式的阶数相同时,δ-M法的图像始终在真实相函数附近波动,始终无法逼近真实的相函数的值,而δ-fit法符合较好。     

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

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