<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>
<div>
<div>Hi Ravi,</div>
<div><br>
</div>
<div>Inspection of the source routine navier5.f will reveal most of the answers to</div>
<div>your questions. That is where drag and torque are calculated. (I recommend</div>
<div>torque_calc, which gives the drag as well and is somewhat superior to drag_calc.)</div>
<div><br>
</div>
<div>hcode = history code --- there are several codes used (the convention</div>
<div>follows the Nekton 2.0 convention established in the 80s). Of relevance to</div>
<div>drag/lift/torque are hcode(10,.), which must be 'I' for Integral</div>
<div>quantity, and 1,2,3, which must be 'F' for force. Just follow the</div>
<div>convention in the ext_cyl.usr example.</div>
<div><br>
</div>
<div>Note that the 2nd argument of hcode indicates how many objects</div>
<div>you have. An Object is a collection of faces. So, you might have</div>
<div>more than one cylinder on which to compute drag, in which case you</div>
<div>would have more than one object.</div>
<div><br>
</div>
<div>Each object comprises "nmember(iobj)" faces.</div>
<div><br>
</div>
<div>nhis is the number of history elements you are tracking.</div>
<div>It is normally zero, plus the number of objects. It is</div>
<div>a measure of the number of lines of information dumped to</div>
<div>the .sch file per timestep.</div>
<div><br>
</div>
<div>cbc is the character array holding the boundary condition</div>
<div>specifier - one per face per field (velocity, temperature, etc.)</div>
<div><br>
</div>
<div>object is a pointer to the members of an object. A member of</div>
<div>an object is a given face (1-6) of a given element. The</div>
<div>element id is global -- (i.e., not the local element number</div>
<div>which is specific to a given mpi rank).</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
</div>
</div>
<div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF785415" 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> Monday, May 26, 2014 4:39 PM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Object definition for surface integrals<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">Dean NEK Users,<br clear="all">
<div><br>
</div>
<div>I am looking for some sort of help or documentation explaining the object definitions for evaluating surface integrals.
<br>
<br>
</div>
<div>My eventual application will be to evaluate the drag and lift forces for the case of a 2D flow past a cylinder, supported on elastic foundation, i.e., vortex-induced vibration simulations using NEK5000.<br>
<br>
</div>
<div>So far, the examples available in the repository and the NEK mailing list have enabled me to simulate the motionless cylinder and the case of cylinder moving with a prescribed motion.
<br>
<br>
However, I have hard time decoding the surface integral evaluation in subroutine "set_obj" available in the "ext_cyl.usr". Overall, I understand that there is a conditional statement to find all the elements (or element edges) with "W" boundary condition and
then an object is created for each of those elements. It would be very helpful is someone could point me to a documentation/help, which gives an explanation for the following functions or variables which have been used in the "set_obj" subroutine.<br>
<br>
</div>
<div>1. variable "hcode"<br>
</div>
<div>2. variable "lochis"<br>
</div>
<div>3. variable "nhis"<br>
</div>
<div>4. function or variable "cbc" <br>
</div>
<div>5. variable "nmember"<br>
</div>
<div>6. function or variable "object"<br>
<br>
</div>
<div>Any hints to understand this subroutine "set_obj" of the example model "ext_cyl" would be of great help for me to understand and develop the model for vortex-induced vibration.
<br>
<br>
</div>
<div>In addition, if there is already an example model to perform vortex-induced vibration, please let me know.<br>
<br>
</div>
<div>Thanks in advance. Any hints and tips are much appreciated. <br>
<br>
</div>
<div>Regards,<br>
</div>
<div>Ravi<br>
<br>
</div>
<div><br>
</div>
<div>-- <br>
Ravi Kumar T R<br>
<br>
University of Illinois at Urbana-Champaign<br>
ph : <a href="tel:%2B1%20217%20778%207538" value="+12177787538" target="_blank">+1 217 778 7538</a><br>
<a href="http://lndvl.mechse.illinois.edu/" target="_blank">http://lndvl.mechse.illinois.edu/</a>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>