Premium
Stabilisation of AMG solvers for saddle‐point stokes problems
Author(s) -
Webster Ronald
Publication year - 2015
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.4199
Subject(s) - scalar (mathematics) , mathematics , saddle point , grid , diagonal , block (permutation group theory) , smoothing , saddle , certificate , block matrix , factorization , mathematical optimization , algorithm , geometry , physics , statistics , eigenvalues and eigenvectors , quantum mechanics
Summary When a block factorisation is used to precondition the saddle‐point equations of the discrete Stokes problem, the stability that this gives for the relaxation of residual errors may not be conserved in the coarse‐grid approximations (CGA) of algebraic multi‐grid (AMG) solvers. If the same first‐order interpolation is used in the inter‐grid transfer operators for the scalar and the vector fields, the conditioning degrades with each coarsening step until eventually a critical coarsening is reached beyond which residual errors are no longer damped and will become divergent with any further coarsening. It is shown that by introducing the same block pre‐conditioner as an integral part of the coarsening algorithm, stable smoothing can be maintained at all levels of the CGA. The pre‐conditioning need only be applied at preselected grid levels, one immediately before the critical threshold and others beyond that level if required. Excessive complexity in the CGA is thereby avoided. The method is purely algebraic and may be used for both classical AMG solvers and for smoothed‐aggregation AMG solvers. It should be applicable to other coupled vector and scalar fields in science and engineering that involve second‐order (block‐diagonal) and first‐order (block‐off‐diagonal) discrete difference operators. Copyright © 2015 John Wiley & Sons, Ltd.