仿真在线:有限元分析咨询与培训中心——Ansys|Fluent|Abaqus|HFSS|SolidWorks|课程培训

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗

2016-12-25  by:CAE仿真在線  來源:互聯(lián)網(wǎng)


者:陳升國

第一作者單位:重慶燃氣集團股份有限公司

摘自《煤氣與熱力》2015年8月刊


1 概述

近年來,天然氣燃燒爆炸(以下簡稱燃爆)事故頻繁發(fā)生,給人民生命安全和國家財產(chǎn)安全造成了極大威脅,嚴重影響天然氣行業(yè)的健康發(fā)展,也給政府主管部門、城鎮(zhèn)燃氣經(jīng)營企業(yè)造成了嚴重的社會負面影響。天然氣安全一直是政府、企業(yè)、社會高度關注的焦點問題。

通常情況下天然氣輸送管道是安全的,因為在管道內(nèi)輸送的天然氣濃度很高,不在天然氣爆炸極限范圍內(nèi),且管道內(nèi)不含氧氣或其他氧化劑,火焰很難在管道內(nèi)產(chǎn)生和傳播。但在特殊條件下,如施工不當、第三方破壞、地震或其他不可預見的意外事故,可能造成天然氣管道某處破裂,引起外界空氣進入管道,形成天然氣與空氣的混合氣。當該混合氣達到爆炸極限時,一旦遇上火源就具備了天然氣燃爆的必要條件。因此,如何預防和減少管道天然氣燃爆事故,已成為燃氣經(jīng)營企業(yè)安全管理亟待解決的問題。

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent分析圖片1

2 天然氣燃爆傳播理論分析

天然氣的主要成分為甲烷,具有易燃易爆特性,其燃爆的實質(zhì)是甲烷和空氣組成的爆炸性混合氣體在高溫誘導下發(fā)生的一種劇烈的化學反應并伴有大量的熱量生成,發(fā)生燃爆時化學反應方程式如下:

CH4+2O2CO2+2H2O

根據(jù)文獻可知:管道天然氣燃爆傳播實際上是以沖擊波方式傳播的,隨傳播時間和空間的推移,沖擊波結構發(fā)生了變化。在起始階段,以爆燃波(爆轟波)方式傳播,隨著甲烷和空氣混合氣體反應結束,最后演變?yōu)閴毫Σㄒ援數(shù)芈曀僭趩渭兛諝饨橘|(zhì)中傳播。

天然氣在管道中燃爆時,火焰面開始呈球面向外界擴張,并隨之向管道內(nèi)氣流方向傳播。由于管道壁面粗糙誘導湍流,使火焰發(fā)生皺褶,加快火焰的傳播,其運動速度將急劇增大,直到火焰竄出管道之后才會逐漸轉化為緩慢燃燒?;鹧娓浇鼩怏w因火焰加熱壓力升高,到一定距離(與天然氣和空氣的混合氣體量有關)形成壓力波,壓力波在傳播過程中,碰到管道壁面后發(fā)生反射,并在一定距離上形成平面壓力波。由此可以看出,天然氣燃爆傳播是火焰和壓力波共同傳播的過程。

3 管道天然氣燃爆數(shù)值模擬

FLUENT是目前國際上比較流行的商用CFD軟件包,只要涉及流體、熱傳遞及化學反應的工程問題都可以用它來進行模擬,因此選用FLUENT軟件對管道天然氣燃爆過程進行模擬。

3.1 建模和劃分網(wǎng)格

對管道天然氣燃爆過程的模擬計算采用流場模擬的方法,因此爆炸發(fā)生場所的幾何建模和計算網(wǎng)格劃分是必不可少的前處理工作。筆者使用Gambit軟件進行天然氣燃爆場所的幾何建模及網(wǎng)格劃分。

3.1.1 管道天然氣燃爆的幾何建模

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent結果圖片2

考慮到天然氣燃爆發(fā)生場所的復雜性,本次模擬采用結構網(wǎng)格和非結構網(wǎng)格形成的混合網(wǎng)格,實際計算中采用了不規(guī)則四邊形網(wǎng)格處理點火區(qū)域,傳播區(qū)域則采用矩形網(wǎng)格。結構網(wǎng)格可以提供更好的結構邊界相容性,因此管道的邊界區(qū)域使用結構網(wǎng)格劃分。局部網(wǎng)格劃分見圖1。

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent結果圖片3

1 局部網(wǎng)格劃分

3.1.2 模型初始條件和邊界條件設置

①初始壓力條件:發(fā)生燃爆時,火源位置產(chǎn)生的壓力相當大。鑒于此,設置點火區(qū)域壓力為1 kPa,其他區(qū)域壓力為0 Pa。

②初始溫度條件:點火區(qū)域為1 600 K,其他區(qū)域為300 K

