<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Hi Neks, <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am trying to run a simulation case similar to the turbulent channel flow case. I initially ran a simulation on the normal turbulent channel flow domain (with the geometry like the fig on the left). I used the parameter p55 in the rea
file to ensure that the bulk velocity of the flow is 1 in the domain, and used the parameter p02 to set the
<b>turbulent</b> Reynolds number, Re_tau to my required value of 180. I used the sub-routine drag_calc to find the mean shear on the wall, from which I calculated the shear stress and then the Re_tau and u_tau values. However, recently I was trying the same
thing on a slightly different domain (fig to the right). <o:p></o:p></p>
<p class="MsoNormal">--------------------------------------------------------------------- ---------------------------------------------------------------------<o:p></o:p></p>
<p class="MsoNormal">| | /
/
<sub>^</sub> Y<o:p></o:p></p>
<p class="MsoNormal">| | /
/ |<o:p></o:p></p>
<p class="MsoNormal">| | /
/ |_____<sub>></sub> X<o:p></o:p></p>
<p class="MsoNormal">| | /
/<o:p></o:p></p>
<p class="MsoNormal">| | /
/<o:p></o:p></p>
<p class="MsoNormal">--------------------------------------------------------------------- ---------------------------------------------------------------------<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">This domain is almost the same as the first except now there is an angle of 30 degrees in the mesh grid about the z-axis. Please note that the perpendicular distance between the walls normal to the y-direction
<b><i>remains the same</i></b> for both the domains. In this domain if I use the same p55 and p02 parameters that I use for the orthogonal domain, then Re_tau, calculated as described above, is increasing to higher values (~220). If I reduce the parameter p55
to lesser values, then the Re_tau value is restored to 180. I could not figure out why exactly this was happening, and whether this difference is now because my way of calculating the Re_tau is not valid anymore in the angled domain. Or is the change because
the mass flow rate is actually higher in the angled domain? If the case is latter, then why is the mass flow rate increasing just due to the angling of the mesh?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I went through the sub-routines drag_calc, vol_flow and set_obj (to see if the integral object needs to be changed) in moderate detail, but could not figure out the reasons for this. Any help in this might be appreciated, as I want to establish
the reason, so that I can trust the accuracy of the results of this simulation. <o:p>
</o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks and regards,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Pradyumna M K Prasad<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">TU Delft,<o:p></o:p></p>
<p class="MsoNormal">Student: Masters in Aerospace Engineering<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<br><br></body>
</html>