[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