转载自:http://blog.sina.com.cn/s/blog_a0b4201d0102v7jt.html
用习惯了FLUENT的操作界面,再使用OpenFOAM就会觉得非常繁琐。遇到的第一个问题就是边界条件的设定。在FLUENT中,例如给定了一个进口边界(质量或者压力或者速度)和一个出口边界(一般是压力),并在进口处设定湍流强度就可以了,不必关心这三个量在求解的各个方程中的分配。而在OpenFOAM中,这个问题就要显得复杂的多,以k-epsl模型求解可压缩流为例为例,要求解质量守恒方程,动量守恒方程和能量守恒方程,湍动能方程以及耗散率方程,分别对应P方程,U方程,T方程,k方程和epsl方程,其中U方程是矢量方程,包含三个标量方程。在求解之前要为每个方程指定边界条件。现在的问题就是:只知道压力进口条件和压力出口回流条件,以及进出口都是室温条件,也就是说只能显式的给定P方程的进口条件和出口回流条件,T方程的进口和出口条件,而U方程的边界条件需要由已知的P方程和T方程的边界条件间接的得到。对于这个问题,OpenFOAM提供了各种各样的边界条件给定方式,多的足够让你眼花缭乱,所以非常有必要深入的扒一扒。先从OpenFOAM在tutorial中给出的例子入手。
以rhoSimplecFoam中的squareBend为例,给出了alphat,epsilon, k, mut, p, T, U 等7个量的初始条件,当然也包括了边界条件。其中alphat指的是turbulence thermal
diffusivity, 例子中给定的壁面边界是
compressible::alphatWallFunction; 进口边界是calculated, value uniform 0;
出口边界是 calculated, value uniform 0; epsilon指的是turbulence
kinetic energy dissipation
rate, k指的是turbulence
kinetic energy,mut指的是turbulence
viscosity,同时alphat, epsilon, k, mut
这个几个参数都是要给定壁面条件模型的。像其他的tutorial例子中给出的也都是速度进口边界条件,跟我的算例所需要的压力进口边界条件不太吻合。
好在OpenFoam给出了各种不同类型的边界条件类型供大家挑选使用,(当然是在不行就只能自己写了,这才是OF最强大的的地方嘛)。具体的不同内容可以查询页面: http://openfoam.org/docs/cpp/
最终我选择的边界条件是: