<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>Hey Stefan,</P>
<P> </P>
<P>I do not know C, and I am still learning the in & outs of Fortran. I certainly think it is possible though. I do use some built in functions that matlab offers to make my life easier, but I'm sure they can be replaced with if needed. </P>
<P> </P>
<P>Michael B,</P>
<P> </P>
<P>I'm still testing various geometries, becoming more complex each time. Our research group is using Nek with applications to Gas turbines, so our motivation to pursue this was indeed because of complex geometries. From Solidworks, you can export as a "ParaSolid". You can then Import this ParaSolid into Gambit. From Gambit, I export a generic Mesh which is a .NEU (neutral file). This is what the matLab code uses, the .neu file. There are many middle steps, but the process seems to work.</P>
<P> </P>
<P>I think there is plenty of room for improvement, but at the moment I'm still testing new geometries.</P>
<P> </P>
<P>- Michael M.</P>
<P><BR>----- Original Message -----<BR>From: nek5000-users@lists.mcs.anl.gov<BR>To: nek5000-users@lists.mcs.anl.gov<BR>Sent: Wednesday, April 21, 2010 10:52:48 AM GMT -06:00 US/Canada Central<BR>Subject: Re: [Nek5000-users] Vanishing Jacobian error description<BR><BR>Great!</P>
<DIV>Any chance to change your MATLAB code into C or Fortran?</DIV>
<DIV><BR></DIV>
<DIV>Stefan</DIV>
<DIV><BR></DIV>
<DIV> </DIV>
<DIV><BR>
<DIV>
<DIV>On Apr 21, 2010, at 5:48 PM, <A href="mailto:nek5000-users@lists.mcs.anl.gov" target=_blank>nek5000-users@lists.mcs.anl.gov</A> wrote:</DIV><BR class=Apple-interchange-newline>
<BLOCKQUOTE><SPAN style="WIDOWS: 2; TEXT-TRANSFORM: none; TEXT-INDENT: 0px; BORDER-COLLAPSE: separate; FONT: medium Helvetica; WHITE-SPACE: normal; ORPHANS: 2; LETTER-SPACING: normal; WORD-SPACING: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px" class=Apple-style-span>
<DIV>
<DIV style="FONT-FAMILY: 'Times New Roman'; COLOR: rgb(0,0,0); FONT-SIZE: 12pt">Paul,<BR><BR>I figured out the issue, and it was a mistake in the routine I was using to put the elements midside nodes in the correct locations. I fixed the issue and you can see the attached result case and image if you (or anyone) would like.<BR><BR>I should note that this geometry, although simple, was created by:<BR><BR>1. Modeled in Solidworks <SPAN class=Apple-converted-space> </SPAN><BR>2. Imported into Gambit<BR>3. Meshed in Gambit<BR>4. Run through our matLab code to generate rea file<BR><BR>This is the result, a proof of concept. If you think there would be significant interest I can post a separate email to the user list, so it is not buried in this "vanishing Jacobian" topic.<BR><BR>Regards,<BR>Michael<BR><BR>----- Original Message -----<BR>From:<SPAN class=Apple-converted-space> </SPAN><A href="mailto:nek5000-users@lists.mcs.anl.gov" target=_blank>nek5000-users@lists.mcs.anl.gov</A><BR>To:<SPAN class=Apple-converted-space> </SPAN><A href="mailto:nek5000-users@lists.mcs.anl.gov" target=_blank>nek5000-users@lists.mcs.anl.gov</A><BR>Sent: Wednesday, April 21, 2010 10:07:52 AM GMT -06:00 Guadalajara / Mexico City / Monterrey<BR>Subject: Re: [Nek5000-users] Vanishing Jacobian error description<BR><BR>
<DIV style="FONT-FAMILY: 'Times New Roman'; COLOR: rgb(0,0,0); FONT-SIZE: 12pt">Hi Paul,<BR><BR>Yes I am checking it right now. The midside nodes on the top layer and bottom layer are working, however the elements in the layers between are not. They are all generated the same way, so at the moment I'm not sure.<SPAN class=Apple-converted-space> </SPAN><BR><BR>When I remove the curved side section in the rea, it works both for the fine and coarse lofts. I'm not sure why it doesn't like the midpoints; it's almost as if it is placing them in the wrong place and fitting a parabola that interferes with surrounding geometry. I'll look into it more here and post again when I have an update.<BR> <BR>Thanks Paul<BR><BR>- Michael<BR><BR>----- Original Message -----<BR>From:<SPAN class=Apple-converted-space> </SPAN><A href="mailto:nek5000-users@lists.mcs.anl.gov" target=_blank>nek5000-users@lists.mcs.anl.gov</A><BR>To:<SPAN class=Apple-converted-space> </SPAN><A href="mailto:nek5000-users@lists.mcs.anl.gov" target=_blank>nek5000-users@lists.mcs.anl.gov</A><BR>Sent: Wednesday, April 21, 2010 10:00:29 AM GMT -06:00 Guadalajara / Mexico City / Monterrey<BR>Subject: Re: [Nek5000-users] Vanishing Jacobian error description<BR><BR>Michael, I looked at this for a couple of different N. nek/postx think that the geometry is messed up --- see<SPAN class=Apple-converted-space> </SPAN><A href="http://www.mcs.anl.gov/~fischer/zzz.jpg" target=_blank>www.mcs.anl.gov/~fischer/zzz.jpg</A><SPAN class=Apple-converted-space> </SPAN>That being said, it could easily be a problem w/ the midside node support that I added last month. (NOTE: to see the geometry I had to fix something in the nek code -- normally when you get vanishing jacobian the code will output "xyzblah.fld01" with the geometry, so that you can then look at this in postx or visit. For pn/pn, this was broken but is now working in the just-updated repo.) Do you have a way to verify that your lofting is working? Paul On Wed, 21 Apr 2010,<SPAN class=Apple-converted-space> </SPAN><A href="mailto:nek5000-users@lists.mcs.anl.gov" target=_blank>nek5000-users@lists.mcs.anl.gov</A><SPAN class=Apple-converted-space> </SPAN>wrote: > Hi, > > Sorry, yes that was for without the curve side data. Here is the case that includes the curved side data. > > - Michael > > ----- Original Message ----- > From:<SPAN class=Apple-converted-space> </SPAN><A href="mailto:nek5000-users@lists.mcs.anl.gov" target=_blank>nek5000-users@lists.mcs.anl.gov</A><SPAN class=Apple-converted-space> </SPAN>> To:<SPAN class=Apple-converted-space> </SPAN><A href="mailto:nek5000-users@lists.mcs.anl.gov" target=_blank>nek5000-users@lists.mcs.anl.gov</A><SPAN class=Apple-converted-space> </SPAN>> Sent: Wednesday, April 21, 2010 8:39:38 AM GMT -06:00 Guadalajara / Mexico City / Monterrey > Subject: Re: [Nek5000-users] Vanishing Jacobian error description > > The case you sent runs fine.... is it w/ or w/o the curved side info? Paul On Wed, 21 Apr 2010,<A href="mailto:nek5000-users@lists.mcs.anl.gov" target=_blank>nek5000-users@lists.mcs.anl.gov</A><SPAN class=Apple-converted-space> </SPAN>wrote: > > > Hi Paul, > > > > Here is the rea, usr, & SIZE file I am using. I ran it again without the curved side information and it worked, so I think the problem is somewhere in the curved side section of the rea. > > > > The geometry is a square base lofted into a circle in the +z direction. The mesh was generated in gambit. I've developed a matlab code that will take a mesh created in gambit and generate the required REA formatting. I have created some basic shapes to test the matlab portion and the curved side portion, and it worked in most cases (not this one yet). I'm not sure what is going on with the curved side data in this case however. > > > > Thanks for looking at this, > > Michael > > > ----- Original Message ----- > From:<SPAN class=Apple-converted-space> </SPAN><A href="mailto:nek5000-users@lists.mcs.anl.gov" target=_blank>nek5000-users@lists.mcs.anl.gov</A><SPAN class=Apple-converted-space> </SPAN>> To: "Nekton User List" > Sent: Tuesday, April 20, 2010 3:36:44 PM GMT -06:00 U S/Canada Central > Subject: Re: [Nek5000-users] Vanishing Jacobian error description > > > > Michael, > > The code is trying to check to ensure that the mesh is > valid. Sometimes, however, there are other causes for > the problem. > > Can you send a gzip'd rea/usr file, and I can take a look > later tonight. > > Paul > > > > On Tue, 20 Apr 2010,<SPAN class=Apple-converted-space> </SPAN><A href="mailto:nek5000-users@lists.mcs.anl.gov" target=_blank>nek5000-users@lists.mcs.anl.gov</A><SPAN class=Apple-converted-space> </SPAN>wrote: > >> Hi, >> >> I am running into an issue with 'Vanishing Jacobian'. To help understand what the problem >> is, could someone explain what is meant in the error? Specifically, what it means by "Xth node" of the element. >> Also, I ran into this issue when I was deforming geometries in the usr file, where one of the elements >> would overlap another by mistake. However in this case I know there are no overlaps. >> >> Here is a portion: >> >> 0 ERROR: Vanishing Jacobian near 91th node of element 69. >> 1.40692804605498586E-002 -7.79619559165757245E-003 >> 0 xyz: 2.65874E-01 3.54498E-01 2.87050E+0 0 >> 0 xyz: 2.07855E-01 1.77249E-01 2.87050E+00 >> >> >> 0 ERROR: Vanishing Jacobian near 76th node of element 70. >> 2.59472856644434001E-002 -2.12442422718029170E-002 >> 0 xyz: 2.11655E-01 4.23310E-01 2.62500E+00 >> 0 xyz: 3.54498E-01 1.77249E-01 2.87050E+00 >> >> Thanks for any help deciphering this error! >> >> - Michael >> >> >> >> >> >> > _______________________________________________ > Nek5000-users mailing list ><A href="mailto:Nek5000-users@lists.mcs.anl.gov" target=_blank>Nek5000-users@lists.mcs.anl.gov</A><SPAN class=Apple-converted-space> </SPAN>><SPAN class=Apple-converted-space> </SPAN><A href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target=_blank>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</A><SPAN class=Apple-converted-space> </SPAN>> > _______________________________________________ Nek5000-users mailing list<SPAN class=Apple-converted-space> </SPAN><A href="mailto:Nek5000-users@lists.mcs.anl.gov" target=_blank>Nek5000-users@lists.mcs.anl.gov</A><SPAN class=Apple-converted-space> </SPAN><A href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target=_blank>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</A><BR>_______________________________________________ Nek5000-users mailing list<SPAN class=Apple-converted-space> </SPAN><A href="mailto:Nek5000-users@lists.mcs.anl.gov" target=_blank>Nek5000-users@lists.mcs.anl.gov</A><SPAN class=Apple-converted-space> </SPAN><A href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target=_blank>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</A></DIV><BR>_______________________________________________ Nek5000-users mailing list<SPAN class=Apple-converted-space> </SPAN><A href="mailto:Nek5000-users@lists.mcs.anl.gov" target=_blank>Nek5000-users@lists.mcs.anl.gov</A><SPAN class=Apple-converted-space> </SPAN><A href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target=_blank>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</A></DIV><SPAN><test4.tar></SPAN><SPAN><visit0001.tif></SPAN>_______________________________________________<BR>Nek5000-users mailing list<BR><A href="mailto:Nek5000-users@lists.mcs.anl.gov" target=_blank>Nek5000-users@lists.mcs.anl.gov</A><BR><A href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target=_blank>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</A><BR></DIV></SPAN></BLOCKQUOTE></DIV><BR></DIV><BR>_______________________________________________ Nek5000-users mailing list Nek5000-users@lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users </div></body></html>