[vmtk-users] cylinder mesh with boundary layer

[Jan H. <jan...@ho...>]

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