Hi Laura,

I cannot reproduce the crash. Can you try with the latest snaphot (not the 
stable release, which is quite old)?

Two remarks:

- your mesh size is too coarse and leads to a self-intersecting mesh
- you probably want to use BooleanFragments{ Volume{1,2,3}; Delete; }{} instead 
of the 2 difference operations to get a conformal mesh at the internal 
interfaces?

Christophe


> On 19 Jul 2018, at 18:00, Laura Holewa <[email protected]> wrote:
> 
> Hello all,
> 
> I have three concentric cylinders. The first physical volume that I want to 
> define is the volume between the outer most cylinder and the middle cylinder. 
> The second physical volume that I want to define is the volume between the 
> middle cylinder and the inner most cylinder. When I run this script gmsh 
> tells me that there are subfaces intersecting and crashes. How can I fix this 
> script? Could this be a possible bug?
> 
> ------------
> 
> // Gmsh project created on Wed Jul 18 13:52:12 2018
> SetFactory("OpenCASCADE");
> 
> Cylinder(1) = {0, 0, 0, 0, 0, 100, 10, 2*Pi};
> Cylinder(2) = {0, 0, 10, 0, 0, 80, 5, 2*Pi};
> Cylinder(3) = {0, 0, 15, 0, 0, 60, 2, 2*Pi};
> 
> BooleanDifference(6) = { Volume{1};Delete; }{ Volume{2}; };
> BooleanDifference(7) = { Volume{2};Delete; }{ Volume{3}; };
> 
> Physical Volume ("outer") = {6};
> Physical Volume ("inner") = {7};
> 
> -------------
> 
> Thanks,
> 
> Laura
> _______________________________________________
> gmsh mailing list
> [email protected] <mailto:[email protected]>
> http://onelab.info/mailman/listinfo/gmsh 
> <http://onelab.info/mailman/listinfo/gmsh>
— 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science 
http://www.montefiore.ulg.ac.be/~geuzaine

Free software: http://gmsh.info | http://getdp.info | http://onelab.info

_______________________________________________
gmsh mailing list
[email protected]
http://onelab.info/mailman/listinfo/gmsh

Reply via email to