Openfoam，multiphaseInterfoam，非恒定入口答案

【问题标题】：Openfoam, multiphaseInterfoam, non-constant inletOpenfoam，multiphaseInterfoam，非恒定入口
【发布时间】：2018-02-20 08:15:51
【问题描述】：

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

这是我的 alpha 文件\

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

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

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

internalField 统一 0;

边界字段 { //- 为约束补丁设置补丁组 #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++ -----------------------------------\

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

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

internalField 统一 0;

边界字段 { //- 为约束补丁设置补丁组 #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

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

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

internalField 统一 0;

边界字段 { //- 为约束补丁设置补丁组 #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++ -----------------------------------\

泡沫文件 { 2.0版； ASCII格式；类 volScalarField；位置“0”；对象阿尔法； } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField 统一 0;

边界字段 { //- 为约束补丁设置补丁组 #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 函数。我尝试过使用以下 alpha.air

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

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

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

internalField 统一 0;

边界字段 { //- 为约束补丁设置补丁组 #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 链接的内容，因此我不允许在此处发布它们。可以在this link to the CDF-online forum 上看到错误消息，我也曾在此提出过这个问题。

是否有人知道如何修改上述文件以使入口处的 alpha 分布保持不变？

【问题讨论】：

标签： openfoam

【解决方案1】：

在 alpa.air 的最后一个版本中，我从“internalField uniform 0”更改为“internalField uniform 1”，现在入口是不变的。我认为 internalFIeld 值必须等于 setFieldsDict 中设置的默认值。

【讨论】：