Fluent提供了众多的近壁面处理方法,这里简单扒一扒。 在临近壁面位置,法向速度非常大的梯度。在非常小的壁面法向距离内,速度从相对较大的值下降到与壁面速度相同。因此对于该区域内流场的计算,通常采用两种方式:(1)利用壁面函数法;(2)加密网格,利用壁面模型法。对于这两类方法的选取,可以通过对于这两类方法的选取,可以通过y+来体现。 如图所示为近壁面位置无量纲速度分布情况。 图中横坐标所表示为无量纲壁面距离y+,纵坐标为无量纲速度u+。其中: 从图中看出,近壁面区域可分为3个区间:
对于近壁区域求解,主要集中在粘性子层的求解上,主要有两种方式: 1、利用壁面模型直接求解粘性子层 若想要求解粘性子层,则需要保证y+值小于1(建议接近1)。由于y+直接影响第一层网格节点位置,因此对于求解粘性子层的情况,需要非常细密的网格,通常要求有10~20层边界层网格。 2、利用壁面函数近似处理。对于湍流模型,需要选择低雷诺数湍流模型(如k-omega模型)。通常来说,若壁面对于仿真结果非常重要(如气动阻力计算、旋转机械叶片性能等),则需要采用此类方法。 壁面函数要求第一层网格尺寸满足条件30<y+<300,当尺寸过小时,壁面函数不可用;当尺寸超出该范围时,无法求解粘性子层。通常使用高雷诺数湍流模型(如标准k-epsilon模型、Realizable k-epsilon模型、RNG k-epsilon模型等)。一般来说,在粘性子层数据不是特别重要的时候可以选用壁面函数进行求解。 如下图所示分别为壁面函数与壁面模型在近壁面区域对网格的需求。 Fluent中壁面函数主要应用于k-epsilon及雷诺应力模型。 如下图所示为k-epsilon湍流模型中的壁面函数。 雷诺应力模型中的壁面函数如下图所示。 其他湍流模型均不支持壁面函数法,因此需要划分更加细密的近壁面网格,通常要求满足y+接近于1。 这些近壁面处理方法包括:标准壁面函数(Standard Wall Functions)、可缩放壁面函数(Scalable Wall Functions)、非平衡壁面函数(Non-Equilibrium Wall Functions)、增强壁面处理(Enhanced Wall Treatment)、Menter-Lechner壁面函数以及自定义壁面函数(User-Defined Wall Functions)。 这些壁面处理方式中,标准壁面函数、可缩放壁面函数以及非平衡壁面函数均为壁面函数法,适合于高雷诺数湍流模型(k-epsilon模型以及雷诺应力模型),其要求第一层网格节点处于湍流核心区域,即y+值处于 而近壁面处理则并非壁面函数法,其适合于低雷诺数湍流模型(k-omega模型),需要在近壁区域划分足够细密的网格,其要求第一层网格节点位于粘性子层内,即y+<5,且要求边界网格层数至少为10~15层。 虽然壁面函数法是一种近似处理方法,然而其在工业流动问题计算中仍然应用非常广泛。对于简单的剪切流动问题,利用标准壁面函数法可以很好的得到解决,而使用非平衡壁面函数法可以对于强压力梯度及分离流动计算进行改善。而可缩放的壁面函数法则可改善第一层网格节点在计算迭代过程中处于粘性子层与核心层之间摇摆从而导致计算不稳定的问题。增强壁面处理通常用于无法应用对数律的复杂流动问题。 近壁面建模的一些推荐策略:
|
|