1. Home
  2. OpenFOAM
  3. interFoam下VOF 跟 LS 的耦合问题

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供参考，不知道是否有相通之处。

    M 影 3 Replies
  • M

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

  • 影

    @李东岳 在 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)

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

  • 李东岳

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

Log in to reply