On Tue, Oct 10, 2023 at 9:34 PM Brandon Denton via petsc-users <petsc-users@mcs.anl.gov> wrote:
Good Evening,

I am looking to implement a form of Navier-Stokes with SUPG Stabilization and shock capturing using PETSc's FEM infrastructure. In this implementation, I need access to the cell's shape function gradients and natural coordinate gradients for calculations within the point-wise residual calculations. How do I get these quantities at the quadrature points? The signatures for fo and f1 don't seem to contain this information.

Are you sure you need those? Darsh and I implemented SUPG without that. You would need local second derivative information, which you can get using -dm_ds_jet_degree 2. If you check in an example, I can go over it.

  Thanks,

     Matt
 
Thank you in advance for your time.
Brandon


--
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