Openfoam，多相Interfoam，非常量入口

我运行 multiphaseInterfoam，但入口不恒定（我希望它是恒定的）有问题。

这是我的 alpha 文件\

/ -------------------------------- - C++ - -------------- -------------------- \

FoamFile { 2.0 版；ASCII格式；类 volScalarField; 位置“0”；对象 alpha.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

尺寸 [0 0 0 0 0 0 0];

internalField 统一 0;

boundaryField { //- 为约束补丁设置补丁组 #includeEtc "caseDicts/setConstraintTypes"

inlet
{
    type            alphaContactAngle;
    thetaProperties
    (
        ( freshWater air ) 90 0 0 0
        ( saltWater air ) 90 0 0 0
        ( freshWater saltWater ) 90 0 0 0
    );
    value           uniform 0;
}

outlet

{
    type            alphaContactAngle;
    thetaProperties
    (
        ( freshWater air ) 90 0 0 0
        ( saltWater air ) 90 0 0 0
        ( freshWater saltWater ) 90 0 0 0
    );
    value           uniform 0;

}

    atmosphere
{
    type            inletOutlet;
    inletValue      uniform 1;
    value           uniform 1;
}


barge
{
    type            alphaContactAngle;
    thetaProperties
    (
        ( freshWater air ) 90 0 0 0
        ( saltWater air ) 90 0 0 0
        ( freshWater saltWater ) 90 0 0 0
    );
    value           uniform 0;
} }

阿尔法淡水：

/ -------------------------------- - C++ - -------------- -------------------- \

FoamFile { 2.0 版；ASCII格式；类 volScalarField; 位置“0”；对象 alpha.freshWater; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

尺寸 [0 0 0 0 0 0 0];

internalField 统一 0;

boundaryField { //- 为约束补丁设置补丁组 #includeEtc "caseDicts/setConstraintTypes"

inlet
{
    type            fixedValue;
    value           $internalField;
}

outlet
{
    type            variableHeightFlowRate;
    lowerBound      0;
    upperBound      1;
    value           $internalField;
}

atmosphere
{
    type            inletOutlet;
    inletValue      $internalField;
    value           $internalField;
}

barge
{
    type            zeroGradient;
} }

阿尔法盐水

FoamFile { 2.0 版；ASCII格式；类 volScalarField; 位置“0”；对象 alpha.saltWater; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

尺寸 [0 0 0 0 0 0 0];

internalField 统一 0;

boundaryField { //- 为约束补丁设置补丁组 #includeEtc "caseDicts/setConstraintTypes"

inlet
{
    type            fixedValue;
    value           $internalField;
}

outlet
{
    type            variableHeightFlowRate;
    lowerBound      0;
    upperBound      1;
    value           $internalField;
}

atmosphere
{
    type            inletOutlet;
    inletValue      $internalField;
    value           $internalField;
}

barge
{
    type            zeroGradient;
} }

// **************************************************** **************************** //

阿尔法

/ -------------------------------- - C++ - -------------- -------------------- \

FoamFile { 2.0 版；ASCII格式；类 volScalarField; 位置“0”；对象阿尔法；} // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

尺寸 [0 0 0 0 0 0 0];

internalField 统一 0;

boundaryField { //- 为约束补丁设置补丁组 #includeEtc "caseDicts/setConstraintTypes"

inlet
{
    type            fixedValue;
    value           $internalField;
}

outlet
{
    type            variableHeightFlowRate;
    lowerBound      0;
    upperBound      1;
    value           $internalField;
}

atmosphere
{
    type            inletOutlet;
    inletValue      $internalField;
    value           $internalField;
}

barge
{
    type            zeroGradient;
} }

// **************************************************** **************************** //

以上给出了时间步长 1 所需的以下流体分布

在此处输入图片说明

但是，经过多次上述更改后，同样在入口处：

在此处输入图片说明

我真的不明白上面 alpha.air 中使用的 contactAngle functino。我曾尝试使用以下 alpha.air

/ -------------------------------- - C++ - -------------- -------------------- \

FoamFile { 2.0 版；ASCII格式；类 volScalarField; 位置“0”；对象 alpha.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

尺寸 [0 0 0 0 0 0 0];

internalField 统一 0;

boundaryField { //- 为约束补丁设置补丁组 #includeEtc "caseDicts/setConstraintTypes"

inlet
{
    type            fixedValue;
    value           $internalField;
}

outlet
{
    type            variableHeightFlowRate;
    lowerBound      0;
    upperBound      1;
    value           $internalField;
}

atmosphere
{
    type            inletOutlet;
    inletValue      $internalField;
    value           $internalField;
}

barge
{
    type            zeroGradient;
} }

// **************************************************** **************************** //

使用上面的 alpha.air，我收到一条很长的错误消息，其中包含被解释为 Internet 链接的内容，因此我不允许在此处发布它们。错误消息可以在指向 CDF 在线论坛的链接上看到，我也在那里问过这个问题。

有没有人知道如何修改上述文件以使 alpha 在入口常数处的分布？