是否可以同时定义边界alpha和接触角？

[yang]

各位大神，
我想用interFoam模拟一个固体表面上的一个静止液滴。想同时固定液滴和表面接触线以及接触角是否能够实现？需要对标准算法或边界做如何修改？

万分感激！

[李东岳]

参考一下interFoam求解器下面的damBreak下的设置。

[yang]

@东岳 多谢东岳大神回复。之前可能我没表述清楚我的问题：我想同时fix液滴与固体表面的接触线和接触角，即壁面处的volume fraction和theta都始终保持t=0时的值，并不随时间而改变。相当于fixed value和constantAlphaContactAngle两个边界条件的结合。如果边界条件设置为constantAlphaContactAngle(zeroGradient也相当于对应theta0=90的constantAlphaContactAngle边界）, 在迭代的过程中会计算当前的curvature, 并通过与设置的theta0比较来更新壁面处的受力情况，从而重新计算边界处的volume fraction, 以使计算的theta趋近于所定义的theta0（如后面摘自interfaceProperties中的代码所示），这并不能保证接触线（即volume fraction）的固定。
我的理解是这样的，不知道对不对：接触角theta与壁面处volume fraction的梯度是相对应的。能否让壁面处的第一层和第二层网格中的volume fraction值都固定，即volume fraction方程在该区域不进行求解而是始终采用t=0时刻的初始值，这样就能实现我前面所说的目的了。不知在openfoam中能否实现。
期待您答复，谢谢！

   forAll(boundary, patchi)
    {
        if (isA<alphaContactAngleFvPatchScalarField>(abf[patchi]))
        {
            alphaContactAngleFvPatchScalarField& acap =
                const_cast<alphaContactAngleFvPatchScalarField&>
                (
                    refCast<const alphaContactAngleFvPatchScalarField>
                    (
                        abf[patchi]
                    )
                );

            fvsPatchVectorField& nHatp = nHatb[patchi];
            const scalarField theta
            (
                convertToRad*acap.theta(U_.boundaryField()[patchi], nHatp)
            );

            const vectorField nf
            (
                boundary[patchi].nf()
            );

            // Reset nHatp to correspond to the contact angle

            const scalarField a12(nHatp & nf);
            const scalarField b1(cos(theta));

            scalarField b2(nHatp.size());
            forAll(b2, facei)
            {
                b2[facei] = cos(acos(a12[facei]) - theta[facei]);
            }

            const scalarField det(1.0 - a12*a12);

            scalarField a((b1 - a12*b2)/det);
            scalarField b((b2 - a12*b1)/det);

            nHatp = a*nf + b*nHatp;
            nHatp /= (mag(nHatp) + deltaN_.value());

            acap.gradient() = (nf & nHatp)*mag(gradAlphaf[patchi]);
            acap.evaluate();
        }
    }