[petsc-users] Continuous FE with Discontinuous Pressure Field
Matthew Knepley
knepley at gmail.com
Thu Mar 19 16:41:00 CDT 2020
On Thu, Mar 19, 2020 at 5:05 PM Jeremy L Thompson <jeth8984 at colorado.edu>
wrote:
> Good Afternoon,
>
> I have a working high-order 3D elasticity code built on PETSc + libCEED.
> My DM has a three component displacement FE field. I'd like to add a
> discontinuous pressure field (P_-1).
>
By P_{-1} do you mean a discontinuous linear basis? I only ask because FV
is normally P_0.
Anyway, PetscFE has discontinuous bases if you want to use that. You use
the flag
https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DUALSPACE/PetscDualSpaceLagrangeSetContinuity.html
Thanks,
Matt
> Am I missing an obvious way to set up a discontinuous FE field? Or is a
> FV field a more natural choice for this sort of thing? Or a different
> approach?
>
> I am managing action of the FDE operator myself via libCEED, so I don't
> need to get too into the weeds about anything but the DM. I just need to
> find the best way to manage the discontinuous pressure values on each
> element in my DM.
>
> Thanks,
>
> Jeremy L Thompson
>
>
--
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
https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200319/76f40681/attachment.html>
More information about the petsc-users
mailing list