<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;">Hi Rif,
<div><br>
</div>
<div>You want to use ux, uy and uz (if 3D).</div>
<div><br>
</div>
<div>For example, if Omega = (0 , 0 , omega), you would set</div>
<div><br>
</div>
<div> omega = whatever</div>
<div><br>
</div>
<div> ffx = 2*uy*omega</div>
<div> ffy = -2*ux*omega</div>
<div> ffz = 0</div>
<div><br>
</div>
<div>in userf.</div>
<div><br>
</div>
<div>ux,uy,uz are user interface values. In userbc and useric, they are set by the user.</div>
<div>They may also be _read_ by the user, since they are assigned the vx(i,j,k,e), vy(i,j,k,e),</div>
<div>and vz(i,j,k,e) values prior to the call, point by point. The same holds for userq and </div>
<div>userf --- values are set prior and read after. If the routine/call under consideration is</div>
<div>looking for values of vx (say, on an inflow boundary in the case of userbc), then the user</div>
<div>sets ux and this value is loaded into vx by Nek.</div>
<div><br>
</div>
<div>hth,</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF451359" 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, February 23, 2016 2:28 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Can nek do this?<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div>
<div>Hi developers and users of nek5000,<br>
</div>
I have been trying for a long time now to rotate the whole computational domain to simulate a rotating wing. I have scoured the nek forums but was unable to get any information on this. Is this even possible? I initially wanted to add a Coriolis force but coriolis
in the Navier stokes is a cross product between angular velocity vector and velocity vector and in the userf, when tried this I received a compilation error saying that the ranks are different. This implied that i must use ux, uy and uz instead of vx,vy
and vz. But as i understand, ux uy and uz are given values and not the velocity field taken from the calculations. Also for a rotating wing, it don't think it is advisable to apply a coriolis force to simulate rotation as the geometry is fixed and the angle
of attack is not changing. This is why I'm trying to rotate the whole computational domain. In that sense, i can avoid using the ALE formulation for deforming meshes and therefore avoid the vanishing Jacobian errors which I have encountered. Can anybody guide
me in the right direction?<br>
<br>
<br>
</div>
Rif <br>
</div>
</div>
</div>
</div>
</div>
</body>
</html>