[Nek5000-users] Vanishing Jacobian error description
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Wed Apr 21 10:02:16 CDT 2010
One other question -
Have you sorted out how the 12 edges are numbered?
On Wed, 21 Apr 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
> Michael,
>
> I looked at this for a couple of different N.
>
> nek/postx think that the geometry is messed up ---
>
> see www.mcs.anl.gov/~fischer/zzz.jpg
>
> 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, nek5000-users at lists.mcs.anl.gov 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: nek5000-users at lists.mcs.anl.gov
>> To: nek5000-users at lists.mcs.anl.gov
>> 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, nek5000-users at lists.mcs.anl.gov 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:
>> nek5000-users at lists.mcs.anl.gov > 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,
> nek5000-users at lists.mcs.anl.gov 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
>> 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