论文部分内容阅读
二维三温辐射热传导方程是激光约束核聚变（Laser-Driven Inertial Confinement Fusion，ICF）数值仿真研究中的一个简化模型，它近似地描述了ICF内爆动力学过程中辐射能量在静止介质中的非线性传播过程。本文从空间离散方法、预条件子构造和网格自适应等几个方面开展了二维三温辐射热传导方程的自适应高效算法研究。论文包含了五个方面的工作。首先，论文在三角形网格的基础上，设计了一类保对称的有限体元方法（SFVEM），并从理论上证明了在一定条件下的二维三温辐射热传导方程SFVEM离散系统的对称性。其次，论文在有限元方法（FEM）预条件理论的基础上，通过分析两类有限体元方法（经典有限体元方法（FVEM）及SFVEM）的系数矩阵和FEM系数矩阵的关系，提出了利用线性有限元的预条件子预处理线性有限体元离散系统的新思路，并从理论上证明了此类预条件子的高效性。第三，论文针对二维三温辐射热传导方程时间半离散后的非线性偏微分方程，设计了一种高效的可并行的两网格离散算法，并给出了该方法在||·||1下的误差分析。第四，本文从基于Hessian矩阵的网格自适应理论出发，设计了一种适用于二维三温辐射热传导方程的网格自适应算法。最后，本文以知名数值仿真软件包为蓝本，将常用的偏微分方程数值方法与论文中新的研究结果相结合，制作了一套专门用于求解二维三温辐射热传导方程的程序包RH3T-2D。论文在RH3T-2D程序包的基础上进行了大量的数值实验，数值实验的结果表明：（1）、SFVEM方法具有比FVEM、FEM更快的离散速度和更小的能量守恒误差；（2）、在使用有限体元方法求解线性椭圆方程和二维三温辐射热传导方程时，随着网格规模的增大，以AMG为预条件的GMRES方法和PCG方法，在求解线性代数系统时，其迭代次数趋于稳定；（3）、两网格离散方法在||·||1下的误差阶得到了验证，相关算法适用于二维三温辐射热传导方程；（4）、基于Hessian矩阵的网格自适应方法能准确追踪二维三温辐射热传导方程的温度突变区域，提高系统的计算速度，相对于基于梯度的网格自适应方法，本文所设计的算法能够更加有效地降低能量守恒误差。