二维三温能量方程的部分Newton迭代法

符尚武等对二维三温能量方程提出了一种9点差分格式，它适用于任意三维网格。但他们用非线性块Gauss—Seidel方法求解所得的非线性代数方程组，收敛得非常慢且经常不得不因为迭代某些次数后仍不收敛而缩小时间步长。     

二维三温能量方程组离散求解的两个新预处理技术

二维三温能量方程离散后得到的稀疏线性代数方程组中,系数矩阵各行的对角占优性相差十分悬殊,矩阵元素相差也十分大.针对前一问题,提出了改善对角占优性的一个新比例化方法.针对后一问题,利用每次舍弃前计算多个行的技术提出了多行ILUT预条件方法.最后,将对角占优性改善技术、多行ILUT与对角元比例化技术、RCM排序联合使用于实际的能量方程离散求解中,取得了较好的加速效果.     

求解二维三温辐射扩散方程组的一种代数两层迭代方法

在二维三温辐射扩散方程离散代数方程组的求解中,由于光子、电子和离子温度之间存在耦合关系,而且三个温度在同种介质中有不同的扩散性质,使得经典的代数多重网格(AMG)方法难以直接应用.基于特殊粗化策略,在粗网格层解除了这种耦合关系,得到一种代数两层网格方法,而粗网格方程由经典AMG方法求解.将这一算法具体应用于JFNK(Jacobian自由的Newton-Krylov)框架中预处理方程的求解,并基于该框架求解二维三温辐射扩散方程组.数值结果显示了算法的可扩展性和健壮性.     

二维三温辐射扩散方程组两层预条件子的自适应求解

针对实际应用中若干典型三温线性系统,分析求解二维三温辐射扩散方程离散线性系统的代数两层预条件子(PCTL)的算法效率.结果表明,PCTL的算法效率与三个温度之间的耦合强度以及单温子系统对角占优性强弱程度有很大关系.为此,通过刻画三温线性系统的耦合强度和单温子系统对角占优性特征,提出一种PCTL中子系统的自适应求解算法.数值结果表明,可以显著改善PCTL的算法效率.对于实际数值模拟应用中37个典型三温线性系统,相对于经典AMG算法,算法整体加速2.5倍.数值实验表明算法具有很强的鲁棒性.