CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    compressibleInterFoam求解器怎么加入空化模型

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

      http://cfd-china.com/topic/95/compressibletwophasemixture库编写问题 链接是一年前有前辈在论坛里问过的问题。我现在也面临着相似的问题,怎么编写compressibleTwoPhaseMixture库来调用空化模型?

      S 1 条回复 最后回复 回复 引用
      • 影
        影川风 最后由 李东岳 编辑

        我想在compressibleInterFoam求解器的基础上加入空化模型,控制方程的编写基本上没啥问题,不过参数调用这块我不怎么会,我请教过别人也看过曾经在这论坛里问过相似问题的前辈问的问题,不过还是不怎么理解。

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

          Okay。我看了一下,interPhaseChangeFoam只能用与不可压缩,所以你们打算添加可压缩VOF的方程;

          1. 一种方式是基于compressibleInterFoam附加空化模型;
          2. 一种方式是基于interPhaseChangeFoam更改速度压力方程适用于可压缩模型;

          怎么编写compressibleTwoPhaseMixture库来调用空化模型?

          这是一个非常大的问题,恐怕需要你自己完成。如果你选择第二种方法,你需要进行2步:

          1. 确认空化模型是否需要改动来可压缩适配,若不需要修改可跳过,若需要修改可修改phaseChangeTwoPhaseMixture.C公式内的函数;
          2. incompressibleTwoPhaseMixture.C这个文件只是计算相界面的粘度、运动粘度等,个人认为不需要修改,可跳过; 不过最好确认下;
          3. 修改VOF模型中的速度和压力方程,可参考compressibleInterFoam中的处理方式,主要在压力方程中附加源项;

          在做改动前,最后要非常熟悉interPhaseChangeFoam的流程,然后逐步的添加自己的模型。

          2023年线下CFD课 报名正式启动 http://dyfluid.com/class.html
          CFD高性能服务器 http://dyfluid.com/servers.html

          影 1 条回复 最后回复 回复 引用
          • 影
            影川风 @李东岳 最后由 编辑

            @李东岳 非常感谢东岳老师耐心的解答,我最近尝试着解决下这个问题

            Y 小 2 条回复 最后回复 回复 引用
            • Y
              yuan_neu @影川风 最后由 编辑

              第二种方案可行性更强
              可以参考这个

              Numerical Modelling of Evaporation and
              Condensation Phenomena

              论文里已经有了现成的compressible 的公式推导

              1 条回复 最后回复 回复 引用
              • 小
                小考拉 @影川风 最后由 编辑

                @影川风 前辈您好,我目前也尝试在compressibleInterFoam里面添加空化模型,可以和前辈交流一下吗?

                我 1 条回复 最后回复 回复 引用
                • 我
                  我爱梨涡 @小考拉 最后由 编辑

                  @小考拉 你好,加个好友?1195143370,我也是interfoam加空化

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

                    不知道这篇论文能否给点帮助(同学的论文):
                    "叶秉晟. 针对非定常空化流动的多相可压缩求解器开发及应用[D]. 中国科学院大学,2019."

                    我 1 条回复 最后回复 回复 引用
                    • 我
                      我爱梨涡 @上级 最后由 编辑

                      @上级 老哥,俺下载不了啊:chouchou: 如有时间,麻烦老哥发一份可否?1195143370@qq.com,拜谢!!!

                      1 条回复 最后回复 回复 引用
                      • S
                        shrine @影川风 最后由 李东岳 编辑

                        @影川风 我也在学习这个模型,请教:

                        1. 看到tutorials里面只有层流的例子,这个不能计算紊流吗?

                        2. thermophysicalProperties三个文件中

                        phases (water air);
                        
                          pMin 10000;
                        
                        sigma 0.07;
                        

                        这两个变量不知道什么意思

                        1. 还有水相的文件中
                        equationOfState
                        {
                        C (0.001278 -2.1055e-06 3.9689e-09 4.3772e-13 -2.0225e-16);
                        }
                        

                        这个也不知道是啥。

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

                          水相前面写的是:

                          equationOfState rPolynomial;
                          

                          所以是多项式状态方程,在rPolynomial.H里面写的是:

                          Reciprocal polynomial equation of state for liquids and solids
                          
                              \f[
                                  1/\rho = C_0 + C_1 T + C_2 T^2 - C_3 p - C_4 p T
                              \f]
                          
                              This polynomial for the reciprocal of the density provides a much better fit
                              than the equivalent polynomial for the density and has the advantage that it
                              support coefficient mixing to support liquid and solid mixtures in an
                              efficient manner.
                          

                          所以就是:
                          \begin{equation}
                          1/\rho = C_0 + C_1 T + C_2 T^2 - C_3 p - C_4 p T
                          \end{equation}

                          2023年线下CFD课 报名正式启动 http://dyfluid.com/class.html
                          CFD高性能服务器 http://dyfluid.com/servers.html

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

                            @李东岳 在 compressibleInterFoam求解器怎么加入空化模型 中说:

                            多项式状态方程

                            这是固液两项?如果算气液两相这么设?

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

                              @shrine 什么都可以这么设定。主要是密度算的对就行。

                              2023年线下CFD课 报名正式启动 http://dyfluid.com/class.html
                              CFD高性能服务器 http://dyfluid.com/servers.html

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