HVAC空调CFD
HVAC空调CFD的理论基础
概要
老师！HVAC空调的CFD解析在什么场景下使用？
针对办公室、商业设施、数据中心、医院等室内环境，通过CFD预测气流分布、温度分布、换气效率、舒适性（PMV/PPD）。从空调设备的设计阶段性能验证到既有建筑的改造规划，应用范围广泛。
支配方程
请讲解室内气流的支配方程。
低速流（马赫数 < 0.3）的不可压缩Navier-Stokes方程 + 能量方程是基础。浮力效应通过Boussinesq近似引入。
此项作为Navier-Stokes方程的源项加入。$\beta$ 是体膨胀系数（空气情况下 $\beta \approx 1/T_0$ [1/K]）。
Boussinesq近似是在温度差较小时的近似，对吧。那室内空调的温度差可以允许多大？
条件是 $\beta \Delta T \ll 1$，所以空气情况下 $\Delta T < 30$ K左右就没问题。室内空调通常在10 K以内，所以无需担心。
舒适性指标
PMV（预测平均评价）在CFD中怎样计算？
PMV是ISO 7730定义的热舒适性指标，由6个参数计算。
| 参数 | 记号 | 从CFD的获取方法 |
|---|---|---|
| 代谢率 | M | 从活动水平设定（办公作业: 1.2 met） |
| 着衣量 | I_cl | 从季节、用途设定（夏装: 0.5 clo） |
| 空气温度 | T_a | 直接从CFD结果获取 |
| 平均辐射温度 | T_r | 从CFD + 辐射模型计算 |
| 气流速度 | v_a | 直接从CFD结果获取 |
| 相对湿度 | RH | 通过Species Transport计算，或假设均匀 |
既然需要平均辐射温度，就必须打开辐射模型，对吧。
是的。使用S2S（Surface-to-Surface）模型或DO（Discrete Ordinates）模型计算壁面间的辐射传热，求得各点的T_r。
换气效率的评价指标
换气效率有哪些指标？
代表性指标汇总如下。
| 指标 | 定义 | 含义 |
|---|---|---|
| ACH (Air Changes per Hour) | Q / V_room | 换气次数 |
| AE (Air Exchange Efficiency) | τ_n / (2τ_mean) | 新鲜空气的分配效率 |
| SVE (Contaminant Removal Effectiveness) | C_e / C_mean | 污染物的除去效率 |
| Local Mean Age of Air | τ(x) | 各点空气的滞留时间 |
Local Mean Age of Air是如何在CFD中计算的？
添加标量输运方程，计算空气的平均滞留时间。在Fluent中使用UDS（User Defined Scalar）。
源项是 $\rho$（均匀的），所以空气停留在室内越久，$\tau$ 就越大，对吧。
空调CFD的源流——1970年代的建筑能源问题所产生的室内气流解析
室内气流的CFD解析（Room CFD）开始广泛应用是在1970年代石油危机之后，省能建筑设计需求上升的时代。初期以区域法（Zone Model）的简单混合空气模型为主流，但随着计算机普及，Navier-Stokes解法开始被应用于室内气流。Nielsen(1974)首次发表了使用k-ε模型的室内气流CFD解析，为吹出口、吸入口的设计提供了科学依据。这是现今CFD-HVAC解析的起点。40年后的现代，数百万网格的非定常LES解析成为标准，实现了对换气效率（Ventilation Efficiency）和在室人员热舒适性（PMV指标）的实时预测。
HVAC空调CFD的数值计算方法
数值方法的详细
请讲解HVAC CFD的具体实现。
湍流模型的选择
室内气流推荐使用SST k-omega 或 RNG k-epsilon。特别是在天井吹出的混合换气中，壁面喷流的行为很重要，所以SST k-omega更受欢迎。
| 换气方式 | 推荐湍流模型 | 原因 |
|---|---|---|
| 混合换气（天井吹出） | SST k-omega | 捕捉壁面喷流的分离 |
| 置换换气（地板吹出） | RNG k-epsilon | 层流-湍流过渡域的处理 |
| 个人换气 | SST k-omega | 低速喷流的精度 |
| 自然换气 | SST k-omega + Boussinesq | 适应浮力驱动流 |
吹出口的建模
空调吹出口形状复杂，是否所有内部都要网格化？
对扩散器的内部形状进行全部网格化效率很低。改为使用Simplified Diffuser Model（SDM），直接在吹出口面施加速度剖面（风速、角度、湍流量）。
常用扩散器类型的设置:
| 扩散器 | 吹出角度 | 有效面积比 | 湍流强度 |
|---|---|---|---|
| 4向天井盒 | 水平～15°下 | 50～70% | 10～15% |
| 风口 | 放射状、45° | 60～80% | 15～20% |
| 线性扩散器 | 水平 | 70～90% | 10% |
| 百叶通风口 | 可变（0～60°） | 80～95% | 5～10% |
| 地板吹出口 | 竖直向上 | 20～40% | 20～30% |
什么是有效面积比？
是扩散器喉部面积与实际吹出面积的比值。百叶和翅片遮挡的部分不算。应根据制造商目录中的Coanda效应到达距离与CFD结果对照验证合理性。
网格策略
室内空间的网格数量有什么参考？
一般的办公室（10m x 15m x 3m）的参考网格数是200万～1000万个。
| 区域 | 单元大小 |
|---|---|
| 吹出口/吸入口周围 | 10～30 mm |
| 人体、家具周围 | 20～50 mm |
| 居住域（FL+0.1m～FL+1.8m） | 30～80 mm |
| 天井附近（喷流域） | 20～50 mm |
| 其他（空间中央） | 80～200 mm |
辐射模型的设置
应该使用哪个辐射模型？
室内环境推荐使用S2S（Surface-to-Surface）模型。预先计算壁面间的View Factor进行辐射传热评估。DO模型也可用，但对于仅有不透明壁面的室内环境，S2S已经足够。
壁面放射率的设置:
| 表面 | 放射率 |
|---|---|
| 混凝土墙 | 0.90～0.95 |
| 玻璃窗 | 0.84～0.90 |
| 金属（涂装） | 0.85～0.95 |
| 金属（未涂装） | 0.05～0.20 |
| 人体表面 | 0.95～0.97 |
玻璃窗的日射怎样处理？
启用Solar Load Model，从太阳位置（纬度、经度、日期、时刻）和窗户的SHGC（Solar Heat Gain Coefficient）计算日射负荷。Fluent标配此功能。
HVAC CFD的湍流模型选择——低雷诺数环境下标准k-ε的局限
室内气流（HVAC）的CFD解析是低速（Re=10³～10⁵）且浮力影响大的特殊环境，不同于外部流体解析。标准k-ε模型在此低雷诺数域易过高估计浮力驱动流，导致温度分层（Thermal Stratification）预测精度下降。更恰当的选择有①RNG k-ε（均衡性好）②Low-Re k-ε模型（Launder-Sharma等）③LES（精度最高、成本高）。日本建筑学会的验证报告表明，在吹出流速0.5m/s以下的自然对流主导域，RNG k-ε与LES结果偏差在10%以内，但在喷流型吹出（高速域）则出现25%以上的偏离。应根据分析目标的支配流动模式选择合适模型才能把握精度。
HVAC空调CFD的实务应用
实践指南
请讲解HVAC CFD的实务案例。
案例1: 数据中心的热管理
数据中心的服务器机架冷却是HVAC CFD最重要的应用领域之一。评估热通道/冷通道配置的温度分布，确定局部热点。
典型的边界条件:
- CRAC（Computer Room Air Conditioning）: 吹出温度15℃、设计风量
- 服务器机架: 发热量5～30 kW/机架
- 地板下防烟通道: 活动地板（瓷砖开口率25%）
听说有个RCI（机架冷却指标）的指标。
RCI是ASHRAE TC 9.9定义的指标。
$T_{in}$ 是机架吸入温度，$T_{rec,max}$ 是ASHRAE推荐上限（27℃）。RCI > 95%是目标。
案例2: 医院的感染控制
COVID-19之后，医院的通风设计备受关注。
是的。飞沫、气溶胶的扩散预测变得重要。用DPM追踪咳嗽喷嚏的飞沫轨迹，评估感染风险。
咳嗽的建模参数:
- 初速: 10～20 m/s
- 飞沫直径: 1～500 um（Rosin-Rammler分布）
- 飞沫数: 约3000个/咳
- 温度: 34℃
- 相对湿度: 95%
飞沫的蒸发也要建模吗？
是的。大飞沫因重力下落，但蒸发后变成飞沫核（droplet nuclei < 5 um）就能长时间悬浮在空中。需将Evaporation model与DPM耦合。
案例3: 大空间的温度分层评估
在中庭、工厂等大空间，会形成温度分层。暖空气积累在天井附近，与居住域的温度差成为问题。
对于置换换气中的温度分层高度预测，Archimedes数是重要参数。
Ar数越大，浮力越强，温度分层越稳定，对吧。
常见失败与对策
| 失败模式 | 原因 | 对策 |
|---|---|---|
| 室温整体过高 | 壁面隔热条件不当 | 设置外墙的热透过率（U值） |
| 吹出喷流迅速消散 | 数值扩散 | 使用二阶及以上的格式 + 网格细分 |
| PMV一致趋于零 | 辐射模型未启用 | 启动S2S模型计算T_r |
| 气流模式与实测不符 | 吹出口建模不当 | 用制造商目录数据验证CFD结果 |
医院空调CFD——手术室的气流设计使感染风险接近零
医院手术室是对空调设计要求最严格的空间之一。为防止细菌、微粒落入手术区域，需维持天井层流下吹（风速0.25～0.45m/s），同时控制人体和设备产生的对流。CFD（HVAC解析）会在空间内配置外科医生、护士的人体模型，预测手术灯的热湍流对层流的影响。实际设计中，手术室污染指标"CFD-OI（Operating room Infection Index）"目标为0.01以下，通过优化空调吹出口的位置和面积，既能维持微粒计数（ISO Class 5相当），还能实现省能的案例已由国内医疗设备厂商报告。
HVAC空调CFD的软件比较
商用工具比较
哪些工具适合HVAC CFD？
既有建筑环境专用CFD工具，也有汎用CFD工具。
| 工具 | 开发商 | 特点 | 目标用户 |
|---|---|---|---|
| Ansys Fluent | Ansys | 太阳负荷、PMV/PPD、DPM | CFD工程师 |
| STAR-CCM+ | Siemens | 自动网格、人体热舒适模型 | CFD工程师 |
| FlowDesigner | Advanced Knowledge Lab | 空调设备设计UI | 空调设计者 |
| Stream (scSTREAM) | MSC Software | 建筑CFD、正交网格 | 空调设计者 |
| SimScale | SimScale GmbH | 云基CFD、有免费额度 | 全体工程师 |
| Autodesk CFD | Autodesk | 与Revit/Inventor联动 | 建筑设计者 |
| OpenFOAM | 开源软件 | buoyantSimpleFoam | 研究者、高级用户 |
| 6SigmaRoom | Future Facilities | 数据中心专用 | DC设施管理者 |
6SigmaRoom是数据中心专用的吗？
是的。服务器机架、瓷砖、CRAC单元的库很齐全，针对数据中心的工作流做了优化。不需要CFD知识的DC工程师也能使用，这是优势。
用途别推荐工具
| 用途 | 推荐 | 原因 |
|---|---|---|
| 办公室舒适性评估 | FlowDesigner / Fluent | PMV/PPD计算、太阳负荷 |
| 数据中心热管理 | 6SigmaRoom / Fluent | DC特有库 / 高精度 |
| 医院感染控制 | Fluent / STAR-CCM+ | DPM飞沫追踪功能完善 |
| 与BIM联动 | Autodesk CFD / SimScale | Revit/IFC数据直接读入 |
| 大规模设施（站、机场） | Fluent / STAR-CCM+ | HPC对应、大规模网格 |
| 自然换气研究 | OpenFOAM | LES、长时间非定常 |
与BIM的联动很重要吗？
建筑设计采用BIM（Revit、ArchiCAD）作为标准，能直接导入BIM模型到CFD会大幅提高工作流效率。但BIM模型包含CFD不需要的细节（布线、金物），简化形状的工作仍然存在。
OpenFOAM进行室内环境解析
用OpenFOAM解析室内环境时用什么求解器？
使用buoyantSimpleFoam（定常、Boussinesq浮力）或buoyantPimpleFoam（非定常）。辐射模型可用viewFactor（S2S相当）。
设置要点:
- turbulenceProperties: kOmegaSST
- RASModel: kOmegaSST
- radiationModel: viewFactor
- g: (0 0 -9.81)
- Boussinesq: 在thermophysicalProperties中设置beta值
HVAC专用CFD工具——IDA与DesignBuilder及汎用CFD的区分
建筑空调的CFD解析工具分"HVAC专用"和"汎用CFD"两个体系。IDA ICE（EQUA）、DesignBuilder（EnergyPlus后端）专门针对建筑能量模拟（BES），可在一年8760小时内计算多个房间、多个区域。但各房间内的详细气流分布用简单的区域模型近似。需要详细室内气流分布（吹风感、污染物浓度分布）时使用汎用CFD（Fluent、StarCCM+、PHOENICS）。近年SimScale（云CFD）为HVAC推出了简易室内气流解析模板，降低了无专业CFD工程师情况下建筑设计者进行概略解析的门槛。
HVAC空调CFD的先端研究
先端话题与研究动向
请讲解HVAC CFD的最新趋势。
1. LES进行室内气流的非定常评估
RANS只能获得时间平均的气流模式，但LES能评估瞬间气流波动（吹风感）。对ISO 7730的吹风风险（DR）预测精度提高。
Tu是湍流强度，LES可以直接计算。
是的。RANS的Tu从定义上是时间平均值，但LES能从瞬间速度波动计算局部Tu，所以DR的空间分布更准确。
2. 人体温热模型的高度化
传统PMV是定常状态的舒适指标，最近使用UCB（UC Berkeley）Comfort Model和Fiala model等，动态计算人体各部位温度的多节点模型。
有研究用人体模型而不是圆柱来建模吧。
CFD人体模型（Computational Thermal Manikin）将头部、躯干、四肢作为独立发热面，同时考虑衣服的隔热效应。STAR-CCM+标配Thermal Comfort model。
3. 基于CO2浓度的DCV（Demand-Controlled Ventilation）
用CO2作为示踪气体，根据在室人数进行换气量控制（DCV）的CFD模拟。用Species Transport方程追踪CO2扩散。
人均CO2排放量大约是多少？
静坐时约0.005 L/s（18 L/hr），办公作业约0.006 L/s，轻作业约0.009 L/s。ASHRAE 62.1推荐室内CO2保持在1000 ppm以下。
4. 数字孪生与BMS联动
建筑管理系统（BMS）的传感器数据与CFD模型耦合的数字孪生逐渐普及。
不是实时运行CFD，而是用ROM和代理模型，对吧。
是的。完整CFD单个案例需要数小时，不能用于实时控制。提前进行大量CFD，构建ROM并整合到BMS中。
数字孪生HVAC——传感器融合与CFD-ROM的实时控制
最先端的智能建筑管理采用"数字孪生HVAC"系统，将室内CO₂浓度、温湿度传感器的数据与CFD-ROM（正交本征分解的低次元模型）相结合。通过POD（Proper Orthogonal Decomposition）事先压缩CFD解空间，用实时传感器数据更新权重系数，可在完整CFD计算时间的1/1000内估计室内气流。新加坡大型购物中心实施的案例，通过该系统实时考虑在室人数和外气条件动态优化AHU（空气处理机组）风量，空调能耗削减了15%的成绩已有报道。
HVAC空调CFD的故障排除
故障排除
HVAC CFD有什么常见问题？
1. 室温整体与设计值不符
检查项:
- 确认壁面热条件。外墙是否使用了Convection BC（外气温度+外侧换热系数）
- 内部发热（照明: 10～20 W/m²、人体: 75 W/人、办公设备: 50～150 W/台）是否正确输入
- 空调吹出温度和风量是否与设计值一致
- 日射负荷（Solar Load）是否被考虑（有窗房间则为主要负荷源）
2. 吹出喷流到达距离很短
天井吹出的喷流很快就消散了。
对策:
- 确认离散格式（First Order Upwind会因数值扩散导致喷流扩散）
- 对吹出口周围网格进行细分（吹出口外10倍距离以上逐步粗化）
- 确认吹出风速和湍流强度是否正确（使用扩散器目录值）
- 为再现Coanda效应（喷流贴附天井面），天井面附近网格要足够细
3. Boussinesq近似下残差振荡
对策:
- 将压力插值改为PRESTO!或Body Force Weighted
- 设置Under-Relaxation Factor中的BodyForce为0.8
- 确认Operating Density（Boussinesq参考密度）是否为合适值
- 确认Operating Pressure和Reference Pressure的设置
4. 辐射模型导致计算时间极长
原因: S2S模型的View Factor计算是瓶颈。面的数量太多。
对策:
- 启用S2S的Face Clustering（面分组化）来减少面数
- 调整Cluster数量（精度和成本的权衡）
- 在辐射影响小的区域（壁面温差3 K以内）忽略辐射
5. 实测的风速分布不吻合
检查清单:
| 确认项 | 常见问题 |
|---|---|
| 吹出口有效面积比 | 未使用目录值 |
| 家具、人体配置 | 作为气流障碍物影响大 |
| 门的开闭状态 | 邻室压力差使气流改变 |
| 外部风压 | 窗、外墙开口的漏气 |
| 吹出温度漂移 | 空调机的控制滞后 |
家具配置也有这么大影响吗？隔板这样的也要建模？
隔板（高度1.2m以上）对气流模式影响大，必须建模。书桌、柜子高度0.7m左右的可以简化成直方体。
空调CFD的"计算合格、现场不快"——吹风感预测的陷阱
空调CFD显示温度分布符合设计值，但现场出现"吹风冷"投诉的情况屡见不鲜。问题在"吹风感"的预测精度上。ISO 7730将人体不快的吹风定义为DR（Draft Rate）指标，在风速0.2m/s以上且湍流强度高的区域易发生。CFD的RANS解析虽然时间平均风速准确，但往往低估局部湍流波动（u'），导致DR预测过于乐观。解决方案有两个：①对时间平均风速加上乱流强度×0.37的修正使用ISO公式、②通过URANS或LES直接计算瞬间风速波动。但实务中能用LES的环境仍然有限。
相关话题
价值
更详细
错误