Implementing simple box fields
Hi Christophe,
I took a look at the tutorials but can't seem to find this. Maybe you can point me to the issue.
I wanted to mesh a simple box, with a small portion to have a refined mesh. I use a simple Box
field, namely
import gmsh
import sys
gmsh.initialize(sys.argv)
min_mesh_size = 0.01
max_mesh_size = 0.1
gmsh.model.add("boxwithRefinements")
surf1 = 1
L = 1
H = 1
eps=1.e-3
c = 0.5*L
thick = 0.05 # dim/20 just for check
box = gmsh.model.occ.add_box(0, 0, 0, L, H, thick, tag=1)
eps_box_field = 0.1*H
gmsh.model.occ.synchronize()
gmsh.model.mesh.field.add("Box", 1)
gmsh.model.mesh.field.setNumber(1, "VIn", min_mesh_size/100.)
gmsh.model.mesh.field.setNumber(1, "VOut", min_mesh_size)
gmsh.model.mesh.field.setNumber(1, "XMin", c-eps)
gmsh.model.mesh.field.setNumber(1, "XMax", L+eps)
gmsh.model.mesh.field.setNumber(1, "YMin", 0.4*H)
gmsh.model.mesh.field.setNumber(1, "YMax", 0.6*H)
gmsh.model.mesh.field.setNumber(1, "ZMin", -eps)
gmsh.model.mesh.field.setNumber(1, "ZMax", thick + eps)
gmsh.model.mesh.field.setNumber(1, "Thickness", 0.5)
# gmsh.option.setNumber("Mesh.MeshSizeExtendFromBoundary", 0)
# gmsh.option.setNumber("Mesh.MeshSizeFromPoints", 0)
# gmsh.option.setNumber("Mesh.MeshSizeFromCurvature", 0)
gmsh.model.mesh.generate(3)
gmsh.write("boxwithRefinements.msh")
if '-nopopup' not in sys.argv:
gmsh.fltk.run()
gmsh.finalize()
but this doesn't result as expected. I expected a finer mesh in the box [XMin, XMax] x [YMin, YMax] x [ZMin, ZMax]