一种求解三维热辐射输运方程的整体预处理迭代方法及并行计算

为三维灰体热辐射输运方程的隐式离散纵标方法发展一个整体预处理迭代方法并研制并行程序.该方法采用组装线性代数方程组策略,同时求出所有离散方向上的辐射强度.借助预处理的Krylov子空间迭代法,避免复杂网格上扫描算法可能遇到的死锁问题,能够提高健壮性和计算效率.空间离散上采用一阶迎风有限体积格式.数值实验测试变形六面体网格...

A Monolithic Preconditioned Iterative Solver and Parallel Computing for Three-dimensional Thermal Radiation Transport Equation

We present a monolithic preconditioned iterative solver for implicit discrete ordinate equations of three-dimensional grey thermal radiation transport and parallel codes are developed. A strategy of assembling linear algebraic systems is used to obtain radiation intensity in all discrete directions simultaneously. With preconditioned Krylov subspace iterative methods, the solver avoids possible mesh cycles in complex grids associated with sweep algorithms, which improves robustness and computational efficiency. First order upwind finite volume scheme is used for space discretization. Numerical experiments verify convergence rate on distorted hexahedral grids and assess performance of preconditioned iterative methods. Problems with coupled radiation and matter are simulated. Simulation results of three-dimensional crooked pipe and hohlraum problems are shown. It shows validity of the codes and flexibility of the method.