-1

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

这是我的 alpha 文件\

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

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

尺寸[0 0 0 0 0 0 0]；

内部字段统一 0；

boundaryField { //- 为约束补丁设置 patchGroups #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]；

内部字段统一 0；

boundaryField { //- 为约束补丁设置 patchGroups #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格式；类 volScalarField；位置“0”；对象 alpha.saltWater; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

尺寸[0 0 0 0 0 0 0]；

内部字段统一 0；

boundaryField { //- 为约束补丁设置 patchGroups #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]；

内部字段统一 0；

boundaryField { //- 为约束补丁设置 patchGroups #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 函数。我尝试过以下 alpha.air

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

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

尺寸[0 0 0 0 0 0 0]；

内部字段统一 0；

boundaryField { //- 为约束补丁设置 patchGroups #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-online forum 的这个链接上看到，我也在那里问过这个问题。

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

4

1 回答 1

0

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

于 2018-02-20T09:46:18.960 回答