[petsc-dev] what's q in snes ex62
Matthew Knepley
knepley at gmail.com
Fri Nov 15 17:21:25 CST 2013
On Fri, Nov 15, 2013 at 5:09 PM, Geoffrey Irving <irving at naml.us> wrote:
> 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
>
u, v \in V and p, q \in P, and hopefully \div V = P for stability.
Matt
> Thanks,
> Geoffrey
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20131115/4f6e5dbd/attachment.html>
More information about the petsc-dev
mailing list