Question

我运行multiphaseInterfoam，我的入口非常不稳定（我希望它不变。）

以下是我的alpha文件\

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

FoamFile {       版本2.0;       格式ascii;       class volScalarField;       位置“0”;       对象alpha.air; } * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField {       // - 为约束补丁设置patchGroup       #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;
} }

alpha.freshwater：

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

FoamFile {       版本2.0;       格式ascii;       class volScalarField;       位置“0”;       object alpha.freshWater; } * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField {       // - 为约束补丁设置patchGroup       #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.saltWater

FoamFile {       版本2.0;       格式ascii;       class volScalarField;       位置“0”;       object alpha.saltWater; } * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField {       // - 为约束补丁设置patchGroup       #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;       class volScalarField;       位置“0”;       对象alphas; } * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField {       // - 为约束补丁设置patchGroup       #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给出了以下所需的流体分布

enter image description here

然而，经过几次上述变化，也在入口处：

enter image description here

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

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

FoamFile {       版本2.0;       格式ascii;       class volScalarField;       位置“0”;       对象alpha.air; } * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField {       // - 为约束补丁设置patchGroup       #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，我得到一条很长的错误消息，其中包含被解释为互联网链接的内容，因此我不允许在这里删除它们。可以在this link to the CDF-online forum上看到错误消息，我也在这里提出了这个问题。

anybvody知道如何修改上述文件以使入口处的alpha分布不变吗？

Answer 1

在alpa.air的最新版本中，我改变了＆＃34; internalField uniform 0＆＃34; to＆＃34; internalField uniform 1＆＃34;，现在入口是恒定的。我认为internalFIeld值必须等于setFieldsDict中设置的默认值。

Openfoam，multiphaseInterfoam，非恒定入口

1 个答案: