[Nek5000-users] 'comp_vort3' in 2D

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

Thanks for the clarification.

Saikat

Saikat Mukherjee,
PhD Student,
Paul Research Group - http://www.me.vt.edu/mpaul/
Engineering Science and Mechanics,
Virginia Tech.

On Sun, Mar 25, 2018 at 4:40 AM, <nek5000-users at lists.mcs.anl.gov> wrote:

> 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
>>
>> _______________________________________________
> 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/20180325/7208793a/attachment.html>