[Gmsh] boolean with surface loop

[Christian Kluge]

Dear Christophe,

thank you very much for your help. We were able to calculate the difference between the cylinder and the fracture and to automate the process. The "BooleanDifference" operation applied to the cylinder and fracture, results in partially very small element sizes along the fracture side (see attached figure). This causes an incomplete cylinder (upper left) or too fine refinements (cylinder bottom). 

Is there a possibility to refine the mesh of the cylinder locally? 

Thank you very much.

All the best,
Christian



-----Ursprüngliche Nachricht-----
Von: Christophe Geuzaine [mailto:cgeuzaine at uliege.be] 
Gesendet: Freitag, 16. November 2018 22:58
An: Guido Blöcher <bloech at gfz-potsdam.de>
Cc: gmsh at geuz.org; Christian Kluge <ckluge at gfz-potsdam.de>
Betreff: Re: [Gmsh] boolean with surface loop



> On 16 Nov 2018, at 16:43, Guido Blöcher <bloech at gfz-potsdam.de> wrote:
> 
> Dear Christoph,
> 
> 
> 
> We try to mesh a three dimensional fracture embedded in a cylindrically shaped rock volume.
> 
> We can mesh the fracture and the cylinder separately but for later use as input for  a numerical simulation, we have to subtract or intersect the two volumes.
> 
> 
> 
> Unfortunately, gmsh provides an error message regarding our surface loop.
> 
> Gmsh complains that our surface loop is not a valid openCascade entity.
> 

Indeed: in your script you SetFactory("OpenCASCADE") after having created the first volume - which was thus created with the default "Built-In" CAD kernel.

Here's a working script:

-------------- next part --------------
A non-text attachment was scrubbed...
Name: 16x32.png
Type: image/png
Size: 1298188 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20181119/a2e25478/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 16x32.geo
Type: application/octet-stream
Size: 266830 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20181119/a2e25478/attachment-0001.geo>