void Foam::interfaceProperties::calculateK()
{
    const fvMesh& mesh = alpha1_.mesh();
    const surfaceVectorField& Sf = mesh.Sf();

    // Cell gradient of alpha
    const volVectorField gradAlpha(fvc::grad(alpha1_, "nHat"));

    // Interpolated face-gradient of alpha
    surfaceVectorField gradAlphaf(fvc::interpolate(gradAlpha));

    // gradAlphaf -=
    //    (mesh.Sf()/mesh.magSf())
    //   *(fvc::snGrad(alpha1_) - (mesh.Sf() & gradAlphaf)/mesh.magSf());

    // Face unit interface normal
    surfaceVectorField nHatfv(gradAlphaf/(mag(gradAlphaf) + deltaN_));
    // surfaceVectorField nHatfv
    // (
    //     (gradAlphaf + deltaN_*vector(0, 0, 1)
    //    *sign(gradAlphaf.component(vector::Z)))/(mag(gradAlphaf) + deltaN_)
    // );
    correctContactAngle(nHatfv.boundaryFieldRef(), gradAlphaf.boundaryField());

    // Face unit interface normal flux
    nHatf_ = nHatfv & Sf;

    // Simple expression for curvature
    K_ = -fvc::div(nHatf_);

    // Complex expression for curvature.
    // Correction is formally zero but numerically non-zero.
    /*
    volVectorField nHat(gradAlpha/(mag(gradAlpha) + deltaN_));
    forAll(nHat.boundaryField(), patchi)
    {
        nHat.boundaryField()[patchi] = nHatfv.boundaryField()[patchi];
    }

    K_ = -fvc::div(nHatf_) + (nHat & fvc::grad(nHatfv) & nHat);
    */
}

interFoam在计算曲率的时候应该是调用了这个函数， 我想问一下这里被注释掉的这段complex expression什么时候可以用？ 有没有原始文献可以参考。

你现在解决虚拟流动的问题了吗？

@东岳
我现在是对所有网格做遍历，遍历每个网格上的面，判断是内部的面还是边界的面，如果是边界的面的话，通量值就会和另一块processor上对应的网格上的值相关，目前就是不知道怎么去找这个网格。

@东岳
谢谢李老师回复，我的问题是这样的：我想要遍历每个面修改通量，然后发现没有考虑网格内部的由于并行出现的每个processor上的交界面，现在不知道怎么去遍历这些边界上的面，以及找到这个边界的某个面的上下游的cell。有没有现成的例子可以参考。
不知道我有没有把问题描述清楚。

我想请教一个问题，在并行下，不同processor数据交换位置的面，patchi和facei一样吗？
还有就是如何去找一个processor上一个面对应的另外一个processor上的面。

@东岳
第一个大概明白了，就是用迎风格式的时候
-fvc::flux(-phir, alpha2, alpharScheme) 和 fvc::flux(phir, alpha2, alpharScheme) 会不一样对吧？

第二个的话

tmp<surfaceScalarField> talphaPhi1Corr(talphaPhi1Un() - alphaPhi10);

那alphaPhi10是高阶的吗？我之前以为是talphaPhi1Un是高阶的，alphaPhi10是低阶的。

         tmp<surfaceScalarField> talphaPhi1Un
         (
             fvc::flux
             (
                 phiCN(),
                 cnCoeff*alpha1 + (1.0 - cnCoeff)*alpha1.oldTime(),
                 alphaScheme
             )
           + fvc::flux
             (
                -fvc::flux(-phir, alpha2, alpharScheme),
                 alpha1,
                 alpharScheme
             )
         );

第二项通量中的-fvc::flux(-phir, alpha2, alpharScheme)为什么要用两个负号？

             MULES::correct
             (
                 geometricOneField(),
                 alpha1,
                 talphaPhi1Un(),
                 talphaPhi1Corr.ref(),
                 Sp,
                 (-Sp*alpha1)(),
                 1,
                 0
             );

调用mules的correct函数的时候，这个talphaPhi1Un()有什么用吗？ 我看了一下好像correct函数本身没有用到这个参数？

你换个并行方式试试

@mr-shan 跟库朗特数有关， 越小的网格精度越高，但是要保证一样的库朗特数就要用越小的时间步长。

我在fluent里面用CLSVOF的时候也碰到类似的情况， 原因我感觉是网格太小或者时间步长太大，你可以试试缩小时间步长。