[Nek5000-users] lambda2 in turbChannel

[nek5000-users at lists.mcs.anl.gov]

Hi Tannoy,

Is there an exist subroutine to compute Q-criterion?

Thanks!

Wei

2016-06-13 17:01 GMT+02:00 <nek5000-users at lists.mcs.anl.gov>:

> You can run nek5000 in postprocessing mode (set NSTEPS = 0). Then in
> userchk, you can do the following
>
>      INCLUDE ’SIZE’
>      INCLUDE ’TOTAL’
>      INCLUDE ’RESTART’
>
>      character*80 filename(9999)
>      real l2(lx1,ly1,lz1,lelv)
>
>      ntot   = nx1*ny1*nz1*nelv
>
>         call load_ fld(filename)  ! select a filename fld file of your
> choice.
>         call lambdat2(l2)
>
>         call outpost(.vx,vy,vz,pr,l2,'lam')
>      enddo
>
> If your file is say turbChannel.fld015, the postprocessing will generate,
> lamturbChannel.fld01. You can type in visnek lamturbChannel in the terminal
> and then open visit and visualize the temperature variable. That will your
> l2.
>
> Best Regards,
> Tanmoy
>
> On Thu, Jun 9, 2016 at 5:13 PM, <nek5000-users at lists.mcs.anl.gov> wrote:
>
>> Hello,
>>
>>
>> In the userchk subroutine for the turbChannel example, lambda2 was
>> computed. How can I visualize the isosurfaces (streamwise structures) in
>> visit?
>>
>>
>>       if(mod(istep,iostep).eq.0 .and. istep.gt.0) then
>>         ! compute lambda2 vortex
>>         call lambda2(t(1,1,1,1,2))
>>         ! compute vorticity --- 3 components returned in t(...)
>>         if(ldimt.ge.5) call comp_vort3(t(1,1,1,1,3),wo1,wo2,vx,vy,vz)
>>       endif
>>
>>
>> Thank you,
>>
>> Virgilio A. Gomez
>>
>>
>>
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>
>>
>
> _______________________________________________
> Nek5000-users mailing list
> 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/20160613/11d34c88/attachment.html>