OpenFOAM中的幾種邊界條件
阿新 • • 發佈:2021-08-09
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為特徵尺寸)。