[Nek5000-users] Representing Curved Side in any plane

Sat Apr 3 20:37:14 CDT 2010

Hi Paul, 

Were you able to figure out the issue with Markus's midpoint case?  Thank you for looking into this! 

Kindly, 

Michael 

----- Original Message ----- 
From: nek5000-users at lists.mcs.anl.gov 
To: nek5000-users at lists.mcs.anl.gov 
Sent: Monday, March 29, 2010 11:37:43 AM GMT -06:00 US/Canada Central 
Subject: Re: [Nek5000-users] Representing Curved Side in any plane 

Hi Markus, 

Thanks -- I'll check into it. 

Paul 

On Mon, 29 Mar 2010, nek5000-users at lists.mcs.anl.gov wrote: 

> Hi, 
> 
> just double-checked, that's what I did. I only ran it for one time step without 
> any meaningful physics, though. 
> Attached are all case files. The nek version I am using is revision 456. 
> 
> Markus 
> 
> 
> 
> Quoting nek5000-users at lists.mcs.anl.gov: 
> 
>> 
>> Hi Markus, 
>> 
>> Did you visualize this with VisIt and with the geometry 
>> put out into (at least) the first .fld or .f file ? 
>> 
>> Paul 
>> 
>> 
>> On Mon, 29 Mar 2010, nek5000-users at lists.mcs.anl.gov wrote: 
>> 
>>> Hi, 
>>> 
>>> I was looking into the midpoint feature and generated a cube (1 element, x, 
>>> y, z from 0 to 1) in prenek, then converted all edges to the midpoint 
>>> notation with prenek, and then manipulated one edge in the rea file. This 
>> is 
>>> the resulting curved side section: 
>>>  ***** CURVED SIDE DATA ***** 
>>>        12 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE 
>>>  1  1  0.500000       0.00000       0.00000       0.00000 0.00000     m 
>>>  2  1   1.10000      0.500000       0.00000       0.00000 0.00000     m 
>>> . 
>>> . 
>>> . 
>>> 10  1   1.00000     -0.500000      0.500000       0.00000 0.00000     m 
>>> 11  1   1.00000       1.00000      0.500000       0.00000 0.00000     m 
>>> 12  1   0.00000       1.00000      0.500000       0.00000 0.00000     m 
>>> 
>>> where edge 10 is supposedly not a straight line any more. 
>>> 
>>> When I run this in nek, however, the cube still comes out with straight 
>>> edges. 
>>> 
>>> Are there any other parameters I need to set? 
>>> 
>>> I checked out the most recent nek version and overcame compilation issues 
>>> with gcc-gfortran 4.4.3-4.fc12 from the fedora 12 repository. 
>>> 
>>> Thanks, 
>>> Markus 
>>> 
>>> 
>>> nek5000-users at lists.mcs.anl.gov wrote: 
>>>> 
>>>> 
>>>> Michael, 
>>>> 
>>>> I've only recently added general purpose midside-node support, 
>>>> which puts a point in 3-space for any one of the 12 edges and 
>>>> nek then fits a parabola to this. 
>>>> 
>>>> I'll set up an example that demos this.  The feature has limited 
>>>> support at the moment -- but does work in nek5000 and generates 
>>>> correct geometry. 
>>>> 
>>>> Depending on what you are after, there may be other ways to 
>>>> generate the geometry.   One of my favorite techniques is to 
>>>> project a given input geometry onto the desired surface. 
>>>> 
>>>> 
>>>> Paul 
>>>> 
>>>> 
>>>> 
>>>> On Sat, 27 Mar 2010, nek5000-users at lists.mcs.anl.gov wrote: 
>>>> 
>>>> 
>>>> 
>>>>  Hello Users, 
>>>> 
>>>> 
>>>> 
>>>> 
>>>>  My question is in regard to the curved side data section in the REA 
>>>> file.  I know that the first 3 terms describe the side, element, and 
>>>> radius, but there are several other numbers that are currently zero, and 
>>>> then the letter C. 
>>>> 
>>>> 
>>>> 
>>>>  Background for the question: I have a situation where I would like to 
>>>> curve an element edge in any plane. I know this is possible from 
>>>> previous posts regarding a sphere.rea. 
>>>> 
>>>> 
>>>> 
>>>>  I know that with two points and a radius, that is enough to describe a 
>>>> circle in a plane.  But lets say I have a side where 2 pts lie on the 
>>>> XY plane, but the center of the circle is located on the YZ plane for 
>>>> example... the question is how to represent this in the REA. In the REA 
>>>> you just give the element side ( 2 pts ) and the radius, which does fix 
>>>> the circle center but the plane that contains the center is not 
>>>> fixed.  So what I am wondering is what the other numbers in this 
>>>> section do...if you could say, give the center instead of the radius, or 
>>>> give a third point with the other 2 pts in the element side to fully 
>>>> define the circle and plane. 
>>>> 
>>>> 
>>>> 
>>>>  Also, Is it possible to curve the remaining 4 sides of the element 
>>>> (edges 9, 10, 11, 12 that would be in the "z" direction) in the same 
>>>> manner as edges 1-8? Thanks for any input on this matter! 
>>>> 
>>>> 
>>>> 
>>>>  Regards, 
>>>> 
>>>>  Michael Meador 
>>>> 
>>>> 
>>>> ------------------------------------------------------------------------ 
>>>> 
>>>> _______________________________________________ 
>>>> 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 
>>> 
>> _______________________________________________ 
>> 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 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20100403/717e1b79/attachment.html>