③初始速度條件:整個區(qū)域初始速度為0 m/s。

④初始組成條件:為了簡化問題,空氣組成定為氧氣體積分數(shù)為22%,氮氣體積分數(shù)為78%

⑤為簡化初始條件,天然氣用CH4替代。

模擬管道壁面按典型的無滑移、無滲透邊界設定,模型左端設置為封閉端,右端設置為泄壓口,使之與實驗裝置盡量保持一致。

3.2 壓力波發(fā)展傳播過程分析

天然氣燃爆壓力波發(fā)展傳播過程見圖2,顯示了從天然氣被點火開始,壓力波以球形波開始傳播,經(jīng)過多次疊加之后逐漸形成平面波的全過程。

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent結果圖片4

圖2 天然氣燃爆壓力波發(fā)展傳播過程

從圖2可以看出,壓力波在點火源處以球形向四周傳播;在碰到壁面反射后,反射波與傳播過來的波相互疊加,疊加后壓力明顯升高,如ef中的紅色區(qū)域。隨著球形波不斷疊加,最終形成了平面波并向管道兩端傳播。左端壓力波在傳播到端頭后產(chǎn)生了反射,這樣反射波對原來向右端傳播的壓力波會起到一定的推動作用,最終導致右邊的壓力波逐漸增強。

3.3 火焰?zhèn)鞑ミ^程分析

3為從點火開始,火焰?zhèn)鞑サ某跏及l(fā)展變化情況。在天然氣被點火的初始階段,火焰以點火區(qū)域為中心向四周蔓延開來。其燃燒擴展是無規(guī)律性的,但是由于管道右端是開口,又受湍流機制影響,火焰向開口端發(fā)展的趨勢更強烈一些,因此右端的火焰區(qū)域要更大些。

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent培訓課程圖片5

圖3 初始階段火焰變化趨勢

4 數(shù)值模擬結果與實驗結果的對比分析

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent培訓課程圖片6

實驗和模擬一樣,管道內(nèi)天然氣體積分數(shù)調(diào)至9.5%,選用DN 700 mm的圓形管道,長度為93 m,共布置9個壓力傳感器和9個火焰?zhèn)鞲衅?壓力傳感器和火焰?zhèn)鞲衅鞣謩e在管道直徑的兩端相對布置,具體布置見圖4。圖4中相應的測點位置的尺寸標注是指該測點距封閉端的距離。

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent培訓課程圖片7

圖4 傳感器布置

4.1 壓力波傳播模擬結果和實驗結果的對比分析

天然氣體積分數(shù)為9.5%的混合氣在管道內(nèi)燃爆時,各測點實驗最大壓力和模擬最大壓力的對比見圖5。從圖5可以看出,壓力波的模擬結果與實驗結果存在一定差異,模擬值在距離管道最左端10~50 m范圍內(nèi)以及55~65 m范圍內(nèi)要明顯大于實驗值,可能與模擬時壁面條件的設定有關,數(shù)值模擬中壁面設定為光滑,而實際管道的壁面粗糙,因此天然氣燃爆傳播過程中損耗了部分能量,使得壓力有所降低。而在出口附近模擬值逐漸降低,這與實驗值逐漸升高的趨勢正好相反,這是由于實驗中在出口處進行了封膜處理,這就對出口附近的壓力波造成擾動,使得壓力升高,而模擬中出口處沒有做任何處理,因此泄壓后壓力是降低的。從曲線的變化趨勢來看,管道天然氣燃爆時最大壓力都是在爆源點附近先降低,然后上升到某一峰值后再逐漸衰減。

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent培訓課程圖片8

圖5 天然氣體積分數(shù)為9.5%的混合氣爆燃時各測點實驗最大壓力和模擬最大壓力的對比

4.2 火焰?zhèn)鞑ツM結果和實驗結果的對比分析

模擬和實驗時火焰在各測點的呈現(xiàn)時間對比見圖6。從圖6可以看出,模擬時和實驗時火焰呈現(xiàn)時間與測點距離的關系是一致的,都是隨著距離的增大,火焰呈現(xiàn)時間單調(diào)遞增,但是模擬的各測點火焰呈現(xiàn)時間要比實驗時火焰呈現(xiàn)時間長,這可能與管壁的粗糙程度有關系。

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent分析案例圖片9

圖6 模擬和實驗時火焰在各測點的呈現(xiàn)時間對比

管道天然氣燃爆時火焰?zhèn)鞑ニ俣扔嬎愎綖?

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent分析案例圖片10

