[Nek5000-users] Representing Curved Side in any plane

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Mon Apr 5 21:29:33 CDT 2010 

Hi,

I've just added the midside-node support to the current svn repo
for nek.

Sorry - I thought I had upgraded this in October but apparently
had not committed to the repo then.

Hopefully all should work -- I retested my original benchmark w/
this new version and it is functioning.

Please let me know if this now works for you.  You should be
able to modify any one of the 12 edges.

Paul


On Sat, 3 Apr 2010, nek5000-users at lists.mcs.anl.gov wrote:

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

More information about the Nek5000-users mailing list