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

Thu Jun 9 07:25:57 CDT 2011

Dear Stefan and Aleks;

With this E-mail I just wanted to say thanks. The IFCYCLIC was the option that I actually were looking for. It works very well, and by "2D" I presume you meant the 2D transformation or planner transformation, since I think it works as well for the 3D cases that transformed in plane.
Also as Stefan mentioned I was interested in transforming the mesh instead of vertexes. Though it was very interesting that I could change the position of vertexes in the 'userdat'.

Many thanks;
Azad

-------------------------------------------------------------------------------------------------------------------

From:     nek5000-users-request at lists.mcs.anl.gov
Reply-to:     nek5000-users at lists.mcs.anl.gov
To:     nek5000-users at lists.mcs.anl.gov
Subject:     Nek5000-users Digest, Vol 28, Issue 6
Date:     Tue, 07 Jun 2011 12:00:13 -0500 (06/07/2011 07:00:13 PM)

Send Nek5000-users mailing list submissions to
        nek5000-users at lists.mcs.anl.gov

To subscribe or unsubscribe via the World Wide Web, visit
        https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
or, via email, send a message with subject or body 'help' to
        nek5000-users-request at lists.mcs.anl.gov

You can reach the person managing the list at
        nek5000-users-owner at lists.mcs.anl.gov

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Nek5000-users digest..."

Today's Topics:

   1. Periodic bc for bent 2D channel (nek5000-users at lists.mcs.anl.gov)
   2. Re: Periodic bc for bent 2D channel
      (nek5000-users at lists.mcs.anl.gov)
   3. Re: Periodic bc for bent 2D channel
      (nek5000-users at lists.mcs.anl.gov)
   4. Re: Periodic bc for bent 2D channel
      (nek5000-users at lists.mcs.anl.gov)

----------------------------------------------------------------------

Message: 1
Date: Tue, 7 Jun 2011 16:14:27 +0000
From: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Periodic bc for bent 2D channel
To: "nek5000-users at lists.mcs.anl.gov"
        <nek5000-users at lists.mcs.anl.gov>
Message-ID: <F035F1DD99A16A45AC22A084B64B6B8F0514A8E3 at EXDB4.ug.kth.se>
Content-Type: text/plain; charset="iso-8859-1"

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20110607/cc2e0e8c/attachment-0001.htm>

------------------------------

Message: 2
Date: Tue, 7 Jun 2011 18:23:38 +0200
From: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Periodic bc for bent 2D channel
To: nek5000-users at lists.mcs.anl.gov
Message-ID: <BANLkTindyX6AU2YKWZuq+1_4aZ05Pfoftw at mail.gmail.com>
Content-Type: text/plain; charset=ISO-8859-1

Hi Azad

I guess you're looking for a cyclic boundary condition. We implemented
something like this for 2D geometries some time ago.Just set the
logical flag IFCYCLIC to true in your rea-file and give it a whirl.

-Stefan

On 6/7/11, nek5000-users at lists.mcs.anl.gov
<nek5000-users at lists.mcs.anl.gov> wrote:
> 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
>

------------------------------

Message: 3
Date: Tue, 7 Jun 2011 11:35:21 -0500 (CDT)
From: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Periodic bc for bent 2D channel
To: nek5000-users at lists.mcs.anl.gov
Message-ID:
        <842151599.12710.1307464521780.JavaMail.root at zimbra.anl.gov>
Content-Type: text/plain; charset=utf-8

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

------------------------------

Message: 4
Date: Tue, 7 Jun 2011 18:46:42 +0200
From: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Periodic bc for bent 2D channel
To: nek5000-users at lists.mcs.anl.gov
Message-ID: <BANLkTinqw3XaWXezDJt_H2gNFxHu84dSKA at mail.gmail.com>
Content-Type: text/plain; charset=ISO-8859-1

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
> _______________________________________________
> 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

End of Nek5000-users Digest, Vol 28, Issue 6
********************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20110609/66e24705/attachment.html>