[Gmsh] Trouble in 3d meshing from a surface in both directions with Netgen

Hans Minnebo hans.minnebo at cenaero.be
Thu Aug 28 16:52:15 CEST 2008


Hello,

I have found the problem.
The transfinite surface 28 has to be oriented "Right" to solve the problem.
I tried to change the orientation of this surface before trying the 
"Right" orientation, but with no results.

Yours,

Hans MINNEBO

Hans Minnebo wrote:
> Hello,
>
> I am trying to mesh a cube, with a refined transfinite volume inside, 
> I have four interface surfaces.
> When using Netgen, I have the nodal continuity on the inner surfaces, 
> but for one (n°28) I do not have the continuity of the elements.
> Have you already encountered that kind of problem ? I tried versions 
> 2.1.1 and 2.2.4.
>
> I join the geo, msh, picture and options.
>
> Yours,
>
> Hans MINNEBO
> ------------------------------------------------------------------------
>
> // Gmsh Option File
> //
> // This file contains configuration options (preferences) that
> // are loaded each time Gmsh is launched. You can create this
> // file by hand, or let Gmsh generate it for you (with
> // 'Tools->Options->Save as defaults'). This file can also be
> // automatically saved every time you quit Gmsh if the option
> // 'General.SaveOptions' is set.
> //
> General.MouseHoverMeshes = 1; // Enable mouse hover on meshes
> General.Color.BackgroundGradient = {255,255,255}; // Background gradient color
> Geometry.LineType = 1; // Display lines as solid color segments (0), 3D cylinders (1) or tapered cylinders (2)
> Geometry.LineWidth = 1; // Display width of lines (in pixels)
> Geometry.PointSelectSize = 4.1; // Display size of selected points (in pixels)
> Geometry.PointSize = 2.8; // Display size of points (in pixels)
> Geometry.Surfaces = 1; // Display geometry surfaces?
> Geometry.SurfaceNumbers = 1; // Display surface numbers?
> Mesh.Algorithm3D = 4; // 3D mesh algorithm (1=Tetgen+Delaunay, 4=Netgen)
> Mesh.SurfaceFaces = 1; // Display faces of surface mesh?
> Mesh.VolumeFaces = 1; // Display faces of volume mesh?
> View.DisplacementFactor = 0.018150461; // Displacement amplification
> View.IntervalsType = 3; // Type of interval display (1=iso, 2=continuous, 3=discrete, 4=numeric)
> View.NbIso = 2; // Number of intervals
> View.RangeType = 2; // Value scale range type (1=default, 2=custom, 3=per time step)
> View.ShowScale = 0; // Show value scale?
> View.VectorType = 5; // Vector display type (1=segment, 2=arrow, 3=pyramid, 4=3D arrow, 5=displacement)
>   
>
> ------------------------------------------------------------------------
>
> ------------------------------------------------------------------------
>
> t1=0.2;
> Point(1) = {-1,-1,-1,t1};
> Point(2) = {-1,-1,1,t1};
> Point(3) = {-1,1,1,t1};
> Point(4) = {-1,1,-1,t1};
> Point(5) = {1,-1,-1,t1};
> Point(6) = {1,-1,1,t1};
> Point(7) = {1,1,1,t1};
> Point(8) = {1,1,-1,t1};
> yref=0.2;
> href=0.4;
> t2=0.1;
> Point(9) = {-1,-href,-href,t2};
> Point(10) = {-1,-href,href,t2};
> Point(11) = {-1,href+yref,href,t2};
> Point(12) = {-1,href+yref,-href,t2};
> Point(13) = {1,-href,-href,t2};
> Point(14) = {1,-href,href,t2};
> Point(15) = {1,href+yref,href,t2};
> Point(16) = {1,href+yref,-href,t2};
> Line(1) = {6,7};
> Line(2) = {7,8};
> Line(3) = {8,5};
> Line(4) = {5,6};
> Line(5) = {2,3};
> Line(6) = {3,4};
> Line(7) = {4,8};
> Line(8) = {5,1};
> Line(9) = {1,4};
> Line(10) = {1,2};
> Line(11) = {6,2};
> Line(12) = {7,3};
> Line(13) = {14,10};
> Line(14) = {15,11};
> Line(15) = {16,12};
> Line(16) = {13,9};
> Line(17) = {13,14};
> Line(18) = {14,15};
> Line(19) = {15,16};
> Line(20) = {16,13};
> Line(21) = {9,10};
> Line(22) = {10,11};
> Line(23) = {11,12};
> Line(24) = {12,9};
> Transfinite Line {17,19,21,23} = 10 Using Bump 2.0 ;//unchanged
> Transfinite Line {18,20,22,24} = 10 Using Bump 3.0 ;
> Transfinite Line {13,14,15,16} = 10 Using Progression 1.0 ;
>
> Line Loop(25) = {18,14,-22,-13};
> Plane Surface(26) = {25};
> Transfinite Surface {26} = {10,11,14,15} ;
> Line Loop(27) = {19,15,-23,-14};
> Plane Surface(28) = {27};
> Transfinite Surface {28} = {11,12,15,16} ;
> Line Loop(29) = {21,-13,-17,16};
> Plane Surface(30) = {29};
> Transfinite Surface {30} = {9,10,13,14} ;
> Line Loop(31) = {24,21,22,23};
> Plane Surface(32) = {31};
> Transfinite Surface {32} = {9,10,11,12} ;
> Line Loop(33) = {19,20,17,18};
> Plane Surface(34) = {33};
> Transfinite Surface {34} = {13,14,15,16} ;
> Line Loop(35) = {15,24,-16,-20};
> Plane Surface(36) = {35};
> Transfinite Surface {36} = {9,12,13,16} ;
> Line Loop(37) = {1,12,-5,-11};
> Plane Surface(38) = {37};
> Line Loop(39) = {10,-11,-4,8};
> Plane Surface(40) = {39};
> Line Loop(41) = {3,8,9,7};
> Plane Surface(42) = {41};
> Line Loop(43) = {2,-7,-6,-12};
> Plane Surface(44) = {43};
> Line Loop(45) = {3,4,1,2};
> Plane Surface(46) = {45,33};
> Line Loop(47) = {6,-9,10,5};
> Plane Surface(48) = {47,31};
> Physical Point(51) = {4};
> Physical Point(52) = {8};
> Physical Point(53) = {3};
> Physical Surface(1) = {38};
> Physical Surface(2) = {42};
> Surface Loop(56) = {30,26,28,36,32,34};
> Volume(57) = {56};
> Transfinite Volume {57} = {9,10,11,12,13,14,15,16};
> Surface Loop(58) = {48,44,46,42,40,38,26,28,36,30};
> Volume(59) = {58};
> Physical Volume(3) = {57,59};
>