空调客车内气流组织与污染物浓度场数值模拟

向立平 ，王汉青

(1. 中南大学 能源科学与工程学院，湖南 长沙，410083；2. 湖南科技大学 能源安全与工程学院，湖南 湘潭，411201；3. 湖南工业大学 土木工程学院，湖南 株洲，412008)

客车是人们出行的主要交通工具之一，人们在车内度过的时间越来越长。目前，要求高速行驶的客车车厢内有很好的密闭性，车内环境主要依靠空调、通风系统维持。据调查统计，约65%的驾驶员在驾车时会出现头晕、困倦和咳嗽等现象[1-2]。车内污染可能是现代社会中人体健康的主要威胁之一，因此，创造舒适的乘车环境、降低能耗和运行成本具有重要意义。人们对此进行了大量研究[3-12]，但大多是集中于对轿车的空气温度场及车内和外界热交换热辐射等进行研究，而对空调客车气流研究较少，对车室内CO2质量分数分布的研究更少。在此，本文作者在考虑人体散热和太阳辐射对流场影响的基础上，研究人体呼出的CO2作为污染源对车内环境的影响，以期指导车室内气流组织优化设计和改善车室内空气质量。

1 模型的建立

1.1 空调客车物理模型的简化

研究对象是有49个座位的豪华空调客车。空调客车车室内部结构复杂，为此建立典型的简化客车结构，如图1所示。其几何参数为：车室内长为11.05 m，宽为2.50 m，高为2.79 m；玻璃窗为普通玻璃，厚度为5.00 mm；左右对称布置12排乘员座椅。对车室内空气流动进行数值模拟时，将蒸发室隔离出来；风道位于车顶左右两侧，送风口沿风道长度方向均匀布置；回风口位于车顶中央，通过走廊集中回风。

图1 空调客车车室简化模型Fig.1 Diagram of air conditioning vehicle

1.2 数学模型

空调客车室内的气流运动属于紊流运动，在紊流的工程计算中，k-ε湍流模型的应用最广泛，并取得了较好的结果[13]。这里采用的 k-ε湍流模型在车壁上符合对数定律，进行如下假设和简化：(1) 车厢内流动和传热视为稳态过程；(2) 车厢内空气视为不可压缩流，而且密度符合 Boussinesq假设；(3) 忽略固体壁面间热辐射，车内空气为辐射透明介质，不参与辐射；(4) 空调车室密闭良好，除送风口和回风口外没有空气泄露；(5) 忽略送风气流中 CO2含量，只考虑人体呼出的CO2含量；(6) 假定车内CO2做跟随运动，不影响空气的流动。在假设和简化的基础上，建立含连续性方程、动量方程、紊流动能方程(k方程)、紊流动能耗散方程(ε方程)、能量方程和组分输运方程等在内的一套封闭的方程组，其通用形式如下：

式中：各参量的意义和模型建立过程见文献[10, 13]。

2 边界条件和数值求解方法

在数值计算中，建立空调车室的边界条件特别是热边界条件是关键。整个车厢按定员考虑，人体散热量按 105 W/人计算，将这些热量均匀分布在座位上方，对于太阳辐射影响，采用Monte Carlo法[14]对太阳透射辐射能在车室内表面的分配比例进行计算，从而得出太阳透射辐射在空调车室内各表面引起的附加热流变化，以此2项作为能量方程的附加源项。

送风口边界条件如下：送风口垂直方向的气流速度w=2 m/s；u和v为相互垂直的两水平方向的气流速度，u=v=0；tin=18 ℃，k=0.02，ε=0.008。

出风口边界条件如下：边界为压力出口边界，出口压力为环境压力；k和ε为自由滑移参数。

对车室内气固耦合问题采用整体求解，固体区域黏性系数取大值，车顶及车两侧壁面、车室地板、行李架符合无滑移边界条件，即ui=0。车顶、地板及两侧壁面取定热流边界条件，给定车厢外部环境温度，根据不同的墙结构确定其换热系数。

污染源边界条件为：车内乘客呼出 CO2含量按0.014 4 m3/(h·人)计算[15]，将这些污染物作为CO2质量分数方程的附加源项。

为较好地逼近车室内复杂的结构，采用贴体坐标划分计算网格。为节省计算机内存和计算时间，在壁面处采用壁面函数法处理。应用有限差分法和交错网格离散控制微分方程，采用压力修正法即SIMPLE算法[13]求解离散控制方程：(1) 对流项采用迎风差分格式，将扩散项与对流项的影响系数分离，使方程绝对稳定；(2) 把相邻节点的影响系数表示成对流分量与扩散分量之和，将对流部分归入源项。

3 计算结果分析

