Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    interFoam下VOF 跟 LS 的耦合问题

    OpenFOAM
    4
    7
    3267
    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.
    • 影
      影川风 last edited by

      参照网上的一个教程修改Coupled Level-Set with VOF interFoam

      其中有一个点是:

      dimensionedScalar deltaX
          (
              twoPhaseProperties.lookup("deltaX")
          );
      

      所给的案列中:
      网格是:

      convertToMeters 1;
      18
      19 vertices
      20 (
      21 (0 0 0)
      22 (1 0 0)
      23 (1 2 0)
      24 (0 2 0)
      25 (0 0 0.1)
      26 (1 0 0.1)
      27 (1 2 0.1)
      28 (0 2 0.1)
      29 );
      30
      31 blocks
      32 (
      33 hex (0 1 2 3 4 5 6 7) (160 320 1) simpleGrading (1 1 1)
      34 );
      

      传输特征transportProperties ;

      deltaX  deltaX[ 0 0 0 0 0 0 0 ] 0.00625; // 0.006667;
      epsilon epsilon [ 0 0 0 0 0 0 0 ] 0.009375; // 1.5* deltaX ;
      

      能够发现:

      deltaX=1/160=0.00625.
      

      我的问题是这个CLSVOF耦合似乎只针对均匀网格的,而我想运用在不均匀的网格里,我能怎么修改这部分呢?

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

        不是很熟悉你提到的这个CLSVOF。

        但是我在高阶格式里面MUSCL里面也需要类似的操作。在我这个里面,你那个

        deltaX=1/160=0.00625.
        

        在笛卡尔网格里面比较好处理,如果非均匀网格,deltaX就是网格中心点到面的距离,比如

        deltaX=mag(mesh.C(cellI).component(0) - mesh.Cf(faceI).component(0)) 
        

        要是非笛卡尔网格,需要转换一下。目前这块我还没处理。我的MUSCL供参考,不知道是否有相通之处。

        CFD课程 改成线上了 http://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        M 影 3 Replies Last reply Reply Quote
        • M
          mohui @李东岳 last edited by

          @李东岳 应该是可以处理的。

          1 Reply Last reply Reply Quote
          • 影
            影川风 @李东岳 last edited by

            @李东岳 在 interFoam下VOF 跟 LS 的耦合问题 中说:

            deltaX

            谢谢东岳老师,我目前在按照你所给思路往下弄。

            1 Reply Last reply Reply Quote
            • 影
              影川风 @李东岳 last edited by 李东岳

              @李东岳 我简单画了一个计算单元,如图所示

              0_1516353981487_单元图.png

              输出体心(CX)与面心值(CX1)如下:

              mesh = CX 4(0.0025 0.0075 0.0025 0.0075)
              
              mesh = CX1 4(0.005 0.0025 0.0075 0.005)
              

              体心还明白是怎么回事,面心有点不明白为什么是这个值?

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

                体心就是网格单元体中心的值..
                面心就是网格单元面中心的值..

                CFD课程 改成线上了 http://dyfluid.com/class.html
                CFD高性能服务器 http://dyfluid.com/servers.html

                1 Reply Last reply Reply Quote
                • S
                  sungda @影川风 last edited by

                  @影川风 您好,我现在也在做这个耦合的内容,您这边这个deltaX最后有什么好的处理办法吗

                  1 Reply Last reply Reply Quote
                  • First post
                    Last post

                  CFD中文网 | 东岳流体 | 京ICP备15017992号-2