Modeling of radiation-diffusion processes has traditionally been accomplished through simulations based on decoupling and linearizing the basic physics equations. By applying these techniques, physicists have simplified their model enough that problems of moderate sizes could be solved. However, new applications demand the simulation of larger problems for which the inaccuracies and nonscalability of current algorithms prevent solution. Recent work in iterative methods has provided computational scientists with new tools for solving these problems.

In this talk, we present an algorithm for the implicit solution of the diffusion approximation coupled to an electron temperature equation. This algorithm uses a stiff ODE solver coupled with Newton's method for solving the implicit equations arising at each time step. The Jacobian systems are solved by applying GMRES preconditioned with a semicoarsening multigrid algorithm. By combining the nonlinear Newton iteration with a multigrid preconditioner, we hope to take advantage of the fast, robust nonlinear convergence of Newton's method and the scalability of the linear multigrid method. Numerical results will be shown indicating that the method is accurate. Furthermore, both algorithm and parallel scalabilities will be explored.

MINISIMPOSIUM SESSION: 9. Solving linear and nonlinear systems in differential equations