<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,Arial,Helvetica,sans-serif;" dir="ltr">
<p>Hi Stefan and Mattias,</p>
<p><br>
</p>
<p>Thank you both for your kind advice.</p>
<p><br>
</p>
<p>Stefan, I had a look at the new 'gfldr' subroutine. It seems that the spectral interpolation procedure is the same as before, however, the latest subroutine breaks down the number of fields to interpolate in one go. I still don't see why the old g2gi hangs
at 'findpts'. Is it mainly because of memory size? (I have been using my own copy of the legacy nek, so I'm quite keen on making it work with g2gi.)</p>
<p><br>
</p>
<p>Mattias, I will give it a go with your suggestion. By reducing the polynomial order, did you mean the order of the old field or the order of the new geometry field, or both?</p>
<p><br>
</p>
<p>Thank you again for both your help.</p>
<p><br>
</p>
<p>Best regards,</p>
<p>Tony<br>
</p>
<br>
<br>
<div style="color: rgb(19, 19, 18);">
<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-bounces@lists.mcs.anl.gov <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users-request@lists.mcs.anl.gov <nek5000-users-request@lists.mcs.anl.gov><br>
<b>Sent:</b> 27 March 2017 12:57<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Nek5000-users Digest, Vol 97, Issue 27</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Send Nek5000-users mailing list submissions to<br>
nek5000-users@lists.mcs.anl.gov<br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a previewremoved="true" id="LPlnk58952" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
or, via email, send a message with subject or body 'help' to<br>
nek5000-users-request@lists.mcs.anl.gov<br>
<br>
You can reach the person managing the list at<br>
nek5000-users-owner@lists.mcs.anl.gov<br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Nek5000-users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Interpolation between two meshes (nek5000-users@lists.mcs.anl.gov)<br>
2. Re: Interpolation between two meshes<br>
(nek5000-users@lists.mcs.anl.gov)<br>
3. Re: Copy velocity slice from one boundary to another<br>
(nek5000-users@lists.mcs.anl.gov)<br>
4. Robin Boundary Condition -- Navier's Slip Length<br>
(nek5000-users@lists.mcs.anl.gov)<br>
5. Re: Interpolation between two meshes<br>
(nek5000-users@lists.mcs.anl.gov)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Mon, 27 Mar 2017 04:37:46 +0000<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: Nek5000 <nek5000-users@lists.mcs.anl.gov><br>
Subject: [Nek5000-users] Interpolation between two meshes<br>
Message-ID:<br>
<mailman.3211.1490589471.2967.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Dear Neks,<br>
<br>
I have been using 'g2gi' to do interpolation between different meshes. Recently, I have noticed that when I interpolate between two big meshes (say 300000 elements in each mesh), the interpolation hangs at 'call findpts' and never moves on. I tried increasing
number of processors and changing 'lpart', but none of these worked.<br>
<br>
I'm wondering what could cause the interpolation to hang at 'findpts'. Has anyone encountered similar issue?<br>
<br>
Hope anyone could give me some advice. Thank you very much in advance.<br>
<br>
Best regards,<br>
Tony<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a previewremoved="true" id="LPlnk227998" href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/b238e605/attachment-0001.html">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/b238e605/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Mon, 27 Mar 2017 09:47:37 +0200<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
Message-ID:<br>
<mailman.3212.1490600871.2967.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset=utf-8<br>
<br>
Please update to the latest master. Note, g2gi was replaced by gfldr. <br>
<br>
Stefan<br>
<br>
-----Original message-----<br>
> From:nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
> Sent: Monday 27th March 2017 6:38<br>
> To: Nek5000 <nek5000-users@lists.mcs.anl.gov><br>
> Subject: [Nek5000-users] Interpolation between two meshes<br>
> <br>
> Dear Neks,<br>
> <br>
> I have been using 'g2gi' to do interpolation between different meshes. Recently, I have?noticed that when I interpolate between two big meshes (say 300000 elements in each mesh), the interpolation hangs at 'call findpts' and never moves on. I tried increasing<br>
> number of processors and changing 'lpart', but none of these worked.<br>
> <br>
> I'm wondering what could cause the interpolation to hang at 'findpts'. Has anyone encountered similar issue?<br>
> <br>
> Hope anyone could give me some advice. Thank you very much in advance.<br>
> <br>
> Best regards,<br>
> Tony<br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> Nek5000-users@lists.mcs.anl.gov<br>
> <a previewremoved="true" id="LPlnk803367" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Mon, 27 Mar 2017 16:47:00 +0700<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Subject: Re: [Nek5000-users] Copy velocity slice from one boundary to<br>
another<br>
Message-ID:<br>
<mailman.3216.1490608022.2967.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Thanks a lot for replying, your answer was very useful.<br>
<br>
I look on multimesh.f file.<br>
As i understood, i need set subroutine 'set_inflag' to look for 'P ' bc.<br>
After this i need to call 'get_values(field)' subroutine, which interpolate<br>
and put the value in valint(ix,iy,iz,ie,ifld) array.<br>
And then this array will be available in userbc in .usr file, to set it on<br>
inflow boundary.<br>
<br>
c-----------------------------------------------------------------------------------------------<br>
<br>
Vatslav,<br>
you will have to use the neknek script that allows you to exchange data<br>
between two Nek5000 sessions that are run simultaneously. You should have a<br>
look at the neknek examples to understand how to set up your usr file.<br>
Some of the subroutines that will be useful to you are located in the<br>
following link: <a previewremoved="true" id="LPlnk944812" href="http://www.mcs.anl.gov/research/projects/nek5/docs/">
http://www.mcs.anl.gov/research/projects/nek5/docs/</a><br>
html/multimesh_8f_source.html<br>
You will have to modify some of them to fit your needs since for example<br>
the subroutine ?set_inflag? is written to look for ?in? boundary conditions<br>
and replace them with ?v ?. You will need a subroutine that looks for ?in?<br>
and substitute it with either ?v ? or ?P ? depending on the which pipe<br>
you are at.<br>
Marco<br>
<br>
Hi, Neks!<br>
I have two pipes in my simulation. They are placed near each other on one<br>
axis with empty space between.<br>
One of them with periodic boundary conditions on z direction, another one -<br>
inflow -> outflow.<br>
How can I copy velocity slice from one of boundary on first pipe with<br>
periodic bc to the inflow boundary on second pipe?<br>
Regards, Vatslav<br>
<br>
c-----------------------------------------------------------------------------------------------<br>
<br>
Regards, Vatslav<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a previewremoved="true" id="LPlnk273894" href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/92c78de8/attachment-0001.html">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/92c78de8/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Mon, 27 Mar 2017 11:31:53 +0200<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Subject: [Nek5000-users] Robin Boundary Condition -- Navier's Slip<br>
Length<br>
Message-ID:<br>
<mailman.3222.1490614366.2967.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Neks,<br>
I wish to simulate a transitional channel flow whose wall respect a<br>
"slip-length" boundary condition, such as:<br>
<br>
$ [u + \delta \frac{\partial u}{\partial n}]_{wall}=0$<br>
<br>
Where \delta is the slip length<br>
<br>
^n<br>
| u u<br>
| |-->/ |----->/<br>
| | / | /<br>
| | / | /<br>
| | / | /<br>
|___|/_______| /_____<br>
| / ^<br>
|/ v \delta<br>
<br>
No-slip Slip-length<br>
<br>
(a better sketch can be found here:<br>
<a previewremoved="true" id="LPlnk669792" href="http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html">http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html</a> )
<div style="margin-bottom: 20px; overflow: auto; width: 100%; text-indent: 0px;" id="LPBorder_GT_14906290202490.2691863624740807">
<table role="presentation" style="width: 90%; background-color: rgb(255, 255, 255); position: relative; overflow: auto; padding-top: 20px; padding-bottom: 20px; margin-top: 20px; border-top: 1px dotted rgb(200, 200, 200); border-bottom: 1px dotted rgb(200, 200, 200);" id="LPContainer_14906290202170.413755950840868" cellspacing="0">
<tbody>
<tr style="border-spacing: 0px;" valign="top">
<td colspan="1" style="width: 250px; position: relative; display: table-cell; padding-right: 20px;" id="ImageCell_14906290202440.8854553625719841">
<div style="background-color: rgb(255, 255, 255); height: 190px; position: relative; margin: auto; display: table; width: 250px;" id="LPImageContainer_14906290202440.19879353770535368">
<a target="_blank" href="http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html" style="display: table-cell; text-align: center;" id="LPImageAnchor_14906290202460.45028870137887456"><img id="LPThumbnailImageID_14906290202460.41945734288710457" aria-label="Preview image with link selected. Double-tap to open the link." style="display: inline-block; max-width: 250px; max-height: 250px; height: 190px; width: 250px; border-width: 0px; vertical-align: bottom;" height="190" width="250" src="http://www.nature.com/nmat/journal/v2/n4/images/nmat857-f1.gif"></a></div>
</td>
<td colspan="2" style="vertical-align: top; position: relative; padding: 0px; display: table-cell;" id="TextCell_14906290202460.8930763948857062">
<div id="LPRemovePreviewContainer_14906290202460.8479600950416156"></div>
<div style="top: 0px; color: rgb(54, 159, 71); font-weight: 400; font-size: 21px; font-family: "wf_segoe-ui_light","Segoe UI Light","Segoe WP Light","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif; line-height: 21px;" id="LPTitle_14906290202460.8013806101660419">
<a target="_blank" href="http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html" style="text-decoration: none;" id="LPUrlAnchor_14906290202480.42011701503029397">Figure 1 : Low-friction flows of liquid at nanopatterned ...</a></div>
<div style="margin: 10px 0px 16px; color: rgb(102, 102, 102); font-weight: 400; font-family: "wf_segoe-ui_normal","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif; font-size: 14px; line-height: 14px;" id="LPMetadata_14906290202480.34224180957589667">
www.nature.com</div>
<div style="display: block; color: rgb(102, 102, 102); font-weight: 400; font-family: "wf_segoe-ui_normal","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif; font-size: 14px; line-height: 20px; max-height: 100px; overflow: hidden;" id="LPDescription_14906290202490.36178553513270073">
Nature Materials journal covers a range of topics within materials science, from materials engineering and structural materials (metals, alloys, ceramics, composites ...</div>
</td>
</tr>
</tbody>
</table>
</div>
<br>
<br>
<br>
I would like to know if this Robin boundary condition for velocity field is<br>
already implemented in the current code.<br>
<br>
Thank you<br>
<br>
Ciccio<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a previewremoved="true" id="LPlnk5611" href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/9922e2e7/attachment-0001.html">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/9922e2e7/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Mon, 27 Mar 2017 13:56:02 +0200<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
Message-ID:<br>
<mailman.3230.1490615875.2967.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="windows-1252"; Format="flowed"<br>
<br>
Hi Tony,<br>
I have encountered the same issue with g2gi when interpolating very <br>
large meshes. I don't know what is causing the issue but a workaround is <br>
to reduce the polynomial order. You obviously loose accuracy by doing <br>
this, but I still find it better than starting from scratch with the new <br>
mesh. Hope this helps.<br>
<br>
Best,<br>
Mattias<br>
<br>
<br>
On 03/27/2017 06:37 AM, nek5000-users@lists.mcs.anl.gov wrote:<br>
> Dear Neks,<br>
><br>
> I have been using 'g2gi' to do interpolation between different meshes. <br>
> Recently, I have noticed that when I interpolate between two big <br>
> meshes (say 300000 elements in each mesh), the interpolation hangs at <br>
> 'call findpts' and never moves on. I tried increasing number of <br>
> processors and changing 'lpart', but none of these worked.<br>
><br>
> I'm wondering what could cause the interpolation to hang at 'findpts'. <br>
> Has anyone encountered similar issue?<br>
><br>
> Hope anyone could give me some advice. Thank you very much in advance.<br>
><br>
> Best regards,<br>
> Tony<br>
><br>
><br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> Nek5000-users@lists.mcs.anl.gov<br>
> <a previewremoved="true" id="LPlnk286223" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a previewremoved="true" id="LPlnk580767" href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/031d792f/attachment.html">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/031d792f/attachment.html</a>><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<br>
<a previewremoved="true" id="LPlnk603251" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
End of Nek5000-users Digest, Vol 97, Issue 27<br>
*********************************************<br>
</div>
</span></font></div>
</div>
</body>
</html>