TMS BOOK · ACADEMY 讲义

除霜除雾风窗仿真

大纲的完整展开版——讲师授课蓝本 / 学员自学材料

J4-03 除霜除雾风窗仿真

课程代码 J4-03 · 板块 J 仿真分析 CAE / J4 乘员舱与舒适性仿真 时长 约 4.0 小时(6 讲 + 1 次除雾仿真实操) 适合对象 仿真工程师,除霜除雾 / 座舱环境专项方向(选修,进阶);HVAC 控制标定工程师读第 6 讲 前置 J4-01 乘员舱气流组织与温度场 CFD 乘员舱气流组织与温度场 CFD;C5-03 低温起雾机理与主动防雾控制 低温起雾机理与主动防雾控制 本讲义定位 讲师授课蓝本 / 学员自学讲义,是本课题的完整展开版


引言:除霜和除雾长得像,可它们是两道不同的题——而仿真最贵的一类错,就发生在把它们当成一道题的时候

一份除雾仿真交出来的样子通常很干净:座舱云图有了,风窗内表面温度场有了,一张红黄绿的「结雾风险图」也有了,报告结论写着「A 区无结雾风险,除雾性能满足要求」。三个月后环境舱里做出来的结果是:满载四人、雨天上车,风窗下沿在第 4 分钟起雾,第 12 分钟才吹清。仿真侧于是开始返工,返工的第一个动作几乎是固定的——加密网格、换湍流模型、把风窗那块面的网格再切细一层。

这条路走到底也不会对,因为错的地方压根不在网格上。它在三个更早的地方:

第一,湿度这一头用的是全舱平均值。结雾是一个逐点成立的判据——某一点的玻璃内表面温度低于那一点的空气露点,那一点就结雾。而四个人的呼出气流会在风窗下沿的回流区聚成一团高湿空气,那团空气的露点比全舱平均值高出一两开尔文。一两开尔文听起来很小,可判据本身是一个减法,减法的两边差不多大的时候,一两开尔文足以让符号翻过来。

第二,凝结这件事根本没算。很多所谓「除雾仿真」只做到「判有没有结雾风险」为止,没有解液膜。可控制标定要问的从来不是「会不会起雾」——冬天上车必然起雾,这不需要仿真——而是「吹多久能吹清」。这个问题只有解了液膜的生长与蒸发才答得出来。

第三,判据的分母取错了。「A 区无结雾风险」这句话里的 A 区如果是拿整块风窗算的,那它系统性地把结论往好的方向偏——而风窗上最先结雾的那一块(下沿、A 柱根部)恰恰是面积最小、最容易被大分母稀释掉的。

这三件事有一个共同点:它们都不会让求解器报错。残差照样掉下去,守恒量照样闭合,云图照样光滑。这正是本课要处理的那一类问题——数值上完全健康、物理上系统性偏掉的仿真。

本课要钉的那一根钉子

本课交付的不是一张结雾风险云图,是一张能被 HVAC 控制标定直接查的判据表;而这张表能不能用,取决于中间那三次收窄做得对不对。

三次收窄是这样一条链:

湿空气场(组分输运解出的 Y_v 场)   → 裕度场(ΔT_margin = T_g,i − T_dew,逐点)   → 膜厚场(δ_f,由壁面液膜模型解出)   → 视区指标(η_vis,A,分母是法规视区 A 区)   → 阈值表(供 ECU 实时查)

每往右走一步,信息量都在减少,而每一步的减法都必须是可交代的:减掉了什么、减掉之后结论还站不站得住、误差往哪个方向偏。走完这五步而不能逐段交代,交出去的表就是一张看起来很确定、其实没人知道它的置信度是多少的表——控制标定会照着它标,标出来的策略在环境舱里对不上,然后没人说得清是哪一段出的问题。

与相邻课的分工,一次说清

本课在 CAE 侧,只讲除雾的数值实现。以下几件是相邻课的地盘,本课一句不重讲,只写去向:

这件事 属主课 本课怎么用它
除霜 / 除雾的法规判据结构、视区 A / A′ / B 怎么画 C5-01 前风窗除霜除雾法规与性能设计 直接取它的视区面积作分母,不重画
除霜的气流组织与共轭瞬态 CFD 方法学 C5-02 风窗气流组织与除霜风道 CFD 本课只讲「除雾与它差在哪」,除霜建模细节全去那边
起雾机理、露点裕度这个被控量、防雾控制闭环 C5-03 低温起雾机理与主动防雾控制 本课全盘沿用它的符号与算例工况,做 CFD 侧的实现
座舱 CFD 的地基(计算域、初场、辐射、湍流、网格) J4-01 乘员舱气流组织与温度场 CFD 本课默认它已经建对,只加水汽这一条方程
加热风窗、电加热镀层这类硬件手柄 C5-04 视野件电加热:加热风窗、后风窗除霜格栅与镜面/雨刮区加热的设计、控制与诊断 本课只把它当作一个内表面热流边界
湿空气物性、焓湿图、Magnus 式的来历 A2-04 湿空气物理与焓湿图:空调负荷计算基石 本课直接用,不推导

