alphaInitializerFoam
-
发布一个VOF计算的前处理工具
alphaInitializerFoam
,用于精确初始化alpha
场。详见 https://github.com/daidezhi/alphaInitializerFoam -
http://dx.doi.org/10.17632/wg4sx7sc57.2 目前正在审核,大概24小时后可以访问。
-
我去,很屌的样子。。多面体网格怎么做的老铁?
-
@东岳 网格用STAR-CCM+做的,就是导出为OpenFOAM的2D网格有点麻烦,有时间出个教程。
-
久仰STAR-CCM+的多面体网格,感觉溜爆了
-
@东岳 他们的算法类似
polyDualMesh
,都是从四面体网格按照格点格式转换过来的,但是CCM会尽量减小转换后的多面体网格face
处的翘曲程度,对附面层网格的支持也很好。polyDualMesh
处理不好附面层网格,转换后的网格第二层和第一层网格厚度比会变大。 -
牛逼了老铁,很实用的功能。
-
@linhan-ge 希望能帮到有需要的人
-
@队长别开枪 一定早点出。
多问您一句,star-ccm+画网格好用吗?
-
@Hungryandfool 好用!教程得八月底了,最近在忙答辩。
-
@Hungryandfool ccm+的网格我觉得真的不错,尤其是对于复杂的结构,自动生成网格,就是不知道移植到OF中收敛性如何。
-
这个是我好像2015年做的
上面用的
polyDualMesh
-
@队长别开枪 OK。
随便问您一句,您那个工具,VOF的可以造波么? -
@mohui 有没有办法,跟您多交流下。
-
很棒的工具,想请问一下您,该工具能在有造波边界的的算例中使用吗?正在使用您的工具,但是报错了,可能是我没正确使用。
-
@guohuiqun 应该可以的。你的算例是啥?报错的信息是怎样的?有时候
Cork
布尔库会报错,不太稳定但是比CGAL
库快,下一版本会添加额外选项让用户指定使用Cork
还是CGAL
。 -
@Hungryandfool 不能,造波器是一类动边界,这个只用于初始化相场。但是如果要初始化类似正弦函数形状的界面还是可以的。
-
@队长别开枪 !!
不能识别我的造波边界 -
@guohuiqun 应该是编译的时候没有加载
waveModels
库,今天我开完组会加一下。 -
@队长别开枪 酷,真期待你的完整版与相关指南
-
@guohuiqun 搞好了,更新到Github页面了。
-
@队长别开枪 谢谢,我已经下载更新并使用了,想请教您该工具是用来处理六边形的网格吗,对于四边形的网格可以吗?我这边使用了您的工具出现下面情况
第一张图是我setFields得到的液面图
第二张是我使用alphaInitializerFoam waterShape.stl
不知道是不是我使用错误?麻烦您指点一下
-
@guohuiqun 工具是应用于任意多面体网格的。你是不是先运行
setFields
,然后直接运行alphaInitializerFoam
的?alphaInitializerFoam
不会清除alpha.water
已有的数据。用一个没有初始化过的alpha.water
文件试试。 -
@队长别开枪
我只使用alphaInitializerFoam waterShape.stl -
@guohuiqun 方便的话你把网格和stl文件发我邮箱(dezhi.dai@mavs.uta.edu)吧,不涉密的话。
-
@队长别开枪 好的,我已经发您邮箱了
-
@guohuiqun 通过邮箱回复你了。
-
@队长别开枪
立马尝试了一下,问题已解决(如上图所示),这前处理工具太棒了。再询问您一个比较愚蠢的问题,我使用alphaInitializerFoam waterShape.stl前处理之后,使用interFoam求解计算波浪问题,计算过程中液面捕捉利用的还是VOF,对吗?
-
@guohuiqun 是的。
-
@队长别开枪 好的,谢谢
-
@队长别开枪 期待您的教程,最近也用star-ccm+画网格给OF使用,想学习一下如何导出OpenFOAM的2D网格,自己尝试过,但是都失败了,如果时间方便的话能推荐一些学习教程就太感谢了。
-
@王金成 最近在准备答辩,有点忙,我会尽量抽时间做一下的。
-
@东岳
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中的验证算例。多面体网格还有很长的路要走。 -
你好,在v1812上编译该工具报错
请问如何解决 -
@夜阑烟寒 编译器?Linux发行版?看你的提示是没找到boost_system库,因为工具建立在第三方几何库上,建议使用Ubuntu或者Mint发行版,以及GNU C++。
-
@队长别开枪 红帽6.6系统 ,gcc版本是5.1.0
-
@队长别开枪 很奇怪,因为安装openfoam时候第三方包里是有个boost库的,应该在安装过程中编译进去了,应该是你的脚本默认的是系统的boost库,而非安装TheTirtyPart包的boost的库
-
@队长别开枪 毕竟我用的是v1812的环境,所以我的第三方编译出来
的是boost_1_64_0
我觉得是要把boost_system库改为boost_1_64_0才合理,但不知道要去改动哪部分 -
@guohuiqun 您好,看到您的这个示意图。我想问您,用setFields怎么设定给三角形结构啊?在region里不能用boxToCell,那么该用什么?我现在在试zoneToCell,但是过不去。希望您能给解惑一下。谢谢🙏
-
@zhe 三角形结构?我没有设置三角形的网格结构哈?
-
@guohuiqun 您图中从左向右,网格变化很像一个爬坡的结构,所以,我以为是类似三角形结构。不好意思,我搞错了。顺便问一句,您知道怎么给一个带坡度的结构使用setFields吗?谢谢
-
@zhe 直接设置水体区域,没有另外处理