摘要
今日学习了钝体绕流模拟的网格处理、四种湍流模型性能对比,以及通过分离泡尺寸和壁面摩擦系数验证模型精度的方法。
Abstract
Today, I learned about grid processing in bluff body flow simulation, performance comparisons of four turbulence models, and methods for validating model accuracy through separation bubble size and wall friction coefficient.
仿真流体Fluent:
介绍了近壁面网格尺寸的理论预估方法,特别是针对平板流动的第一层网格高度计算。它首先给出了湍流摩擦系数的经验公式,并基于该公式推导出了从壁面到第一层网格中心距离的计算公式,该公式与壁面雷诺数密切相关。通过一个具体的工程实例进行演示,在已知空气流速、密度、粘度和平板长度的条件下,计算得出若要达到目标无量纲距离y+为50,第一层网格的高度应设定为约1毫米。这为准确定义前处理中的边界层网格提供了直接的理论依据和实用计算手段。

讲解了CFD仿真中进出口边界湍流参数的设置方法,具体内容如下:
-
设置必要性:当湍流通过入口或出口进入流体域时,必须定义湍流参数,具体设定取决于选用的湍流模型。
-
四种设置方法:
-
直接输入k、ε、ω或雷诺应力分量(唯一支持分布输入的方法);
-
湍流强度和长度尺度(边界层流动取0.4699L,下游流动取开口尺寸);
-
湍流强度和粘性比(默认方法,适合外流和内流);
-
湍流强度和水力直径(主要适合内流)。
-
-
参数推荐值:
-
湍流强度通常为1%~5%,其中5%适用于大多数问题;
-
粘性比推荐使用1~100,默认为10;
-
内流场中,粘性比通常为10~100。
-
-
软件操作示例:在Ansys界面中,速度入口的湍流设置下拉菜单展示了四种可选方法,当前高亮的是"Intensity and Hydraulic Diameter"。
梳理了进出口湍流参数的设置逻辑、选项与推荐值,为仿真前处理提供了明确的操作指导。
。针对Realizable k-ε模型,指出在所有壁面位置保持wall Y+在30-300范围内较困难,因此推荐使用可扩展壁面函数(scalable wall function),该模型可无视Y+数值直接处理近壁面速度梯度,减轻前处理负担。根据网格疏密给出曾经的建议:网格较大时采用k-ε配合标准壁面函数(对应y+约30-300的对数律区域),网格较密时采用SST配合k-ε增强壁面函数(对应y+≈1的粘性底层)。图中还展示了软件界面中湍流模型的层级选择路径,以及新旧版本建议的对比提示。


该案例为钝体平板绕流模拟,主要内容如下:
-
案例设置:采用四种不同湍流模型进行模拟,划分8700个四边形网格,在回流再附着区和前缘附近局部加密,并采用非平衡边界层处理;来流雷诺数Re_D=50,000,研究对象为矩形钝体平板绕流,关注其尾部分离泡与再附现象。
-
分离泡预测结果:标准k-ε模型(SKE)严重低估分离泡尺寸,而可实现的k-ε模型(RKE)能准确预测分离泡大小,显示其对逆压梯度和分离现象有更好的捕捉能力。
-
壁面摩擦系数验证:通过对比不同模型预测的壁面摩擦系数C_f分布与实验数据,发现标准k-ε预测结果与实验偏差较大(再附着点位置预测不准),而RNG k-ε、雷诺应力模型和可实现的k-ε与实验数据吻合较好;实验观测到的再附着点位于x/D=4.7处,表明改进型k-ε模型在近壁面和分离流模拟中具有更高精度。