一句话分工:C5-03 低温起雾机理与主动防雾控制 回答「为什么起雾、怎么控」,C5-01 前风窗除霜除雾法规与性能设计 回答「合格线画在哪」,本课回答「这两件事怎么在 CFD 里算出来」。

符号与口径:本课沿用 C5 系列,只新增四个

跨课符号不统一是这条链上最贵的一种错——两门课各自都没写错,读者互代一次就全错了。所以先把口径钉死。⚠ 下表前六行沿用 C5-01 前风窗除霜除雾法规与性能设计 与 C5-03 低温起雾机理与主动防雾控制 已经定下的口径,本课一个字不改;后四行是本课新增。

符号 含义与单位 出处
T_g,i / T_g,o 玻璃内 / 外表面温度 [℃]。⛔ 全课不出现不带下标的 T_g、T_glass C5-01 前风窗除霜除雾法规与性能设计 / C5-03 低温起雾机理与主动防雾控制
T_dew 空气露点温度 [℃]。本课起它是一个场,写 T_dew(x) C5-03 低温起雾机理与主动防雾控制
ΔT_margin 露点裕度 = T_g,i − T_dew [K],恒为「离起雾还有几 K」 C5-03 低温起雾机理与主动防雾控制
d 含湿量 [g/kg 干空气] C5-03 低温起雾机理与主动防雾控制
S 面积 [m²]。⛔ 不用字母 A——A 在本系列专指视区 A 区 C5-01 前风窗除霜除雾法规与性能设计
δ_g / λ_g 玻璃厚度 [m] / 玻璃导热系数 [W/(m·K)] C5-01 前风窗除霜除雾法规与性能设计
Y_v 水蒸气质量分数 [kg 水汽 / kg 湿空气],⚠ 与 d 不是同一个量,换算见 2.2 节 本课新增
δ_f 壁面液膜等效厚度 [µm]。⚠ 下标 f = film,⛔ 与 δ_g(glass)不是一回事 本课新增
η_vis,z 视区 z 的可视面积比 = S_clear,z / S_zone,z × 100% 本课新增(分母口径承 C5-01 前风窗除霜除雾法规与性能设计)
ṁ_cond / ṁ_evap 单位面积凝结 / 蒸发质量流 [kg/(m²·s)] 本课新增

⚠ 一处跨课口径提醒,必须读:本课题的公式表里把结雾判据写作「T_glass,in ≤ T_dp」。这个写法在物理上没有错,但它与 C5-03 低温起雾机理与主动防雾控制 已经立住的符号不一致——那门课明确规定「正文一律带下标,⛔ 不裸写 T_glass」,露点写 T_dew。本课统一采用 C5-03 低温起雾机理与主动防雾控制 的写法:结雾判据写作 T_g,i ≤ T_dew,等价地 ΔT_margin ≤ 0。⇒ 读者在别处见到 T_glass,in / T_dp 这一对时,把它读成 T_g,i / T_dew 即可,两者指的是同一个量、同一个位置(玻璃内表面、局部空气露点),⛔ 不存在口径差。

⚠ 另一处跨课口径提醒,方向相反、必须区分:本课的 η_vis,z 与 C5-01 前风窗除霜除雾法规与性能设计 的已除霜面积占比形状一样(都是「某区里合格的面积 ÷ 该区总面积」),而分子的定义完全不同——除霜那边的分子是「霜已化掉」的面积,本课的分子是「没有起雾 / 雾已吹清」的面积,判读手段也不同(除霜看透光,除雾看散射)。⇒ 两个数不可互代,也不可放进同一张曲线图。本课全篇写 η_vis,z 并永远带下标 z(本课以 A 区为主,写 η_vis,A),就是为了在版面上就把它和除霜那条曲线分开。

本课的统一算例工况(全部教学假设值,⛔ 不对应任何平台)

为了让六讲的数字互相咬得住,本课从头到尾只用一组工况,且这组工况直接沿用 C5-03 低温起雾机理与主动防雾控制 的那一组——两门课的算例数字因此可以逐个对上,读者从那门课过来不用换脑子。

量 值 说明
舱内空气 22 ℃ / 50% RH ⇒ d = 8.2 g/kg,T_dew = 11.1 ℃
车外空气 2 ℃ / 95% RH ⇒ d_amb = 4.1 g/kg,T_dew,amb = 1.3 ℃
风窗玻璃 δ_g = 5 mm,λ_g = 1.0 W/(m·K) ρ = 2500 kg/m³,c = 840 J/(kg·K)
换热系数 h_i = 10、h_o = 30 W/(m²·K) 未开除雾风时的基准值
由上四行现算 T_g,i = 7.5 ℃ ⇒ ΔT_margin = −3.6 K 上车那一刻就已经在起雾
舱内干空气质量 3.6 kg(≈ 3.0 m³) 用于湿量平衡
乘员产湿 4 人共 200 g/h 静坐档
风窗面积 S_ws = 1.0 m²,其中 A 区 S_zone,A = 0.35 m² 用于视区指标
大气压 p = 101.325 kPa Magnus 式与含湿量换算的基准

