[Nek5000-users] Nek5000-users Digest, Vol 111, Issue 3

Sat May 5 10:25:36 CDT 2018

Hello Paul and Ketan,

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?

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.

       1 0.1617415E+03  1.56161E-02  1.22442E-07  7.85399E-01  7.85399E-01
volflow Z
       2 0.1617415E+03  1.56215E-02  9.51058E-08  7.85399E-01  7.85399E-01
volflow Z
       3 0.1617415E+03  1.56186E-02  9.08498E-08  7.85399E-01  7.85399E-01
volflow Z
       4 0.1617415E+03  1.56207E-02  1.09035E-07  7.85399E-01  7.85399E-01
volflow Z
       5 0.1617416E+03  1.56204E-02  1.30839E-07  7.85399E-01  7.85399E-01
volflow Z
       6 0.1617416E+03  1.56200E-02  1.57003E-07  7.85399E-01  7.85399E-01
volflow Z
       7 0.1617416E+03  1.56198E-02  1.88400E-07  7.85399E-01  7.85399E-01
volflow Z
       8 0.1617417E+03  1.56196E-02  2.26078E-07  7.85399E-01  7.85399E-01
volflow Z

Kind regards,

Jian

On 4 May 2018 at 18:00, <nek5000-users-request at lists.mcs.anl.gov> wrote:

