豎直射流火撞擊障礙管道數(shù)值模擬分析

廖珊珊，張少剛，陶駿駿，劉家豪，汪金輝

（1 上海海事大學(xué)海洋科學(xué)與工程學(xué)院，上海 201306； 2 常熟理工學(xué)院材料工程學(xué)院，江蘇 常熟 215500）

引 言

隨著現(xiàn)代工業(yè)的發(fā)展，各種氣體燃料被廣泛運(yùn)用。美國(guó)能源局發(fā)布的《年度能源展望》中預(yù)測(cè)，在未來(lái)幾十年里，天然氣能源消耗將保持最高的年增長(zhǎng)速率[1-2]，高熱值的燃料資源提高了生產(chǎn)和工作效率。由于射流火的熱輻射通量相對(duì)較小，且隨著距離的增加而迅速降低，所以當(dāng)氣體燃料發(fā)生泄漏時(shí)，噴射火不會(huì)造成嚴(yán)重的后果。但在實(shí)際發(fā)生射流火災(zāi)時(shí)，管道、容器等設(shè)備，如果受到火焰沖擊，設(shè)備上的熱通量可能會(huì)非常高，可在短時(shí)間內(nèi)導(dǎo)致管道、設(shè)備等出現(xiàn)高溫、融化以及管道內(nèi)物質(zhì)的汽化與分解，從而引發(fā)多米諾骨牌效應(yīng)。如2016 年4 月，廣西防城港一化工廠管道泄漏引發(fā)火災(zāi)，而后火羽流撞擊鄰近管道，導(dǎo)致事故迅速蔓延和擴(kuò)大[3]。學(xué)者們已針對(duì)氣體燃料泄漏造成的射流火展開(kāi)大量研究[4]，其中射流火焰的形態(tài)特征是重中之重。氣體射流火的火焰長(zhǎng)度是研究火羽流行為演化特征的一個(gè)重要參數(shù)，它不僅給出了火焰發(fā)展的程度，而且可以用于確定輻射分?jǐn)?shù)，影響火焰的熱輻射以及傳播等特性[5-6]。

無(wú)量綱熱釋放速率Q*最早由Zokosiki[7]通過(guò)分析層流擴(kuò)散火焰羽流得到，Heskestad[8]引入與Q*有關(guān)的無(wú)量綱參數(shù)N，建立了火焰長(zhǎng)度與無(wú)量綱參數(shù)N的關(guān)系，該式能夠很好地描述平均火焰長(zhǎng)度。而后，Heskestad[9]對(duì)浮力主控的湍流射流火焰的火焰高度進(jìn)行了無(wú)量綱化，建立了火焰高度、燃燒器尺寸和Q*之間的關(guān)系。然而該火焰高度模型僅適用于湍流強(qiáng)度較弱的浮力主控射流火焰，并不能預(yù)測(cè)湍流強(qiáng)度較強(qiáng)時(shí)的火焰高度。Delichatsios 等[10-11]分析了浮力主控和動(dòng)量主控的擴(kuò)散射流火的演化行為特征，提出了基于Froude 數(shù)的無(wú)量綱火焰長(zhǎng)度計(jì)算公式，進(jìn)一步將火焰長(zhǎng)度預(yù)測(cè)模型擴(kuò)展到了動(dòng)量控制區(qū)域，并且發(fā)現(xiàn)隨著Froude 數(shù)的增加，熱輻射分?jǐn)?shù)前期幾乎不變，隨后減小，最后又幾乎不變[12]。實(shí)際上，射流火焰形態(tài)特征還會(huì)受到燃料類型[13]、環(huán)境氣壓[14-15]、環(huán)境風(fēng)[16-17]、多火源耦合[18-19]及泄漏口形狀[20]等多種因素的影響，而且射流火焰的最大寬度與長(zhǎng)度存在一定的比例關(guān)系，因燃料、環(huán)境氣壓等因素呈現(xiàn)不同比值，主要在0.189～0.119之間浮動(dòng)[4,21]。

