<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<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">PS - also, please have a look at the .usr file in the previous post.</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, February 9, 2018 9:15:15 AM<br>
<b>To:</b> nek5000-users<br>
<b>Subject:</b> Re: [Nek5000-users] angular velocity to the mesh velocity</font>
<div> </div>
</div>
<div>
<div style="font-family:arial,helvetica,sans-serif; font-size:12pt; color:#000000">
<div><style>
<!--
@font-face
{font-family:"Cambria Math"}
@font-face
{font-family:Calibri}
p.x_x_MsoNormal, li.x_x_MsoNormal, div.x_x_MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif}
a:x_x_link, span.x_x_MsoHyperlink
{color:#0563C1;
text-decoration:underline}
a:x_x_visited, span.x_x_MsoHyperlinkFollowed
{color:#954F72;
text-decoration:underline}
p.x_x_msonormal0, li.x_x_msonormal0, div.x_x_msonormal0
{margin-right:0in;
margin-left:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif}
span.x_x_EmailStyle18
{font-family:"Calibri",sans-serif;
color:windowtext}
.x_x_MsoChpDefault
{font-size:10.0pt}
@page WordSection1
{margin:1.0in 1.0in 1.0in 1.0in}
-->
</style><style style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style></div>
<div>Thank you very much Paul for your detailed replay and I am very glad about it.<br>
</div>
<div><br>
</div>
<div>But I am sorry if i wrongly said my question because the answer you said was not the one I was expecting. The thing which you have said I had already taken into account (the blending function ). In fact we already had a talk about this previously and
I also had developed another blending function with another equation. </div>
<div><br>
</div>
<div>my question was is apart from wx,wy,wz. Is there any way to give the angular velocity to the mesh (like wr,,wtheta,wz). Since my proile pitch about an axis I will have an angular velocity and I am bit confused how should I give this angular velocity as
an input to the mesh velocity.<br>
</div>
<div><br>
</div>
<div>Thanks with regards<br>
</div>
<div>Sijo GEORGE<br>
</div>
<div><br>
</div>
<hr id="x_zwchr">
<div><b>De: </b>"nek5000-users" <nek5000-users@lists.mcs.anl.gov><br>
<b>À: </b>"nek5000-users" <nek5000-users@lists.mcs.anl.gov><br>
<b>Envoyé: </b>Vendredi 9 Février 2018 16:02:34<br>
<b>Objet: </b>Re: [Nek5000-users] angular velocity to the mesh velocity<br>
</div>
<div><br>
</div>
<div>
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Sorry for the slow response but was on travel yesterday.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">For moving meshes with<span style="font-size:12pt"> rigid-body motion of pitch+translation, say,</span></p>
<p style="margin-top:0; margin-bottom:0">you essentially need a smooth function that is 1 at the moving boundary and 0 at the far field boundaries. Call this function h(X,t)=h(x,y,z,t).</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Then your mesh velocity (W=[wx,wy,wz]), which is required by Nek's ALE formulation, is given by</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0"> Translational Velocity(t)*h + Rotational Velocity*h</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Where h==1 (i.e., on the boundary) you will have the correct velocity.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Where h==0 (on the far field), you will have no mesh velocity. </p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">h is in effect a blending function between these regions. </p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">My favorite way of generating h is to solve a Laplace equation with a variable diffusivity. A high diffusivity near the object ensures h=1 near</p>
<p style="margin-top:0; margin-bottom:0">the object and thus elements near the object move _with_ the object, which means that they do not get excessively deformed.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">I attach a .usr file from a case I ran several years ago -- I've not tested it recently so please take it as a guide and not a solution to your problem.</p>
<p style="margin-top:0; margin-bottom:0">Look for "Greens Function" to see how to compute h.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Here is a paper describing the approach</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0"><a href="http://www.mcs.anl.gov/~fischer/pubhtml/ccfd.pdf" id="LPlnk675412" target="_blank">http://www.mcs.anl.gov/~fischer/pubhtml/ccfd.pdf</a><br>
</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>
<p style="margin-top:0; margin-bottom:0">PS - Note - Nek requires _only_ the mesh velocity -- mesh position</p>
<p style="margin-top:0; margin-bottom:0">is updated by Nek.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Note that you may, or may not, need to recompute the Green's function.... It is often a smooth function even after the mesh is deformed, so you can keep using it.... I suggest doing some experiments for your particular
case.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
</div>
<hr style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><span color="#000000" face="Calibri, sans-serif" style="font-size:11pt; color:#000000; font-family:Calibri,sans-serif"><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, February 9, 2018 8:36:39 AM<br>
<b>To:</b> nek5000-users<br>
<b>Subject:</b> Re: [Nek5000-users] angular velocity to the mesh velocity</span>
<div> </div>
</div>
<div>
<div style="font-family:arial,helvetica,sans-serif; font-size:12pt; color:#000000">
<div></div>
<div>Thank you for your replay Marco.<br>
</div>
<div><br>
</div>
<div>I will explain from the begening clearly. My current task is to validate Fluid structre interaction on aero foil(NACA profile). The equations which i have to solve I have attached for aero foil. Before this case, I had succesfully done a case of vortex
induced vibration of cylinder in NEK5000. And the code below shows a part of the forced oscillating cylinder. From this code I learned that the body motion prescribed through the mesh velocities but I also think that I can only give translational velocity
components in NEK5000 to presribe a body motion (mesh as well).<br>
</div>
<div><br>
</div>
<div>Since my aero foil has pitching (it rotates up and down about an axis at angle theta). Since it has to be a FSI problem my plan is by using " call torque_calc(scale,x0,.true.,.true.)" I will get the M (moment about the elastic axis (x0)) and it gives
as an input to the Equation of angular motion (attached with this email). Then using Range Kutta method I get the angular velocity of the profile (r, theta,z co ordinates). Then I am planning to conver this angular velocity into linear velocity components
(x, y, z direction) then impose into the wx,wy,wz components and I hope it leads to the movement of mesh thus the geometry.
<br>
</div>
<div><br>
</div>
<div>My question is could you think that its possible? Or is there any other simple way to give this angular velocity to my NACA profile about the elastic axis (x0)? And I feel difficulty to convert the angular velocity into linear velocty and gives its to
each nodes which is on "mv" B.C . <br>
</div>
<div>Could you please help me with this situation by an example if you have.<br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>Thanks in advance<br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div> c Cylinder velocity components<br>
<br>
ucx = 0.0<br>
ucy = amplitude * omega * cos(omega*time)<br>
ucz = 0.0<br>
<br>
</div>
<div><br>
</div>
<div> n = nx1*ny1*nz1*nelv<br>
</div>
<div> do i=1,n ! Translational velocity<br>
wx(i,1,1,1) = basev(i)*ucx ! component.<br>
wy(i,1,1,1) = basev(i)*ucy<br>
wz(i,1,1,1) = basev(i)*ucz<br>
<br>
umeshx(i,1,1,1) = wx(i,1,1,1)<br>
umeshy(i,1,1,1) = wy(i,1,1,1)<br>
umeshz(i,1,1,1) = wz(i,1,1,1)<br>
<br>
enddo<br>
<br>
</div>
<div><br>
</div>
<hr id="x_x_zwchr">
<div><b>De: </b>"nek5000-users" <nek5000-users@lists.mcs.anl.gov><br>
<b>À: </b>"nek5000-users" <Nek5000-users@lists.mcs.anl.gov><br>
<b>Envoyé: </b>Vendredi 9 Février 2018 15:03:26<br>
<b>Objet: </b>Re: [Nek5000-users] angular velocity to the mesh velocity<br>
</div>
<div><br>
</div>
<div>
<div class="x_x_WordSection1">
<p class="x_x_MsoNormal">Sijo,</p>
<p class="x_x_MsoNormal"> </p>
<p class="x_x_MsoNormal">The mv_cyl example is a 2-D moving body example. As for the mesh velocity, can you please be more specific?</p>
<p class="x_x_MsoNormal"> </p>
<p class="x_x_MsoNormal">Marco</p>
<p class="x_x_MsoNormal"> </p>
<div style="border:none; border-top:solid #B5C4DF 1.0pt; padding:3.0pt 0in 0in 0in">
<p class="x_x_MsoNormal"><b><span style="font-size:12.0pt; color:black">From: </span>
</b><span style="font-size:12.0pt; color:black">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>Reply-To: </b>"nek5000-users@lists.mcs.anl.gov" <Nek5000-users@lists.mcs.anl.gov><br>
<b>Date: </b>Friday, February 9, 2018 at 8:56 AM<br>
<b>To: </b>nek5000-users <Nek5000-users@lists.mcs.anl.gov><br>
<b>Subject: </b>[Nek5000-users] angular velocity to the mesh velocity</span></p>
</div>
<div>
<p class="x_x_MsoNormal"> </p>
</div>
<div>
<div>
<p class="x_x_MsoNormal"><a name="x_x__MailOriginalBody"><span style="font-size:12.0pt; font-family:"Arial",sans-serif; color:black"> </span></a><br>
</p>
</div>
<div>
<div>
<p class="x_x_MsoNormal"><span style=""><span style="font-size:12.0pt; font-family:"Arial",sans-serif; color:black">Hi Neks,</span></span></p>
</div>
<div>
<p class="x_x_MsoNormal"><span style=""><span style="font-size:12.0pt; font-family:"Arial",sans-serif; color:black"> </span></span></p>
</div>
<div>
<p class="x_x_MsoNormal"><span style=""><span style="font-size:12.0pt; font-family:"Arial",sans-serif; color:black">Its a long time that i came up with something. So here we go,</span></span></p>
</div>
<div>
<p class="x_x_MsoNormal"><span style=""><span style="font-size:12.0pt; font-family:"Arial",sans-serif; color:black">SInce I was doinf the FSI of aero foil, I have to deal with 2 degrees of freedom (its a 2D case) one in vertical and another in angular direction
about the elastic axis of the profile. I could code the vertical motion of the foil but before starting the angular motion I have some doubts. Is there anyway to give the angular velocities directly to the mesh velocities? Apart from the translational velocity
components Or Do you want me to covert the angluar velocities into translational components and give to mesh velocities? Do you have any example with a body which is rotating?</span></span></p>
</div>
<div>
<p class="x_x_MsoNormal"><span style=""><span style="font-size:12.0pt; font-family:"Arial",sans-serif; color:black"> </span></span></p>
</div>
<div>
<p class="x_x_MsoNormal"><span style=""><span style="font-size:12.0pt; font-family:"Arial",sans-serif; color:black">Thanks in advance</span></span></p>
</div>
<div>
<p class="x_x_MsoNormal"><span style=""><span style="font-size:12.0pt; font-family:"Arial",sans-serif; color:black"> </span></span></p>
</div>
<div>
<p class="x_x_MsoNormal"><span style=""><span style="font-size:12.0pt; font-family:"Arial",sans-serif; color:black">Sijo GEORGE</span></span></p>
</div>
</div>
</div>
</div>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<br>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br>
</div>
</div>
</div>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<br>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br>
</div>
</div>
</div>
</body>
</html>