Implicit solution of large-scale radiation - material energy transfer problems
- LLNL
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 paper, we present an algorithm for the implicit solution of the multi- group 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 take advantage of the fast, robust nonlinear convergence of Newton's method and the scalability of the linear multigrid method. Numerical results show that the method is accurate and scalable.
- Research Organization:
- Lawrence Livermore National Laboratory, Livermore, CA
- Sponsoring Organization:
- USDOE Office of Defense Programs (DP)
- DOE Contract Number:
- W-7405-ENG-48
- OSTI ID:
- 3416
- Report Number(s):
- UCRL-JC-132831; DP0101031; ON: DE00003416
- Country of Publication:
- United States
- Language:
- English
Similar Records
Physics-based preconditioning and the Newton-Krylov method for non-equilibrium radiation diffusion
A multigrid preconditioner and automatic differentiation for non-equilibrium radiation diffusion problems