In order to make conventional implicit algorithm to be applicable in large scale parallel computers,an interface correction method is introduced to solve the particle transport S_N equations.Domain decomposition is adopted in the computational domain.On the interface,an upwind explicit scheme is applied to give an incident boundary condition,which enables the subdomain problem iterated independently.The interface values are updated by an implicit scheme concurrently in iteration.The scheme shows good precision,parallelism and simplicity in numerical experiments.

We developed a high performance algebraic solver for nonlinear systems discretized from two-dimensional energy equations with three temperatures by a nine point scheme.The main idea is to solve the system by an inexact Newton method and preconditioned Krylov subspace methods in the frame of PNK and JFNK methods.Numerical experiments show the efficiency of the algebraic solvers.It is shown that our PNK method is 6 times faster than the nonlinear block Gauss-Seidel method. The JFNK and PNK methods are also compared.

A new type of hybrid iterative method is presented, which is competent in solving the large scale sparse linear systems derived from the 2-dimensional 3-temperature radiation dynamic energy equations. Numerical results show that the new method is as 4 times fast as the old ones.Especially on the cases the old one doesn't converge, the new method can easily get the solution to the precision required. It can successfully complete the simulation and the final physical parameters of the simulation coincide with the theory.