TMS BOOK · ACADEMY 课题
风窗气流组织与除霜风道 CFD
除霜风道的活儿,是把一层热风毯贴着玻璃铺满 A 区。风吹得响不等于霜化得快——决定快慢的是近壁换热系数。本课从贴壁射流与射流冲击讲到共轭瞬态 CFD,教你在数模阶段就把 A 柱根部那块死角找出来并改掉。
C5-02 · 乘员舱空调与 HVAC / 除霜除雾与视野安全
待认领listed
→
众筹中recruiting
→
已开课open · Course 售卖
一课题可并存多讲师版本
课程大纲
学完能做什么
- 说清贴壁射流覆盖与射流冲击换热在除霜上的分工与各自判据
- 用开口面积、导流叶片角度、风道走向三个手柄把流量分到 A / A' 区
- 建出含风窗共轭壁面的瞬态除霜 CFD 模型并读懂近壁换热系数分布
- 用 DOE 把 CFD 结论收敛成可制造的风道改进方案
- 判断什么交给 CFD、什么必须回环境舱实测
内容大纲
- 除霜气流组织的物理
- 射流冲击:格栅正对区换热强,但覆盖范围窄,只够啃一小块
- 贴壁射流与 Coanda 效应:靠玻璃壁面把射流「粘住」向上铺开,A 区覆盖主要靠它
- 脱流的两个诱因:玻璃倾角突变(A 柱根部)与逆压梯度;一脱流,换热系数掉一个量级
- 对撞与回流:左右格栅射流在中央对撞会抬起边界层,看似风大实则贴不住
- 除霜格栅与风道结构设计
- 格栅位置与出流角:太靠前打到玻璃中段、太靠后贴不住,出流角需与玻璃倾角配合
- 开口面积分布与局部阻力:流量分配由 A_i 与 ζ_i 共同决定,只改面积不改阻力常常失效
- 导流叶片:把流量往两侧 A' 区搬,是最便宜的整改手柄
- 风道走向与热损:仪表台下的长风道既掉温又掉压,还容易内漏到吹脚
- 可制造性约束:脱模方向、拼合面、公差对实际出流角的侵蚀
- 除霜风量的系统侧约束与模式门耦合
- 除霜模式下的风量占比与鼓风机工作点上移,需与 C2-04 蜗壳能力对齐
- 模式门半开/密封不良导致的吹脚泄漏,把除霜风白白漏掉
- 混合风门位置与出风温度、风量的耦合:温度上去了风量往往下来
- 整车级约束:除霜风量提高带来的风噪与 NVH 代价(见 C2-11)
- 除霜 CFD 建模
- 网格:格栅出口与近壁区加密,y+ 必须匹配所选壁面处理,否则 h 系统性失真
- 湍流模型:射流附壁与分离对模型敏感,realizable k-ε / SST k-ω 需做对比而非默认
- 风窗必须建成共轭壁面(CHT):玻璃导热 + 外侧低温迎风对流缺一不可。绝热壁与定壁温是方向相反的两个理想化——绝热壁零热流、无外散失,把内表面温度算成上界(偏高、偏乐观);定壁温则把壁温当输入而不是输出,壁面成了无限热汇/热源,偏差方向随所取壁温值(T_amb 还是 0℃)与评价判据而翻,本身不可预判。只有 CHT 让内表面温度是解出来的
- 瞬态而非稳态:融霜本质是有相变的瞬态过程,稳态云图会给出「全绿」的错觉
- 边界条件的真值来源:格栅出口温度/风量应取自 1D 风道模型或实测,别直接用 HVAC 出风口值;除霜支路的进口温度/风量也可由 J4-04「HVAC 空调箱总成级 CFD」的全箱+风道联算给出——热损与内漏正发生在箱内到仪表台下风道那一段,那一段无人算,进口边界就先错了
- 从 CFD 结果到设计闭环与试验对齐
- 评价量的排序:近壁换热系数与内表面温度 > 融霜面积-时间曲线 > 速度云图
- 找末位死角:按分区统计 h 的最低分位区间,定位成因(脱流 / 打偏 / 流量不足)
- DOE:以叶片角、开口面积比、格栅位置为因子,用融霜面积-时间曲线做响应
- 与试验相关性:红外热像看内表面温度场、PIV 看射流轨迹,先对齐再信 CFD(见 J7-01)
关键公式
Nu = C·Re^m·Pr^n
贴壁/冲击射流换热经验关联式的通用形式;除霜局部 h 估算的骨架,系数需按具体构型的文献或实测标定
q(x,y) = h(x,y)·[T_jet(x,y) − T_g(x,y)]
内表面局部热流——除霜快慢的真正决定量。速度大而 h 小(射流脱离玻璃)照样不化霜
U = 1/(1/h_i + δ_g/λ_g + 1/h_o)
风窗共轭壁面的串联热阻;CFD 里必须建成 CHT。绝热壁相当于令 1/h_o→∞(外散失归零),内表面温度算成上界、结论偏乐观;定壁温相当于把 T_g 直接指定掉,这条串联链压根不参与求解,偏差方向随取值与判据而变。两者都不能当「保守估计」用
ṁ_i/ṁ_total = f(A_i, ζ_i)
各格栅流量分配由开口面积与局部阻力系数共同决定;只放大面积而不降阻力,流量未必跟着走
y⁺ = ρ·u_τ·Δy/μ
近壁网格判据;y+ 落在标准壁面函数与低雷诺模型的适用区之间(缓冲层 5~30)时换热系数会系统性偏差。第三条路是 all-y+/增强壁面处理,但缓冲层仍是精度最低段
Q_melt = m_ice·(c_p,ice·ΔT + L_f),L_f ≈ 334 kJ/kg
瞬态融霜的能量口径,与 C5-01 一致;CFD 后处理把它折成融霜面积-时间曲线
关键概念
贴壁射流Coanda 效应射流冲击除霜格栅导流叶片出流角局部阻力系数 ζ共轭传热 CHTy+壁面函数回流死角融霜面积-时间曲线
推荐工具与标准
STAR-CCM+ / Fluent(瞬态共轭除霜 CFD) GT-SUITE / Flowmaster(1D 风道压降与流量分配,给 CFD 提供边界) HEEDS / Optimus / Isight(DOE 与形状优化) 红外热像仪(内表面温度场,CFD 相关性) PIV / 烟流显示(射流轨迹与脱流位置验证)
GB 11555-2025《汽车风窗玻璃除霜和除雾系统技术规范》(2026-01-01 实施,全部代替 2009 版)——CFD 的判据来源 FMVSS 103(风窗除霜与除雾系统) GB 11562-2014《汽车驾驶员前方视野要求及测量方法》(视区划分依据;**GB 11562-2025 已于 2025-12-31 发布、2027-01-01 实施**,全部代替 2014 版,适用范围扩至 M1+N1——现在按 2014 版写仍正确,交付前须按实施日复核) ASME V&V 20(计算流体力学与传热的验证与确认)
工程案例
某纯电轿车除霜时 A 区靠 A 柱两侧的角落迟迟不化,中间早已通透。CFD 一看:中央格栅出流强、两侧射流在玻璃倾角突变处提前脱离,A 柱根部形成低速回流死角——速度云图还挺好看,近壁换热系数却掉了一个量级。改导流叶片角度并重排开口面积分布,把流量往两侧搬,融霜面积-时间曲线整体左移。
动手做
交付物 · 拿一版风窗 + 除霜风道数模:① 建含玻璃共轭壁面的瞬态除霜 CFD,跑到法规时间节点;② 提取 A / A' / B 区内表面换热系数与温度分布,找出末位分位的死角并归因(脱流 / 打偏 / 流量不足);③ 出 2 版风道改进并用 DOE 对比,交付改进前后的融霜面积-时间曲线与一页可制造性说明。
常见误区
- 用稳态 CFD + 定壁温评价本质瞬态的融霜过程:稳态抹掉了融霜的时间维,定壁温又把内表面温度从待求量降成输入,云图全绿而实车 A 柱根部照样死角——这里的毛病不是「一定偏乐观」,而是这套设置根本给不出可判读的壁温历程
- 把绝热壁与定壁温当成同一类错误:绝热壁是无热损、T_g 算高、除霜时间偏短(偏乐观);定壁温是把 T_g 钉死成输入,方向取决于所取壁温与评价判据——取 T_amb 而按「内表面升到 0℃」判读则永远不达标(偏保守),换成按 q=h·(T_jet−T_g) 积分到 Q_melt 后处理,ΔT 与 q 双双高估、融霜时间反而算短(偏乐观)。拿错方向去解释自家老模型与实测的偏差,会把整改方向带反;只有共轭两头都对
- 只看速度云图不看近壁换热系数:射流打偏或离玻璃太远,风速再高也不换热
- 近壁网格 y+ 不匹配所选壁面处理,h 系统性偏差,做再多 DOE 也是在错的响应面上找最优
- 用 HVAC 出风口温度当格栅出口温度,忽略仪表台下风道的热损与内漏,进口边界就先错了
相关课题