Premium
Error bounds on block Gauss–Seidel solutions of coupled multiphysics problems
Author(s) -
Whiteley J. P.,
Gillow K.,
Tavener S. J.,
Walter A. C.
Publication year - 2011
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.3217
Subject(s) - gauss–seidel method , mathematics , block (permutation group theory) , multiphysics , linear algebra , linear system , discretization , gauss , partial differential equation , iterative method , mathematical optimization , finite element method , mathematical analysis , physics , geometry , quantum mechanics , thermodynamics
Mathematical models in many fields often consist of coupled sub‐models, each of which describes a different physical process. For many applications, the quantity of interest from these models may be written as a linear functional of the solution to the governing equations. Mature numerical solution techniques for the individual sub‐models often exist. Rather than derive a numerical solution technique for the full coupled model, it is therefore natural to investigate whether these techniques may be used by coupling in a block Gauss–Seidel fashion. In this study, we derive two a posteriori bounds for such linear functionals. These bounds may be used on each Gauss–Seidel iteration to estimate the error in the linear functional computed using the single physics solvers, without actually solving the full, coupled problem. We demonstrate the use of the bound first by using a model problem from linear algebra, and then a linear ordinary differential equation example. We then investigate the effectiveness of the bound using a non‐linear coupled fluid‐temperature problem. One of the bounds derived is very sharp for most linear functionals considered, allowing us to predict very accurately when to terminate our block Gauss–Seidel iteration. Copyright © 2011 John Wiley & Sons, Ltd.