[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