[Gmsh] Transfinite interpolation cuts through the boundary.

Maka Mohu maka.mohu at gmail.com
Thu Jan 31 20:07:16 CET 2008


If you try to run the following example and then zoom in near the highest
point in the curved boundary you will see that the faces of the element cuts
through the boundary. Is there a way to control that. Thanks.

ex:

---------------------------------------------------------------------------------


lc = 0.1;

Point(1) = {0, 0, 0, lc};
Point(2) = {.1, 0,  0, lc} ;
Point(3) = {.1, .3, 0, lc} ;
Point(4) = {0,  .3, 0, lc} ;
Point(5) = {.05, .15, 0, lc};

Spline(1) = {1,5,2} ;
Line(2) = {3,2} ;
Line(3) = {3,4} ;
Line(4) = {4,1} ;


Line Loop(5) = {4,1,-2,3} ;

Plane Surface(6) = {5} ;


Transfinite Line{4} = 20 ;
Transfinite Line{2} = 20 ;//Using Bump 0.05;
Transfinite Line{-1,3} = 20 ;// Using Progression 1.2;
Transfinite Surface{6} = {1,2,3,4};

Mesh.ElementOrder = 1;

// Recombine the triangles into quads
Recombine Surface{6};

// Apply an elliptic smoother to the grid
Mesh.Smoothing = 100;
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20080131/81045ed5/attachment.html>