针对实际应用中若干典型三温线性系统,分析求解二维三温辐射扩散方程离散线性系统的代数两层预条件子(PCTL)的算法效率.结果表明,PCTL的算法效率与三个温度之间的耦合强度以及单温子系统对角占优性强弱程度有很大关系.为此,通过刻画三温线性系统的耦合强度和单温子系统对角占优性特征,提出一种PCTL中子系统的自适应求解算法.数值结果表明,可以显著改善PCTL的算法效率.对于实际数值模拟应用中37个典型三温线性系统,相对于经典AMG算法,算法整体加速2.5倍.数值实验表明算法具有很强的鲁棒性.
Performance of physical-variable based coarsening two-level (PCTL) preconditioner is analyzed for typical linear systems discretizated from two-dimensional (2-D) radiative diffusion equations with photon, electron, ion temperatures (3-T). It reveals that performance of PCTL strongly depends on both coupling of three temperatures and diagonally dominance of three diagonal sub-matrices of coefficient matrix. An adaptive algorithm for sub linear systems in PCTL is proposed. Numerical results show efficiency and robustness of the method. For 37 2-D 3-T linear systems in simulations, PCTL based on the algorithm speeds up 2.5 times compared with classical algebraic muhigrid (AMG) preconditioners.