[Gmsh] difference between gmsh svn nightly build

Benedikt Oswald benedikt.oswald at lspr.swiss
Mon Mar 13 12:11:04 CET 2017


Dear gmsh developers,

there seems to be quite a different behavior with respect to certain
.geo files between

the svn nightly builds of 27.2.2017 and 13.3.2017. I have attached the
.geo files.

Which version is more relevant ?

Thanks & greetings, Benedikt


 
----------------------------------------------------------------------------------------------------------------------
Dr. sc. techn. Benedikt Oswald - first engineer - LSPR AG - phone +41 43 366 90 74
Grubenstrasse 9, CH-8045 Zürich, benedikt.oswald at lspr.swiss
----------------------------------------------------------------------------------------------------------------------

-------------- next part --------------
// objective - gmsh opencascade script for creating a filleted cuboid embedded in free space
// creator - benedikt oswald, all rights reserved, 2017 march 3
// characteristic length scane - nanometer [nm]


SetFactory("OpenCASCADE");

Mesh.CharacteristicLengthMin = 10.0;
Mesh.CharacteristicLengthMax = 10.0;


SyncModel;

//////////////////////////////////////////////////////////////////////////////////
//specify dimensions
//////////////////////////////////////////////////////////////////////////////////
R_OUTER_SPHERE      =   300.0;

LENGTH_A            =   100.0;
LENGTH_B            =   100.0;
LENGTH_D            =   40.0;

CUBOOID_FILLET_LEN  =   0.1*(LENGTH_A+LENGTH_B+LENGTH_A)/3.0;   // choose 10% of the average edge length as the filleting radius


//////////////////////////////////////////////////////////////////////////////////
// specify physical tags
//////////////////////////////////////////////////////////////////////////////////
VACUUM_M                =   100;
CUBOID_M                =   400;
TRUNCATION_BOUNDARY_M   =   800;
FILLETED_CUBOID_M       =   500;





//////////////////////////////////////////////////////////////////////////////////
// specify internal volume id
//////////////////////////////////////////////////////////////////////////////////
OUTER_SPHERE_G      =   20;

CUBOID_G            =   50;
CUBOID_FILLETED_G   =   50;

SPHERICAL_SHELL_G   =   80;



//////////////////////////////////////////////////////////////////////////////////
// create the cuboid
//////////////////////////////////////////////////////////////////////////////////
xlo     =   -0.5*LENGTH_A;
xhi     =    0.5*LENGTH_A;

ylo     =   -0.5*LENGTH_B;
yhi     =    0.5*LENGTH_B;

zlo     =   -0.5*LENGTH_D;
zhi     =    0.5*LENGTH_D;

Block(CUBOID_G) = {xlo,ylo,zlo,xhi,yhi,zhi};


f() = Boundary{ Volume{CUBOID_G}; };
e() = Unique( Boundary{ Surface{f()}; } );

blob = Fillet{CUBOID_G}{e()}{CUBOOID_FILLET_LEN};

// delete the base of the filleted structure
Delete{ Volume{CUBOID_G}; Surface{f()}; Line{e()}; }


// create outer sphere
Sphere(OUTER_SPHERE_G) = {0,0,0,R_OUTER_SPHERE};


// in order to assign the physical surface tag TRUNCATION_BOUNDARY to surface SPHERICAL_SHELL, extract the surface from the outer sphere
b() = Boundary{ Volume{OUTER_SPHERE_G}; };
u() = Unique( Boundary{ Surface{b()}; } );
Physical Surface(TRUNCATION_BOUNDARY_M) = {u};




// make sure, the sphere & the cuboid are meshed simultaneously
BooleanDifference(SPHERICAL_SHELL_G) = { Volume{OUTER_SPHERE_G}; Delete; }{ Volume{blob}; };



// define the physical tags for the volumes
//Physical Volume(FILLETED_CUBOID_M) = {51};
Physical Volume(FILLETED_CUBOID_M) = {blob};
Physical Volume(VACUUM_M) = {SPHERICAL_SHELL_G};


SyncModel;


-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 882 bytes
Desc: OpenPGP digital signature
URL: <http://onelab.info/pipermail/gmsh/attachments/20170313/a7eb36b3/attachment.asc>


More information about the gmsh mailing list