<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">Dear Valerio,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Just modify the routine I sent to compute the average, but now weight the average by unx to get the x-component.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Do the same with uny for y, and unz for z.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Best,</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> Saturday, May 12, 2018 5:44:50 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Re: [Nek5000-users] Integral on the wall</font>
<div> </div>
</div>
<meta content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
p
{margin-top:0;
margin-bottom:0}
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Dear Paul, </p>
<p>using the subroutine in zfunc.f I got the value of the pressure averaged on the surface. However I would like to compute the force on the wall as the line integral of the pressure times the normal (maybe my first question was unclear). Do you think this
can be achieved modifying a bit the subroutine in zfunc.f or should I do it from scratch?</p>
<p>Thanks. Best,</p>
<p><br>
</p>
<p>Valerio<br>
</p>
<div dir="ltr" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>Da:</b> Nek5000-users <nek5000-users-bounces@lists.mcs.anl.gov> per conto di nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
<b>Inviato:</b> mercoledì 9 maggio 2018 15:11<br>
<b>A:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Oggetto:</b> Re: [Nek5000-users] Integral on the wall</font>
<div> </div>
</div>
<div>
<p>Dear Paul,</p>
<p><br>
</p>
<p>Thank you for the help, I will test the routine zfunc.</p>
<p><br>
</p>
<p>Best,</p>
<p>Valerio<br>
</p>
<div dir="ltr" style="color:rgb(33,33,33)">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>Da:</b> Nek5000-users <nek5000-users-bounces@lists.mcs.anl.gov> per conto di nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
<b>Inviato:</b> mercoledì 9 maggio 2018 15:05<br>
<b>A:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Oggetto:</b> Re: [Nek5000-users] Integral on the wall</font>
<div> </div>
</div>
<div>
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,'Apple Color Emoji','Segoe UI Emoji',NotoColorEmoji,'Segoe UI Symbol','Android Emoji',EmojiSymbols">
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Dear Valerio,</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">I've not tested or even compiled the attached code...</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">But it should be fairly close to being functional.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">The idea would be to pass in the pressure on mesh 1, which</p>
<p style="margin-top:0; margin-bottom:0">means you need to interpolate from pr to pm1 ---- there are</p>
<p style="margin-top:0; margin-bottom:0">several routines floating around to do this.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Please let me know if you have any questions and if you</p>
<p style="margin-top:0; margin-bottom:0">are able to make this routine work.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Best,</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 that you should be able to send in "zm1" to this</p>
<p style="margin-top:0; margin-bottom:0">routine and get back the correct z distribution. This is what</p>
<p style="margin-top:0; margin-bottom:0">I would do if I wished to plot pr vs z. Also, you might want</p>
<p style="margin-top:0; margin-bottom:0">to dsavg() pm1 if you want a continuous pressure estimate.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><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> Wednesday, May 9, 2018 5:55:29 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Re: [Nek5000-users] Integral on the wall</font>
<div> </div>
</div>
<meta content="text/html; charset=us-ascii">
<div dir="ltr">
<div id="x_x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,'Apple Color Emoji','Segoe UI Emoji',NotoColorEmoji,'Segoe UI Symbol','Android Emoji',EmojiSymbols">
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Dear Valerio,</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">You would need to modify torque calc to do this,</p>
<p style="margin-top:0; margin-bottom:0">assuming that your objective is to compute a line</p>
<p style="margin-top:0; margin-bottom:0">integral of the pressure on the GLL points.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Assuming that your pipe mesh was made by </p>
<p style="margin-top:0; margin-bottom:0">sweeping the geometry in the z direction with </p>
<p style="margin-top:0; margin-bottom:0">n2to3, the modification should be quite easy.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">The routine I would start with would be drgtrq</p>
<p style="margin-top:0; margin-bottom:0">in navier5.f, which I see uses the two-index</p>
<p style="margin-top:0; margin-bottom:0">loop to traverse faces based on jskip1 and jskip2.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">For your case, I would use the facind() 3 index</p>
<p style="margin-top:0; margin-bottom:0">loop approach (e.g., as used in norcmp2 routine).</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Then you track the outer "k" loop and accumulate one</p>
<p style="margin-top:0; margin-bottom:0">value for each k value and for each z level of the elements.</p>
<p style="margin-top:0; margin-bottom:0">You will know the z level because it will be:</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0"> klevel_e = 1 + (eg-1)/nelxy</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">where nelxy is the number of elements in your pipe</p>
<p style="margin-top:0; margin-bottom:0">cross section.</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Hopefully this will get you pointed in the right direction...</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0">Best,</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"></p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
<p style="margin-top:0; margin-bottom:0"><br>
</p>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><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> Wednesday, May 9, 2018 3:13:19 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Integral on the wall</font>
<div> </div>
</div>
<meta content="text/html; charset=iso-8859-1">
<style type="text/css" style="">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Dear all,</p>
<p>I want to compute the integral of the pressure on pipe wall at different cross-sections. I tried using the routine set_obj, setting as many objects as the number of elements along the pipe length and then using the routine torque_calc. However I would like
to compute it on each circumference where GLL nodes are defined, for example if I use lx1 = 8 I would like to compute it on 8 cross-sections per element. Is there an efficient way to do this using the set_obj routine?</p>
<p><br>
</p>
<p>Thanks. Best regards.</p>
<p>Valerio<br>
</p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>