Hi,I am trying to create a cylinder mesh with boundary layer cells. The created mesh looks fine except the fact that there is no boundary layer created. As a starting point I use an open, non-capped cylinder surface created with DOLFIN (attached to the mail). Then I use the following commands.
vmtk = []
outfile = "cylinder.xml"
outfiles.append(outfile)
vmtk.append("vmtksurfacetransform -ifile cylinder.vtp -rotation 0.0
0.0 90.0 -translation 0.0 %f 0.0" %float(height/2))
#vmtk.append("vmtksurfaceviewer")
vmtk.append("vmtksurfacetriangle")
# Mesh generator (comment out unwanted lines)
cmd = "vmtkmeshgenerator"
cmd += " -edgelength %f" %0.11
cmd += " -maxedgelength %f" %(1.2*0.12)
cmd += " -minedgelength %f" %(0.8*0.12)
cmd += " -boundarylayer 1"
cmd += " -thicknessfactor 0.1"
cmd += " -tetrahedralize 1"
vmtk.append(cmd)
vmtk.append("vmtkmeshwriter -entityidsarray CellEntityIds")
vmtk = " --pipe ".join(vmtk)
vmtk = vmtk+" -ofile %s" %outfile
print vmtk
subprocess.call(vmtk, shell=True)
Does anybody understand why there is no boundary layer created?
Maybe it has something to do that the surface consists out of very
stretched triangles after the command "vmtksurfacetriangle".
Thanks a lot in advance, Jan
cylinder.vtp
Description: Binary data
------------------------------------------------------------------------------ Everyone hates slow websites. So do we. Make your web apps faster with AppDynamics Download AppDynamics Lite for free today: http://p.sf.net/sfu/appdyn_d2d_feb
_______________________________________________ vmtk-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/vmtk-users
