[Gmsh] structured mesh in a wedge

CHARCHI AGHDAM, ALI alic at email.sc.edu
Tue May 17 17:40:02 CEST 2016


Hi Deepak,
Still not working. I guess the issue is not from geometry definition section but from mesh options.

Ali

From: Deepak Garg [mailto:deepak.garg at ingv.it]
Sent: Monday, May 16, 2016 5:43 PM
To: CHARCHI AGHDAM, ALI <alic at email.sc.edu>
Cc: Felix Salazar <felix.a.salazar at gmail.com>; gmsh at onelab.info
Subject: Re: [Gmsh] structured mesh in a wedge

Here is your file. Try with this.

Deepak

On 16 May 2016 at 23:37, Deepak Garg <deepak.garg at ingv.it<mailto:deepak.garg at ingv.it>> wrote:
Hi Ali
In place of     Plane Surface(6) = {5};
you should write   Ruled Surface(6) = {5};
Best
Deepak

On 16 May 2016 at 20:17, CHARCHI AGHDAM, ALI <alic at email.sc.edu<mailto:alic at email.sc.edu>> wrote:
Hi Felix,
Thanks for your reply. I made the changes you suggested but I am still getting some errors:

Info    : Meshing surface 14 (transfinite)
.Error   : Surface 14 cannot be meshed using the transfinite algo 1 != 19
Info    : Meshing surface 14 (extruded)
Info    : Meshing surface 18 (transfinite)
Info    : Meshing surface 21 (transfinite)
.Error   : Surface 21 cannot be meshed using the transfinite algo 19 != 1

My new geo file is:

n1=10;
n2=20;

Point(1) =  {0.000000000000e+00 , 0.000000000000e+00 , 0.000000000000e+00 , 1e-5 };
Point(2) =  {1.998096443164e-05 , 0.000000000000e+00 , -8.723877473067e-07 , 1e-5 };
Point(3) =  {1.198857865898e-04 , 0.000000000000e+00 , -5.234326483840e-06 , 1e-5 };
Point(4) =  {0.000000000000e+00 , 1.400000000000e-04 , 0.000000000000e+00 , 1e-5 };
Point(5) =  {0.000000000000e+00 , 8.000000000000e-05 , 0.000000000000e+00 , 1e-5 };
Rotate {{0, 1, 0}, {0, 0, 0}, -0.04363323129985824} {
  Point{5, 4, 2, 3};
}

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

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

Extrude {{0, 1, 0},{0,0,0},0.08726646259971647}{
   Surface{6};
   Layers{1};
   Recombine;
}

Surface Loop(24) = {6, 14, 18, 21, 23};
Volume(25) = {24};

Transfinite Line {1,2} = n1;
Transfinite Line {3,4} = n2;

Transfinite Surface {6} = {};
Transfinite Surface {21} = {3,4,13};
Transfinite Surface {14} = {2,9,5};
Transfinite Surface {18} = {};
Transfinite Surface {23} = {};

Recombine Surface {6,21,14,18,23};

Transfinite Volume{25} = {2,3,4,5,9,13,4,5};

Physical Surface("front") = {6};
Physical Surface("back") = {23};
Physical Surface("right") = {21};
Physical Surface("down") = {18};
Physical Surface("electrode") = {14};

Note that 14 and 21 surfaces are bounded by 3 lines so, I have specified the corners.

Ali

From: Felix Salazar [mailto:felix.a.salazar at gmail.com<mailto:felix.a.salazar at gmail.com>]
Sent: Monday, May 16, 2016 12:09 PM
To: CHARCHI AGHDAM, ALI <alic at email.sc.edu<mailto:alic at email.sc.edu>>
Cc: gmsh at onelab.info<mailto:gmsh at onelab.info>
Subject: Re: [Gmsh] structured mesh in a wedge

Hello,

All your surfaces should be Transfinite. And if they are bounded by 4 lines, you don't need to specify the corners. You might try:

Transfinite Surface {6} = {};
Transfinite Surface {21} = {};
Transfinite Surface {14} = {};
Transfinite Surface {18} = {};
Transfinite Surface {23} = {};

Recombine Surface {6,21,14,18,23};

Although, you might just need to increase the number of layers in your extrusion. This PC does not have GMSH installed so I can't check your GEO file

-------------------------
Félix Salazar
felix.a.salazar at gmail.com<mailto:felix.a.salazar at gmail.com>
-------------------------

On Thu, May 12, 2016 at 4:50 PM, CHARCHI AGHDAM, ALI <alic at email.sc.edu<mailto:alic at email.sc.edu>> wrote:
Hi everyone,
I would like to mesh a wedge using gmsh. I have already been able to mesh all the surfaces but the 3D mesh module fails to extrude the 2D mesh and generate 3D mesh. I am getting the following error:

Incompatible surface 14 in transfinite volume 25.

And here is my geo file:

n1=20;
n2=40;

Point(1) =  {0.000000000000e+00 , 0.000000000000e+00 , 0.000000000000e+00 , 1e-5 };
Point(2) =  {1.998096443164e-05 , 0.000000000000e+00 , -8.723877473067e-07 , 1e-5 };
Point(3) =  {1.198857865898e-04 , 0.000000000000e+00 , -5.234326483840e-06 , 1e-5 };
Point(4) =  {0.000000000000e+00 , 1.400000000000e-04 , 0.000000000000e+00 , 1e-5 };
Point(5) =  {0.000000000000e+00 , 8.000000000000e-05 , 0.000000000000e+00 , 1e-5 };
Rotate {{0, 1, 0}, {0, 0, 0}, -0.04363323129985824} {
  Point{5, 4, 2, 3};
}

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

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

Extrude {{0, 1, 0},{0,0,0},0.08726646259971647}{
   Surface{6};
   Layers{1};
   Recombine;
}

Surface Loop(24) = {6, 14, 18, 21, 23};
Volume(25) = {24};

Transfinite Line {1,2} = n1 Using Progression 1;
Transfinite Line {3,4} = n2 Using Progression 1;

Transfinite Surface {6} = {2,3,4,5};

Recombine Surface {6,21,14,18,23};

Transfinite Volume{25} = {2,3,4,5,9,13,4,5};

Would you please help?

Best regards,

Ali Charchi

PhD Student of Mechanical Engineering
University of South Carolina
Columbia, SC



_______________________________________________
gmsh mailing list
gmsh at onelab.info<mailto:gmsh at onelab.info>
http://onelab.info/mailman/listinfo/gmsh


_______________________________________________
gmsh mailing list
gmsh at onelab.info<mailto:gmsh at onelab.info>
http://onelab.info/mailman/listinfo/gmsh


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20160517/7bf3bbb0/attachment-0001.html>


More information about the gmsh mailing list