上述研究主要集中在自由射流火焰上。在某些場(chǎng)景下，噴射火路徑中可能會(huì)出現(xiàn)一些障礙物，如管道、燃料儲(chǔ)罐、隔離墻或其他一些設(shè)備，導(dǎo)致噴射火撞擊障礙物[22-23]，火焰發(fā)展受壁面阻力、浮力和初始動(dòng)量的控制。因此，射流火焰的傳播機(jī)理和空氣卷吸行為比自由射流火焰更為復(fù)雜。Zhang 等[24-26]測(cè)量了矩形火源射流火撞擊到傾斜頂棚后火焰從撞擊點(diǎn)向上或向下的延伸長(zhǎng)度。Wang 等[27]通過(guò)實(shí)驗(yàn)研究了水平射流火焰撞擊壁面時(shí)火焰形態(tài)隨噴嘴-壁面間距的變化規(guī)律，發(fā)現(xiàn)隨著噴嘴-壁面間距的增大，火焰形態(tài)由“U”形轉(zhuǎn)變?yōu)椤癘”形。周夢(mèng)雅等[28]針對(duì)坑道限制條件下的水平射流火火焰長(zhǎng)度建立了基于Froude 數(shù)的預(yù)測(cè)模型。Kashi 等[29]通過(guò)設(shè)置不同直徑及擺放位置的障礙管道，發(fā)現(xiàn)火焰路徑中較小的設(shè)備對(duì)火焰形狀沒(méi)有太大影響，當(dāng)存在大型設(shè)備時(shí)，火焰會(huì)蔓延并改變其整體形狀。李玉星等[30]發(fā)現(xiàn)在其他條件一定時(shí)，隨著障礙物寬度的增加，火焰高度先減小后增大；隨著障礙物與泄露孔間距的增大，火焰高度先增大后減小。吳月瓊等[31]通過(guò)實(shí)驗(yàn)發(fā)現(xiàn)，與自由射流火相比，管道壁面的阻塞效應(yīng)會(huì)降低火焰的推舉速度。Foroughi等[32]通過(guò)實(shí)驗(yàn)發(fā)現(xiàn)火焰撞擊管道，會(huì)使管壁中產(chǎn)生極高的溫度，影響該管道的安全性，但僅關(guān)注了該過(guò)程的傳熱方面，未對(duì)火羽流的形態(tài)進(jìn)一步研究。Wang 等[33]通過(guò)水平射流火撞擊障礙管道實(shí)驗(yàn)，探究了泄漏速度、噴嘴直徑及管壁-火源間距對(duì)推舉距離的影響，發(fā)現(xiàn)隨著射流速度的增加，射流火焰撞擊圓柱的推舉距離比自由射流火焰的推舉距離增加得更明顯。

總地來(lái)說(shuō)，當(dāng)前關(guān)于火焰長(zhǎng)度[34-35]研究較為完善，但主要集中于受限空間屏障射流火，且屏障多為天花板或豎直壁面，有關(guān)氣體燃料泄漏誘發(fā)的射流火撞擊管道的研究則較少。然而在工業(yè)車間內(nèi)，各種管道縱橫交錯(cuò)，射流火焰撞擊管道后，羽流上升軌跡發(fā)生轉(zhuǎn)變，對(duì)空氣的卷吸強(qiáng)度發(fā)生變化，同時(shí)障礙管道將助長(zhǎng)火焰的蔓延，不但影響障礙管道的性能，還會(huì)使事故迅速蔓延和擴(kuò)大，造成更嚴(yán)重的后果。因此，本文擬探討射流火焰撞擊管道的演化特征規(guī)律，揭示不同障礙管道直徑和管壁-火源間距射流火形態(tài)特征演化規(guī)律。為現(xiàn)代工廠車間等場(chǎng)所中重要管道設(shè)備與燃?xì)夤艿腊踩g距的確定和火災(zāi)防治提供理論依據(jù)，對(duì)設(shè)施設(shè)備重點(diǎn)部位的耐高溫防護(hù)設(shè)計(jì)提供理論指導(dǎo)；同時(shí)，拓展射流火撞擊障礙物的火焰形態(tài)特征研究，進(jìn)一步豐富射流火災(zāi)研究理論。

1 數(shù)值模擬

