[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