⚠ 上表除大气压、水的相变潜热、玻璃物性区间之外,全部是教学假设值——它们的作用是让六讲的算例能互相校核,⛔ 不对应任何车型,不得作为任何设计输入或验收依据。真值须由整车实测与平台数据确定,拿到它的路见第 6 讲的环境舱标定一节。

左右两栏对照图,中间画一块共用的风窗玻璃剖面(厚度方向放大),玻璃两侧各标出内表面与外表面。左栏标题为除霜,画一条相变链:水汽在低于零度的外表面直接凝华成霜(气到固),再由内侧热流融化成水(固到液),链上两个箭头分别标注凝华与融化,并注明相变面钉在玻璃外表面、相变期间锋面温度恒为零度。左栏下方一个热阻串联条,自左到右三段分别是内侧射流膜阻、玻璃导热阻、外侧对流膜阻,标注净热流必须穿透前两段、并被第三段持续削减。右栏标题为除雾,画另一条相变链:水汽在低于露点的内表面凝结成液膜(气到液),再由除雾风蒸发回气相(液到气),链上两个箭头分别标注凝结与蒸发,并注明相变面在玻璃内表面、与内侧对流换热同在一侧。右栏下方同样画一个热阻条,但只有内侧那一段被高亮,标注凝结潜热直接释放在内表面、与对流同侧。图中央下方一个红框写出把霜层相变源项误挂到内表面的量化后果:正确挂法净热流七百九十三点三瓦每平方米,误挂内表面一千六百瓦每平方米,比值二点零,除霜时间由七点一分钟被算成三点五分钟。框内注明这两个数由本课统一教学假设工况现算。图右下角一句独立标注写明本图为机理与挂法对照示意,两栏热阻条按算例现算的相对长度作出,绝对值仅用于教学对照,不得据图读取任何车型的热阻、热流或除霜时间。
图1 除霜与除雾是两条不同的相变链,而且相变面分居玻璃两侧。要读出三件事:一,除霜是「气→固→液」两次相变、除雾是「气→液→气」一去一回,驱动力与潜热方向都不同;二,霜结在外表面、雾结在内表面,所以除霜的净热流要穿透玻璃并被外侧对流持续削减,而除雾的凝结潜热与内侧对流同在一侧;三,把霜层相变源项误挂到内表面,本课工况下净热流被算成 2.0 倍、除霜时间被算成一半(7.1 min → 3.5 min),而求解器不会报任何错。本图为机理与挂法对照示意,两栏热阻条按算例现算的相对长度作出,绝对值仅用于教学对照,⛔ 不得据图读取任何车型的热阻、热流或除霜时间。

按任务导航

手上正在做哪件事,就从哪一行进来。⛔ 本课不建议顺序通读之外再乱跳——每一讲都用到前一讲现算出来的数。

你手上的任务 直接读 你会拿到什么
别人给了我一份除霜模型,让我「改成除雾」 第 1 讲 一张必须逐项换掉的清单,以及漏换哪一项会错多少
我要建座舱湿空气 CFD,从零开始 第 2 讲 组分输运方程怎么加、源汇项怎么给、Y_v 与 d 怎么换算
报告说「A 区无结雾风险」,我要复核这句话 第 2 讲 2.4 节 + 第 6 讲 6.3 节 风险图的正确读法(看分位不看平均)与分母的正确取法
我要拿这套模型顺手算一下车内 PM2.5 / CO₂ 第 3 讲 能复用哪一部分、必须换掉哪两处,以及照抄会怎么错
我要算「吹多久能吹清」,不只是「会不会起雾」 第 4 讲 壁面液膜模型怎么选、薄膜假设什么时候失效
我要搭一个除霜 / 除雾共用几何的一体化算例 第 5 讲 共轭壁面怎么建、两侧源项怎么挂、时间步怎么定
我要把仿真结果交给控制标定,让他们能用 第 6 讲 三次收窄的完整链条与一张可查的阈值表模板
我要判断这份仿真值不值得信 常见误区一节 32 条可当场对照的检查项

⛔ 本课不覆盖——⚠ 本课已有一张「与相邻课的分工,一次说清」表(三列:这件事/属主课/本课怎么用它),⇒ 那张表是唯一事实来源,⛔ 本节不另写一份。要点:除霜/除雾的法规判据结构与视区 A/A′/B 怎么画 → C5-01 前风窗除霜除雾法规与性能设计(本课直接取它的视区面积作分母,⛔ 不重画);除霜的气流组织与共轭瞬态 CFD 方法学 → C5-02 风窗气流组织与除霜风道 CFD(本课只讲「除雾与它差在哪」);起雾机理、露点裕度与防雾控制闭环 → C5 板块。⇒ 本课在 CAE 侧只讲除雾的数值实现。


第 1 讲 除霜与除雾:两套不同的仿真问题

