boundary layer
i'm tryng to make a 2D hibryd mesh in .occ model and i get it in the following way..
import gmsh
import math
import os
import sys
from pathlib import Path
gmsh.initialize()
path = os.path.dirname(os.path.abspath(__file__))
v = gmsh.model.occ.importShapes(os.path.join(path, '2D_simple_engine_big_domain.step'))
gmsh.model.occ.synchronize()
N=53 #number of layer
r =1.2 # ratio
d = [0.00001] # thickness of first layer
for i in range(1, N): d.append((d[0]) * r**i+d[-1])
d.append(1)
gmsh.model.occ.extrude([(1,19),(1,3)],0,200,0,[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1],d,True)
gmsh.model.occ.cut([(2,1)],[(2,2),(2,3)],-1,True,False)
gmsh.model.occ.synchronize()
gmsh.model.mesh.setTransfiniteCurve(3,200)
gmsh.model.mesh.setTransfiniteCurve(19,200)
gmsh.model.mesh.setTransfiniteCurve(27,30)
gmsh.model.mesh.setTransfiniteCurve(44,50,"Beta",coef=1.1)
gmsh.model.mesh.setTransfiniteCurve(28,50,"Beta",coef=-1.1)
gmsh.model.occ.synchronize()
if '-nopopup' not in sys.argv:
gmsh.fltk.run()
gmsh.finalize()
Is it possible to trasform the command gmsh.model.geo.extrudeBoundaryLayer into a .occ command and so trasform the line separating the two typs of mesh into a discrete line (as in the naca_boundary_layer example ). Otherwise how can i call up a specific entity from .occ to .geo in order to use that command?
Edited by Guido Vallifuoco