[Nek5000-users] Vanishing Jacobians while using moving mesh

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Wed Feb 15 08:39:10 CST 2012


Pradeep,

Which version of Visit are you using?  We have experienced similar
issues with version 2.4 when viewing 2D meshes.  If this is the case,
try using 2.3.2 and see if you have the same problem.

Josh

On Tue, Feb 14, 2012 at 8:08 AM,  <nek5000-users at lists.mcs.anl.gov> wrote:
> Hi Neks,
>
> I'm trying to simulate 2D flow over an oscillating cylinder using moving
> mesh, and tried to set up a test case by looking at the peris example. I
> however am getting vanishing Jacobians after the first few timesteps.
>
> I am modifying the the boundary velocity of the cylinder in userbc instead
> of userchk (as in the peris example). Would it be possible for you to look
> at my .usr file and tell me what I am doing wrong?
>
> I have attached the .rea, .usr, and SIZE with this mail.
>
> Further, when I try to view the mesh of the output (or the xyz fld
> outputfile) in visit, visit crashes with the following error:
> VisIt's viewer exited abnormally! Aborting the Graphical User Interface.
>
> Lastly, would it be possible to simulate moving walls with large
> deformations (i.e. where multiple elements get deformed), where the
> deformation is known ahead of time? (eg. a flapping aerofoil)
>
> Thank you for your help!
>
> Regards,
> Pradeep
>
> Extract from logfile:
>
>         60 U-PRES gmres:    100  7.5345E+42  1.0000E-02  4.0062E+57
> 6.7556E-03  1.5093E-02
>  ERROR:  alphad .le. 0 in econjp -1.31175831444004495E+179          16
>  ERROR:  alphad .le. 0 in econjp -1.31175831444004495E+179          16
>  ERROR:  alphad .le. 0 in econjp -1.31175831444004495E+179          16
>  ERROR:  alphad .le. 0 in econjp -1.31175831444004495E+179          16
>           60  DNORM, DIVEX  1.18120867781751263E+043
> 7.53448962619089522E+042
>          60   6.0000E-02  1.6697E-02 Fluid done
>      DT/DTCFL/DTFS/DTINIT   0.635E-45   0.635E-45   0.000E+00   0.100E-02
> Step     61, t= 6.0800000E-02, DT= 8.0000000E-04, C=******* 2.2667E-01
> 1.8287E-02
>              Solving for fluid
>
>
>     3  ERROR:  Vanishing Jacobian near      2th node of element         3.
>
>
>     2  ERROR:  Vanishing Jacobian near      2th node of element         1.
>   2.27348348709230109E+081 -5.12478988977734591E+080
>  -2.60483379627720531E+079  3.21465902039155941E+039
>     3 xyz:  -5.30056E+39  -5.84642E+39
>     3 xyz:  -5.24850E+39  -6.93109E+39
>
>
>     3  ERROR:  Vanishing Jacobian near      2th node of element         4.
>     2 xyz:   1.25061E+38   4.64957E+39
>   9.79493533554062390E+080 -3.01548938509156252E+080
>     2 xyz:  -2.46589E-08  -5.64092E-01
>
>
>     3 xyz:  -5.30056E+39  -5.84642E+39
>     3 xyz:  -5.24850E+39  -6.93109E+39
>     2  ERROR:  Vanishing Jacobian near      3th node of element         2.
>
>
>     3  ERROR:  Vanishing Jacobian near      2th node of element         5.
>   1.43125817410005932E+039 -6.03607926692375119E+038
>   5.54918579647991880E+039 -1.67302572815141788E+039
>     3 xyz:  -5.30056E+39  -5.84642E+39
>     2 xyz:  -2.46589E-08  -5.64092E-01
>     3 xyz:  -5.24850E+39  -6.93109E+39
>
>
>     3  ERROR:  Vanishing Jacobian near      2th node of element         6.
>     2 xyz:  -3.07794E-08  -7.04112E-01
>
>
>     2  ERROR:  Vanishing Jacobian near      7th node of element         7.
>   8.35923896949543198E+038 -3.55455289325888722E+038
>     3 xyz:  -5.30056E+39  -5.84642E+39
>   1.21747337192549314E+039 -1.79454406120891328E+038
>     3 xyz:  -5.24850E+39  -6.93109E+39
>     2 xyz:  -5.66434E-08  -1.29581E+00
>     2 xyz:  -6.27639E-08  -1.43583E+00
>
>
>     2  ERROR:  Vanishing Jacobian near      9th node of element         8.
>   1.17802538305091326E+039 -4.21489851051369387E+038
>     2 xyz:  -6.55671E-08  -1.49996E+00
>     2 xyz:  -3.25016E+38   4.93026E+39
>  call outfld: ifpsco: F
>
>        61  6.0800E-02 Write checkpoint:
>
>        61  6.0800E-02 OPEN:
> xyzcyl2.fld01
>  Jac error 1, setting p66=4, ifxyo=t
>
> call exitt: dying ...
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>



-- 
Josh Camp

"All that is necessary for the triumph of evil is that good men do
nothing" -- Edmund Burke



More information about the Nek5000-users mailing list