> Send Nek5000-users mailing list submissions to
>         nek5000-users at lists.mcs.anl.gov
>
> To subscribe or unsubscribe via the World Wide Web, visit
>         https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> or, via email, send a message with subject or body 'help' to
>         nek5000-users-request at lists.mcs.anl.gov
>
> You can reach the person managing the list at
>         nek5000-users-owner at lists.mcs.anl.gov
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Nek5000-users digest..."
>
>
> Today's Topics:
>
>    1. Re: ffx ffy ffz from current run (nek5000-users at lists.mcs.anl.gov)
>    2. Re: ffx ffy ffz from current run (nek5000-users at lists.mcs.anl.gov)
>    3. Postprocessing: Append data to output file
>       (nek5000-users at lists.mcs.anl.gov)
>    4. Re: Postprocessing: Append data to output file
>       (nek5000-users at lists.mcs.anl.gov)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Thu, 3 May 2018 17:36:06 +0000
> From: nek5000-users at lists.mcs.anl.gov
> To: "nek5000-users at lists.mcs.anl.gov"
>         <nek5000-users at lists.mcs.anl.gov>
> Subject: Re: [Nek5000-users] ffx ffy ffz from current run
> Message-ID:
>         <mailman.841.1525368969.86639.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="utf-8"
>
> Hi Jian,
>
> 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.
>
> 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.
>
>
>       if (nio.eq.0) write(6,1) istep,chv(icvflow)
>      $   ,scale,delta_flow,current_flow,flow_rate
>     1    format(i11,'  volflow ',a1,11x,1p4e13.4)
>
> c     Note, this scale factor corresponds to FFX, provided FFX has
> c     not also been specified in userf.   If ffx is also specified
> c     in userf then the true FFX is given by ffx_userf + scale.
>
>
> Ketan
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/
> attachments/20180503/19d76cf4/attachment-0001.html>
>
> ------------------------------
>
> Message: 2
> Date: Thu, 3 May 2018 23:03:18 +0000
> From: nek5000-users at lists.mcs.anl.gov
> To: "nek5000-users at lists.mcs.anl.gov"
>         <nek5000-users at lists.mcs.anl.gov>
> Subject: Re: [Nek5000-users] ffx ffy ffz from current run
> Message-ID:
>         <mailman.860.1525388608.86639.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="us-ascii"
>
>
> Dear Zhai,
>
>
> Assuming you have a periodic domain in direction "d" (d=1,2,3 or x,y,z),
> and
>
> that you have ffx=ffy=ffz=0 in your .usr file, then
>
>
> grep volf logfile
>
>
> will show a list of time step, time, and "ffd" + other columns.
>
>
> Column 3 will be "ffd"
>
>
> where ffd is the force in the direction "d" and d = |p54|.
>
>
> p54 = +/-  1  ==>   d=x direction
>
>            +/-  2 ==>         y direction
>
>            +/-  3 ==>         z direction
>
>
> If p54 < 0, then p55 determines the mean velocity in direction d
>
>
> If p54 > 0, then p55 determines the volumetric flux in direction d
>
>
> (Note that p55 < 0 would be the flow is in the reverse direction.)
>
>
> hth,
>
>
> Paul
>
>
> ________________________________
> From: Nek5000-users <nek5000-users-bounces at lists.mcs.anl.gov> on behalf
> of nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
> Sent: Thursday, May 3, 2018 11:47:44 AM
> To: nek5000-users at lists.mcs.anl.gov
> Subject: [Nek5000-users] ffx ffy ffz from current run
>
> Hi Nek users,
>
> How can I get the ffx ffy and ffz from the simulation with flowrate fixed,
> where ffx ffy and ffz are in usr file?
>
> Kind regards,
>
> Zhai
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/
> attachments/20180503/ecd674f3/attachment-0001.html>
>
> ------------------------------
>
> Message: 3
> Date: Fri, 4 May 2018 06:39:25 +0000
> From: nek5000-users at lists.mcs.anl.gov
> To: "nek5000-users at lists.mcs.anl.gov"
>         <nek5000-users at lists.mcs.anl.gov>
> Subject: [Nek5000-users] Postprocessing: Append data to output file
> Message-ID:
>         <mailman.885.1525415969.86639.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi,
>
>
> I am postprocessing output files in order to compute the vorticity.
>
> So far I have followed the instructions in https://lists.mcs.anl.gov/
> pipermail/nek5000-users/2010-March/000171.html,
>
> where it is explained how to output the vorticity components as velocity
> components.
>
>
> I however would like to make a simultaneous visualization of the vorticity
> field, the velocity field and some of the scalar fields.
>
> If postprocessing a file "mycase0.f000XX" - is it then possible to append
> the vorticity to this file?
>
> If this is not possible, how would I to in order to write the velocity and
> the vorticity in the same file?
>
>
> The interface of the subroutine outpost is outpost(vx,vy,vz,pr,t,nam), but
> I have 3 scalars in addition to the temperature.
>
> How can I write a file containing both the vorticity and the four scalar
> fields?
>
>
> Best,
>
>
> Johan
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/
> attachments/20180504/c8b262ce/attachment-0001.html>
>
> ------------------------------
>
> Message: 4
> Date: Fri, 4 May 2018 11:30:41 +0200
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov     <nek5000-users at lists.mcs.anl.gov>
> Subject: Re: [Nek5000-users] Postprocessing: Append data to output
>         file
> Message-ID:
>         <mailman.888.1525426277.86639.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="utf-8"
>
> Assuming you are using par, define a scalar:
>
> [SCALAR01] # vorticity
> solver = none
>
> Then, just store the vorticity in t(1,1,1,1,iscal +1) where iscal is your
> passive scalar index (in this example 1).
>
> On 4 May 2018, at 08:39, "nek5000-users at lists.mcs.anl.gov <mailto:
> nek5000-users at lists.mcs.anl.gov> " <nek5000-users at lists.mcs.anl.gov
> <mailto:nek5000-users at lists.mcs.anl.gov> > wrote:
>
> Hi,
>
>
>
> I am postprocessing output files in order to compute the vorticity.
>
> So far I have followed the instructions in https://lists.mcs.anl.gov/
> pipermail/nek5000-users/2010-March/000171.html,
>
> where it is explained how to output the vorticity components as velocity
> components.
>
>
>
> I however would like to make a simultaneous visualization of the vorticity
> field, the velocity field and some of the scalar fields.
>
> If postprocessing a file "mycase0.f000XX" - is it then possible to append
> the vorticity to this file?
>
> If this is not possible, how would I to in order to write the velocity and
> the vorticity in the same file?
>
>
>
> The interface of the subroutine outpost is outpost(vx,vy,vz,pr,t,nam), but
> I have 3 scalars in addition to the temperature.
>
> How can I write a file containing both the vorticity and the four scalar
> fields?
>
>
>
> Best,
>
>
>
> Johan
>
>
> _______________________________________________
>
> Nek5000-users mailing list
>
> Nek5000-users at lists.mcs.anl.gov <mailto:Nek5000-users at lists.mcs.anl.gov>
>
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/
> attachments/20180504/e0c1a83b/attachment-0001.html>
>
> ------------------------------
>
> Subject: Digest Footer
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
> ------------------------------
>
> End of Nek5000-users Digest, Vol 111, Issue 3
> *********************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180505/dd348295/attachment.html>