1. 程式人生 > 其它 >OpenFOAM中的幾種邊界條件

OpenFOAM中的幾種邊界條件

1.turbulentIntensityKineticEnergyInlet

邊界名稱
{
   type              turbulentIntensityKineticEnergyInlet;
   intensity         數值;      
   value             uniform 數值/$internalField;
 }

說明:

intensity指的是湍流強度,如果不知道怎麼計算,便可以指定為0.05。關於為什麼指定為5%,可以參考文獻《The numerical computation of turbulent flows》。可通過湍流強度來k的值,value關鍵字下可填寫任意數值或$internalField(僅是將字串拷貝到本地,等於內部場的值),僅起到佔位的作用,並不對計算造成影響。OpenFOAM自帶tutorial中,基本上都用在了k的進口邊界上。很顯然在壁面處,k應該為0。

2.omegaWallFunction/epsilonWallFunction

邊界名稱
{
    type               omegaWallFunction;//或者是epsilonWallFunction
    value              uniform 數值;
}

說明:

該邊界條件僅對壁面設定壁面函式,應用在k-e或者k-w湍流模型中,從而對方程進行求解。value關鍵字下可填寫任意數值或$internalField(僅是將字串拷貝到本地,等於內部場的值),僅起到佔位的作用,並不對計算造成影響。

3.turbulentMixingLengthFrequencyInlet/turbulentMixingLengthDissipationRateInlet

邊界名稱
{
    type                  turbulentMixingLengthFrequencyInlet;//或者turbulentMixingLengthDissipationRateInlet
    mixingLength          數值;
    value                 uniform 數值;
}

說明:

可通過混合長度來計算omega或者epsilon的值,value關鍵字下可填寫任意數值或$internalField(僅是將字串拷貝到本地,等於內部場的值),僅起到佔位的作用,並不對計算造成影響。mixingLength為湍流尺度,其計算公式為:I=0.07L(其中I為mixingLength的值,L為特徵尺寸)。