多块结构网格上的Kershaw扩散格式

Kershaw格式是在四边形结构网格上求解扩散方程的-种经典格式.基于对Kershaw格式中"流"的深入理解,将其拓展到包含非结构点的多块结构网格,分别推导退化非结构点和强化非结构点情况下的Kershaw格式,拓展的Kershaw格式满足流连续条件.三个数值算例的计算结果与精确解吻合得很好,表明将Kershaw格式拓展到多块结构网格的正确性和有效性.     

二维三温能量方程的差分格式及平面靶数值模拟

通过统一二维直角坐标与柱坐标下三温能量方程的差分格式,实现了LARED-H(laser radiationelectron dynamic holehurm)程序既可以模拟正入射激光平面靶过程,也可模拟斜入射平面靶问题,从而扩展了LARED-H程序的功能.使用发展后的LARED-H程序进行了直角坐标系下线聚焦打靶耦合过程数值模拟,给出了不同入射角度下靶等离子体的膨胀过程.     

积分梯度法的两种格式

在二维柱坐标系下Lagrange流体力学的计算中,积分梯度法是动量方程的一种有效离散方法.积分梯度法中,IGT(Integral Gradient Total)格式不能保持柱几何下一维球对称性;IGA(Integral Gradient Average)格式可以保持一维球对称性,但当相邻网格质量相差比较大时,会得到远远脱离真实物理现象的加速度.深入研究IGA和IGT格式发现,当相邻网格边界压力取为质量加权时,即使相邻网格质量相差较大,对于一维平面和一维柱问题,IGT与IGA等价;在二维情形下,可以缩小IGT和IGA之间的差异.理论证明,IGA格式不能保持系统的动量守恒,IGT格式能保持系统的动量守恒性.数值模拟结果进一步显示了这两个格式的优缺点.     

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.     

空腔靶大变形问题的网格构造

针对ICF(Inertial Confinement Fusion)空腔靶大变形问题,发展了九点网格重构方法.拓展后的九点网格重构方法,保留了九点网格重构方法的基本特征(逼近椭圆型或抛物型方程,并具有一维球对称性),而且包含更多形式的九点网格重构方法,扩展了实际问题的使用需要.最后LARED-H程序采用ALE(ArbitraryLagrangian-Eulerian)方法,对空腔靶进行数值模拟,给出了合理的物理图像.     

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

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

基于神经网络的粒子输运问题高效计算方法

蒙特卡罗方法是求解粒子输运问题的有力工具之一,其局限性在于为达到精度要求需模拟大量粒子,计算耗时长,这阻碍了该方法的进一步应用,尤其在需快速响应的情形.本文结合神经网络和若干蒙特卡罗方法基本原理发展了一种计算方法,能够实现源分布可变,几何、材料和目标计数不变的中子输运问题的快速准确求解.首先,为高效生成用于神经网络训练的数据,利用重要性原理实现在同样模拟次数基础上有效扩充训练数据集容量,在一定程度上克服了使用蒙特卡罗计算获取训练数据耗时长的缺点.进而,基于目标计数是源分布与重要性函数乘积积分的事实,设计了利用神经网络实现快速输运计算的策略.该网络的输入是中子源项,输出是目标计数,在几何、材料和目标计数固定的情况下,该神经网络可重复使用,根据新的源项快速准确得到目标计数.本文所提出方法的原理和框架同样适用于其他种类粒子的同类型输运问题.基于若干基准模型的验证表明,训练得到的神经网络能在不到1 s的时间内得到目标计数,且与蒙特卡罗大样本模拟得到基准结果的平均相对偏差均低于5%.