<div dir="ltr">Hello Paul and Ketan,<div><br></div><div>I really appreciate your reply. Do you mean the third column which is 1.56E-02, is the ffz force if para(54) and para(55) are 0?</div><div><br></div><div>To be honest, I tried a case that 1.56E-02 was used as ffz in the usr file while setting para(54) and (55) are 0. The case was set to start from where keeping flowrate stopped which have the right mean velocities. But the mean streamwise velocity decreased a lot. So I am wondering how much ffz should be. </div><div><br></div><div><div> 1 0.1617415E+03 1.56161E-02 1.22442E-07 7.85399E-01 7.85399E-01 volflow Z</div><div> 2 0.1617415E+03 1.56215E-02 9.51058E-08 7.85399E-01 7.85399E-01 volflow Z</div><div> 3 0.1617415E+03 1.56186E-02 9.08498E-08 7.85399E-01 7.85399E-01 volflow Z</div><div> 4 0.1617415E+03 1.56207E-02 1.09035E-07 7.85399E-01 7.85399E-01 volflow Z</div><div> 5 0.1617416E+03 1.56204E-02 1.30839E-07 7.85399E-01 7.85399E-01 volflow Z</div><div> 6 0.1617416E+03 1.56200E-02 1.57003E-07 7.85399E-01 7.85399E-01 volflow Z</div><div> 7 0.1617416E+03 1.56198E-02 1.88400E-07 7.85399E-01 7.85399E-01 volflow Z</div><div> 8 0.1617417E+03 1.56196E-02 2.26078E-07 7.85399E-01 7.85399E-01 volflow Z</div></div><div><br></div><div>Kind regards,</div><div><br></div><div>Jian </div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 4 May 2018 at 18:00, <span dir="ltr"><<a href="mailto:nek5000-users-request@lists.mcs.anl.gov" target="_blank">nek5000-users-request@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send Nek5000-users mailing list submissions to<br>
<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>mailman/listinfo/nek5000-users</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:nek5000-users-request@lists.mcs.anl.gov">nek5000-users-request@lists.<wbr>mcs.anl.gov</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:nek5000-users-owner@lists.mcs.anl.gov">nek5000-users-owner@lists.mcs.<wbr>anl.gov</a><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: ffx ffy ffz from current run (<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>)<br>
2. Re: ffx ffy ffz from current run (<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>)<br>
3. Postprocessing: Append data to output file<br>
(<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>)<br>
4. Re: Postprocessing: Append data to output file<br>
(<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>)<br>
<br>
<br>
------------------------------<wbr>------------------------------<wbr>----------<br>
<br>
Message: 1<br>
Date: Thu, 3 May 2018 17:36:06 +0000<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
To: "<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>"<br>
<<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>><br>
Subject: Re: [Nek5000-users] ffx ffy ffz from current run<br>
Message-ID:<br>
<<a href="mailto:mailman.841.1525368969.86639.nek5000-users@lists.mcs.anl.gov">mailman.841.1525368969.86639.<wbr>nek5000-users@lists.mcs.anl.<wbr>gov</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Jian,<br>
<br>
Not exactly sure what you mean, but if you are using p54, p55 to set your flowrate, then you can use your logfile to get what the corresponding forcing should be using ffx/ffy/ffz.<br>
<br>
As per drive2.f, if you are using param(54),param(55), then in logfile you can grep for ?volflow?, and ?scale? represents the ffx/ffy/ffz.<br>
<br>
<br>
if (nio.eq.0) write(6,1) istep,chv(icvflow)<br>
$ ,scale,delta_flow,current_<wbr>flow,flow_rate<br>
1 format(i11,' volflow ',a1,11x,1p4e13.4)<br>
<br>
c Note, this scale factor corresponds to FFX, provided FFX has<br>
c not also been specified in userf. If ffx is also specified<br>
c in userf then the true FFX is given by ffx_userf + scale.<br>
<br>
<br>
Ketan<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180503/19d76cf4/attachment-0001.html" rel="noreferrer" target="_blank">http://lists.mcs.anl.gov/<wbr>pipermail/nek5000-users/<wbr>attachments/20180503/19d76cf4/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Thu, 3 May 2018 23:03:18 +0000<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
To: "<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>"<br>
<<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>><br>
Subject: Re: [Nek5000-users] ffx ffy ffz from current run<br>
Message-ID:<br>
<<a href="mailto:mailman.860.1525388608.86639.nek5000-users@lists.mcs.anl.gov">mailman.860.1525388608.86639.<wbr>nek5000-users@lists.mcs.anl.<wbr>gov</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
<br>
Dear Zhai,<br>
<br>
<br>
Assuming you have a periodic domain in direction "d" (d=1,2,3 or x,y,z), and<br>
<br>
that you have ffx=ffy=ffz=0 in your .usr file, then<br>
<br>
<br>
grep volf logfile<br>
<br>
<br>
will show a list of time step, time, and "ffd" + other columns.<br>
<br>
<br>
Column 3 will be "ffd"<br>
<br>
<br>
where ffd is the force in the direction "d" and d = |p54|.<br>
<br>
<br>
p54 = +/- 1 ==> d=x direction<br>
<br>
+/- 2 ==> y direction<br>
<br>
+/- 3 ==> z direction<br>
<br>
<br>
If p54 < 0, then p55 determines the mean velocity in direction d<br>
<br>
<br>
If p54 > 0, then p55 determines the volumetric flux in direction d<br>
<br>
<br>
(Note that p55 < 0 would be the flow is in the reverse direction.)<br>
<br>
<br>
hth,<br>
<br>
<br>
Paul<br>
<br>
<br>
______________________________<wbr>__<br>
From: Nek5000-users <<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov">nek5000-users-bounces@lists.<wbr>mcs.anl.gov</a>> on behalf of <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>><br>
Sent: Thursday, May 3, 2018 11:47:44 AM<br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
Subject: [Nek5000-users] ffx ffy ffz from current run<br>
<br>
Hi Nek users,<br>
<br>
How can I get the ffx ffy and ffz from the simulation with flowrate fixed, where ffx ffy and ffz are in usr file?<br>
<br>
Kind regards,<br>
<br>
Zhai<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180503/ecd674f3/attachment-0001.html" rel="noreferrer" target="_blank">http://lists.mcs.anl.gov/<wbr>pipermail/nek5000-users/<wbr>attachments/20180503/ecd674f3/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Fri, 4 May 2018 06:39:25 +0000<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
To: "<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>"<br>
<<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>><br>
Subject: [Nek5000-users] Postprocessing: Append data to output file<br>
Message-ID:<br>
<<a href="mailto:mailman.885.1525415969.86639.nek5000-users@lists.mcs.anl.gov">mailman.885.1525415969.86639.<wbr>nek5000-users@lists.mcs.anl.<wbr>gov</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi,<br>
<br>
<br>
I am postprocessing output files in order to compute the vorticity.<br>
<br>
So far I have followed the instructions in <a href="https://lists.mcs.anl.gov/pipermail/nek5000-users/2010-March/000171.html" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>pipermail/nek5000-users/2010-<wbr>March/000171.html</a>,<br>
<br>
where it is explained how to output the vorticity components as velocity components.<br>
<br>
<br>
I however would like to make a simultaneous visualization of the vorticity field, the velocity field and some of the scalar fields.<br>
<br>
If postprocessing a file "mycase0.f000XX" - is it then possible to append the vorticity to this file?<br>
<br>
If this is not possible, how would I to in order to write the velocity and the vorticity in the same file?<br>
<br>
<br>
The interface of the subroutine outpost is outpost(vx,vy,vz,pr,t,nam), but I have 3 scalars in addition to the temperature.<br>
<br>
How can I write a file containing both the vorticity and the four scalar fields?<br>
<br>
<br>
Best,<br>
<br>
<br>
Johan<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180504/c8b262ce/attachment-0001.html" rel="noreferrer" target="_blank">http://lists.mcs.anl.gov/<wbr>pipermail/nek5000-users/<wbr>attachments/20180504/c8b262ce/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Fri, 4 May 2018 11:30:41 +0200<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>><br>
Subject: Re: [Nek5000-users] Postprocessing: Append data to output<br>
file<br>
Message-ID:<br>
<<a href="mailto:mailman.888.1525426277.86639.nek5000-users@lists.mcs.anl.gov">mailman.888.1525426277.86639.<wbr>nek5000-users@lists.mcs.anl.<wbr>gov</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Assuming you are using par, define a scalar:<br>
<br>
[SCALAR01] # vorticity<br>
solver = none<br>
<br>
Then, just store the vorticity in t(1,1,1,1,iscal +1) where iscal is your passive scalar index (in this example 1).<br>
<br>
On 4 May 2018, at 08:39, "<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a> <mailto:<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.<wbr>mcs.anl.gov</a>> " <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a> <mailto:<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.<wbr>mcs.anl.gov</a>> > wrote:<br>
<br>
Hi,<br>
<br>
<br>
<br>
I am postprocessing output files in order to compute the vorticity.<br>
<br>
So far I have followed the instructions in <a href="https://lists.mcs.anl.gov/pipermail/nek5000-users/2010-March/000171.html" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>pipermail/nek5000-users/2010-<wbr>March/000171.html</a>,<br>
<br>
where it is explained how to output the vorticity components as velocity components.<br>
<br>
<br>
<br>
I however would like to make a simultaneous visualization of the vorticity field, the velocity field and some of the scalar fields.<br>
<br>
If postprocessing a file "mycase0.f000XX" - is it then possible to append the vorticity to this file?<br>
<br>
If this is not possible, how would I to in order to write the velocity and the vorticity in the same file?<br>
<br>
<br>
<br>
The interface of the subroutine outpost is outpost(vx,vy,vz,pr,t,nam), but I have 3 scalars in addition to the temperature.<br>
<br>
How can I write a file containing both the vorticity and the four scalar fields?<br>
<br>
<br>
<br>
Best,<br>
<br>
<br>
<br>
Johan<br>
<br>
<br>
______________________________<wbr>_________________<br>
<br>
Nek5000-users mailing list<br>
<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.<wbr>gov</a> <mailto:<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.<wbr>mcs.anl.gov</a>> <br>
<br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>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/20180504/e0c1a83b/attachment-0001.html" rel="noreferrer" target="_blank">http://lists.mcs.anl.gov/<wbr>pipermail/nek5000-users/<wbr>attachments/20180504/e0c1a83b/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>mailman/listinfo/nek5000-users</a><br>
<br>
<br>
------------------------------<br>
<br>
End of Nek5000-users Digest, Vol 111, Issue 3<br>
******************************<wbr>***************<br>
</blockquote></div><br></div>