interFoam下VOF 跟 LS 的耦合问题



  • 参照网上的一个教程修改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耦合似乎只针对均匀网格的,而我想运用在不均匀的网格里,我能怎么修改这部分呢?


  • 管理员

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

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

    deltaX=1/160=0.00625.
    

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

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

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



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



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

    deltaX

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



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

    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)
    

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


  • 管理员

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


登录后回复
 

与 CFD中文网 的连接断开,我们正在尝试重连,请耐心等待