[Nek5000-users] Periodic bc for bent 2D channel
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Tue Jun 7 11:35:21 CDT 2011
Hi Azad,
I think that a correct way of implementing this (non-affine) transformation would be in userdat instead of userdat2 -- there you can just adjust element vertex coordinate arrays xc,yc,zc
INPUT:69: COMMON /INPUT5/ XC(8,LELT),YC(8,LELT),ZC(8,LELT)
which are a part of TOTAL
The piece of the code should be something like
subroutine usrdat
include 'SIZE'
include 'TOTAL'
ntot = 8*nelt
do i=1,ntot
xc(i,1) = ...
...
Best,
Aleks
----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Sent: Tuesday, June 7, 2011 11:14:27 AM
Subject: [Nek5000-users] Periodic bc for bent 2D channel
Hello;
I have a question regarding the Periodic BC in Nek5000. I already setup a simulation, it is a 2D channel generated by genbox. Using userdat2 I would be able able to bend it to (0,pi/2) as Half Arch with flow direction from up-left to right-bottom; however when I tried to run it, it seems the periodic bc do not take care of the change in the direction of the flow at the boundary. Now since the left-side periodic domain is in the same position as before, but the right-side one face downward with a 90 deg rotation, I would like to define the Periodic bc as (u_left,v_left)=(-v_right,u_right). Unfortunately after searching a lot in the source code I did not find any relevant information regarding this adjustment or at least re-define the Periodic bc. I really appreciate if you could help me by this.
Many thanks
Azad
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
More information about the Nek5000-users
mailing list