火焰?zhèn)鞑ニ俣妊毓艿赖淖兓闆r見圖7。由圖7可以看出,由于點火位置不是在封閉端口中央而是在管壁靠近封閉端1.3 m處,造成火焰?zhèn)鞑ブ练忾]端口后形成了反射火焰,因此剛開始火焰?zhèn)鞑ニ俣扔兴鶞p小,但是隨著傳播距離的增大,火焰?zhèn)鞑ニ俣仍龃?達到峰值后速度降低,在靠近出口處,實驗值又有所上升,而模擬值降低。這可能是由于實驗時有封膜處理并且管壁粗糙以及測試處的極度凹凸造成了對火焰的擾動,而模擬時假設管壁光滑且無封膜處理,因此導致火焰?zhèn)鞑ニ俣葘嶒炛稻笥谀M值。出口處模擬值下降而實驗值上升是由于模擬時無封膜處理所致。因此整體來說,火焰?zhèn)鞑ニ俣饶M結果和實驗結果是基本吻合的。

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent分析案例圖片11

圖7 模擬和實驗所得的各測點火焰?zhèn)鞑ニ俣?/span>

5 天然氣燃爆的微觀解釋

Fluent軟件模擬出了隨著管道內(nèi)天然氣的逐漸消耗,其對應的燃燒反應速率、火焰波陣面以及壓力波陣面的變化趨勢,這是實驗無法觀測到的,見圖8。圖8a中藍色區(qū)域為天然氣已經(jīng)消耗的區(qū)域,紅色區(qū)域為天然氣未消耗的區(qū)域。

基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗fluent圖片12

圖8 燃燒反應速率、火焰波陣面及壓力波陣面與天然氣消耗量之間的關系

在模擬管道內(nèi)隨機截取了天然氣燃爆過程中5個不同時刻下各參數(shù)的云圖,時間分別為點火后第263、305、342、361、402 ms。由圖8可以看出,在某個時刻,當管道內(nèi)燃燒的天然氣消耗量達到某位置時,相對應的天然氣燃燒反應速率和火焰波陣面也到達同一位置,然而此時壓力波陣面是超前此位置的。這正好從模擬的角度合理解釋了管道天然氣燃爆時傳播中的兩波三區(qū)結構,即火焰波陣面和壓力波陣面形成了三個區(qū)域,分別為天然氣和空氣混合氣體的初始狀態(tài)、前驅(qū)沖擊波通過后的狀態(tài)、爆燃波陣面通過后的狀態(tài),且壓力波陣面超前火焰波陣面。

6 結語

隨著天然氣市場的不斷開拓,龐大的管網(wǎng)系統(tǒng)和多樣化的用氣環(huán)境給安全管理工作提出了更高的要求。筆者采用FLUENT軟件模擬出了管道天然氣燃爆的全過程,并結合模擬參數(shù)下的實驗數(shù)據(jù)對比分析了天然氣燃爆時壓力和火焰的傳播規(guī)律,給出了實驗無法觀測到的一些物理化學變化規(guī)律,為預防和減少天然氣燃爆事故提供了理論依據(jù),對天然氣行業(yè)安全技術研究具有一定的借鑒意義。


開放分享:優(yōu)質(zhì)有限元技術文章,助你自學成才

相關標簽搜索:基于FLUENT的天然氣燃燒爆炸數(shù)值模擬與實驗 Fluent培訓 Fluent流體培訓 Fluent軟件培訓 fluent技術教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學反應 fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析 

編輯
在線報名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進行在線報名,我們會及時給您回復電話,謝謝!
驗證碼

全國服務熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓QQ咨詢:點擊咨詢 點擊咨詢
項目QQ咨詢:點擊咨詢
email:kf@1cae.com