这一讲要交付的是一句听起来简单、做起来常错的判断:拿一份除霜 CFD 模型改成除雾模型,需要换掉的东西比大多数人以为的多。而漏换任何一项,模型都照样跑得完、收敛得好、云图漂亮。

1.1 两条相变链:类型不同、方向不同、潜热的量级也不同

先把物理摆清楚。霜与雾在中文里都是「玻璃上看不清了」,可它们的形成路径完全不同:

  • 结霜:水汽在温度低于冰点的表面上直接凝华成冰晶(气 → 固,跳过液相)。驱动力是水汽分压高于该温度下冰面的饱和分压。
  • 除霜:把已经存在的冰层融化(固 → 液),再由雨刮或气流带走。驱动力是热流。
  • 结雾:水汽在温度低于露点的表面上凝结成液(气 → 液)。驱动力是水汽分压高于该温度下水面的饱和分压。
  • 除雾:把已经存在的液膜蒸发回气相(液 → 气)。驱动力是壁面处的水汽分压差。

四个动作,三种相变,两个方向。把它们并排放到一张表上,仿真侧要建的东西立刻就不一样了:

结霜 除霜 结雾 除雾
相变类型 凝华 气→固 融化 固→液 凝结 气→液 蒸发 液→气
潜热 [kJ/kg] ≈ 2838(升华) ≈ 334(熔化) ≈ 2470(@10 ℃) ≈ 2470(@10 ℃)
相变面在哪一侧 玻璃外表面 玻璃外表面 玻璃内表面 玻璃内表面
驱动力 水汽分压差(对冰面) 热流 水汽分压差(对水面) 水汽分压差(对水面)
判据的量 — 已除霜面积占比 ΔT_margin ≤ 0 δ_f → 0 / η_vis,A
CFD 侧要解的 — 共轭导热 + 融化锋面 组分输运 组分输运 + 壁面液膜

⚠ 上表的三个潜热是物性值:水的熔化潜热 334 kJ/kg、10 ℃ 下的汽化潜热约 2470 kJ/kg、冰的升华潜热约 2838 kJ/kg(≈ 334 + 2504,即熔化加汽化),出自标准水物性表。⛔ 其余各行是本课整理的对照口径,不是任何标准的要求。

这张表上最容易被忽略的一行是「潜热」。汽化潜热是熔化潜热的 7.4 倍(2470/334)。这个倍数直接决定了两件事的难度不在一个量级上:融掉 1 g 冰要 334 J,而蒸发掉 1 g 水要 2470 J。所以「把雾吹干」在能量上比「把霜化掉」贵得多——只不过雾的量(微米级液膜)比霜的量(毫米级冰层)小得多,两者相乘之后总量才回到同一个量级。这也是为什么除雾的时间尺度是分钟,而除霜是数十分钟。

1.2 算例 1-A:相变源项挂错侧,除霜时间会被算成一半

这是本讲最重要的一个数,因为它量化了「两条链的相变面分居玻璃两侧」这句话的代价。

工况(教学假设值,与 C5-01 前风窗除霜除雾法规与性能设计 的除霜算例对齐):环境 T_amb = −18 ℃,除霜射流 T_jet = 40 ℃,除霜模式下贴附区内侧换热系数 h_i,defrost = 40 W/(m²·K),外侧 h_o = 30 W/(m²·K)(法规试验规定的迎面风速下),玻璃 δ_g = 5 mm、λ_g = 1.0 W/(m·K)。霜层结在外表面,相变期间融霜锋面温度 T_f = 0 ℃(恒定)。

正确挂法——相变源项挂在外表面。到达锋面的净热流要先穿过「内侧膜阻 + 玻璃导热阻」,到了之后还要扣掉外侧对流带走的那一份:

q_net = (T_jet − T_f) / (1/h_i,defrost + δ_g/λ_g) − h_o·(T_f − T_amb)

代入:

1/h_i,defrost = 1/40 = 0.0250 m²·K/W   δ_g/λ_g = 0.005/1.0 = 0.0050 m²·K/W   ⇒ 串联阻 = 0.0300 m²·K/W

导入项 = (40 − 0)/0.0300 = 1333.3 W/m²   外侧散失 = 30 × (0 − (−18)) = 30 × 18 = 540.0 W/m²   q_net = 1333.3 − 540.0 = 793.3 W/m²

错误挂法——把相变源项挂在内表面(这正是「拿除雾模板改除霜」时最常犯的一次操作,因为除雾的相变面确实在内表面)。此时玻璃导热阻与外侧散失都不在热流路径上:

q_wrong = h_i,defrost × (T_jet − T_f) = 40 × 40 = 1600.0 W/m²

后果。倍数 = 1600.0 / 793.3 = 2.02。取 C5-01 前风窗除霜除雾法规与性能设计 现算的除霜能量底账 Q_A = 337.4 kJ/m²(含霜显热、霜潜热、玻璃显热三笔,其中霜面密度 0.4 kg/m² 是教学假设值):

