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