Re: [vmtk-users] cylinder mesh with boundary layer

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

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...@ho... 
>>> <mailto:jan...@ho...>> 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
>>>     vmt...@li...
>>>     <mailto:vmt...@li...>
>>>     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
>> vmt...@li...
>> https://lists.sourceforge.net/lists/listinfo/vmtk-users
>