王佳文,江 平,王 旭
(沈陽(yáng)航空航天大學(xué) 安全工程學(xué)院,遼寧 沈陽(yáng) 110136)
火旋風(fēng)是一種特殊的火焰行為,它是由燃燒產(chǎn)生的火羽流與周圍環(huán)境渦量場(chǎng)相互作用而誘發(fā)的劇烈燃燒的旋轉(zhuǎn)火焰?;鹦L(fēng)形成需具備3個(gè)條件,即有1個(gè)生成渦存在,并在生成渦內(nèi)存在1個(gè)流體匯,并存在某種摩擦力或拖拽力為地表面的空氣運(yùn)動(dòng)提供能量,以形成渦的底部邊界層[1]?;鹦L(fēng)火焰不僅具有很大的軸向上升速度,還有較大的切向旋轉(zhuǎn)速度,是類似于龍卷風(fēng)的強(qiáng)對(duì)流燃燒現(xiàn)象[2]?;鹦L(fēng)常出現(xiàn)在森林、草原與城市建筑火災(zāi)、一些中小規(guī)模的人工焚燒和大面積可燃液體池火中[3-4],如1923年發(fā)生在日本關(guān)東大地震震后的火災(zāi)中,在環(huán)境風(fēng)作用下,L形狀的被服廠著火誘發(fā)了火旋風(fēng),旋轉(zhuǎn)火焰脫離被服廠襲擊了在廣場(chǎng)上避難人群,分鐘內(nèi)造成38 000人死亡;2010年巴西圣保羅州南部阿拉薩圖巴地區(qū)麥田地發(fā)生火災(zāi),長(zhǎng)期干旱助長(zhǎng)了火勢(shì),并在環(huán)境風(fēng)的作用下向前移動(dòng),阻斷了道路;2012年澳大利亞愛(ài)麗斯斯普林斯地區(qū)出現(xiàn)了一個(gè)高達(dá)30 m的火焰龍卷風(fēng),其持續(xù)時(shí)間長(zhǎng)達(dá)40多分鐘[5-10]。當(dāng)火旋風(fēng)發(fā)生后,其燃燒速率,火焰高度,火焰溫度以及輻射強(qiáng)度都會(huì)顯著增加,會(huì)加快甚至改變火災(zāi)轉(zhuǎn)播和蔓延方式[11]。在實(shí)際多火點(diǎn)的火災(zāi)中,可能會(huì)出現(xiàn)由于多個(gè)火源相互影響從而使多個(gè)火旋風(fēng)同時(shí)出現(xiàn)在火場(chǎng)中的情況。因此,開展對(duì)多火源火旋風(fēng)的實(shí)驗(yàn)研究具有重要意義。
本文應(yīng)用自主研發(fā)的四面邊墻夾縫式多火源火旋風(fēng)實(shí)驗(yàn)裝置,研究實(shí)驗(yàn)室環(huán)境下以自然誘導(dǎo)方式形成的火旋風(fēng)。采用直徑為150,200,250 mm 3種尺寸的油盤,分別以80,100,120 mL的航空煤油作為燃料,在油盤中預(yù)先加入一定高度的水來(lái)保持液面上表面與油盤壁面基本持平。在實(shí)驗(yàn)中對(duì)3種尺寸的油盤進(jìn)行排列組合,全過(guò)程共進(jìn)行20次實(shí)驗(yàn),實(shí)驗(yàn)工況中油盤直徑大小和擺放位置如圖1所示。實(shí)驗(yàn)時(shí)采用攝像機(jī)監(jiān)測(cè)火旋風(fēng)的誘發(fā)過(guò)程,捕捉火焰實(shí)時(shí)動(dòng)態(tài)畫面以獲取火焰高度,并運(yùn)用Zukoski火焰高度模型來(lái)計(jì)算其特征火焰的火焰高度;用電子天平記錄油盤中剩余的燃料質(zhì)量,在其上方放置防火板,并使防火板上方的油盤與天平中心重合,為后續(xù)計(jì)算質(zhì)量燃燒速率提供數(shù)據(jù);同時(shí)采用風(fēng)速計(jì)來(lái)測(cè)量3個(gè)狹縫中心處的風(fēng)速,以獲得相應(yīng)的外加環(huán)量。
圖1 3油盤火旋風(fēng)裝置斜視圖Fig.1 Oblique view of the device of the fire whirls of three oil pan
實(shí)驗(yàn)裝置具體包括4塊規(guī)格為2 000 mm×500 mm×10 mm的防火玻璃、2塊規(guī)格為2 000 mm×600 mm×2 mm的防火板、1塊規(guī)格為2 000 mm×1 000 mm×2 mm的防火板、400 mm×400 mm×200 mm的電子天平支架、數(shù)個(gè)油盤墊板、自制的卡槽固定器以及電子天平和風(fēng)速計(jì)。具體實(shí)驗(yàn)裝置圖和效果圖如圖2~3。
圖2 3油盤火旋風(fēng)實(shí)驗(yàn)效果Fig.2 Experimental renderings of the fire whirls of three oil pan
圖3 3油盤火旋風(fēng)實(shí)驗(yàn)裝置尺寸Fig.3 Dimensions of the experiment device of the fire whirls of three oil pan
熱釋放速率是表征火災(zāi)危險(xiǎn)性的一個(gè)非常重要的參數(shù),它被定義為燃料的質(zhì)量燃燒速率與有效燃燒熱之積。m′為燃料的質(zhì)量燃燒速率(g/s),即單位時(shí)間內(nèi)通過(guò)燃燒反應(yīng)消耗的燃料質(zhì)量。圖4顯示的是第015次實(shí)驗(yàn)中1號(hào)盤燃料質(zhì)量隨時(shí)間變化的曲線。
圖4 油盤燃料質(zhì)量隨時(shí)間變化曲線Fig.4 The fuel quality of oil pan changes with time
圖4顯示了燃料燃燒過(guò)程中的質(zhì)量變化趨勢(shì),取準(zhǔn)穩(wěn)態(tài)時(shí)間段的數(shù)據(jù)進(jìn)行擬合運(yùn)算,得到第015次實(shí)驗(yàn)1號(hào)盤的火旋風(fēng)質(zhì)量燃燒速率絕對(duì)值為m′=1.28 g/s,同理可得各個(gè)工況下的質(zhì)量燃燒速率的絕對(duì)值,如表1所示。
首先對(duì)3個(gè)油盤形成的火旋風(fēng)內(nèi)部進(jìn)行分析,對(duì)比001,002,003次實(shí)驗(yàn)的3號(hào)油盤,如圖5所示??梢园l(fā)現(xiàn)在1,2號(hào)油盤直徑不發(fā)生改變時(shí),隨著3號(hào)油盤直徑的增加,其質(zhì)量燃燒速率不斷增加。同理對(duì)于005,006次實(shí)驗(yàn)的3號(hào)油盤,004,006次實(shí)驗(yàn)以及007,008次實(shí)驗(yàn)的1號(hào)油盤也可以發(fā)現(xiàn),對(duì)于3油盤形成的火旋風(fēng),在其他2個(gè)油盤直徑相同的時(shí)候,第3個(gè)直徑大的油盤所形成火旋風(fēng)的質(zhì)量燃燒速率明顯大于直徑小的油盤,說(shuō)明在3油盤形成的火旋風(fēng)之間,當(dāng)其他2個(gè)油盤直徑相同時(shí),另一油盤直徑的增加對(duì)其形成火旋風(fēng)的質(zhì)量燃燒速率有增大的作用。
表1 實(shí)驗(yàn)工況、質(zhì)量燃燒速率和風(fēng)速記錄
圖5 質(zhì)量燃燒速率與油盤號(hào)數(shù)關(guān)系Fig.5 Relation diagram of mass combustion rate and oil pan number
接著對(duì)比分析直徑相同油盤個(gè)數(shù)不同時(shí)的質(zhì)量燃燒速率差異,如圖6所示。
圖6 相同直徑不同油盤個(gè)數(shù)之間的質(zhì)量燃燒速率Fig.6 The mass combustion rate of the same diameters with different numbers of oil pans
對(duì)比3組數(shù)據(jù)發(fā)現(xiàn)隨著油盤數(shù)量的增加,明顯發(fā)現(xiàn)油盤的質(zhì)量燃燒速率依次遞增。說(shuō)明火源個(gè)數(shù)對(duì)于質(zhì)量燃燒速率是有相互作用的:火源個(gè)數(shù)越多,相互作用越明顯。
為了對(duì)實(shí)驗(yàn)數(shù)據(jù)做進(jìn)一步定量分析,引進(jìn)單位面積的質(zhì)量燃燒速率m″和速度環(huán)量Γ。
m″=m′/πr2
(1)
式中:m′為燃料的質(zhì)量燃燒速率,g/s;r為油盤的半徑,m。
本實(shí)驗(yàn)中火旋風(fēng)的質(zhì)量燃燒速率是燃燒準(zhǔn)穩(wěn)態(tài)階段數(shù)據(jù)的擬合直線斜率,因此在結(jié)合燃燒錄像視頻和質(zhì)量隨時(shí)間變化曲線的基礎(chǔ)上,確定了50~100 s的準(zhǔn)穩(wěn)態(tài)階段,觀察這段時(shí)間風(fēng)速計(jì)的變化,每5 s統(tǒng)計(jì)一個(gè)數(shù)值之后取風(fēng)速的平均值作為火旋風(fēng)自誘導(dǎo)渦量場(chǎng)的風(fēng)速。
本文將火旋風(fēng)看成是半無(wú)限長(zhǎng)的直線渦,根據(jù)“畢奧-薩伐”公式,速度環(huán)量如下。
Γ=2VP
(2)
式中:Γ為自誘導(dǎo)速度環(huán)量,m2/s;V為環(huán)境氣流在裝置入口處的風(fēng)速,m/s;P為形成單個(gè)火旋風(fēng)裝置的周長(zhǎng),2.06 m;單位面積的質(zhì)量燃燒速率和自誘導(dǎo)速度環(huán)量計(jì)算結(jié)果如表2所示。
表2 單位面積的質(zhì)量燃燒速率和速度環(huán)量記錄
火旋風(fēng)的研究中,周魁斌等人建立了火旋風(fēng)質(zhì)量燃燒速率公式和火焰高度公式,成功解除了燃燒速率與旋轉(zhuǎn)環(huán)量的耦合,并將結(jié)果與雷佼等人的進(jìn)行比較,可以很好擬合實(shí)驗(yàn)數(shù)據(jù)。本文以周魁斌的無(wú)量綱質(zhì)量燃燒速率與無(wú)量綱環(huán)量關(guān)系模型為參考進(jìn)行多火源火旋風(fēng)的質(zhì)量燃燒速率半經(jīng)驗(yàn)關(guān)系式的確定,分別用單油盤、雙油盤和3油盤形成的火旋風(fēng)質(zhì)量燃燒速率的數(shù)據(jù)與雷佼、周魁斌的數(shù)據(jù)進(jìn)行對(duì)比,數(shù)據(jù)分析見(jiàn)圖7。
圖7 質(zhì)量燃燒速率半經(jīng)驗(yàn)公式模型Fig.7 Model diagram of semi-empirical formula of mass combustion rate
如圖7所示,本實(shí)驗(yàn)中3組實(shí)驗(yàn)數(shù)據(jù)均呈線性擬合,說(shuō)明無(wú)量綱質(zhì)量燃燒速率模型適用于本次實(shí)驗(yàn),同時(shí)得出單油盤質(zhì)量燃燒速率半經(jīng)驗(yàn)公式,見(jiàn)式(3)。
lnm″*=-0.81+0.59lnΓ*
(3)
雙油盤質(zhì)量燃燒速率半經(jīng)驗(yàn)公式,見(jiàn)式(4)。
lnm″*=-0.16+0.41lnΓ*
(4)
3油盤質(zhì)量燃燒速率半經(jīng)驗(yàn)公式,見(jiàn)式(5)。
lnm″*=0.42+0.2lnΓ*
(5)
與雷佼、周魁斌的公式對(duì)比發(fā)現(xiàn),其半經(jīng)驗(yàn)公式斜率均大于本實(shí)驗(yàn)3組質(zhì)量燃燒速率半經(jīng)驗(yàn)公式斜率,這是由實(shí)驗(yàn)所用燃料類型的不同造成的。對(duì)比式(3),(4),(5)可知,隨著火源個(gè)數(shù)的增加,擬合線性的斜率也隨之降低。這說(shuō)明火旋風(fēng)之間的相互作用對(duì)于質(zhì)量燃燒速率有增強(qiáng)的作用。同時(shí)在進(jìn)行線性模擬時(shí),去掉了5個(gè)不符合準(zhǔn)穩(wěn)態(tài)階段質(zhì)量燃燒的數(shù)據(jù)。分別為第006次3號(hào)盤、007次2號(hào)盤、008次2號(hào)盤、008次3號(hào)盤和010次3號(hào)盤,其質(zhì)量與時(shí)間變化曲線與單油盤明顯不符,并且發(fā)現(xiàn)其油盤直徑均為250 mm,如圖8所示為燃料的質(zhì)量隨時(shí)間變化。
圖8 燃料的質(zhì)量隨時(shí)間變化Fig.8 The change of fuel mass with time
將火焰形態(tài)演化過(guò)程由點(diǎn)燃到熄滅劃分為4個(gè)階段:醞釀階段、生長(zhǎng)階段、準(zhǔn)穩(wěn)態(tài)階段和衰減階段。相比較于正常的單油盤火旋風(fēng)燃燒過(guò)程,在準(zhǔn)穩(wěn)態(tài)階段,其質(zhì)量都會(huì)有一段迅速均勻的下降過(guò)程,接著在衰減階段會(huì)有一次幅度很小的升高過(guò)程。在這5次實(shí)驗(yàn)過(guò)程中,衰減階段被大大延長(zhǎng),并且出現(xiàn)天平顯示的質(zhì)量不斷往復(fù)升高下降的現(xiàn)象。在這個(gè)階段火旋風(fēng)的形狀也會(huì)出現(xiàn)改變,火焰直徑會(huì)不斷變化,出現(xiàn)不穩(wěn)定的火焰形態(tài),直到火旋風(fēng)最終消失。
火旋風(fēng)為火羽流與環(huán)境渦量場(chǎng)相互作用所產(chǎn)生。在單油盤火旋風(fēng)形成后,拖拽力將圍繞一個(gè)穩(wěn)定值小幅度波動(dòng),這個(gè)階段即為準(zhǔn)穩(wěn)態(tài)階段。當(dāng)燃料逐漸耗盡時(shí),為旋轉(zhuǎn)提供的動(dòng)力減弱,則作用在液面上的拖拽力逐漸減弱,會(huì)使天平質(zhì)量上升,即衰減階段。但是在多火源火旋風(fēng)中,當(dāng)其中一個(gè)油盤燃料減少,進(jìn)入衰減階段時(shí),由于其它火源的渦量和拖拽力相互作用,會(huì)使這一衰減階段變得延長(zhǎng),質(zhì)量不斷往復(fù),使其火焰形狀發(fā)生改變。同時(shí)發(fā)現(xiàn)這5次實(shí)驗(yàn)中都是250 mm大直徑油盤質(zhì)量發(fā)生變化,說(shuō)明在火旋風(fēng)相互作用時(shí),大直徑火旋風(fēng)更容易收到周圍拖拽力和渦量的作用,直徑小的火旋風(fēng)反而會(huì)比較穩(wěn)定。
平均火焰高度的定義有多種,Zukoski等引入間歇率的概念來(lái)定義火焰高度,平均火焰高度為間歇率為0.5時(shí)對(duì)應(yīng)的H值。截取火焰穩(wěn)定期連續(xù)火焰圖像,通過(guò)等比例計(jì)算出火焰高度,并運(yùn)用Origin軟件計(jì)算其累積概率分布。累積概率為0.5時(shí)所對(duì)應(yīng)的火焰高度即為該火旋風(fēng)火焰的平均高度。
對(duì)D=250 mm的單油盤形成的火旋風(fēng)火焰高度進(jìn)行計(jì)算,結(jié)果如圖9所示。當(dāng)D=250 mm時(shí)火旋風(fēng)的平均火焰高度約為1 384 mm。
圖9 火焰高度的累積概率分布Fig.9 The cumulative probability distribution of flame height
同理運(yùn)用累積概率計(jì)算得出各個(gè)實(shí)驗(yàn)狀況下形成的火旋風(fēng)火焰高度,結(jié)果如表3所示。
表3 各個(gè)實(shí)驗(yàn)工況下的火焰高度
首先對(duì)油盤尺寸進(jìn)行分析,如圖10所示。
圖10 火焰高度與油盤號(hào)數(shù)關(guān)系Fig.10 Diagram of flame height with oil pan number
對(duì)比001,002,003次實(shí)驗(yàn),1,2號(hào)位置油盤尺寸相同,3號(hào)位置油盤分別放置150,200,250 mm油盤,發(fā)現(xiàn)隨著3號(hào)位置油盤尺寸增加,火焰高度有增加的趨勢(shì)。同時(shí)對(duì)比005,006次實(shí)驗(yàn)3號(hào)位置油盤及007,008次實(shí)驗(yàn)1號(hào)位置油盤,也發(fā)現(xiàn)火焰高度有明顯增加的趨勢(shì)。說(shuō)明多源火旋風(fēng)之間,油盤直徑的增大對(duì)其形成的火焰高度有增強(qiáng)的趨勢(shì)。
同時(shí)對(duì)尺寸相同時(shí)油盤個(gè)數(shù)不同的火焰高度進(jìn)行分析,如圖11所示。
圖11 相同直徑不同油盤個(gè)數(shù)之間的火焰高度Fig.11 Flame height between different oil pan Numbers in the same diameter
通過(guò)對(duì)比008,009,010次、008,009,016次和001,013,014次實(shí)驗(yàn)的火焰高度數(shù)據(jù),發(fā)現(xiàn)相同尺寸油盤,隨油盤個(gè)數(shù)增加,火焰高度呈上升趨勢(shì),說(shuō)明油盤個(gè)數(shù)對(duì)火焰高度有增強(qiáng)作用。
為了對(duì)數(shù)據(jù)做進(jìn)一步的定量分析,引入速度環(huán)量,并對(duì)其無(wú)量綱化。由Thomas火焰高度公式對(duì)火焰高度數(shù)據(jù)進(jìn)行處理,以周魁斌的無(wú)量綱火焰高度與無(wú)量綱環(huán)量關(guān)系模型為參考進(jìn)行多火源火旋風(fēng)的火焰高度半經(jīng)驗(yàn)關(guān)系式的確定,分別用單油盤、雙油盤和3油盤形成的火旋風(fēng)火焰高度的數(shù)據(jù)與雷佼、周魁斌的數(shù)據(jù)進(jìn)行對(duì)比,數(shù)據(jù)分析見(jiàn)圖12。
圖12 無(wú)量綱高度與無(wú)量綱環(huán)量半經(jīng)驗(yàn)公式模型Fig.12 The model diagram of semi-empirical formula with dimensionless height and dimensionless circulation
如圖12所示,本實(shí)驗(yàn)中3組實(shí)驗(yàn)數(shù)據(jù)均呈線性擬合,說(shuō)明無(wú)量綱火焰高度模型適用于本次實(shí)驗(yàn),同時(shí)得出單油盤火焰高度半經(jīng)驗(yàn)公式,見(jiàn)式(6)。
lnH*=0.05+0.35lnΓ*
(6)
雙油盤火焰高度半經(jīng)驗(yàn)公式,見(jiàn)式(7)。
lnH*=0.95+0.04lnΓ*
(7)
3油盤火焰高度半經(jīng)驗(yàn)公式,見(jiàn)式(8)。
lnH*=1.00+0.07lnΓ*
(8)
與雷佼、周魁斌的公式對(duì)比發(fā)現(xiàn),其半經(jīng)驗(yàn)公式斜率均大于本實(shí)驗(yàn)3組質(zhì)量燃燒速率半經(jīng)驗(yàn)公式斜率,這是由實(shí)驗(yàn)條件和燃料類型不同造成的。對(duì)比式(6),(7),(8)可知,隨著火源個(gè)數(shù)的增加,擬合線性的斜率也隨之降低。這說(shuō)明火旋風(fēng)之間的相互作用對(duì)于火焰高度有增強(qiáng)的作用。同時(shí)發(fā)現(xiàn)雙油盤、3油盤的火焰高度半經(jīng)驗(yàn)公式的斜率遠(yuǎn)遠(yuǎn)小于單油盤的斜率,這主要的原因是多火源火旋風(fēng)之間渦量和拖拽力的相互作用會(huì)對(duì)火旋風(fēng)的火焰高度帶來(lái)影響,這還需要做進(jìn)一步的理論分析。
1)單油盤火旋風(fēng)的無(wú)量綱質(zhì)量燃燒速率模型適用于本次實(shí)驗(yàn),同時(shí)得出多火源火旋風(fēng)質(zhì)量燃燒速率的半經(jīng)驗(yàn)公式,與實(shí)驗(yàn)數(shù)據(jù)有著較高的吻合度。且多源火旋風(fēng)質(zhì)量燃燒速率數(shù)據(jù)擬合直線的斜率小于相應(yīng)的單火源火旋風(fēng)數(shù)據(jù)擬合直線,擴(kuò)展了理論模型的應(yīng)用范圍。
2)隨著火源個(gè)數(shù)的增加,其質(zhì)量燃燒速率的斜率越低。
3)在多源火旋風(fēng)準(zhǔn)穩(wěn)態(tài)階段會(huì)出現(xiàn)明顯的相互增強(qiáng)現(xiàn)象,并且在衰減階段會(huì)出現(xiàn)劇烈的燃燒現(xiàn)象,天平質(zhì)量大幅度升降,并往復(fù)幾次,相比于單油盤火旋風(fēng)短暫的衰減期,衰減期的時(shí)間被大大加長(zhǎng)。
4)在火旋風(fēng)相互作用時(shí),大直徑火旋風(fēng)更容易收到周圍拖拽力和渦量的作用,直徑小的火旋風(fēng)反而會(huì)比較穩(wěn)定。
5)多油盤的火焰高度半經(jīng)驗(yàn)公式的斜率遠(yuǎn)遠(yuǎn)小于單油盤,這是由于多火源火旋風(fēng)之間渦量與拖拽力的互相作用的影響,需要對(duì)火旋風(fēng)之間渦量和拖拽力對(duì)火焰高度的作用做進(jìn)一步的理論分析。
[1]周魁斌.火旋風(fēng)的燃燒規(guī)律及其火焰移動(dòng)機(jī)制研究[D].合肥:中國(guó)科學(xué)技術(shù)大學(xué),2013.
[2]雷佼.火旋風(fēng)燃燒動(dòng)力學(xué)的實(shí)驗(yàn)與理論研究[D].合肥:中國(guó)科學(xué)技術(shù)大學(xué),2012.
[3]江平,王爽,劉穎,等.油盤直徑和高度對(duì)火旋風(fēng)特征影響的試驗(yàn)研究[J].中國(guó)安全科學(xué)學(xué)報(bào),2015,25(9):66-72.
JIANG Ping, WANG Shuang, LIU Ying, et al. Experimental research on influence of fire source heights and diameters on characteristics of fire whirls[J]. China Safety Science Journal,2015,25(9):66-72.
[4]余罕元.火場(chǎng)進(jìn)風(fēng)口情況對(duì)火旋風(fēng)特性影響的研究[D].南京:南京理工大學(xué)大學(xué),2014.
[5]王鵬飛.火旋風(fēng)的火焰與流動(dòng)特性研究[D].合肥:中國(guó)科學(xué)技術(shù)大學(xué),2015.
[6]Kuwana K, Sekimoto K, Minami T, Tashiro T, Saito K. Scale-model experiments of moving fire whirl over a line fire[J]. Proceedings of the Combustion Institute, 2013, 34(2): 2625-2631.
[7]Umscheid ME, Monteverdi JP, Davies JM, Photographs and analysis of an unusually large and long-lived fire whirl[J]. E-Journal of Severe Storms Meteorology, 2006, 1(2): 1-13.
[8]Soma S, Saito K. Reconstruction of fire whirls using scale models[J]. Combustion and Flame, 1991,86(3): 269-84.
[9]Kuwana K, Sekimoto K, Saio K, et al. Can we predict the occurrence of extreme fire whirls[J]. AIAA JOURNAL, 2007, 45(1): 9-16.
[10]Kuwana K, Sekimoto K, Saito K, et al. Scaling fire whirls[J]. Fire Safety Journal, 2008,43(4): 7-252.
[11]王爽,王振峰,張衛(wèi),等.外加環(huán)量作用下單油盤火旋風(fēng)實(shí)驗(yàn)研究[J].消防科學(xué)與技術(shù),2016,35(1):18-21.
WANG Shuang, WANG Zhenfeng, ZHANG Wei, et al. Experimental study on fire whirl of single oil pan under the plus circulation function[J]. Fire Science and Technology, 2016,35(1):18-21.