Is the second equation the weak form of incompressibility? If so, why
isn't q u? It's used as q several places throughout the code, so it
seems more likely to be my misunderstanding than a typo.
< \nabla v, \nabla u + {\nabla u}^T > - < \nabla\cdot v, p > + < v, f > = 0
< q, \nabla\cdot v > = 0
Thanks,
Geoffrey