雪峰山隧道通风三维CFD局部仿真优化.doc

雪峰山隧道通风三维CFD局部仿真优化 摘要：介绍了三维CFD仿真软件 三维CFD仿真软件Computational Fluid Dynamics)在近40年得到了迅猛发展。一方面是解决实际问题的需要；另一方面得益于计算机技术的飞跃发展和巨型计算机的出现。计算流体动力学是多种领域的交叉学科，它所涉及的学科有流体力学、偏微分方程的数学理论、计算几何、数值分析、计算科学等等。其基本思想是采用一定的计算方法（如差分方法等）求解流体力学N—S方程。基于计算流体动力学的基本思想，国际上出现了许多计算流体动力学通用软件，包括FLUENT、PHOENCIS 、CFX以及CFDesign等。课题组采用了国际上通用的计算流体动力学软件CFDesign对雪峰山特长公路隧道进行了仿真模拟分析计算。CFDesign是英国著名的BLUE BRIDGE NUMERICS 公司的主导产品。CFDesign软件主要是用来解决与热传导、流体运动相关的工程设计问题。 CFDesign3.1主要包含三个模块：CFDesign Modeler、CFDesign Solver和CFDisplay。 CFDesign具有解决与流体运动相关的以下一系列问题的功能和特征，包括： A. 二维平面、轴对称和三维空间流体问题 B. 可压缩和不可压缩流体的计算问题包括次音速、接近音速和超音速流体运动的问题 C. 稳定状态和短暂状态流体运动问题包括边界条件随时间变化的流体问题 D. 紊流和层流问题 E. 可以采用笛卡尔坐标和圆柱坐标进行分析 F. 牛顿流和非牛顿流的分析 G. 外部流体和内部流体 任何一个软件在正式运用到实际科学研究前，必须对该软件可靠性进行严格的测试计算。课题组在采用CFDesign进行具体的公路隧道通风仿真模拟前，对CFDesign软件进行了大量的测试和验算。通过建立一定的计算模型，同时采用理论计算和CFDesign数值模拟的方法进行分析。结果证明采用CFDesign进行仿真计算是可行的，同时其计算结果也是可信的。 几何参数： 圆管段断面半径：1m 圆管段长度：100m 计算参数： 圆管段入口流速：30m/s 圆管段出口压强：0Pa 计算模型如下： 计算结果如下表： 指标 断面 流量 平均速率 相对静压 绝对静压 绝对总压 入口 72.8229 23.5375 358.88 101682.8 102000.2 出口 70.0016 23.344 0 101325 101745 在计算过程中取管段壁面平均粗糙度为3mm。则根据现行《规范》可以有如下计算公式： 式中： △——平均壁面粗糙度（mm） ——隧道断面当量直径 通过以上公式可以计算出直线型、断面不变的隧道沿程阻力损失为： 式中： ——隧道长度 ——为隧道内平均速度 ——隧道断面当量直径 则可以计算如下： 则 而通过CFDesign计算到压力损失为：358.88Pa. 比较以上计算两者的相对误差为2.915%。 由此可见采用CFDesign进行公路隧道通风仿真模拟是可行的，同时其结果也是可信的。 模型侧视图 模型俯视图 建立的计算流体动力学（CFD）有限元分析模型中，为了能够尽量再现实际工程情况，模型采用的尺寸均参照雪峰山特长公路隧道通风系统设计所给定值。同时为了保证有限元分析能够收敛，对模型容易导致计算发散的局部区域进行适当调整。进行了2种工况的计算。 模型几何参数： 通风道长度：50m 通风道渐变段长度：4m 连接段长度：7.35 m 轴流风机断面积：3.10833m2 工况Ⅰ： 轴流风机流量：114m3/s 可得以下计算结果： 工况Ⅰ各断面计算结果表 指标 断面 流量 平均速率 相对静压 绝对静压 绝对总压 风道入口 114.856 5.17531 0 101325 101344 风道出口 116.382 7.16657 0 101325 101364 风机入口 112.361 29.4201 -1097.55 100227 100848 风机出口 114.481 29.4591 112.2258 101437.2257 102008.3 工况Ⅱ： 轴流风机流量：164m3/s 可得以下计算结果： 工况Ⅱ各断面计算结果表 指标 断面 流量 平均速率 相对静压 绝对静压 绝对总压 风道入口 164.236 7.39835 0 101325 101364 风道出口 166.258 10.2392 0 101325 101404 风机入口 160.516 42.0287 -2224.51 99100.5