[Gmsh] MathEval field in a circle with Open Cascade

[Christophe Geuzaine]

> On 22 Jun 2018, at 00:55, Hector Gabriel Espinoza Román <hespinoza at utb.edu.co> wrote:
> 
> Hello:
> 
> This simple script should draw a circle (2D surface) and when meshed it should be refined on the boundary, but I am not getting that. I have done the same without Open Cascade and it works.
> 
> cx = 0.0 ;
> cy = 0.0 ;
> cz = 0.0 ;
> D = 1.0 ;
> R = D/2;
> SetFactory("OpenCASCADE");
> Disk(1) = {cx,cy,cz,R} ; // creates a circle surface
> Field[1] = MathEval ;
> //Field[1].F = "0.15 - 0.13*((x^2)/(0.5^2))" ; // works 2D, 617 elements: 0.02 to 0.15
> Field[1].F = "0.15 - 0.13*(x^2 + y^2)/(0.5^2)" ; // doesnt work 2D
> Background Field = 1 ;
> 
> What I am doing wrong?
> 

Add

Mesh.CharacteristicLengthExtendFromBoundary = 0;

Christophe


> Sincerely,
> 
> Héctor
> _______________________________________________
> gmsh mailing list
> gmsh at onelab.info <mailto:gmsh at onelab.info>
> 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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20180622/eedafc24/attachment.html>