<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: Times New Roman; font-size: 12pt; color: #000000'><P>Hi Paul,</P>
<P> </P>
<P>Were you able to figure out the issue with Markus's midpoint case?  Thank you for looking into this!</P>
<P> </P>
<P>Kindly,</P>
<P> </P>
<P>Michael</P>
<P><BR>----- Original Message -----<BR>From: nek5000-users@lists.mcs.anl.gov<BR>To: nek5000-users@lists.mcs.anl.gov<BR>Sent: Monday, March 29, 2010 11:37:43 AM GMT -06:00 US/Canada Central<BR>Subject: Re: [Nek5000-users] Representing Curved Side in any plane<BR><BR><BR>Hi Markus,<BR><BR>Thanks -- I'll check into it.<BR><BR>Paul<BR><BR><BR>On Mon, 29 Mar 2010, nek5000-users@lists.mcs.anl.gov wrote:<BR><BR>> Hi,<BR>><BR>> just double-checked, that's what I did. I only ran it for one time step without<BR>> any meaningful physics, though.<BR>> Attached are all case files. The nek version I am using is revision 456.<BR>><BR>> Markus<BR>><BR>><BR>><BR>> Quoting nek5000-users@lists.mcs.anl.gov:<BR>><BR>>><BR>>> Hi Markus,<BR>>><BR>>> Did you visualize this with VisIt and with the geometry<BR>>> put out into (at least) the first .fld or .f file ?<BR>>><BR>>> Paul<BR>>><BR>>><BR>>> On Mon, 29 Mar 2010, nek5000-users@lists.mcs.anl.gov wrote:<BR>>><BR>>>> Hi,<BR>>>><BR>>>> I was looking into the midpoint feature and generated a cube (1 element, x,<BR>>>> y, z from 0 to 1) in prenek, then converted all edges to the midpoint<BR>>>> notation with prenek, and then manipulated one edge in the rea file. This<BR>>> is<BR>>>> the resulting curved side section:<BR>>>>  ***** CURVED SIDE DATA *****<BR>>>>        12 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE<BR>>>>  1  1  0.500000       0.00000       0.00000       0.00000 0.00000     m<BR>>>>  2  1   1.10000      0.500000       0.00000       0.00000 0.00000     m<BR>>>> .<BR>>>> .<BR>>>> .<BR>>>> 10  1   1.00000     -0.500000      0.500000       0.00000 0.00000     m<BR>>>> 11  1   1.00000       1.00000      0.500000       0.00000 0.00000     m<BR>>>> 12  1   0.00000       1.00000      0.500000       0.00000 0.00000     m<BR>>>><BR>>>> where edge 10 is supposedly not a straight line any more.<BR>>>><BR>>>> When I run this in nek, however, the cube still comes out with straight<BR>>>> edges.<BR>>>><BR>>>> Are there any other parameters I need to set?<BR>>>><BR>>>> I checked out the most recent nek version and overcame compilation issues<BR>>>> with gcc-gfortran 4.4.3-4.fc12 from the fedora 12 repository.<BR>>>><BR>>>> Thanks,<BR>>>> Markus<BR>>>><BR>>>><BR>>>> nek5000-users@lists.mcs.anl.gov wrote:<BR>>>>><BR>>>>><BR>>>>> Michael,<BR>>>>><BR>>>>> I've only recently added general purpose midside-node support,<BR>>>>> which puts a point in 3-space for any one of the 12 edges and<BR>>>>> nek then fits a parabola to this.<BR>>>>><BR>>>>> I'll set up an example that demos this.  The feature has limited<BR>>>>> support at the moment -- but does work in nek5000 and generates<BR>>>>> correct geometry.<BR>>>>><BR>>>>> Depending on what you are after, there may be other ways to<BR>>>>> generate the geometry.   One of my favorite techniques is to<BR>>>>> project a given input geometry onto the desired surface.<BR>>>>><BR>>>>><BR>>>>> Paul<BR>>>>><BR>>>>><BR>>>>><BR>>>>> On Sat, 27 Mar 2010, nek5000-users@lists.mcs.anl.gov wrote:<BR>>>>><BR>>>>><BR>>>>><BR>>>>>  Hello Users,<BR>>>>><BR>>>>><BR>>>>><BR>>>>><BR>>>>>  My question is in regard to the curved side data section in the REA<BR>>>>> file.  I know that the first 3 terms describe the side, element, and<BR>>>>> radius, but there are several other numbers that are currently zero, and<BR>>>>> then the letter C.<BR>>>>><BR>>>>><BR>>>>><BR>>>>>  Background for the question: I have a situation where I would like to<BR>>>>> curve an element edge in any plane. I know this is possible from<BR>>>>> previous posts regarding a sphere.rea.<BR>>>>><BR>>>>><BR>>>>><BR>>>>>  I know that with two points and a radius, that is enough to describe a<BR>>>>> circle in a plane.  But lets say I have a side where 2 pts lie on the<BR>>>>> XY plane, but the center of the circle is located on the YZ plane for<BR>>>>> example... the question is how to represent this in the REA. In the REA<BR>>>>> you just give the element side ( 2 pts ) and the radius, which does fix<BR>>>>> the circle center but the plane that contains the center is not<BR>>>>> fixed.  So what I am wondering is what the other numbers in this<BR>>>>> section do...if you could say, give the center instead of the radius, or<BR>>>>> give a third point with the other 2 pts in the element side to fully<BR>>>>> define the circle and plane.<BR>>>>><BR>>>>><BR>>>>><BR>>>>>  Also, Is it possible to curve the remaining 4 sides of the element<BR>>>>> (edges 9, 10, 11, 12 that would be in the "z" direction) in the same<BR>>>>> manner as edges 1-8? Thanks for any input on this matter!<BR>>>>><BR>>>>><BR>>>>><BR>>>>>  Regards,<BR>>>>><BR>>>>>  Michael Meador<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>>>> Nek5000-users mailing list<BR>>>> Nek5000-users@lists.mcs.anl.gov<BR>>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<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>><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></P></div></body></html>