C5-02 风窗气流组织与除霜风道 CFD
课程代码 C5-02 · 板块 C 乘员舱空调与 HVAC / 除霜除雾与视野安全 时长 约 4.0 小时(7 讲) 前置 C5-01 除霜除雾法规与视区划分;J2-01 CFD 基础、湍流模型与网格划分;A2-05 空气侧压降、风量与风阻 本讲义定位 讲师授课蓝本 / 学员自学讲义,是 C5-02 大纲的完整展开版
引言:风吹得响,不等于霜化得快
冬天清晨那几分钟,几乎每个人都经历过:除霜档拨到底,鼓风机的声音大到说话要提高嗓门,风窗中间那一块很快化开、能看见前车了,而靠 A 柱两侧的两个角落还蒙着白,迟迟不退。副驾会问一句"风都开这么大了怎么还不化"——这句话里藏着一个很自然、也很难自己察觉的判据错误:把"风量"当成了除霜的驱动量。
工程侧的版本更体面,但结构一模一样。方案评审上摆一张速度云图,玻璃前面一片饱满的高速区,会上没人提出异议;实车进环境舱,A 柱根部照旧是最后化开的那一块。于是下一轮动作往往是"再加点风量"。加风量要付鼓风机工作点上移、风噪与整车功率的代价,这些代价确定会发生,而收益可能全部落在早就化开的中央区域。方案没错,判据错了。
这门课要交付的,就是把判据换过来,并且把换过来之后该怎么算、怎么读、怎么改,一步一步走通。
钉子一:能动的只有换热系数的空间分布
除霜快慢由玻璃内表面的局部净热流决定:
其中 h(x,y) 是该点的内侧对流换热系数(W/(m²·K)),T_jet(x,y) 是打到该点的射流温度,T_g,i(x,y) 是该点的玻璃内表面温度。这个式子是逐点成立的:两个因子都是场量,不是标量。而气流组织这件事,能动的只有两项——h 的空间分布,以及 T_jet 沿风道、沿贴附行程的衰减。
请注意式子里没有风速。这不是笔误,是本课要敲的第一根钉子:本课全部手柄——格栅位置、出流角、开口面积 A_i、局部阻力 ζ_i、导流叶片、风道走向——最终都只通过一件事起作用,就是改变 h(x,y) 的分布(⚠ 这条"把手段收敛到同一把尺子上"的口径是本课归纳的,⛔ 不是大纲或任何标准里的既有条款;立它是为了让度数、平方毫米、布置走向这些量纲完全不同的手柄能放在一起比较)。
反过来说,⛔ 风速不是判据。射流一旦离开玻璃,主流可以继续以很高的速度掠过,而近壁速度梯度已经塌掉,近壁换热系数随之出现量级跌落(⚠ 这是本课引述的一句定性量级判断,只用来说明"速度不是判据"这一件事,⛔ 不给具体倍数,也⛔ 不得据它反算脱流之后的任何绝对值)。也就是说,"速度云图很漂亮"与"A 柱根部化不开",是可以同时成立的两件事——它们之间没有矛盾,只是被两个不同的场量分别描述了。
钉子二:CFD 交付的是空间分布与相对排序
第二根钉子是方法学侧的:在本课里,CFD 的交付物是"空间分布 + 相对排序",⛔ 不是合规结论、⛔ 不是一个分钟数。
一个除霜 CFD 模型有三处会让 h 场系统性地错,而云图、残差曲线和守恒量全都看不出来:① 进口边界取错——拿 HVAC 出风口的温度与风量当格栅出口值,中间那段仪表台下风道的热损与内漏无人计算;② 壁面理想化取错——风窗没有建成共轭壁面;③ 近壁网格与所选壁面处理不配对——第一层网格落进了缓冲层。三处里任何一处错,整个 h 场都是系统性偏的。
而 DOE 恰恰建在这个 h 场上。系统性偏差对每一个算例都一致地作用,反而让响应面拟合得更平滑、置信区间更窄——偏差越系统,报告看起来越可信。于是最优点落在了错的地方,而报告里的统计检验一切正常。这就是本课为什么要花三讲(第 3、4、5 讲)去讲那三处设置,而不是直接跳到"怎么优化风道"。
反钉子:⛔ "简化总是往安全侧偏",这句话在这里是假的
读者最容易读反的一条是这样的:"壁面边界条件本来就是个简化,简化总是往安全侧偏——所以我这套(绝热壁 / 定壁温)算出来的除霜时间是上界,实车只会更快。"
它之所以有说服力,是因为它从一句真话(理想化必然带来偏差)推出了一句假话(偏差的方向是保守的)。很多人手上已经有一版跑得动、也跟实测大致对得上的老模型,"简化=保守"给了它一个不必再查的理由。
两种理想化的真实方向是这样的:
- 绝热壁相当于令 1/h_o → ∞,把外侧低温迎风的散失整项归零。结果是内表面温度被算成上界(偏高)、越过融点更早、锋面净得热被高估 ⇒ 除霜时间偏短,是偏乐观的那一侧。方向可预判,但恰好与直觉相反。
- 定壁温相当于把玻璃温度从待求量降成输入,串联热阻链压根不参与求解,壁面成了无限热汇或热源。它的偏差方向随"取什么壁温 × 按什么判据读"而翻:同样取环境温度,按"内表面升到 0 ℃"判读则永远不达标(偏保守,而且是不可用的那种保守);改成按 q = h_i·(T_jet − T_g,i) 积分到融霜所需能量后处理,驱动温差被钉在最大值不衰减,融霜时间反而算短(偏乐观)。同一套设置、换一个后处理口径,方向就翻过来 ⇒ 方向不可预判。
读反之后,会做出三个很具体的错误动作:
- 拿"简化是保守的"去解释老模型与实测的偏差。模型算得比实测快,于是说"仿真偏乐观是正常的,留个系数就行"。而真因可能是进口边界拿了 HVAC 出风口温度,或者第一层网格落进了缓冲层——方向判错,就会去查错了地方,整改方向整条带反。
- 在偏乐观的绝热壁模型上做 DOE 并直接采信排序。模型系统性抬高内表面温度,把"勉强够"的方案判成"有裕量",A 柱根部的死角在模型里比实车浅;挑出来的最优方案在实车上不达标,而 DOE 报告上全是绿的。
- 把定壁温模型的"永远不达标"当成设计缺陷去整改。取环境温度作定壁温、又按内表面升到 0 ℃ 判读,任何方案都不达标,于是把风量与热源功率一路往上堆——那是模型的产物,不是设计的问题,而堆上去的风量还要付风噪与鼓风机工作点的代价。
正确做法只有一条:只有共轭壁面让内表面温度成为解出来的量——内侧对流、玻璃固体导热、外侧低温迎风对流三段串联同时在场,两头都对。⛔ 绝热壁与定壁温都不能当"保守估计"用:一个方向已知偏乐观,一个方向不可预判。要声明偏差方向,必须把"哪种理想化 + 取什么值 + 按什么判据读"三件一起声明,缺任何一件方向都不成立(⚠ 这条三件套交付口径同样是本课立的,⛔ 不是某个标准的要求)。
本课的射程边界与去向
本课只讲两件事:除霜的气流组织,以及共轭瞬态除霜 CFD 的方法学。以下几件相邻的事,本课一句不重讲,只写去向:
- 法规判据结构、视区作图链、能量底账(含玻璃显热)、以及"风量与温升数学上是乘积、工程上不可互换"这条结论——在 C5-01 前风窗除霜除雾法规与性能设计,本课引用其结论;
- 除雾侧的湿空气传输与凝结成膜数值方法——归 J4-03 除霜除雾风窗仿真;
- 空调箱内到仪表台下风道那一段的建模与内漏标定——归 J4-04 HVAC 空调箱总成级 CFD:箱内建模、开度扫描与内漏标定,本课第 3 讲的进口边界正取自与它的接界断面(⚠ 接界要说清的四件事是本课自己给的口径);
- 风道与出风分配的通用设计——在 C2-07 风道设计与出风分配(吹面/吹脚/除霜);鼓风机工作点在 C2-04 鼓风机与蜗壳设计:风量—噪声平衡,风噪与再生噪声在 C2-11 空调箱 NVH:风噪/风门异响/鼓风机噪声/再生噪声;
- CFD 通用判据(湍流模型选型、y⁺ 三档、网格质量、网格收敛指数、假收敛识别)——在 J2-01 CFD 基础、湍流模型与网格划分(Fluent/STAR-CCM+),本课只讲这些判据在除霜构型上怎么取值;
- 电加热风窗——在 C5-04 视野件电加热:加热风窗、后风窗除霜格栅与镜面/雨刮区加热的设计、控制与诊断(它在法规试验中作为除霜手段能不能计入那一层属法规判据结构,在 C5-01 前风窗除霜除雾法规与性能设计;⚠ 这一支的去向是两门,与 5.2、小结同口径),本课只在第 5 讲提一句"热源位置改变热阻路径 ⇒ 判据随之改变"。
标准方面,本课一律只写标准号与去向,⛔ 不写任何条款内容与限值——本课没有独立核对过它们的原文。已知的版次关系是:GB 11555-2025《汽车风窗玻璃除霜和除雾系统技术规范》全部代替 2009 版;GB 11562-2014 现行,GB 11562-2025 已发布并将自其实施日起代替 2014 版。⚠ 版次与实施日期以现行目录为准,引用前须核。凡是时间节点、分区面积门槛、试验环境温度与外侧气流条件这类量,本课一概不给,写成"须按现行版本核"。
同样地,本课对平台相关量一律不给数:h_i 与 h_o 及其空间分布、贴附行程与铺展宽度、脱流与对撞线位置、各支路的 A_i 与 ζ_i、除霜档风量与分区流量占比、换热关联式的系数与指数——它们全部取决于具体构型,须由本项目的仿真、1D 风道模型或实测给出。这不是本课偷懒:本课能交付的是"怎么把它算出来、怎么读它、怎么改它",而不是一组可以抄走的数。
七讲的走法
- 第 1 讲 除霜气流组织的物理:判据是贴不贴得住,不是吹得响不响。射流冲击与贴壁射流是分工不是二选一,脱流的两个诱因,以及中央对撞为什么在速度场与换热系数场上给出相反的结论。
- 第 2 讲 五个几何手柄:把流量分到 A / A′ 区。手柄是五个而不是三个(格栅位置与出流角、开口面积分布、局部阻力、导流叶片、风道走向),第六项"可制造性"不是手柄而是对前五项的侵蚀,单列。这一讲收在一条排序上:改分布优于加总量。
- 第 3 讲 进口边界的真值从哪来:除霜档的风量受系统约束、混合风门把温度与风量耦合在一起、模式门半开会把风漏到吹脚;以及 1D 与 3D 的分工——分配达标而分布不达标,是本课最常见的形态。
- 第 4 讲 数值设置:换热系数的系统性偏差从网格与壁面处理的配对来。这一讲有一个可以自己动手复算的算例,把 y⁺ 判据落到毫米上;也会讲一条反直觉的结论——把网格无脑加密可能让结果更差。
- 第 5 讲 物理理想化:三段串联热阻、为什么必须建成共轭壁面、为什么必须瞬态而不是稳态。上面那条反钉子在这一讲被当场点破。
- 第 6 讲 从换热系数场到设计动作:评价量的排序、按分区最低分位(⛔ 不是平均值)找末位死角、按本课归纳的八类成因归因(大纲给的三类都在气动线上,另外五类不在),以及 DOE 的三条前置检查与最优点的可制造性复核。
- 第 7 讲 凭什么信它:红外热像与粒子图像测速各校准证据链的哪一端、验证与确认是两件事、以及"什么交给 CFD、什么必须回环境舱"的三条互不相干的判据——能不能算、算了信不信、谁是法定证据。⛔ 三条不得互相推导。
最后是一次动手做:拿一版风窗与除霜风道数模,建含玻璃共轭壁面的瞬态模型、找出末位分位的死角并归因、出两版改进用 DOE 对比,交付改进前后的融霜面积-时间曲线与一页可制造性说明。走完这一遍,本课七讲的判据就全部用过一次了。
第 1 讲 除霜气流组织的物理:判据是贴不贴得住,不是吹得响不响
这一讲不讲任何一个可以填进设计输入的数,也不讲怎么建模型——它只做一件事:把「除霜风道设计得好不好」这个含糊的问句,收敛成一个逐点成立的式子和几条能判真假的判据。之所以要先做这一步,是因为除霜整改里最常见的僵局不是算错了,而是整改的人盯着的量,和决定成败的量根本不是同一个:会上看的是风量够不够、出风温度上没上来、速度云图漂不漂亮,而实际决定 A 柱根部那块霜化不化得开的,是那一点上的近壁换热系数。这两组量之间没有可以互相替代的关系,甚至可以给出方向相反的结论。
本讲的走法是:先把驱动量写成式子、看清气流组织到底能动哪两项(1.1);再把玻璃上同时存在的两种换热机制拆开——射流冲击(1.2)与贴壁射流(1.3),并给出它们各自的判据;然后说明两者是分工而不是二选一,以及为什么两套换热关联式绝对不能互套(1.4);接着把关联式这个骨架本身讲清楚,说明它为什么不是答案(1.5);再进入两处让覆盖失效的机理——脱流(1.6)与对撞回流(1.7),它们正是「速度云图好看而角落化不开」的来源;最后把全讲收口成一条本课自己立的换算规则(1.8),后面几讲的手柄、判据与实验设计都要挂在这条规则上比。
1.1 驱动量是一个逐点成立的式子:q(x,y) = h(x,y)·[T_jet(x,y) − T_g,i(x,y)]
是什么。玻璃内表面上任意一点 (x, y) 处,从空气侧进入玻璃的局部热流密度是:
q(x,y) = h(x,y)·[T_jet(x,y) − T_g,i(x,y)]
式中 q 为局部热流密度 [W/m²];h 为该点的内侧对流换热系数 [W/(m²·K)];T_jet 为该点近旁的射流当地温度 [℃];T_g,i 为该点的玻璃内表面温度 [℃](符号口径承 C5-01 前风窗除霜除雾法规与性能设计:玻璃温度一律拆成内表面 T_g,i 与外表面 T_g,o,全篇不出现不带下标的玻璃温度符号;q 表示热流密度 [W/m²],与热功率 Q̇ [W] 和按单位风窗面积记账的能量 Q [kJ/m²] 三者不混用——⚠ 除霜整门课的能量与热阻一样按 m² 记账,⛔ 与不带面积的总量 [kJ] 混算前必须先乘除面积因子,混用不会报错、只会静默错一个面积因子)。⚠ 这里的 (x, y) 是风窗曲面上的两个面内坐标,⛔ 不是壁面法向距离;本课的 y 分两义并在此分列:主式 q(x,y) 里的 y 是面内坐标,而第 4 讲 y⁺ 与 Δy 里的 y 是壁面法向距离,⛔ 两者不得互代。
这个式子有一个容易被滑过去的性质:它是逐点成立的,而且式子右边的两项都是场量,不是标量。h 是一个分布在整块风窗上的场,T_jet 也是——空气从格栅出来之后,一边沿风道往上走一边被玻璃夺热,越往前走温度越低。于是 q 也是一个场。每一点的融霜推进速度由该点自己的 h 与该点自己的驱动温差共同决定,与隔壁那一点算出来的值没有任何关系。
为什么这一步必须先做。因为读者手上的直觉判据是「风量」与「出风温度」这两个总量,而它们在这个式子里只决定上界:风量决定射流能带多少动量和多少焓进来,出风温度决定 T_jet 在源头能有多高。上界够高是必要的,但一点也不充分——真正被考核的那条「已除霜面积—时间」曲线(判据体系承 C5-01 前风窗除霜除雾法规与性能设计,本课不重讲),本质上问的是「区内有多大比例的面积,其局部 q 已经把该点的霜化完了」,这是一个分布问题。总量达标而分布不达标,是本课要处理的主要形态。
顺着这个式子往下看,还能读出一条对后面全部内容都有约束力的结论:气流组织这件事,真正能动的只有两项——
- h 的空间分布:把高 h 的区域挪到需要它的地方去;
- T_jet 的沿程衰减:包括风道内那一段(热损)与玻璃表面那一段(沿贴附行程被玻璃夺热)。
风速本身不在这个式子里。速度不是不重要,而是它只能通过 h 起作用;一旦速度与 h 之间的那条通路断掉(1.6 与 1.7 会给出两种断法),速度再高也不会兑换成 q。这一句就是本课的第一根钉子,1.7 末尾会正式收束它。
工程量级。h(x, y) 与 T_jet(x, y) 的具体分布是平台相关量,取决于格栅位置、出流角、除霜档风量与贴附状态,⛔ 本课一律不给数值、不给区间、也不给「典型值」。本课能交付的是另外三件事:怎么把这个场算出来(第 3~5 讲)、怎么读它(第 6 讲)、怎么改它(第 2 讲与第 6 讲)。
易错点。最常见的一处是把 q 写成 ṁ_air·c_p·ΔT_air(ΔT_air 是空气侧温升;⚠ 本课的 ΔT 另有专用含义——它只指「霜从初始温度升到融点的温差」,见 5.8 与 F6,所以这里必须带下标写全)这个总量式,然后据此在风量与温升之间做「等焓流互换」——鼓风机已经到顶了,那就把出风温度拉高一点补回来。⛔ 这条推理是错的,而且它错在一个很不显眼的地方:那个乘积式在数学上确实成立,C5-01 前风窗除霜除雾法规与性能设计 已经把它原样写出来过,也把它当成该课的反钉子当场点破了。本课直接继承它的结论——数学上是乘积,工程上不可互换;对本课尤其致命的是它的后半句:降风量首先塌掉的就是贴附覆盖,而贴附覆盖正是 1.3 要讲的那件事,也正是 A 柱两侧角落唯一的供给来源。换句话说,用降风量去换温升,换掉的恰好是最不该动的那一项。
1.2 射流冲击:驻点区的 h 是本构型的峰值,但它的判据是覆盖尺度
是什么。当射流以垂直或较大的角度撞上一块壁面时,在正对喷口的那一小片区域(驻点区)会发生一件特殊的事:主流不断把外部的新鲜空气推到壁面上,把刚被加热或冷却过的近壁流体挤走,热边界层被反复更新而始终维持在很薄的状态。边界层薄意味着壁面法向的温度梯度大,于是局部 h 在这里达到本构型能给出的最高值。离开驻点之后,流体转为沿壁面径向铺开,边界层重新开始增长,h 沿径向迅速衰减。
它的判据是什么。⚠ 不是「驻点 h 有多高」,而是驻点区的尺度 vs 要覆盖的区域的尺度。理由在 1.1 里已经埋好了:除霜的考核量是面积占比随时间,不是某一点最快化开。一个 h 极高但作用半径很小的驻点,对面积—时间曲线的贡献,完全可能小于一片 h 只有中等水平、却铺得开的贴附区。把峰值当成判据,会挑出一个「化出来几个亮点、周围一圈仍然是白的」的方案。
工程量级。驻点 h 的绝对值与它的覆盖半径都是构型相关量——它们取决于喷口的水力直径、冲击距离(喷口到壁面的法向距离)与冲击角,⛔ 本课不给任何数值。可以交付的只有一条量级上的定性关系:覆盖半径随喷口特征尺寸增大而增大,而以除霜格栅那样的开口尺度去比整块 A 区的尺度,前者通常小得多。这条量级关系的用途只有一个——说明「靠冲击啃满整块 A 区」这条路在几何上就走不通,⛔ 不得据它去反算任何覆盖面积或方案的可行性。
易错点。两个方向的错都常见:
- 拿平板或圆管的 Nu 关联式去估驻点 h。这是射程错误而不是精度问题——外掠平板与管内充分发展流的关联式,其特征长度、特征速度与流动结构与冲击射流完全不同(这一族关联式的选型自检见 A1-02 对流换热:自然对流、强制对流与换热系数估算)。代进去会得到一个数,而且这个数不会有任何东西提示你它是错的。
- 看到某方案驻点 h 很高就判它更好,而没有看它铺不铺得开。这正是把峰值当判据的形态。
1.3 贴壁射流与 Coandă 效应:A 区的覆盖主要靠它,判据是贴附行程与铺展宽度
是什么。当射流以较小的角度贴着壁面出流时,射流两侧的卷吸条件是不对称的:靠壁面那一侧的空气被壁面挡住、补充受限,形成局部低压;另一侧则可以自由卷吸。这个压差把射流拉向壁面并按住它,使射流沿壁面继续延伸——这就是 Coandă 效应。除霜格栅要的正是这个效果:让热风变成一层贴着玻璃铺开的「热风毯」。
但这个「按住」不是无限的。射流沿程一边卷吸周围较冷的空气、一边被壁面摩擦消耗动量,速度与温度都在衰减;走到某处动量不足以维持贴附,或者遇到 1.6 要讲的两种诱因,它就会离开壁面。因此贴壁射流的能力必须用两个几何量来度量:
- 贴附行程——射流能贴着玻璃走多远;
- 铺展宽度——射流沿横向铺开多宽。
它的判据是这两个几何量,⛔ 不是出口速度。出口速度高只说明源头动量大,说明不了这股动量能不能一直被壁面按住;1.6 会给出一个反例,说明速度和贴附可以完全脱钩。
为什么这一条对除霜特别关键。因为 A 区不是仪表台正前方那一小块。C5-01 前风窗除霜除雾法规与性能设计 已经把这条几何链讲透了:视区是从 V 点发出的角度锥投影到大倾角风窗曲面上,倾角越大投影拉伸越强,A 区一直铺到 A 柱根部。本课直接承接这个结论,并把它翻译成对气流组织的一条要求:射流必须贴得足够远、铺得足够宽,才够得着 A 区的外缘。这一步是把「造型定的玻璃倾角」这个纯外形参数,换算成「对射流的一条硬性要求」的那一步——倾角在造型冻结时就定了,而买单的是风道。
工程量级。贴附行程与铺展宽度都是平台相关量,取决于出口几何、出流角、玻璃倾角与除霜档风量,须由 CFD 或粒子图像测速 / 烟流显示实测给出,⛔ 本课不给数值。能给的量级关系只有一条,而且它是承 C5-01 前风窗除霜除雾法规与性能设计 的:射流动量随质量流量增长,降风量时贴附行程与铺展宽度会同时塌缩——注意是「同时」,这两个量不会一个塌一个不塌,所以降风量对覆盖的伤害是双重的。
易错点。把 Coandă 当成一句「只要贴着玻璃吹就一定能贴到底」的保证。⚠ 贴附是可以被打断的——壁面倾角突变与逆压梯度都能打断它,而且打断之后的后果不是「稍微差一点」,见 1.6。
1.4 两者是分工不是二选一:谁负责哪一块,以及⛔ 两套关联式一律不得互套
是什么。射流冲击与贴壁射流在同一块风窗上同时存在,只是各自主导不同的区域:冲击啃驻点附近那一小块局部硬骨头,贴壁负责把整片 A 区铺满。所以设计上要判的从来不是「我这套方案该走冲击路线还是贴壁路线」,而是哪一块该由谁负责,以及负责那一块的机制够不够用。
为什么必须把两者的关联式隔开。这不是一条讲究,而是一条硬禁令:⛔ 两套换热关联式的公式与取值规程一律不得互套。原因在于它们的特征量定义根本不同——
- 冲击射流关联式:特征长度取喷口水力直径,特征速度取喷口出口速度;
- 贴壁射流关联式:特征长度取沿壁行程,特征速度取当地射流最大速度。
把一套的 C、m、n 代进另一套的特征量里去,Re 本身就已经差了一到两个量级,再经过幂次放大,出来的 h 与真值不是差百分之几的关系。⚠ 这是量级级别的错误,不是精度问题——而它在计算书上看起来完全正常,因为每一步都是「按关联式算的」。
工程量级。⛔ 本课不给任何一组 C / m / n 的数值。它们必须按具体构型的公开文献或本项目实测标定,而且必须与它自己的特征长度、特征速度、物性取值温度成套引用——只抄三个系数不抄定义,等于抄了一个不知道自变量是什么的函数。
易错点。
- 抄系数不抄定义:这是上面那条禁令最常见的犯法形态,而且往往发生在转手一次之后——第一个人是从原文抄的、知道特征量,第二个人只从计算书里抄了三个数。
- 在同一张图上把两种机制的 h 曲线画在一起,却不标各自横轴的含义:两条曲线的横轴一个是径向距离、一个是沿壁行程,画在同一根轴上就是在暗示它们可以直接比大小。本课在图 1 上用分区着色与两套独立标签处理这件事,⛔ 刻意不画成一条连续曲线。
1.5 Nu = C·Re^m·Pr^n 是估算 h 的骨架,不是答案
是什么。上一节反复提到的「关联式」,指的是这样一族经验式:
Nu = C·Re^m·Pr^n
它把局部换热折算成 Re 与 Pr 的幂函数。三个无量纲数的定义是:
- Nu = h·L/k(无量纲)——L 为该关联式规定的特征长度 [m],k 为空气导热系数 [W/(m·K)];
- Re = u·L/ν(无量纲)——u 为该关联式规定的特征速度 [m/s],ν 为运动黏度 [m²/s];
- Pr = ν/α(无量纲)——α 为热扩散率 [m²/s]。
注意两处「该关联式规定的」——这不是修辞。L 与 u 不是通用的几何量和速度,它们是关联式定义的一部分,换一组定义,同一个物理构型算出来的 Re 就是另一个数。
适用条件与失效边界。每一组 C / m / n 只在它标定时的构型、Re 区间与边界条件下成立。换构型、换特征量定义、或者把 Re 外推出标定区间,误差都不是百分之几,而是倍数。这就是 1.4 那条禁令的定量根据。
工程量级。本课只在这里给出一个可查的物性口径:空气侧的 Pr ≈ 0.708(干空气、20 ℃、101.325 kPa,来源 CoolProp)。⚠ 这个口径必须连同两条限定一起读:① 真实除霜工况下,空气从 T_jet 一路被冷却到接近 T_g,i,物性严格说应取膜温 (T_jet + T_g,i)/2 现算,本课把它固定在 20 ℃ 是为了让后面的手算算例能被读者一步步复算——这是显示口径,不是工程取法,换温度整个算例会平移;② Pr 在本课的用途只有一条——说明空气侧 Pr 在这个工况附近近似为常数,所以 Nu 的工况依赖主要走 Re 这一项。⛔ 不得把它代入任何 Nu 算例去求 h。
至于 C / m / n 本身:⛔ 本课一个数都不给,也不代入这个式子算任何 h。理由很直接——C / m / n 未经本构型标定时代进去,会得到一个「看起来很专业的假 h」,而它后面每一步计算都会把这个假值原样传下去,没有任何一个环节会报错。
易错点。
- 抄系数不抄定义(同 1.4,这里给出它的定量后果:特征量换一套,Re 差一到两个量级)。
- 把关联式给出的平均 Nu 当成局部 Nu 用。这一条对除霜格外致命:除霜要的恰恰是局部分布,而平均值这个运算正好抹掉的就是死角——A 柱根部那一小块拖后腿的面积,在面平均里几乎不占权重。局部 h 与平均 h 的分野,A1-02 对流换热:自然对流、强制对流与换热系数估算 已有专节,本课不重讲,只强调它在除霜上的后果。
1.6 脱流:两个诱因,以及「速度云图好看而 A 柱根部化不开」的机理来源
是什么。贴壁射流沿着一块曲面前进时,如果壁面向背离射流的方向急转,或者沿程压力升高,近壁那一薄层流体的动量就可能不足以继续跟着壁面走。此时边界层从壁面分离,射流的主体离开玻璃继续向前,壁面附近留下的是一团低速滞止或回流的空气。这就是脱流。
两个诱因,⛔ 不许只数一个。
- 诱因一:玻璃倾角突变。典型位置就是 A 柱根部——风窗曲面在这里向两侧急转、曲率变化最剧烈,而它恰恰是 A 区外缘所在(1.3 已说明这条几何链承 C5-01 前风窗除霜除雾法规与性能设计)。最难贴的地方和最必须贴到的地方是同一个地方,这是除霜气流组织最不讨巧的一处。
- 诱因二:逆压梯度。射流沿程扩散、动量衰减时,或者下游存在滞止区、对撞区抬高了当地静压时,射流会顶着一个升高的压力往前走。这一支与几何无关,改倾角改不掉它。
为什么脱流之后 h 会塌。关键在于一句常被跳过的话:h 由近壁速度梯度决定,而不是由主流速度决定。贴附时,高速主流紧挨着壁面,壁面法向的速度梯度 ∂u/∂n 很大(⚠ n 是壁面法向坐标;⛔ 这里刻意不写成 ∂u/∂y——1.1 的 (x, y) 是玻璃面内坐标,与法向不是同一个方向,第 4 讲 y⁺ 与 Δy 里的 y 才是法向距离),边界层薄、换热强;脱流之后,主流可以继续以很高的速度掠过,但它已经不贴着玻璃了,玻璃表面那一层是低速回流,∂u/∂n 塌掉,h 随之塌掉。
⚠ 这就是本课那句「速度云图漂亮与 A 柱根部化不开可以同时成立」的完整机理。它不是一个玄学现象,也不是仿真精度不够——它是两个量之间的通路在脱流点被切断了。
工程量级,以及这里的一条严格自限。脱流的发生位置、以及脱流之后 h 的绝对值,都是平台相关量,⛔ 本课不给数。至于跌落幅度,本课引述的表述是「近壁换热系数出现量级跌落」——⚠ 这是一个定性的量级判断,本课对它加三条限制:① ⛔ 不给任何具体倍数;② ⛔ 不得据它反算脱流后的绝对 h、也不得反算任何融霜速率或时间;③ 它的用途被限定为唯一一条——说明「速度不是判据」这一件事。之所以要写得这么死,是因为 h 在本课是明确不给单点值的量,而一个不给单点值的量参与倍数比较,本身就踩在边界上;把用途锁死、把倍数拿掉,是这里能取的最保守写法。
易错点。最贵的一处是把脱流当成「风量不够」去整改。这个误判很自然:报告上写着「A 柱根部换热系数偏低」,最顺手的动作就是把风量往上加。但脱流位置主要由几何与压力梯度决定——对称地加风量,弱侧的相对劣势一点没变,两侧按同一比例变强,最慢的那一支还是最慢的那一支,达标节点仍由它决定;而付出去的是鼓风机工作点、风噪与整车功率的代价。正确的第一动作是改分布:导流叶片、开口面积与局部阻力的组合、格栅位置与出流角的配合——这些手柄第 2 讲逐个展开。
1.7 对撞与回流:速度场上的高值区,可能正是 h 场上的低值区
是什么。大多数布置里,除霜格栅是左右分布的,两侧射流沿玻璃向上铺开时会在风窗中央迎面相遇。两股流对撞之后,切向动量互相抵消,法向动量叠加,结果是整片边界层被抬离壁面,形成一股向上的自由射流;而被抬起来之后,壁面附近剩下的是低速滞止或回流。
于是中央区域出现一个很反直觉的组合:速度场上这里是高值区(两股流叠加,甚至比单侧还高),h 场上这里可能是低值区。「看似风大实则贴不住」说的就是它。
为什么要把它单独拎出来讲。因为它是本课第一根钉子最直观的一个反例:同一块区域,两种场图给出方向相反的结论。如果验收只看速度云图,对撞区会被判成「风量充足、没问题」,而它在 h 场上恰恰是需要处置的一块。1.6 的脱流讲的是「速度还在但 h 没了」,本节讲的是「速度更高了但 h 反而没了」——两者的机制不同,但落到判据上是同一句话。
工程量级。对撞线的位置与边界层被抬升的高度,都是平台相关量,取决于左右格栅的位置、出流角与两侧的流量分配比例,⛔ 本课不给数。
易错点。最典型的一处误动作是:为了「加强中央覆盖」,把左右格栅的出流角进一步向内偏。这个动作方向上就是错的——它让两股射流更正面地撞在一起,把对撞加剧、把边界层抬得更高。正确的处置方向是让两股射流错开(调整出流角的水平分量、重排开口分布,使两侧射流的主轴不在同一条线上迎头相遇),或者干脆接受中央那一带由驾驶员正前方的另一支射流负责覆盖,而不指望两侧射流在中央会合。
⚠ 这里还有一条判据上的纪律,必须显式写出来:中央对撞与 A 柱根部脱流是两个独立成因。它们各自成立、各有各的证据形态、各有各的处置动作,⛔ 不得互相推导——「中央对撞的问题解决了」推不出「A 柱根部也会跟着好」,反过来也不成立;图上同样 ⛔ 不得用一条引线把这两处连起来当因果链。把两个独立成因画成一条链,会让读者只处置其中一个就宣布收工。
钉子①,到这里正式收束。把 1.1 的式子和 1.6、1.7 两处机理放在一起,得到本课的第一根钉子:除霜快慢由玻璃内表面的局部净热流 q(x,y) = h(x,y)·[T_jet(x,y) − T_g,i(x,y)] 决定,而气流组织唯一能动的是 h 的空间分布与 T_jet 的沿程衰减。⛔ 风速不是判据——射流离开玻璃之后,速度可以照旧很高而 h 出现量级跌落;速度云图漂亮与 A 柱根部化不开,是可以同时成立的两件事。这条钉子在第 6 讲会以「评价量排序」的形态回来接受验收:近壁 h 与内表面温度 > 融霜面积—时间曲线 > 速度云图,速度场排在最后不是因为它不重要,而是因为它会给出相反的结论。
1.8 本讲收口:全部整改动作最终只通过一件事起作用——改变 h 的空间分布
到这里可以把整讲收成一条换算规则。
规则本身。本课后面要展开的全部手柄——格栅位置、出流角、开口面积分布 A_i、局部阻力系数 ζ_i、导流叶片、风道走向,以及第 3 讲的系统侧约束与第 6 讲的试验设计因子——最终都只通过同一件事起作用:改变 h 的空间分布(风道走向另有一条通路,它同时压低 T_jet,第 2 讲会单独处理)。于是每拿到一个手柄,都按同一对问句去过:
- 它改的是 h 的哪一块分布?
- 那一块是不是当前的末位分位区?
第二问不能省。改到了 h 的分布,但改的是本来就已经达标的那一块,等于没改——甚至更差,因为你为它付了成本,还占掉了一轮改型的机会。
为什么要立这条规则。因为不这么收敛,手柄之间根本没法比较:叶片角的单位是度,开口面积的单位是 mm²,风道走向压根不是一个数而是一种布置。三者放在一张表里,排序只能靠经验和嗓门。收敛到「改的是 h 的哪一块分布」之后,它们的收益就有了同一个度量——同一条融霜面积—时间曲线上的位移。这也是第 6 讲能把它们放进同一个试验设计里做因子的前提。
工程量级。⛔ 本课不给任何手柄的收益系数、不给任何「叶片角每变一度带来多少改善」这类换算。收益是构型相关的,而且它正是第 6 讲的试验设计要回答的问题——如果这里就能给出系数,那一整讲就不必存在了。
⚠ 一条出处声明,必须写明。上面这条「全部整改动作最终只通过改变 h 的空间分布起作用」是本课自己归纳出来的收敛口径,用途是把量纲各异的手柄放到同一把尺子上比。⛔ 它不是任何标准里的条款,也不是任何既有规范里已有的规定;跨课引用它时须说明这一点。本课后面还会出现若干条同样性质的自建判据,凡属这一类的,都会在当处写明。
与下一讲的接口。第 1 讲交付的是判据(什么算好、什么不算好),第 2 讲交付的是手柄(拿什么去改)。⚠ 这里先钉一句免得读错:能改流量分配与落点的几何手柄是五项而不是三项,另外还有一项不是手柄、而是对全部五项的侵蚀——它们的完整清单与各自的换算方式,是第 2 讲开场要做的第一件事。
后面还有 6 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做