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 讲 | 一张必须逐项换掉的清单,以及漏换哪一项会错多少 |
| 我要建座舱湿空气 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。⇒ 越是法规关心的低温工况,这个错越严重。
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 ≥ 阈值且总除湿率 ≥ 产湿率」——三个字母都要有数,而这三个数只能由本课的仿真链算出来。
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 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做