【CFD学习】yplus，y+，边界层理论，壁面第一层网格

最新推荐文章于 2025-03-18 22:23:25 发布

BLADE_XuYiBing

最新推荐文章于 2025-03-18 22:23:25 发布

阅读量9.4w

点赞数 132

分类专栏： CFD学习笔记文章标签：流体力学 cfd

本文链接：https://blog.csdn.net/dabing1212/article/details/107637215

版权

CFD学习笔记专栏收录该内容

2 篇文章

订阅专栏

谈谈本人对计算流体力学中y+的理解，欢迎批评指正！

思考几个问题：

1. Q:什么时候需要考虑y+？

A:粘性计算（求解N-S方程）才涉及y+，无粘计算（求解Euler方程）不考虑边界层，不涉及y+的讨论

2.Q:为什么需要考虑y+?（理论基础）

A：粘性计算时，壁面附近是流场中可以预先知道的高梯度区域，即边界层区域。即使边界层不是流场中速度梯度最大的区域，但至少是已知的大速度梯度区域。因此，应当预先考虑：需要用多少个网格，才能分辨出边界层呢？

要回答这个问题，不仅需要预估边界层厚度，还需要大致了解边界层内的速度分布规律。

冯·卡门在普朗特混合长度理论的基础上，用量纲分析的方法推导出了对数律（log law），后续的理论与实验研究，界定了对数律的成立范围，并进一步将边界层细分为4个区域：粘性（线性）底层、缓冲区、对数律区，外层，如下图：

由上图可知，在对数坐标系下，两个区域具有解析的函数关系式。分别是：

粘性（线性）底层：惯性力相对粘性力较小，粘性力与速度梯度呈线性关系。

对数率层：无量纲速度 $U^{+}$ 与无量纲壁面距离 $y^{+}$ 呈对数关系，即卡门对数律。

可以利用这两个特殊区域中的关系，建立所谓Wall function，替代解析求解边界层内流动的过程。这样做的优势是，允许我们用更大的壁面第一层网格，来捕捉（毫米，微米量级的）边界层。极大的减少壁面法向的网格数量，避免了壁面处网格长宽比（aspect ratio）过大造成的数值误差和数值不稳定，化解了网格数和网格质量的矛盾。

$y^{+}$ 作为区分边界层内部结构的特征参数，在网格划分时需要重点关注。

3. Q:求解中如何利用wall function?

A:首先需要明确的是，流动并不一定具有如上的边界层结构。

对于低雷诺数流动（Low-Re）：

边界层较厚，主流和边界层的掺混，会使得对数律区的范围很小，可以认为没有对数律区的存在，不存在Lindgren图中所示的内部结构。此外，对于分离流动，一样不存在对数律区，如，存在分离的翼型、叶轮机内的流动等。这时候，对数律区的logarithmic function不能使用，因此，必须将壁面第一层网格 $y^{wall}$ ，划分到粘性底层，对应的 $y_{1}^{+}< 10$ ，并且能接近于1（这里的 $y_{1}^{+}$ 是壁面第一层网格的无量纲壁面距离）。此方案，称为低雷诺方案（Low-Re mothod）。

对于高雷诺数流动（High-Re）：

边界层较薄，存在Lindgren图中所示的内部结构，也存在对数律区。这时候，既可以把壁面第一层网格划分到粘性底层，使得 $y_{1}^{+}< 10$ ，即，采用低雷诺数方案；也可将其划分到对数律区域，使得 $20< y_{1}^{+}< 50$ ，利用对数律wall function接近，即，采用高雷诺数方案（High-Re mothod）。

需要说明的是：这里的所谓高低雷诺数，并非以过流装置的尺寸和平均流速为特征参数计算得到的装置雷诺数（device Reynolds number），而是以壁面局部区域，边界层发展长度和壁面附近自由流为特征参数，计算得到的局部的湍流雷诺数（turbulence Reynolds number）。

