[Nek5000-users] Periodic bc for bent 2D channel

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

It depends ...

If you modify the element vertices you should do it in usrdat(), as
suggested by Aleks.
If you want modify the SEM mesh (in our case the GLL points) you're ok
with usrdat2().

-Stefan

On 6/7/11, nek5000-users at lists.mcs.anl.gov
<nek5000-users at lists.mcs.anl.gov> wrote:
> 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
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>