<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">
Thank you for the discussion.</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">
Are we agreed then that the derivatives of the natural coordinates are required for the described approach? If so, is this something PETSc can currently do within the point-wise residual functions?</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">
Matt - Thank you for the command line option for the 2<span><sup>nd</sup> derivatives. Those will be needed to implement the discussed approach. Specifically in the stabilization and shock capture parameters. (Ref.: B. Kirk's Thesis). What is a good reference
for the usual SUPG method you are referencing? I've been looking through my textbooks but haven't found a good reference.</span></div>
<div style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span><br>
</span></div>
<div style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span>Jed - Thank you for the link. I will review the information on it. </span></div>
<div style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span><br class="Apple-interchange-newline ContentPasted0">
<span style="display: inline !important; background-color: rgb(255, 255, 255);" class="ContentPasted0">Sorry about the attachment. I will upload it to this thread later (I'm at work right now and I can't do it from here).</span><br>
</span></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> Wednesday, October 11, 2023 1:38 PM<br>
<b>To:</b> Matthew Knepley <knepley@gmail.com><br>
<b>Cc:</b> Brandon Denton <bldenton@buffalo.edu>; petsc-users <petsc-users@mcs.anl.gov><br>
<b>Subject:</b> Re: [petsc-users] FEM Implementation of NS with SUPG Stabilization</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">Matthew Knepley <knepley@gmail.com> writes:<br>
<br>
> On Wed, Oct 11, 2023 at 1:03 PM Jed Brown <jed@jedbrown.org> wrote:<br>
><br>
>> I don't see an attachment, but his thesis used conservative variables and<br>
>> defined an effective length scale in a way that seemed to assume constant<br>
>> shape function gradients. I'm not aware of systematic literature comparing<br>
>> the covariant and contravariant length measures on anisotropic meshes, but<br>
>> I believe most people working in the Shakib/Hughes approach use the<br>
>> covariant measure. Our docs have a brief discussion of this choice.<br>
>><br>
>> <a href="https://libceed.org/en/latest/examples/fluids/#equation-eq-peclet">https://nam12.safelinks.protection.outlook.com/?url=https%3A%2F%2Flibceed.org%2Fen%2Flatest%2Fexamples%2Ffluids%2F%23equation-eq-peclet&data=05%7C01%7Cbldenton%40buffalo.edu%7Cd9372f934b26455371a708dbca80dc8e%7C96464a8af8ed40b199e25f6b50a20250%7C0%7C0%7C638326427028053956%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=skMsKDmpBxiaXtBSqhsyckvVpTOkGqDsNJIYo22Ywps%3D&reserved=0</a><br>
>><br>
>> Matt, I don't understand how the second derivative comes into play as a<br>
>> length measure on anistropic meshes -- the second derivatives can be<br>
>> uniformly zero and yet you still need a length measure.<br>
>><br>
><br>
> I was talking about the usual SUPG where we just penalize the true residual.<br>
<br>
I think you're focused on computing the strong diffusive flux (which can be done using second derivatives or by a projection; the latter produces somewhat better results). But you still need a length scale and that's most naturally computed using the derivative
of reference coordinates with respect to physical (or equivalently, the associated metric tensor).<br>
</div>
</span></font></div>
</body>
</html>