[petsc-users] Using DMPlex for Darcys law

[Justin Chang]

Hi,

So in ex62 of SNES it solves the stokes equation with P2-P1 elements. I
want to solve the darcy equation with P2-P1 elements in a similar fashion.

The difference between stokes and darcy is that in the weak form of darcy,
the velocity-velocity term is v*u whereas in stokes it's grad[v]*grad[u]

In ex62, the function g3_uu (lines 172-183) formulates the grad[v]*grad[u]
term but how would I do it for Darcy's equation where its simply v*u
(assuming drag coefficient is scalar quantity of 1.0)?

Thanks,
Justin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140806/7470be12/attachment.html>