<div dir="ltr"><div>Hi Paul and Phillip,<br><br></div><div>I have validated my model at constant AOA and am now returning to my pitch and hold scenario to do some work. I definitely think the ALE formulation applied to the entire domain is the way to go, but I'm feeling a little confused about the "mv" boundary condition. In the mv_wall example, the "mv" boundary condition is applied to the section of the wall that is moving, so in my mind it follows logically that to have the entire domain pitch about the z axis, the entire domain would need the "mv" condition. This doesn't seem possible while preserving my inlet and outlet BCs. I understand that for ALE to work appropriately, we need to match the normal component of fluid velocity to that of the mesh velocity at any boundary. Does this apply also to the inlet and outlet, or only internal fluid boundaries such as my airfoil?<br></div><div><br></div><div>I am also looking at the mv_cyl and ocyl/ocyl2 examples, but these are in rea/re2 format and I cannot look at the element boundary conditions (forgive me, I'm still learning terminology so some of what I say is not 100% accurate). While looking through the example files, I see that the logical switch "ifusermv" is assigned in "usrchk" in some of the examples, and in one of the "usrdat" subroutines in others. I don't really appreciate the difference, is there one?<br></div><div><br></div><div>I am building my mesh in ANSYS Fluent, and then converting the mesh to ASCII .rea format with mshconvert2, and then extruding to 3D and converting to rea/re2 format with n2to3. Before, I would scroll to the bottom of the .msh file and add Nek5000 BC flags to each boundary (i.e. adding an _W to wall, _O to outlet, and _v to inlet) before converting to .rea format, and this would successfully apply boundary conditions. However, when attempting to insert the _mv boundary condition on my wall, the boundary conditions in the .rea file are unchanged. Is there a better/more appropriate way to assign the "mv" BC when converting a mesh?<br></div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Aug 21, 2018 at 11:30 AM <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div id="m_6723969809221433119m_-1656579852477447871divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Hi Harry,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">One simple approach is to use ALE but to have the entire domain pitch,</p>
<p style="margin-top:0;margin-bottom:0">thus preserving your mesh structure.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">That would properly account for all the accelerations in your non inertial</p>
<p style="margin-top:0;margin-bottom:0">domain. </p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">The slightly tricky part would be to determine what BCs you need. You</p>
<p style="margin-top:0;margin-bottom:0">need to be careful about changing the character of the BCs (e.g., from inflow</p>
<p style="margin-top:0;margin-bottom:0">to outflow) in the middle of a computation because that changes which quantities</p>
<p style="margin-top:0;margin-bottom:0">are prescribed and which are unknown. Perhaps one approach would be</p>
<p style="margin-top:0;margin-bottom:0">to have a circular domain with 3/4 of it Dirichlet and 1/4 being outflow (i.e.,</p>
<p style="margin-top:0;margin-bottom:0">Neumann for velocity).</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">hth,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Paul</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_6723969809221433119m_-1656579852477447871divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Nek5000-users <<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.mcs.anl.gov</a>> on behalf of <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>><br>
<b>Sent:</b> Tuesday, August 21, 2018 9:43:56 AM<br>
<b>To:</b> <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a><br>
<b>Subject:</b> [Nek5000-users] Rotating computational domain</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div>Hello Neks,<br>
<br>
</div>
<div>I am an undergraduate student new to Nek and I am working on a simulation studying dynamic stall in a dynamically pitching airfoil. I want to do a pitch and hold simulation at low Reynold's number (Re = 12000) where I rotate the entire computational domain
up to about 50 degrees AOA and then hold at 50 degrees, with steady inflow ux = 1. I would like to rotate the domain about the z axis, which passes through the quarter chord location of the airfoil profile (modified NACA 0012).<br>
</div>
<div><br>
</div>
<div>Because of the steep AOA, I would like to stay away from pitching the airfoil within the stationary domain with ALE to avoid too much mesh deformation. To this end, I have been scouring documentation to try and find a method of applying a rotational mesh
velocity to the entire domain about the z axis while keeping the flow within the domain at ux = 1. I am feeling very stuck. Others have suggested applying a Coriolis force to the flow, but I don't really see how this would simulate rotation of the entire
domain. It seems to me that simulation of a Coriolis force would be acceptable to correct the flow field to always be ux = 1 inside a rotating domain, but I don't see how that would apply a mesh velocity anywhere.</div>
<div><br>
</div>
<div>Are there any examples where an entire computational domain is rotated that I can look at, or has anybody done similar work who can point me in the right direction?</div>
<div><br>
</div>
<div>Regards,</div>
<div>-Harry</div>
</div>
</div>
</div>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="m_6723969809221433119gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Harry E. Werner IV<div>Secretary</div><div>Clarkson University SCUBA Club</div><div><u>Cell:</u> (716) 570-4023</div></div></div>