<!DOCTYPE html><html><head><title></title><style type="text/css">
p.MsoNormal,p.MsoNoSpacing{margin:0}</style></head><body><div>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.)<br></div><div><br></div><div><a href="https://libceed.org/en/latest/examples/fluids/">https://libceed.org/en/latest/examples/fluids/</a><br></div><div><a href="https://github.com/CEED/libCEED/blob/main/examples/fluids/qfunctions/stabilization.h#L76">https://github.com/CEED/libCEED/blob/main/examples/fluids/qfunctions/stabilization.h#L76</a><br></div><div><br></div><div><br></div><div>On Tue, Oct 10, 2023, at 7:34 PM, Brandon Denton via petsc-users wrote:<br></div><blockquote type="cite" id="qt" style=""><div style="font-family:Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif;font-size:12pt;color:rgb(0, 0, 0);" class="qt-elementToProof">Good Evening,<br></div><div style="font-family:Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif;font-size:12pt;color:rgb(0, 0, 0);" class="qt-elementToProof"><br></div><div style="font-family:Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif;font-size:12pt;color:rgb(0, 0, 0);" class="qt-elementToProof">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.<br></div><div style="font-family:Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif;font-size:12pt;color:rgb(0, 0, 0);" class="qt-elementToProof"><br></div><div style="font-family:Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif;font-size:12pt;color:rgb(0, 0, 0);" class="qt-elementToProof">Thank you in advance for your time.<br></div><div style="font-family:Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif;font-size:12pt;color:rgb(0, 0, 0);" class="qt-elementToProof">Brandon<br></div></blockquote><div><br></div></body></html>