因此，对于高雷诺数下的设备，也有可能必须采用低雷诺数方案（Low-Re mothod），如存在分离的叶轮机内的流动。并且，即使是局部雷诺数依然很大的请况，一样可以采用Low-Re mothod解析到粘性底层，只是耗费了更多网格而已。

总结一下：

当 $y_{1}^{+}< 10$ 时，我们认为PDE的求解已经数值积分到壁面，边界层被完整求解。

当 $20< y_{1}^{+}< 50$ 时（一般在此区间，取决于具体流动类型的buffer layer），壁面第一层网格落在log law layer，此时我们在求解过程中引入壁面函数wall function来近似求解第一层网格上的速度，此时，PDE的求解没有积分到壁面。

当 $10< y_{1}^{+}< 20$ 或 $y_{1}^{+}>50$ 时，不被允许用这样的网格求解含边界层的流动。

注：工程计算中认为粘性底层的边界是 $y_{1}^{+}< 10$ ，理论研究计算中认为 $y_{1}^{+}< 5$

4. Q:怎样确定y+的?（壁面第一层网格和y+的关系）

A: $y_{1}^{+}$ 是考虑了第一层网格厚度 $y^{wall}$ 和具体流动特征参数的无量纲壁面距离，表征了第一层网格在Lindgren图边界层分区结构中的位置。

正问题求 $y_{1}^{+}$ ：

根据已知流场和壁面第一层网格的厚度 $y^{wall}$ ，可以直接求解 $y_{1}^{+}$ ，公式如下：

$y_{1}^{+}=\frac{\rho u_{\tau }y_{wall}}{ \mu }$

其中 $u_{\tau }$ 是壁面摩擦速度（friction velocity）：

$u_{\tau }=\sqrt{\frac{\tau_{wall}}{\rho }} =\sqrt{\frac{1}{2} (V_{ref})^{2}C_{f}}$

反问题预估 $y^{wall}$ ：

由于流场也未知，因此，预估壁面第一层网格比较难。

一种方法是，通过Blasius 方程的截断级数解，给出 $y^{wall}$

$y^{wall}=6(\frac{V_{ref}}{\upsilon })^{-\frac{7}{8}}(\frac{L_{ref}}{2 })^{\frac{1}{8}}y_{1}^{+}...$

推荐网站，yplus计算工具：

http://www.pointwise.com/yplus/

https://www.cfd-online.com/Tools/yplus.php

一些网格绘制的建议：

湍流模型与y+的匹配

低雷诺数模型（如k-w模型，SA模型等），需要满足 $y_{1}^{+}< 10$ ， $y_{1}^{+}\sim 1$ 最佳

高雷诺数模型（如k-Epsilon模型、雷诺应力模型等），需要满足 $20< y_{1}^{+}< 50$ ， $y_{1}^{+}\sim 30$ 最佳

自带壁面函数型模型（如SST，RNG模型），需要满足 $20< y_{1}^{+}< 50$ ， $y_{1}^{+}\sim 30$ 最佳

大分离流动相关模型（LES，DES，SAS） $y_{1}^{+}\sim 1$ ，趋向正方型网格最佳

另外，所有的模型均可以强行添加wall function，添加wall function后原则上，需要满足 $20< y_{1}^{+}< 50$ ， $y_{1}^{+}\sim 30$ 最佳

在商用软件中y+重要么

需要说明的是，很多商用软件，可以实现自动壁面处理，可以适应于任意精细的网格，如CFX中的scalable wall function，自动壁面处理可以实现（Low-Re mothod）和（High-Re mothod）的光滑切换，因此对 $y_{1}^{+}$ 不必严苛要求。

尽管如此，仍然建议，采用Low-Re mothod时将第一层网格画到 $y_{1}^{+}\sim 1$ ，采用High-Re mothod时画到 $y_{1}^{+}\sim 30$ 。因为，本人在做网格无关性检验时发现，网格数保持不变，仅改变壁面第一层网格从 $y_{1}^{+}\approx10$ （很少的局部区域大于10）到 $y_{1}^{+}\sim 5$ ，流动参数的定量影响仍相当可观（离心叶轮，效率变化约0.5个百分点）。