停暖后房间为什么拖三天才凉：暖气片毕渥数0.00013，墙体1.29，指数曲线在这里根本不成立
停暖后房间为什么拖三天才凉：暖气片毕渥数0.00013，墙体1.29，指数曲线在这里根本不成立
几乎所有讲"房间多久凉透"的文章，用的都是同一条曲线：温差按指数衰减，时间常数 τ = 蓄热 ÷ 漏热。本站用时间常数 τ 算地暖该一直开还是出门就关那篇把这套算式做得非常完整，而且在第十一节老老实实写了一句话：严格说毕渥数 Bi = hL/λ ≈ 0.53 > 0.1，集总假设不成立，τ 只是工程近似，量级对就行。
这句话是对的，但它停在了最难的地方。本文接着往下问三件事：0.1 这个门槛是怎么定出来的、不成立的时候真实曲线长什么样、以及你手里只有一支温度计的时候怎么分辨。答案有点意外——墙那一半根本没有时间常数。
一、两个"凉"是两种数学
停暖之后，房间里同时发生着两件降温的事：
- 暖气片自己凉下来：钢板薄、导热快，内部温度处处相等，是教科书意义上的一阶惯性。
- 墙体把存的热慢慢吐出来：热要一层一层往外挪，是扩散过程。
前者有 τ，后者没有。把两者混在一条指数曲线里，就是"算出来三天该冻透、实测还停在 6 ℃"这个矛盾的来源。
二、毕渥数到底在比什么
Bi = h·L/λ，分子是表面把热交给空气的能力，分母是材料内部把热传开的能力。它比的其实是两个温差的比值：
- 外部温差：表面温度减流体温度，等于 q″/h；
- 内部温差：物体中心减表面，等于 q″·L/(2λ)（平板、长圆柱）或 q″·R/(3λ)（球）。
两式相除，热流 q″ 被约掉，得到一个干净的闭式：
内部最大温差 ÷ 表面-流体温差 = Bi/(n+1)，平板与长圆柱 n=1，球 n=2。
这就是 0.1 的全部来历：把 Bi 压到 0.1 以下，内部温差就只剩外部温差的 5%（球是 3.3%），此时"整个物体只有一个温度"这个假设才勉强站得住。门槛值不是拍出来的，是一道除法算出来的。
三、同一个房间里，Bi 差 9877 倍
取 h = 8.7 W/(m²·K)（GB 50736 内表面换热系数），特征长度按集总参数定义取 Lc = V/A：
| 对象 | λ W/(m·K) | Lc | Bi | 内部不均匀度 | 集总假设 |
|---|---|---|---|---|---|
| 暖气片钢板 1.5 mm | 50 | 0.75 mm | 0.0001305 | 0.0065% | 成立 |
| 铸铁片 3 mm | 52 | 1.5 mm | 0.0002510 | 0.0126% | 成立 |
| 100 mm 混凝土楼板 | 1.74 | 50 mm | 0.250 | 12.5% | 失效 |
| 3 cm 有效渗透层 | 0.50 | 30 mm | 0.522 | 26.1% | 失效 |
| 240 mm 黏土砖墙 | 0.81 | 120 mm | 1.289 | 64.4% | 失效 |
| 50 mm EPS 保温层 | 0.041 | 25 mm | 5.30 | 265% | 失效 |
钢板与砖墙之间差 9877 倍。暖气片本体是真正的集总体，所以用 τ 描述"停暖后暖气片多久变凉"完全合法（这与钢制暖气片停暖后余温时间常数那篇的算法一致）；而墙体内部温差高达外部温差的 64.4%，用单一温度描述它，等于默认墙心和墙面一样热——这显然不对。
EPS 那一行的 265% 也别当成错误：保温层本来就该把温差全扛在自己身上，它不作为蓄热体参与这场计算。
还有人会问 h 该取哪个。上表用的是内表面换热系数 8.7；如果按整栋外围折算（总漏热 82.69 W/K 摊到约 60 ㎡ 围护面积），h 只有 1.38，砖墙的 Bi 会降到 0.204。仍然大于 0.1。h 从 8.7 换到 1.38 缩小 6.3 倍，Bi 从 1.289 掉到 0.204，跨越门槛却没跨过去——这个判据在两种极端口径下给出同一个答案，比它在某一种口径下的具体数值更有说服力。
四、判据对 L 怎么取并不敏感
有人会问：L 取 3 cm（几小时尺度下的有效渗透深度）算出来 Bi = 0.522，取 12 cm（整墙半厚）算出来是 1.289，差 2.5 倍，结论会不会摇摆？
不会。两档都在 0.1 以上，而且都高出数倍。判据的失效区离门槛很远，L 的口径之争改变的是"失效多少"，不是"失不失效"。这一点很重要——它说明"墙不能用集总"不是参数取值造成的假象。
五、失效之后该用什么：误差函数
墙相对几小时的降温尺度足够厚，可以按半无限体处理。取 λ = 0.81、ρ = 1800、c = 880，得热扩散率：
α = λ/(ρc) = 5.114 × 10−7 m²/s
这个数与本站外墙内保温与外保温对比那篇表里黏土砖砌体的 a = 5.11 × 10−7 只差 0.071%，两条独立路径对得上。
停暖相当于在墙表面加一个阶跃，解有两种写法，取决于边界条件：
- 墙内部某深度 x（定壁温）：剩余比例 = erf(z)，z = x/(2√(αt))；
- 墙表面（对流边界，也就是你真正测到的地方）：剩余比例 = exp(H²)·erfc(H)，H = h√(αt)/λ。
两个式子里都没有 e−t/τ，也就不存在"经过一个 τ 走完 63.2%"这回事。63.2% 是指数函数的独有性质，换到误差函数上它就消失了。
六、最硬的一块砖：三个普适数
既然没有 τ，就用一个不依赖 τ 的办法衡量曲线形状：记下"降到一半"的时刻 t50 和"降到一成"的时刻 t10，取比值。
| 模型 | 剩余比例表达式 | t10 / t50 | 与材料有关吗 |
|---|---|---|---|
| 指数（一阶集总） | e−t/τ | 3.322 | 无关 |
| 扩散·墙内部 | erf(z) | 28.810 | 无关 |
| 扩散·墙表面 | exp(H²)·erfc(H) | 52.163 | 无关 |
这三个数值得单独说清楚。
3.322 就是 ln10/ln2，τ 被完全约掉，任何一阶系统都是它。
28.810 来自 t ∝ 1/z²，解 erf(z) = 0.5 得 z = 0.4769、erf(z) = 0.1 得 z = 0.08886，比值平方正好 28.810。α、x、λ 全部约掉了。
52.163 同理，解 exp(H²)·erfc(H) = 0.5 得 H = 0.7691、等于 0.1 得 H = 5.5546，(5.5546/0.7691)² = 52.163。有意思的是它对 h 也不敏感——h 从 0.5 取到 8.7（17 倍范围），比值恒为 52.163，因为 h 只改变时间轴的整体缩放。
三个数、三种曲线、与材料厚度换热系数统统无关。指数的 3.322 是三者的下限。
七、为什么 τ 在实测里经常"看起来是对的"
既然 Bi 早就超标，为什么工程上用 τ 几十年也没出大问题？两个原因，都不构成对 τ 的辩护。
第一，观察窗口太窄。只在开头 6~8 小时内取样时，erf 曲线和指数曲线的差别小于测温误差——前者在 t/t50 = 2 时剩 25%，后者剩 25%，几乎重合；分歧全在 t/t50 > 4 之后才拉开（18.38 对 3.00）。而绝大多数人测的是"关了两小时降了几度"，正好落在重合区。
第二，τ 是被拟合出来的，不是被算出来的。拿一条扩散曲线去套指数，总能解出一个 τ，代价是这个 τ 只在拟合区间内有效，换个区间它自己就变了。本站τ 那篇第十一节把"有效热容取 5 cm 还是 10 cm 渗透深度，τ 可差 2 倍"列成边界条件，其实就是这个症状——真正的物理常数不会因为你换了观测尺度就变。
八、等差与等比：一眼能看出的差别
把"剩余量每减一半需要多久"排成一行，两种模型的骨架就露出来了：
| 剩余 | 指数 t/τ | 指数 t/t50 | 扩散（墙内部）t/t50 | 扩散（墙表面）t/t50 |
|---|---|---|---|---|
| 50% | 0.693 | 1.000 | 1.00 | 1.00 |
| 25% | 1.386 | 2.000 | 4.48 | 7.12 |
| 12.5% | 2.079 | 3.000 | 18.38 | 32.81 |
| 6.25% | 2.773 | 4.000 | 73.99 | 136.09 |
| 3.125% | 3.466 | 5.000 | 296.42 | 549.38 |
指数这一列是 1、2、3、4、5，等差；扩散这两列越往后越接近 ×4，等比（表面版倍率 7.12 → 4.61 → 4.15 → 4.04 → 4.01，收敛到 4）。原因是小剩余时 erf(z) ≈ 1.128 z，于是 t ∝ 1/f²，剩余量减半就翻 4 倍。
这条"时间翻 4 倍、多跌一半"的规律，是重质房尾巴拖三天的真正机理，跟 τ 大不大没有直接关系。
九、能量账：一个会饱和，一个不会
从墙表面流出的累积热量也有闭式：
Q(0→t) = 2·ΔT·√(λρc)·√(t/π)，即正比于 √t。
其中 √(λρc) 叫吸热系数，黏土砖砌体算出来是 1132.7 J/(m²·K·s0.5)。把它换成规范里的蓄热系数：S24 = √(2π/T)·√(λρc) = 9.66 W/(m²·K)，与 GB 50176 表中黏土砖砌体的取值一致——又一条独立验证。
对照之下，指数模型的累积放热是 C·ΔT·(1−e−t/τ)，会饱和在 C·ΔT 上。扩散模型不会：24 小时放热系数 375.7 kJ/(m²·K)，96 小时（时间翻 4 倍）是 751.4，放热翻 2 倍。墙里的热是取之不尽的，因为它一直在往更深处挖。
十、具体到 240 砖墙的数字
| 内表面下深度 | 降到一半 | 降到一成 | 比值 |
|---|---|---|---|
| 3 cm | 0.54 h | 15.48 h | 28.81 |
| 5 cm | 1.49 h | 43.00 h | 28.81 |
| 10 cm | 5.97 h | 172.00 h | 28.81 |
| 20 cm | 23.88 h | 688.01 h | 28.81 |
比值恒为 28.81，深度只改变绝对时间。再看停暖 72 小时后各深度还剩多少余热：5 cm 处 7.74%、10 cm 处 15.40%、20 cm 处 30.23%、24 cm 处 35.89%。表层早凉了，深处还压着三分之一。
同一时刻，单一 τ 模型（取 τ = 29.8 h，本站停暖后房间多久凉透那篇的中质房取值）给出剩余 exp(−72/29.8) = 8.93%，折算室温 −5 + 25 × 0.0893 = −2.77 ℃，早该冻透了。可实测重质房三天后还有 6 ℃左右。差的不是快慢，是曲线的形状。
十一、落到你家：掐表测两个数就够了
找一个停暖日（或出门前关阀），室内 20 ℃、室外 −5 ℃ 时：
- 记下室温跌到 12.5 ℃（25 K 温差跌掉一半）的时刻 t50；
- 记下室温跌到 −2.5 ℃ 之前、也就是跌掉九成（剩 2.5 K）的时刻 t10；
- 相除。
| t10 / t50 | 说明 | 该用什么外推 |
|---|---|---|
| ≈ 3.3 | 换气/空气主导，真一阶 | τ 法完全可用 |
| 5 ~ 15 | 混合区，两者都在起作用 | τ 法只估量级 |
| 20 以上 | 墙体扩散主导，没有时间常数 | 用"时间翻 4 倍、多跌一半" |
大多数有混凝土或砖墙的房子会落在第三档。这时候别再用 τ 去外推"三天后多少度"，按 √t 律估：已经跌了 5 K 用了 10 小时，再跌 5 K 大约还要 30 小时，再跌 5 K 又要 120 小时。
反过来，测暖气片本体（用红外枪打散热片表面）会得到接近 3.3 的比值——它 Bi = 0.00013，是货真价实的指数衰减。同一个房间、同一个小时，两件东西服从两条不同的曲线，这就是 Bi 判据的意义。
一个算例
某砖混房实测：停暖后 3 小时室温从 20.0 ℃ 跌到 12.5 ℃（t50 = 3 h），跌到 −2.5 ℃（剩一成）用了 86 小时。比值 86/3 = 28.7，落在墙内部的 28.81 上。
用 τ 法会怎么做？由 t50 = 3 h 得 τ = 3/0.693 = 4.33 h，于是预测 24 小时后剩 exp(−24/4.33) = 0.39%，室温 −4.9 ℃——可实测 24 小时还有 5 ℃ 上下。错得离谱，而且错的方向永远是"低估尾巴"。改用 √t 律：3 小时跌一半，再跌一半（到剩 25%）需要 3 × 4.48 = 13.4 小时，再跌一半要 55 小时，正好对上那条 86 小时的曲线。
十二、本篇不算什么
- 不算推翻 τ 那套算法。它在估量级、比方案（出门几小时该不该关）时仍然好用，本文只是给出它的适用边界和边界外的替代品。
- 不算谐波响应。室外温度按 24 小时周期起伏的那套（衰减倍数 ν、延迟时间 ξ、渗透深度 δ = √(aT/π)）由外保温那篇讲完了，本文只处理停暖这种单次阶跃，不处理周期波动。
- 不含太阳辐射、户间传热。南向房间白天得热 20~40 W/m²，中间户降不到室外温度，这两项会让实测比值偏大，判读时按上表区间取。
- 空气本身不参与 Bi 判定：它不是导热体，靠换气混合，其时间常数就是 1/换气次数（0.5 次/h 对应 2 h），这一项永远是指数的。
所以完整答案是三段拼起来的：空气是指数（分钟到小时）、暖气片是指数（小时）、墙体是平方根（天）。房间降温曲线之所以永远拟合不好，是因为它本来就不是一条曲线。
相关产品推荐
结合本文场景，瑞华特可提供以下采暖产品和配套选型建议：
- 铜铝复合暖气片 TAL-800：铜铝复合暖气片TAL-800适合集中供暖、自采暖和水质复杂场景，兼顾耐腐蚀、快散热和长期稳定。核心参数为材质铜管+铝合金，高度800mm，工作压力1.2MPa。
- PPR采暖主管 DN25/DN32：用于壁挂炉、分水器和采暖主管连接，适合地暖与暖气片混合系统配套。
- 钢制暖气片 GZ-600：钢制暖气片GZ-600适合家庭集中供暖、自采暖和工程采暖场景，兼顾散热效率、耐用性和安装适配。核心参数为材质低碳钢，高度600mm，工作压力1.0MPa。瑞华特
需要确认规格、数量或工程采购清单，可联系瑞华特暖气片厂家获取免费选型指导。