因车厢几何结构左右对称，取车室室内左边结构各断面的速度场、温度场和 CO2质量分数分布组合，如图2～7所示，其中：图2～4所示为x=1.31 m和x=10.20 m时断面的速度、温度、CO2质量分数的分布图；图5～7所示为x=5.20 m和x=9.14 m时断面的速度、温度和CO2质量分数的分布图；图8所示为z=1.46 m断面的温度和CO2质量分数分布图。x=10.20 m时断面为驾驶员区，x=5.20 m时断面为车厢中部，且经过回风口。

图2 x为10.20 m和1.31 m时的断面速度场Fig.2 Velocity fields at x=10.20 m and 1.31 m

图3 x为10.20 m和1.31 m时的断面温度场Fig.3 Temperature fields at x=10.20 m and 1.31 m

图4 x为10.20 m和1.31 m时的断面CO2质量分数分布Fig.4 Concentration fields at x=10.20 m and 1.31 m

3.1 气流速度分布

由图2和图5可知：流场具有很强的对称性，不同断面的气流流线也几乎相似；送出的冷射流轴心速度垂直向下，扩散角度很小，在座椅处产生较大的流动涡旋，空气流动强度较大；回风口对空调车内流场影响较大，距回风口越近，流动强度越强，速度流线偏向回风口，风速逐渐趋近于回风口速度，回风口处风速大约为4.64 m/s；回风速度较大，从而导致风机噪声增加，影响乘客舒适性。

3.2 温度分布

由图3、图6和图8可知：由于空调客车车厢较长，驾驶区因太阳辐射温度较高，司机区的温度高达35.4 ℃；车厢尾部因太阳辐射和发动机传热温度稍高，达26.5 ℃。显然，驾驶区温度较高，应该增加驾驶区的送风量。由于送风口靠近窗侧，邻窗区的温度比邻过道区温度稍高；冷射流在座位区形成涡旋区，使座位区的温度比座位区上部和座椅下部温度高。由图 8可知：在乘坐区头部高度水平面内，除靠近驾驶区温度偏高外，温度分布均匀，比较理想，平均温度为26 ℃。

图6 x为5.20 m和9.14 m时的断面温度场Fig.6 Temperature fields at x=5.20 m and 9.14 m

图7 x为5.20 m和9.14 m时的断面CO2质量分数分布Fig.7 Concentration fields at x=5.20 m and 9.14 m

图8 z为1.46 m时断面温度场和CO2质量分数分布Fig.8 Temperature field and concentration field at z=1.46 m

3.3 浓度分布

根据国家标准《长途客车内空气质量要求》(GB/T 17729—1999)规定车内CO2质量分数不得高于0.15%，换算成质量分数应不得高于0.22%。由图4、图 7及图8所示的CO2质量分数分布图可知：在驾驶区域，CO2质量分数为 0.07%，车厢前、后部乘客区域的质量分数平均约为0.24%，而车厢中部CO2质量分数高达0.41%。车厢内CO2质量分数超出标准值，尤其是在车厢中部区域CO2质量分数严重超标。其原因是乘客区域人员密集，而回风口在车厢的车顶中部位置，车厢前后区域的污染物通过中部回风口排走，因此，中部区域污染物CO2质量分数大大高于前后区域的污染物CO2质量分数(见图8)；座位区是乘客的呼吸区，同时，由于空调系统送出的冷射流在座位区形成回旋气流，污染物受座椅的影响不容易排出，从而污染物更容易在座椅区聚集；车厢的上部和下部CO2质量分数较低，且上部区域CO2质量分数随高度增加而逐渐降低，但中部区是乘客主要呼吸区，污染物CO2质量分数聚集，CO2质量分数最高达0.40%(见图4和图7)，使车厢内乘客主要活动区空气质量变差，人体呼吸区空气比较浑浊，CO2质量分数偏高，严重影响了乘客的健康和舒适性。因此，应优化气流组织设计，科学利用通风换气，降低车内座位区污染物CO2质量分数。全面权衡车内空气质量、能量消耗和噪声等诸多因素，据控制车内空气质量达标、经济可行的原则，可以适当增加新风量。新风量的增加会使送风速度增大，从而引起能耗增加，噪声增大，同时，使车内微风速增大，引起乘客有吹风感。为了改善车厢中部CO2质量分数高于前后部CO2质量分数的措施，可以在车厢的前后部各设置1个回风口，减少污染物在车内的排出时间，以尽快置换人员区的污染空气，提高乘客的舒适性，防止人体健康受到危害。

4 结论

(1) 对空调客车室内气流组织、传热和污染物的数值分析，可以预测流场、温度场和CO2质量分数分布情况，筛选空调车室内气流组织设计最佳方案，从而缩短设计周期，节省人力、物力与财力。可以优化控制，为乘客提供舒适、健康的乘车环境。

