[Gmsh] How to set element size on the surface of a cylinder in OpenCascade?

Kópházi József j.kophazi at imperial.ac.uk
Thu Aug 1 20:07:13 CEST 2019


Dear Gmsh community,

I'm fighting with Gmsh to set up a model of a quarter cylinder using the 
new OpenCascade based facilities. The following script is considered:

SetFactory("OpenCASCADE");

cylinder = newv;
Cylinder(cylinder)         = {  0.0  ,  0.0  ,  0.0  , 0.0  ,  0.0  ,  
60.0,  175.0  };
box = newv;
Box(box)                 = {  0.0  ,  0.0  , -50.0 ,200.0  , 200.0  
,200.0  };
lst() = BooleanIntersection { Volume{ cylinder }; Delete; }
                             { Volume{ box }; Delete; };
pt_volume = lst[0];

Characteristic Length{ PointsOf{ Volume{:}; } } = 6.0;

However, the surface mesh on the (remaining part of the) cylinder is 
rather odd, it contains one single layer of elements axially. (Please 
find attached the image).

How could I get a proper, isotropic mesh?

Thank you for your help in advance,

Jozsef


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20190801/75222449/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: bghpnghkeenmmnod.png
Type: image/png
Size: 25651 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20190801/75222449/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: p.png
Type: image/png
Size: 25651 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20190801/75222449/attachment-0003.png>


More information about the gmsh mailing list