正确:t = 337 400 / 793.3 = 425 s = 7.1 min   错挂:t = 337 400 / 1600.0 = 211 s = 3.5 min

这个错误让除霜时间少了一半,而它在数值上完全健康。残差正常、能量守恒(在错误的那本账里守恒)、锋面推进曲线光滑。要发现它,只有一条路:在后处理里把外侧对流那一项单独积分出来看它是不是零。是零,就说明源项挂错了侧——玻璃外表面根本没参与传热。

⚠ 两个方向性的补充。第一,误差永远偏乐观:漏掉的两项(玻璃导热阻、外侧散失)都是减少净热流的,漏掉它们必然把热流算大。第二,外温越低这个倍数越大:把 T_amb 换成 −8 ℃,外侧散失变成 240 W/m²,q_net = 1093.3,倍数掉到 1.46;换成 −25 ℃,散失 750 W/m²,q_net = 583.3,倍数升到 2.74。⇒ 越是法规关心的低温工况,这个错越严重。

一张两栏手柄方向对照表的图形化版本。左栏标题为除霜工况,工况条件写明约零下十八摄氏度、压缩机被低环温闭锁、外气几乎无湿可除。右栏标题为除雾工况,工况条件写明零到十摄氏度、压缩机可运行、舱内湿源活跃。表体共六行,每行一个手柄,行首写手柄名,左右两栏各画一支箭头表示该手柄在该工况下应往哪个方向拧,箭头旁写一句理由。六个手柄依次为:热源功率,左栏箭头向上标注热量全由正温度系数加热器与暖风芯体提供、是唯一驱动力,右栏箭头向上但标注受限于抬高玻璃温度与压低再热温升之间的权衡;风量,左栏箭头向上标注加大射流换热系数直接加大净热流,右栏画一支带问号的双向箭头标注方向不唯一、见第一讲算例,并加一个感叹号标记;内外循环,左栏标注不敏感因为外气无湿,右栏箭头指向外循环并标注低温外气本身就是除湿手段;空调压缩机,左栏画一个禁止符号标注被低环温闭锁、无法除湿,右栏箭头向上标注蒸发器表面除湿是主力;出风分配,左右两栏箭头都指向风窗但标注除霜要覆盖到玻璃两侧下角、除雾要优先覆盖驾驶员视区;再热温升,左栏箭头向上标注越高越好,右栏画一支带问号的双向箭头标注与风量此消彼长。表格底部一个红色横条写出本讲的结论句:六个手柄里有两个在两种工况下方向相反或不唯一,把除霜的边界条件原样拿来跑除雾必然错,两个不唯一的手柄旁各引一条线指向这句话。图右下角一句独立标注写明本图为手柄方向的定性对照,箭头只表示方向不表示幅度,图上无数值轴无刻度,不得据图读取任何风量、温升或功率数值。
图2 除霜与除雾的手柄方向不是同一套,其中两个手柄的方向甚至不唯一。要读出三件事:一,压缩机这一行在两栏里是「禁止」对「主力」——低温除霜工况下压缩机多被低环温闭锁,除湿这条路整个不存在;二,风量与再热温升这两行在除雾栏里是带问号的双向箭头,因为抬风量同时干三件方向相反的事(详见 1.3 节算例);三,正因为有这两个不唯一,「拿除霜的边界条件跑一遍就当除雾做过了」必然错。本图为手柄方向的定性对照,箭头只表示方向、不表示幅度,图上无数值轴、无刻度,⛔ 不得据图读取任何风量、温升或功率数值。

1.3 算例 1-B:除雾时抬风量,同时干三件方向相反的事

「除雾风量开大一点总没错」——这句话在工程上非常常见,而它是错的,或者更准确地说,它的对错取决于三笔账相加的结果,而这三笔账方向不一致。

抬高除雾风量,同时发生三件事:

  • ① 内侧换热系数 h_i 上升 ⇒ 玻璃内表面温度 T_g,i 被拉向送风温度 ⇒ 对 ΔT_margin 有利;
  • ② 再热温升 ΔT_reheat 下降(热源功率给定时,Q̇ = ṁ·c_p·ΔT,ṁ 涨则 ΔT 跌)⇒ 送风温度降低 ⇒ 对 ΔT_margin 不利;
  • ③ 蒸发器表面温度抬高(同样蒸发功率摊到更大风量上)⇒ 单位质量除湿量下降 ⇒ 对 T_dew 不利。

把三笔账算出来。工况(教学假设值):热源功率上限 Q̇_heat = 2.0 kW,空气 c_p = 1.006 kJ/(kg·K)、ρ = 1.2 kg/m³,外温 2 ℃,玻璃与外侧同统一工况(δ_g/λ_g = 0.0050 m²·K/W、h_o = 30)。基准送风温度取舱温 22 ℃ 之上再加再热温升。

工况 A:风量 300 m³/h

ṁ = 300 × 1.2 / 3600 = 0.1000 kg/s   ΔT_reheat = 2000 / (0.1000 × 1006) = 19.9 K ⇒ T_jet = 22.0 + 19.9 = 41.9 ℃

