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]

牛逼了老铁，很实用的功能。

[队长别开枪]

@linhan-ge 希望能帮到有需要的人

[Hungryandfool]

@队长别开枪 一定早点出。

多问您一句，star-ccm+画网格好用吗？

[队长别开枪]

@Hungryandfool 好用！教程得八月底了，最近在忙答辩。

[mohui]

@Hungryandfool ccm+的网格我觉得真的不错，尤其是对于复杂的结构，自动生成网格，就是不知道移植到OF中收敛性如何。

[李东岳]

这个是我好像2015年做的 上面用的polyDualMesh

[Hungryandfool]

@队长别开枪 OK。
随便问您一句，您那个工具，VOF的可以造波么？

[Hungryandfool]

@mohui 有没有办法，跟您多交流下。

[guohuiqun]

很棒的工具，想请问一下您，该工具能在有造波边界的的算例中使用吗？正在使用您的工具，但是报错了，可能是我没正确使用。

[队长别开枪]

@guohuiqun 应该可以的。你的算例是啥？报错的信息是怎样的？有时候Cork布尔库会报错，不太稳定但是比CGAL库快，下一版本会添加额外选项让用户指定使用Cork还是CGAL。

[队长别开枪]

@Hungryandfool 不能，造波器是一类动边界，这个只用于初始化相场。但是如果要初始化类似正弦函数形状的界面还是可以的。

[guohuiqun]

@队长别开枪 !!

不能识别我的造波边界

[队长别开枪]

@guohuiqun 应该是编译的时候没有加载waveModels库，今天我开完组会加一下。

[guohuiqun]

@队长别开枪 酷，真期待你的完整版与相关指南

[队长别开枪]

@guohuiqun 搞好了，更新到Github页面了。

[guohuiqun]

@队长别开枪 谢谢，我已经下载更新并使用了，想请教您该工具是用来处理六边形的网格吗，对于四边形的网格可以吗？我这边使用了您的工具出现下面情况
第一张图是我setFields得到的液面图

第二张是我使用alphaInitializerFoam waterShape.stl

不知道是不是我使用错误？麻烦您指点一下

[队长别开枪]

@guohuiqun 工具是应用于任意多面体网格的。你是不是先运行setFields，然后直接运行alphaInitializerFoam 的？alphaInitializerFoam 不会清除alpha.water已有的数据。用一个没有初始化过的alpha.water文件试试。

[guohuiqun]

@队长别开枪
我只使用alphaInitializerFoam waterShape.stl

[队长别开枪]

@guohuiqun 方便的话你把网格和stl文件发我邮箱(dezhi.dai@mavs.uta.edu)吧，不涉密的话。

[guohuiqun]

@队长别开枪 好的，我已经发您邮箱了

[队长别开枪]

@guohuiqun 通过邮箱回复你了。

[guohuiqun]

@队长别开枪
立马尝试了一下，问题已解决（如上图所示），这前处理工具太棒了 。再询问您一个比较愚蠢的问题，我使用alphaInitializerFoam waterShape.stl前处理之后，使用interFoam求解计算波浪问题，计算过程中液面捕捉利用的还是VOF，对吗？

[队长别开枪]

@guohuiqun 是的。

[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才合理，但不知道要去改动哪部分

[zhe]

@guohuiqun 您好，看到您的这个示意图。我想问您，用setFields怎么设定给三角形结构啊？在region里不能用boxToCell，那么该用什么？我现在在试zoneToCell,但是过不去。希望您能给解惑一下。谢谢🙏

[guohuiqun]

@zhe 三角形结构？我没有设置三角形的网格结构哈？

[zhe]

@guohuiqun 您图中从左向右，网格变化很像一个爬坡的结构，所以，我以为是类似三角形结构。不好意思，我搞错了。顺便问一句，您知道怎么给一个带坡度的结构使用setFields吗？谢谢

[guohuiqun]

@zhe 直接设置水体区域，没有另外处理