[Nek5000-users] 'comp_vort3' in 2D

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

Hi,
if you look into the code you see that for 2D cases the result it stored 
in vort(:,1), even though it is of course the z-vorticity.

00791       if (if3d) then
00792 c        work1=dw/dy ; work2=dv/dz
00793            call dudxyz(work1,w,rym1,sym1,tym1,jacm1,1,2)
00794            call dudxyz(work2,v,rzm1,szm1,tzm1,jacm1,1,3)
00795            call sub3(vort(1,1),work1,work2,ntot)
00796 c        work1=du/dz ; work2=dw/dx
00797            call dudxyz(work1,u,rzm1,szm1,tzm1,jacm1,1,3)
00798            call dudxyz(work2,w,rxm1,sxm1,txm1,jacm1,1,1)
00799            call sub3(vort(1,2),work1,work2,ntot)
00800 c        work1=dv/dx ; work2=du/dy
00801            call dudxyz(work1,v,rxm1,sxm1,txm1,jacm1,1,1)
00802            call dudxyz(work2,u,rym1,sym1,tym1,jacm1,1,2)
00803            call sub3(vort(1,3),work1,work2,ntot)
00804       else
00805 c        work1=dv/dx ; work2=du/dy
00806            call dudxyz(work1,v,rxm1,sxm1,txm1,jacm1,1,1)
00807            call dudxyz(work2,u,rym1,sym1,tym1,jacm1,1,2)
00808            call sub3(vort,work1,work2,ntot)
00809       endif

Philipp


On 2018-03-25 00:49, nek5000-users at lists.mcs.anl.gov wrote:
> Hello neks,
> 
> I am using 'comp_vort3' for a 2D flow, with a forcing in x direction,
> 
> call comp_vort3 (vort, w1, w2, vx, vy, vz)
> 
> call outpost (vort(1,1), vort(1,2), vort(1,3), pr, t, 'vrt')
> 
> 
> In this case, the x and y component of vorticity should be trivial, and 
> only z component should remain. However I am getting a non-zero 
> component of x-component, vort(1,1) while vort(1,2) and vort(1,3) are zero.
> 
> However if I add a z-direction, comp_vort3 returns a non-zero value of 
> vort(1,3) and the x and y-components are zero, which should be the case.
> 
> So my question is, in 2D, does comp_vort3 return only one value which is 
> in vort(1,1)? The vorticity fields in both cases look sort of similar.
> 
> Thanks,
> Saikat
> 
> 
> Saikat Mukherjee,
> PhD Student,
> Paul Research Group - http://www.me.vt.edu/mpaul/
> Engineering Science and Mechanics,
> Virginia Tech.
> 
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>