工况 B:风量 500 m³/h

ṁ = 500 × 1.2 / 3600 = 0.1667 kg/s   ΔT_reheat = 2000 / (0.1667 × 1006) = 11.9 K ⇒ T_jet = 22.0 + 11.9 = 33.9 ℃

换热系数按湍流冲击射流的经验标度 h_i ∝ V^0.8(行业常用指数,冲击射流一般在 0.7~0.8,⛔ 具体值须由本平台风道 CFD 或台架标定):

h_i,B = 10 × (500/300)^0.8 = 10 × 1.473 = 14.7 W/(m²·K)

算 T_g,i(三段串联,驱动温差取 T_jet 到 T_amb):

工况 A:R_i = 1/10 = 0.1000,R_g = 0.0050,R_o = 0.0333 ⇒ R_tot = 0.1383

T_g,i = T_jet − (T_jet − T_amb)·R_i/R_tot = 41.9 − 39.9 × 0.7229 = 41.9 − 28.8 = 13.1 ℃

工况 B:R_i = 1/14.7 = 0.0680,R_g = 0.0050,R_o = 0.0333 ⇒ R_tot = 0.1064

T_g,i = 33.9 − 31.9 × (0.0680/0.1064) = 33.9 − 31.9 × 0.6396 = 33.9 − 20.4 = 13.5 ℃

第一个结论:风量抬了 67%,T_g,i 只涨了 0.4 K。因为 ① 和 ② 几乎正好抵消——换热系数涨的那部分,被送风温度跌的那部分吃掉了。

再算第三笔账(除湿)。设蒸发器出口空气接近饱和,风量 300 m³/h 时表面温度 3 ℃、500 m³/h 时抬到 7 ℃(教学假设值;实际抬升量由蒸发器 UA 与制冷量决定,⛔ 须由系统仿真给)。查本课统一算式(Magnus + 含湿量换算,p = 101.325 kPa):

d_sat(3 ℃) = 4.7 g/kg  d_sat(7 ℃) = 6.2 g/kg

内循环进风 d = 8.2 g/kg,则单位质量除湿量:

工况 A:8.2 − 4.7 = 3.5 g/kg   工况 B:8.2 − 6.2 = 2.0 g/kg ⇒ 掉了 43%

总除湿率:

工况 A:0.1000 kg/s × 3.5 g/kg = 0.350 g/s = 1.26 kg/h   工况 B:0.1667 kg/s × 2.0 g/kg = 0.333 g/s = 1.20 kg/h

第二个结论:风量抬了 67%,总除湿率反而跌了 5%。

合起来:从 300 抬到 500 m³/h,玻璃那一头只赚了 0.4 K,空气那一头亏了 5% 的除湿率——而付出的是 67% 的风机功耗与相应的噪声。⇒ 「除雾风量越大越好」不成立。⚠ 同样要注意,本讲也没有证明「风量越小越好」——从 300 往下降,h_i 掉得很快,T_g,i 会跟着掉。正确的说法只有一句:除雾风量存在一个由「热源功率上限 × 蒸发器 UA × 风道换热特性」共同决定的最优点,它必须靠寻优找出来,而寻优的目标函数是 ΔT_margin 与总除湿率的组合,⛔ 不是风量本身。这也正是本课要交付的那张判据表存在的理由。

⛔ 特别提醒:把这一节的结论写成「适度的风量」是没有信息量的。工程上可交付的说法必须是:「在本平台 PTC 功率 X kW、蒸发器 UA = Y 的条件下,除雾风量的最优区间是 Z1~Z2 m³/h,判据是 ΔT_margin ≥ 阈值且总除湿率 ≥ 产湿率」——三个字母都要有数,而这三个数只能由本课的仿真链算出来。

