2.2 网格划分实战：从物理场控制到用户控制

2.2 网格划分实战：从物理场控制到用户控制

本节摘要：网格是把连续空间翻译成离散代数的协议，翻译质量决定计算成败。本节先讲清四类网格单元的分工与单元阶次的含义，再给散热器实际划网：物理场控制网格快速出活，用户控制网格用尺寸节点与扫掠把单元数压下来，最后用统计数字验收。

网格选型第一课：单元形状决定命运

COMSOL 里能用的网格形状就四种，各自的地盘划分得很清楚。自由四面体是默认全能选手，任意复杂几何都能填满，代价是单元数多、精度密度一般；六面体（配合映射与扫掠生成）是精度之王，同样节点数下误差最小，但只肯住在规则的几何里；棱柱（三棱柱）专攻一件事——边界层，把薄层单元贴着壁面铺出去，流体与共轭传热的刚需；金字塔是过渡件，专职给六面体与四面体的交界处补缝。选型逻辑因此简单：几何规则用六面体系，几何不规则用四面体，算流动必加棱柱边界层，三者混用靠金字塔过渡。

单元"阶次"是另一个独立维度。传热与结构力学接口默认用二阶拉格朗日单元——每个单元边上还套着中间节点，两单元就能弯曲逼近二次的温度分布；层流接口默认一阶速度、一阶压力。同样的几何，二阶单元数少而自由度多，精度更高但内存更大。新手常见误区是把"加密网格"和"提高阶次"混为一谈：网格细不细说的是单元尺寸，阶次说的是单元内部的插值多项式，两条路都能提精度，成本结构不同。

图：四类网格单元与适用场景

给散热器划网：两条路线对比

**路线一，物理场控制网格。**右键"网格 1"→ 编辑物理场引导网格，单元尺寸选"常规"，全部构建。几十秒后完成，右键"统计"看家底：单元数约 5 万出头，最小单元质量 0.18，平均单元质量 0.66。能用，但翅片根部有不少"瘦长"单元——翅片只有 2 mm 厚，全局尺寸约 3.5 mm 的四面体为了塞进缝隙，被迫畸变。

**路线二，用户控制网格。**点"网格 1"设为用户控制，先添一个"自由四面体 1"，再右键它添加"尺寸"子节点：几何实体层选"域"，点选全部 7 片翅片，定制里勾选"最大单元大小"填 1.5[mm]。全部构建后统计变为：单元数约 4.6 万，最小单元质量升到 0.41，平均 0.78。翅片内的单元横竖均匀了，梯度大的翅根区域也顺带被加密——因为单元小的地方自然能看清梯度。

**再进一步，扫掠。**基板与翅片都是标准的"挤出体"，横截面规整，正适合扫掠：右键"网格 1"先加"映射 1"，选散热器底面（z=0 的面）作为铺面目标，构建得到规则四边形面网格；再加"扫掠 1"，源面自动取映射结果，目标选剩余的顶面与翅片顶，分布设"单元数"填 4 层，构建。统计对比路线一：单元数约 1.1 万，平均单元质量 0.93。同样的几何，六面体网格用五分之一的单元拿到了更高的质量——这就是选型的价值，不是玄学，是账面数字。

网格统计对比（散热器，同一几何）
路线  物理场控制-常规    用户控制+尺寸域定制   映射+扫掠
单元数     约 52 000           约 46 000          约 11 000
平均质量      0.66                0.78               0.93
最小质量      0.18                0.41               0.71
​

单元尺寸控制的三个旋钮

用户控制网格里最常用的是"尺寸"节点的定制项。最大单元大小管绝对尺寸，梯度大的区域按"特征长度内至少给五到六个单元"的经验取值——散热器翅片高 20 mm、温度沿高度近似线性，其实不用很细；而第 4 章流道的壁面边界层，就要靠"边界层属性"节点铺棱柱了。曲率因子默认 0.3，数值越小弯曲处单元越密，倒角圆孔多的模型可收到 0.2。狭窄区域分辨率默认 1，控制薄缝隙里至少铺几层单元，薄翅片、薄垫片类几何建议保持 1 以上，否则薄板方向只有一层单元，厚度方向的梯度直接失真。

⚠️ 网格划分报错"几何中存在比单元尺寸小得多的边"，先回 2.1 的几何检查找短边，而不是把网格调到极细去迁就脏几何。

本节要点回顾

• 四类单元四种地盘：四面体通用、六面体精度高但挑几何、棱柱管边界层、金字塔管过渡。
• 阶次与尺寸是两码事：传热默认二阶拉格朗日；加密是缩小尺寸，升阶是提高单元内部插值。
• 散热器实测账目：物理场控制 5.2 万单元、映射加扫掠 1.1 万单元且平均质量 0.93，规则几何优先扫掠。
• 三个旋钮：最大单元大小按特征长度给五到六个单元，曲率因子默认 0.3，狭窄区域分辨率保住薄壁层数。
• 网格统计（单元数、平均与最小单元质量）从本节起作为每次划网的验收指标，诊断细则在下一节。

补一个网格密度分配的实用原则：不是均匀加密，而是按物理梯度分配。温度梯度大的区域（热源附近、材料界面处、几何突变处）需要细网格，温度均匀的区域用粗网格就够了——COMSOL 的"物理场控制网格"会自动做这件事，但"用户控制网格"让你有更精确的控制。实际操作建议：先用物理场控制网格跑一遍，看温度分布图上的梯度集中区域；如果关键位置的等温线太稀疏（说明网格太粗），就在那个区域加尺寸控制（Size 子节点）加密。另一个容易忽略的细节：薄壁结构的厚度方向至少要有两层单元，否则厚度方向的温度梯度会被完全忽略——这对于散热器的翅片和 PCBA 板尤其重要。网格划分不是一次到位的——先跑粗网格看趋势、再加密验证收敛——比一开始就建极细网格然后等半天求解效率高得多。

再补一个网格尺寸的估算方法。COMSOL 的物理场控制网格默认分了从极粗到极细多个档位，但每档对应的单元尺寸并没有在界面上明示。一个实用的估算方法：看求解后的日志里"自由度数"——稳态传热的自由度约等于节点数，节点数与单元数的关系因单元阶次而异（二阶拉格朗日单元的自由度约等于节点数）。对于散热器这种三维多域问题，XS 档可能在 5000 自由度、S 档 20000、M 档 80000 的量级。自由度数直接决定求解时间和内存消耗：100 万自由度约需 8-16 GB 内存、求解时间数分钟——如果日志显示自由度超出你的内存预算，就要从物理场控制切换到用户控制网格并适当放粗。总之，自由度数是连接"网格设置"与"计算资源"的桥梁——每做一步都要看一眼日志里的这个数。

网格划分的最后一个实战技巧是用"统计"面板做全局体检：右键网格节点选统计，查看总单元数、最小单元质量、平均单元质量与单元类型分布。最小单元质量低于 0.2 的网格在求解时可能产生数值振荡或假解——定位方法是在图形窗口按质量着色显示网格，红色单元就是问题单元。如果红色单元集中在某个几何尖角处，用虚拟操作（删除面、忽略边）或局部尺寸控制修正即可。如果红色单元分散在多个位置，说明整体网格策略需要调整——比如从自由四面体切换到以六面体为主的扫掠网格。网格统计面板是你网格质量的第一道体检，每次修改网格后都应该看一眼，养成习惯后你会对"这个网格能不能用"产生直觉判断——看到总自由度数和平均质量就能预判求解的时间和精度，这是资深仿真工程师的核心竞争力之一。