<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" 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 Sijo,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Torque should simply compute the torque in each direction about the 3 axes emanating from the point prescribed by the coordinates of X0.</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%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Nek5000-users <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> Friday, June 22, 2018 6:20:05 PM<br>
<b>To:</b> nek5000-users<br>
<b>Subject:</b> [Nek5000-users] torque calculation on airfoil</font>
<div> </div>
</div>
<meta content="text/html; charset=utf-8">
<div>
<div style="font-family:arial,helvetica,sans-serif; font-size:12pt; color:#000000">
<div>Hello Neks,</div>
<div><br>
</div>
<div>Finally I was successfully able to import mesh from ICEM to Nek5000 with periodic BCs. Now I have to set up a fluid structure interaction code in NEK5000.</div>
<div><br>
</div>
<div>The test case is called self sustained pitching of an airfoil. For that I had to calculate torque about the Elastic Axis(EA). I know that using toque_calc function I can calculate torque in x and y direction.</div>
<div><br>
</div>
<div>But I am bit confused with the input parameters for this function.</div>
<div><br>
</div>
<div>Here I show the code,</div>
<div></div>
<div>
<div><br>
</div>
<div> real x0(3)</div>
<div> save x0</div>
<div> data x0 /3*0/</div>
<div> <br>
</div>
<div> </div>
<div> x0(1) = 0.29016 ! X- Direction (Chord wise direction where chord length = 1.56)</div>
<div> x0(2) = 0.0 ! Y -Direction</div>
<div> x0(3) = 0.2 ! Z- Direction. (Span wise direction where total span width = 0.4)</div>
<div><br>
</div>
<div> n = nx1*ny1*nz1*nelv</div>
<div></div>
<div><br>
</div>
<div> if (istep.eq.0) call set_obj ! define objects for surface integrals</div>
<div><br>
</div>
<div> scale = 2. ! Cd = F/(.5 rho U^2 ) = 2*F</div>
<div><br>
</div>
<div> call torque_calc(scale,x0,.true.,.true.)</div>
<div><br>
</div>
<div>I just have to confirm that,</div>
<div><br>
</div>
<div>1) Here x0 is the Elastic Axis point right?</div>
<div><br>
</div>
<div>2) Am I correct with the Axis notations? (for example: x0(1) represents X direction?)</div>
<div><br>
</div>
<div>3)In fact, the simulation is 3D, but the thickness is just 0.4 and there is no motion in Z direction. So, to calculate torque for pitching do I have to give the value for x0(3)? If yes is that okay if I give the midpoint in z direction?</div>
<div></div>
<div><br>
</div>
<div>Thanks in advance</div>
<div><br>
</div>
<div>Sijo GEORGE</div>
<div><br>
</div>
<div><br>
</div>
<div></div>
</div>
</div>
</div>
</body>
</html>