说起来有点讽刺搞通风空调的人天天被甲方问“你这新风量够不够”但真正到了现场一测风口轰轰响工位上的人还是觉得闷二氧化碳浓度该超标还是超标。问题出在哪绝大多数项目在做气流组织设计时只算了“送进去多少风”没算“风到了哪里、花了多长时间才到”。这时候光靠手算换气次数是完全不够的必须上CFD。我手头正好做完一个办公室气流组织评估项目用Fluent把房间里的速度场、温度场和空气龄Mean Age of Air, MAA全部跑了一遍。这套“流场分析空气龄仿真”的组合基本是复杂室内环境通风评估的标配思路。这篇文章不打算写教科书就直接把我踩过的坑、试过的设置、能直接套用的UDF和边界条件都放出来给做暖通设计、洁净室评估、室内环境研究的同行一个能落地的参考。1. 室内流场分析为什么离不开空气龄1.1 速度场只能告诉你“风从哪里来”空气龄才能告诉你“空气走了多久”很多刚接触CFD的同学会有个误区以为把速度云图跑出来看哪里有风、哪里没风就够了。但实际上速度场解决的是“空气怎么运动”的问题它回答不了“新鲜空气到达某个位置的效率有多高”。说得直白一点一个房间的换气次数如果是10次/小时理论名义时间常数τn6分钟听起来很美好。但真实房间存在短路流、回流区、死区有些区域的空气可能在里面待了40分钟还没被换掉。速度云图上看那些角落可能只是低速区你无法从速度大小直接判断空气的“新鲜程度”。空气龄的定义很朴素从送风口进入房间的空气微团到达房间内某个点所经历的时间。这个值越小说明该点被新鲜空气“洗”得越勤值越大说明该点空气陈旧污染物容易积聚。在Fluent里我们通常用用户自定义标量方程UDS来求解它后面我会详细讲UDS设置。流场分析和空气龄仿真必须“搭伙”的原因也在这里空气龄本质上是流场输运的结果流场不收敛空气龄必错流场收敛了但流型本身有问题比如送排风短路空气龄云图会非常直观地暴露问题。它是连接“送风量”和“实际通风效果”之间的桥梁。1.2 空气龄仿真的两条技术路线UDS标量法与示踪气体法做空气龄仿真业内基本有两条路一种是稳态标量输运法另一种是瞬态示踪气体法。稳态UDS标量法的原理是在Fluent中增加一个额外的标量输运方程方程的非稳态项设为零对流项照常扩散项使用有效扩散系数源项设为1代表单位体积内空气老化速率为1即每立方米空气每秒钟老1秒。入口边界给固定值0出口边界用法向梯度为零。这样解出来的标量场就是“空气龄”场单位是秒。这种方法计算量小工程上最常用尤其适合做方案比选。瞬态示踪气体法则更接近实验在送风口注入一定浓度的示踪气体比如SF6通过监测室内各点浓度随时间的变化曲线来推算空气龄。这个方法可以用Fluent的组分输运模型瞬态求解也可以用实验实测。它的好处是能同时分析污染物传播过程坏处是计算量成倍增加而且实验注入气体的操作门槛高。两种方法怎么选我一般这么判断对比项UDS标量法空气龄示踪气体法浓度衰减求解方式稳态多算一个UDS方程瞬态组分输运/多组分模型典型用途快速评估气流组织、方案比选与实验对标、分析污染物动态传播边界条件入口0出口零梯度入口给示踪气体浓度曲线计算成本低高精度工程上足够更贴近实验但对网格和时间步长敏感如果你只是想知道“哪个区域通风差”UDS法完全够用如果你要复现实验数据或者做污染物扩散全过程分析就得老老实实上瞬态示踪法。1.3 空气龄衍生指标与评价口径光有一张空气龄云图还不够工程上不能指着云图跟甲方说“你看这里红”。我们需要把空气龄转换成几个能写进报告、能跟规范对标的评价指标名义时间常数τn V/QV是房间体积Q是送风量。它代表理想完全混合条件下的平均换气时间全室平均空气龄——把空间所有点的空气龄做体积平均换气效率ε τn / (2 × 平均空气龄)——这个值越接近1说明流动越接近“活塞式”理想推流一般混合通风的ε在0.3~0.5之间置换通风可以做到0.5~0.7局部空气龄——直接看房间内某个具体人员呼吸区的空气龄值。我在实际项目中习惯先看局部空气龄云图找“红区”再用全室平均空气龄做不同方案之间的量化对比。比如方案A比方案B平均空气龄降低了15%这比说“A看起来更通畅”要有说服力得多。2. 复杂室内环境的前处理几何、网格与边界条件2.1 几何模型简化到什么程度才算“够用”复杂室内环境建模的第一个拦路虎不是软件操作而是“模型画到多细”。把办公桌、椅子、电脑、人体、吊灯全部画出来网格量轻松破亿算个空气龄要等三天而且很多小特征还会引发收敛问题。我的原则是“抓大放小、以流场影响论英雄”必须保留送风口、排风口、房间围护结构、大的隔断、主要障碍物柱子、大柜子、发热量大的设备服务器机柜、打印机、大型设备、人员用简化体代替圆柱方体的组合就够可以简化文件柜上的小把手、圆桌变成正多边形都不影响必须剔除墙面上的插座、开关、小孔洞、踢脚线、螺丝钉之类这些特征对气流组织的影响可以忽略却会严重浪费网格。有人问我人体模型做不做我的建议是如果分析的是工位呼吸区空气龄或污染物暴露最好做简化人体模型模拟人体的热羽流对气流有显著影响如果只是分析整体通风效率人体可以不做用热流密度替代就行。CAD清理阶段还要注意“长细比”陷阱——只要有很薄的狭长结构网格划分时就会产生大量高扭曲率单元收敛困难。遇到这种直接压缩厚度或删除。2.2 网格划分的核心策略局部加密、边界层与“孤儿网格”室内CFD的网格划分核心就一句话在流动变化剧烈的区域加密在均匀流动区域放松。具体拆开来是三件事第一全局尺寸先给一个基准值。比如会议室模型约50㎡、层高3m基准体网格尺寸可以设在80~120mm送风口附近、排风口附近、热源表面加密到10~20mm。注意不要全局都加密那会让网格量爆炸。第二壁面边界层必须加。送风口、排风口、热源表面、人员模型表面这些涉及传热或有明显速度梯度的壁面建议加3~5层棱柱层首层厚度按照y要求估算。对于室内通风这类低雷诺数问题壁面附近不需要追求y≈1那是LES或者低雷诺数模型的要求用增强壁面处理时y控制在30~300都行但如果你用SST k-ω就得控制得严一些。第三复杂几何连接处最容易出现“孤儿网格”。这个词在中文社区里经常出现其实指的就是孤立网格区域——几何上本该连通的体区域因为缝隙、重叠、导入错误变成了互不相连的几个“孤岛”。具体表现是网格导入Fluent后提示存在多个volume regions或者计算时某些区域流速始终为0。排查方法很简单先执行Mesh → Check看是否有负体积再在Display里按region显示看是否存在没和主流域联通的部分。修复手段一般是回到几何软件里做布尔运算或gap修复或者在Fluent里用Mesh → Repair偶尔有效。2.3 边界条件的合理设定送风口、排风口、热源怎么给边界条件决定流场像不像话。室内通风仿真最常见的一批问题几乎都是送风口边界给错的。送风口最稳妥的是用速度入口velocity-inlet直接给送风风速和送风温度湍流强度按5%~10%估算水力直径按风口当量直径算。如果是散流器这种出流方向复杂的风口直接建模太费网格工程上常用“等效风口法”把散流器简化成一个平面这个平面给法向速度加切向偏转速度或者在风口后方加一个多孔介质区来等效阻力。这也就是为什么网上很多人问“多孔介质fluent参数设定”怎么算——风口处的均流板、滤网、百叶都可以用多孔介质模型替代避免画大量小孔网格。多孔介质参数怎么定通常需要两个系数粘性阻力系数1/α和惯性阻力系数C2。如果厂家能给压降-流速测试数据直接拟合ΔP/L (μ/α)v (0.5ρC2)v²没有数据时可以参考经验值比如常见金属滤网1/α约1e6~1e8 m⁻²C2约10~100 m⁻¹但做了之后最好和设计手册上的阻力特性核对一下别差一个数量级。排风口我一般用压力出口pressure-outlet表压设为0回流总温设为房间温度。如果排风口直接对着室外也可以用outflow但要保证出口流动充分发展。热源人体、设备、照明、外窗传热。最省事的方法是把所有热源简化为壁面热流密度W/m²也可以做带体积热源比如服务器机柜内部热源体。人体热羽流在不做人体模型时可以近似为座椅位置附近的一个热流面做了人体模型的话皮肤表面给33℃等温或按人体显热散热量折算热流。3. Fluent中的流场与空气龄实操全流程3.1 模型与求解器设置别一上来就挑大梁Fluent版本不同菜单略有区别但核心逻辑一样。我先说求解器整体设置再说空气龄UDS的细节。室内通风流场属于低速不可压流动压力基求解器Pressure-Based是标准选择。速度入口到压力出口之间的流动稳态求解通常就能收敛但如果研究的是人员走动干扰、门窗开关后的瞬态响应或者送风波动才需要切瞬态。湍流模型标准k-ε在室内通风领域用了二十年但它对旋流、强浮力流场的预报能力偏弱。我个人的习惯是普通混合通风用Realizable k-ε配合增强壁面处理有明显的自然对流主导或置换通风时RNG k-ε或者SST k-ω更稳一些。置换通风的核心特征是热羽流驱动浮力项必须正确开启开启能量方程并在Operating Conditions里设置重力加速度密度模型如果温差不大可以用Boussinesq近似密度只随温度线性变化温差超过40℃就直接用理想气体不可压密度。压力速度耦合用SIMPLEC收敛速度比SIMPLE快动量、能量、k、ε离散格式全部用二阶迎风First Order解析解精度差容易出现假扩散空气龄这种输运量对假扩散尤其敏感。浮力驱动的流场收敛残差经常卡在1e-3下不去。这时候不要硬降残差标准重点看监测点的速度和温度曲线是否平稳、出入口质量流量是否守恒。3.2 空气龄UDS设置全记录这一步是全文的重头戏我按菜单顺序把操作写清楚。第一步启用UDS方程。路径是Define → User-Defined → Scalars或者在Models面板里把User-Defined Scalars的个数设为1。这个UDS标量没有量纲但我们心里清楚它代表“秒”。第二步设置边界条件。入口边界在送风口边界条件的UDS标签页里把标量值设为0入口空气龄为0这是物理上最合理的设定排风口保持默认的零通量即法向梯度为零出口边界的标量会由流动自然带走。第三步也是最容易踩坑的一步扩散系数和源项必须通过UDF定义因为在Fluent的面板里UDS的扩散系数默认是0源项默认也是0你不定义算出来的空气龄就是纯对流保形输运数值完全不对。UDF写法我直接贴出来#include udf.h /* 空气龄UDS的源项单位体积内空气老化速率为1即每秒老1秒 */ DEFINE_SOURCE(age_source, c, t, dS, eqn) { dS[eqn] 0.0; return 1.0; } /* 空气龄UDS的扩散系数分子扩散湍流扩散 */ DEFINE_DIFFUSIVITY(age_diffusivity, c, t, i) { real mu_t C_MU_T(c, t); /* 湍流粘度 */ real rho C_R(c, t); real sc_t 0.7; /* 湍流施密特数经验值 */ real d_eff 2.0e-5 mu_t / (rho * sc_t); return d_eff; }扩散系数为什么这么写空气龄标量的输运包括分子扩散和湍流扩散两部分。分子扩散系数对于空气中的示踪气体大概是2e-5 m²/s量级但实际房间里湍流扩散远大于分子扩散用湍流粘度除以湍流施密特数Sc_t在0.7~1.0之间取值来估算。这个有效扩散系数决定了空气龄在空间中的“抹平”程度取太小会让空气龄分布呈清晰的“细丝”状取太大会让空气龄分布过于均匀、失去局部差异。一般0.7是稳妥值。第四步加载UDF。这就顺便回答一个很多人问的高频问题“Fluent里UDF文件在哪里编辑”UDF本质上就是一个C语言文本文件用什么编辑器都行——记事本、VS Code、Notepad都可以。写好以.c后缀保存后在Fluent里通过Define → User-Defined → Functions → Compile编译选好源文件然后Build简单函数也可以用Interpreted解释型不用装C编译器但语法限制多一些。如果你是2021R2之后的版本Fluent自带了UDF编辑器界面点Define → User-Defined → UDF Editor就能直接写。加载成功后回到User-Defined Scalars面板给UDS指定Diffusivity选择age_diffusivitySource Term选择age_source。第五步求解策略。我强烈建议两步走先把没有UDS的纯流场算收敛再把UDS方程加进去继续算。流场不收敛就加UDS很容易把空气龄的残差搞得很乱且难以判断到底是流场没收敛还是UDS设置有问题。UDS加入后残差曲线里会多一条uds-0一般收敛到1e-5以下就说明标量场已经稳定。3.3 迭代策略、保存在后处理前的关键操作很多人算完流场直接File → Write → Data结果后处理时发现梯度信息或者某些自定义变量没了这就要说“fluent计算保存梯度”这件事的性质了。在Fluent里梯度gradient本质上是在求解过程中为每个网格单元计算出来的量它被放在内存和case/data文件里。默认保存的data文件会包含梯度和单元数据至少在正常流程下通过Write Data生成的dat文件重新Read进后处理比如CFD-Post或Tecplot时是可以读出这些变量的梯度的。如果你发现打开后云图里没有梯度显示多半是在后处理软件里没勾选Compute Derivatives或者导出数据时格式选成了ASCII的简化格式。对于空气龄这种自定义标量在后处理前最好做几件事在Fluent的File → Write → Data保存之前确保UDS残差已经收敛如果要导出场数据给第三方后处理软件Tecplot、Paraview建议用File → Export → Solution Data选择需要的变量速度、温度、UDS-0即空气龄避免dat文件转格式时把UDS变量漏掉。如果只是想看云图那直接在Fluent的Graphics → Contours里选UDS-0画等值线就行不需要额外保存梯度。真正需要保存梯度的场景是你要沿某条路径计算速度梯度、壁面剪应力等派生量那时就用标准的Write Data流程保存然后在CFD-Post里用Calculator做梯度运算。4. 真实案例会议室置换通风的流场与空气龄分析4.1 案例基本设置我拿一个典型会议室开刀房间长12m、宽8m、高3.5m中间一张长桌桌两侧坐16个人简化人体模型每个约0.4m×0.35m×1.2m坐姿每人显热散热量按75W估算。送风方式采用置换通风——地板两侧各设一个置换通风末端送风风速0.25m/s送风温度19℃送风量合计约1500m³/h排风口设在房间顶部西侧两个400×400mm的排风口压力出口边界。这套系统的换气次数大约是1500/(12×8×3.5)≈4.5次/小时标称值不错。但从置换通风的角度看决定通风效果的关键不是换气次数而是送风能否“推”着室内污浊空气向上走并且人员呼吸区能否始终处于新风的覆盖范围内。网格策略全局基准尺寸100mm送风口附近加密到15mm人体及桌面表面加了边界层首层2mm增长比1.2共3层总网格量约620万。湍流模型用RNG k-ε流动为混合对流浮力强迫对流开启能量方程和重力密度用Boussinesq近似参考温度设为26℃。4.2 典型结果解读速度分布、空气龄分布、通风死角算完以后最主要的几张图里第一张是呼吸面高度z1.1m坐姿人员呼吸带高度的速度云图。置换通风在这个高度的速度特征应该是“低速、均匀”实测模拟结果大约0.1~0.15m/s但有两点异常靠近送风末端的一侧速度较高房间远端速度衰减明显人体热羽流在头部位置出现局部的上升气流速度约0.2m/s这使得呼吸区不再是单一的水平推送气流主导而是水平送风与人体热羽流的交汇。第二张是同一剖面高度的空气龄云图。结果很有意思送风口正对着的区域空气龄只有80~120秒但会议室远端角落的空气龄高达900秒以上。对比名义时间常数τn V/Q 336/0.417 ≈ 806秒也就是说远端角落的实际空气龄已经超过了“理论全场平均换气一次”的时间这在通风评估里是明显的死角信号。第三张是纵向剖面的空气龄云图y4m能看到置换通风经典的“分层”下部空气龄低200秒左右顶部靠近排风口处空气龄偏高600~800秒这是正常的——污染热空气上浮到顶部排出。但如果顶部出现局部高空气龄但排风口流速又低就说明排风短路或者部分气流在顶部循环。4.3 案例里的两个改进方向针对这个结果我做了两个对比工况一个是把总送风量从1500m³/h提高到1800m³/h另一个是保持1500m³/h不变但把远端送风末端的送风角度向角落调整同时增大末端出风面积。结果是教科书般典型的单纯加风量远端角落空气龄从900秒降到780秒左右改善有限而且近风口处平均速度升高到0.2m/s开始有明显的吹风感调整送风口之后空气龄分布均匀性大幅改善远端角落空气龄降到500秒左右人员区平均空气龄从420秒降到310秒换气效率从0.48提高到0.55。这个案例想说明的道理很简单空气龄仿真不是用来“算一个漂亮数字交差”的它是用来定位问题、指导风口布置和送风参数调整的。在复杂室内环境里风口的空间布置对通风效率的影响往往比单纯增加风量更大。5. 常见问题与排查技巧实录5.1 残差不收敛或震荡的排查室内通风流场最常见的收敛问题是残差在某个水平上反复震荡尤其是开启浮力之后。我遇到这种情况习惯按顺序排查第一看网格质量。负体积是致命的扭曲率大于0.9的单元如果数量多必须回炉重画。在Fluent里用Mesh → Quality查看最小正交质量低于0.1就说明网格有问题了。第二看浮力处理。开启Boussinesq近似时参考温度必须接近全场平均温度。参考温度给错了浮力驱动就会异常流场会像水烧开一样到处翻滚。另一个常见问题是Operating Pressure室内通风一般保持默认101325 Pa不需要改。第三调欠松弛因子。能量方程和k方程不太好收敛时把能量欠松弛从默认1.0降到0.9湍动能欠松弛从0.8降到0.6压力欠松弛保持0.3~0.5。尤其是加入UDS后如果uds-0残差乱跳给UDS单独设置欠松弛因子0.7左右效果立竿见影。第四流场本身存在物理上的不稳定比如自然对流占主导时本来就是周期性震荡的这种情况硬调松弛因子没有用需要判断是否应该转瞬态计算。5.2 空气龄负值、值不合理怎么处理UDS算出来的空气龄出现负值这个问题几乎每位做空气龄的工程师都遇到过。原因无非几个最常见的是出口边界条件给错。如果排风口边界把UDS的specified value设成了除0以外的数或者设了finite flux导致标量从外部往内流入负值就会出现。排风口正确的设置是流出通量自动计算不要强制指定标量值。第二个原因是UDS的残差没有收敛到足够低。空气龄标量是一个全场输运量在低速区收敛尤其慢。如果残差只降到1e-3就停了有些区域的值可能还在缓慢变化。在Continue迭代几百步后看单点空气龄监测值是否基本平了这个比只看残差更靠谱。第三个原因比较隐蔽扩散系数设置不合理。如果扩散系数给的过大UDS方程会表现出非常强的椭圆性使得某些区域出现“源头在上游却从下游反算”的非物理值。有效扩散系数一般控制在分子扩散量级到0.1 m²/s之间不要超过太多。5.3 高频疑问速查UDF、多孔介质、梯度保存、孤儿网格我把这些年被问到最多的几个问题的处理方法整理成了一张表问题快速处理办法UDF在哪里写、在哪加载用记事本/VS Code写.c文件Fluent中Define → User-Defined → Functions → Compile/Interpreted加载新版可用内置UDF Editor多孔介质参数怎么定粘性阻力系数1/α和惯性阻力系数C2由压降-流速曲线拟合无数据时参考经验值送风口滤网和均流板均可简化算完保存后云图丢了梯度后处理前用File → Write Data导出到Tecplot/Paraview时用Export → Solution Data勾选变量读入后需手动计算梯度Fluent里提示孤立区域或表面不连通回到CAD修布尔运算Import时尝试边界网格修复在Fluent里Display→Region确认连通性空气龄入口设了值还是不对检查边界类型速度入口的UDS指定值要为0若用的是mass-flow-inlet要额外设置UDS压力出口不要给UDS值模型明明有送风口但流场无风检查是否把风口设成了wallShow Mesh检查风口面元是否被其他面遮挡看进出口质量流量是否守恒做CFD本身就很考验细节空气龄这类用户自定义标量的计算更是如此。把UDS方程、扩散系数、边界条件每一个环节都抠清楚出来的结果才能经得起推敲。我在实际项目里最后还会做一件事把空气龄结果跟二氧化碳实测浓度做对照。空气龄低的位置实测CO₂浓度通常也低空气龄高的位置CO₂浓度对应偏高。数值上虽然没有严格一对一的关系但趋势对上了报告的可信度就上来了。如果后续还要扩展我建议可以把这套UDS思路延伸到污染物浓度预测、颗粒物沉降、飞沫传播风险分析上去思路完全一样只是把UDS标量换成组分方程或者DPM模型。基础打好了换场景只是改个方程的事。