一张三栏并列的折线图,三栏共用同一条横轴,横轴为除雾风量,范围三百到五百立方米每小时,标出三百与五百两个刻度点。左栏纵轴为再热温升,单位开尔文,画一条随风量单调下降的曲线,标出两个算例点,三百对应十九点九开尔文、五百对应十一点九开尔文,曲线旁一句标注写明热源功率给定时温升与质量流量成反比,并写出公式热功率等于质量流量乘比热乘温差。中栏纵轴为玻璃内表面温度,单位摄氏度,纵轴范围刻意放窄到十二到十四摄氏度以便看清变化,画一条几乎水平的曲线,标出两个算例点,三百对应十三点一摄氏度、五百对应十三点五摄氏度,两点之间画一支双箭头标注只涨了零点四开尔文,曲线上方一个文字框写明换热系数上升带来的好处几乎被送风温度下降完全抵消,并列出两栏的三段热阻数值。右栏纵轴为总除湿率,单位千克每小时,画一条先平后微降的曲线,标出两个算例点,三百对应一点二六千克每小时、五百对应一点二零千克每小时,两点之间画一支向下的短箭头标注反而跌了百分之五,曲线旁一个文字框列出单位质量除湿量由三点五克每千克跌到二点零克每千克、跌幅百分之四十三,以及蒸发器表面温度由三摄氏度抬到七摄氏度这个假设。三栏顶部一条横贯的标题带写明抬风量同时干三件方向相反的事。三栏底部一条横贯的结论带写明风量越大越好不成立,风量越小越好同样不成立,最优点必须寻优。图右下角一句独立标注写明三条曲线的两个端点由本课统一教学假设工况现算,曲线中段为示意插值,蒸发器表面温度的抬升幅度是教学假设值,须由系统仿真确定,不得据图读取任何车型的风量最优点、温升或除湿率。
图3 抬高除雾风量同时干三件方向相反的事,净结果不是「更好」。要读出三件事:一,左栏——热源功率给定时再热温升与风量成反比,300→500 m³/h 让温升从 19.9 K 掉到 11.9 K;二,中栏——换热系数涨的好处被送风温度跌的坏处几乎完全抵消,T_g,i 只从 13.1 涨到 13.5 ℃(+0.4 K,注意纵轴是放窄过的);三,右栏——蒸发器表面温度被抬高后单位质量除湿量掉 43%,总除湿率反而跌 5%。⇒ 结论只有一句:最优风量必须寻优,⛔ 不能靠「开大一点」。三条曲线的两个端点由本课统一教学假设工况现算,曲线中段为示意插值,蒸发器表面温度的抬升幅度是教学假设值、须由系统仿真确定,⛔ 不得据图读取任何车型的风量最优点、温升或除湿率。

1.4 判据的分母不同:一个看「化开了多少」,一个看「看得清多少」

第 1 讲的第三处差异在判据侧,它比前两处更容易被忽略,因为两个判据长得几乎一样。

  • 除霜的判据是一条已除霜面积—时间曲线:多个时间节点 × 分区(A / A′ / B)× 各自的面积百分比门槛,逐节点全过才算过。判读手段是看玻璃透不透光(霜化了就透)。判据结构与门槛的属主是 C5-01 前风窗除霜除雾法规与性能设计。
  • 除雾的判据是可视面积比 η_vis,z:某个视区里「看得清」的面积占该视区面积的百分比。判读手段是看散射(有雾就散射,透光率反而可能还不低)。

三处必须分开:

① 分母都是视区,不是整块风窗——这一点两者相同,而且都是最常错的一处。把风窗四周的黑边、下沿被仪表台遮挡的部分、以及 B 区那些「看得见但不用于主要驾驶视野」的面积算进分母,会把最先失效的那块小面积稀释掉。⚠ 视区怎么画(座椅 R 点 → V1/V2 点 → 视区角度锥 → 与风窗曲面求交)是 C5-01 前风窗除霜除雾法规与性能设计 的地盘,本课直接用它算出来的 S_zone。

⛔ 这里有一处高频错,本课点名:眼椭圆(SAE J941 那一族)不投影进除霜除雾视区。它和 V 点常在同一套视野数模里建、图层看起来平级,但它服务的是后视镜视野、抬头显示、遮阳板遮蔽与 A 柱障碍角。拿眼椭圆去截风窗得到的那块区域,与法规视区不是同一块——而拿错了包络,仿真与试验都不会报错,它们会老老实实地把一块错的区域算得很精确。

② 除雾以 A 区为主,⛔ 不要照搬除霜的 A′/B 分时段门槛。除霜的考核是「多节点 × 多区域」的一张表,除雾的考核口径与它不同,本课以 A 区为主。⚠ 具体的门槛数值、时间节点、各区适用范围,以现行版本条文为准——GB 11555《汽车风窗玻璃除霜和除雾系统技术规范》与 FMVSS 103(Windshield Defrosting and Defogging Systems)都在其中规定了考核方法,版次与实施日期以现行标准目录为准。⛔ 本课不复述任何具体数值门槛,因为一旦复述,下游就会当限值抄进需求书,而不会有人再回去核它的出处。

③ 分子的定义不可互代。除霜的「已除霜」是一个相态判断(冰化没化),除雾的「已除清」是一个光学判断(还散不散射)。而光学判断这一头不是膜厚的单值函数——同样多的水,摊成连续水膜和聚成一片小液滴,可见度差得很远。这一条是第 6 讲的主题,本讲只先把话放在这里:δ_f 是凝结量的代理量,⛔ 不是可见度本身。

1.5 一张「除霜模型改除雾」的换件清单

把前四节压成一张可以逐行打勾的表。⚠ 这张表是本课整理的工程口径,⛔ 不是任何标准的要求。

