[petsc-users] Continuous FE with Discontinuous Pressure Field

[Jeremy L Thompson]

Thanks Matt, that's what I need.

Jeremy

On 3/19/20 3:41 PM, Matthew Knepley wrote:
> On Thu, Mar 19, 2020 at 5:05 PM Jeremy L Thompson 
> <jeth8984 at colorado.edu <mailto: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/080e34b7/attachment.html>