泄压口比率对管道内爆炸压力与火焰速度的影响.docxVIP

泄压口比率对管道内爆炸压力与火焰速度的影响 燃料爆炸一直是石油、化工、采矿和其他行业的主要灾难之一。因为管道具有长期的封闭空间特征。如果在其中发生局部燃烧，火焰将加快传播，并很快发展为爆炸，造成巨大的生命和财产损失。泄爆装置是防止可燃气体爆炸的手段之一,泄爆是指通过泄爆装置将腔室或容器内部已燃高压混合气体导入到外部环境中,使内部压力迅速降低以防爆炸的技术,开口泄爆是可燃介质流动与热化学反应耦合的非定常过程。泄放装置设计的关键是设计合理的泄放面积,以便在发生爆炸时能及时将爆炸超压泄放。 国内外学者的研究主要集中在内流场,在不同泄爆面积和不同泄爆压力条件下的容器内压力发展过程,讨论了不同泄爆条件对容器内爆炸发展产生的影响。尤明伟等用Fluent软件对球形容器泄爆过程中流场进行数值分析,对最终爆炸压力和爆炸强度的影响进行深入研究。 本文通过数值模拟的方法比较泄压口比率不同时得到火焰速度和压力波峰的变化,为泄爆装置的设计提供理论依据。 1 模型和数值方法 1.1 扩散速率的测量 可燃性气体爆炸是气体动力学流动与化学反应的综合作用过程。在建立数学模型时,该过程可以由流体力学基本方程组描述。 燃烧过程中的湍流过程,是气体爆炸机理的关键因素。国外广泛运用的描述模型为k-ε模型,该模型为湍流动能k和湍流动能的扩散速率ε的守恒方程: ??t(ρk)+??xj(ρujk)=??xj(Γk??xj)+τij?uj?xj?ρε(1)??t(ρε)+??xj(ρujε)=??xj(Γε??xj)+c1τij?uj?xj?c2ρε2k(2)??t(ρk)+??xj(ρujk)=??xj(Γk??xj)+τij?uj?xj-ρε(1)??t(ρε)+??xj(ρujε)=??xj(Γε??xj)+c1τij?uj?xj-c2ρε2k(2) 式中,ρ为密度;u为速度;τ为湍流切应力;Γ为湍流扩散因子;c1、c2为模型常数。 湍流强度u′和特征长度Lt分别为 u′=23k???√(3)Lt=C0.75μk3/2ε(4)u′=23k(3)Lt=Cμ0.75k3/2ε(4) 式中,Cμ为模型常数。 1.2 燃烧模型的建立 在建模时燃烧过程被看作一个由燃料-空气混合物转变成燃烧产物的简单的单向转变过程,以CH4为爆炸介质,其爆炸过程中的单步化学反应为 CH4+O2=CO2+2H20+Q(Q=880kJ/mol)CΗ4+Ο2=CΟ2+2Η20+Q(Q=880kJ/mol) 燃烧模型选用EBU-Arrhenius燃烧模型: Rfu=min(|Rfu,A|,|Rfu,T|)(5)Rfu,A=Bρ2w1w2exp(?ERT)(6)Rfu,T=CEBUρεkmin(w1,w2,w3)(7)Rfu=min(|Rfu,A|,|Rfu,Τ|)(5)Rfu,A=Bρ2w1w2exp(-ERΤ)(6)Rfu,Τ=CEBUρεkmin(w1,w2,w3)(7) 式中,Rfu为燃烧速率;B为阿仑尼乌斯因子;w1、w2、w3为燃料的质量分数;E为表观活化能;R为气体常数,kJ/(kg·K);T为气体温度,K;CEBU为计算常数。 1.3 网格划分的正交试验 实验管道物理模型为长方形管道,长为100 cm,截面尺寸为6 cm×6 cm,管道等效直径为10 cm,整个空间充满了甲烷-空气混合物,在左侧封闭端点火,在管道右侧设置泄压口,其它侧面均封闭,在下侧封闭端设置10个测点,测点均匀布置在管道中央,各测点间距为10 cm,采用与坐标轴平行的正交网格对整个计算域进行网格划分,网格数:100×10×10,如图1所示。 1.3.1 边境条件 壁面采用无滑移、无传热的边界条件,采用壁面函数法补充壁面区域的流动条件,采用非平衡壁面函数法处理壁面边界条件。 1.3.2 模拟的湍流系数 在左侧封闭端中心处全球点火,甲烷-空气混合浓度为9.5%;初始压力为101 325 Pa;初始温度为288 K;层流火焰拉伸系数Fs=0.15;无量纲湍流系数Ct=45;k-ε模型常数:Cμ=0.09,C1=1.44,C2=1.79,σk=1.0,σε=1.3,Fk=0.1(湍流动能转换系数)。 2 模拟结果分析 笔者分别模拟了泄压口比率为10%、20%、25%、30%、40%、60%、80%和100%(泄压口面积/管道截面积)的情况。 2.1 超压和火焰传播的阶段性关系 不同泄压口比率条件下爆炸压力与火焰速度的变化历程如图2所示。从图2中可以看出,超压的变化大致分为4个阶段:第1阶段(0~0.010 s),该测点处超压有规则地逐渐上升,是因为已燃气体释放出的热量使气体的温度骤然升高导致体积膨胀,进而诱导形成前驱压力波,前驱压力波使未燃气体发生扰动并推动气体向管道的出口端流动,测点处的压力值随前驱压力波的到来而上升;第2阶