[Gmsh] Transfinite Volume command - Example of use -

Paolo Tricerri paolotricerri at gmail.com
Tue Oct 23 13:56:36 CEST 2012


Dear all,

I am having problem at using the command Transfinite Volume. I am
trying to give to the definition of the volume the list of the eight
points at the corners of the geometry but the volume does not get
meshed.

Hereafter you can see the code I have using. Btw, I would like to ask
if any of you could give an example of use of the command Transfinite
Volume {expression-list} | "*" < = {expression-list} since I do not
understand if I have to pass to is Physical Point or Points.

Many thanks,

I wish you all a nice day,

Paolo Tricerri

Code:

lc = 0.1;
H= 0.44;
L = 0.48;
pH = 0.16;
thickness = pH/10;

// ############# POINTS ##############
//Base
Point(0) = {0, 0, 0, lc};
Point(1) = {0, H, 0, lc};
Point(2) = {L, H+pH, 0, lc};
Point(3) = {L, H, 0, lc};

Point(4) = {0, 0, thickness, lc};
Point(5) = {0, H, thickness, lc};
Point(6) = {L, H+pH, thickness, lc};
Point(7) = {L, H, thickness, lc};

Physical Point(1000) = {0};
Physical Point(2000) = {1};
Physical Point(3000) = {2};
Physical Point(4000) = {3};
Physical Point(5000) = {4};
Physical Point(6000) = {5};
Physical Point(7000) = {6};
Physical Point(8000) = {7};


//############### LINES ##############
//Base`
Line (1) = {0,1};
Line (2) = {1,2};
Line (3) = {2,3};
Line (4) = {3,0};

Line (5) = {4,5};
Line (6) = {5,6};
Line (7) = {6,7};
Line (8) = {7,4};

Line (9) = {0,4};
Line (10) = {1,5};
Line (11) = {3,7};
Line (12) = {2,6};

// ############# LINE LOOPS ############
Line Loop(1) = {1,2,3,4}; //on lines
Line Loop(2) = {5,6,7,8}; //on lines
Line Loop(3) = {1,10,-5,-9}; //on lines
Line Loop(4) = {3,11,-7,-12}; //on lines
Line Loop(5) = {9,-8,-11,4}; //on lines
Line Loop(6) = {10,6,-12,-2}; //on lines

// ############ SURFACES ############

Plane Surface(1) = {1}; //on line loops
Plane Surface(2) = {2}; //on line loops
Plane Surface(3) = {3}; //on line loops
Plane Surface(4) = {4}; //on line loops
Plane Surface(5) = {5}; //on line loops
Plane Surface(6) = {6}; //on line loops

// ############ TRANSFINITE LINES ############
Transfinite Line{1,3,5,7}=23; //Number of points on a 1/4 of circle
Transfinite Line{2,4,6,8}=15;       //Number of structure layers
Transfinite Line{9,10,11,12}=4;
// ############ TRANSFINITE SURFACES ############

Transfinite Surface { 1 } = { 0,1,2,3} ; //on points
Transfinite Surface { 2 } = { 4,5,6,7} ; //on points
Transfinite Surface { 3 } = { 0,4,5,1} ; //on points
Transfinite Surface { 4 } = { 3,7,6,2} ; //on points
Transfinite Surface { 5 } = { 0,4,7,3} ; //on points
Transfinite Surface { 6 } = { 1,5,6,2} ; //on points
Physical Surface(20)={1};
Physical Surface(30)={2};
Physical Surface(40)={3};
Physical Surface(50)={4};
Physical Surface(60)={5};
Physical Surface(70)={6};

//Surface Loop(100) = {0,1,2,3,4,5,6,7};

Volume(1) = {0,1,2,3,4,5,6,7};
//Transfinite Volume(1) = {0,1,2,3,4,5,6,7};
Physical Volume(1000) = {1};
//############# EXTRUSION TO 3D ###############