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> 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
>> 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