本文運(yùn)用Fluent 軟件進(jìn)行模擬計(jì)算，F(xiàn)luent 軟件根據(jù)控制流體流動(dòng)的三大守恒定律推導(dǎo)出能量、動(dòng)量以及流體流動(dòng)所必須遵循的N-S 方程組計(jì)算。本文數(shù)值模擬采用三維穩(wěn)態(tài)求解器進(jìn)行計(jì)算，湍流模型選用標(biāo)準(zhǔn)k-ε湍流模型計(jì)算火場(chǎng)特征，該湍流模型通過(guò)求解湍流動(dòng)能(k)方程和湍流耗散率(ε)方程，得到k和ε的解，進(jìn)而計(jì)算湍流強(qiáng)度，最終通過(guò)Boussinesq 假設(shè)得到雷諾應(yīng)力的解。大量模擬驗(yàn)證表明，該模型能較好地計(jì)算射流火羽流有關(guān)特征?；瘜W(xué)反應(yīng)選用渦耗散模型，并使用標(biāo)準(zhǔn)壁面函數(shù)，輻射采用P-1 輻射模型，壓力與速度的耦合采用SIMPLE算法。

1.1 網(wǎng)格劃分與網(wǎng)格獨(dú)立性測(cè)試

在CFD 數(shù)值模擬中，為了能夠得到與實(shí)際情況更相符的模擬結(jié)果，需要有足夠精確的網(wǎng)格尺寸及較高的網(wǎng)格質(zhì)量。當(dāng)網(wǎng)格太大時(shí)，模擬計(jì)算結(jié)果粗糙；當(dāng)網(wǎng)格太小時(shí)，計(jì)算時(shí)間較長(zhǎng)且耗費(fèi)計(jì)算機(jī)資源。因此，選用合適的網(wǎng)格大小格外重要。本文采用六面體結(jié)構(gòu)化網(wǎng)格對(duì)模型進(jìn)行幾何劃分，對(duì)噴口及管道撞擊處進(jìn)行局部加密，并保證網(wǎng)格質(zhì)量均在0.4以上。網(wǎng)格大小尺寸以管道-火源距離1.0 m、管道直徑0.4 m為例，分別取網(wǎng)格數(shù)量為179萬(wàn)、209萬(wàn)、248萬(wàn)個(gè)，射流速度為300 m/s，得到如圖1所示溫度分布圖。相比179 萬(wàn)個(gè)的網(wǎng)格數(shù)量，另外兩組模擬結(jié)果更為可靠。綜上，為節(jié)約計(jì)算成本，本文選擇網(wǎng)格數(shù)量為209萬(wàn)個(gè)的劃分方式進(jìn)行網(wǎng)格劃分及計(jì)算模擬。

圖1 不同網(wǎng)格數(shù)量下障礙管道下方切面溫度分布Fig.1 Temperature distribution of the section below the obstacle pipeline under different grid numbers

1.2 模型準(zhǔn)確性驗(yàn)證

為驗(yàn)證計(jì)算模型的可靠性，本文對(duì)孔祥曉[36]研究的實(shí)驗(yàn)工況進(jìn)行模擬，噴嘴直徑設(shè)置為0.004～0.01 m，燃料豎直向上噴射。對(duì)于數(shù)值模擬射流火焰，前人用600、750 和800 K 為界限劃分火焰形狀，Huang 等[37]發(fā)現(xiàn)在數(shù)值模擬中用800 K 為界限來(lái)確定火焰形狀相對(duì)于實(shí)驗(yàn)值偏小，這可能是Fluent 數(shù)值模擬軟件中，重力對(duì)火焰形態(tài)的影響大于實(shí)際影響所導(dǎo)致的[38]?；诖?，本文對(duì)模擬結(jié)果分別以600、750、800 K 為界限來(lái)確定火焰形狀。模擬結(jié)果與實(shí)驗(yàn)結(jié)果相比發(fā)現(xiàn)，以800 K 為界限獲得的火焰形狀更符合實(shí)際情況，誤差均控制在9%以內(nèi)，在可接受范圍之內(nèi)。因此，本文可用上述模型來(lái)模擬射流火羽流演化行為特征。

圖2 射流模擬與實(shí)驗(yàn)結(jié)果對(duì)比Fig.2 Comparison between simulation and experimental results of jet fire flame length

1.3 模型工況設(shè)置及邊界條件

