<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
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?</div>
<div style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="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="elementToProof">
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.</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Jed Brown <jed@jedbrown.org><br>
<b>Sent:</b> Tuesday, October 10, 2023 10:18 PM<br>
<b>To:</b> Brandon Denton <bldenton@buffalo.edu>; petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
<b>Subject:</b> Re: [petsc-users] FEM Implementation of NS with SUPG Stabilization</font>
<div> </div>
</div>
<style type="text/css">
<!--
p.x_MsoNormal, p.x_MsoNoSpacing
{margin:0}
-->
</style>
<div>
<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/" originalsrc="https://libceed.org/en/latest/examples/fluids/" shash="Y5bMoEgFjVAidRE+oF0YK4L+CTjHA3ZaUL2tjwlmABAHsZRY4TEDEgSiGLBgIQQEGForJt13+cclPlU1rpHTEXihVgQP0QHUkqrH9IYASrPKScP65qAVsMtvD3nvJOnQ/MAmgUhGCalUunnUob+2rvIHmAVamGXmq78VR9wX9bw=">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" originalsrc="https://github.com/CEED/libCEED/blob/main/examples/fluids/qfunctions/stabilization.h#L76" shash="Ye7x6TGQpg/7p5vwTdllUSaf4uG8LxKin7Pe1OlWRirNeiN67Ov8YbFkZhMAsfTKmp6t7eo4IdgNlHjhLpcbIsVkDGCO3rELRLxW8nAZlNPlg8xq9Tz9SyOhd5qcPbPWjMiZ+IsofhP6o4rLR49F3lf+5MPPsPPZUalw2Vz/nNY=">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="x_qt" style="">
<div class="x_qt-elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Good Evening,<br>
</div>
<div class="x_qt-elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div class="x_qt-elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
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 class="x_qt-elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div class="x_qt-elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Thank you in advance for your time.<br>
</div>
<div class="x_qt-elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Brandon<br>
</div>
</blockquote>
<div><br>
</div>
</div>
</body>
</html>