# 要换的东西 除霜里是什么 除雾里必须换成什么 漏换的后果
1 相变源项挂在哪一侧 玻璃外表面 玻璃内表面 净热流算成 2.0 倍(算例 1-A)
2 要不要解组分输运 不用(外气几乎无湿) 必须解 Y_v 根本算不出结雾风险
3 要不要解壁面液膜 融化锋面(另一类模型) 必须解 δ_f 答不出「吹多久能清」
4 环境温度 约 −18 ℃ 量级 0~10 ℃ 量级 外侧膜阻与凝结驱动力全错
5 压缩机状态 多被低环温闭锁 运行,除湿主力 湿量平衡整个记反
6 内外循环位 不敏感 敏感,是主手柄 漏掉最有效的一支除湿手段
7 舱内湿源 可忽略 必须建(呼出 + 湿衣物) 露点场整体偏低,风险图偏乐观
8 玻璃内表面涂层状态 影响小 影响可见度阈值(第 6 讲) 判据阈值差数倍
9 判据的分子 已除霜面积(相态) 可视面积(光学) 两条曲线画进一张图,读者互代
10 时间尺度 数十分钟 分钟 时间步与总时长都定错

十项里,前三项是模型结构级的改动(要开新的物理模型),中间四项是边界条件级的,后三项是后处理与判据级的。⇒ 「把除霜模型改成除雾模型」不是改几个数,是重建一半的模型。

1.6 本讲的易错点

易错点 1-①:把霜层相变源项挂到内表面。 形态:从除雾模板起手改除霜(或反过来),源项位置没跟着换侧。 为什么会错:两个模型的几何完全一样、面组名也一样(都是「风窗玻璃」那一组面),前处理界面上「内表面」和「外表面」只差一个字,而选错之后没有任何一个求解器设置项会变红——两侧都是合法的壁面,都能挂源项。 怎么当场发现:在后处理里对外表面单独做一次热流面积分。除霜工况下它必须是一个显著的负值(向环境散热,本课算例是 −540 W/m²);如果它是零或接近零,源项一定挂错了侧。

易错点 1-②:拿除霜的边界条件跑一遍就当除雾做过了。 形态:一份报告里两个工况用同一套 HVAC 边界(同样的风量、送风温度、内外循环位)。 为什么会错:除霜工况下压缩机被低环温闭锁、外气几乎无湿可除,所以「内外循环」这个手柄在除霜里根本不敏感,工程师习惯性地不去动它;换到除雾工况,这个手柄变成了最有效的一支除湿手段(低温外气的含湿量只有舱内的一半——本课工况 4.1 对 8.2 g/kg),不动它等于把主力扔掉了。 怎么当场发现:翻边界条件表,看「内外循环位」「压缩机开关」这两行在两个工况里是不是同一个值。是同一个值,就一定有一个是错的。

易错点 1-③:把「除雾风量」当成越大越好的单调手柄。 形态:整改建议写「加大除雾风量」,没有给上界,也没给判据。 为什么会错:抬风量同时干三件事(算例 1-B),其中两件是不利的——热源功率给定时再热温升按 Q̇ = ṁ·c_p·ΔT 反比下降,蒸发器表面温度被抬高使单位质量除湿量下降。三件相加之后净收益可能是负的。 怎么当场发现:要求提出者给出两个数:抬风量后的送风温度、抬风量后的蒸发器出口露点。给不出这两个数,这条建议就只算了三笔账里的一笔。

易错点 1-④:拿眼椭圆截风窗当成法规视区。 形态:视区面积由视野数模里的「眼椭圆图层」投影得到。 为什么会错:眼椭圆与 V 点同属驾驶员视野包络工具族,常建在同一套数模里、图层平级,可它服务的对象完全不同(后视镜视野、抬头显示、遮阳板、A 柱障碍角),不投影进除霜除雾视区。而拿错包络之后,仿真与试验都会把一块错的区域算得非常精确,不会报任何错。 怎么当场发现:问一句「V1/V2 点的坐标是多少、从哪个 R 点推出来的」。答不上来,用的多半就是眼椭圆。

易错点 1-⑤:把已除霜面积曲线与可视面积比画进同一张图。 形态:一页报告上两条百分比曲线并列,共用一根「面积占比 %」纵轴。 为什么会错:两条曲线的分子定义不同——一个是相态判断(冰化没化),一个是光学判断(还散不散射)。而它们的分母口径又恰好相同(都是视区面积),于是在版面上长得一模一样,读者会不自觉地互相比较甚至互相代入。 怎么当场发现:看图例里有没有写明「分子是什么、怎么判读的」。只写「A 区面积占比」而没写判读手段的,就是这个毛病。

易错点 1-⑥:把汽化潜热和熔化潜热记混,导致除雾能耗被低估。 形态:估算「吹干液膜需要多少能量」时按 334 kJ/kg 算。 为什么会错:334 是熔化潜热(冰→水),而吹干液膜是汽化(水→气),10 ℃ 下约 2470 kJ/kg,差 7.4 倍。这个混淆特别容易发生在同一份报告同时处理除霜与除雾的时候——两个数就在相邻的表格里。 怎么当场发现:看单位面积能量的量级。本课工况下 10 min 凝结的液膜约 11.5 µm(算例 4-A),即 11.5 g/m²,汽化需 11.5 × 2470 = 28.4 kJ/m²;若算成 3.8 kJ/m²,就是记混了。

后面还有 5 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做

会员专属

后续为会员深水区内容——四库数据与深度拆解。

查看会员方案