<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div><br>
</div>
Hi Andrew,
<div><br>
</div>
<div>I checked into this... I think what is happening is that the 'O ' bc for the stress formulation</div>
<div>means stress-free, which is not guaranteed to yield a parabolic profile at the outlet.</div>
<div><br>
</div>
<div>If you increase your Reynolds number I'm guessing that you'll recover the parabolic</div>
<div>profile because the viscous stresses will diminish -- this is what I observed.</div>
<div><br>
</div>
<div>From my perspective, the outlet boundary is not a region where I would count on accurate</div>
<div>physics --- it is, after all, a truncated domain, so I don't generally worry too much about</div>
<div>the behavior there.</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF241228" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> nek5000-users-bounces@lists.mcs.anl.gov [nek5000-users-bounces@lists.mcs.anl.gov] on behalf of nek5000-users@lists.mcs.anl.gov [nek5000-users@lists.mcs.anl.gov]<br>
<b>Sent:</b> Tuesday, December 01, 2015 2:44 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Stress formulation<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div><span lang="en"><span title="Andrew"></span></span>Hi Neks,<br>
<br>
I have one question about stress formulation in NEK5000.<br>
I tried to make a simulation for Poiseuil flow with stress formulation. In inlet boundary condition use was made of parabolic profile. If IFSTRS=false the outlet profile is exactly the same as inlet one but in case of IFSTRS=true the outlet profile changes
significantly (both the maximum value and profile shape). Could you tell me how I can get the parabolic profile at outlet, please?<br>
The same thing occurs with pressure values. The maximum pressure value with IFSTRS=true is three times higher then the one with IFSTRS=false.<br>
<br>
I found that the difference comes from subroutines where the stiffness matrix is calculated (axehlm for IFSTRS = false and axhmsf for IFSTRS = true) but I don't undestand what is happening exactly.<br>
<br>
Thanks in advance for your help,<br>
Best regards,<br>
Andrew<br>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>