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

mohui

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

影川风

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

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

李东岳

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

sungda

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