OpenFOAM 中边界条件的设定【转载】

转载自：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/

 

最终我选择的边界条件是：