CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    MPPICInterFoam中静置情况下水流沸腾

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

      如题,在图中底下的凹槽里填满里粒子,红色为水域,蓝色为空气。在静置的状态下,水流会像开锅一样沸腾。有遇到过类似问题的同学吗?希望能指出问题所在,谢谢。

      1.png

      2.png

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

        @zhe 啥也不设置自己沸腾?压力?

        m.sui20@foxmail.com

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

          @星星星星晴 对,什么设置都没有。观察了流速和压强状况,在粒子堆与纯水的交界处最强烈。那个地方的不同就是由于粒子的存在,alphac会出现梯度。但是在动量方程里又找不到这一项东西。有没有想法,有几个猜想也是好的。先感谢!

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

            @zhe 是不是粒子到了packing limit了

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

            1 条回复 最后回复 回复 引用
            • T
              Tens 讲师 最后由 编辑

              颗粒和水的密度比大吗,是不是有可能是颗粒在浮力作用下运动引起的波动呢

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

                @星星星星晴 @李东岳 @Tens 找到原因了。就像之前观测发现的现象,在粒子群和纯水的界限处流速和压强极强。问题是出在压力泊松方程上。MPPICInterFoam应该按着interFoam的连续方程,考虑alphac的影响。但目前MPPICInterFoam是按照单项流的DPMFoam的方程。ddt(alphac)这一项是造成水沸腾的原因。把现有的

                    surfaceScalarField Dp("Dp", alphacf*rAUcf);
                
                    fvScalarMatrix p_rghEqn
                    (
                        fvm::laplacian(Dp, p_rgh)
                     ==
                        fvc::ddt(alphac) + fvc::div(alphacf*phiHbyA)
                    );
                

                改成:

                    surfaceScalarField Dp("Dp", alphacf*rAUcf);
                
                    fvScalarMatrix p_rghEqn
                    (
                        fvm::laplacian(alphacf*Dp, p_rgh) == fvc::div(alphacf*phiHbyA)
                    );
                

                这样就不会出现“无风自动”的现象了。是不是对别的地方会造成影响,还需要再研究。(其中出现的alphac*alphac是因为我的粒子浓度太高,不能忽略压力项。李东岳老师在其他帖子里也讨论过这个问题。根据自己的粒子粒径或占网格比来定。

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

                  @zhe :146: :146:

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

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