主站蜘蛛池模板: 芜湖液压配件-液压过滤器滤芯-水过滤器-芜湖新俊液压设备有限公司 | 日本国际高中_上海日本国际高中学校排名_日本国际高中留学课程_上海日语国际高中学校学费-上海工程技术大学国际多语种特色高中课程【官网】 | 密集架,密集柜,智能密集架,档案密集架,智慧档案馆(室)一体化建设 | 泰安兴润建材有限公司,泰安井盖定做,泰安警示桩定做,泰安雨水篦子定做,泰安操场篦子定做,泰安标志牌定做 | 上海物流公司_上海冷链运输_空运_电商仓储配送电话-迈泽物流 | 汽油发电机,柴油发电机,小型汽油发电机,小型柴油发电机,家用发电机生产厂家——上海东明动力设备有限公司 | 粮食加工设备_玉米_大米_面粉_燕麦_豆类杂粮加工设备-华豫万通 | 世博磁力泵,耐腐蚀磁力泵,耐酸碱磁力泵,氟塑料磁力泵-苏州冠裕机电科技有限公司 | 活性氧化铝球-活性氧化铝干燥剂厂家价格-3A-5A分子筛-巩义市亿洋水处理材料有限公司 | 耐磨焊条_高硬度堆焊焊条_碳化钨合金耐磨焊丝_北京耐默 | 实木中药柜,实木中药斗,木制中药柜,木制中药柜的价格,实木中草药柜,安国美佳中药柜厂家 | 钕铁硼-强力磁铁-永磁铁「聚盛磁铁」10年定制加工生产厂家 | 品牌T恤零售/批发厂商/供应商_男女式T恤生产/设计-圣达信T恤 | 普利卡管|普利卡管接头|普利卡接头-上海闵彬管业有限公司 | 选金银铅锌钽铌锡铂钨铑铱镐钯宝石&金属尾矿重选-青州冠诚重工机械 | 纸袋机|多层纸袋机|高速纸袋机|无锡市天天友情机械有限公司 | 蓝禹太阳能蓄电池 风能储能胶体铅酸电池-扬州东泰电源有限公司 | 乌海市腐植酸盐加工,内蒙古腐植酸钠,内蒙古型煤粘合剂生产厂家|创联腐植酸盐加工有限公司 | 螺带混合机|卧式螺带混合机|双动力混合机-无锡鑫海干燥粉体设备有限公司 | 手板机箱_亚克力机箱_医疗机箱_美容机箱_医疗设备外壳厂家_鸿堪医疗器械有限公司 | 上海国际餐饮博览会|餐饮供应|餐饮服务|餐饮加盟 | 旋风铣_数控旋风铣整机_旋风铣头|定制生产厂家【腾创机械厂官网】 | 营口新北方制糖有限公司| 郑州课桌椅|学生课桌椅|升降课桌椅批发|厂家|价格-新科教育用品 郑州井盖雨水篦子厂家-建联建材 | 平邑衡器厂--临沂衡器--电子汽车衡 -- 十佳知名品牌企业 | 上海钧尚电器有限公司 - Faulhaber电机 AMETEK pittman电机 AMETEK ROTRON军用航空风机 Exlar电动缸 MAE电机 MCG电机 CP电动工具 马头工具 AMCI驱动器 直流电机 减速箱 直流伺服电机,无刷电机,直线电机 直流防爆电机 防爆电机 汽车助力转向电机 EPS电机 faulhaber motor faulhaber gearbox NANOTEC电机 ELWOOD电机 PHYTRON电机 EXLAR伺服电动缸 高力矩、高性能直流电机,音圈电机,风机,直流风机,航空风机 | 手板模型-温州手板模-快速成型厂家-温州星科模具加工厂 | 易众拍卖行-事故车拍卖,残值车拍卖,水淹车拍卖,全损车拍卖,修复车拍卖,碰橦车拍卖,瑕疵车拍卖,报废车拍卖,泡水车拍卖,拆车件拍卖,配件拍卖,火烧车拍卖,二手车拍卖专业线上平台 | 专业色素炭黑生产厂家,提供各种用途色素炭黑价格-枣庄鑫源化工 | 轻型防化服|重型防化服|全封闭防化服|济南三安安全防护设备有限公司 | 眉山净源居环保科技有限公司,眉山除甲醛公司,眉山甲醛治理,眉山保洁服务,眉山家政保洁,眉山家电维修 - 眉山净源居环保科技有限公司,眉山除甲醛公司,眉山甲醛治理,眉山保洁服务,眉山家政保洁,眉山家电维修 | 圆锯机-滚牙机-滚丝轮-滚丝机-滚牙轮-切断机-东莞市溪远泰五金机械有限公司 | 生物安全柜检测,GMP设备确认,仪器性能确认,洁净厂房检测,仓储温湿度检测-上海熙迈 | 立式/卧式/暖气/管道泵,管道离心泵选型,管道增压泵型号 - 安平鼎千泵业制造厂 | 聚丙烯酰胺,聚合氯化铝,重金属捕捉剂,污泥调理剂,活性氧化铝,生石灰,反渗透阻垢剂,工业葡萄糖,硫酸铝,果壳活性炭,柱状活性炭,蜂窝活性炭,石英砂,锰砂-北京雁归来环保科技有限公司-以真诚为立足之本,以质量为生存之本,愿与海内外同仁共创双赢。雁归来人一路走来,气贯长虹,勇锐盖过怯弱,进取压倒苟安!我们紧扣时代脉搏,专注水处理、继往开来! | 江西省屹诚建筑装饰工程有限公司 | 南京自学考试网-南京自考报名| 深圳注册公司-工商注册代理-深圳注册公司流程及费用-记帐报税-公司注销-[深圳市悟空企业管理(深圳)有限公司] | 葫芦岛市鹏翔农药化工科技有限公司-苦参碱可溶液剂-氯氰菊酯乳油-阿维菌素乳油 | 宣传片拍摄-宣传片制作-三维动画-企业视频广告报价-北京环宇影音 | 商用车之家——观点有态度 热点不缺席 靠谱的商用车门户网站 |