<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;"><br>
Dear Stefan and Aleks;<br>
<br>
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. <br>
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'.<br>
<br>
Many thanks;<br>
Azad<br>
<br>
-------------------------------------------------------------------------------------------------------------------<br>
<br>
From: nek5000-users-request@lists.mcs.anl.gov<br>
Reply-to: nek5000-users@lists.mcs.anl.gov<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Subject: Nek5000-users Digest, Vol 28, Issue 6<br>
Date: Tue, 07 Jun 2011 12:00:13 -0500 (06/07/2011 07:00:13 PM)<br>
<br>
<br>
Send Nek5000-users mailing list submissions to<br>
nek5000-users@lists.mcs.anl.gov<br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br>
or, via email, send a message with subject or body 'help' to<br>
nek5000-users-request@lists.mcs.anl.gov<br>
<br>
You can reach the person managing the list at<br>
nek5000-users-owner@lists.mcs.anl.gov<br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Nek5000-users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Periodic bc for bent 2D channel (nek5000-users@lists.mcs.anl.gov)<br>
2. Re: Periodic bc for bent 2D channel<br>
(nek5000-users@lists.mcs.anl.gov)<br>
3. Re: Periodic bc for bent 2D channel<br>
(nek5000-users@lists.mcs.anl.gov)<br>
4. Re: Periodic bc for bent 2D channel<br>
(nek5000-users@lists.mcs.anl.gov)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Tue, 7 Jun 2011 16:14:27 +0000<br>
From: nek5000-users@lists.mcs.anl.gov<br>
Subject: [Nek5000-users] Periodic bc for bent 2D channel<br>
To: "nek5000-users@lists.mcs.anl.gov"<br>
<nek5000-users@lists.mcs.anl.gov><br>
Message-ID: <F035F1DD99A16A45AC22A084B64B6B8F0514A8E3@EXDB4.ug.kth.se><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hello;<br>
<br>
<br>
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.<br>
<br>
<br>
Many thanks<br>
<br>
Azad<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20110607/cc2e0e8c/attachment-0001.htm><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Tue, 7 Jun 2011 18:23:38 +0200<br>
From: nek5000-users@lists.mcs.anl.gov<br>
Subject: Re: [Nek5000-users] Periodic bc for bent 2D channel<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Message-ID: <BANLkTindyX6AU2YKWZuq+1_4aZ05Pfoftw@mail.gmail.com><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
Hi Azad<br>
<br>
I guess you're looking for a cyclic boundary condition. We implemented<br>
something like this for 2D geometries some time ago.Just set the<br>
logical flag IFCYCLIC to true in your rea-file and give it a whirl.<br>
<br>
-Stefan<br>
<br>
On 6/7/11, nek5000-users@lists.mcs.anl.gov<br>
<nek5000-users@lists.mcs.anl.gov> wrote:<br>
> Hello;<br>
><br>
><br>
> I have a question regarding the Periodic BC in Nek5000. I already setup a<br>
> simulation, it is a 2D channel generated by genbox. Using userdat2 I would<br>
> be able able to bend it to (0,pi/2) as Half Arch with flow direction from<br>
> up-left to right-bottom; however when I tried to run it, it seems the<br>
> periodic bc do not take care of the change in the direction of the flow at<br>
> the boundary. Now since the left-side periodic domain is in the same<br>
> position as before, but the right-side one face downward with a 90 deg<br>
> rotation, I would like to define the Periodic bc as<br>
> (u_left,v_left)=(-v_right,u_right). Unfortunately after searching a lot in<br>
> the source code I did not find any relevant information regarding this<br>
> adjustment or at least re-define the Periodic bc. I really appreciate if you<br>
> could help me by this.<br>
><br>
><br>
> Many thanks<br>
><br>
> Azad<br>
><br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Tue, 7 Jun 2011 11:35:21 -0500 (CDT)<br>
From: nek5000-users@lists.mcs.anl.gov<br>
Subject: Re: [Nek5000-users] Periodic bc for bent 2D channel<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Message-ID:<br>
<842151599.12710.1307464521780.JavaMail.root@zimbra.anl.gov><br>
Content-Type: text/plain; charset=utf-8<br>
<br>
Hi Azad,<br>
<br>
<br>
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<br>
<br>
INPUT:69: COMMON /INPUT5/ XC(8,LELT),YC(8,LELT),ZC(8,LELT)<br>
<br>
which are a part of TOTAL<br>
<br>
<br>
The piece of the code should be something like<br>
<br>
subroutine usrdat<br>
include 'SIZE'<br>
include 'TOTAL'<br>
<br>
ntot = 8*nelt <br>
<br>
do i=1,ntot<br>
<br>
xc(i,1) = ...<br>
<br>
...<br>
<br>
<br>
Best,<br>
Aleks<br>
<br>
<br>
<br>
<br>
----- Original Message -----<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Sent: Tuesday, June 7, 2011 11:14:27 AM<br>
Subject: [Nek5000-users] Periodic bc for bent 2D channel<br>
<br>
<br>
<br>
<br>
<br>
Hello; <br>
<br>
<br>
<br>
<br>
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. <br>
<br>
<br>
<br>
<br>
Many thanks <br>
<br>
Azad <br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<br>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br>
<br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Tue, 7 Jun 2011 18:46:42 +0200<br>
From: nek5000-users@lists.mcs.anl.gov<br>
Subject: Re: [Nek5000-users] Periodic bc for bent 2D channel<br>
To: nek5000-users@lists.mcs.anl.gov<br>
Message-ID: <BANLkTinqw3XaWXezDJt_H2gNFxHu84dSKA@mail.gmail.com><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
It depends ...<br>
<br>
If you modify the element vertices you should do it in usrdat(), as<br>
suggested by Aleks.<br>
If you want modify the SEM mesh (in our case the GLL points) you're ok<br>
with usrdat2().<br>
<br>
-Stefan<br>
<br>
On 6/7/11, nek5000-users@lists.mcs.anl.gov<br>
<nek5000-users@lists.mcs.anl.gov> wrote:<br>
> Hi Azad,<br>
><br>
><br>
> I think that a correct way of implementing this (non-affine) transformation<br>
> would be in userdat instead of userdat2 -- there you can just adjust element<br>
> vertex coordinate arrays xc,yc,zc<br>
><br>
> INPUT:69: COMMON /INPUT5/ XC(8,LELT),YC(8,LELT),ZC(8,LELT)<br>
><br>
> which are a part of TOTAL<br>
><br>
><br>
> The piece of the code should be something like<br>
><br>
> subroutine usrdat<br>
> include 'SIZE'<br>
> include 'TOTAL'<br>
><br>
> ntot = 8*nelt<br>
><br>
> do i=1,ntot<br>
><br>
> xc(i,1) = ...<br>
><br>
> ...<br>
><br>
><br>
> Best,<br>
> Aleks<br>
><br>
><br>
><br>
><br>
> ----- Original Message -----<br>
> From: nek5000-users@lists.mcs.anl.gov<br>
> To: nek5000-users@lists.mcs.anl.gov<br>
> Sent: Tuesday, June 7, 2011 11:14:27 AM<br>
> Subject: [Nek5000-users] Periodic bc for bent 2D channel<br>
><br>
><br>
><br>
><br>
><br>
> Hello;<br>
><br>
><br>
><br>
><br>
> I have a question regarding the Periodic BC in Nek5000. I already setup a<br>
> simulation, it is a 2D channel generated by genbox. Using userdat2 I would<br>
> be able able to bend it to (0,pi/2) as Half Arch with flow direction from<br>
> up-left to right-bottom; however when I tried to run it, it seems the<br>
> periodic bc do not take care of the change in the direction of the flow at<br>
> the boundary. Now since the left-side periodic domain is in the same<br>
> position as before, but the right-side one face downward with a 90 deg<br>
> rotation, I would like to define the Periodic bc as<br>
> (u_left,v_left)=(-v_right,u_right). Unfortunately after searching a lot in<br>
> the source code I did not find any relevant information regarding this<br>
> adjustment or at least re-define the Periodic bc. I really appreciate if you<br>
> could help me by this.<br>
><br>
><br>
><br>
><br>
> Many thanks<br>
><br>
> Azad<br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> Nek5000-users@lists.mcs.anl.gov<br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> Nek5000-users@lists.mcs.anl.gov<br>
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br>
><br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<br>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br>
<br>
<br>
End of Nek5000-users Digest, Vol 28, Issue 6<br>
********************************************<br>
</div>
</body>
</html>