PMV/PPD温热舒适性模拟
PMV/PPD温热舒适性的理论基础
PMV模型的概况
老师，我听说PMV能将"舒适性数字化"，真的吗？舒适度不是主观的吗？
问得好。PMV（预测平均投票）是丹麦工程师Ole Fanger在1970年提出的指标，本质上是"用数学模型统计预测大多数人的主观感受"。
怎样才能把主观感受变成公式呢？
Fanger进行了一项规模庞大的实验——1396名被测者在各种温度、湿度、风速条件下，用-3（冷）到+3（热）的7级量表评价他们的感受。然后他把这些统计数据与人体热收支方程结合起来。简单说，如果人体产生的热量与环境散发的热量平衡，就是PMV=0（中立＝舒适）；如果热量过剩就是PMV>0（热），热量不足就是PMV<0（冷）。
这个模型在什么场景下应用呢？
最常见的是办公楼和商业建筑的空调设计。在设计阶段，用CFD（计算流体力学）计算室内温度、气流分布，然后在各点计算PMV并画成等高线图，一眼就能看出"这个座位太冷""窗边太热"。ISO 7730标准规定PMV在±0.5以内才算舒适。
6因子的物理意义
我听说PMV有6个影响因子。各自代表什么意思？
6个因子分为环境方面4个和人体方面2个。这是很重要的区分。
| 因子 | 符号 | 单位 | 典型值（办公） | 物理意义 |
|---|---|---|---|---|
| 空气温度 | $T_a$ | °C | 22~26 | 直接影响对流传热 |
| 平均辐射温度 | $\bar{T}_r$ | °C | 20~28 | 周围表面红外辐射的合成效果 |
| 相对风速 | $v_a$ | m/s | 0.1~0.3 | 气流带走体表热量的速率 |
| 水蒸气分压 | $p_a$ | Pa | 1000~1500 | 出汗蒸发的驱动力（湿度影响） |
| 代谢率 | $M$ | W/m² | 58~70（坐位） | 人体内部产热量。1 met = 58.15 W/m² |
| 着装量 | $I_{cl}$ | clo | 0.5~1.0 | 衣服的隔热性能。1 clo = 0.155 m²K/W |
"met"和"clo"是特殊单位。具体数值是多少？
1 met对应"安静坐着"的代谢量，等于58.15 W/m²。走路的话约2 met，轻度运动约3 met。着装量的1 clo大约是"标准商务套装一套"。夏季清凉着装约0.5 clo，冬季穿外套约1.5 clo。这两个参数对结果影响很大，必须在设计条件中明确指定。
PMV式的公式化
PMV的公式是什么样的？
PMV的核心是求"人体热负荷$L$"，然后把它转换成主观评分尺度。整体形式是这样的：
其中$L$是人体热负荷（内部产热减去各种散热途径），展开为：
哇，这是个很长的式子…各项分别代表什么？
按顺序讲解：
- $(M - W)$：代谢率减外部工作，等于纯粹的内部产热。办公工作中$W \approx 0$
- 第2项：从皮肤不感蒸发（你察觉不到的蒸发）散失的热
- 第3项：出汗蒸发散失的热。当代谢超过58.15 W/m²（1 met）时开始出汗
- 第4项：呼吸潜热。吸入的空气和呼出的空气中水蒸气含量差异导致的散热
- 第5项：呼吸显热。吸入的冷空气被加热到体温需要消耗的热量
- 第6项：衣服外表面的辐射热交换。遵循Stefan-Boltzmann定律
- 第7项：衣服外表面的对流热交换
其中$f_{cl}$是衣服面积比（衣服增加多少表面积），$T_{cl}$是衣服表面温度，$h_c$是对流传热系数。$T_{cl}$本身需要通过迭代法从$T_a$、$\bar{T}_r$、$I_{cl}$计算，这是程序实现中最麻烦的地方。
PPD（预测不满意率）
从PMV还能预测"不满意的人的比例"？PPD是什么？
PPD（预测不满意率）就是根据PMV值预测"在这个环境中感到不适的人口比例"的指标。公式是：
诶，即使PMV=0（完全中立），PPD也不会是0？
没错，PMV=0时PPD=5%。也就是说再完美的环境，也至少有5%的人会感到不适——这反映了人体温感的个体差异。在实际工程中，目标通常设定为PPD<10%（对应PMV±0.5）。
| PMV | PPD (%) | 主观评价 | ISO 7730类别 |
|---|---|---|---|
| 0 | 5 | 中立（舒适） | A（PPD<6%） |
| ±0.5 | 10 | 略微温暖/略微凉爽 | B（PPD<10%） |
| ±0.7 | 15 | — | C（PPD<15%） |
| ±1.0 | 26 | 温暖/凉爽 | 标准外 |
| ±2.0 | 77 | 炎热/寒冷 | 标准外 |
ISO 7730与ASHRAE 55
关于PMV/PPD有什么国际标准吗？
主要有两个标准：
- ISO 7730（国际标准）：规定了PMV/PPD的计算方法和舒适等级A/B/C的定义。在欧洲空调设计中是标准
- ASHRAE标准55（美国标准）：除了PMV法，还包括自然通风建筑的自适应舒适模型。在国际项目中应用更广泛
日本有SHASE-S102（日本空调卫生工程学会标准），以ISO 7730为基础。实际工程中通常这样说："设计目标是ISO 7730的B级以上"。
Fanger与1396人的实验
Ole Fanger（1934~2006）是丹麦工科大学教授，在1967年的博士论文中发表了PMV模型的原型。他的实验规模在当时是罕见的，使用了精密控制温度、湿度和风速的气候室，让被试者在各种着装和活动条件下用7级量表评价自己的温感。他把统计结果与人体热平衡方程结合，这个思路在当时是革命性的。1984年，这个模型被标准化为ISO 7730。Fanger被誉为"舒适工程之父"，他的公式50多年来一直是全球空调设计的基础。
数值求解与CFD耦合
CFD-PMV耦合的基本策略
PMV本身是个代数式，为什么必须与CFD耦合呢？只要测量房间的温度和湿度，不就能算出来了吗？
如果只用一点的温度和湿度来算，确实可以。但实际的室内空间里温度和气流分布完全不同。比如空调出风口下面风速0.5 m/s很凉爽，但3米外的窗边因为有日光辐射会很热。要显示这种空间分布，CFD是必不可少的。
原来如此。所以就是用CFD算出温度场和速度场，然后代入PMV公式？
正是。耦合策略主要有两种：
- 后处理法（单向）：先用CFD求稳态解，得到各网格点的$T_a$、$v_a$、湿度。从湿度算出$p_a$，从辐射模型得到$\bar{T}_r$，设计条件给出$M$和$I_{cl}$，然后计算PMV。这是最常用的方法
- 耦合法（双向）：在CFD网格中放入人体假体模型，把人体的发热和发汗作为热源项。精度更高，但建模工作量大
实际上，后处理法大多数情况下已经足够。人体产生的热量相对于空调容量来说比例很小，对整体的温度场影响不大。
人体热收支的离散化
从CFD结果计算PMV时，具体怎么写代码呢？
在Ansys Fluent中用UDF（用户自定义函数），在STAR-CCM+中用Field Function，对每个网格单元执行下列步骤：
- 从CFD结果读取单元的$T_a$（空气温度）、$v_a$（速度向量的模）、$p_a$（水蒸气分压）
- 从辐射模型或表面对表面辐射计算$\bar{T}_r$（平均辐射温度）
- 从设计条件读入$M$、$W$、$I_{cl}$（通常为常数）
- 用Newton-Raphson法迭代求解$T_{cl}$（初始值35°C，收敛条件$10^{-3}$°C）
- 计算热负荷$L$
- 计算PMV = [0.303·exp(-0.036M) + 0.028]·$L$
- 计算PPD = 100 - 95·exp(-0.03353·PMV⁴ - 0.2179·PMV²)
计算量本身很小——问题的关键在于CFD的温度场、速度场、辐射场的精度。特别是辐射温度，对结果影响很大。
辐射温度的CFD计算
平均辐射温度$\bar{T}_r$怎样在CFD中计算？普通CFD不是只给出温度和速度吗？
这是个很关键的问题。$\bar{T}_r$是周围所有表面红外辐射的合成"等效温度"，原理上用黑球温度计测量。CFD中有几个办法：
- 表面对表面（S2S）辐射模型：计算view factor（视角系数），求解表面间的辐射热交换。精度高但计算量大
- 离散纵坐标（DO）模型：把辐射强度按方向离散化求解。参与媒介（烟雾、气体）存在时需要
- 简化公式：用墙面温度和view factor的加权平均近似$\bar{T}_r$。计算快，壁面温度差异小时有效
简化公式的形式是：
其中$F_{p-i}$是评估点到表面$i$的view factor，$T_{s,i}$是表面$i$的绝对温度。
网格策略与居住区的定义
为了评估PMV，CFD网格有什么特别的要求吗？
重点是居住区（Occupied Zone）的定义和网格分辨率。ISO 7730规定居住区是"距地面0.1~1.8m（站立）或1.1m（坐位），距外墙0.5m以上"的范围。PMV评估只在这个区域进行。
网格的关键点：
- 居住区内：网格间距50~100mm，均匀分布。要准确捕捉温度和速度的梯度
- 空调出风口周边：10~20mm以下的细网格。要分辨喷流的初始行为
- 壁面附近：y⁺ ≈ 1（近壁面低Reynolds处理）或y⁺ = 30~300（壁面函数）。直接影响辐射温度精度
- 天花板以上、地板以下等非居住区：可以用粗网格降低计算成本
PMV/PPD温热舒适性的实践应用
分析工作流程
实际工作中怎样开展PMV仿真呢？
典型的工作流程是这样的：
- 导入建筑CAD模型：从BIM（Revit/IFC）或CAD数据导入CFD预处理器。保留墙、窗、天花板、地板各个表面
- 配置空调设备：设定出风口类型（喷淋、直线等）、出风温度、风量、喷出角度等
- 设置边界条件：外墙温度（或传热系数+外气温），窗户日射量，内部热源（照明、设备、人体）
- CFD计算：用RANS（k-ε或SST k-ω）求稳态解。对室内气流，SST k-ω精度通常更好
- PMV/PPD后处理：在居住区水平截面（FL+0.6m、FL+1.1m、FL+1.7m）输出PMV等高线
- 识别不适区并优化：找出PMV超出±0.5的区域，调整出风位置或风量，重新计算
边界条件的设定
设置边界条件时，有什么容易出错的地方？现场常见的误区？
有3个容易踩的坑：
- 忽视窗户日射：有日射的窗附近，$\bar{T}_r$会比$T_a$高5~10°C。夏季南向窗一定要考虑日射取得，不然会低估热负荷
- 遗漏人体发热：每个人约100W的热发出（显热60W+潜热40W）。30人的办公室就有3kW，忽略的话室温会算低几摄氏度
- 设备发热过估：不要直接用设备设计负荷，要考虑实际使用率。一台台式电脑实际发热60~80W左右，不是额定功耗
验证与确认
计算结果对不对，怎样检验？
验证分两个层面：
- CFD层面：拿实测的室温和流速来对比。如果是竣工后的建筑，用温湿度数据记录仪在多个位置测量，CFD预测值的误差应在±1°C以内。流速用风速仪测量，误差应在±0.05 m/s以内
- PMV层面：在同样的测点处计算PMV，和使用者的问卷调查结果比较。但是个体差异很大，需要N≥20的样本量才能有效。实际上，如果温度和风速都对了，PMV一般也就对了
网格收敛性也要验证：网格数加倍后，PMV值的变化应小于0.1。
案例研究：办公室空调设计
能给个具体的例子吗？办公室的PMV分析通常是什么样的？
典型场景：10m × 15m × 2.8m的办公楼层，40人坐着工作，4台天花板嵌入式空调机。夏季设计工况的分析：
- 设计参数：$M$ = 1.2 met（PC工作）、$I_{cl}$ = 0.5 clo（清凉装）
- 空调设置：出风温度16°C、总风量3600 m³/h
- 外部条件：外气35°C、南向窗日射300 W/m²
计算通常会发现：窗边1.5m范围内PMV > +0.5（太热），空调下方PMV < -0.5（吹风感冷）。改进方案：
- 在窗边增加专用周边空调
- 窗户百叶遮光率设定60%以上
- 空调出风角度调整为水平，减轻吹风感
一般要迭代2~3轮，反复调整空调布局，最后才能落实到最优方案。这个过程就是PMV仿真的最大价值——在纸面上发现和解决舒适性问题。
PMV/PPD温热舒适性的软件对比
商用工具对比
做PMV仿真有什么软件可以用？
大概分两类："通用CFD+PMV后处理"和"建筑环境专用工具"：
| 工具 | 开发商 | PMV支持 | 特点 |
|---|---|---|---|
| Ansys Fluent | Ansys Inc. | UDF实现 | 通用CFD龙头。辐射模型完整。PMV需用户自己编写UDF |
| STAR-CCM+ | Siemens | Field Function | 多面体网格、自动化强。有HVAC专用模板 |
| COMSOL Multiphysics | COMSOL AB | 热传递+自定义 | 多物理耦合方便。教学和研究用 |
| IDA ICE | EQUA Simulation AB | 原生支持 | 建筑环保模拟专用。Comfort Analyzer直接输出PMV分布 |
| EnergyPlus + OpenStudio | 美国DOE | 原生支持 | 全年能耗模拟。可直接输出区域PMV |
| DesignBuilder | DesignBuilder Software | 原生支持 | EnergyPlus的可视化前端。有CFD模块可计算室内气流 |
通用CFD和专用工具，哪个更好？
看用途。设计初期的全年能耗评估用EnergyPlus或IDA ICE最快，几分钟就能算出全年每小时的PMV。但如果需要详细的空间分布可视化（"这个座位太冷"），就需要用Fluent或STAR-CCM+这样的CFD。实务中常见的做法是：先用EnergyPlus定空调容量，再用CFD细化出风设计。
开源选项
有免费的工具吗？学生预算有限…
有几个选择：
- OpenFOAM + PMV库：室内CFD可加自定义Function Object实现PMV。GitHub上有公开实现
- CBE热舒适工具（加州伯克利大学）：网页版，可交互式计算单点PMV。没有CFD耦合，但快速检验很方便
- pythermalcomfort：Python包，符合ASHRAE 55/ISO 7730。能读CFD结果CSV文件批量计算PMV
- JOS-3（产业技术总合研究所）：日本开发的人体热模型。比PMV更细致（部位级），适应日本人体型
黑球温度计的原理
平均辐射温度$\bar{T}_r$的实测仪器是"黑球温度计"（Globe Thermometer）。它是一个直径150mm的薄铜球，涂成黑色，中间放温度计。球面吸收辐射（近似黑体），经过对流和辐射平衡后达到平衡温度。球温$T_g$与空气温$T_a$、风速$v_a$的关系式为：$\bar{T}_r = [(T_g+273)^4 + 2.5 \times 10^8 \cdot v_a^{0.6}(T_g - T_a)]^{1/4} - 273$。缺点是响应慢（平衡要15~20分钟），现在逐渐被椭球形辐射计和红外传感器阵列取代。
PMV/PPD温热舒适性的先进研究
自适应舒适性模型（自适应舒适）
我听说PMV模型有局限。是什么意思？
PMV最大的局限是基于稳态假设+空调建筑。对自然通风建筑（可以开窗的办公室）来说，人会根据外界气温变化来调整着装和窗户开度——这就是"适应"。例如，外气30°C那天，室温28°C就感到舒适，但外气20°C时室温28°C就觉得太热了。这个适应效应PMV无法表现。
为此，ASHRAE 55引入了自适应舒适模型，用月平均外气温的函数来定义舒适温度：
其中$T_{out,mean}$是月平均室外温度。这个公式只适用于可自主调节环境（开窗、调衣服）的自然通风建筑。
也就是说，空调楼用PMV，自然通风建筑用自适应模型，这样分开用？
正是。最近的"混合模式建筑"（既有空调又能开窗）会根据季节或运行模式切换两个模型。LEED认证和WELL认证的取得需要这种分层方法的文件说明。
数字孪生与实时PMV
"数字孪生"与PMV怎样结合？我听说能实时控制空调。
把建筑BIM模型与IoT传感器数据（温度、湿度、CO2、在室人数等）实时关联起来，用提前训练的代理模型（神经网络）秒级预测PMV分布。CFD用于生成训练数据，运行期间用更快的机器学习模型替代。大林组和NTT Facilities在试点。
基于PMV预测的空调控制比传统的温度点控制能降低15~25%的能耗。
局部不适感的评估
PMV在±0.5以内，但还是有人说"脚冷"、"脖子吹得难受"。这怎么办？
这是局部不适（局部热不适）问题，PMV/PPD单独无法评估。ISO 7730建议追加4个局部因素的评价：
- 吹风感（Draft Risk）：局部气流冷却。用公式DR(%) = (34-$T_a$)(v-0.05)⁰·⁶²(0.37·$v$·Tu+3.14) 估算，其中Tu是气流湍流度（%）
- 竖向温度差：头部和脚部温差≥3°C时不适感明显
- 地面温度：推荐19~29°C范围（特别是地暖时要注意）
- 辐射不对称：窗户或散热器的单侧辐射。温暖天花的辐射不对称应≤5°C
从CFD结果计算吹风风险相对容易，Fluent/STAR-CCM+的PMV UDF中常会一并输出吹风等级。
PMV/PPD温热舒适性的故障排除
常见错误与对策
初学者在PMV仿真中容易踩什么坑？
现场常见的故障：
1. PMV在全域都是极端值（如全部>+1.5或<-2）
症状：PMV值在所有网格上都非常高或非常低
原因：$M$或$I_{cl}$单位设错。代谢率用了met而不是W/m²（相差58.15倍），或着装量用了m²K/W而不是clo（相差0.155倍）
对策：在单点用手算或pythermalcomfort验证，与UDF/Field Function输出对比，确保单位一致
2. $T_{cl}$迭代计算不收敛
症状：PMV计算子程序死循环
原因：$T_{cl}$初值太极端（如<0°C），或收敛条件过严（如$10^{-4}$°C）。辐射项的4次方在数值上可能不稳定
对策：初值固定为35°C，收敛判定放宽到$10^{-2}$°C，反复次数上限150次，如仍未收敛就用最后的值
3. 辐射温度评估错误
症状：窗边PMV异常高，或辐射效应完全看不出
原因：S2S辐射模型的view factor精度不足，或根本没启用辐射模型
对策：Fluent中启用Surface-to-Surface或DO模型。检查墙面射出率设置（一般内装0.9，铝型0.1~0.3）。窗玻璃对红外近似不透明（射出率≈0.84）
PMV值与实测不符的情况
实测问卷反映"太热"，但仿真PMV=0左右。什么原因？
通常查三个方面：
- 日射效应：CFD中日射常作为"壁面吸收热"来模型化，但透射光直接加热人体的影响容易漏掉。有日射时$\bar{T}_r$能升高10°C以上
- 代谢量：是否假设了"安静坐着"（M=1.0 met）？实际PC工作约1.1~1.2 met，讨论会1.4 met
- 心理因素：PMV是物理模型，无法表现"房间太小""噪音大""灯光冷"等非热因素。这些会强化"太热"的主观感受
先验证1和2，如果还有差异，就要在报告中声明PMV模型的适用范围和局限。
老师，PMV仿真的全体过程我大致理解了。核心还是CFD的温度场、风速场、辐射场的精度？
完全同意。PMV的公式本身是"成熟技术"，50年来没有大改变，计算本身很简单。仿真的质量取决于CFD的精度和边界条件的合理性——这是工程师的功力。我建议从单点的手算开始，然后和CFD结果对标，这样能深刻理解整个流程。