> On 1 Aug 2019, at 20:07, Kópházi József <[email protected]> wrote: > > 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). >
It's a bug in the MeshAdapt algorithm. Using the frontal algo (Mesh.Algorithm = 6) gives the expected result. Can you open an issue in our gitlab so we can track this problem? Thanks, Christophe > <bghpnghkeenmmnod.png> > > How could I get a proper, isotropic mesh? > > Thank you for your help in advance, > > Jozsef > > > > <p.png>_______________________________________________ > gmsh mailing list > [email protected] > http://onelab.info/mailman/listinfo/gmsh — Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine _______________________________________________ gmsh mailing list [email protected] http://onelab.info/mailman/listinfo/gmsh
