CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    OF有欧拉拉格朗日求解器吗

    OpenFOAM
    4
    27
    3572
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • S
      Samuel-Tu 最后由 编辑

      铁铁们,看见东岳老师在白话多相流里面说欧拉拉格朗日法已经被广泛植入OF中,现在OF有这个求解器吗?

      1 条回复 最后回复 回复 引用
      • S
        Samuel-Tu 最后由 编辑

        好像搜索按钮不见了。。。

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

          才发现 我看看

          CFD高性能服务器 http://dyfluid.com/servers.html
          2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

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

            欧拉拉格朗日就是DPMFoam啊

            论坛最近迁移服务器,才发现搜索功能丢了,我去找找,铁铁

            CFD高性能服务器 http://dyfluid.com/servers.html
            2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

            S 1 条回复 最后回复 回复 引用
            • S
              Samuel-Tu @李东岳 最后由 编辑

              @东岳 好像还是用不了,我搜欧拉,拉格朗日等词都都不出来。

              1 条回复 最后回复 回复 引用
              • S
                Samuel-Tu 最后由 编辑

                最近有这么一个问题,就是给定了一个稳态的流场了,然后呢我想根据这个稳态的流场驱动颗粒的运动。这个不知道DPM能做到吗
                另外我想考虑的是这个流场中的颗粒的粒径是不同的。

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

                  icoUncoupleKinematicParcelFoam

                  CFD高性能服务器 http://dyfluid.com/servers.html
                  2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

                  1 条回复 最后回复 回复 引用
                  • 星
                    星星星星晴 @Samuel-Tu 最后由 编辑

                    @Samuel-Tu
                    可以的,具体教程见下
                    https://www.foamacademy.com/wp-content/uploads/2018/03/particles_slides.pdf

                    m.sui20@foxmail.com

                    S 1 条回复 最后回复 回复 引用
                    • S
                      Samuel-Tu @星星星星晴 最后由 编辑

                      @星星星星晴 非常感谢,太有用了,还有编译过程。。这个是哪个大学的教材吗。有没有官网呀,看看还有没有其他的教程

                      星 1 条回复 最后回复 回复 引用
                      • 星
                        星星星星晴 @Samuel-Tu 最后由 编辑

                        @Samuel-Tu google 大学
                        哈哈哈 搜的

                        m.sui20@foxmail.com

                        1 条回复 最后回复 回复 引用
                        • Q
                          qingdong_wang 最后由 编辑

                          @Samuel-Tu 这是德国罗斯托克大学的, 我也找这个做的, 不过我做的是气泡流, 加上粒子的影响。 你的颗粒是一开始静止的吗? 我写了一个python脚本, 来注射粒子。你需要的话, 我可以贴上来,但是我写的是均匀直径的, 只要用random做一下就可以变直径了。

                          李东岳 S 2 条回复 最后回复 回复 引用
                          • 李东岳
                            李东岳 管理员 @qingdong_wang 最后由 编辑

                            @qingdong_wang 气泡流欧拉拉格朗日模拟?结果怎么样?我现在也在做,发现dispersion模型需要改进。

                            CFD高性能服务器 http://dyfluid.com/servers.html
                            2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

                            星 Q 2 条回复 最后回复 回复 引用
                            • 星
                              星星星星晴 @李东岳 最后由 编辑

                              @东岳

                              我们这边同事用的这个dispersion模型

                              J. Lipowsky, M. Sommerfeld, “LES-Simulation of the formation of particle
                              strands in swirling flows using an unsteady Euler-Lagrange approach,” in
                              6th International Conference on Multiphase Flow ICMF 2007, Leipzig, Germany,
                              2007.

                              m.sui20@foxmail.com

                              1 条回复 最后回复 回复 引用
                              • S
                                Samuel-Tu @qingdong_wang 最后由 编辑

                                @qingdong_wang 太感谢了,麻烦贴出来学习下。我确实是想做变直径的。我是做稳态风场中雪颗粒的运动。。你是按这个教程里面开发的吗,不知道OF7能编译通过吗。。

                                Q 1 条回复 最后回复 回复 引用
                                • Q
                                  qingdong_wang @李东岳 最后由 编辑

                                  @东岳 这个我不是很清楚啊, 我用laminar。关于后处理想请教您一个问题, 如何用preview将粒子显示在alpha.water 的上面。 我现在只会用extract block来显示粒子以及粒子的glyph, 但是当我把alpha.water调成visible之后,他就覆盖粒子的分布和粒子的速度矢量。 我搜了搜网上没找到怎么弄:134:

                                  1 条回复 最后回复 回复 引用
                                  • Q
                                    qingdong_wang @Samuel-Tu 最后由 编辑

                                    @Samuel-Tu 可以通的过,我本想把我照着他做的不可压缩单相流贴在这边, 但是文件大小有限制。 那个python我也顺便贴一下。这个是配合manuelInjection model一起使用。你到时候使用的时候改一下用random函数改一下直径就好了。如果有问题, 也烦请你再联系我一下, 也让我改一下自己的错误。manualInjection.zip 在这个压缩包里,因为Python文件不能直接上传

                                    S 2 条回复 最后回复 回复 引用
                                    • S
                                      Samuel-Tu @qingdong_wang 最后由 编辑

                                      @qingdong_wang 感谢了,有任何进展我都会及时和你联系的:baobao:

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

                                        稳态风场中雪颗粒的运动

                                        @Samuel-Tu 为什么不直接用icoUncoupleKinematicParcelFoam,你不需要做任何编译

                                        @qingdong_wang 你可以打开两个case,一个显示粒子,一个现实alpha :142:

                                        CFD高性能服务器 http://dyfluid.com/servers.html
                                        2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

                                        S 2 条回复 最后回复 回复 引用
                                        • S
                                          Samuel-Tu @李东岳 最后由 编辑

                                          @东岳 嗯嗯,我正在看icoUncoupleKinematicParcelFoam的源代码。另外我看了下射入颗粒的直径分布,有那种随机直径分布函数,normalDistribution。但是我想做的可能是PantchInjection,然后呢,入射的颗粒是跟patch的z坐标有关的。就是说这个patch,低处颗粒直径大,高处颗粒直径小。这个好像没有这个功能。
                                          另外一个功能是想,颗粒和某一个壁面碰撞后,反射角度和入射角度是不一样的,并且根据和壁面碰撞的能量,这个壁面还会激发颗粒从这个壁面入射。不知道OF有没有这个功能。。
                                          感觉可能要研究下Inject的代码了,可能要自己修改

                                          1 条回复 最后回复 回复 引用
                                          • S
                                            Samuel-Tu @qingdong_wang 最后由 编辑

                                            @qingdong_wang 非常感谢。写的代码非常清晰。我看了下tutorials的mannualInjection,但是发现PositonFile里面没有记录直径的,只记录了位置信息,不知道mannualInjection是否支持变直径这种功能。
                                            感觉可能要深入这个源代码,看他如何初始化颗粒信息的

                                            Q 1 条回复 最后回复 回复 引用
                                            • S
                                              Samuel-Tu @李东岳 最后由 编辑

                                              @东岳 另外,我这个是湍流算出来的流场,我看ico这个solver初始化的时候,是用的laminar的mu,不知道是否能用。。

                                              1 条回复 最后回复 回复 引用
                                              • Q
                                                qingdong_wang @Samuel-Tu 最后由 编辑

                                                @Samuel-Tu 你可以看一下这个kinematicproperties

                                                injectionModels
                                                   {
                                                       model1
                                                       {
                                                           type            patchInjection;
                                                           massTotal       140e-3;
                                                           SOI             0;
                                                           parcelBasisType fixed;//mass;
                                                           nParticle       1;
                                                           patch           inlet;
                                                           duration        1;
                                                           parcelsPerSecond 1e5;
                                                           U0              (0 0 0.1);
                                                           flowRateProfile constant 1;
                                                           sizeDistribution
                                                           {
                                                               type        normal;
                                                               normalDistribution
                                                               {
                                                                   expectation 550e-6;
                                                                   variance 50e-6;
                                                                   minValue 400e-6;
                                                                   maxValue 800e-6;
                                                               }
                                                           }
                                                       }
                                                   }
                                                
                                                

                                                这个sizedistribution就是设置粒子直径, 这个是正态分布, 你可以看看其他设置方法, 看看有没有这种达到你想要的方法。
                                                @东岳 这个方法我试验过, 他的图层是按照次序自然叠加的吗? 比如说第一个case显示粒子, 第二个case显示alpha, 我就是这么做的, 但是alpha还是覆盖粒子。 是不是我应该第一个case显示alpha, 然后第二个case显示粒子, 这样就可以显示?因为我就是想要粒子显示在alpha上方。

                                                S 李东岳 2 条回复 最后回复 回复 引用
                                                • S
                                                  Samuel-Tu @qingdong_wang 最后由 编辑

                                                  @qingdong_wang 谢谢,我也发现了这个设置。说明在初始化粒子的时候,是可以对每个粒子分别设置的,我需要去研究一下粒子的构造函数是如何初始化和调用的。

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

                                                    @qingdong_wang 你加一些transparent就可以了 你试试

                                                    CFD高性能服务器 http://dyfluid.com/servers.html
                                                    2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

                                                    1 条回复 最后回复 回复 引用
                                                    • Q
                                                      qingdong_wang 最后由 编辑

                                                      @东岳 我试过, 但是我觉得比较丑,哈哈。 我看网上有人用python的paraview module进行图层的叠加。 我再看看那个。

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

                                                        捕获.JPG

                                                        这是我们放在文章里面的图,我们觉得还行 :140: :140: 我们把粒子显示在了U上面

                                                        CFD高性能服务器 http://dyfluid.com/servers.html
                                                        2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

                                                        1 条回复 最后回复 回复 引用
                                                        • S
                                                          Samuel-Tu 最后由 编辑

                                                          研究了下manualInjection的构造函数,直径和速度分别有构造函数,构造了List,说明可以分别设置,有了更系统的发现再跟大家汇报

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