Newton-Krylov methods applied to nonequilibrium radiation diffusion
The authors present results of applying a matrix-free Newton-Krylov method to a nonequilibrium radiation diffusion problem. Here, there is no use of operator splitting, and Newton`s method is used to convert the nonlinearities within a time step. Since the nonlinear residual is formed, it is used to monitor convergence. It is demonstrated that a simple Picard-based linearization produces a sufficient preconditioning matrix for the Krylov method, thus elevating the need to form or store a Jacobian matrix for Newton`s method. They discuss the possibility that the Newton-Krylov approach may allow larger time steps, without loss of accuracy, as compared to an operator split approach where nonlinearities are not converged within a time step.
- Research Organization:
- Los Alamos National Lab., Applied Theoretical and Computational Physics Div., NM (United States)
- Sponsoring Organization:
- USDOE Assistant Secretary for Management and Administration, Washington, DC (United States)
- DOE Contract Number:
- W-7405-ENG-36
- OSTI ID:
- 666037
- Report Number(s):
- LA-UR--98-1076; CONF-980393--; ON: DE98006340
- Country of Publication:
- United States
- Language:
- English
Similar Records
Physics-based preconditioning and the Newton-Krylov method for non-equilibrium radiation diffusion
A multigrid preconditioned Newton-Krylov method