為探究撞擊管道火羽流演化行為特征，如圖3所示，本文設(shè)置了6 m×6 m×4 m 的全尺寸受限空間模型，模型上方設(shè)置為壓力出口，四周設(shè)置為墻面。為了便于模擬計(jì)算和分析，將模型內(nèi)簡(jiǎn)化設(shè)置，只保留噴口和被撞擊管道。環(huán)境溫度設(shè)置為300 K，壓強(qiáng)設(shè)置為一個(gè)標(biāo)準(zhǔn)大氣壓，甲烷沿豎直向上方向泄漏，考慮重力的影響且設(shè)置重力加速度為9.8 m/s2，噴嘴直徑設(shè)置為0.006 m。由于地面會(huì)影響火焰對(duì)空氣的卷吸情況[39]，為降低壁面邊界對(duì)火焰形態(tài)的影響，將噴口到地面的距離設(shè)置為0.1 m。

圖3 射流火模擬模型Fig.3 Simulation model of jet fire

表1 模擬工況設(shè)置情況Table1 Information of simulation tests

2 模擬結(jié)果及分析

2.1 障礙管道直徑的影響

2.1.1 火焰寬度 模擬結(jié)果顯示，障礙管道有效抑制了豎直方向的射流速度，同時(shí)增大了撞擊后火焰的橫向?qū)挾取?duì)于管道上方火焰寬度而言，如圖4所示為管壁-火源間距H為1.5 m，Q為334.3 kW，障礙管道直徑分別為0.2、0.4、0.6 m 時(shí)的溫度切片云圖。當(dāng)豎直向上射流火撞擊管道后，管道下方火焰寬度不變，管道上方火焰寬度明顯增大，這是由于火羽流撞擊管道后沿著管道橫向蔓延，增大了空氣的接觸面積，從而加劇燃燒，增強(qiáng)了對(duì)空氣的卷吸作用。而隨著障礙管道直徑的增加，障礙管道上方火焰寬度也明顯增大。

圖4 模擬溫度切片云圖Fig.4 Contour of temperature in different simulation tests

由于本研究中關(guān)注圓管障礙物對(duì)豎直射流火焰形態(tài)的影響，因此，將豎直射流火焰向上撞擊圓管后在圓管壁面伸展的最大寬度定義為火焰寬度，圖5為火羽流撞擊障礙管道后沿著管道蔓延的火焰寬度隨不同管道直徑的變化。不難發(fā)現(xiàn)，在相同熱釋放速率條件下，隨著障礙管道直徑D1增大，豎直向上射流火撞擊管道后沿管道蔓延的行進(jìn)距離也隨之增大。一方面，火焰表面積增大使得對(duì)空氣的卷吸效果加強(qiáng)，促進(jìn)了燃料與氧氣的混合；另一方面，障礙管道直徑的增大對(duì)火羽流豎直向上蔓延的阻擋效果增強(qiáng)，使火焰沿著管道水平蔓延，導(dǎo)致火焰寬度增大。當(dāng)Q為445.8 kW 時(shí)，火焰寬度隨障礙管道直徑的增加而增加得最多；當(dāng)Q為167.2 kW時(shí)，火焰寬度隨障礙管道直徑的增加而增加，但變化趨勢(shì)最為緩慢。當(dāng)障礙管道直徑一定時(shí)，不同熱釋放速率下的火焰寬度也不一致，這是因?yàn)楫?dāng)管道障礙直徑較小時(shí)，對(duì)火羽流的阻礙也較小，增加射流速度，火焰寬度變化也較??；而當(dāng)障礙管道直徑較大時(shí)，對(duì)火羽流的阻礙較大，增加火源熱釋放速率，火焰寬度變化也較大。

圖5 不同D1和Q下火焰沿管道蔓延的水平寬度變化Fig.5 Variation of horizontal width of flame spread along the pipe under different D1 and Q

2.1.2 火焰高度 在泄漏孔徑和管壁-火源間距不變的情況下，將管壁-火源間距設(shè)置為1.5 m 來(lái)模擬分析障礙管道直徑對(duì)射流火的影響。如圖6 所示，泄漏速度一定時(shí)，隨著障礙管道直徑的增大火焰高度逐漸減小。此外，與Wang 等[24]開(kāi)展的實(shí)驗(yàn)現(xiàn)象一致，火羽流始終包裹著管道蔓延，并未發(fā)生火羽流撞擊障礙管道而分離的情況，如圖7 所示。這是由于障礙管道表面邊界層具有黏性力作用，可燃?xì)怏w撞擊障礙管道后，一部分沿著管道水平蔓延，一部分包絡(luò)管道向上蔓延?？傮w而言，火羽流豎直向上的動(dòng)能分量減小，直徑越大的障礙管道火羽流豎直高度下降得越多，且火羽流水平寬度增加越明顯。

