本系列来看看 OpenFOAM 中的壁面函数。壁面函数的本质,是边界条件。这里主要来看看壁面函数的基本原理,OpenFOAM 中实现了的壁面函数,以及选择壁面函数的一些参考依据。
1. 壁面函数的基本原理
湍流模拟中,需要对近壁区域进行处理。一般来讲,壁面处理方法包含两类,一类是使用很细的网格,使靠近壁面的第一层网格在粘性层内($y^+ 30$),然后用经验公式来将粘性层和对数区关联起来。下图是一个典型的壁面附近的 $U^+ text{-} y^+$ 关系图。
图片来自 Wikipedia:Law of the wall。
在粘性层,满足如下关系
$$
u^+ = y^+
$$
而在对数区,则满足
$$
U^+ = frac{1}{kappa}ln(Ey^+)
$$
其中 $U^+ = U/u_tau$, $y^+ = yu_tau/nu$, $u_tau = sqrt{tau_w/rho}$,$kappaapprox 0.41$,$E approx 9.8$,$y$ 表示与壁面的距离。1$),然后里可以直接解析到粘性层的低雷诺湍流模型;另一类,不直接解析粘性层,而是将第一层网格设置在对数区($y^+>
本篇以标准壁面函数法来讨论一下壁面函数方法的基本原理,以及壁面函数在 OpenFOAM 中的实现。下面的讨论,先局限在 $k-varepsilon$ 模型,且第一层网格在对数区的情形。
先来看一下壁面函数方法需要解决什么问题。
有限体积方法中,扩散项的离散可以表示如下:
$$
nabla cdot (nu nabla U) = sum_f left [nu_f cdot (nabla U)_f right]
$$
当 $f$ 表示的是壁面边界单元时