最近总有网友问关于Y+的问题,这里来做一个总结。参考自cfd-online
wiki。
1、为什么会存在第一层网格的问题
关于这个问题,实际上要从边界层说起。实验表明,边界层内根据流动状态的不同可以分为三层,自壁面向流动核心区分别为:粘性子层、过渡层和湍流核心层。边界层很薄,一般都是毫米~微米级,因此,若采用划分网格进而利用数值方法求解的话,势必会大大增加计算网格的数量,从而急剧增加计算工作量。又有实验发现,在粘性子层和过渡层内,主要是粘性力在起主导作用,惯性力的作用几乎可以忽略。在该区域内,粘性力与速度梯度成线性关系,因此在于核心层为高雷诺数湍流流动的情况下,过渡层与粘性子层内的速度分布可以通过经验公式直接计算得到,而无需划分网格,换句话说,在这种情况下,可以将计算节点的第一层网格节点放置在湍流核心区内,而过渡层与粘性子层中则无需要任何网格。这部分区域中的物理量分布采用壁面函数(wall function)来计算完成。需要用到壁面函数的湍流模型包括:k-epsilon模型,雷诺应力模型。
另外一种低雷诺数湍流模型的情况则与之不同,其不采用壁面函数来求解粘性子层与过渡层中的流动物理量分布,而是采用NS方程离散求解,与核心区域求解方式一样,如K-W模型,SA模型等。
2、Y+的问题
1、为什么会存在第一层网格的问题
关于这个问题,实际上要从边界层说起。实验表明,边界层内根据流动状态的不同可以分为三层,自壁面向流动核心区分别为:粘性子层、过渡层和湍流核心层。边界层很薄,一般都是毫米~微米级,因此,若采用划分网格进而利用数值方法求解的话,势必会大大增加计算网格的数量,从而急剧增加计算工作量。又有实验发现,在粘性子层和过渡层内,主要是粘性力在起主导作用,惯性力的作用几乎可以忽略。在该区域内,粘性力与速度梯度成线性关系,因此在于核心层为高雷诺数湍流流动的情况下,过渡层与粘性子层内的速度分布可以通过经验公式直接计算得到,而无需划分网格,换句话说,在这种情况下,可以将计算节点的第一层网格节点放置在湍流核心区内,而过渡层与粘性子层中则无需要任何网格。这部分区域中的物理量分布采用壁面函数(wall function)来计算完成。需要用到壁面函数的湍流模型包括:k-epsilon模型,雷诺应力模型。
另外一种低雷诺数湍流模型的情况则与之不同,其不采用壁面函数来求解粘性子层与过渡层中的流动物理量分布,而是采用NS方程离散求解,与核心区域求解方式一样,如K-W模型,SA模型等。
2、Y+的问题
