multiphaseInterfoam を実行していますが、インレットが一定でないことに問題があります (一定にしたいのです)。
ここに私のアルファファイルがあります\
/ -------------------------------- - C++ - -------------- -------------------- \
FoamFile { バージョン 2.0; ascii をフォーマットします。クラスvolScalarField; ロケーション "0"; オブジェクト alpha.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
次元 [0 0 0 0 0 0 0];
internalField 均一 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; } }
alpha.淡水:
/ -------------------------------- - C++ - -------------- -------------------- \
FoamFile { バージョン 2.0; ascii をフォーマットします。クラスvolScalarField; ロケーション "0"; オブジェクト alpha.freshWater; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
次元 [0 0 0 0 0 0 0];
internalField 均一 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];
internalField 均一 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];
internalField 均一 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];
internalField 均一 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 を使用すると、インターネット リンクとして解釈されるものを含む長いエラー メッセージが表示されるため、ここで公開することはできません。エラー メッセージは、CDF オンライン フォーラムへのこのリンクで見ることができます。ここでもこの質問をしました。
上記のファイルを変更して、入口でのアルファの分布を一定にする方法を知っている人はいますか?