<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: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;" dir="ltr">
<p>Hi Paul,</p>
<p><br>
</p>
<p>Thanks for that, so you're saying that I can set up my box file by individually prescribing the coordinates of each element? I think I am able to do this, I still find, however, that when I post-process my results in Python and use the element map in the
0.f##### files that the elements are not mapped to the correct places. Are you saying that there's a way to pre-emptively fix this problem in the box file before I do my simulation?</p>
<p><br>
</p>
<p>Cheers,</p>
<p><br>
</p>
<p>Daniel<br>
</p>
<br>
<br>
<div style="color: rgb(49, 55, 57);">
<div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> Nek5000-users <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> 18 September 2017 13:04<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Nek5000-users Digest, Vol 103, Issue 17</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 href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" id="LPlnk701874" previewremoved="true">
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. Re: Nek5000-users Digest, Vol 103, Issue 12<br>
(nek5000-users@lists.mcs.anl.gov)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Mon, 18 Sep 2017 12:03:59 +0000<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: "nek5000-users@lists.mcs.anl.gov"<br>
<nek5000-users@lists.mcs.anl.gov><br>
Subject: Re: [Nek5000-users] Nek5000-users Digest, Vol 103, Issue 12<br>
Message-ID:<br>
<mailman.953.1505736244.10190.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
<br>
Daniel,<br>
<br>
<br>
If you want a custom distribution in any of the directions, you just list<br>
<br>
a positive value, nel, for the number of elements in that direction, and then<br>
<br>
provide a list with nel+1 entries. Note, the list of entries must be shorter than<br>
<br>
80 characters, so you have to wrap the line if it goes over 80 characters. This<br>
<br>
is because I wrote a parser to allow for comments in the box files and as such<br>
<br>
it only parses 80 character strings. Just so long as the sequence of strings has<br>
<br>
nel+1 values, it will keep reading the strings till satisfied. It will then go on to<br>
<br>
the next required read (e.g., the next list of element locations for the next<br>
<br>
direction).<br>
<br>
<br>
2D example:<br>
<br>
<br>
-2 9 # 2 elements in x, 9 in y<br>
<br>
0 1 0.5 x0,x1, ratio: element 2 is half the size of element 1<br>
<br>
0 .1 .2 .3<br>
<br>
.4 .5 .6 .7<br>
<br>
.8 .9 coordinates for y (comments only on last string)<br>
<br>
<br>
etc<br>
<br>
<br>
Paul<br>
<br>
<br>
________________________________<br>
From: 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>
Sent: Monday, September 18, 2017 4:02:28 AM<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Subject: Re: [Nek5000-users] Nek5000-users Digest, Vol 103, Issue 12<br>
<br>
<br>
Hi Paul,<br>
<br>
<br>
Sorry for not replying earlier, I didn't see your reply last week unfortunately. So in my case I'm creating a partition as follows:<br>
<br>
<br>
Box_1<br>
-20 -50 nel_1, nel_2<br>
0. 50. 0.8 x0, x1<br>
0. 100. 1.1 y0, y1<br>
W ,E ,W ,W (bc's) velocity (west, east, south, north)<br>
T ,E ,t ,W (bc's) temperature (west, east, south, north)<br>
Box_2<br>
-20 -50 nel_1, nel_2<br>
50. 100. 1.2 x0, x1<br>
0. 100. 1.1 y0, y1<br>
E ,W ,W ,W (bc's) velocity (west, east, south, north)<br>
E ,T ,t ,W (bc's) temperature (west, east, south, north)<br>
<br>
<br>
This makes two boxes, separated in the x-direction. Is what you're saying, then, that my two 20x50 boxes are stacked one on top of the other to, in effect, create one 20x100 box? If that's the case, how do I take my element map (which is just a list of 2000
numbers mapping local to global elements) and make sure that my two boxes are concatenated in the x-direction, instead of the z?<br>
<br>
<br>
Thanks in advance,<br>
<br>
<br>
Daniel<br>
<br>
________________________________<br>
From: Nek5000-users <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>
Sent: 15 September 2017 13:05:25<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Subject: Nek5000-users Digest, Vol 103, Issue 12<br>
<br>
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 href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" id="LPlnk14285" previewremoved="true">
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. Re: time integration scheme for temperature/scalar equations<br>
(nek5000-users@lists.mcs.anl.gov)<br>
2. Re: Nek5000-users Digest, Vol 103, Issue 9<br>
(nek5000-users@lists.mcs.anl.gov)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Fri, 15 Sep 2017 13:02:20 +0200<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Subject: Re: [Nek5000-users] time integration scheme for<br>
temperature/scalar equations<br>
Message-ID:<br>
<mailman.839.1505473384.10190.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="windows-1252"; Format="flowed"<br>
<br>
Hi,<br>
To be more precise. Routines "makeabq" and "makebdq" are called if you<br>
use hemholtz solver for passive scalars. The other option is cvode,<br>
which uses slightly different scheme. I guess hemholtz is default and to<br>
use cvode you have to add CVODE to preprocessing flag. About the link to<br>
doxygen. This one is outdated so it could be dangerous to use it. There<br>
is doxygen showing current beta2 version under<br>
<a href="https://adampep.github.io/KTH_Framework/files.html" id="LPlnk363085" previewremoved="true">https://adampep.github.io/KTH_Framework/files.html</a><br>
Regards<br>
Adam<br>
<br>
On 15/09/17 12:11, nek5000-users@lists.mcs.anl.gov wrote:<br>
><br>
> Hi Bolun Xu,<br>
><br>
> Yes,<br>
><br>
> you can see that in the conduct.f file( in which convection-diffusion<br>
> equation for temperature is solved):<br>
> <a href="http://www.mcs.anl.gov/research/projects/nek5/docs/html/conduct_8f_source.html" id="LPlnk993676" previewremoved="true">
http://www.mcs.anl.gov/research/projects/nek5/docs/html/conduct_8f_source.html</a><br>
><br>
> In the subroutine "makeabq" the same Adam-Bashforth coefficients AB()<br>
> are called and in the subroutine "makebdq", the BD() coefficients are<br>
> used again.<br>
><br>
> best,<br>
> Pierre-Emmanuel<br>
><br>
><br>
> On 2017-09-15 11:47, nek5000-users@lists.mcs.anl.gov wrote:<br>
>> Hi! All,<br>
>><br>
>> I learned from the subroutines makeabf and makebdf in navier1.f that the<br>
>> time integration scheme for momentum equation is the AB-CN scheme, that<br>
>> explicit schemes for nonlinear and forcing term, and implicit ones for<br>
>> viscous term. The prefactors of those schemes are stored in AB() and BD()<br>
>> arrays. Now my question is, whether the temperature/scalar equations also<br>
>> use the same schemes as momentum ones and also share the same prefactors in<br>
>> AB() / BD() arrays.<br>
>><br>
>> Any help is appreciated.<br>
>><br>
>> Regards.<br>
>><br>
>> ----------------------------------------------------------------<br>
>><br>
>> Bolun Xu<br>
>><br>
>> University of Science and Technology of China<br>
>><br>
>> Department of Modern Mechanics<br>
>><br>
>> Hefei, Anhui, China<br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>> _______________________________________________<br>
>> Nek5000-users mailing list<br>
>> Nek5000-users@lists.mcs.anl.gov<br>
>> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" id="LPlnk442620" previewremoved="true">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
><br>
><br>
><br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> Nek5000-users@lists.mcs.anl.gov<br>
> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" id="LPlnk810651" previewremoved="true">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170915/b385ae23/attachment-0001.html" id="LPlnk367876" previewremoved="true">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170915/b385ae23/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Fri, 15 Sep 2017 07:51:15 +0000<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: "nek5000-users@lists.mcs.anl.gov"<br>
<nek5000-users@lists.mcs.anl.gov><br>
Subject: Re: [Nek5000-users] Nek5000-users Digest, Vol 103, Issue 9<br>
Message-ID:<br>
<mailman.845.1505477125.10190.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Hi Goran,<br>
<br>
<br>
First of all thanks for having a look, my box file is as follows:<br>
<br>
<br>
plume_v2_largeIni.rea<br>
2 spatial dimension<br>
2 number of fields<br>
#<br>
#<br>
#<br>
#========================================================<br>
#<br>
Box_1<br>
-20 -50 nel_1, nel_2<br>
0. 50. 0.8 x0, x1<br>
0. 100. 1.1 y0, y1<br>
W ,E ,W ,W (bc's) velocity (west, east, south, north)<br>
T ,E ,t ,W (bc's) temperature (west, east, south, north)<br>
Box_2<br>
-20 -50 nel_1, nel_2<br>
50. 100. 1.2 x0, x1<br>
0. 100. 1.1 y0, y1<br>
E ,W ,W ,W (bc's) velocity (west, east, south, north)<br>
E ,T ,t ,W (bc's) temperature (west, east, south, north)<br>
<br>
I think where I'm struggling is in working out how the output (0.f#####) files are structured, or how their structure relates to the actual solution. I've managed to get the (2000, in this case) element mappings out of the files, I just can't then get this
mapping to actually map the elements correctly.<br>
<br>
Thanks again,<br>
<br>
Daniel<br>
<br>
<br>
________________________________<br>
From: Nek5000-users <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>
Sent: 14 September 2017 17:29:13<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Subject: Nek5000-users Digest, Vol 103, Issue 9<br>
<br>
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 href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" id="LPlnk654204" previewremoved="true">
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. Problems with post-processing multiple-boxes in MATLAB and<br>
Python (nek5000-users@lists.mcs.anl.gov)<br>
2. Re: Problems with post-processing multiple-boxes in MATLAB<br>
and Python (nek5000-users@lists.mcs.anl.gov)<br>
3. Re: Problems with post-processing multiple-boxes in MATLAB<br>
and Python (nek5000-users@lists.mcs.anl.gov)<br>
4. Re: proj_ortho error in pipe flow<br>
(nek5000-users@lists.mcs.anl.gov)<br>
5. problem with applying boundary conditions<br>
(nek5000-users@lists.mcs.anl.gov)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Thu, 14 Sep 2017 14:57:27 +0000<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: "nek5000-users@lists.mcs.anl.gov"<br>
<nek5000-users@lists.mcs.anl.gov><br>
Subject: [Nek5000-users] Problems with post-processing multiple-boxes<br>
in MATLAB and Python<br>
Message-ID:<br>
<mailman.763.1505402018.10190.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi All,<br>
<br>
<br>
I am using Nek5000 to simulate a buoyant plume rising from the centre of the floor of a square 2D domain. I am trying to use MATLAB and/or Python to post-process my data. Everything works fine when I post-process a simulation that was initialised with only
one box in my .box file. When, however, I run a simulation with multiple boxes (i.e. to cluster nodes), the local to global mapping of elements appears to go wrong and puts data into the wrong places. When plotting a contour plot of coordinates for example,
where I would expect a constantly changing colour gradient, I instead get back colour gradients in four quarters.<br>
<br>
<br>
I haven't used this list before so if I need to provide examples of images or code I can do. If anybody else has had this problem and has found a way to overcome it by finding the correct mapping some help would be very much appreciated.<br>
<br>
<br>
Best regards,<br>
<br>
<br>
Daniel Ward<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170914/94fe07e8/attachment-0001.html" id="LPlnk425037" previewremoved="true">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170914/94fe07e8/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Thu, 14 Sep 2017 15:20:27 +0000<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: shivang agarwal <nek5000-users@lists.mcs.anl.gov><br>
Subject: Re: [Nek5000-users] Problems with post-processing<br>
multiple-boxes in MATLAB and Python<br>
Message-ID:<br>
<mailman.765.1505402478.10190.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Daniel,<br>
<br>
I've had a similar problem before. Can you send your .box file so we can have a look?<br>
<br>
Best regards,<br>
<br>
Goran<br>
<br>
On Sep 14, 2017, at 11:13, "nek5000-users@lists.mcs.anl.gov<mailto:nek5000-users@lists.mcs.anl.gov>" <nek5000-users@lists.mcs.anl.gov<mailto:nek5000-users@lists.mcs.anl.gov>> wrote:<br>
<br>
<br>
Hi All,<br>
<br>
<br>
I am using Nek5000 to simulate a buoyant plume rising from the centre of the floor of a square 2D domain. I am trying to use MATLAB and/or Python to post-process my data. Everything works fine when I post-process a simulation that was initialised with only
one box in my .box file. When, however, I run a simulation with multiple boxes (i.e. to cluster nodes), the local to global mapping of elements appears to go wrong and puts data into the wrong places. When plotting a contour plot of coordinates for example,
where I would expect a constantly changing colour gradient, I instead get back colour gradients in four quarters.<br>
<br>
<br>
I haven't used this list before so if I need to provide examples of images or code I can do. If anybody else has had this problem and has found a way to overcome it by finding the correct mapping some help would be very much appreciated.<br>
<br>
<br>
Best regards,<br>
<br>
<br>
Daniel Ward<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<mailto:Nek5000-users@lists.mcs.anl.gov><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" id="LPlnk855779" previewremoved="true">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170914/e8a271d0/attachment-0001.html" id="LPlnk432709" previewremoved="true">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170914/e8a271d0/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Thu, 14 Sep 2017 15:33:07 +0000<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: "nek5000-users@lists.mcs.anl.gov"<br>
<nek5000-users@lists.mcs.anl.gov><br>
Subject: Re: [Nek5000-users] Problems with post-processing<br>
multiple-boxes in MATLAB and Python<br>
Message-ID:<br>
<mailman.767.1505403191.10190.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
<br>
Hi All,<br>
<br>
<br>
Genbox generates a tensor-product array of elements with layout:<br>
<br>
<br>
(ix , iy, iz) (ix fastest moving, iy 2nd fastest, iz 3rd fastest)<br>
<br>
<br>
If you want to have multiple boxes that continue in this ordering, you could<br>
<br>
do the following, say:<br>
<br>
<br>
box (1)<br>
<br>
-5 -5 -5<br>
<br>
-1 1 1.0<br>
<br>
-1 1 1.0<br>
<br>
0 2 1.0<br>
<br>
<br>
box(2)<br>
<br>
-5 -5 -10<br>
<br>
-1 1 1<br>
<br>
-1 1 1<br>
<br>
1.0 9.0<br>
<br>
<br>
etc<br>
<br>
<br>
This would stack a 5x5x10 box on top of a 5x5x5 box<br>
<br>
<br>
and topologically is the same as a 5x5 x (15) box<br>
<br>
<br>
... i.e., things get concatenated in the z direction.<br>
<br>
<br>
<br>
Also, note that you can use custom distributions in each direction, e.g.,<br>
<br>
<br>
box<br>
<br>
3 3 -5<br>
<br>
0 0.1 0.3 10<br>
<br>
0 .9 22 33<br>
<br>
-5 9 1.5<br>
<br>
<br>
which would give spacing [ 0 .1 .3 10] in the x direction<br>
<br>
<br>
[0 .9 22 33] in y<br>
<br>
<br>
and geometric partition of 5 elements in z on [-5:9]<br>
<br>
with size ratio 1.5 between each adjacent z level.<br>
<br>
<br>
hth<br>
<br>
<br>
Paul<br>
<br>
<br>
________________________________<br>
From: 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>
Sent: Thursday, September 14, 2017 10:20:27 AM<br>
To: shivang agarwal<br>
Subject: Re: [Nek5000-users] Problems with post-processing multiple-boxes in MATLAB and Python<br>
<br>
Hi Daniel,<br>
<br>
I've had a similar problem before. Can you send your .box file so we can have a look?<br>
<br>
Best regards,<br>
<br>
Goran<br>
<br>
On Sep 14, 2017, at 11:13, "nek5000-users@lists.mcs.anl.gov<mailto:nek5000-users@lists.mcs.anl.gov>" <nek5000-users@lists.mcs.anl.gov<mailto:nek5000-users@lists.mcs.anl.gov>> wrote:<br>
<br>
<br>
Hi All,<br>
<br>
<br>
I am using Nek5000 to simulate a buoyant plume rising from the centre of the floor of a square 2D domain. I am trying to use MATLAB and/or Python to post-process my data. Everything works fine when I post-process a simulation that was initialised with only
one box in my .box file. When, however, I run a simulation with multiple boxes (i.e. to cluster nodes), the local to global mapping of elements appears to go wrong and puts data into the wrong places. When plotting a contour plot of coordinates for example,
where I would expect a constantly changing colour gradient, I instead get back colour gradients in four quarters.<br>
<br>
<br>
I haven't used this list before so if I need to provide examples of images or code I can do. If anybody else has had this problem and has found a way to overcome it by finding the correct mapping some help would be very much appreciated.<br>
<br>
<br>
Best regards,<br>
<br>
<br>
Daniel Ward<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<mailto:Nek5000-users@lists.mcs.anl.gov><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" id="LPlnk692471" previewremoved="true">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170914/4ab9769c/attachment-0001.html" id="LPlnk225478" previewremoved="true">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170914/4ab9769c/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Thu, 14 Sep 2017 15:34:39 +0000<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: "nek5000-users@lists.mcs.anl.gov"<br>
<nek5000-users@lists.mcs.anl.gov><br>
Subject: Re: [Nek5000-users] proj_ortho error in pipe flow<br>
Message-ID:<br>
<mailman.770.1505403282.10190.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="windows-1252"<br>
<br>
<br>
Hi Steffen,<br>
<br>
<br>
Thanks for bringing this to our attention.<br>
<br>
<br>
Adam and I have been looking into this... (no doubt, it's something<br>
<br>
stupid done by me awhile back).<br>
<br>
<br>
Hope to have it resolved by next week.<br>
<br>
<br>
Best,<br>
<br>
<br>
Paul<br>
<br>
<br>
________________________________<br>
From: 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>
Sent: Thursday, September 14, 2017 8:18:27 AM<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Subject: Re: [Nek5000-users] proj_ortho error in pipe flow<br>
<br>
Hi,<br>
<br>
I have made some more tests to find the cause of this issue.<br>
Attached is a laminar pipe flow simulation with temperature and one<br>
passive scalar. The original setup is from Paul's "nusselt.tgz", where<br>
he answered another question of mine.<br>
<br>
I found that I can reproduce the same error I see in my turbulent<br>
simulations when there is an additional passive scalar with a different<br>
type of boundary condition, i.e. Neumann boundary condition at the pipe<br>
wall for temperature and Dirichlet boundary condition for the passive<br>
scalar.<br>
<br>
Can someone please check my setup and tell me why this causes the<br>
projection to fail?<br>
<br>
<br>
I hope to save some computational resources by using the projection.<br>
However, if it is not possible, I will keep param(94) = 0.<br>
<br>
Kind Regards,<br>
Steffen Straub<br>
<br>
--<br>
Karlsruhe Institute of Technology (KIT)<br>
Institute of Fluid Mechanics<br>
<br>
M.Sc. Steffen Straub<br>
Doctoral Researcher<br>
<br>
Kaiserstra?e 10<br>
Building 10.23<br>
76131 Karlsruhe, Germany<br>
<br>
Phone: +49 721 608-43027<br>
E-mail: steffen.straub@kit.edu<br>
Web: <a href="http://www.istm.kit.edu" id="LPlnk738293" previewremoved="true">http://www.istm.kit.edu</a>
<div id="LPBorder_GT_15057414194940.15931319070604988" style="margin-bottom: 20px; overflow: auto; width: 100%; text-indent: 0px;">
<table id="LPContainer_15057414194840.08258323161507208" 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);" role="presentation" cellspacing="0">
<tbody>
<tr style="border-spacing: 0px;" valign="top">
<td id="TextCell_15057414194880.7496341569732489" style="vertical-align: top; position: relative; padding: 0px; display: table-cell;" colspan="2">
<div id="LPRemovePreviewContainer_15057414194880.8207059982816429"></div>
<div id="LPTitle_15057414194880.10707463068426337" style="top: 0px; color: rgb(0, 120, 215); 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;">
<a id="LPUrlAnchor_15057414194900.132676917235796" style="text-decoration: none;" href="http://www.istm.kit.edu/" target="_blank">KIT - Institut für StrömungsmechanikStartseite</a></div>
<div id="LPMetadata_15057414194920.13752182250303113" 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;">
www.istm.kit.edu</div>
<div id="LPDescription_15057414194930.6617985319694651" 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;">
Aktuelle Stellenangebote - Tutoren gesucht! Wenn Sie als Student Interesse an einer Tutorentätigkeit haben, melden Sie sich bitte bei sekretariat@istm.kit.edu</div>
</td>
</tr>
</tbody>
</table>
</div>
<br>
<br>
KIT ? The Research University in the Helmholtz Association<br>
<br>
Since 2010, the KIT has been certified as a family-friendly university.<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170914/3306082d/attachment-0001.html" id="LPlnk969070" previewremoved="true">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170914/3306082d/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Thu, 14 Sep 2017 12:29:10 -0400<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Subject: [Nek5000-users] problem with applying boundary conditions<br>
Message-ID:<br>
<mailman.778.1505406553.10190.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear Neks,<br>
<br>
I am trying to run a simulation, replicating the Nek5000 run from this<br>
youtube video:<br>
<a href="https://www.youtube.com/watch?v=RmmN-KFSvAY" id="LPlnk842615" previewremoved="true">https://www.youtube.com/watch?v=RmmN-KFSvAY</a>
<div id="LPBorder_BVTaHR0cHM6Ly93d3cueW91dHViZS5jb20vd2F0Y2g:dj1SbW1OLUtGU3ZBWQ.._15057414226590.3661356573067165" style="margin-bottom: 20px; overflow: auto; width: 100%; text-indent: 0px;">
<table id="LPContainer_15057414226370.11833918968539603" 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);" role="presentation" cellspacing="0">
<tbody>
<tr style="border-spacing: 0px;" valign="top">
<td id="ImageCell_15057414226390.6435275563013065" style="width: 250px; position: relative; display: table-cell; padding-right: 20px;" colspan="1">
<div id="LPImageContainer_15057414226400.4416967528323523" style="background-color: rgb(255, 255, 255); height: 126px; position: relative; margin: auto; width: 250px;">
<a id="LPImageAnchor_15057414226430.45281044033885753" style="display: table-cell; text-align: center;" href="https://www.youtube.com/watch?v=RmmN-KFSvAY" target="_blank"><img style="display: inline-block; max-width: 250px; max-height: 250px; height: 126px; width: 250px; border-width: 0px; vertical-align: bottom;" id="LPThumbnailImageID_15057414226480.49121627633005993" height="126" width="250" src="https://www.bing.com/th?id=OVP.XZlUydh0y5CDPktWEircHgEsCY&pid=Api"><video style="display: none; max-width: 250px; max-height: 250px; height: 128px; width: 250px; border-width: 0px; vertical-align: bottom;" aria-label="Preview image with link selected. Double-tap to open the link."><source src="https://www.bing.com/th?id=OM2.36BwthA6vvuNLQ_1499566665&pid=Api"></video></a></div>
</td>
<td id="TextCell_15057414226530.015460458821632939" style="vertical-align: top; position: relative; padding: 0px; display: table-cell;" colspan="2">
<div id="LPRemovePreviewContainer_15057414226530.7610334282678067"></div>
<div id="LPTitle_15057414226530.5431119746575268" style="top: 0px; color: rgb(0, 120, 215); 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;">
<a id="LPUrlAnchor_15057414226550.6293445411749663" style="text-decoration: none;" href="https://www.youtube.com/watch?v=RmmN-KFSvAY" target="_blank">nek5000: Fluid flow through orifice</a></div>
<div id="LPMetadata_15057414226550.6592983453215557" 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;">
www.youtube.com</div>
<div id="LPDescription_15057414226570.6584386298458718" 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;">
Showing the magnitude of the velocity Orifice width 40% of channel width Orifice thickness 10% of channel width Incompressible Newtonian fluid Reynolds number 1000</div>
</td>
</tr>
</tbody>
</table>
</div>
<br>
<br>
The inlet BC is constant velocity (which is set to 1.0), and as can be seen<br>
from the video, the velocity magnitude reasonably away from the inlet is<br>
zero at starting times. However, in all of my simulations, I see the<br>
velocity magnitude as one at all the interior points away from the inlet<br>
from the beginning (the inlet BC ux=1.0 I am applying is somehow propagated<br>
to all interior elements). Am I doing something wrong somewhere, or is this<br>
the expected behavior? Please see my userbc() and useric() from my .usr<br>
file below:<br>
<br>
subroutine userbc (ix,iy,iz,iside,eg)<br>
include 'SIZE'<br>
include 'TOTAL'<br>
include 'NEKUSE'<br>
integer ix,iy,iz,eg,e,i,n<br>
<br>
CC "My inlet is at x=0"<br>
<br>
if (x.le.1e-4) then<br>
ux=1.0<br>
endif<br>
<br>
return<br>
end<br>
<br>
subroutine useric (ix,iy,iz,eg)<br>
include 'SIZE'<br>
include 'TOTAL'<br>
include 'NEKUSE'<br>
<br>
ux=0<br>
uy=0<br>
<br>
return<br>
end<br>
<br>
I believe userbc() is executed for only 'v' in my case. For the other<br>
constraints on the elements (like 'W', 'O', 'E'), I think userbc() would<br>
not run (please correct me if I am wrong, I am learning here). Since the<br>
above functions wouldn't work as expected, I changed userchk() and removed<br>
the above if statements in userbc(). My userchk() is as follows:<br>
<br>
subroutine userchk<br>
include 'SIZE'<br>
include 'TOTAL'<br>
integer i,n<br>
<br>
n = nx1*ny1*nz1*nelv<br>
<br>
if (istep.eq.0) then<br>
do i=1,n<br>
if (xm1(i,1,1,1).le.1e-4) then<br>
vx(i,1,1,1)=1<br>
endif<br>
enddo<br>
endif<br>
<br>
return<br>
end<br>
<br>
This didn't work either. Can somebody point me in the right direction?<br>
<br>
Thanks,<br>
Ajay<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170914/bb155242/attachment.html" id="LPlnk74900" previewremoved="true">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170914/bb155242/attachment.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" id="LPlnk22617" previewremoved="true">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
------------------------------<br>
<br>
End of Nek5000-users Digest, Vol 103, Issue 9<br>
*********************************************<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170915/322d13e9/attachment.html" id="LPlnk63266" previewremoved="true">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170915/322d13e9/attachment.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" id="LPlnk388760" previewremoved="true">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
------------------------------<br>
<br>
End of Nek5000-users Digest, Vol 103, Issue 12<br>
**********************************************<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170918/3edb650e/attachment.html" id="LPlnk42701" previewremoved="true">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170918/3edb650e/attachment.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" id="LPlnk583401" previewremoved="true">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
------------------------------<br>
<br>
End of Nek5000-users Digest, Vol 103, Issue 17<br>
**********************************************<br>
</div>
</span></font></div>
</div>
</body>
</html>