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

    建筑體型廓線對建筑物表面粒子沉積特性的影響

    2017-03-13 08:16:18華鳳皎亢燕銘
    關(guān)鍵詞:背風(fēng)面體型通量

    華鳳皎, 亢燕銘, 鐘 珂

    (東華大學(xué) 環(huán)境科學(xué)與工程學(xué)院,上海 201620)

    建筑體型廓線對建筑物表面粒子沉積特性的影響

    華鳳皎, 亢燕銘, 鐘 珂

    (東華大學(xué) 環(huán)境科學(xué)與工程學(xué)院,上海 201620)

    選取典型建筑體型,利用計算流體力學(xué)方法研究了建筑體型對氣溶膠粒子在建筑表面沉積速度和沉積通量的影響.結(jié)果表明:迎風(fēng)面寬度是影響氣溶膠粒子在建筑表面沉積速度和沉積通量的重要因素,且迎風(fēng)面和背風(fēng)面對其的影響規(guī)律不同;建筑表面上氣溶膠粒子的沉積通量隨著其粒徑的變化而變化,且不同建筑表面最小沉積通量對應(yīng)的粒徑不同.

    建筑體型; 建筑外表面;慣性沉積; 沉積速度; 沉積通量; 氣溶膠粒子

    目前,城市大氣霾污染已成為我國重要的環(huán)境問題.工業(yè)過程、機動車尾氣以及其他人類活動產(chǎn)生的氣溶膠粒子,不僅對人體健康產(chǎn)生威脅,也將對建筑物、構(gòu)筑物以及其他各種市政設(shè)施產(chǎn)生嚴重污染腐蝕和損害,如建筑表面的粒子沉積會影響其圍護結(jié)構(gòu)的熱工性能及美觀要求.近幾年,隨著霾天氣在我國大范圍地區(qū)的頻繁出現(xiàn),這種氣溶膠污染效應(yīng)也變得越發(fā)嚴重.即使大氣中有害氣溶膠已被有效清除,沉積在建筑物表面的粒子仍然會長時間地影響建筑物附近的空氣質(zhì)量.為了減少氣溶膠粒子造成的危害,須對這些粒子沉積特征有清楚的認識.此外,古建筑的保護也需要了解和預(yù)測建筑各表面的粒子沉積情況,通過對建筑物墻體進行及時清洗,減弱氣溶膠對文物和古跡的腐蝕.

    粒子在固體表面的沉積特性研究報道較多.如文獻[1]給出了封閉空間內(nèi)均勻湍流粒子沉積的解析解.文獻[2]對粒子的沉積速度進行了參數(shù)化分析,給出了計算粒子干沉積速度的簡化公式.文獻[3]采用多層箱體模型模擬了包含多種地表形態(tài)的大尺度區(qū)域粒子的沉積通量,但并沒有求解區(qū)域內(nèi)的流動情況.文獻[4]采用離散軌跡模型研究了不同通風(fēng)房間內(nèi)粒子的質(zhì)量濃度分布和沉積特性.文獻[5]研究了城市環(huán)境下墊面、建筑物墻體和屋頂?shù)臍馊苣z粒子沉積速度,指出僅用一個或幾個物理量來表示城市氣溶膠粒子沉積速度的方法過于簡化了對問題的分析.

    上述文獻側(cè)重于研究粒子沉積速度與氣流分布的關(guān)系以及室內(nèi)粒子的沉積,均沒有考慮建筑體型對室外粒子在墻面沉積的影響.通常情況下,粒子的沉積效應(yīng)主要受擴散和慣性的共同作用[6].由于不同建筑物體型對氣流的阻擋效應(yīng)和其周圍流場特征不同,導(dǎo)致了粒子在不同體型建筑周圍隨氣流運動時的動力學(xué)行為的差異.

    為考察不同的建筑體型廓線下粒子在建筑物表面的沉積特征,本文將針對城市中常見的點式建筑和條形建筑體型,采用計算流體力學(xué)方法構(gòu)建粒子沉積模型,分析建筑體型對其圍護結(jié)構(gòu)外表面上的粒子沉積特性的影響.

    1 數(shù)值計算模型和方法

    1.1 數(shù)值計算模型

    表1給出了典型建筑體型的特征參數(shù),其中,建筑體型尺寸長為L,寬為W,高度為H.計算區(qū)域的大小對于計算量和模擬結(jié)果的準確性會造成影響,當(dāng)計算區(qū)域長度和高度不小于建筑物高度的5~6倍時,可消除邊界效應(yīng)對計算結(jié)果的影響[7-8].因此,本研究中計算區(qū)域高度取為6H,入口距建筑物迎風(fēng)面的距離為5H,出口距建筑物背風(fēng)面距離為10H.本文建筑體型case 2 L的計算區(qū)域如圖1所示.

    表1 建筑體型尺寸

    圖1 模型計算域及網(wǎng)格示意圖Fig.1 Schematic diagram of the computational domain and enlarged mesh

    采用Gambit軟件對求解區(qū)域進行非結(jié)構(gòu)化網(wǎng)格劃分,并對建筑外表面附近區(qū)域進行網(wǎng)格局部加密處理,而其他區(qū)域則采用正常尺寸網(wǎng)格來優(yōu)化整體網(wǎng)格數(shù)量并節(jié)約數(shù)值計算成本.模擬計算前,進行了網(wǎng)格獨立性檢驗,即增加網(wǎng)格數(shù),直至網(wǎng)格數(shù)的增加對計算結(jié)果沒有影響為止[9].經(jīng)檢驗,最終網(wǎng)格數(shù)量為163~291萬.

    1.2 計算方法和邊界條件的設(shè)定

    湍流流場可以采用不同的湍流模型進行模擬,為節(jié)約計算成本并獲得簡潔合理的數(shù)值結(jié)果,根據(jù)文獻[10-13]可知,標準κ-ε模型能很好地模擬本文條件下的湍流流場,其中壁面湍流特征用標準壁面函數(shù)描述,采用有限容積法離散控制方程,并對動量、湍流動能κ、能量耗散率ε的輸運方程采用二階迎風(fēng)格式進行差分,以提高計算精度.

    污染物輸送質(zhì)量濃度變化的對流擴散方程為

    (1)

    式中:Ci為污染物質(zhì)量濃度;K為粒子擴散系數(shù);Si為污染源項;uj為j方向的速度分量;xj為笛卡爾坐標分量.

    計算域的進口邊界條件設(shè)為速度入口,風(fēng)速大小采用速度廓線形式.風(fēng)速的垂直特性與地形及地面粗糙度有關(guān),速度廓線可由指數(shù)律分布[14-16]表示,如式(2)所示.

    (2)

    式中:z為任一高度, m;zs為參考高度, m;u(zs)為zs處的參考風(fēng)速, m/s, 考慮到城市區(qū)域的穩(wěn)定氣象條件,本文取20 m高度處的風(fēng)速為2 m/s;α為地面粗糙系數(shù)(風(fēng)剖面冪指數(shù)),本文α取0.25[17-18].速度廓線采用C++語言用戶編程并與Fluent程序提供的UDF(user-defined functions)對接.

    計算區(qū)域出口處可視為充分發(fā)展的湍流,采用出流邊界條件(outflow),即該流出區(qū)域的質(zhì)量流率和流入?yún)^(qū)域的質(zhì)量流率相等.建筑表面及地面采用無滑移條件,其余各面則采用滑移速度邊界.

    假設(shè)城市背景氣流中粒子質(zhì)量濃度均勻,因此,設(shè)每種粒徑粒子在計算域進口處空氣中的質(zhì)量濃度為10 μg/m3.城市大氣中,氣溶膠粒子的粒徑范圍很廣,生物有害性粒子尺寸(空氣動力學(xué)直徑)通常在1~10 μm范圍,其中2 μm左右及以下的粒子占絕大多數(shù)[3],粒子密度約為1 300 kg/m3[19].故本文選取粒子密度取ρp= 1 500 kg/m3,對粒徑dp=0.5, 1.0, 1.5, 2.5, 3.5, 5.0, 7.0和10.0 μm分別進行計算和分析.

    1.3 粒子輸運理論基礎(chǔ)和計算方法

    粒子在空氣中的輸運模擬采用歐拉法,選用混合模型求解第二相粒子的動量和體積分數(shù).大氣環(huán)境中粒子可視為稀薄相,只考慮流場對顆粒物的單向耦合作用.壁面處的粒子沉積通量J可用式(3)計算.

    (3)

    式中:CVOF-dA為靠近壁面控制體中粒子相體積分數(shù);vd-dA為某一點粒子的沉積速度; dA為控制體壁面的面積;Awall為建筑物表面面積.

    建筑物不同外表面上的粒子沉積速度vd為垂直面上的沉積速度vdv

    (4)

    朝上的水平面沉積速度vdu

    (5)

    式中:vs為粒子的終端沉降速度, m/s;u*為摩擦速度, m/s; 參量I的具體表達式見文獻[20];u*為通過流體的湍流動能k, m2/s,由式(6)計算得到.

    (6)

    1.4 數(shù)值模型的合理性驗證

    為了保證數(shù)值模擬方法分析建筑流場分布特征的可靠性,需要對所采用數(shù)學(xué)模型進行驗證.本文利用漢堡大學(xué)氣象研究所提供的BLASIUS風(fēng)洞試驗數(shù)據(jù)庫[21]對速度場進行驗證.模擬參數(shù)與試驗保持一致,模型采用邊長為125 mm的立方體,入口風(fēng)速采用速度廓線,參考高度500 mm處的參考風(fēng)速為6 m/s,風(fēng)速廓線指數(shù)為0.22.采用文獻[22]的試驗結(jié)果驗證濃度場,具體如下:風(fēng)洞中的模型尺寸為60 mm×60 mm,其長度與風(fēng)洞長度相等,充分發(fā)展區(qū)的風(fēng)速廓線指數(shù)為0.20,參考高度650 mm處的參考風(fēng)速為3 m/s.污染物質(zhì)量濃度采用無因次濃度C*來表示,如式(7)所示.

    (7)

    式中:Ci為i點處的污染物質(zhì)量濃度, kg/m3;q0為污染物散發(fā)強度, kg/s;Ls為污染源長度,取0.9 m.

    建筑迎風(fēng)區(qū)和背風(fēng)區(qū)的數(shù)值計算結(jié)果與風(fēng)洞試驗結(jié)果如圖2所示,其中u(z)為沿?zé)o量綱高度Z(Z=z/H)方向的氣流速度分量.由圖2(a)可以看出,本文模擬結(jié)果在建筑物高度范圍內(nèi)與風(fēng)洞試驗結(jié)果吻合得比較好.圖2(b)給出的污染物質(zhì)量濃度分布驗證說明,迎風(fēng)區(qū)質(zhì)量濃度計算值與實測值符合得較好,背風(fēng)區(qū)幾乎完全吻合.因此,本文采用的數(shù)值模型可以很好地預(yù)測城市中單體建筑物附近的氣流流動狀況和粒子質(zhì)量濃度分布情況.

    (a) 氣流速度

    (b) 質(zhì)量濃度

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

    2.1 建筑體型對建筑附近氣流速度的影響

    對表1給出的各建筑體型逐一進行了模擬.在case 1L和case 3L兩種典型情況下,在z=5 m高度處建筑周圍的氣流流線分布如圖3所示.

    圖3 典型建筑z=5 m平面上的氣流流線分布Fig.3 Streamline distributions on typical planes of the buildings at z=5 m

    由圖3可以看出,由于建筑物對來流的阻擋,導(dǎo)致在建筑側(cè)面及背風(fēng)面區(qū)域均會形成渦流.由于case 3L的迎風(fēng)面寬度較大,其背風(fēng)面形成的回漩渦流區(qū)域也明顯大于case 1L. case 3L側(cè)面近壁區(qū)完全處于渦流狀態(tài),而case 1L側(cè)面較長,渦流主要分布在擾流區(qū)前端.二次流的分布不同可能會導(dǎo)致粒子沉積到建筑物表面的概率不同,從而造成建筑表面粒子沉積通量的差異.圖3還表明,不同體型的建筑物附近的氣流流動速度不同,可能會導(dǎo)致建筑物不同表面的粒子沉積速度有所不同.

    為定量分析建筑各表面附近氣流速度的差異,圖4給出了距離建筑物背風(fēng)面、迎風(fēng)面和側(cè)面0.25 m處投影平面的垂直中心線上氣流速度沿?zé)o量綱高度Z的變化情況.

    (a) 迎風(fēng)面

    (b) 背風(fēng)面

    (c) 側(cè)面

    由圖4(a)可見,在迎風(fēng)面附近,高度相同、體型不同的建筑物附近氣流速度變化曲線相似,僅在Z<0.3的近地區(qū)域略有不同.高度不同的建筑物附近氣流速度變化趨勢也基本一致,這表明建筑體型和高度對迎風(fēng)面中心附近氣流速度的影響很?。畯膱D4(b)和4(c)中可看出,建筑物背風(fēng)面和側(cè)面的氣流速度受建筑體型和迎風(fēng)面寬度影響很大.不同體型的建筑物背風(fēng)面在Z<0.3的近地區(qū)域氣流速度基本相等,上部流速隨著迎風(fēng)面寬度增大而增大,且建筑物越高,規(guī)律越明顯.另外,較高建筑的背風(fēng)區(qū)氣流速度也高于相同平面結(jié)構(gòu)的較低建筑.

    由圖4(c)還可以看到,在case 2L/case 2H和case 3L/case 3H兩種情況下,建筑側(cè)面近地區(qū)域的氣流速度很接近,且遠大于case 1L/case 1H的氣流速度.這是因為前者建筑側(cè)面尺寸較小,建筑側(cè)面完全位于渦流區(qū)內(nèi),流速受渦流影響而明顯增大,而后者側(cè)面部分的前端區(qū)域位于渦流區(qū)內(nèi),而所考察的速度樣本點位置已在渦流區(qū)外,故氣流速度較?。@一差別將會對粒子在建筑物表面的沉積特性產(chǎn)生影響.

    2.2 建筑各表面上的粒子沉積通量

    由數(shù)值模擬結(jié)果和式(3)沉積通量的計算可以得到表1中6種情形的迎風(fēng)面、側(cè)面、背風(fēng)面和屋頂上粒子沉積通量J與粒徑dp的關(guān)系,結(jié)果如圖5所示.

    由圖5(a)可知,不論建筑高低,均表現(xiàn)為迎風(fēng)面越寬,迎風(fēng)面上的沉積通量越小,此外,粒子粒徑越小,沉積通量受建筑體型影響越小.這是因為小粒子的沉積主要受擴散機理控制,而大粒子的沉積主要受慣性力控制.由于建筑體型是通過影響周圍流場來影響粒子沉積效果的,因此只能影響到受慣性力控制的大粒子的沉積效果.

    (a) 迎風(fēng)面和側(cè)面

    (b) 背風(fēng)面

    (c) 屋頂

    (d) 建筑物總表面

    另外,比較不同建筑體型中迎風(fēng)面最小沉積通量對應(yīng)的粒子直徑可以看到,除了迎風(fēng)面面積較小的case 1H迎風(fēng)面上的最小沉積通量對應(yīng)粒徑為2.5 μm外,其他體型建筑對應(yīng)粒徑均為3.5 μm.這是因為case 1H對氣流的阻擋作用小,迎風(fēng)面附近氣流平均速度大于其他建筑體型的情況,粒子受到的慣性作用大,故導(dǎo)致最小沉積通量對應(yīng)的粒徑小于其他建筑體型.由圖5(a)還可以看到,雖然側(cè)面粒子沉積通量略小于迎風(fēng)面,但沉積通量隨粒徑的變化曲線與迎風(fēng)面非常相似,最小沉積通量對應(yīng)的粒徑也是3.5 μm.

    由圖5(b)可見,背風(fēng)面上沉積通量隨粒徑的變化規(guī)律為,迎風(fēng)面越寬,背風(fēng)面粒子沉積通量越大.此外,對比圖5(a)和5(b)還可以看到,背風(fēng)面上的沉積通量明顯小于迎風(fēng)面,同時最小沉積通量對應(yīng)粒徑較大,約為5.0~7.0 μm.這是因為背風(fēng)面附近氣流速度遠低于迎風(fēng)面,故粒子在這里受到的慣性力較小,于是最小沉積通量對應(yīng)的dp必然增大.由圖5(c)可看出,在重力沉降的作用下,屋頂上的沉積通量隨著粒徑的增加單調(diào)上升.由于屋頂表面附近的局地風(fēng)速受建筑體型的影響較小,故不同建筑體型的沉積通量差別很小.由圖5(d)可知,由于屋頂水平面上的沉積通量遠大于其他垂直表面的沉積通量,建筑物表面總的沉積通量主要取決于屋頂上的沉積通量,因此,其與屋頂上呈現(xiàn)出相同的沉積規(guī)律.

    2.3 建筑體型對沉積速度的影響

    根據(jù)氣溶膠動力學(xué)原理[23-24],粒子的沉積速度與建筑物附近氣流速度密切相關(guān).為了對建筑物表面沉積情況進行總體評價,在距離建筑物表面0.25 m的投影平面上均勻取點,每個點間距1 m,取值點總數(shù)為442~2 502,提取這些取樣點處的數(shù)值結(jié)果來討論氣溶膠粒子在外墻面的沉積規(guī)律. 根據(jù)流場的數(shù)值模擬結(jié)果和式(4)和(5),在不同建筑體型下,dp=0.5和10.0 μm的粒子在建筑迎風(fēng)面、背風(fēng)面、側(cè)面和屋頂上的沉積速度統(tǒng)計結(jié)果如圖6所示.

    (a) 迎風(fēng)面

    (b) 背風(fēng)面

    (c) 側(cè)面

    (d) 屋頂

    由圖6(a)可見,在兩種建筑高度下,均表現(xiàn)為迎風(fēng)面越窄,粒子沉積速度越大的特征.由圖6(b)可知,背風(fēng)面與迎風(fēng)面表現(xiàn)出的規(guī)律則恰好相反.這一規(guī)律在高建筑(即H=40 m)上表現(xiàn)得更為明顯.這可能是由于迎風(fēng)面尺度越大,背風(fēng)區(qū)的湍動能亦越大,故式(6)中的摩擦速度將增加,從而使得沉積速度增大.

    圖6(c)中表明,粒子在建筑側(cè)面上沉積速度與迎風(fēng)面沉積速度在量級上相當(dāng),但與建筑體型和迎風(fēng)面寬度的關(guān)系并沒有非常明確的相關(guān)性.這很可能是因為,不同體型的側(cè)面寬度不同,而建筑側(cè)面空氣流動情況受到迎風(fēng)面和側(cè)面尺寸的共同影響.由圖6(d)可見,dp=10.0 μm的粒子在屋頂上的水平沉積速度明顯高于dp=0.5 μm粒子的沉積速度.這是因為大粒子還會受到重力沉積的作用,從而強化了大粒子的沉積.此外,屋頂上的沉積速度與建筑體型和迎風(fēng)面寬度的相關(guān)性與背風(fēng)面相同,即沉積速度隨著迎風(fēng)面寬度的增大而增大.綜合以上分析,可以認為迎風(fēng)面寬度W與沉積速度vdv存在明確的相關(guān)性.

    為評價建筑物迎風(fēng)面寬度W對粒子沉積效果的影響,現(xiàn)以迎風(fēng)面寬度最小的case 1L和case 1H為基準,定義其他建筑模型表面不同dp粒子的沉積速度和沉積通量變化率η,如式(8)所示.

    (8)

    式中:Y代表沉積速度vdv或沉積通量J;x=2, 3,分別為表1所示的case 2和case 3類型的建筑;i=L, H,分別表示低建筑和高建筑.

    根據(jù)以上定義,圖7給出了不同建筑體型迎風(fēng)面上的粒子沉積通量變化率η與其粒徑的關(guān)系.

    圖7 不同建筑迎風(fēng)面上不同粒子沉積通量 變化率隨其粒徑的變化Fig.7 Changed ratios of deposition flux with dpon windward of different buildings

    由圖7可以看到,建筑迎風(fēng)面寬度的增大將使得所有粒徑的粒子在迎風(fēng)面的沉積通量均產(chǎn)生不同程度的減小,這與上述迎風(fēng)面寬度對粒子沉積速度的影響的結(jié)論相吻合.由圖7還可以看到,迎風(fēng)面寬度對粒子沉積通量的影響程度隨著粒徑和建筑高度的增大而增加.

    為了比較建筑體型對粒子沉積速度和沉積通量的影響程度,圖8分別給出了dp=0.5和10.0 μm的粒子在建筑迎風(fēng)面和背風(fēng)面上沉積速度與沉積通量變化率η.

    (a) 迎風(fēng)面

    (b) 背風(fēng)面

    由圖8(a)可以看到,迎風(fēng)面上的粒子沉積速度與沉積通量的變化率η的變化規(guī)律保持一致,隨著迎風(fēng)面寬度和建筑高度的增加,迎風(fēng)面上的粒子沉積速度和沉積通量都減小,且粒徑越大越明顯.由圖8(b)可以看到:對于較高的建筑,迎風(fēng)面寬度較大的體型使得背風(fēng)面上大粒子和小粒子沉積速度和沉積通量都增大,但沉積通量的增大幅度明顯低于沉積速度;對于較低建筑,背風(fēng)面沉積速度隨著迎風(fēng)面寬度增大而增加,但沉積通量卻出現(xiàn)了減小現(xiàn)象.這是因為沉積通量由沉積速度和粒子質(zhì)量濃度共同決定,沉積通量沒有完全隨著沉積速度的增大而增加,這可能是建筑表面粒子質(zhì)量濃度隨迎風(fēng)面寬度變化的規(guī)律與沉積速度不一致造成的.為此,圖9給出了不同建筑體型下建筑中心剖面上dp=0.5 μm粒子的質(zhì)量濃度分布.

    圖9 不同建筑體型中心剖面上的粒子質(zhì)量濃度分布Fig.9 Particle concentration distributions on centre section of different buildings

    由圖9可以看到,由于沉積效應(yīng)對粒子的清除作用,隨著迎風(fēng)面寬度的增大,背風(fēng)面附近的粒子質(zhì)量濃度下降,因此,使得粒子沉積通量隨迎風(fēng)面寬度增加的幅度減小(對于較高建筑).盡管隨著建筑迎風(fēng)面寬度增加,粒子沉積速度增大,但當(dāng)沉積速度的增大量不足以抵消粒子質(zhì)量濃度的減小量時,沉積通量隨迎風(fēng)面寬度增大的變化趨勢就不再與沉積速度相同.

    3 結(jié) 語

    建筑物表面的粒子沉積不僅嚴重影響建筑外觀,而且在雨雪沖刷和太陽輻射的共同作用下,還會因外墻材料的腐蝕而改變圍護結(jié)構(gòu)的熱工特性和室內(nèi)熱環(huán)境.由于建筑物的體型對其周圍流場有重要的畸變作用,進而影響到隨背景空氣輸送的大氣氣溶膠在建筑物表面上的沉積特性,但相關(guān)的研究還很缺乏.因此,本文應(yīng)用粒子沉積模型和計算流體力學(xué)方法,對幾種典型體型廓線的建筑表面上粒子沉積特性進行了數(shù)值計算和分析,討論了不同表面上發(fā)生粒子沉積的特征.主要結(jié)論如下:

    (1) 由于不同體型建筑對來流的阻擋作用不同,形成了不同的流場特征.本文模擬結(jié)果表明,建筑迎風(fēng)面尺寸對流場的影響作用大于建筑高度的影響效果.因此,建筑高度增大雖然增大了迎風(fēng)面面積,但對粒子沉積速度的影響沒有迎風(fēng)面寬度變化的作用明顯.

    (2) 隨著建筑物迎風(fēng)面寬度的增大,建筑迎風(fēng)面上粒子沉積速度將會下降,但在背風(fēng)面上的沉積速度則會升高.建筑物迎風(fēng)面上的粒子沉積通量隨迎風(fēng)面寬度的變化規(guī)律與沉積速度一致,但對于背風(fēng)面,由于墻面附近粒子質(zhì)量濃度受迎風(fēng)面寬度的影響趨勢與沉積速度不同,沉積通量與沉積速度變化規(guī)律不完全一致.

    (3) 由于不同建筑表面附近氣流慣性作用不同,故不同表面上最小沉積通量對應(yīng)的粒徑也不同.由于背風(fēng)面處在建筑尾流區(qū),氣流平均速度低,故粒子的慣性和湍流擴散效應(yīng)均較弱,因此,背風(fēng)面最小沉積通量對應(yīng)的粒徑將大于迎風(fēng)面和側(cè)面.

    [1] CORNER J, PENDLEBURY E D. The coagulation and deposition of a stirred aerosol[J]. Proceeding of the Physical Society, 1951,64(8):645-654.

    [2] ZHANG L, GONG S, PADRO J, et al. A size-segregated particle dry deposition scheme for an atmospheric aerosol module[J]. Atmospheric Environment, 2001,35(3):549-560.

    [3] MONFORTI F, BELLASIO R, CLAI G, et al. An evaluation of particle deposition fluxes to cultural heritage sites in Florence, Italy[J]. Science of the Total Environment, 2004,334/335(1):61-72.

    [4] ZHAO B, ZHANG Y, LI X T, et al. Comparison of indoor aerosol particle concentration and deposition in different ventilated rooms by numerical method[J]. Building and Environment, 2004,39(1):1-8.

    [5] JONSSON L, KARLSSON E, JONSSON P. Aspects of particulate dry deposition in the urban environment[J]. Journal of Hazardous Materials, 2008,153(1):229-243.

    [6] GREENFIELD S M. Rain scavenging of radioactive particulate matter from the atmosphere[J]. Journal of Atmospheric Sciences, 1957,14(2):115-125.

    [7] TOMINAGA Y, MOCHIDA A, YOSHIE R, et al. AIJ guidelines for practical applications of CFD to pedestrian wind environment around buildings[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2008,96(10):1749-1761.

    [8] NAKAMURA Y,OKE T R. Wind temperature and stability conditions in an east-west oriented urban canyon[J]. Atmospheric Environment, 1988,22(12):2691-2700.

    [9] LEI Y G, HE Y L. Effects of baffle inclination angle on flow and heat transfer of a heat exchanger with helical baffles[J]. Chemical Engineering and Processing: Process Intensification, 2008,47(12):2336-2345.

    [10] LAUNDER B E, SPALDING D B. The numerical computation of turbulent flows[J]. Computer Methods in Applied Mechanics and Engineering, 1974,3(2):269-289.

    [11] NAZRIDOUST K, AHMADI G. Airflow and pollutant transport in street canyons[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2006,94(6):491-522.

    [12] DI SABATINO S, BUCCOLIERI R, PULVIRENTI B, et al. Flow and pollutant dispersion in street canyons using FLUENT and ADMS-Urban[J]. Environ Mental Modeling & Assessment, 2007,13(3):369-381.

    [13] NG W Y, CHAU C K. A modeling investigation of the impact of street and building configurations on personal air pollutant exposure in isolated deep urban canyons[J]. Science of the Total Environment, 2014,468/469:429-448.

    [14] CHEUNG J O P, LIU C H. CFD simulations of natural ventilation behaviour in high-rise buildings in regular and staggered arrangements at various spacings[J]. Energy and Buildings, 2011,43(5):1149-1158.

    [15] LEUNG K K, LIU C H, WONG C C C, et al. On the study of ventilation and pollutant removel over idealized two-dimensional urban street canyons[J]. Building Simulation, 2012,5(4):359-369.

    [16] RAJAPAKSHA I, NAGAI H, OKUMIYA M. A ventilated courtyard as a passive cooling strategy in the warm humid tropics[J]. Renewable Energy, 2003,28(11):1755-1778.

    [17] MOHAMED F Y.A wind tunnel study on the effect of thermal stability on flow and dispersion of rooftop stack emissions in the near wake of a building[J]. Atmospheric Environment, 2013,65:89-100.

    [18] WANG H J, CHEN Q Y. A new empirical model for predicting single-sided, wind-driven natural ventilation in buildings[J]. Energy and Buildings, 2012,54:386-394.

    [19] TISA L S, KOSHIKAWA T, GERHARDT P.Wet and dry bacterial spore densities determined by buoyant sedimentation[J]. Applied and Environmental Microbiology, 1982,43(6):1307-1310.

    [20] LAI A C K, NAZAROFF W W. Modeling indoor particle deposition from turbulent flow onto smooth surfaces[J]. Journal of Aerosol Science, 2000,31(4):463-476.

    [21] SCHATZMANN M. Compilation of experimental data for validation of microscale dispersion models[DB]. (1999-05-15)[2015-03-20]. http://www.mi.zmaw.de/index.php?id=628.

    [22] MERONEY R N, PAVAGEAU M, RAFADALIS S, et al. Study of line source characteristics for 2- D physical modelling of pollutant dispersion in street canyon[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1996,62(1):37-56.

    [23] FRIEDLANDER S K, JOHNSTONE H F. Deposition of suspended particles from turbulent gas streams[J]. Journal of Industrial and Engineering Chemistry, 1957,49(7):1151-1156.

    [24] FUCHS N A. The mechanics of aerosols[M]. New York: Dover Publications Inc,1964:159-170.

    (責(zé)任編輯:劉園園)

    Effects of Building Configuration on the Particle Deposition Features on Building Surfaces

    HUAFengjiao,KANGYanming,ZHONGKe

    (School of Environmental Science and Engineering, Donghua University, Shanghai 201620, China)

    The influences of building configuration on the deposition velocity and deposition flux of aerosol particles on the building surface were investigated by using computational fluid dynamics for typical building configurations. The numerical results show that the width of windward side of a building is a key factor to influence deposition velocity and deposition flux. Deposition velocities and deposition fluxes of aerosol particles are different on windward and leeward sides of the building. Deposition fluxes of aerosol particles on the building surfaces vary with particle diameters. Particle diameters corresponding to minimum deposition flux are also different for different building surfaces.

    building configuration; building surface; inertia deposition; deposition velocity; deposition flux; aerosol particles

    1671-0444(2017)01-0115-08

    2016-03-23

    國家自然科學(xué)基金資助項目(40975093);上海市教委科研創(chuàng)新重點資助項目(14ZZ073)

    華鳳皎(1989—),女,山東煙臺人,博士研究生,研究方向為城市大氣環(huán)境與室內(nèi)空氣品質(zhì). E-mail: 2111173@mail.dhu.edu.cn 鐘 珂(聯(lián)系人),女,教授,E-mail:zhongkeyx@dhu.edu.cn

    X 513; P 426.6

    A

    猜你喜歡
    背風(fēng)面體型通量
    體型增大,同樣靈巧
    車主之友(2022年5期)2022-11-23 07:23:28
    冬小麥田N2O通量研究
    體型消瘦,如何科學(xué)變強壯
    體型特別的勇士
    幽默大師(2019年10期)2019-10-17 02:09:12
    體型最小的貓,竟然是個恐怖的大胃王
    文苑(2019年14期)2019-08-09 02:14:30
    非均勻等離子體Ka-Band傳輸性能中繼法優(yōu)化研究
    高超聲速風(fēng)洞子母彈大迎角拋殼投放試驗
    高壓輸電鐵塔塔身背風(fēng)面風(fēng)荷載遮擋效應(yīng)研究
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    化工進展(2015年6期)2015-11-13 00:26:29
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    国产一区二区三区在线臀色熟女| 18禁裸乳无遮挡免费网站照片 | 免费无遮挡裸体视频| netflix在线观看网站| 高潮久久久久久久久久久不卡| 精品第一国产精品| 人妻丰满熟妇av一区二区三区| 国产精品一区二区免费欧美| 久久 成人 亚洲| 免费一级毛片在线播放高清视频 | 欧美日本中文国产一区发布| 亚洲狠狠婷婷综合久久图片| 黄色视频,在线免费观看| 亚洲精品粉嫩美女一区| 欧美日韩亚洲综合一区二区三区_| 久久中文字幕一级| 中亚洲国语对白在线视频| 真人一进一出gif抽搐免费| 国产av一区在线观看免费| 韩国av一区二区三区四区| 国产精品综合久久久久久久免费 | 精品日产1卡2卡| 在线观看舔阴道视频| 亚洲专区中文字幕在线| 亚洲av电影在线进入| 精品欧美国产一区二区三| 正在播放国产对白刺激| 夜夜看夜夜爽夜夜摸| 欧美在线黄色| 亚洲精品中文字幕一二三四区| aaaaa片日本免费| 日日摸夜夜添夜夜添小说| 亚洲一区二区三区色噜噜| 宅男免费午夜| 男人的好看免费观看在线视频 | 久久久精品国产亚洲av高清涩受| 欧美最黄视频在线播放免费| 久久久久久人人人人人| 日韩三级视频一区二区三区| 美女 人体艺术 gogo| 国产在线精品亚洲第一网站| 丝袜在线中文字幕| 国产伦人伦偷精品视频| 成年女人毛片免费观看观看9| 正在播放国产对白刺激| 日韩欧美在线二视频| 亚洲aⅴ乱码一区二区在线播放 | 久久婷婷成人综合色麻豆| 999久久久国产精品视频| 日韩 欧美 亚洲 中文字幕| 一二三四社区在线视频社区8| 两个人视频免费观看高清| 国产麻豆69| 九色国产91popny在线| 欧美性长视频在线观看| 黄色 视频免费看| 精品人妻在线不人妻| 亚洲欧美一区二区三区黑人| 欧美日韩福利视频一区二区| 久久久久久久久免费视频了| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品aⅴ一区二区三区四区| 天堂动漫精品| 久久久久久久精品吃奶| 国产麻豆成人av免费视频| 亚洲天堂国产精品一区在线| 天天一区二区日本电影三级 | 欧美成人一区二区免费高清观看 | 老司机福利观看| 欧美乱码精品一区二区三区| 一级a爱片免费观看的视频| 欧美激情高清一区二区三区| 日本在线视频免费播放| av免费在线观看网站| 在线永久观看黄色视频| 国产欧美日韩一区二区精品| 9191精品国产免费久久| av欧美777| 啪啪无遮挡十八禁网站| 日韩欧美国产在线观看| 精品电影一区二区在线| 嫩草影院精品99| 日本vs欧美在线观看视频| 韩国av一区二区三区四区| 国产野战对白在线观看| 久久精品国产亚洲av高清一级| 女人被狂操c到高潮| 久久精品亚洲熟妇少妇任你| 免费av毛片视频| a在线观看视频网站| 久久国产精品人妻蜜桃| 午夜福利高清视频| 搡老妇女老女人老熟妇| 欧美日韩亚洲国产一区二区在线观看| 亚洲熟女毛片儿| 88av欧美| av片东京热男人的天堂| 欧美日韩精品网址| 啦啦啦 在线观看视频| 女警被强在线播放| 色婷婷久久久亚洲欧美| 国产真人三级小视频在线观看| 最近最新中文字幕大全免费视频| 黄片播放在线免费| 久久人人97超碰香蕉20202| 精品久久久久久久毛片微露脸| 亚洲一区高清亚洲精品| 在线观看免费视频日本深夜| 欧美激情极品国产一区二区三区| 国产精品秋霞免费鲁丝片| 9191精品国产免费久久| 精品国产美女av久久久久小说| 亚洲五月婷婷丁香| 欧美乱码精品一区二区三区| 久久天堂一区二区三区四区| 非洲黑人性xxxx精品又粗又长| 精品国产国语对白av| 久久人人精品亚洲av| 一夜夜www| 精品久久久久久成人av| 一进一出抽搐gif免费好疼| 成人免费观看视频高清| 精品久久蜜臀av无| 不卡一级毛片| 免费观看人在逋| 成熟少妇高潮喷水视频| 每晚都被弄得嗷嗷叫到高潮| 91成年电影在线观看| 亚洲成国产人片在线观看| 精品福利观看| 美女高潮喷水抽搐中文字幕| 亚洲一区高清亚洲精品| 动漫黄色视频在线观看| 极品教师在线免费播放| 法律面前人人平等表现在哪些方面| 黑人巨大精品欧美一区二区mp4| 99热只有精品国产| 桃色一区二区三区在线观看| 欧美性长视频在线观看| 人人妻,人人澡人人爽秒播| 亚洲精品美女久久久久99蜜臀| 在线观看日韩欧美| 色老头精品视频在线观看| 久久天堂一区二区三区四区| 97超级碰碰碰精品色视频在线观看| 国产精品乱码一区二三区的特点 | 国产精品电影一区二区三区| 日韩 欧美 亚洲 中文字幕| 亚洲一码二码三码区别大吗| 成人亚洲精品av一区二区| 一级a爱视频在线免费观看| 欧美黑人欧美精品刺激| 可以免费在线观看a视频的电影网站| 精品福利观看| 精品久久久久久久毛片微露脸| 精品日产1卡2卡| 黄网站色视频无遮挡免费观看| 国产三级黄色录像| 国产成人啪精品午夜网站| 亚洲成人免费电影在线观看| 午夜亚洲福利在线播放| 夜夜躁狠狠躁天天躁| 嫩草影视91久久| 亚洲欧美一区二区三区黑人| 高潮久久久久久久久久久不卡| 婷婷精品国产亚洲av在线| 日日干狠狠操夜夜爽| 丰满人妻熟妇乱又伦精品不卡| 精品一区二区三区视频在线观看免费| 日本a在线网址| 久久人妻av系列| 精品电影一区二区在线| 超碰成人久久| 日本欧美视频一区| 精品免费久久久久久久清纯| 国产欧美日韩综合在线一区二区| 看免费av毛片| 9色porny在线观看| 亚洲精品国产一区二区精华液| 波多野结衣巨乳人妻| 国产亚洲av高清不卡| 国产成人系列免费观看| 久久精品国产亚洲av香蕉五月| 宅男免费午夜| 久热这里只有精品99| 久久欧美精品欧美久久欧美| 女人被躁到高潮嗷嗷叫费观| 免费在线观看影片大全网站| 妹子高潮喷水视频| 精品不卡国产一区二区三区| 又紧又爽又黄一区二区| 欧美精品亚洲一区二区| 亚洲国产高清在线一区二区三 | 久久久久久人人人人人| 色综合亚洲欧美另类图片| 亚洲精华国产精华精| 人人澡人人妻人| 国产亚洲精品第一综合不卡| 波多野结衣一区麻豆| 亚洲自偷自拍图片 自拍| 国产精品免费一区二区三区在线| 亚洲欧美精品综合久久99| 精品国产超薄肉色丝袜足j| 黑丝袜美女国产一区| 高清在线国产一区| 日韩成人在线观看一区二区三区| 久99久视频精品免费| 岛国在线观看网站| 免费高清在线观看日韩| 脱女人内裤的视频| 日韩欧美一区视频在线观看| 成人国语在线视频| 日本免费一区二区三区高清不卡 | 国产一区二区三区视频了| 免费看美女性在线毛片视频| 成人欧美大片| 男女做爰动态图高潮gif福利片 | 黑人巨大精品欧美一区二区蜜桃| 午夜免费激情av| av电影中文网址| 欧美黄色片欧美黄色片| 精品久久久久久成人av| 亚洲精品在线观看二区| 18禁裸乳无遮挡免费网站照片 | 欧美成人午夜精品| 最近最新中文字幕大全免费视频| 黑人操中国人逼视频| 在线永久观看黄色视频| 成熟少妇高潮喷水视频| 午夜福利18| 国产色视频综合| 无限看片的www在线观看| 国产精品一区二区三区四区久久 | 美女 人体艺术 gogo| 一级,二级,三级黄色视频| 女警被强在线播放| 九色国产91popny在线| 国产99白浆流出| 老汉色av国产亚洲站长工具| 大香蕉久久成人网| 夜夜躁狠狠躁天天躁| 日韩成人在线观看一区二区三区| 男女下面插进去视频免费观看| 欧美丝袜亚洲另类 | 99国产精品一区二区蜜桃av| 搡老妇女老女人老熟妇| 欧美激情高清一区二区三区| 一夜夜www| 亚洲五月天丁香| 久久久久国内视频| 色精品久久人妻99蜜桃| 高清毛片免费观看视频网站| 一本久久中文字幕| 久久精品成人免费网站| 日韩av在线大香蕉| 久久精品91无色码中文字幕| 午夜精品在线福利| 熟女少妇亚洲综合色aaa.| 日韩精品中文字幕看吧| 在线观看免费视频网站a站| 母亲3免费完整高清在线观看| 国产亚洲欧美精品永久| 99国产极品粉嫩在线观看| 亚洲精品美女久久久久99蜜臀| 真人一进一出gif抽搐免费| 老司机深夜福利视频在线观看| 日韩三级视频一区二区三区| 9色porny在线观看| 精品国产美女av久久久久小说| 在线av久久热| 神马国产精品三级电影在线观看 | 亚洲一码二码三码区别大吗| 欧美乱色亚洲激情| 亚洲人成77777在线视频| 亚洲情色 制服丝袜| 久久人人精品亚洲av| 国产亚洲av嫩草精品影院| 老司机深夜福利视频在线观看| 欧美激情极品国产一区二区三区| 国产人伦9x9x在线观看| 一区二区三区高清视频在线| 午夜影院日韩av| 成人国产综合亚洲| 婷婷精品国产亚洲av在线| 国产欧美日韩精品亚洲av| 人成视频在线观看免费观看| 国产激情久久老熟女| 在线观看66精品国产| 黄色女人牲交| 我的亚洲天堂| 12—13女人毛片做爰片一| 成人亚洲精品一区在线观看| 欧美黄色淫秽网站| 亚洲最大成人中文| 国产av一区二区精品久久| 日韩大码丰满熟妇| 看黄色毛片网站| 国产野战对白在线观看| 电影成人av| 国内精品久久久久精免费| 亚洲一区中文字幕在线| 此物有八面人人有两片| www.999成人在线观看| 桃红色精品国产亚洲av| 国产人伦9x9x在线观看| 麻豆久久精品国产亚洲av| 天堂√8在线中文| 日韩有码中文字幕| 久久国产亚洲av麻豆专区| 久久精品91蜜桃| 19禁男女啪啪无遮挡网站| 热99re8久久精品国产| 欧美乱码精品一区二区三区| 国产又爽黄色视频| 亚洲人成电影观看| 免费在线观看日本一区| 伊人久久大香线蕉亚洲五| 国产成人精品在线电影| 人妻久久中文字幕网| 成人国产一区最新在线观看| 黄色片一级片一级黄色片| 国产99白浆流出| 视频在线观看一区二区三区| 欧美精品亚洲一区二区| 亚洲五月天丁香| 日韩精品中文字幕看吧| 久久久久久久精品吃奶| 成人国产综合亚洲| av网站免费在线观看视频| 99国产精品一区二区三区| 极品教师在线免费播放| 亚洲一区二区三区不卡视频| 久久中文字幕一级| 无限看片的www在线观看| 91成年电影在线观看| 窝窝影院91人妻| 国产成年人精品一区二区| 天堂影院成人在线观看| 国产精品免费一区二区三区在线| 久久国产精品人妻蜜桃| 亚洲成av片中文字幕在线观看| 久久久久精品国产欧美久久久| 久久狼人影院| 欧洲精品卡2卡3卡4卡5卡区| 精品人妻在线不人妻| 午夜福利成人在线免费观看| 精品国产一区二区三区四区第35| 亚洲精品一区av在线观看| av有码第一页| √禁漫天堂资源中文www| 成人特级黄色片久久久久久久| 国产精品爽爽va在线观看网站 | 国产一区二区三区在线臀色熟女| 国产高清videossex| 男人舔女人的私密视频| 久久天躁狠狠躁夜夜2o2o| 窝窝影院91人妻| 啦啦啦免费观看视频1| 一区二区三区精品91| 香蕉久久夜色| 啦啦啦 在线观看视频| 老司机福利观看| 每晚都被弄得嗷嗷叫到高潮| 一级,二级,三级黄色视频| 波多野结衣巨乳人妻| 久久精品aⅴ一区二区三区四区| 欧美日韩一级在线毛片| 一级黄色大片毛片| 免费看美女性在线毛片视频| 啦啦啦韩国在线观看视频| 国产真人三级小视频在线观看| 麻豆国产av国片精品| 一边摸一边做爽爽视频免费| 国产精华一区二区三区| 亚洲国产欧美日韩在线播放| 午夜老司机福利片| 久久久国产欧美日韩av| 亚洲欧美精品综合久久99| 99re在线观看精品视频| 亚洲第一电影网av| 亚洲,欧美精品.| 最新美女视频免费是黄的| 亚洲男人天堂网一区| 亚洲情色 制服丝袜| 中文字幕人妻丝袜一区二区| 久久 成人 亚洲| videosex国产| 欧美乱色亚洲激情| 免费搜索国产男女视频| 午夜两性在线视频| 久久香蕉激情| 亚洲精品在线观看二区| 中出人妻视频一区二区| 国产主播在线观看一区二区| 91大片在线观看| 国产精品二区激情视频| 久久 成人 亚洲| 国产成人精品在线电影| 免费在线观看视频国产中文字幕亚洲| 99国产精品一区二区蜜桃av| 女同久久另类99精品国产91| 亚洲成av人片免费观看| 午夜福利高清视频| 亚洲一区高清亚洲精品| 一区福利在线观看| 久久人妻av系列| 亚洲熟妇中文字幕五十中出| 国产成人精品久久二区二区免费| 香蕉丝袜av| 欧美日韩瑟瑟在线播放| 18美女黄网站色大片免费观看| 侵犯人妻中文字幕一二三四区| 男人的好看免费观看在线视频 | 精品少妇一区二区三区视频日本电影| 如日韩欧美国产精品一区二区三区| 黄色成人免费大全| 亚洲精品中文字幕一二三四区| 看黄色毛片网站| 免费av毛片视频| 午夜精品在线福利| 国产伦一二天堂av在线观看| 久久天堂一区二区三区四区| 首页视频小说图片口味搜索| 亚洲中文字幕一区二区三区有码在线看 | 国产精品,欧美在线| 女人精品久久久久毛片| 999久久久精品免费观看国产| 精品免费久久久久久久清纯| 久久国产亚洲av麻豆专区| 国产主播在线观看一区二区| 99久久精品国产亚洲精品| 12—13女人毛片做爰片一| 啪啪无遮挡十八禁网站| 在线观看舔阴道视频| 这个男人来自地球电影免费观看| 国产熟女xx| 老汉色av国产亚洲站长工具| 亚洲av电影不卡..在线观看| 大型黄色视频在线免费观看| 日韩欧美免费精品| 午夜免费鲁丝| 国产av一区二区精品久久| 国内精品久久久久久久电影| 精品国产亚洲在线| 天天添夜夜摸| 91精品国产国语对白视频| 国产精品久久久久久精品电影 | 成人亚洲精品av一区二区| 人人妻,人人澡人人爽秒播| 啪啪无遮挡十八禁网站| 日韩国内少妇激情av| 久久久久久久午夜电影| 亚洲在线自拍视频| 国产高清激情床上av| 成人国语在线视频| 欧美一区二区精品小视频在线| 一本久久中文字幕| 久久国产精品影院| 一级,二级,三级黄色视频| 午夜福利影视在线免费观看| 久久久国产精品麻豆| 成人18禁在线播放| ponron亚洲| 日韩av在线大香蕉| 国产亚洲精品综合一区在线观看 | 久久精品亚洲精品国产色婷小说| 女性生殖器流出的白浆| 欧美久久黑人一区二区| 国产不卡一卡二| av片东京热男人的天堂| 亚洲一区二区三区色噜噜| 一区二区三区精品91| 黄色成人免费大全| 日韩av在线大香蕉| 色婷婷久久久亚洲欧美| 免费看美女性在线毛片视频| 99久久99久久久精品蜜桃| 成人免费观看视频高清| 久久久久久久久中文| 精品免费久久久久久久清纯| 日韩中文字幕欧美一区二区| 免费在线观看亚洲国产| 精品一品国产午夜福利视频| 精品国产乱子伦一区二区三区| 大陆偷拍与自拍| 欧美精品亚洲一区二区| 可以在线观看的亚洲视频| 欧美av亚洲av综合av国产av| 日本三级黄在线观看| 变态另类丝袜制服| 久久久久国产精品人妻aⅴ院| 国产精华一区二区三区| 欧美日本亚洲视频在线播放| 天堂影院成人在线观看| 亚洲人成77777在线视频| 亚洲色图av天堂| 欧美色视频一区免费| 91成人精品电影| 国产高清videossex| 欧美日韩中文字幕国产精品一区二区三区 | 老司机靠b影院| 亚洲免费av在线视频| 精品人妻在线不人妻| 丁香六月欧美| 免费人成视频x8x8入口观看| 久久婷婷成人综合色麻豆| 久久久国产成人精品二区| 国产精品久久久久久精品电影 | 亚洲伊人色综图| 国产亚洲欧美精品永久| 久久久水蜜桃国产精品网| 欧美性长视频在线观看| 国产av又大| 天天躁狠狠躁夜夜躁狠狠躁| 国产aⅴ精品一区二区三区波| or卡值多少钱| 精品少妇一区二区三区视频日本电影| 精品无人区乱码1区二区| 精品人妻在线不人妻| 91精品国产国语对白视频| 国产精品久久久久久人妻精品电影| 黄色 视频免费看| 久久精品国产99精品国产亚洲性色 | 免费看十八禁软件| 久久久久亚洲av毛片大全| 国内久久婷婷六月综合欲色啪| 久久久久久国产a免费观看| 9热在线视频观看99| 可以在线观看的亚洲视频| 日本 av在线| 丁香六月欧美| a级毛片在线看网站| 天天躁夜夜躁狠狠躁躁| 久久伊人香网站| 亚洲av日韩精品久久久久久密| 国语自产精品视频在线第100页| 女警被强在线播放| av视频免费观看在线观看| 国产精品野战在线观看| 久久久久久久久久久久大奶| 嫩草影视91久久| 国产片内射在线| 亚洲av电影在线进入| 国产亚洲av高清不卡| 两人在一起打扑克的视频| 18禁国产床啪视频网站| 欧美乱色亚洲激情| 男人操女人黄网站| 亚洲国产看品久久| 99国产精品99久久久久| 一边摸一边抽搐一进一小说| 国产男靠女视频免费网站| 制服人妻中文乱码| 日本免费a在线| 妹子高潮喷水视频| 精品日产1卡2卡| 久久香蕉国产精品| 在线观看66精品国产| 两人在一起打扑克的视频| 88av欧美| 成人国产综合亚洲| 久久久精品国产亚洲av高清涩受| 久久精品国产清高在天天线| 无限看片的www在线观看| 午夜精品在线福利| 一a级毛片在线观看| 日日爽夜夜爽网站| 亚洲五月天丁香| 免费看a级黄色片| 91老司机精品| 亚洲aⅴ乱码一区二区在线播放 | 老汉色∧v一级毛片| 精品午夜福利视频在线观看一区| 国产精品九九99| 亚洲欧美精品综合一区二区三区| 国产麻豆69| 岛国视频午夜一区免费看| 欧美日韩瑟瑟在线播放| 人人妻人人澡人人看| 免费av毛片视频| 又黄又粗又硬又大视频| 中文字幕人妻熟女乱码| 男女做爰动态图高潮gif福利片 | 女生性感内裤真人,穿戴方法视频| 国产精品1区2区在线观看.| 老司机午夜十八禁免费视频| 国产精品免费一区二区三区在线| 亚洲va日本ⅴa欧美va伊人久久| 欧美人与性动交α欧美精品济南到| 亚洲片人在线观看| 欧美乱码精品一区二区三区| 真人一进一出gif抽搐免费| 亚洲精品久久成人aⅴ小说| 国产亚洲精品av在线| 午夜福利视频1000在线观看 | 91老司机精品| 啪啪无遮挡十八禁网站| 最好的美女福利视频网| 欧美最黄视频在线播放免费| 1024视频免费在线观看| 美女大奶头视频| 免费一级毛片在线播放高清视频 | 婷婷精品国产亚洲av在线| 淫妇啪啪啪对白视频| 夜夜爽天天搞| 久久久国产欧美日韩av| 日日夜夜操网爽| 一边摸一边抽搐一进一小说| 97人妻精品一区二区三区麻豆 | 曰老女人黄片| 精品无人区乱码1区二区| 91精品三级在线观看| 一级a爱视频在线免费观看| www国产在线视频色| 久久精品人人爽人人爽视色|