Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    reactingTwoPhaseEulerFoam求解器的组分输运方程中添加质量源相导致质量分数无界

    OpenFOAM
    3
    15
    2442
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • Z
      zwl last edited by

      reactingTwoPhaseEulerFoam求解器的组分输运方程中添加质量源相导致质量分数无界,这是什么原因,为什么没添加源相之前有界呢?如何处理才能避免质量分数远大于1?

      1 Reply Last reply Reply Quote
      • 李东岳
        李东岳 管理员 last edited by

        是的,添加源项会导致越界。$\frac{\p \phi}{\p t}+\bfU\cdot\nabla\phi=S$,只有在$S=0$的时候才有界。目前只能帮你到这里了

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

        Z 1 Reply Last reply Reply Quote
        • Z
          zwl @李东岳 last edited by

          @东岳 非常感谢东岳老师的回答,我有个疑问,这个求解器是可以算蒸发的,而且用的就是界面组分,他们是怎么处理的呢?就避免了越界。

          李东岳 1 Reply Last reply Reply Quote
          • 李东岳
            李东岳 管理员 @zwl last edited by

            @zwl 你是说植入的Y方程么

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

            Z 1 Reply Last reply Reply Quote
            • Z
              zwl @李东岳 last edited by

              @东岳 对,就是自带的Y方程。

              1 Reply Last reply Reply Quote
              • 李东岳
                李东岳 管理员 last edited by

                Y方程没有详细看过。粗略看一眼并没有用MULES求解。另外,带源项要有界可以参考可压缩的相方程

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

                Z 1 Reply Last reply Reply Quote
                • Y
                  yfclark 讲师 last edited by

                  @zwl 会不会是源项处理有问题,求解不稳定,Interface Composition搞了个fvm::Sp:
                  c86d36d2-987e-433f-b3ad-2d3f527a8772-image.png
                  此外可以考虑更换对流项格式,Laneer,Gamma01等

                  Z 1 Reply Last reply Reply Quote
                  • Z
                    zwl @yfclark last edited by

                    @yfclark 非常感谢你的建议。源相我确实是直接添加的,没有再做处理,这个Sp会有什么作用(虽然有时我也会模仿着这样处理)?越界的主要原因是什么?

                    1 Reply Last reply Reply Quote
                    • Z
                      zwl @李东岳 last edited by

                      @东岳 那么下面带拉普拉斯项的组分输运方程如何采用MULES求解呢?谢谢。

                              fvm::ddt(alpha, rho, Yi)
                            + fvm::div(alphaRhoPhi, Yi, "div(" + alphaRhoPhi.name() + ",Yi)")
                      
                            - fvm::laplacian
                              (
                                  fvc::interpolate(alpha)
                                 *fvc::interpolate(this->muEff()/Sct_),
                                  Yi
                              )
                           ==
                              alpha*this->R(Yi)
                             + Source
                            + fvc::ddt(residualAlpha_*rho, Yi)
                            - fvm::ddt(residualAlpha_*rho, Yi)
                      
                      
                      1 Reply Last reply Reply Quote
                      • 李东岳
                        李东岳 管理员 last edited by

                        扩散项+对流项用高阶/低阶格式处理,通过反扩散算法做出通量,源项进Sp或者Su就行

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

                        Z 1 Reply Last reply Reply Quote
                        • Z
                          zwl @李东岳 last edited by

                          @东岳 我看MULES代码,感觉他就是为求解相方程定制的,没有扩散相的处理,所以我不知道这里的扩散相放在哪个位置,怎么处理。

                          李东岳 2 Replies Last reply Reply Quote
                          • 李东岳
                            李东岳 管理员 @zwl last edited by

                            @zwl 扩散项、梯度项、对流项都是通过高斯定律化为$\sum$加和的形式,都是一样的,以梯度项举例:

                            捕获.JPG

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

                            1 Reply Last reply Reply Quote
                            • 李东岳
                              李东岳 管理员 last edited by

                              对了,扩散项不会引起越界问题,或许你都不用处理了,扩散项是好朋友

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

                              1 Reply Last reply Reply Quote
                              • 李东岳
                                李东岳 管理员 @zwl last edited by

                                @zwl 你在研究E-E模型下面的组分传输么

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

                                1 Reply Last reply Reply Quote
                                • Z
                                  zwl last edited by

                                  是的,求解器自带的是蒸发,我想计算空化,二者的计算原理大不相同,我更改了质量源相,导致质量分数越界。

                                  1 Reply Last reply Reply Quote
                                  • First post
                                    Last post

                                  CFD中文网 | 东岳流体 | 京ICP备15017992号-2
                                  论坛登录问题反馈可联系 li.dy@dyfluid.com