alphaInitializerFoam
-
@东岳
polyDualMesh
转换的多面体网格在锐边会有碎面,原先我们按照 https://openfoamwiki.net/index.php/Polyhedral_mesh_generation 里的方法修复二维的多边形网格的,在对三维的进行修复的时候发现很难保持锐边的形状,这个工具不知道最新版本的如何,原先对边界层网格也不友好,之后我们就转向STAR CCM+
了。现在主流多面体网格生成软件,商用的,开源的,基本都是先生成传统四面体网格,然后按照格心格式转换过来的,带来的问题是生成的多面体网格的面单元有时翘曲比较严重(STAR CCM+
对这个有优化,但也不能百分百保证没有翘曲),对PLIC
一类的几何VOF
方法不友好,所以一般会对翘曲面进行三角化,但是这个三角化面中心选的不好又会形成凹(体)单元 (所以最近Dr. Lopez等人开发出了任意凹单元界面重构快速算法 https://www.sciencedirect.com/science/article/pii/S0021999119303080 ,我写好了接口,最近正在测试,可以让interPlicFoam
调用他们组的VOFTools库(Fortran)),所以有时候蛮让人头大的。第二类多面体网格生成方法就比较直接了,先生成规则六面体网格(不是结构网格,结构非结构不能靠单元形状判断,根本区别在数据结构),然后每个六面体单元分解为24个四面体单元,然后再按照格心格式转成多面体网格,这样面单元就没有翘曲了,但是第二种方法很难应用在包含复杂几何形状的实际工程算例中,多用于paper中的验证算例。多面体网格还有很长的路要走。