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为特征尺寸)。