(2) 空调通风系统布置对空调车室内微环境有较大影响，送、回风口的布置对空调车室内空气流场、温度场和CO2质量分数分布影响较大。送风气流在车内两侧形成比较大的回旋区，各对称面气流流型几乎相同；靠近回风口区域的流速较大，车内气流流线偏向回风口，且距回风口越近，流线偏转程度越大。

(3) 太阳辐射和人体散热对车厢内流场、CO2质量分数分布影响不大，但对车内温度场影响较大。由于受人体散热和太阳辐射的影响，车内座椅区温度高于座椅上部和下部温度，存在温度梯度；由于车厢前玻璃面积大，受到太阳辐射强烈，导致驾驶区和前部乘客区温度高于其他区温度。

(4) 车室内驾驶区的 CO2质量分数较低，车室内后部区CO2质量分数低于中部区域CO2质量分数。但座位区由于乘客密集，人体呼吸区污染物聚集，CO2质量分数偏高；因回风口位于车厢中部顶端，车室内中部污染物CO2质量分数严重偏高。

[1] 朱天乐. 微环境空气质量控制[M]. 北京: 北京航空航天大学出版社, 2006: 339-360.ZHU Tian-le. Microclimate air quality control[M]. Beijing:Beihang University Press, 2006: 339-360.

[2] 邓大跃, 陈双基. 汽车内空气污染研究综述[J]. 北京联合大学学报: 自然科学版, 2004, 18(2): 52-59.DENG Da-yue, CHEN Shuang-ji. Study summary of pollutants inside vehicles[J]. Journal of Beijing Union University: Natural Science, 2004, 18(2): 52-59.

[3] Tohru K. Analysis of vehicle passenger compartment ventilation using experimental and numerical models[J]. SAE Trans, 1997,89(6): 392-400.

[4] Han T. Three-dimensional Navier-stokes simulation for passenger compartment cooling[J]. International Journal of Vehicle Design, 1989(2): 175-186.

[5] Hara J, Fujitanj K. Computer simulation of passenger compartment air-flow[J]. SAE Trans, 1988(6): 1375-1382.

[6] China H, Yamgimoto K. Prediction of internal air flow analysis of passenger car by three dimensional numerical flow analysis[J].JSAE Review, 1987(12): 1376-1381.

[7] Currle J. Numerical simulation of the flow in a passenger compartment and evaluation of the thermal comfort of the occupants[J]. SAE Trans, 1997, 106(1): 806-816.

[8] Lin C H, Han T, Koromilas C A. Effects of HVAC design parameters on passenger thermal comfort[J]. SAE Trans, 1992,101(55): 209-231.

[9] Fujita A, Kanemaru J, Nakagawa H, et al. Numerical simulation method to predict the thermal environment inside a car cabin[J].JSAE Review, 2001, 22(1): 39-47.

[10] 靳谊勇, 郁永章. 铁路空调客车车内气流组织的数值模拟[J].制冷学报, 2002, 33(2): 30-34.JIN Yi-yong, YU Yong-zhang. Numerical simulation on airflow organizing inside the railway air-conditioning passenger car[J].Journal of Refrigeration, 2002, 33(2): 30-34.

[11] 连之伟, 张桂荣, 叶晓江. 铁路空调客车气流组织评价[J]. 上海交通大学学报, 2004, 38(6): 961-966.LIAN Zhi-wei, ZHANG Gui-rong, YE Xiao-jiang. Evaluation of air distribution in an air conditioned railway carriage[J]. Journal of Shanghai Jiaotong University, 2004, 38(6): 961-966.

[12] 杨培志, 顾小松, 傅俊萍. YW25G型空调硬卧列车车厢内气流数值计算[J]. 中南大学学报: 自然科学版, 2005, 36(5):888-891.YANG Pei-zhi, GU Xiao-song, FU Jun-ping. Numerical simulation of air distribution in YW25G train compartment[J].Journal of Central South University: Science and Technology,2005, 36(5): 888-891.

[13] 陶文铨. 数值传热学[M]. 西安: 西安交通大学出版社, 2001:430-447.TAO Wen-quan. Numerical heat transfer[M]. Xi’an: Xi’an Jiaotong University Press, 2001: 430-447.

[14] 卞伯绘. 辐射换热的分析与计算[M]. 北京: 清华大学出版社,1988: 79-88.BIAN Bo-hui. Analysis and simulation of radiation and heat transfer[M]. Beijing: Tsinghua University Press, 1988: 79-88.

[15] 陈焕新, 杨培志, 张登春. 列车车厢中CO2浓度控制系统仿真[J]. 湘潭矿业学院学报, 2002, 23(1): 132-136.CHEN Huan-xin, YANG Pei-zhi, ZHANG Deng-chun.Emulation about the control system of CO2concentration in train compartment[J]. Journal of Xiangtan Mining Institute, 2002,23(1):132-136.