引水隧洞独头掘进工作面风流组织与CO扩散的模拟 Simulation on ventilation air flow and CO diffusion in leading face of excavation in tunnel.pdfVIP

水 利 学 报 SHUIU 2008年1月 XUEBAO 第39卷第1期 文章编号：0559．9350(2008)01．0121．07 引水隧洞独头掘进工作面风流组织与CO扩散的模拟 王晓玲1，陈红超1，刘雪朋1，张静2 (1．天津大学环境科学与工程学院。天津300072；2．天津大学建筑工程学院，天津300072) 摘要：建立了独头引水隧洞压入式通风紊流三维高雷诺数k-E数学模型，在动量方程中考虑了浮力的作用，采用 PISO算法求解非稳态的风流组织和CO浓度场。应用计算流体力学软件STAR—CD对云南南汀河引水隧洞独头掘 进工作面通风进行了模拟。结果表明独头引水隧洞掘进工作面风流流场可分为贴壁射流区，冲击射流区，回流区 和涡流区；工作面附近涡流区内、风管周围以及隧洞出口边壁处存在风流停滞区；分析了不同的通风时刻隧洞内 CO迁移和分布规律，并以Nakayama等的矿井通风甲烷浓度分布实验结果进行了验证，模拟结果与实验基本吻合。 关键词：引水隧洞独头掘进工作面；压入式风流组织；CO浓度分布；三维高雷诺数“e紊流模型；CFD模拟 中图分类号：TV554 文献标识码：A。 目前国内外对于隧洞通风进行了大量研究¨9。，对水电站引水隧洞施工期三维风流流动与CO扩散 分布的研究目前报道较少。本文根据流体力学和射流理论，利用三维缸e紊流数学模型来描述独头引 水隧洞掘进压入式通风不可压缩气体流动过程，在求解非稳态的风流组织和CO浓度场过程中采用了 with of PISO(Pressure-ImplicitOperators)算法，在非线性方程组迭代求解过程中采用了低松弛法来 Splitting 防止发散。与Nakayama…矿井实测甲烷分布进行对比，验证了CFD应用于隧洞通风和污染物扩散的可 行性。 1数学模型及求解方法 1．1风流模型描述黏性流体流动的风流模型由连续性方程和动量方程组成。 连续性方程： 誓+暑(1Dq)=0 、‘7 t)a1。a善；、广’7一 ( 动量方程： 1／11¨ (2) 警+壶(P∽)=一罄+鑫∽(誊+堕axl 参考密度，kg，m3。 1．2 CO输移模型CO输移模型描述了CO与其他组分在流体内质量传递的关系。 未(川+鑫(P吩匕)=杀(O丝y蔷)_形 (3) 收稿日期：2007．01．18 基金项目：国家自然科学基金；国家自然科学基金委员会、二滩水电开发有限责任公司雅砻江水电开发联合研究基金项 目；国家杰出青年科学基金资助项目 作者简介：王晓玲(1968一)，女，重庆人，副教授，主要从事环境数值模拟。E．rail：wmgxl@tin．edu．∞ 一12l一 万方数据 式中：k为质量分率；d，为紊流Prandtl数(一般取1．0)，E为时均组分的反应率，kg／m3·s。 1．3紊流模型紊流封闭通过两方程模型来实现。 后——输移模型：