CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    alphaInitializerFoam

    OpenFOAM
    9
    42
    8848
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 队长别开枪
      队长别开枪 教授 最后由 编辑

      发布一个VOF计算的前处理工具alphaInitializerFoam,用于精确初始化alpha场。详见 https://github.com/daidezhi/alphaInitializerFoam

      1 条回复 最后回复 回复 引用
      • 队长别开枪
        队长别开枪 教授 最后由 编辑

        http://dx.doi.org/10.17632/wg4sx7sc57.2 目前正在审核,大概24小时后可以访问。

        1 条回复 最后回复 回复 引用
        • 李东岳
          李东岳 管理员 最后由 编辑

          我去,很屌的样子。。多面体网格怎么做的老铁?

          替代文字

          线上CFD课程开始报名:http://www.dyfluid.com/class.html

          CFD高性能服务器 http://dyfluid.com/servers.html

          队长别开枪 1 条回复 最后回复 回复 引用
          • 队长别开枪
            队长别开枪 教授 @李东岳 最后由 编辑

            @东岳 网格用STAR-CCM+做的,就是导出为OpenFOAM的2D网格有点麻烦,有时间出个教程。

            H 1 条回复 最后回复 回复 引用
            • 李东岳
              李东岳 管理员 最后由 编辑

              久仰STAR-CCM+的多面体网格,感觉溜爆了

              线上CFD课程开始报名:http://www.dyfluid.com/class.html

              CFD高性能服务器 http://dyfluid.com/servers.html

              队长别开枪 1 条回复 最后回复 回复 引用
              • 队长别开枪
                队长别开枪 教授 @李东岳 最后由 编辑

                @东岳 他们的算法类似polyDualMesh,都是从四面体网格按照格点格式转换过来的,但是CCM会尽量减小转换后的多面体网格face处的翘曲程度,对附面层网格的支持也很好。polyDualMesh处理不好附面层网格,转换后的网格第二层和第一层网格厚度比会变大。

                1 条回复 最后回复 回复 引用
                • linhan.ge
                  linhan.ge 最后由 编辑

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

                  Working on CFDEM project.

                  队长别开枪 1 条回复 最后回复 回复 引用
                  • 队长别开枪
                    队长别开枪 教授 @linhan.ge 最后由 编辑

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

                    1 条回复 最后回复 回复 引用
                    • H
                      Hungryandfool @队长别开枪 最后由 编辑

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

                      多问您一句,star-ccm+画网格好用吗?

                      队长别开枪 M 2 条回复 最后回复 回复 引用
                      • 队长别开枪
                        队长别开枪 教授 @Hungryandfool 最后由 编辑

                        @Hungryandfool 好用!教程得八月底了,最近在忙答辩。:chouchou:

                        H 王 2 条回复 最后回复 回复 引用
                        • M
                          mohui @Hungryandfool 最后由 编辑

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

                          H 1 条回复 最后回复 回复 引用
                          • 李东岳
                            李东岳 管理员 最后由 编辑

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

                            替代文字

                            线上CFD课程开始报名:http://www.dyfluid.com/class.html

                            CFD高性能服务器 http://dyfluid.com/servers.html

                            队长别开枪 1 条回复 最后回复 回复 引用
                            • H
                              Hungryandfool @队长别开枪 最后由 编辑

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

                              队长别开枪 1 条回复 最后回复 回复 引用
                              • H
                                Hungryandfool @mohui 最后由 编辑

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

                                1 条回复 最后回复 回复 引用
                                • G
                                  guohuiqun 最后由 编辑

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

                                  队长别开枪 1 条回复 最后回复 回复 引用
                                  • 队长别开枪
                                    队长别开枪 教授 @guohuiqun 最后由 编辑

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

                                    G 1 条回复 最后回复 回复 引用
                                    • 队长别开枪
                                      队长别开枪 教授 @Hungryandfool 最后由 编辑

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

                                      1 条回复 最后回复 回复 引用
                                      • G
                                        guohuiqun @队长别开枪 最后由 编辑

                                        @队长别开枪 !!1.png
                                        2.png

                                        3.png
                                        不能识别我的造波边界

                                        队长别开枪 1 条回复 最后回复 回复 引用
                                        • 队长别开枪
                                          队长别开枪 教授 @guohuiqun 最后由 编辑

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

                                          G 1 条回复 最后回复 回复 引用
                                          • G
                                            guohuiqun @队长别开枪 最后由 编辑

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

                                            队长别开枪 1 条回复 最后回复 回复 引用
                                            • 队长别开枪
                                              队长别开枪 教授 @guohuiqun 最后由 编辑

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

                                              G 1 条回复 最后回复 回复 引用
                                              • G
                                                guohuiqun @队长别开枪 最后由 编辑

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

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

                                                队长别开枪 zhe 2 条回复 最后回复 回复 引用
                                                • 队长别开枪
                                                  队长别开枪 教授 @guohuiqun 最后由 编辑

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

                                                  G 1 条回复 最后回复 回复 引用
                                                  • G
                                                    guohuiqun @队长别开枪 最后由 编辑

                                                    @队长别开枪 微信图片_20190620190457.png
                                                    我只使用alphaInitializerFoam waterShape.stl

                                                    队长别开枪 1 条回复 最后回复 回复 引用
                                                    • 队长别开枪
                                                      队长别开枪 教授 @guohuiqun 最后由 编辑

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

                                                      G 1 条回复 最后回复 回复 引用
                                                      • G
                                                        guohuiqun @队长别开枪 最后由 编辑

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

                                                        队长别开枪 1 条回复 最后回复 回复 引用
                                                        • 队长别开枪
                                                          队长别开枪 教授 @guohuiqun 最后由 编辑

                                                          @guohuiqun 通过邮箱回复你了。

                                                          G 1 条回复 最后回复 回复 引用
                                                          • G
                                                            guohuiqun @队长别开枪 最后由 编辑

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

                                                            队长别开枪 1 条回复 最后回复 回复 引用
                                                            • 队长别开枪
                                                              队长别开枪 教授 @guohuiqun 最后由 编辑

                                                              @guohuiqun 是的。

                                                              G 1 条回复 最后回复 回复 引用
                                                              • G
                                                                guohuiqun @队长别开枪 最后由 编辑

                                                                @队长别开枪 好的,谢谢

                                                                1 条回复 最后回复 回复 引用
                                                                • 王
                                                                  王金成 @队长别开枪 最后由 编辑

                                                                  @队长别开枪 期待您的教程,最近也用star-ccm+画网格给OF使用,想学习一下如何导出OpenFOAM的2D网格,自己尝试过,但是都失败了,如果时间方便的话能推荐一些学习教程就太感谢了。

                                                                  队长别开枪 1 条回复 最后回复 回复 引用
                                                                  • 队长别开枪
                                                                    队长别开枪 教授 @王金成 最后由 编辑

                                                                    @王金成 最近在准备答辩,有点忙,我会尽量抽时间做一下的。

                                                                    1 条回复 最后回复 回复 引用
                                                                    • 队长别开枪
                                                                      队长别开枪 教授 @李东岳 最后由 编辑

                                                                      @东岳 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中的验证算例。多面体网格还有很长的路要走。

                                                                      1 条回复 最后回复 回复 引用
                                                                      • 夜
                                                                        夜阑烟寒 最后由 编辑

                                                                        你好,在v1812上编译该工具报错
                                                                        ec64997b-6b03-4130-9eb3-013393be9aaf-image.png
                                                                        请问如何解决

                                                                        队长别开枪 1 条回复 最后回复 回复 引用
                                                                        • 队长别开枪
                                                                          队长别开枪 教授 @夜阑烟寒 最后由 编辑

                                                                          @夜阑烟寒 编译器?Linux发行版?看你的提示是没找到boost_system库,因为工具建立在第三方几何库上,建议使用Ubuntu或者Mint发行版,以及GNU C++。

                                                                          夜 3 条回复 最后回复 回复 引用
                                                                          • 夜
                                                                            夜阑烟寒 @队长别开枪 最后由 编辑

                                                                            @队长别开枪 红帽6.6系统 ,gcc版本是5.1.0

                                                                            1 条回复 最后回复 回复 引用
                                                                            • 夜
                                                                              夜阑烟寒 @队长别开枪 最后由 编辑

                                                                              @队长别开枪 很奇怪,因为安装openfoam时候第三方包里是有个boost库的,应该在安装过程中编译进去了,应该是你的脚本默认的是系统的boost库,而非安装TheTirtyPart包的boost的库

                                                                              1 条回复 最后回复 回复 引用
                                                                              • 夜
                                                                                夜阑烟寒 @队长别开枪 最后由 编辑

                                                                                @队长别开枪 毕竟我用的是v1812的环境,所以我的第三方编译出来
                                                                                的是boost_1_64_0
                                                                                26bd7400-0a42-4818-8c0a-d22f15c6f42f-image.png
                                                                                54c087f4-3c7d-4192-b90a-c0ac1c7a0c56-image.png
                                                                                我觉得是要把boost_system库改为boost_1_64_0才合理,但不知道要去改动哪部分

                                                                                1 条回复 最后回复 回复 引用
                                                                                • zhe
                                                                                  zhe @guohuiqun 最后由 编辑

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

                                                                                  G 1 条回复 最后回复 回复 引用
                                                                                  • G
                                                                                    guohuiqun @zhe 最后由 编辑

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

                                                                                    zhe 1 条回复 最后回复 回复 引用
                                                                                    • zhe
                                                                                      zhe @guohuiqun 最后由 编辑

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

                                                                                      G 1 条回复 最后回复 回复 引用
                                                                                      • G
                                                                                        guohuiqun @zhe 最后由 编辑

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

                                                                                        1 条回复 最后回复 回复 引用
                                                                                        • First post
                                                                                          Last post