Hi Dave, My understanding is that stabilization by pressure weighted interpolation in FVM is in essence a fourth order pressure derivative with proper scaling. But in FVM it is usually implemented in defect correction form: only part of the stabilization is in the A11 block of the coupled matrix, the rest is in de rhs vector. This makes the coupled system diagonally dominant and solvable by AMG. Early work was done at Waterloo in the late eighties and early nineties, see http://dx.doi.org/10.2514/6.1996-297
Chris dr. ir. Christiaan Klaij CFD Researcher Research & Development E mailto:c.kl...@marin.nl T +31 317 49 33 44 MARIN 2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands T +31 317 49 39 11, F +31 317 49 32 45, I www.marin.nl