圖6 泄漏速度為300m/s時(shí)火羽流撞擊不同直徑管道溫度分布Fig.6 Temperature distribution of fire plume hitting pipes with different diameters when u=300m/s

圖7 水平射流火羽流撞擊障礙管道[24]Fig.7 Horizontal jet fire plume hiting obstacle pipe[24]

圖8 為管壁-火源間距為1.5 m 時(shí)，不同障礙管道直徑與不同泄漏速度對(duì)火焰長(zhǎng)度的影響，加入豎直向上的自由射流火作為對(duì)照組。在無(wú)障礙管道的情形下，燃燒火焰的高度最高，且火焰高度隨泄漏速度的增加而增大。在相同泄漏速度及管壁-火源間距的同一障礙物管道直徑條件下，增加障礙物后，火焰高度隨泄漏速度的增加而呈現(xiàn)出增長(zhǎng)的趨勢(shì)。從模擬結(jié)果來(lái)看，當(dāng)障礙管道直徑為0.2 m，Q為167.2 kW 時(shí)，障礙管道對(duì)火羽流的高度影響很小。而在相同火源功率下，障礙管道直徑為0.4～0.8 m 區(qū)間內(nèi)，障礙管道對(duì)火羽流高度影響則較明顯。這是因?yàn)樵诖嘶鹪垂β氏?，管道壁面距離火源1.5 m時(shí)，障礙管道處于火羽流后端，燃燒程度降低，導(dǎo)致管道對(duì)火羽流高度的影響也降低。當(dāng)障礙管道直徑為0.8 m，火源功率為334.3 kW 時(shí)，火焰高度反而出現(xiàn)降低的情況。從模擬結(jié)果看，此時(shí)火羽流在障礙管道的作用下向上的動(dòng)量撞擊管道，大部分沿著管道水平方向蔓延，剩下的少部分則沿管道壁面向上蔓延。整體來(lái)說(shuō)，當(dāng)管道直徑為0.8 m 時(shí)，障礙物有效阻擋了火焰豎直向上的蔓延程度，這將引起射流火燃燒導(dǎo)致周圍空氣的湍流強(qiáng)度加大，而空氣的湍流運(yùn)動(dòng)又將會(huì)帶動(dòng)燃燒火焰形態(tài)進(jìn)行無(wú)規(guī)則快速變化。因此，火焰高度在火源功率為334.3 kW時(shí)反而出現(xiàn)降低的情況。

圖8 管壁-火源間距為1.5m時(shí)火焰高度與熱釋放速率關(guān)系Fig.8 Relationship between the Lf/D and Q at H=1.5m

2.2 障礙管道壁面-火源間距的影響

2.2.1 火焰寬度 控制障礙管道直徑為0.4 m，隨著管壁-火源間距的增加，火焰撞擊管道處的水平寬度減小，即管壁-火源間距越小，火焰水平寬度也就越大。這種現(xiàn)象的產(chǎn)生是由于火焰前端的燃料未能充分燃燒，而障礙物的出現(xiàn)導(dǎo)致火羽流沿著管道向水平方向蔓延，增大了火羽流與空氣的接觸面積，從而加劇燃燒。位于火焰后端的燃料大部分已消耗殆盡，障礙管道雖也導(dǎo)致火羽流沿著管道在水平方向蔓延，但相對(duì)于較小的管壁-火源間距的屏障射流而言，其寬度增加的幅度較小。

圖9 不同管壁-火源間距溫度切片圖Fig.9 Temperature slice under different spacing between pipewall and vent of fire source

