My initial plan was to write a new code using only PETSc. However, I don't see how to do what I want within the point-wise residual function. Am I missing something?
Yes. I would be interested in collaborating on the ceed-fluids. I took a quick look at the links you provided and it looks interesting. I'll warn you though. I'm a Mechanical Engineer by trade/training. The calculus and programming sometimes take me a little
while to wrap my head around. Let me know how I can help. In the meantime, I'll continue to review the information you sent over.
Do you want to write a new code using only PETSc or would you be up for collaborating on ceed-fluids, which is a high-performance compressible SUPG solver based on DMPlex with good GPU support? It uses the metric to compute covariant length for stabilization.
We have YZƁ shock capturing, though it hasn't been tested much beyond shock tube experiments. (Most of our work has been low Mach.)
On Tue, Oct 10, 2023, at 7:34 PM, Brandon Denton via petsc-users 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.
Thank you in advance for your time.
Brandon