CFD网格质量计算器 返回
流体仿真

CFD网格质量指标计算器

实时计算宽高比、偏度、正交性。与OpenFOAM、ANSYS Fluent、Star-CCM+的允许范围对比,即刻评估网格质量。

拉普拉斯平滑改善网格质量(实时)

暂停时动画将静止。点击播放可恢复自动动画。

单元数
平均偏度
最大宽高比
最小正交质量
低于阈值的单元占比
平滑迭代
理论·主要公式

该动画对一个被扭曲的2D网格反复施加拉普拉斯平滑,并按质量为每个单元着色(绿=良好,红=不良)。可以观察到平均偏度随迭代不断下降。

等角偏度:

$$S = \frac{\theta_{max}- \theta_{equi}}{180 - \theta_{equi}}$$

$\theta_{equi}$:理想等角(四边形=90°)。$S\in[0,1]$(0=理想)

宽高比: $AR = L_{max}/ L_{min}$

正交质量: $Q_{orth} = 1 - S$

拉普拉斯平滑: $\mathbf{x}_i \leftarrow \mathbf{x}_i + \alpha\,(\bar{\mathbf{x}}_{nb} - \mathbf{x}_i)$(移向相邻节点的平均位置)

CFD网格质量指标说明

🙋
网格质量为什么这么重要?不就是计算的"前期准备"吗?
🎓
简单来说,网格的"形状缺陷"直接决定了计算结果的"可信度"。比如在汽车空气动力学分析中,如果引擎盖上的网格明显扭曲,该区域的压力分布就会发生锯齿波动,完全无法准确再现真实流场。这个计算器上的雷达图表,如果5项指标同时变红,就说明"用这个网格计算基本没有意义"。
🙋
那如果我拖动上面的"偏度"滑块,网格会变成什么样子呢?
🎓
偏度高的网格(接近0.9)中,四边形会变成菱形,三角形会变得又细又长。实际工程中最常见的就是复杂几何形状被强行网格化后产生的"变形单元"。这类单元出现时,Navier-Stokes方程在离散化时的"梯度计算"精度会大幅下降,导致压力和速度计算出现振荡。通过调节计算器中的数值,你能看到雷达图表的突刺变长,直观地反映出"这个网格有问题"。
🙋
"宽高比"和"膨胀比"看起来差不多,都是尺寸差异,有什么区别吗?
🎓
区别很大!宽高比是"单个单元内"最长边和最短边的比值。比如为了表现边界层,工程师会刻意用宽高比=1000以上的细长六面体网格。膨胀比是指"相邻两个单元的体积变化程度"。膨胀比过大时,流场的物理量(温度、压力)在单元边界处会产生不连续,导致计算容易发散。用计算器分别调节这两个参数,看雷达图表的反应不同,你就明白区别了。

常见问题

宽高比是网格单元最长边的长度除以最短边的长度。数值越接近1表示越接近理想的正方形或立方体,数值越大说明单元越细长。OpenFOAM一般推荐小于10,ANSYS Fluent推荐小于5。
偏度高表示单元角度与理想形状偏差大,导致离散化误差增大。特别是梯度计算和扩散项精度会严重下降,容易引起解的发散或产生非物理的数值振荡。一般来说,偏度超过0.8就需要重新调整网格。
工具支持OpenFOAM、ANSYS Fluent、Star-CCM+三个主流求解器。自动将计算的宽高比、偏度、正交性等指标与各求解器推荐范围对比,并用"良好""注意""需改进"三个等级评判网格质量。
正交性差时,需要重新检查网格的贴体处理和层状网格设置。具体方法是:降低边界层网格的增长率、增加光滑迭代次数,或对六面体网格重新调整块分割,使得单元面尽可能与流向垂直。

现实应用案例

汽车与飞机外部气动分析:车体或翼面表面的边界层网格采用宽高比很大的细长单元,但车体棱角等复杂区域容易产生高偏度单元。如果不仔细检查这些区域,阻力系数(Cd值)的预测精度会大幅下降。

电子设备热流体分析:散热片翅片间、芯片周围等几何形状复杂,网格质量容易恶化。膨胀比急速变化的网格会导致热传递系数计算不稳定,产生过度高估或低估的风险。

伴有化学反应的燃烧分析:在火焰面附近,浓度和温度梯度非常陡峭。该区域偏度大、正交性差的网格(偏度接近0.9)会严重扭曲离散化,错误计算反应速率和火焰传播速度。

与CAE软件联动:OpenFOAM的checkMesh、ANSYS Fluent的网格质量报告、Star-CCM+的网格诊断工具都计算类似指标。优质网格不仅能提高求解器收敛性,还能在更短时间内获得可信的结果。

常见误区与注意点

首先,"雷达图全绿就完美"的误解。实际上,根据分析类型和区域特点,有时需要"战略性地恶化"某些指标。例如,飞机翼型周围的边界层分析会故意使用宽高比超过1000的细长网格。这是为了准确捕捉物理现象的"有意设计",而非真正的缺陷。计算器给出的许可范围只是一般参考,最终判断还是要看"要计算什么"。

其次,"工具检查过就万事大吉"的自满。这个工具只看网格的"几何质量",不足以评估完整的仿真适用性。比如在超音速流中会产生激波,即使网格几何质量优秀,如果流向分辨率不够细(这是工具无法检测的),激波位置就会偏移。几何质量和物理分辨率都要兼顾。

最后,忽视局部坏单元的问题。整体平均指标很好,但在流场关键部位(比如汽车后视镜后的涡流区)隐藏着一个"偏度0.9"的单元,就足以从那里放大误差,破坏整个流场预测。计算器给出指标后,一定要将数值差的单元可视化出来,确认它们"是否处在关键位置"。

使用指南

  1. 从网格统计信息中输入宽高比(AR值)。滑块范围为1~100,边界层中更大的值需另行核对
  2. 设置偏度值(0~1)。越接近0越好,OpenFOAM推荐≤0.85,Fluent推荐≤0.90
  3. 指定非正交角θ(0~89°)。cos θ是输出的质量值,角度越小正交性越好
  4. 输入平滑度比和膨胀比ER(增长率),用100分制的综合评分判定网格可用性

具体计算示例

良好的翼型周围网格示例:宽高比AR=20、偏度=0.35、非正交角θ=25°、平滑度比=1.2、膨胀比ER=1.2。按页面权重计算,各子评分为AR 70、偏度88、正交性85、平滑度100、ER 100,综合评分85分,判定为“优良”。

工程实务注意事项