模型不大,网格质量检查全绿,残差曲线稳稳往下走,收敛得很漂亮。云图出来,速度场、压力场看着都顺。可把关键的结果拿去跟实测一比,差得不多不少,就那么百分之十几,怎么调都对不上。
你开始怀疑边界条件给错了,改了一遍,重算,没好。又怀疑湍流模型选得不对,换了一个,重算,还是没好。松弛因子调小,网格加密一轮,机时花出去几宿,那个差值就在那儿待着,不增不减。
我在这个问题上卡过挺久。后来才弄明白,问题不在我改过的任何一个地方,在我从来没动过的地方。贴着壁面的那第一层网格。
流体贴着壁面流的时候,紧挨着壁面那一薄层的流动,跟管道中间的主流完全是两回事。
主流里速度变化是均匀的,你把网格放大一点,算出来的速度分布依然八九不离十。可贴着壁面的那一层,速度从零快速爬升到主流速度,这个爬升的过程高度非线性,中间还夹着粘性底层、过渡层、对数律层这几个脾气不同的区段。你在这儿放一层多高的网格,等于决定了你去分辨这一段里的哪一部分。
网格第一层太厚,这一薄层就被你一个网格跨过去了。软件只能在网格中心取一个代表值,把从零到主流那段变化糊成一格。算出来的壁面剪应力、壁面换热系数,自然就对不上。
y+ 就是干这个用的。它是无量纲的壁面距离,衡量的是你第一层网格中心,落在离壁面多远的位置,用壁面摩擦速度和运动粘度做尺度归一化。y+ 等于 1,说明你的网格已经细到能钻进粘性底层里。y+ 等于 100,说明你的第一层网格还悬在对数律层里,底层压根没碰到。
壁面第一层的高度,跟你算得准不准,是直接挂钩的。
这里最容易出的错,是把 y+ 当成一个孤立指标去追。y+ 本身没有绝对的好坏,它得跟湍流模型配套。
工程上常用的三类处理方式,对 y+ 的要求是分开的。
第一类是壁面函数法。这条路不打算解析粘性底层,而是在壁面和第一层网格之间架一个半经验公式,把底层的影响补上去。用它的时候,第一层网格中心得落在对数律层,y+ 一般在 30 到 300 之间,别往细里加。你要是好心把网格搞到 y+ 等于 5,反而掉进粘性底层里,壁面函数那个公式的前提就不成立了,结果比粗网格还差。
第二类是增强壁面处理,也叫两层模型。它把整个边界层区域都用公式覆盖掉,包括粘性底层,所以要求第一层网格落得足够细,y+ 最好在 1 上下。这条路适合流动分离、壁面换热这类底层影响大的问题,代价是网格量大。
第三类是直接解析到壁面。低雷诺数模型那类走的就是这条路,y+ 必须小于 1,边界层里要铺十几层以上的网格。算得准,也最费机时,通常用在换热、气动噪声这些对壁面极敏感的场合。
三种方式没有谁更高级。做了换热分析,用壁面函数法,结果就可能偏。做了简单的内流压降,硬要解析到壁面,机时全花在了没必要的地方。选哪个,取决于你这个问题里,壁面到底有多重要。
这一条是很多人漏掉的。
y+ 不是一个画完网格就定下来的数。它是在求解之后才算出来的,取决于壁面摩擦速度,而壁面摩擦速度取决于流场本身。也就是说,网格画的时候你觉得 y+ 是 30,算完之后一查,可能是 8,也可能是 120。
更麻烦的是换模型的时候。同一套网格,你从 k-epsilon 换到 k-omega SST,或者打开增强壁面处理,y+ 的适用区间就变了。上一版落在 30 到 300 里,看着挺合适。新版要求 y+ 接近 1,同一套网格瞬间就不合格了。
所以每次动了湍流模型或者壁面处理方式,回头复查一次 y+ 分布,这一步不能省。软件里都能直接画出壁面的 y+ 云图,看一眼数值范围,跟你选的模型对不对得上,心里就有数了。
说点能直接上手的。
第一步,先估一个壁面摩擦速度。手头有经验公式,按来流速度、雷诺数估个量级出来,不用很准,能定个量级就行。
第二步,按你选的湍流模型,定一个目标 y+。壁面函数那类取 30 到 300,增强壁面处理取 1 左右,解析到壁面取小于 1。
第三步,反推第一层网格高度。有了目标 y+ 和估出来的摩擦速度,加上流体的运动粘度,就能算出第一层网格该多厚。
第四步,往回铺边界层网格。光定第一层不够,还得在边界层厚度范围内铺够层数,让网格从第一层平滑地过渡到主流网格尺寸。层数不够,第一层再薄,网格尺寸梯度也会跳。
第五步,算完回查。前面说的 y+ 云图,看实际值有没有落在你定的区间里。落在里面,这一步就算过了。跑出去了,回头调第一层高度,重算。
拿一段内部流道的换热算例说说。
第一版网格,全局尺寸设了 3 毫米,壁面附近没单独处理。算完 y+ 云图一看,壁面大部分区域在 200 到 400 之间。那时候用的是增强壁面处理,要求 y+ 接近 1,差了三个数量级。算出来的换热系数比实测低了两成多。
第二版,把边界层网格打开,第一层设 0.05 毫米,铺了 10 层,增长因子 1.2。重算,y+ 落到 30 上下。这个值对壁面函数法合适,对增强壁面处理还是偏粗。换热系数跟实测的差距缩到了一成。
第三版,第一层压到 0.005 毫米,铺 15 层。重算,y+ 落在 1 附近,换热系数跟实测相差在百分之五以内。相邻两版的第一层再减半,结果几乎不动了,收手。
整个过程三版网格,重点全在壁面那几层,主流网格尺寸一次没动。
回头看这件事,卡住我很久的不是软件不会用,是壁面附近这一层的物理我一直没想明白。网格画得漂亮、残差收敛得漂亮,跟结果对不对,是三件事。
壁面这一段该怎么处理,第一层该多高、铺几层、跟哪个湍流模型配,这些都有现成的判断依据。要是有人从一开始就把这层物理讲透,再带着在真实算例里走一遍,能少花很多几宿的机时。
仿真高研院的 Fluent 流体工程师双证研修班,讲的就是把这类工程判断立起来,从理论到软件操作再到案例交付,一步一步走。边界层和壁面处理,是里面绕不开的一环。
如果你也在自己摸索,有些判断拿不准,可以先领一个《仿真工程师从入门到就业知识库》,看看自己卡在哪个阶段,心里有个数。

扫码即可领取。对着看看自己在哪个阶段,下一步该补哪一块。
了解一下
Fluent 流体工程师双证班
如果你厌倦了无序自学的低效,想要走一条更高效的学习路径,推荐你了解仿真高研院的《fluent流体双证班》。
课程分六个阶段,从 CFD 理论和流体力学基础讲起,再到网格技术、求解器实战、湍流模型深度应用、多相流与动网格进阶,最后落到工程仿真方案设计。每一步都踩在前一步的底子上。
主讲老师长期泡在企业仿真一线,讲的是踩过坑之后的正确做法。案例取材真实企业项目,练的是能写进简历的活儿,不是课本例题。
每个阶段的作业有人逐条批改,错在哪、该怎么改,都会讲清楚。学习支持是三合一的,录播课永久回看,配上直播面授和社群答疑,漏了课不慌,卡住了有人接。
课程结束后参加考试,合格可以拿到工业和信息化部教育与考试中心颁发的证书。
对于培训内容,你还有什么想要知道的,可以扫码咨询。