[Gmsh] structured mesh in a wedge

[Deepak Garg]

Here is your file. Try with this.


Deepak

On 16 May 2016 at 23:37, Deepak Garg <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> 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]
>> *Sent:* Monday, May 16, 2016 12:09 PM
>> *To:* CHARCHI AGHDAM, ALI <alic at email.sc.edu>
>> *Cc:* 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
>>
>> -------------------------
>>
>>
>>
>> On Thu, May 12, 2016 at 4:50 PM, CHARCHI AGHDAM, ALI <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
>> http://onelab.info/mailman/listinfo/gmsh
>>
>>
>>
>> _______________________________________________
>> gmsh mailing list
>> gmsh at onelab.info
>> http://onelab.info/mailman/listinfo/gmsh
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20160516/41e49c12/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wedge.geo
Type: application/octet-stream
Size: 1038 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20160516/41e49c12/attachment-0001.geo>