Hi Luca,
There is indeed a boundary layer that is hiding from me:)
Thanks a lot!
Jan
On 02/27/2013 02:01 PM, Luca Antiga wrote:
Hi Jan,
are you sure that no boundary layer is generated? Note that in the
current release
version you can't see the boundary layer from outside, you have to
inspect the inner
mesh to see it, since a boundary layer is generated at the endcaps as
well.
Try to cut the mesh using, for instance, Paraview and see if you can
see the prismatic layer.
Thanks to David Steinman's support I'm currently working on a boundary
layer
code that avoids the boundary layers at the endcaps and that allows
for much thicker
boundary layers. The code is currently at https://github.com/lantiga/vmtk,
"bl" branch. It will be merged to master as soon as it's stable and
tested.
Hope this helps
Luca
On Feb 27, 2013, at 12:59 PM, Jan Helmig wrote:
Hi Arjan,
thanks, but I have already tried vmtksurfaceremeshing with different
options but it does not work as well.
Jan
On 02/27/2013 11:41 AM, Arjan Geers wrote:
Hi Jan,
Try vmtksurfaceremeshing instead of vmtksurfacetriangle before
generating the volumetric mesh. The large aspect ratio of the
triangles is likely what throws off vmtkmeshgenerator.
Arjan
On Wed, Feb 27, 2013 at 9:31 AM, Jan Helmig <jan_hel...@hotmail.com
<mailto:jan_hel...@hotmail.com>> wrote:
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
------------------------------------------------------------------------------
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
vmtk-users@lists.sourceforge.net
<mailto:vmtk-users@lists.sourceforge.net>
https://lists.sourceforge.net/lists/listinfo/vmtk-users
------------------------------------------------------------------------------
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
vmtk-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/vmtk-users
------------------------------------------------------------------------------
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
vmtk-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/vmtk-users
