吳文亮,袁魁浩,李智
(華南理工大學(xué)土木與交通學(xué)院,廣東廣州510641)
隨著改革開放的發(fā)展,我國城市規(guī)模逐漸擴(kuò)大。但城市熱島效應(yīng)也在逐漸加重,嚴(yán)重影響了城市的生態(tài)環(huán)境。道路等建筑在受到太陽輻射快速提升溫度,從而影響周圍建筑及路面上層空氣溫度,成為了熱島效應(yīng)非常大的影響因素。因此,為了緩解城市熱島效應(yīng)問題,改善道路路面溫度變得尤為重要。
針對(duì)如何改善道路路面溫度,國內(nèi)外學(xué)者大多數(shù)的研究方向分為兩大類:太陽熱反射涂層及透水路面。加藤寬道等[1-6]對(duì)太陽熱反射涂層材料、現(xiàn)狀、問題、未來發(fā)展趨勢展開了重點(diǎn)探究,但路面材料由于雨水沖刷,車輛磨損易使涂層剝落,涂層的耐久性問題一直未能得到解決。同時(shí),反射涂料只是單方面地降低路面的溫度,卻由于反射增加周圍建筑的溫度,在緩解城市熱島效應(yīng)方面未能取得滿意的效果。在透水路面方向上,張玉等[7-11]究了透水瀝青路面的材料組成、結(jié)構(gòu)組合、性能對(duì)比等方向,但透水路面由于本身空隙率大,在易儲(chǔ)存水分的同時(shí)也會(huì)吸附灰塵,使用久后影響水分蒸發(fā)降溫效果。
路面灑水降溫是目前我國最常見的道路降溫方法,在特定的時(shí)間通過灑水車灑水能有效降低路面溫度,同時(shí)清洗路面灰塵,相較于熱反射涂層材料和透水路面顯得更加簡便、有效及環(huán)保。本文擬通過室內(nèi)灑水試驗(yàn)及Fluent有限元模擬研究瀝青路面在風(fēng)速、灑水量、路面輻射吸收率和空氣濕度等影響因素變化下的灑水降溫效果,并根據(jù)影響特性探究實(shí)際一天內(nèi)的最佳灑水次數(shù)和最佳灑水時(shí)間。
影響瀝青路面灑水降溫的因素主要包括風(fēng)速、空氣濕度、單次灑水量和輻射吸收率。由于實(shí)際實(shí)驗(yàn)條件的限制,室內(nèi)試驗(yàn)中只研究分析風(fēng)速和單次灑水量對(duì)路面灑水降溫的影響,其余因素控制為常量。
本次實(shí)驗(yàn)場地選擇在華南理工大學(xué)道路工程研究所內(nèi)進(jìn)行,研究所內(nèi)因環(huán)境較為密閉,且有上面建筑物的遮擋,所以空氣濕度和太陽輻射量能維持在一個(gè)相對(duì)穩(wěn)定的值,利于采用控制變量法去研究。實(shí)驗(yàn)所需要用到的器材和功能見表1。
表1 實(shí)驗(yàn)器材及功能Tab.1 Laboratory equipment and functions
本次試驗(yàn)所采用的試件為瀝青混凝土車轍板,通過控制風(fēng)速和單次灑水量的變化,保持空氣濕度和路面輻射吸收率不變,灑水后記錄車轍板表面的溫度,以研究路面灑水降溫的影響特性。在材料選擇方面,其中瀝青選用SBS改性瀝青,集料采用經(jīng)初步篩選后的10~18、5~10、3~5 mm碎石和0~3 mm機(jī)制砂,填料采用普通礦粉。根據(jù)礦料篩分結(jié)果,算出密集配瀝青混凝土的配合比,并根據(jù)該比值合成最終瀝青混合料的級(jí)配見表2。
表2 瀝青混合料礦料級(jí)配組成及要求Tab.2 Asphalt mixture mineral grade composition and requirements
根據(jù)設(shè)計(jì)規(guī)范[12]要求制備車轍板300 mm×300 mm×50 mm,一共制備6個(gè)。試驗(yàn)開始前在所有車轍板側(cè)邊貼上錫箔紙,以確保側(cè)邊不受外部傳熱的影響。為對(duì)路面溫度進(jìn)行初始化,全部將車轍板置于原先設(shè)置為60 ℃的烘烤箱中。在試驗(yàn)地方架起碘鎢燈,使用太陽輻射測量器測量路面處的輻射值,通過反復(fù)調(diào)試碘鎢燈到車轍板的距離[13]使其達(dá)到一年最大的輻射值962 W/m2。
研究風(fēng)速對(duì)路面灑水降溫的影響時(shí),需控制單次灑水量為固定值,其余影響因素因試驗(yàn)場地原因可維持在一個(gè)固定值。從烘烤箱中取出1個(gè)車轍板,使用橡膠板作為墊面以確保底面不受外部傳熱影響。根據(jù)SASS[14]的研究,瀝青路面灑水后最大的吸附量1 kg/m2,超過1 kg/m2的部分會(huì)發(fā)生徑流。根據(jù)車轍板面積換算后1 kg/m2所對(duì)應(yīng)的灑水重量為90 g。實(shí)驗(yàn)開始前,打開可調(diào)速風(fēng)扇,使用風(fēng)速儀控制路面上方的風(fēng)速為2 m/s。稱取90 g水均勻噴灑到路面上,用測溫儀記錄路面表面隨時(shí)間變化的溫度。待路面恢復(fù)至穩(wěn)定溫度時(shí),結(jié)束本組實(shí)驗(yàn)。再次取出另一車轍板,重復(fù)以上實(shí)驗(yàn),依次記錄車轍板表面在風(fēng)速1 m/s、0 m/s下隨時(shí)間變化的溫度,繪制溫度—時(shí)間關(guān)系如圖1所示。
當(dāng)研究單次灑水量對(duì)路面灑水降溫的影響時(shí),需控制風(fēng)速為固定值。同樣從烘烤箱中取出1個(gè)車轍板,使用橡膠板作為墊面以確保底面不受外部傳熱影響。打開可調(diào)速風(fēng)扇,使用風(fēng)速儀調(diào)試路面上方風(fēng)速,使其保持在恒定值1 m/s,然后根據(jù)單次灑水量0.6 kg/m2換算的質(zhì)量為54 g,稱取54 g水均勻噴灑到路面上,用測溫儀記錄路面表面隨時(shí)間變化的溫度。重復(fù)以上實(shí)驗(yàn),分別記錄車轍板表面在單次灑水量為0.8 kg/m2和1 kg/m2下隨時(shí)間變化的溫度,繪制的溫度—時(shí)間關(guān)系如圖2所示:
通過圖1可以看出,當(dāng)風(fēng)速為0 m/s時(shí),曲線的變化趨勢較為平穩(wěn);而當(dāng)風(fēng)速為1 m/s和2 m/s時(shí),曲線下降速率明顯增快。風(fēng)速越大,路面溫度降溫幅度也越大,從曲線上表現(xiàn)為谷值越低,但路面恢復(fù)后穩(wěn)定值會(huì)越低。這是因?yàn)轱L(fēng)速除了提高水分蒸發(fā)速率外,還會(huì)促進(jìn)空氣的對(duì)流換熱,從而增大路面的降溫幅度。而幅度的增大,會(huì)增加路面溫度恢復(fù)穩(wěn)定值的時(shí)間。其次,在風(fēng)速持續(xù)影響下,空氣的對(duì)流換熱會(huì)使路面長時(shí)間保持在相對(duì)較低的溫度,故最終路面溫度穩(wěn)定值會(huì)隨著風(fēng)速的增大而減小。
圖1 風(fēng)速變化下溫度—時(shí)間關(guān)系圖Fig.1 Temperature-time relationship diagram under wind speed change
圖2 單次灑水量變化下溫度—時(shí)間關(guān)系圖Fig.2 Temperature-time relationship diagram under a single sprinkling volume change
通過圖2可以看出,在前期灑水量足夠的情況下,路面溫度下降速率不會(huì)隨著灑水量的增加而增加。但是隨著灑水量的增大,路面的降溫幅度和降溫時(shí)間會(huì)增加,相應(yīng)路面的恢復(fù)時(shí)間也增加。而從曲線最終的穩(wěn)定值可以看出,灑水量的增加,對(duì)路面溫度的穩(wěn)定值影響甚微。從以上曲線變化中可以看出,灑水量作為外部影響因素不會(huì)改變?yōu)⑺舭l(fā)吸熱的性質(zhì),包括降溫速率和溫度恢復(fù)程度。但是灑水量的增加,能延長水分蒸發(fā)降溫的時(shí)間,進(jìn)一步影響溫度的下降程度。
由于室內(nèi)試驗(yàn)的條件限制,只是對(duì)灑水路面影響因素中的風(fēng)速和單次灑水量進(jìn)行了分析,為更加充分地分析路面輻射吸收率及空氣濕度的影響特性,本文擬采用Fluent有限元模擬灑水蒸發(fā),同時(shí)確定在不同情況下的詳細(xì)灑水方案,包括灑水量、灑水次數(shù)和灑水時(shí)間。
由于Fluent系統(tǒng)內(nèi)未自帶灑水蒸發(fā)系統(tǒng),需選用合適的蒸發(fā)公式建立USD(user subroutine dflux)進(jìn)行分析計(jì)算。根據(jù)路面水分蒸發(fā)的外部條件和原理,采用質(zhì)交換準(zhǔn)則公式來模擬灑水降溫。
質(zhì)交換準(zhǔn)則公式中水蒸發(fā)的質(zhì)通量和熱通量由以下公式確定:
(1)
qe=MAγ,
(2)
式中:MA蒸發(fā)質(zhì)通量;hD為質(zhì)交換系數(shù);RA為空氣摩爾數(shù);T為空氣溫度;pA,1為蒸發(fā)面水蒸氣分壓;pA,2為空氣中水蒸氣分壓;qe為蒸發(fā)熱通量;γ為水的蒸發(fā)潛熱。
由于在蒸發(fā)過程中存在著風(fēng)速和汽車運(yùn)行過程中所帶的空氣流動(dòng),故可視為受迫對(duì)流,其中質(zhì)交換系數(shù)可由以下公式確定:
(3)
(4)
(5)
式中:D為質(zhì)擴(kuò)散系數(shù);D0為T0=273 K時(shí)的擴(kuò)散率;T0為初始空氣溫度;ν為運(yùn)動(dòng)粘度;L為定型尺寸;Sc為施米特?cái)?shù);Sh為宣烏特?cái)?shù);當(dāng)Re<5×105時(shí),Sh=0.664Re1/2Sc1/3;當(dāng)Re>5×105時(shí),Sh=(0.037Re0.8-870)Sc1/3,Re為雷諾數(shù)。
結(jié)合路面實(shí)際組成結(jié)構(gòu),于Ansys系列產(chǎn)品中的ICEM CFD建立300 cm×300 cm×10 cm的三維瀝青路面模型。同時(shí),為模擬空氣和路面的傳熱過程,根據(jù)路面尺寸在原來路面模型上方疊加模型300 cm×300 cm×90 cm,整體的模型尺寸為300 cm×300 cm×100 cm。建立物理模型后,在軟件中設(shè)置映射關(guān)系,并將長寬方向均設(shè)定為30節(jié)點(diǎn),高度方向設(shè)定為10節(jié)點(diǎn),劃分網(wǎng)格。
劃分網(wǎng)格結(jié)束后,將保存的文件使用Read功能導(dǎo)入進(jìn)Fluent中進(jìn)行求解設(shè)置。具體求解設(shè)置流程如圖3所示。根據(jù)殘差曲線中最終的分布,可見各指標(biāo)均趨于穩(wěn)定,因此可判定計(jì)算結(jié)果收斂。
圖3 Fluent各參數(shù)求解設(shè)置流程圖Fig.3 Fluent parameter setting flow chart
當(dāng)研究路面輻射吸收率對(duì)路面灑水降溫的影響時(shí),模擬設(shè)置時(shí)需保持其余變量不變,根據(jù)fluent設(shè)置流程,將風(fēng)速設(shè)置為1 m/s,單次灑水量設(shè)置為1 kg/m2,空氣濕度設(shè)置為20 %。在Fluent路面區(qū)域選定SourceTerms選項(xiàng)卡,在Energy中分別設(shè)置路面輻射吸收率為90 %、85 %、80 %和70 %,以模擬路面溫度在不同輻射吸收率下隨時(shí)間的變化情況。最終根據(jù)模擬計(jì)算結(jié)果繪制的溫度—時(shí)間關(guān)系如圖4所示。
當(dāng)研究空氣濕度對(duì)路面灑水降溫的影響時(shí),模擬設(shè)置時(shí)同樣需保持其余變量不變。其中風(fēng)速設(shè)置為1 m/s,單次灑水量設(shè)置為1 kg/m2,路面輻射吸收率設(shè)置為90 %。在邊界條件設(shè)置中選擇Inlet面,在Species選項(xiàng)卡中分別輸入空氣濕度值,以模擬計(jì)算空氣濕度為20 %、40 %、60 %、80 %時(shí)路面灑水后的溫度變化,經(jīng)匯總后的溫度—時(shí)間關(guān)系如圖5所示。
從圖4得知,路面輻射吸收率的變化對(duì)曲線的下降速率和恢復(fù)速率均影響不大。但是輻射吸收率增大時(shí),會(huì)使得降溫幅度減小,溫度穩(wěn)定值增大。這是因?yàn)檩椛湮章屎吐访嫖盏臒崃肯⑾⑾嚓P(guān),輻射吸收率越小,路面吸收的熱量也越小。而吸收熱量的減小表觀上為路面溫度更低,降溫幅度更大,溫度穩(wěn)定值更低。而降溫幅度越大時(shí),降溫時(shí)間會(huì)越長,整體路面的降溫持續(xù)時(shí)間也越長。
從圖5可得知,空氣濕度的增加,會(huì)降低路面溫度的下降速率和恢復(fù)速率,整體的降溫持續(xù)時(shí)間會(huì)增大。這是因?yàn)樗终舭l(fā)微觀上實(shí)則為空氣和路面的水分相互流動(dòng)的互逆過程,當(dāng)空氣濕度變大時(shí),會(huì)延遲兩過程趨向平衡的時(shí)間,表觀上為增大了降溫時(shí)間和恢復(fù)時(shí)間。但是,由于在太陽輻射的持續(xù)作用下,路面上的水分總會(huì)蒸發(fā)完,故當(dāng)灑水量不變時(shí)溫度的最低值和穩(wěn)定值是不變的。這一現(xiàn)象說明了空氣濕度對(duì)路面灑水的降溫幅度影響甚微。
圖4 輻射吸收率變化下溫度—時(shí)間關(guān)系圖Fig.4 Temperature-time relationship diagram under the change of radiation absorption rate
圖5 空氣濕度變化下溫度—時(shí)間關(guān)系圖Fig.5 Temperature-time relationship diagram under changes in air humidity
上文已經(jīng)對(duì)路面灑水降溫的影響特性進(jìn)行詳細(xì)地研究,為得到具體的灑水方案,還需通過Fluent分情況模擬一天內(nèi)路面的具體灑水時(shí)間和灑水次數(shù)??紤]到瀝青的最低軟化點(diǎn)在50 ℃以下,如AH-90為45.2 ℃[16]。故在灑水時(shí)機(jī)選擇上可設(shè)置為當(dāng)路面溫度達(dá)到50 ℃時(shí)需要灑水,這樣既能保證瀝青路面的結(jié)構(gòu)安全,也能起到降溫的作用。
計(jì)算前,保留之前的計(jì)算模型,并按照參數(shù)要求在具體的選項(xiàng)中設(shè)置。但在初始設(shè)置中,關(guān)掉Patch功能,表示一開始尚未灑水。創(chuàng)建點(diǎn)Point1(150 cm,150 cm, 10 cm),打開Monitor選項(xiàng)卡,將此點(diǎn)設(shè)為溫度監(jiān)控點(diǎn),以判斷灑水時(shí)機(jī)。另為了模擬一天內(nèi)氣溫和輻射量的變化情況,需在計(jì)算過程中按時(shí)暫停計(jì)算,并人為地補(bǔ)充更換設(shè)置。
在計(jì)算過程中,當(dāng)路面每次升溫至50 ℃時(shí)暫停計(jì)算,打開Patch功能設(shè)置往路面灑水再繼續(xù)計(jì)算。路面溫度會(huì)經(jīng)歷先下降后上升過程。待路面再次達(dá)到50 ℃時(shí)重復(fù)以上操作,直至路面不再回升至50 ℃。另外根據(jù)民用建筑熱工設(shè)計(jì)規(guī)范查閱一天每時(shí)段的輻射量和溫度,從10∶00到18∶00每隔1 h暫停計(jì)算更換路面輻射值和空氣溫度,以模擬真實(shí)一天的情況。同時(shí),灑水次數(shù)相對(duì)于灑水量來說成本更高,為探究更加合理的灑水方案,故灑水量均默認(rèn)為瀝青路面的最大吸附量1 kg/m2。
圖6 不同車流量下溫度—時(shí)間關(guān)系圖Fig.6 Temperature-time relationship diagram under different traffic flows
在路面水分蒸發(fā)的過程中,風(fēng)速對(duì)水分蒸發(fā)的影響和車輛在行駛過程中所帶動(dòng)空氣流動(dòng)對(duì)水分蒸發(fā)影響實(shí)質(zhì)上非常相似,均為通過影響空氣對(duì)流換熱來影響水分蒸發(fā),所以在研究不同車流量變化時(shí)可通過設(shè)置風(fēng)速來簡化控制。在實(shí)際變量控制中,通過控制風(fēng)速為2 m/s和0 m/s來簡化模擬車流量的“大”和“小”,其余變量設(shè)置為固定值。為方便統(tǒng)計(jì)灑水次數(shù)和灑水時(shí)間,建立“50 ℃灑水線”。當(dāng)溫度超過該線時(shí)開始灑水,超過該線的時(shí)間即為灑水時(shí)間,最終一天內(nèi)超過該線的頂點(diǎn)個(gè)數(shù)即為灑水次數(shù)。根據(jù)定義繪制的溫度—時(shí)間關(guān)系如圖6所示。
根據(jù)上圖分析,隨著車流量的減小,由于風(fēng)速越小路面升溫速率越快,故更加提前達(dá)到灑水點(diǎn)。其次,雖然風(fēng)速越小降溫時(shí)間越大,但是溫度會(huì)因?yàn)榭諝鈱?duì)流換熱較弱急劇上升,整體的降溫持續(xù)時(shí)間會(huì)更短。從一天內(nèi)比較可以發(fā)現(xiàn)12∶00~14∶00期間灑水持續(xù)時(shí)間更短,降溫幅度更小,這是因?yàn)樘栐诖藭r(shí)段的輻射值最大,溫度最高。最終結(jié)果:當(dāng)車流量較小時(shí),需要灑水6次。當(dāng)車流量較大時(shí),則需要灑水3次。
在研究不同使用年限路面灑水方案時(shí),根據(jù)LI等[17]的研究,瀝青路面的使用年限越長,其路面輻射吸收率會(huì)越小。故在實(shí)際變量控制中,通過控制輻射吸收率為0.9 m/s和0.75 m/s來簡化模擬路面的“新”和“舊”,其余變量設(shè)置為固定值。最終根據(jù)計(jì)算結(jié)果繪制的溫度—時(shí)間關(guān)系如圖7所示。
根據(jù)上圖分析,隨著路面輻射吸收率減小,路面所吸收的太陽輻射量也減小,升溫速度變慢,故較為延遲達(dá)到灑水點(diǎn)。同時(shí)路面使用年限變大時(shí),降溫時(shí)間和恢復(fù)時(shí)間都變長。最終結(jié)果:當(dāng)路面使用年限較小時(shí),需要灑水5次。當(dāng)路面使用年限較大時(shí),需要灑水3次。
在研究不同濕度地區(qū)路面灑水方案時(shí),在實(shí)際變量控制中,通過控制空氣濕度為80 %和20 %來簡化模擬不同地區(qū)濕度的“大”和“小”,其余變量設(shè)置為固定值。最終根據(jù)計(jì)算結(jié)果繪制的溫度—時(shí)間關(guān)系如圖8所示:
圖7 不同使用年限路面下溫度—時(shí)間關(guān)系圖Fig.7 Temperature-time relationship diagram for pavement under different service life
圖8 不同濕度地區(qū)路面下溫度—時(shí)間關(guān)系圖Fig.8 Temperature-time relationship diagram under different roads in different regions
根據(jù)上圖分析,當(dāng)空氣濕度變小時(shí),路面升溫的抑制效果減弱,故較為提前達(dá)到灑水點(diǎn)。同時(shí)隨著空氣濕度的減小,降溫時(shí)間和升溫時(shí)間變短,灑水間隔變短,灑水次數(shù)增多。最終結(jié)果:當(dāng)路面處于小濕度地區(qū)時(shí),需要灑水4次。當(dāng)路面處于大濕度地區(qū)時(shí),則需要灑水3次。
上文已經(jīng)對(duì)路面在不同情況下的灑水方案進(jìn)行了詳細(xì)地研究,但得出的灑水時(shí)機(jī)和灑水次數(shù)仍有些凌亂,需對(duì)上面結(jié)果進(jìn)行匯總。從不同情況下的溫度—時(shí)間圖(圖6~圖8)可以看出,最早的灑水時(shí)間為11∶02,最遲的灑水時(shí)間為15∶57。故可將所有灑水時(shí)機(jī)簡化選為11∶00~16∶00,灑水間隔計(jì)算下:
另外上述的灑水次數(shù)只是考慮了單方面的情況,當(dāng)多情況疊加時(shí),如大濕度地區(qū)車流較大的路段,可考慮按照最不利情況考慮,即按照灑水次數(shù)較大的情況來考慮,各情況疊加后灑水次數(shù)匯總見表3。
表3 多情況組合下灑水方案表/次數(shù)Tab.3 Multi-case combination sprinkling scheme
本文通過瀝青路面室內(nèi)灑水試驗(yàn)和Fluent有限元模擬對(duì)影響路面灑水降溫的各因素進(jìn)行了分析,并結(jié)合分析結(jié)果進(jìn)一步探究不同路面情況下的灑水方案,包括灑水時(shí)機(jī)和灑水次數(shù)。最終得出的結(jié)論如下:
①通過室內(nèi)試驗(yàn)得知,風(fēng)速越大,灑水后路面溫度降幅越大,降溫所持續(xù)的時(shí)間越長;單次灑水量越大,灑水后路面溫度降幅越大,降溫所持續(xù)的時(shí)間越長。
②Fluent有限元軟件結(jié)合質(zhì)交換準(zhǔn)則一定程度上能模擬灑水蒸發(fā)降溫的過程。
③通過Fluent模擬分析得知,路面輻射吸收率越大,灑水后路面溫度降幅越小,降溫所持續(xù)的時(shí)間越短;空氣濕度越大,降溫所持續(xù)的時(shí)間越大,但灑水后路面溫度降幅變化不顯著。
④當(dāng)車流量較小時(shí),需要灑水6次;當(dāng)車流量較大時(shí),需要灑水3次。當(dāng)路面使用年限較小時(shí),需要灑水5次;當(dāng)路面使用年限較大時(shí),則需要灑水3次。當(dāng)路面處于小濕度地區(qū)時(shí),需要灑水4次;當(dāng)路面處于大濕度地區(qū)時(shí),則需要灑水3次。
⑤灑水時(shí)機(jī)統(tǒng)一選為11∶00~16∶00,時(shí)間間隔視灑水次數(shù)而定。
⑥當(dāng)路面多情況疊加時(shí),按照最不利情況來決定灑水次數(shù)。