2.2.2 火焰高度 在控制障礙管道直徑為0.4 m 不變時(shí)，改變熱釋放速率和管壁-火源間距，研究障礙管道的位置對(duì)射流火羽流傳播特性的影響。如圖10 所示，隨著Q值的增大，火焰高度整體保持增加的趨勢(shì)。隨著管壁-火源間距的增大，火焰高度也呈現(xiàn)增長(zhǎng)的趨勢(shì)。當(dāng)管壁-火源間距為0.5 m 時(shí)，火焰高度最低；當(dāng)管壁-火源間距為1.5m 時(shí)，火焰高度明顯大于另外兩障礙射流組。這是因?yàn)榛鹧孀矒艄艿篮?，部分火羽流沿著管道蔓延，增大了與空氣接觸的面積，加劇了火焰燃燒。而管壁-火源間距越小時(shí)，火羽流沿管道蔓延的分量越多，導(dǎo)致大部分燃料在撞擊管道后沿水平方向蔓延燃燒從而被消耗掉，剩余火羽流包絡(luò)管道向上蔓延，火焰高度隨之下降。而管壁-火源間距較大時(shí)，射流火焰撞擊管道時(shí)未燃燃料的量偏小，火羽流具有相對(duì)較大的豎直運(yùn)動(dòng)速度，火焰更容易繞過(guò)管道向上運(yùn)動(dòng)，而不是沿管道水平延展，最終導(dǎo)致管壁-火源間距越大，火焰高度相對(duì)越高。

圖10 障礙管道直徑為0.4 m時(shí)火焰高度與熱釋放速率關(guān)系Fig.10 Relationship between flame height and heat release rate at the diameter of the obstacle pipe 0.4 m

2.3 單障礙管道影響下火焰高度綜合分析

國(guó)內(nèi)外學(xué)者對(duì)自由射流火的火焰長(zhǎng)度做了實(shí)驗(yàn)和預(yù)測(cè)，研究得出射流火無(wú)量綱火焰長(zhǎng)度與Froude數(shù)有關(guān)，如表2所示。

表2 目前關(guān)于豎直射流火焰長(zhǎng)度的部分研究Table 2 Some current studies on flame length of jet fire

本文射流火屬于動(dòng)量主控射流火，故采用基于Froude 數(shù)的無(wú)量綱火焰模型進(jìn)行擬合，其擬合情況如圖11所示，得到了一個(gè)動(dòng)量主導(dǎo)的甲烷氣體射流火焰長(zhǎng)度關(guān)聯(lián)式。

圖11 自由射流火數(shù)據(jù)擬合Fig.11 Fitting of free jet fire data

由此可以得出不同圓管障礙物影響下的射流火火焰高度表征模型為

圖12 障礙管道射流火數(shù)據(jù)擬合Fig.12 Data fitting of jet fire in obstructed pipeline

根據(jù)分析可知，式（4）建立在動(dòng)量主控射流火的基礎(chǔ)上，能很好地描述障礙管道直徑、管壁-火源間距對(duì)豎直向上射流火焰無(wú)量綱長(zhǎng)度的影響。

3 結(jié) 論

本文對(duì)豎直向上射流火焰撞擊障礙管道進(jìn)行了數(shù)值模擬，對(duì)比分析了不同熱釋放速率、不同障礙管道直徑和不同管壁-火源間距下火焰幾何形態(tài)的變化特征，通過(guò)對(duì)火焰寬度及長(zhǎng)度分析得出以下結(jié)論。

（1）障礙管道能夠使火焰沿著管道水平蔓延。在相同釋放速率和管壁-火源間距條件下，障礙管道直徑越大，火焰水平蔓延距離越長(zhǎng)，且管道上方火羽流寬度隨著障礙管道直徑的增大而增大。在相同釋放速率和障礙管道直徑條件下，管壁-火源間距越小，火焰沿管道蔓延距離越長(zhǎng)。

（2）障礙管道能夠有效減緩火焰縱向蔓延，火焰高度隨著障礙管道直徑的增大而降低，隨管壁-火源間距的減小而增加。

（3）通過(guò)加入影響因子A，建立基于Froude 數(shù)的計(jì)算公式，能夠較好地?cái)M合障礙管道作用下無(wú)量綱火焰高度。

本文僅對(duì)障礙單管道對(duì)火焰形態(tài)的影響做了探究，在后續(xù)的工作中會(huì)對(duì)多管道的間隔距離、擺放位置、直徑等因素進(jìn)行進(jìn)一步探究。

符 號(hào) 說(shuō) 明

ΔTf,a—— 火焰平均溫升

u—— 燃料泄漏速度，m/s

φ—— 化學(xué)計(jì)量比

ρe—— 燃料密度，kg/m3

ρ∞—— 環(huán)境空氣密度，kg/m3