[Gmsh] MathEval field in a circle with Open Cascade

Hector Gabriel Espinoza Román hespinoza at utb.edu.co
Fri Jun 22 00:55:15 CEST 2018


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?


Sincerely,

Héctor
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20180621/150d435c/attachment.html>


More information about the gmsh mailing list