• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    壁面綠化及熱效應(yīng)對淺型街谷內(nèi)污染物擴散與轉(zhuǎn)化的影響研究

    2022-09-21 01:04:02黃遠東王可心劉宇辰崔鵬義
    上海理工大學(xué)學(xué)報 2022年4期
    關(guān)鍵詞:背風(fēng)面旋渦順時針

    黃遠東,王可心,劉宇辰,崔鵬義,羅 楊

    (上海理工大學(xué) 環(huán)境與建筑學(xué)院,上海 200093)

    街道峽谷是一種類似自然峽谷的城市結(jié)構(gòu),由兩側(cè)連續(xù)建筑物以及中央的區(qū)域構(gòu)成[1]。較低的天空視野系數(shù),使得日間太陽輻射存儲于兩側(cè)建筑壁面之中,加劇了熱島效應(yīng)和光化學(xué)反應(yīng);半封閉的街谷結(jié)構(gòu)進一步降低機動車尾氣污染物擴散效率,嚴重威脅臨街居民和行人的身體健康[2-3]。

    壁面綠化是近年來發(fā)展起來的一種新型立體綠化模式,研究表明其可以有效改善城市微氣候環(huán)境。在不同天氣下,Cuce[4]對配置綠化的壁面和無綠化壁面進行了溫度實測,結(jié)果表明:在晴天時,壁面綠化能夠使墻體的平均溫度降低約6.1 K;多云天氣時,則能夠使墻體的平均溫度降低約4.0 K。Schettini 等[5]進行的場地實測也得到了類似結(jié)果:在日間溫度最高的時刻,壁面綠化能夠降低壁面的表溫約3.0~4.5 K。Basher 等[6]選用長壽命且能夠承受更多太陽輻射的翼豆植物作為壁面綠化的主體,發(fā)現(xiàn)其能夠降低壁面平均溫度約2.4 K,且壁面最大下降溫度約6.4 K。Jayasooriya等[7]進行的城市尺度環(huán)境影響研究表明,雖然壁面綠化對顆粒和氣態(tài)污染物的去除能力不如傳統(tǒng)樹木,但其較小的占地面積和潛在的凈化能力使其應(yīng)用前景廣闊。Ysebaert 等[8]研究發(fā)現(xiàn),壁面綠化對顆粒污染物的截留效果取決于顆粒物在街谷內(nèi)的停留時間;在繁忙的主干道上,每平方米葉面積的壁面綠化能夠截留約1010~1011個顆粒,可以有效降低街谷內(nèi)的顆粒污染物濃度。Viecco 等[9]使用ENVI-met 模型就壁面綠化對PM2.5的去除效果進行數(shù)值計算,發(fā)現(xiàn)最大去除率可以達到7.5%,其主要受綠化覆蓋率、污染源位置和建筑物高度等影響。

    綜上所述,有關(guān)研究主要聚焦在壁面綠化對街谷內(nèi)顆粒污染物濾除和熱島效應(yīng)的緩解,并且側(cè)重于寬高比為1 的標準街谷,以及單獨考慮壁面綠化或太陽輻射熱效應(yīng)對峽谷內(nèi)流場和污染物擴散的影響。近年來,大氣光化學(xué)反應(yīng)產(chǎn)生的O3污染日漸嚴峻,但有關(guān)街谷內(nèi)綠化對反應(yīng)性NOx-O3動力學(xué)特性影響的研究卻沒有引起足夠重視。鑒于此,本文利用經(jīng)風(fēng)洞實驗驗證的三維CFD 模型,探究壁面綠化與熱效應(yīng)的耦合作用對淺型街道峽谷內(nèi)污染物擴散與NOx-O3反應(yīng)動力學(xué)的影響。該研究可為城市規(guī)劃和壁面綠化種植提供指導(dǎo),以緩解城市熱島效應(yīng),改善局部空氣質(zhì)量。

    1 研究方法

    1.1 物理模型

    如圖1 所示,研究采用三維全尺度街道峽谷。峽谷的高寬比AR和長高比分別為0.5 和20.0,亦稱為淺型長街道峽谷。圖1(a)和圖1(b)展示了街道峽谷的幾何模型:每個建筑物的長、寬、高分別為10H(H=20 m),1.0H和0.5H,兩棟建筑物的間距(即峽谷寬度W)為1.0H。圖1(b)還給出了街道峽谷內(nèi)部壁面綠化和近地面污染源的具體布置情況:壁面綠化被種植在上風(fēng)向建筑物的背風(fēng)壁與下風(fēng)向建筑物的迎風(fēng)壁上,厚度為0.01H(20 cm)。對于壁面綠化覆蓋的流場區(qū)域,考慮其空氣動力學(xué)效應(yīng)(對氣流運動的阻礙效應(yīng))和熱力學(xué)效應(yīng)(包含綠化的蒸騰降溫效應(yīng)以及對太陽輻射的遮蔽而帶來的溫度下降),研究中考慮4 種典型的壁面綠化的葉面積密度LAD=0.5,1.0,1.5,2.0;地面機動車尾氣污染源位于峽谷中央,其尺寸為0.6H×0.02H×10H,考慮嚴峻交通污染情況,污染源均勻釋放CO,NO 以及NO2的速率分別為7.37×10-7,5.98×10-8,9.17×10-9kg/(m3·s)[10]。

    圖1 具有壁面綠化的街道峽谷模型Fig.1 Model of the street canyon with green walls

    在研究中,使用設(shè)置壁面溫度來實現(xiàn)夏季不同時刻太陽輻射所導(dǎo)致的壁面升溫[11-13],包括背風(fēng)面受熱、迎風(fēng)面受熱、地面受熱和全壁面受熱(也稱夜間墻體保溫或放熱)。需要指出的是環(huán)境溫度設(shè)為300 K。本文研究的工況共計21 個:1 個既無壁面受熱又無壁面綠化的空白對照工況、4 個僅有壁面受熱而無壁面綠化的工況,以及16 個既有壁面受熱又有壁面綠化的工況(計及4 種壁面綠化的葉面積密度和4 種不同壁面受熱模式)。研究工況的壁面受熱溫度設(shè)置如表1 所示[14]。

    表1 研究工況的壁面溫度設(shè)置Tab.1 Wall temperature setting for the studied cases

    此外,研究還考慮了NOx-O3光化學(xué)反應(yīng)動力學(xué),用來實現(xiàn)反應(yīng)性污染物的光化學(xué)過程。具體的化學(xué)方程式為[15]

    式中:hv為光照強度;M為第3 種分子。

    式中:[O3]為O3的濃度;[O3]b為O3的背景濃度。

    1.2 數(shù)值模型

    1.2.1 控制方程

    街道峽谷內(nèi)的湍流、熱傳遞以及污染物的擴散和轉(zhuǎn)化采用標準k-ε湍流模型進行求解,控制方程為[16]

    式中:xi和xj為笛卡爾坐標;ui和uj分別為xi和xj方向上的空氣時均速度分量;ρ為空氣密度;p為空氣壓強;μ為分子黏性系數(shù);μt為湍流黏性系數(shù);gi為重力加速度分量;T為溫度;Pr為分子普朗特數(shù),Pr=0.85;σt為湍流普朗特數(shù);ST為溫度源項。

    此外,在動量方程求解時,采用Boussinesq 假定[17]來近似處理由溫度變化引起的峽谷內(nèi)空氣浮升力項。

    各污染物的組分輸運方程為[10]

    式中:[X]代表物質(zhì)X 的濃度;Dm為污染物的分子擴散系數(shù);Dt為湍流擴散系數(shù);SX代表物質(zhì)X 的源項;為NO2的光解速率;k1,k2分別為式(2)和式(3)所代表的反應(yīng)速率,其與溫度和污染物濃度呈正相關(guān)[15]。

    壁面綠化覆蓋區(qū)域?qū)е驴刂品匠淘黾拥脑错椃謩e為[18]

    式中:SE,,Sk,Sε分別為能量源項、動量源項、湍流能k的源項和湍動能耗散率 ε的源項。Pc為單位LAD植被的體積冷卻功率;模型常數(shù){Cd,βp,βd,Cε4,Cε5}={0.2,1.0,5.1,0.9,0.9}[19]。

    1.2.2 邊界條件

    入口邊界條件采用速度入口(velocity-inlet),風(fēng)速U(z)、湍動能k(z)和湍動能耗散率 ε(z)的定義為[19]

    式中:Uref表示高度為Href處的參考風(fēng)速,Uref=3.0 m/s,此處Href=10 m,為建筑物高度;α為風(fēng)剖指數(shù),反映較密集建筑群的城市郊區(qū)狀況[13],α=0.187;u*為摩擦速度,u*=0.54 m/s;Cμ為模型常數(shù),Cμ=0.09;κ 為馮卡門常數(shù),κ=0.4;δ為邊界層高度,δ=960 m。此外,入口處還包含體積分數(shù)為2×10-8的O3與體積分數(shù)為0.21 的O2。

    出口邊界選擇自由出流邊界條件;側(cè)面和頂部邊界選擇對稱邊界條件;壁面選擇無滑移壁面邊界條件。

    1.2.3 求解設(shè)置

    控制方程采用有限體積法(FVM)進行離散;對流離散項選用二階迎風(fēng)差分格式;數(shù)值算法選用壓力耦合方程組的半隱式方法(SIMPLE)。求解過程中,當(dāng)各變量的相對計算殘差均降至10-6以下時即視為收斂。

    街道峽谷中尾氣污染物發(fā)生光化學(xué)反應(yīng)需要合適的氣象條件,為了盡可能地接近現(xiàn)實情況,將數(shù)值模擬分為3 步進行[10]:

    第1 步保持關(guān)閉污染源與體積反應(yīng)模塊,先進行計算,得到穩(wěn)定流場。此時,計算域中充滿體積分數(shù)為2×10-8的O3與體積分數(shù)為0.21 的O2。

    第2 步在第1 步的計算結(jié)果的基礎(chǔ)上,開啟污染源,不開啟體積反應(yīng)模塊。通過此步驟能夠得到無光化學(xué)反應(yīng)條件下穩(wěn)定的污染物分布,為計算后續(xù)的光化學(xué)反應(yīng)提供初始條件。

    第3 步保持污染源開啟,并開啟體積反應(yīng)。最終得到化學(xué)平衡時各污染物的分布情況。

    1.3 模型驗證

    采用德國卡爾斯魯厄大學(xué)風(fēng)洞實驗數(shù)據(jù)對CFD 模型進行驗證和完善,最終確定使用標準k-ε湍流模型進行后續(xù)的工況計算。具體的驗證細節(jié)參見文獻[20]。

    2 結(jié)果與討論

    2.1 街谷內(nèi)流場變化分析

    圖2 為4 種壁面受熱和2 種壁面綠化(LAD分別為1.0 和2.0)配置下街道峽谷豎直中心面上的無量綱速度云圖及流線分布。作為空白對照,圖2也給出既無壁面受熱又無壁面綠化下的街道峽谷豎直中心面處的相關(guān)流場信息。

    在既無壁面綠化又無壁面受熱條件下,從圖2(a)可以看出:在峽谷豎直中心面的背風(fēng)側(cè)處形成一股(單束)上升氣流,該上升氣流流出峽谷與峽谷頂部自由來流匯合;另一股氣流則于峽谷迎風(fēng)面附近形成一個順時針旋渦。

    無論有無壁面綠化,背風(fēng)面受熱所產(chǎn)生的熱浮力會驅(qū)散背風(fēng)側(cè)的單束向上運動的氣流,從而在峽谷垂直中心面上形成一個順時針大旋渦(圖2(b)、圖2(c)和圖2(d)),該大旋渦的外沿氣流從背風(fēng)側(cè)流出峽谷而與峽谷上部自由來流匯集。在背風(fēng)面受熱時,隨著壁面綠化的LAD的增加,峽谷背風(fēng)面底部的低速區(qū)域增大,而當(dāng)LAD增加至2.0時,背風(fēng)墻底部的低速區(qū)域演化成一個風(fēng)速極低的逆時針小角渦。

    圖2 街道峽谷豎直中心面上的無量綱速度云圖及流線Fig.2 Contours of dimensionless velocity and streamlines at the vertical center plane of the street canyon

    相比無壁面綠化及無壁面受熱工況(圖2(a)),在無壁面綠化條件下,迎風(fēng)面受熱會極大削弱迎風(fēng)側(cè)的順時針旋渦,從而間接增強背風(fēng)側(cè)的單束向上氣流:迎風(fēng)面受熱產(chǎn)生的熱浮力方向與旋渦氣流在近迎風(fēng)面的流動方向相反,故會極大削弱順時針旋渦尺度,并于迎風(fēng)面附近形成一塊低速區(qū)域(圖2(e))。當(dāng)有壁面綠化的峽谷迎風(fēng)面受熱時,隨著壁面綠化的LAD的增大,其蒸騰降溫效應(yīng)隨之增強,迎風(fēng)面受熱對流場的影響減小,故迎風(fēng)側(cè)順時針旋渦尺度得以恢復(fù)增大(圖2(e)、圖2(f)和圖2(g))。

    相比無壁面綠化及無壁面受熱工況(圖2(a)),地面受熱也能驅(qū)散背風(fēng)側(cè)的單束上升氣流,使得順時針旋渦進一步得以發(fā)展(圖2(h),2(i)和2(j))。很明顯,在地面受熱條件下,壁面綠化有助于順時針旋渦的形成(比較圖2(b)與2(c)及2(d))。同時得益于更大的受熱面積,地面受熱有助于增加峽谷內(nèi)部的平均氣流速度,增強峽谷通風(fēng)。

    全壁面受熱則是受熱面積最大的受熱模式,在該種情況下,三壁面均產(chǎn)生強烈的向上熱浮升力,使得峽谷通道內(nèi)的橫向氣流均向中心面處匯集,并向上運動流出峽谷(圖2(k))。當(dāng)具有壁面綠化的全壁面受熱時,隨著壁面綠化的LAD上升至1.0,壁面受熱產(chǎn)生的熱浮升力也被削弱,迎風(fēng)面附近的順時針旋渦也初現(xiàn)雛形(圖2(l)),而隨著LAD進一步增加至2.0 時,該旋渦尺度略微增大(圖2(m))。

    2.2 街谷內(nèi)溫度變化分析

    圖3 為不同工況條件下街道峽谷豎直中心面處的溫度分布。圖4 為不同工況條件下街道峽谷內(nèi)部平均溫度。如圖3(a)所示,在無壁面綠化條件下,背風(fēng)面受熱使得背風(fēng)面附加空氣溫度升高至302 K,且其周圍有較為明顯的高溫區(qū)域,而順時針主旋渦則能夠攜帶壁面周圍的熱量于峽谷中循環(huán),提高峽谷內(nèi)部的平均溫度。從圖4 可知,在無壁面綠化條件下,背風(fēng)面受熱能夠提升峽谷內(nèi)部平均溫度約0.35 K(相對于空白對照工況)。圖3(b)顯示壁面綠化的蒸騰降溫效應(yīng)能夠有效降低兩側(cè)壁面周圍的溫度,由于密度差異原因,被降溫的空氣于近壁面處向下沉積,聚集在街道峽谷底部并得以停留,故可以在峽谷左右下角發(fā)現(xiàn)較為明顯的低溫區(qū)域。進一步,從圖4 可知,在LAD=2.0 的情況下,壁面綠化能夠降低峽谷內(nèi)部平均溫度約0.38 K(相對于無綠化的背風(fēng)面受熱)。

    圖3 街道峽谷豎直中心面處的溫度分布Fig.3 Distribution of temperature at vertical center plane of street canyon

    圖4 街道峽谷內(nèi)部平均溫度Fig.4 Averaged temperature within the street canyon

    如圖3(c)所示,在無壁面綠化條件下,迎風(fēng)面受熱使其周圍仍然存在較為明顯的高溫區(qū)域,順時針小旋渦的運動致使大量熱量在此處循環(huán),難以擴散出峽谷。迎風(fēng)面受熱能夠提高街道峽谷內(nèi)部的平均溫度約0.60 K(相對于空白對照工況)。圖3(d)顯示壁面綠化能夠有效降低兩側(cè)壁面周圍的溫度。從圖4 可知,迎風(fēng)面受熱時,壁面綠化能夠降低峽谷內(nèi)部平均溫度約0.40 K(相對于無綠化的迎風(fēng)面受熱)。

    如圖3(e)所示,在無壁面綠化條件下,地面受熱使街道峽谷內(nèi)部的有效受熱面積更大,相較于前兩種受熱方式,峽谷內(nèi)部平均溫度有明顯升高,約0.90 K(相對于空白對照工況)。在該系列工況中,氣流會攜帶近地面的熱空氣向背風(fēng)面運動,隨后再沿入口氣流方向偏轉(zhuǎn)并流出峽谷。這個過程使得大量熱量被運輸至峽谷左下角及近背風(fēng)面處,使得該區(qū)域溫度遠高于迎風(fēng)面周圍區(qū)域,且最高溫度達到了約305 K。從圖4 可知,在該種受熱情況下,壁面綠化對街道峽谷內(nèi)部熱環(huán)境改善的能力十分有限,僅降低內(nèi)部平均溫度約0.32 K(相對于無綠化的地面受熱),這在4 種受熱模式中為最低值(圖3(f))。

    如圖3(g)所示,在全壁面受熱情況下,峽谷內(nèi)部溫度進一步升高,高溫區(qū)由三處壁面向峽谷內(nèi)部擴散,也大幅度提高了峽谷上空的溫度。由于峽谷內(nèi)向上運動的氣流依然向背風(fēng)面偏轉(zhuǎn),故溫度最高的區(qū)域仍然是峽谷的左下角部分。從圖4可知,無壁面綠化的全壁面受熱時峽谷內(nèi)部平均溫度提升也最大,約為1.43 K(相對于空白對照工況);如圖3(h)顯示,有綠化的全壁面受熱時,峽谷內(nèi)溫度分布與無綠化條件的全壁面受熱情況下的峽谷內(nèi)溫度分布式樣相似,但從圖4 可知,壁面綠化的存在則能夠均勻地降低峽谷內(nèi)部的溫度約0.39 K(相對于無綠化的全壁面受熱)。

    2.3 街谷內(nèi)CO 擴散與NOx-O3 反應(yīng)情況分析

    本研究涉及3 種類型的污染物,分別為CO,NO 和NO2。此處以CO 為例(不涉及反應(yīng)),圖5所示為各種工況下街道峽谷豎直中心面的CO 體積分數(shù)分布。在既無壁面綠化又無壁面受熱條件下(空白對照工況),峽谷內(nèi)部的流動結(jié)構(gòu)由單束氣流及順時針旋渦共同主導(dǎo),在兩束氣流交匯的區(qū)域污染物大量堆積,并朝背風(fēng)面方向向上逐漸遞減,使得峽谷背風(fēng)面區(qū)域的污染物濃度遠高于迎風(fēng)面區(qū)域的污染物濃度(如圖5(a)所示)。

    背風(fēng)面受熱會增強順時針旋渦,并增強上升氣流速度,使得峽谷內(nèi)部流動由順時針單旋渦結(jié)構(gòu)主導(dǎo)。近地面污染源釋放的CO 隨順時針旋渦在峽谷中運動,沿背風(fēng)面向上遞減分布,且背風(fēng)面附近的污染物濃度遠高于迎風(fēng)面附近(如圖5(b)、圖5(c)和圖5(d)所示)。隨著LAD的增加,迎風(fēng)面附近的污染物濃度明顯降低,大量污染物被輸送至峽谷的背風(fēng)面附近;在LAD=2.0 的工況中,在旋渦氣流的控制下,在背風(fēng)側(cè)角渦處形成一塊較小的高污染區(qū)域。

    在迎風(fēng)面受熱情況下,近地面污染物CO 受單束氣流影響,濃度沿水平方向遞減分布,且污染物大多聚集于峽谷中段區(qū)域(如圖5(e),5(f)和5(g)所示);而迎風(fēng)側(cè)小旋渦則能夠?qū)⒂L(fēng)面附近污染物向背風(fēng)面方向輸送,使得迎風(fēng)面附近污染物濃度偏低,導(dǎo)致污染物分布在水平方向發(fā)生略微偏移。在該種情況下,由于大量污染物仍集中于街道峽谷中心,故行人呼吸面上平均污染物CO 濃度也略微提高。隨著LAD的上升,迎風(fēng)側(cè)小型旋渦規(guī)模逐漸增大,逐漸恢復(fù)對峽谷內(nèi)部流動結(jié)構(gòu)的控制。

    地面受熱情況下峽谷內(nèi)部的流動結(jié)構(gòu)與背風(fēng)面受熱極為類似,故地面受熱情況下峽谷內(nèi)部污染物的分布規(guī)律(如圖5(h),5(i)和5(j)所示)也與背風(fēng)面受熱情況下峽谷內(nèi)部污染物的分布極為相似。

    圖5 街道峽谷豎直中心面處CO 體積分數(shù)分布Fig.5 Distributions of CO volume fraction at vertical center plane of street canyon

    如圖5(k),5(l)和5(m)所示,在全壁面受熱情況下,兩部分氣流聚集于峽谷中心并向上涌出,橫向氣流攜帶了大量污染物向峽谷中心運動,使得中心面上的污染物濃度極高。大量污染物積聚于峽谷的背風(fēng)側(cè),并向上遞減。隨著LAD的上升,壁面綠化的蒸騰降溫效應(yīng)使得峽谷內(nèi)部氣流得以穩(wěn)定,重新回歸單旋渦結(jié)構(gòu),污染物的分布也轉(zhuǎn)化為沿背風(fēng)面向上遞減分布的結(jié)構(gòu)(如圖5(m)所示)。

    需要指出的是,未涉及反應(yīng)時,NO 和NO2的分布特征與CO 基本相同,且濃度比為10∶1。此處僅列舉在既無壁面綠化又無壁面受熱條件下NO 與NO2前后體積分數(shù)對比與全部工況下的臭氧消耗率。圖6 為反應(yīng)前后NO 和NO2體積分數(shù)變化(無壁面受熱及無綠化條件下)。由圖6 可知,環(huán)境中反應(yīng)后NO 體積分數(shù)略微降低而NO2體積分數(shù)升高,伴隨環(huán)境中背景臭氧的消耗,說明反應(yīng)式(3)在街道峽谷中占據(jù)主導(dǎo)地位。

    圖6 反應(yīng)前后NO 和NO2 體積分數(shù)變化(無壁面受熱及無綠化)Fig.6 Changes of volume fraction of NO and NO2 before and after reactions

    背風(fēng)面受熱改變了峽谷內(nèi)部的流動結(jié)構(gòu),從而使得污染物大部分集中在背風(fēng)面附近,且沿其向上遞減分布,也呈現(xiàn)出類似的趨勢:高區(qū)域均集中于峽谷背風(fēng)面附近,并朝向峽谷外部依次遞減(如圖7(b),7(c)和7(d)所示)。隨著LAD的增加,近迎風(fēng)面污染物濃度降低而近背風(fēng)面則升高,=60%的等值線從迎風(fēng)面附近收縮至峽谷中心處。進一步分析表2,發(fā)現(xiàn)LAD=2.0 的壁面綠化使得峽谷內(nèi)部平均下降了約5.02%(相對于無壁面綠化的背風(fēng)面受熱工況),這歸因于其蒸騰降溫效應(yīng)極大地削弱了兩側(cè)壁面附近的化學(xué)反應(yīng)強度。在這種情況下,較低的背景臭氧濃度成為了達到化學(xué)平衡的限制因素。

    如圖7(e),7(f)和7(g)所示,迎風(fēng)面受熱條件下,高區(qū)域均集中于峽谷底部中段,即污染源所在位置,并向上遞減,峽谷內(nèi)部的幾乎均大于70%。隨著LAD增加至2.0,=70%的區(qū)域也從迎風(fēng)面底部收縮至峽谷中段區(qū)域,峽谷內(nèi)部平均也下降至69.19%,較無綠化工況下降了3.63%,如表2 所示。

    如圖7(h),7(i)和7(j)所示,在地面受熱條件下,分布規(guī)律幾乎與背風(fēng)面受熱情形相同。從表2 可知,地面受熱情況下,LAD=2.0 的壁面綠化工況相較于無綠化的受熱工況,峽谷內(nèi)部平均下降約3.59%。

    如圖7(k),7(l)和7(m)所示,全壁面受熱能夠大幅度提高峽谷內(nèi)部的平均溫度,故峽谷內(nèi)部平均也普遍較高,均大于70%。在無綠化的條件下,橫向氣流使得大量污染物聚集于峽谷左半部分區(qū)域,使得該處有一塊較大的高區(qū)域。隨著LAD的增加,峽谷內(nèi)部的污染物濃度也轉(zhuǎn)變?yōu)檠乇筹L(fēng)面向上遞減分布,此時,分布特點也與背風(fēng)面受熱和地面受熱工況相似,但是,由于峽谷內(nèi)部溫度普遍較高,峽谷內(nèi)部平均也高于上述兩類工況。

    圖7 街道峽谷豎直中心面處分布Fig.7 Distributions of at vertical center plane of the street canyon

    進一步,表2 表明,在每一種壁面受熱模式下,峽谷內(nèi)部平均隨LAD的增大而降低。

    表2 街道峽谷內(nèi)部平均Tab.2 Average within street canyon

    表2 街道峽谷內(nèi)部平均Tab.2 Average within street canyon

    綜上所述,在AR=0.5 的淺型街道峽谷中,壁面綠化的蒸發(fā)降溫效應(yīng)可以有效降低峽谷內(nèi)部的溫度,并緩解峽谷內(nèi)光化學(xué)反應(yīng)的劇烈程度。

    3 結(jié)論

    建立了耦合壁面綠化及熱效應(yīng)的淺型街道峽谷(AR=0.5)內(nèi)污染物擴散與轉(zhuǎn)化的CFD 模型,并利用風(fēng)洞實驗數(shù)據(jù)對選用的標準k-ε湍流模型及相關(guān)模型參數(shù)進行驗證。采用CFD 模型對壁面綠化及熱效應(yīng)影響下的淺型街道峽谷內(nèi)流場、污染物擴散及NOx-O3光化學(xué)轉(zhuǎn)化開展模擬研究。研究中考慮4 種壁面綠化的葉面積密度(LAD=0.5,1.0,1.5,2.0)以及4 種壁面加熱模式(背風(fēng)面加熱、迎風(fēng)面加熱、地面加熱和全壁面加熱)。得到如下主要結(jié)論:

    a.在既無壁面綠化又無壁面加熱條件下,峽谷垂直中心面上的流動由橫向氣流和順時針旋渦共同主導(dǎo)。在無壁面綠化時,迎風(fēng)面受熱會削弱順時針旋渦的尺度,從而使得橫向氣流占據(jù)主導(dǎo)地位;背風(fēng)面受熱和地面受熱會強化順時針旋渦;而全壁面受熱則會使旋渦結(jié)構(gòu)消散。

    b.每種壁面受熱模式下,壁面綠化對峽谷內(nèi)部流場的影響隨LAD的增大而加?。ㄌ貏e對于迎風(fēng)面受熱和全壁面受熱工況尤為如此)。

    c.壁面受熱導(dǎo)致峽谷內(nèi)的平均溫度升高程度為:全壁面受熱最高,其次為地面受熱、迎風(fēng)面受熱和背風(fēng)面受熱。相較于既無壁面綠化又無壁面加熱工況,根據(jù)不同的壁面受熱模式,無壁面綠化的壁面加熱可使峽谷內(nèi)平均氣溫升高0.35~1.43 K;壁面綠化能夠有效改善峽谷內(nèi)的熱環(huán)境,根據(jù)不同的壁面受熱模式,LAD=2.0 的壁面綠化可使峽谷內(nèi)的平均氣溫降低0.38~0.48 K。

    d.壁面受熱會影響污染物的擴散,尤其是迎風(fēng)面受熱會促使大量污染物聚集在峽谷底部中心位置附近。

    猜你喜歡
    背風(fēng)面旋渦順時針
    為什么鐘表順時針轉(zhuǎn)?
    最后才吃梨
    小心,旋渦來啦
    大班科學(xué)活動:神秘的旋渦
    旋渦笑臉
    山間湖
    揚子江(2019年1期)2019-03-08 02:52:34
    為什么表的指針都按照順時針方向轉(zhuǎn)動
    非均勻等離子體Ka-Band傳輸性能中繼法優(yōu)化研究
    高超聲速風(fēng)洞子母彈大迎角拋殼投放試驗
    高壓輸電鐵塔塔身背風(fēng)面風(fēng)荷載遮擋效應(yīng)研究
    麻豆一二三区av精品| 国产熟女xx| 欧美不卡视频在线免费观看| 日韩成人在线观看一区二区三区| 天天添夜夜摸| 97超级碰碰碰精品色视频在线观看| 亚洲人成网站高清观看| av天堂中文字幕网| 久久精品国产自在天天线| 精品久久久久久久毛片微露脸| 免费观看的影片在线观看| 国产成人啪精品午夜网站| 中文字幕av在线有码专区| 精品国产超薄肉色丝袜足j| 草草在线视频免费看| 性色av乱码一区二区三区2| 无人区码免费观看不卡| 国产亚洲精品综合一区在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲av电影在线进入| tocl精华| 一级毛片女人18水好多| 老熟妇仑乱视频hdxx| 黄片小视频在线播放| 老鸭窝网址在线观看| 亚洲欧美精品综合久久99| 国产亚洲精品av在线| 一进一出抽搐gif免费好疼| 国产精品免费一区二区三区在线| 俺也久久电影网| 国产单亲对白刺激| 欧美日韩瑟瑟在线播放| 午夜免费激情av| 一卡2卡三卡四卡精品乱码亚洲| 国产伦在线观看视频一区| 91久久精品电影网| АⅤ资源中文在线天堂| 国产亚洲av嫩草精品影院| 日本在线视频免费播放| 波野结衣二区三区在线 | 综合色av麻豆| 亚洲人成伊人成综合网2020| 看片在线看免费视频| 国内精品美女久久久久久| 一进一出好大好爽视频| 久久九九热精品免费| 精品久久久久久久久久久久久| 成人国产综合亚洲| 亚洲自拍偷在线| 亚洲狠狠婷婷综合久久图片| 很黄的视频免费| 成人一区二区视频在线观看| 久久久久久人人人人人| 国产精品爽爽va在线观看网站| 99精品久久久久人妻精品| 亚洲在线自拍视频| 麻豆久久精品国产亚洲av| 亚洲,欧美精品.| 成人无遮挡网站| 日本在线视频免费播放| 亚洲中文字幕一区二区三区有码在线看| 久久精品国产清高在天天线| 亚洲av免费高清在线观看| 制服人妻中文乱码| 精品国产超薄肉色丝袜足j| 亚洲av免费在线观看| 午夜老司机福利剧场| 在线播放国产精品三级| 不卡一级毛片| 90打野战视频偷拍视频| 97超级碰碰碰精品色视频在线观看| 韩国av一区二区三区四区| 神马国产精品三级电影在线观看| 天美传媒精品一区二区| 色综合欧美亚洲国产小说| 亚洲国产日韩欧美精品在线观看 | 久久久久国内视频| 欧美中文综合在线视频| 亚洲国产高清在线一区二区三| 欧美午夜高清在线| av在线天堂中文字幕| 午夜老司机福利剧场| 一卡2卡三卡四卡精品乱码亚洲| 国产又黄又爽又无遮挡在线| 三级国产精品欧美在线观看| 综合色av麻豆| 国产精品久久电影中文字幕| 成人特级黄色片久久久久久久| 高清日韩中文字幕在线| 精品一区二区三区视频在线 | 国产一区二区在线av高清观看| 亚洲精品美女久久久久99蜜臀| 亚洲精品在线美女| 日韩欧美精品v在线| АⅤ资源中文在线天堂| 亚洲五月天丁香| 一进一出抽搐动态| 最新美女视频免费是黄的| 18+在线观看网站| 日韩国内少妇激情av| 成年女人看的毛片在线观看| av天堂在线播放| 亚洲精品美女久久久久99蜜臀| 亚洲精华国产精华精| 精品熟女少妇八av免费久了| 国产精品久久电影中文字幕| 亚洲国产欧美网| 动漫黄色视频在线观看| 国产精品影院久久| 精品乱码久久久久久99久播| 综合色av麻豆| 女人十人毛片免费观看3o分钟| 国产麻豆成人av免费视频| 少妇的逼水好多| 女人高潮潮喷娇喘18禁视频| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩亚洲国产一区二区在线观看| 动漫黄色视频在线观看| 久久久久久大精品| 蜜桃亚洲精品一区二区三区| 国产精品亚洲一级av第二区| 美女高潮喷水抽搐中文字幕| 亚洲成人久久爱视频| 最近最新中文字幕大全电影3| 亚洲成人免费电影在线观看| 亚洲无线在线观看| 女警被强在线播放| 最近最新中文字幕大全电影3| 国产精品电影一区二区三区| 91久久精品国产一区二区成人 | 老鸭窝网址在线观看| 久久久久久久精品吃奶| 国产亚洲av嫩草精品影院| 亚洲欧美精品综合久久99| 人人妻人人看人人澡| 一级黄色大片毛片| 精品一区二区三区视频在线观看免费| 2021天堂中文幕一二区在线观| 一a级毛片在线观看| 国产亚洲欧美在线一区二区| 亚洲国产精品成人综合色| 欧美色视频一区免费| 久久久国产成人精品二区| 午夜精品一区二区三区免费看| 在线十欧美十亚洲十日本专区| 91久久精品电影网| 男人舔女人下体高潮全视频| 亚洲内射少妇av| 嫩草影院入口| 啦啦啦观看免费观看视频高清| 日韩免费av在线播放| 免费观看人在逋| 久久伊人香网站| 国产亚洲精品久久久久久毛片| 国产蜜桃级精品一区二区三区| 99久久精品一区二区三区| 少妇高潮的动态图| 婷婷精品国产亚洲av在线| 久久久久久久久大av| 色播亚洲综合网| 女人被狂操c到高潮| 窝窝影院91人妻| 在线观看免费视频日本深夜| 国产高清激情床上av| 少妇裸体淫交视频免费看高清| 国产视频内射| 18+在线观看网站| 亚洲精品一区av在线观看| 五月玫瑰六月丁香| 久久久久九九精品影院| 国产久久久一区二区三区| 天天一区二区日本电影三级| 国产真人三级小视频在线观看| 九色成人免费人妻av| 亚洲av免费高清在线观看| 在线播放国产精品三级| 亚洲欧美日韩东京热| 国产男靠女视频免费网站| 久久精品国产清高在天天线| 亚洲av美国av| 露出奶头的视频| 欧美性感艳星| 亚洲成人久久爱视频| 99热这里只有是精品50| 午夜两性在线视频| 国产不卡一卡二| 91在线精品国自产拍蜜月 | 一个人看视频在线观看www免费 | 久99久视频精品免费| 国产又黄又爽又无遮挡在线| av福利片在线观看| 亚洲乱码一区二区免费版| 亚洲av成人av| 亚洲av美国av| 精华霜和精华液先用哪个| 最好的美女福利视频网| 国产免费一级a男人的天堂| av欧美777| 国产亚洲欧美在线一区二区| 一本一本综合久久| av黄色大香蕉| 熟妇人妻久久中文字幕3abv| 国产精品98久久久久久宅男小说| 久久国产乱子伦精品免费另类| 88av欧美| 一级黄片播放器| 午夜日韩欧美国产| 男女午夜视频在线观看| 免费看美女性在线毛片视频| 中文字幕精品亚洲无线码一区| 中文字幕人妻熟人妻熟丝袜美 | 亚洲熟妇中文字幕五十中出| 午夜老司机福利剧场| 亚洲av成人精品一区久久| 悠悠久久av| 老汉色av国产亚洲站长工具| 欧美色视频一区免费| 久久久久久久亚洲中文字幕 | 757午夜福利合集在线观看| 美女免费视频网站| 国产精品女同一区二区软件 | 久久香蕉国产精品| 最近在线观看免费完整版| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 黄色片一级片一级黄色片| 亚洲avbb在线观看| 99国产精品一区二区三区| 精品电影一区二区在线| 亚洲内射少妇av| 别揉我奶头~嗯~啊~动态视频| 国产在线精品亚洲第一网站| 亚洲中文字幕一区二区三区有码在线看| 九九久久精品国产亚洲av麻豆| 国产精品久久久久久人妻精品电影| 人妻久久中文字幕网| 国内精品一区二区在线观看| 一区福利在线观看| 啪啪无遮挡十八禁网站| 亚洲久久久久久中文字幕| 国产精品一区二区三区四区免费观看 | 欧美日韩瑟瑟在线播放| 床上黄色一级片| 亚洲久久久久久中文字幕| 精品电影一区二区在线| 成年人黄色毛片网站| 特级一级黄色大片| 国产精品日韩av在线免费观看| 日本三级黄在线观看| 我的老师免费观看完整版| eeuss影院久久| 国产毛片a区久久久久| 老熟妇仑乱视频hdxx| 免费电影在线观看免费观看| 18禁美女被吸乳视频| 久久久久九九精品影院| 蜜桃久久精品国产亚洲av| 国产视频一区二区在线看| 香蕉av资源在线| 青草久久国产| 亚洲欧美日韩高清在线视频| 老熟妇仑乱视频hdxx| or卡值多少钱| 日本撒尿小便嘘嘘汇集6| 国产亚洲精品一区二区www| 久久国产乱子伦精品免费另类| 狠狠狠狠99中文字幕| 欧美极品一区二区三区四区| 国产精品一区二区三区四区免费观看 | 成人av一区二区三区在线看| 变态另类丝袜制服| 很黄的视频免费| 一个人看视频在线观看www免费 | 欧美国产日韩亚洲一区| 国内精品美女久久久久久| 免费高清视频大片| 桃色一区二区三区在线观看| 成人国产综合亚洲| 国产伦精品一区二区三区视频9 | 国产精品亚洲av一区麻豆| 丁香六月欧美| 最好的美女福利视频网| 99国产精品一区二区蜜桃av| 国产av一区在线观看免费| 熟女电影av网| 国产精品三级大全| 欧美丝袜亚洲另类 | 黄片大片在线免费观看| 中文字幕久久专区| 99国产极品粉嫩在线观看| 一个人免费在线观看电影| 蜜桃久久精品国产亚洲av| 夜夜看夜夜爽夜夜摸| 麻豆国产97在线/欧美| 久久草成人影院| 国产高清视频在线播放一区| 亚洲欧美日韩高清专用| 国产午夜福利久久久久久| 久久亚洲真实| 欧美日韩瑟瑟在线播放| 亚洲av第一区精品v没综合| a级一级毛片免费在线观看| 黄片小视频在线播放| 熟女少妇亚洲综合色aaa.| 亚洲真实伦在线观看| 好男人在线观看高清免费视频| 久久国产精品影院| 国产精品久久视频播放| 亚洲精品一卡2卡三卡4卡5卡| 天美传媒精品一区二区| 又黄又粗又硬又大视频| 亚洲av一区综合| 美女被艹到高潮喷水动态| 欧美午夜高清在线| 日本在线视频免费播放| 九九久久精品国产亚洲av麻豆| 99久国产av精品| 亚洲精品一区av在线观看| 国产精品美女特级片免费视频播放器| 亚洲一区二区三区色噜噜| 老鸭窝网址在线观看| 免费看光身美女| 国产亚洲欧美在线一区二区| 法律面前人人平等表现在哪些方面| 久久久久久人人人人人| 精品熟女少妇八av免费久了| 日韩免费av在线播放| 久久久色成人| 精品一区二区三区人妻视频| 午夜福利视频1000在线观看| 成年女人看的毛片在线观看| 日韩欧美精品免费久久 | 婷婷六月久久综合丁香| svipshipincom国产片| 3wmmmm亚洲av在线观看| 亚洲avbb在线观看| 草草在线视频免费看| 一卡2卡三卡四卡精品乱码亚洲| 日韩成人在线观看一区二区三区| 久久精品国产99精品国产亚洲性色| 久久久久亚洲av毛片大全| 成人av在线播放网站| 国产99白浆流出| 夜夜爽天天搞| 手机成人av网站| 欧美成人性av电影在线观看| 伊人久久精品亚洲午夜| 欧美激情在线99| 成人永久免费在线观看视频| 国产精品爽爽va在线观看网站| 国产伦精品一区二区三区四那| 中文字幕人妻丝袜一区二区| 啦啦啦韩国在线观看视频| 久久中文看片网| 中文字幕人妻熟人妻熟丝袜美 | a级一级毛片免费在线观看| 长腿黑丝高跟| 精品久久久久久久人妻蜜臀av| 亚洲国产欧洲综合997久久,| 精品一区二区三区av网在线观看| 五月玫瑰六月丁香| 窝窝影院91人妻| 少妇的逼好多水| 国内精品久久久久精免费| 亚洲,欧美精品.| 九色国产91popny在线| 美女黄网站色视频| 怎么达到女性高潮| 亚洲精品国产精品久久久不卡| 99精品在免费线老司机午夜| 国产黄a三级三级三级人| 啦啦啦韩国在线观看视频| 亚洲熟妇熟女久久| 午夜视频国产福利| 国产精品女同一区二区软件 | 成人欧美大片| 国产三级黄色录像| 亚洲中文日韩欧美视频| 最近视频中文字幕2019在线8| 小蜜桃在线观看免费完整版高清| 天天添夜夜摸| 伊人久久精品亚洲午夜| 男女视频在线观看网站免费| 欧美日韩综合久久久久久 | 夜夜爽天天搞| 久久久久精品国产欧美久久久| 99视频精品全部免费 在线| 长腿黑丝高跟| 欧美3d第一页| 亚洲国产精品合色在线| 国产精品,欧美在线| 久久久久性生活片| 九九在线视频观看精品| 偷拍熟女少妇极品色| 天堂网av新在线| av欧美777| 免费在线观看成人毛片| 全区人妻精品视频| 亚洲av一区综合| 国产激情欧美一区二区| 国产精品香港三级国产av潘金莲| 精品熟女少妇八av免费久了| 国产精品99久久99久久久不卡| 舔av片在线| 国产成人a区在线观看| 一a级毛片在线观看| 国产中年淑女户外野战色| 精品不卡国产一区二区三区| 国产精品精品国产色婷婷| 国产又黄又爽又无遮挡在线| 午夜福利欧美成人| 别揉我奶头~嗯~啊~动态视频| 日韩精品中文字幕看吧| 哪里可以看免费的av片| 亚洲国产中文字幕在线视频| 欧美+亚洲+日韩+国产| 男插女下体视频免费在线播放| 成人永久免费在线观看视频| 亚洲国产欧美人成| 国产aⅴ精品一区二区三区波| 2021天堂中文幕一二区在线观| 最近在线观看免费完整版| 欧美最新免费一区二区三区 | av国产免费在线观看| 女生性感内裤真人,穿戴方法视频| 国产综合懂色| 色吧在线观看| 天堂av国产一区二区熟女人妻| 噜噜噜噜噜久久久久久91| 无限看片的www在线观看| 久久6这里有精品| 在线a可以看的网站| www.www免费av| 动漫黄色视频在线观看| 三级毛片av免费| 日本五十路高清| 国产精品亚洲av一区麻豆| 一a级毛片在线观看| 国产在视频线在精品| 亚洲国产欧洲综合997久久,| 色综合婷婷激情| 天美传媒精品一区二区| 又紧又爽又黄一区二区| 国产精品嫩草影院av在线观看 | 国产激情偷乱视频一区二区| 黄色视频,在线免费观看| 少妇熟女aⅴ在线视频| 欧美xxxx黑人xx丫x性爽| 全区人妻精品视频| 少妇人妻精品综合一区二区 | 黄片小视频在线播放| 免费人成在线观看视频色| 亚洲狠狠婷婷综合久久图片| 少妇高潮的动态图| 亚洲av免费高清在线观看| 成人国产一区最新在线观看| 又爽又黄无遮挡网站| 婷婷六月久久综合丁香| 免费在线观看日本一区| 黄色视频,在线免费观看| 欧美激情久久久久久爽电影| a在线观看视频网站| 老熟妇仑乱视频hdxx| 最新在线观看一区二区三区| 日韩欧美一区二区三区在线观看| 窝窝影院91人妻| 国内毛片毛片毛片毛片毛片| 欧美最新免费一区二区三区 | 欧美黄色淫秽网站| 日韩欧美一区二区三区在线观看| 丁香六月欧美| 69人妻影院| 99视频精品全部免费 在线| 在线观看免费午夜福利视频| 99国产极品粉嫩在线观看| 又爽又黄无遮挡网站| 亚洲成人精品中文字幕电影| 国产69精品久久久久777片| 性欧美人与动物交配| 色哟哟哟哟哟哟| tocl精华| 亚洲中文字幕一区二区三区有码在线看| 国产精品久久视频播放| 老司机福利观看| 丁香欧美五月| 久久精品国产99精品国产亚洲性色| 伊人久久大香线蕉亚洲五| 日本一二三区视频观看| 又黄又爽又免费观看的视频| 一个人观看的视频www高清免费观看| 亚洲人成网站在线播放欧美日韩| 黄色片一级片一级黄色片| 18禁黄网站禁片免费观看直播| 国产又黄又爽又无遮挡在线| 亚洲狠狠婷婷综合久久图片| 欧美最新免费一区二区三区 | 亚洲最大成人手机在线| 午夜影院日韩av| 久久久久精品国产欧美久久久| 国产一区二区激情短视频| 无限看片的www在线观看| 最后的刺客免费高清国语| 亚洲精品亚洲一区二区| 亚洲国产精品久久男人天堂| 免费高清视频大片| 亚洲最大成人中文| 欧美成狂野欧美在线观看| 国产精品久久久久久亚洲av鲁大| 久久伊人香网站| 亚洲国产欧美人成| 精品福利观看| 国产真实乱freesex| 18禁黄网站禁片免费观看直播| 亚洲人成网站高清观看| 免费观看人在逋| 亚洲中文字幕日韩| 美女高潮的动态| 国产精品,欧美在线| 91在线观看av| 国模一区二区三区四区视频| 此物有八面人人有两片| www.熟女人妻精品国产| 久久久久久久午夜电影| 黄色日韩在线| 可以在线观看的亚洲视频| 色噜噜av男人的天堂激情| 脱女人内裤的视频| 69av精品久久久久久| 日本黄大片高清| 亚洲第一欧美日韩一区二区三区| 十八禁网站免费在线| 亚洲一区二区三区不卡视频| 免费观看的影片在线观看| 久久人人精品亚洲av| 最近最新中文字幕大全电影3| 乱人视频在线观看| 国产精品久久久久久久电影 | 国产欧美日韩一区二区三| 欧美高清成人免费视频www| 欧美xxxx黑人xx丫x性爽| 神马国产精品三级电影在线观看| 色尼玛亚洲综合影院| 久久久久久国产a免费观看| 国产乱人伦免费视频| 国产麻豆成人av免费视频| 草草在线视频免费看| 亚洲中文日韩欧美视频| 日韩欧美精品免费久久 | 我要搜黄色片| 不卡一级毛片| 人妻久久中文字幕网| 床上黄色一级片| 网址你懂的国产日韩在线| 少妇人妻精品综合一区二区 | 久久久久久国产a免费观看| 一个人看的www免费观看视频| 狂野欧美白嫩少妇大欣赏| 国产蜜桃级精品一区二区三区| 免费av观看视频| 午夜免费激情av| 在线视频色国产色| 国产探花在线观看一区二区| 免费搜索国产男女视频| 亚洲精品亚洲一区二区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一个人免费在线观看电影| 亚洲一区二区三区色噜噜| 成人亚洲精品av一区二区| 久久精品综合一区二区三区| 波野结衣二区三区在线 | 一区二区三区高清视频在线| 内地一区二区视频在线| 舔av片在线| 国产午夜福利久久久久久| 精品日产1卡2卡| 19禁男女啪啪无遮挡网站| 国产精品98久久久久久宅男小说| 老司机午夜十八禁免费视频| 88av欧美| 99久久九九国产精品国产免费| 国产精品美女特级片免费视频播放器| 色综合亚洲欧美另类图片| 久久久久免费精品人妻一区二区| 久久久色成人| 久久亚洲精品不卡| 91九色精品人成在线观看| 男人和女人高潮做爰伦理| 午夜a级毛片| 小蜜桃在线观看免费完整版高清| 18禁裸乳无遮挡免费网站照片| 成人三级黄色视频| 免费一级毛片在线播放高清视频| 日韩欧美精品免费久久 | 两个人视频免费观看高清| 亚洲专区中文字幕在线| 一本精品99久久精品77| 在线观看美女被高潮喷水网站 | 国产成人系列免费观看| av专区在线播放| 久久精品91无色码中文字幕| 亚洲真实伦在线观看| 啦啦啦观看免费观看视频高清| 国产精品 国内视频| 国产三级在线视频| 脱女人内裤的视频| 久久久久国内视频| 一级a爱片免费观看的视频| 麻豆一二三区av精品| 午夜福利视频1000在线观看| 一a级毛片在线观看| 国产亚洲精品av在线| 真人一进一出gif抽搐免费| 久久婷婷人人爽人人干人人爱| 欧美乱色亚洲激情| 一进一出抽搐动态| 深夜精品福利|