[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