<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;} p
{margin-top:0;
margin-bottom:0}p
{margin-top:0;
margin-bottom:0}--></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<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="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><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="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<div id="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="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="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> 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_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_divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="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> 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>
</body>
</html>