3. 网格生成策略与热管理模型前处理
在新能源汽车热管理仿真中,几何与网格前处理往往占据整个项目 60% 以上的工作量。电池包液冷板内部流道细长、风冷风道截面突变、电芯与导热垫之间存在大量薄壁接触,这些特征决定了网格不能“一键生成”了事。本节围绕几何清理、流体域抽取、边界命名、多域网格划分、边界层设置与网格质量评价六个环节展开,并以电池包液冷板与风冷风道为实战对象,给出可复现的操作路径与判断依据。
一、几何清理:先让模型“可网格化”
CAD 模型面向制造,而 CFD 模型面向求解,两者诉求不同。导入 Fluent Meshing 或 SpaceClaim 后,首先要处理三类问题:微小特征(倒角、螺纹、工艺孔、Logo 凹槽)、拓扑缺陷(重复面、缝隙、自由边、干涉)和薄壁结构(液冷板上下盖板、电芯外壳)。
清理原则是“该删的删、该补的补、该简化的简化”。例如液冷板进出口的螺纹接头,若研究目标是流道流量分配而非接头局部阻力,可直接用圆柱面替代螺纹,直径取公称通径。对于电芯与导热垫之间 0.2 mm 的装配间隙,若小于最小网格尺寸,应合并为共享拓扑或按接触热阻处理,而不是强行加密网格。清理完成后用“Geometry → Fault Detection”检查自由边与干涉,目标是将自由边数量降为 0,干涉体积小于总体积的 0.1%。
二、流体域抽取与边界命名
热管理仿真通常同时存在固体域(电芯、冷板金属、导热垫)和流体域(冷却液、空气)。流体域抽取的常用方法是 SpaceClaim 的“Volume Extract”或 Fluent Meshing 的“Fluid Volume Extraction”,通过封闭面组生成内部流体体积。抽取后必须确认流体域是单一连通体,若出现多个连通域,说明存在未封闭的缝隙或多余隔板。
边界命名是后续设置边界条件的基础,建议在几何阶段就完成命名,而不是在求解器里逐个点选。命名规范可采用“域类型_物理位置_编号”的形式,例如:
coolant_inlet_01 冷却液进口
coolant_outlet_01 冷却液出口
coldplate_wall 冷板与流体交界面
air_inlet_01 风道进口
air_outlet_01 风道出口
cell_surface 电芯表面
thermal_pad_bottom 导热垫底面
命名时要注意:进口与出口不能混用同一名称;固体域与流体域的交界面应单独命名,便于后续设置为 coupled wall 或 interface。若使用 Fluent Meshing 的 Watertight Geometry 工作流,边界命名会直接映射为 mesh 的 boundary zone,减少求解器中的重复劳动。
三、多域网格划分与边界层设置
热管理模型通常包含流体域与固体域,两者对网格的需求不同。流体域需要边界层来捕捉速度与温度梯度,固体域则更关注导热路径上的网格连续性。推荐采用 Poly-Hexcore 或 CutCell 方法:流体域用多面体网格配合棱柱边界层,固体域用六面体主导网格,两者在交界面处通过共节点或 interface 连接。
边界层设置的核心参数是首层高度、层数和增长率。首层高度由目标 y+ 反推:对于液冷板,冷却液流速约 0.5~1.5 m/s,水力直径 4~8 mm,若采用 SST k-ω 模型,目标 y+ 取 1~5,首层高度约 0.02~0.05 mm;对于风冷风道,空气流速 3~8 m/s,若采用 Realizable k-ε 配合标准壁面函数,目标 y+ 取 30~100,首层高度约 0.1~0.3 mm。层数一般取 5~10 层,增长率 1.2 左右,总边界层厚度不宜超过流道特征尺寸的 1/3,否则会挤压核心区网格。
多域网格划分时,建议先划分流体域并检查边界层质量,再划分固体域。若使用 Fluent Meshing 的“Share Topology”功能,可保证交界面网格连续,避免插值误差。对于液冷板与电芯之间的导热垫,若厚度仅 0.5 mm,可单独划分 2~3 层六面体网格,或将其简化为薄壁热阻边界。
四、网格质量评价与实战判断
网格质量直接决定求解稳定性与精度。Fluent 中常用的评价指标包括正交质量(Orthogonal Quality)、偏斜度(Skewness)、长宽比(Aspect Ratio)和边界层 y+。工程上可参考以下判断标准:
| 指标 | 可接受范围 | 目标值 | 不达标时的处理 |
|---|---|---|---|
| 正交质量 | ≥ 0.15 | > 0.3 | 局部加密、改用 Poly-Hexcore |
| 偏斜度 | < 0.85 | < 0.7 | 调整节点分布、光顺网格 |
| 长宽比 | 边界层内 < 50 | < 30 | 减小增长率或增加层数 |
| 边界层 y+ | 液冷 1~5,风冷 30~100 | 按湍流模型匹配 | 调整首层高度 |
实战中,液冷板蛇形流道转弯处容易出现高偏斜度网格,原因是转弯半径小、边界层挤压。处理方法是:在转弯处设置局部尺寸控制,将面网格尺寸从 1 mm 细化到 0.5 mm,同时将边界层增长率从 1.2 降到 1.1。风冷风道渐缩段则容易出现长宽比过大,可在渐缩段增加横向节点数,或改用扫掠网格。
以某 60 kWh 电池包液冷板为例(示例参数,非实测):流道宽 6 mm、高 3 mm,冷却液为 50% 乙二醇水溶液,进口流量 8 L/min,目标 y+ 取 3,计算得首层高度约 0.035 mm,设置 6 层边界层,增长率 1.2,总厚度约 0.28 mm,占流道高度约 9%,核心区网格仍保持 0.8 mm 的多面体尺寸。风冷风道截面 120 mm × 40 mm,空气流量 0.05 kg/s,目标 y+ 取 50,首层高度约 0.15 mm,设置 8 层,增长率 1.25。按此策略生成的网格,正交质量最低约 0.22,最大偏斜度约 0.78,满足求解要求。
网格生成不是一次性动作,而是“划分—检查—局部调整—再检查”的迭代过程。把边界命名和几何清理做扎实,后续的边界条件设置、材料赋值和求解监控都会顺畅很多。