Navier-Stokes方程在隧道通风改建中的应用.pdfVIP

矗盒^’ ■■一‘ 征庆横～ 玛一翻～ 《·一鲁大●■#±木工疆羲木研究所，●●t20∞”--青扛●古爿市盛碧t理魅。一■·317000) 一l-搴土用№嘶一＆o妇膏鲁柞为■重矗丑的曩事曩曩，采用三戎■值请重点基木■簟曩肯曩．对鼍扛啬爿I“■道冀土■ t童■曩ttz●柞Tt●l担和丹析．身优化■兵育叠量tT饭量． 羹■■，■薯．t^．tt●担 在囊曩许多早期毫成的公路雒道中，由于其车流量与原设计的流量相比较，有了成倍的增加．在这 种蕾况下，隧道内的废气和烟雾已不越够安全依靠自然风和车辆所产生的活塞风捧出隧道外，因此需要 改善鞋童的通风条件．为了合理经挤地增设通风设施，因此有必要对遥风方案进行计算机数值模狠，为 公麝蔓遣叠风方案构造摔和优化提侠定性和定量的分析． 一、Navier—Stokes方程及其CIP解法 在文t[1]中介绍了鞋道通风数值摸拟的三种数掌模噩，其中模型z是将隧道内的风泷视为非定常 牯性不可压的藏体，褥掰了如下矢量形式的连续性方程和N¨h—Stokb方程． l 可·U 0 鼍+(u。审)u掌一VP+壶叫 其中，可·u裹示矢量U的鼓虚，(U·守)U表示肘藏项。v为梯度算子t△为词和算子，■诺鼓Re 、，J =孚．v为空气的运莉桔性系数，V为来藏蘧庄，d为隧道的当量直径． 求■N¨岫一Stok∞方程，一般墨将速度蔫和压力场分开来处理．常用的数值算法有MAC算法、 (Cm)作一筲要舟绍． InterpoIated Method)方法，是由T．Yabe 三次羹韭赞质点方法简称CIP(Cubk P靶udo—plrtick 等人于八十年代提出并发晨的用于求解双曲曩蕾t分方程的一种羲值方法．CIP方法精度高、数值扩散 小、I膏便直曩、面且要求存麓量比较小．当■诺致报大时，Navier—stokes方程的类噩趋寓于双曲塑方 程，因此可以用Cm方法来隶■。用c球方法隶鼻N丑vier—stokes方程，嚣为其主要过程采用的是曼格 式，所以不必束●I蛾性方程组，从两可以加快诗算速度． cIP方法基奉恩越是将Nyier—stohs方程 警+(u·V)￡，葺一VP+忐凹 分■成事移灌硬和穆藏项两部分： ，r 1 ‘ 著2一VP+忐叫 (1) V)c，置O (2) 警+(u· 对于非芬童嘎(1)采用差分法或者有限元盛来采■．因为谈式右■同耐古有压力项和精性项，所以将其 分成两步来处理，取时间步长为△t．即 鼍笋一vP．．掣一恚彬+，其中，凹。一警 对于咎蠢嘎方程(2)，采用cIP方法来求鲁(详觅[4])，量后得旁U“1． 一136一 一————1l_■■≯。耵I_ 二、浙江台州104国道黄土岭隧道通风改建方案的模拟 104一连黄主岭奠遵由于其净膏有限。为此设矗在蠢道外北蛙增t一十弓l舅，采甩毫注羹鲁量风方 式．但选一蠢风方式在国内一首次采用，器耍撵人研究gI景形状和风机位量财奠道内风建风压分布的髟 畸．在这里。将风机嚷出的气藏橇为粘性不可压藏体，人口条件和韧妨条件嵌据风机出口喷速决定，出口 条件采用一阶外推格式，由于薯遭譬面固定，所以连界■子t面无滑移情况．用第一节中奔绍的C疋方 法采■非定1|不可压的N¨ier—Stok∞方程，来分析奠道引涓内的风压风速． 由于此通风方式是在引稠口均匀布重风机，风压风毫在隧道竞度方向上应基本一致，所以我们可根 据不同的引■赢度和长度以及风机位置进行二壤数值蕞擅．下面以无量纲标量形式蛤出基奉方程和韧 始条件、边界条件，即 ． 皇+皇。o