PythonAPI: Fields: BoundaryLayer
I'm having an issue meshing using the BoundaryLayer Field method. It doesn't appear to mesh (below is a simplified example that demonstrates the issue:
import gmsh
gmsh.initialize()
gmsh.clear()
# definitions
gmsh.model.add(nameModel)
geoMod = gmsh.model.geo
meshSize = 5
# 2D geometry
P6 = geoMod.addPoint(0, 0, 0, meshSize)
P5 = geoMod.addPoint(2000, 0, 0, meshSize)
P4 = geoMod.addPoint(2000, 600, 0, meshSize)
P3 = geoMod.addPoint(0, 600, 0, meshSize)
P2 = geoMod.addPoint(0, 310, 0, meshSize)
P1 = geoMod.addPoint(0, 290, 0, meshSize)
LMouth = geoMod.addLine(P1, P2)
L1 = geoMod.addLine(P2, P3)
L2 = geoMod.addLine(P3, P4)
L3 = geoMod.addLine(P4, P5)
L4 = geoMod.addLine(P5, P6)
L5 = geoMod.addLine(P6, P1)
Curve = geoMod.addCurveLoop([LMouth, L1, L2, L3, L4, L5])
Surf = geoMod.addPlaneSurface([Curve])
geoMod.synchronize()
# hex-mesh
gmsh.model.mesh.field.add("BoundaryLayer", 1)
gmsh.model.mesh.field.setNumber(1, "Quads", 1)
gmsh.model.mesh.field.setNumber(1, "Size", 2)
gmsh.model.mesh.field.setNumber(1, "NbLayers", 20)
gmsh.model.mesh.field.setNumber(1, "Thickness", 200)
gmsh.model.mesh.field.setNumber(1, "Ratio", 1.2)
gmsh.model.mesh.field.setNumber(1, "CurvesList", L3)
gmsh.model.mesh.field.setNumber(1, "PointsList", P4)
gmsh.model.mesh.field.setNumber(1, "PointsList", P5)
gmsh.model.mesh.field.setNumber(1, "SizeFar", 5)
gmsh.model.mesh.field.setAsBoundaryLayer(1)
gmsh.option.setNumber("Mesh.Algorithm", 5)
geoMod.synchronize()
No errors are seen in Python nor GMSH's Message Console's output.
Any advice, or am I missing something? Thanks in advance, Nick
Edited by Nick Howlett