🙋
CFD中经常听到"y+",具体是什么意思?模拟器中的"初始网格高度Δy"与它如何关联?
🎓
简单来说,它是决定壁面附近网格细度的无量纲数。y+太小计算成本太高,太大精度下降。这个工具根据雷诺数、粘性系数、流速等参数,实时计算推荐的初始网格高度Δy。比如在汽车外部空气动力学中调高流速,你会看到所需的Δy变小。
🙋
下面图表中看到的"粘性底层"和"对数律区域"有什么区别?为什么k-ε和k-ω SST推荐的y+不同?
🎓
贴着壁面速度线性增长的区域叫粘性底层(y+<5),外侧速度对数增长的叫对数律区域(y+>30)。k-ε模型用壁面函数,所以第一个网格中心点要放在对数律区域(y+≈30〜300)。而k-ω SST是"低Re数模型",直接解粘性底层,所以第一个网格要设y+≈1才能保证精度。切换模型看看,你会清楚地理解y+推荐值和图表的关系。
🙋
在工程实践中怎样使用这些信息?输入"湍流强度"和"长度尺度"是做什么的?
🎓
实际工作中,先用这个工具算出Δy,再根据它设置壁面网格第一层厚度。比如航空翼解析,选SST模型,设y+≈1得到Δy。"湍流强度"和"长度尺度"是用来估算k、ε这些湍流变量初值的。如果有风洞实验数据就用实验值,没有的话就从工具默认值开始,然后对照模拟结果调整。这是标准工业流程。
工具输出的y+是目标值。第1层高度Δy由预期的雷诺数和摩擦速度反算:Δy = y+ × ν / u*。以y+=1为目标时,按此公式得Δy,在建网时作为初值。
k-ε适合高Re数、无旋流(如管道流),需要标准壁面函数(y+≈30〜300)。k-ω SST特别擅长分离、逆压力梯度(如翼周)流动,可在低y+(≈1)也可在高y+(≈30)使用,通用性强。
计算值是根据设定湍流强度(如5%)和特征长度推算的估值。k代表湍流运动能,ε和ω代表散逸率。必须与实验数据或经验对比验证量纲是否合理,再才能输入软件。
LES直接解大涡,适于非定常现象(渡涡释放、噪声预测),但需极细网格(y+≈1)和小时间步。RANS(k-ε、k-ω SST等)主要预测平均场,定常分析和设计初期参数扫描更实用。
汽车空气动力学分析:车身周围流场和阻力系数(Cd)预测广泛应用CFD。通常采用k-ω SST模型配y+≈1网格,精确捕捉后视镜、支柱周围的复杂分离流。
飞机和翼设计:升力、阻力评估,失速特性预测是主要目标。需处理逆压力梯度下的边界层发展和分离,SST模型最合适,翼面网格设y+≈1。
建筑周围风环境分析:用于城市规划和建筑风环评。形状复杂,往往用相对粗网格的k-ε模型配标准壁面函数(y+≈30-300),在成本与精度间取平衡。
涡轮机械(泵、风机)内流:CFD对效率预测不可或缺。受强曲率和旋转影响,推荐SST模型,壁面网格非常细密(y+≈1)。
首先,"只要y+设对就万事大吉"这种误解要警惕。y+确实是第一步,但仅凭它不够。即使y+≈1,如果垂直方向网格层数太少(如只有5层)就无法解析整个边界层速度梯度,分离流预测会出偏差。通常需要10〜15层以上网格来覆盖边界层。
其次,"推荐y+值就是绝对准则"的思想要改正。工具给出"k-ε模型y+≈30"只是起点。如果实际流场很复杂(强压力梯度、分离和再贴附),可能需要用更小y+值(如10〜20)。必须检查壁面剪应力、压力分布结果,必要时重新生成网格,这种灵活度是专业人士的标志。
最后,随意设置入口"湍流强度"和"长度尺度"要避免。这些直接影响计算初值和流动发展。若模拟风洞实验,务必用实验值。数据缺失时,内流路(配管)一般用几%强度,外流(汽车周围)用0.5〜1%左右开始,然后必须做灵敏度分析确认影响,这才是职业做法。