[Gmsh] crash on 3D meshing (AW: Anisiotropic mesh in stacked layers)

[Zenker, Dr. Matthias]

I forgot to say that the crash happens under Windows 7 with gmsh 2.8.0 64bit and with 2.7.1 32bit and 64bit.

Matthias

      _____________________________________________
     Von: Zenker, Dr. Matthias
     Gesendet: Donnerstag, 11. Juli 2013 16:57
     An: 'gmsh at geuz.org'
     Cc: 'Christophe Geuzaine [cgeuzaine at ulg.ac.be] (cgeuzaine at ulg.ac.be)'
     Betreff: crash on 3D meshing (AW: Anisiotropic mesh in stacked layers)


     Hi,

     now that I have successfully created a anisotropic mesh width field for my geometry, I get a crash when I do the 3D meshing. 2D meshing works without problem. I have the impression that the desired mesh width may be too small. Is this so, and what is the smallest mesh width gmsh can cope with?

     Thanks for an answer, this time...

     Matthias

      < Datei: layers.geo >>
           _____________________________________________
           Von: Zenker, Dr. Matthias
           Gesendet: Donnerstag, 11. Juli 2013 16:41
           An: 'gmsh at geuz.org'
           Cc: 'Christophe Geuzaine [cgeuzaine at ulg.ac.be] (cgeuzaine at ulg.ac.be)'
           Betreff: AW: Anisiotropic mesh in stacked layers


           Hi again,

           I have got it to work, but I still don't understand why. It seems that the fields in MathEvalAniso have to be 1/(mesh_width)^2. But why?
           I have read the relevant part of the gmsh manual, but didn't find the answer there.
           A short explanation and/or a pointer to more info would be highly appreciated.

           Matthias

                 _____________________________________________
                 Von: Zenker, Dr. Matthias
                 Gesendet: Mittwoch, 10. Juli 2013 17:09
                 An: gmsh at geuz.org<mailto:gmsh at geuz.org>
                 Cc: Christophe Geuzaine [cgeuzaine at ulg.ac.be] (cgeuzaine at ulg.ac.be<mailto:cgeuzaine at ulg.ac.be>)
                 Betreff: Anisiotropic mesh in stacked layers


                 Hi,

                 I want to mesh a geometry with layers stacked in z direction, where each layer should have a mesh width in the xy plane, and a different one in z direction. I have tried to combine the Box, Min and MathEval Aniso fields to achieve this (see attachment). But the mesh widths are not as I expect, in particular in layer 2 in z direction it should be 1/20 its height. What do I do wrong?
                 How is the mesh width derived from the m11, m22 and m33 fields in MathEvalAniso?

                 Thanks for some explanations and hints,

                 Matthias

                  < Datei: layers.geo >>



_____________________________________________________________________
ERBE Elektromedizin GmbH
Firmensitz: 72072 Tuebingen
Geschaeftsfuehrer: Christian O. Erbe, Reiner Thede
Registergericht: Stuttgart HRB 380137

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20130711/111d49f0/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: layers.geo
Type: application/octet-stream
Size: 3424 bytes
Desc: layers.geo
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20130711/111d49f0/attachment.geo>