Parallel Newton-Krylov methods for the solution of
unstructured finite element reacting flow applications
John N. Shadid, Raymond S. Tuminaro, Homer F. Walker
Sandia National Laboratories
Albuquerque, New Mexico 87185
Abstract
The solution of the governing steady transport equations for
momentum, heat and mass transfer in flowing fluids can be very
difficult. These difficulties arise from the nonlinear, coupled,
nonsymmetric nature of the system of algebraic equations that results
from spatial discretization of the PDEs. In this talk we focus on
evaluating a proposed nonlinear and linear solution method based
on an inexact Newton method with backtracking and on preconditioned
Krylov methods. In this context we use a particular
spatial discretization based on a pressure stabilized Petrov-Galerkin
finite element formulation of the low Mach number Navier-Stokes
equations with heat and mass transport. Our discussion considers
computational efficiency, robustness and some implementational issues
related to the proposed nonlinear solution scheme. Computational results
are presented for several challenging CFD benchmark problems as well
as a large scale 3D flow simulation.