Hi Christophe,
thanks for your quick reply.
The optimization steps (Mesh.Optimize=1 and Mesh.OptimizeNetgen) improve
the mesh, but there is still a sliver in my mesh (see attached picture).
I attached also my GMSH log file.
Is there a possibility to set a minimum element volume in GMSH?
Another interesting point is that the standard 3D GMSH Delauny algorithm
(Mesh.Algorithm3D = 1) does not finish 3D meshing. The tetgen step is
finished, but GMSH optimization step not.
Do you feel that the surface mesh is the problem? Is it of very bad quality?
Best Regards
Alex
On 10/29/2015 10:44 AM, Christophe Geuzaine wrote:
Mesh.Optimize=1;" and/or "Mesh.OptimizeNetgen=1
--
Alexander Tismer
Institut für Strömungsmechanik
und Hydraulische Strömungsmaschinen
Universität Stuttgart
Pfaffenwaldring 10
70550 Stuttgart
Tel.: +49 (0) 711 / 685 63007
Fax: +49 (0) 711 / 685 63255
eMail: [email protected]
http://www.ihs.uni-stuttgart.de
Info : -------------------------------------------------------
Info : Gmsh version : 2.10.1
Info : Build OS : Linux
Info : Build options : 64Bit Ann Bamg Bfgs Blas(Generic) Blossom Cairo
Chaco DIntegration Dlopen Fltk GMP Gmm Jpeg Kbipack Lapack(Generic)
LinuxJoystick MathEx Mesh Metis Mmg3d Mpeg NativeFileChooser Netgen ONELAB
ONELABMetamodel OpenGL OptHom Parser Plugins Png Post Solver Tetgen(1.5) Voro3D
Zlib
Info : Build date : 20151028
Info : Build host : ringkolbenschieber
Info : Packager : swadm
Info : Executable :
/mnt/opt.net/src/gmsh/Linux~13.2~x86_64/bin/gmsh79Cí
Info : Home directory : /mnt/fs2/home/atismer/
Info : Launch date : Thu Oct 29 11:27:12 2015
Info : Command line : gmsh gmsh.geo
Info : -------------------------------------------------------
Info : Reading 'gmsh.geo'...
Info : Reading 'transferToTetgen0.stl'...
Info : 29168 facets in solid 0
Info : Done reading 'transferToTetgen0.stl'
Info : Done reading 'gmsh.geo'
Info : Meshing 1D...
Info : Done meshing 1D (0 s)
Info : Meshing 2D...
Info : Done meshing 2D (5.88894e-05 s)
Info : Meshing 3D...
Info : Delaunay Meshing 1 volumes with 1 connected components
Info : Meshing volume 1 (Delaunay)
Info : 14584 points 73390 edges and 29168 faces in the initial mesh
Info : 14584 points created -- Worst tet radius is 751.249
Info : 19584 points created -- Worst tet radius is 2.29727
Info : 24584 points created -- Worst tet radius is 1.31693
Info : 29584 points created -- Worst tet radius is 1.02252
Info : 34584 points created -- Worst tet radius is 1.2026
Info : 39584 points created -- Worst tet radius is 1.6418
Info : Done meshing 3D (14.961 s)
Info : Optimizing 3D mesh...
Info : Optimizing volume 1
Info : Opti : START with 3.48125E-01 QBAD 1.70755E-04 QAVG 5.46624E-01
Info : Opti : 0.00 < QUAL < 0.10 : 5122 elements
Info : Opti : 0.10 < QUAL < 0.20 : 9143 elements
Info : Opti : 0.20 < QUAL < 0.30 : 9708 elements
Info : Opti : 0.30 < QUAL < 0.40 : 7863 elements
Info : Opti : 0.40 < QUAL < 0.50 : 8404 elements
Info : Opti : 0.50 < QUAL < 0.60 : 12706 elements
Info : Opti : 0.60 < QUAL < 0.70 : 24112 elements
Info : Opti : 0.70 < QUAL < 0.80 : 24827 elements
Info : Opti : 0.80 < QUAL < 0.90 : 10719 elements
Info : Opti : 0.90 < QUAL < 1.00 : 1579 elements
Info : Opti : (10924,9,5408) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.80184E-01 ( 3.875 sec)
Info : Opti : (13272,89,8217) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.85011E-01 ( 7.212 sec)
Info : Opti : (14277,129,9970) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.86735E-01 ( 10.455 sec)
Info : Opti : (14761,150,11091) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.87419E-01 ( 13.669 sec)
Info : Opti : (15071,167,11916) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.87851E-01 ( 16.874 sec)
Info : Opti : (15326,176,12537) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88269E-01 ( 20.069 sec)
Info : Opti : (15481,185,13109) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88439E-01 ( 23.257 sec)
Info : Opti : (15583,187,13623) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88537E-01 ( 26.432 sec)
Info : Opti : (15656,189,14057) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88693E-01 ( 29.594 sec)
Info : Opti : (15720,195,14332) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88806E-01 ( 32.759 sec)
Info : Opti : (15746,199,14516) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88832E-01 ( 35.928 sec)
Info : Opti : (15761,202,14635) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88833E-01 ( 39.082 sec)
Info : Opti : (15768,203,14746) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88822E-01 ( 42.302 sec)
Info : Opti : (15780,206,14812) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88829E-01 ( 45.491 sec)
Info : Opti : (15784,207,14891) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88820E-01 ( 48.652 sec)
Info : Opti : (15789,211,14961) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88814E-01 ( 51.926 sec)
Info : Opti : (15798,213,15000) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88811E-01 ( 55.074 sec)
Info : Opti : (15806,213,15038) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88814E-01 ( 58.196 sec)
Info : Opti : (15809,214,15060) = 3.48125E-01 QBAD 6.92043E-04 QAVG
5.88821E-01 ( 61.333 sec)
Info : Opti : 2679 ill-shaped tets are still in the mesh
Info : Opti : 0.00 < QUAL < 0.10 : 2679 elements
Info : Opti : 0.10 < QUAL < 0.20 : 7732 elements
Info : Opti : 0.20 < QUAL < 0.30 : 7402 elements
Info : Opti : 0.30 < QUAL < 0.40 : 4755 elements
Info : Opti : 0.40 < QUAL < 0.50 : 4667 elements
Info : Opti : 0.50 < QUAL < 0.60 : 13782 elements
Info : Opti : 0.60 < QUAL < 0.70 : 25193 elements
Info : Opti : 0.70 < QUAL < 0.80 : 26717 elements
Info : Opti : 0.80 < QUAL < 0.90 : 12674 elements
Info : Opti : 0.90 < QUAL < 1.00 : 1561 elements
Info : Done optimizing 3D mesh (64.975 s)
Info : Optimizing 3D mesh with Netgen...
Info : Optimizing volume 1
Info : CalcLocalH: 25521 Points 107162 Elements 29168 Surface Elements
Info : Remove Illegal Elements
Info : 2296 illegal tets
Info : SplitImprove
Info : badmax = 9.502e+06
Info : 373 splits performed
Info : SwapImprove
Info : 358 swaps performed
Info : SwapImprove2
Info : 16 swaps performed
Info : 1499 illegal tets
Info : SplitImprove
Info : badmax = 7.23561e+06
Info : 334 splits performed
Info : SwapImprove
Info : 184 swaps performed
Info : SwapImprove2
Info : 42 swaps performed
Info : 727 illegal tets
Info : SplitImprove
Info : badmax = 2.61504e+06
Info : 190 splits performed
Info : SwapImprove
Info : 129 swaps performed
Info : SwapImprove2
Info : 20 swaps performed
Info : 240 illegal tets
Info : SplitImprove
Info : badmax = 1.64169e+06
Info : 76 splits performed
Info : SwapImprove
Info : 57 swaps performed
Info : SwapImprove2
Info : 7 swaps performed
Info : 69 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 17 splits performed
Info : SwapImprove
Info : 21 swaps performed
Info : SwapImprove2
Info : 4 swaps performed
Info : 26 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 4 splits performed
Info : SwapImprove
Info : 11 swaps performed
Info : SwapImprove2
Info : 0 swaps performed
Info : 15 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 0 splits performed
Info : SwapImprove
Info : 2 swaps performed
Info : SwapImprove2
Info : 1 swaps performed
Info : 15 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 0 splits performed
Info : SwapImprove
Info : 2 swaps performed
Info : SwapImprove2
Info : 0 swaps performed
Info : 15 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 0 splits performed
Info : SwapImprove
Info : 2 swaps performed
Info : SwapImprove2
Info : 0 swaps performed
Info : 15 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 0 splits performed
Info : SwapImprove
Info : 2 swaps performed
Info : SwapImprove2
Info : 0 swaps performed
Info : 15 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 0 splits performed
Info : SwapImprove
Info : 0 swaps performed
Info : SwapImprove2
Info : 0 swaps performed
Info : 15 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 0 splits performed
Info : SwapImprove
Info : 0 swaps performed
Info : SwapImprove2
Info : 0 swaps performed
Info : 15 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 0 splits performed
Info : SwapImprove
Info : 0 swaps performed
Info : SwapImprove2
Info : 0 swaps performed
Info : 15 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 0 splits performed
Info : SwapImprove
Info : 0 swaps performed
Info : SwapImprove2
Info : 0 swaps performed
Info : 15 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 0 splits performed
Info : SwapImprove
Info : 0 swaps performed
Info : SwapImprove2
Info : 0 swaps performed
Info : 15 illegal tets
Info : SplitImprove
Info : badmax = 1.61339e+06
Info : 0 splits performed
Info : SwapImprove
Info : 0 swaps performed
Info : SwapImprove2
Info : 0 swaps performed
Info : 15 illegal tets
Info : Volume Optimization
Info : CombineImprove
Info : 1275 elements combined
Info : ImproveMesh
Info : Total badness = 388312
Info : Total badness = 345362
Info : SplitImprove
Info : badmax = 972989
Info : 46 splits performed
Info : ImproveMesh
Info : Total badness = 373083
Info : Total badness = 343857
Info : SwapImprove
Info : 3679 swaps performed
Info : SwapImprove2
Info : 38 swaps performed
Info : ImproveMesh
Info : Total badness = 320961
Info : Total badness = 302251
Info : CombineImprove
Info : 148 elements combined
Info : ImproveMesh
Info : Total badness = 300600
Info : Total badness = 299950
Info : SplitImprove
Info : badmax = 716164
Info : 12 splits performed
Info : ImproveMesh
Info : Total badness = 308568
Info : Total badness = 300803
Info : SwapImprove
Info : 1919 swaps performed
Info : SwapImprove2
Info : 36 swaps performed
Info : ImproveMesh
Info : Total badness = 284519
Info : Total badness = 276936
Info : CombineImprove
Info : 56 elements combined
Info : ImproveMesh
Info : Total badness = 276300
Info : Total badness = 276141
Info : SplitImprove
Info : badmax = 565802
Info : 6 splits performed
Info : ImproveMesh
Info : Total badness = 281251
Info : Total badness = 273455
Info : SwapImprove
Info : 1193 swaps performed
Info : SwapImprove2
Info : 41 swaps performed
Info : ImproveMesh
Info : Total badness = 263824
Info : Total badness = 261879
Info : Done optimizing 3D mesh with Netgen (55.074 s)
Info : 25100 vertices 132723 elements
Info : Writing '/mnt/scratch/atismer/tmp/26_netgenGmsh/gmsh.txt'...
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh
