• <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
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    成人特级黄色片久久久久久久 | 性高湖久久久久久久久免费观看| 一级,二级,三级黄色视频| 欧美黄色淫秽网站| 精品国产国语对白av| 国产无遮挡羞羞视频在线观看| 精品亚洲成a人片在线观看| 亚洲成人免费电影在线观看| 午夜精品久久久久久毛片777| 69av精品久久久久久 | 中文字幕人妻丝袜一区二区| a级片在线免费高清观看视频| 免费少妇av软件| 欧美+亚洲+日韩+国产| 老熟妇仑乱视频hdxx| 丁香六月天网| 巨乳人妻的诱惑在线观看| 亚洲精品久久成人aⅴ小说| 亚洲第一av免费看| 国产一区有黄有色的免费视频| 精品国产一区二区三区四区第35| 亚洲精品一二三| 成年人免费黄色播放视频| 90打野战视频偷拍视频| 成人免费观看视频高清| 久久久精品国产亚洲av高清涩受| 人人妻人人添人人爽欧美一区卜| 久久久久久久国产电影| 久久国产精品大桥未久av| 国产极品粉嫩免费观看在线| 王馨瑶露胸无遮挡在线观看| 狠狠精品人妻久久久久久综合| 天堂动漫精品| 免费在线观看黄色视频的| 午夜免费鲁丝| 亚洲国产成人一精品久久久| 国产一区二区激情短视频| 69精品国产乱码久久久| 久久国产精品人妻蜜桃| 欧美日韩黄片免| 黄片小视频在线播放| 一边摸一边抽搐一进一小说 | 欧美日韩av久久| 国产又色又爽无遮挡免费看| 亚洲国产欧美一区二区综合| 最新美女视频免费是黄的| 久久精品国产综合久久久| 久久精品国产a三级三级三级| 十八禁网站网址无遮挡| 亚洲九九香蕉| 日韩欧美免费精品| 欧美+亚洲+日韩+国产| 99久久人妻综合| av天堂在线播放| 99re在线观看精品视频| videos熟女内射| 国产精品香港三级国产av潘金莲| 精品一区二区三卡| 80岁老熟妇乱子伦牲交| 宅男免费午夜| 久久精品人人爽人人爽视色| 久久精品人人爽人人爽视色| 精品福利永久在线观看| 日本av免费视频播放| 久久久久久亚洲精品国产蜜桃av| 国产片内射在线| 丝袜美足系列| 手机成人av网站| avwww免费| 性色av乱码一区二区三区2| 亚洲精品乱久久久久久| 久久国产精品影院| 亚洲精品美女久久av网站| 免费观看a级毛片全部| av国产精品久久久久影院| 999精品在线视频| 十八禁网站网址无遮挡| 两人在一起打扑克的视频| 午夜福利影视在线免费观看| avwww免费| 久久精品91无色码中文字幕| 日日爽夜夜爽网站| 亚洲欧美色中文字幕在线| 一夜夜www| 国产亚洲欧美精品永久| 亚洲成人免费电影在线观看| 十八禁网站免费在线| 汤姆久久久久久久影院中文字幕| 国精品久久久久久国模美| 男人舔女人的私密视频| 日韩欧美免费精品| 欧美精品av麻豆av| 999精品在线视频| 亚洲精品粉嫩美女一区| 热99re8久久精品国产| 午夜免费成人在线视频| 大型av网站在线播放| 亚洲欧洲日产国产| 国产成人欧美| 精品国产一区二区三区久久久樱花| 女人精品久久久久毛片| 亚洲色图av天堂| 国产不卡av网站在线观看| 欧美乱妇无乱码| 亚洲精品av麻豆狂野| 乱人伦中国视频| 91av网站免费观看| av天堂久久9| 欧美亚洲日本最大视频资源| 性高湖久久久久久久久免费观看| 欧美另类亚洲清纯唯美| 999久久久精品免费观看国产| 精品人妻熟女毛片av久久网站| 99精品久久久久人妻精品| 亚洲成人手机| 一区福利在线观看| 欧美人与性动交α欧美精品济南到| 国产1区2区3区精品| 亚洲精品在线观看二区| 肉色欧美久久久久久久蜜桃| 亚洲情色 制服丝袜| 老司机靠b影院| 午夜福利免费观看在线| 国产精品欧美亚洲77777| 建设人人有责人人尽责人人享有的| 亚洲成av片中文字幕在线观看| 精品国产乱码久久久久久小说| 亚洲午夜理论影院| 国产精品久久久人人做人人爽| 黄片大片在线免费观看| 男男h啪啪无遮挡| 丁香六月天网| 久久九九热精品免费| 国产精品电影一区二区三区 | 亚洲成人国产一区在线观看| aaaaa片日本免费| 在线 av 中文字幕| 国产欧美亚洲国产| 国产激情久久老熟女| 一区二区三区激情视频| 露出奶头的视频| 久久久久视频综合| 在线观看www视频免费| 亚洲精品国产精品久久久不卡| 极品人妻少妇av视频| 久久精品亚洲熟妇少妇任你| 亚洲久久久国产精品| 精品久久久久久电影网| 国产精品99久久99久久久不卡| 一本一本久久a久久精品综合妖精| 日本一区二区免费在线视频| 久热爱精品视频在线9| 亚洲精品美女久久久久99蜜臀| 99精品久久久久人妻精品| 夜夜骑夜夜射夜夜干| 80岁老熟妇乱子伦牲交| 99久久精品国产亚洲精品| 一边摸一边抽搐一进一出视频| 亚洲第一青青草原| 首页视频小说图片口味搜索| 国产一区二区三区视频了| 日韩视频一区二区在线观看| 成人18禁在线播放| 极品少妇高潮喷水抽搐| 欧美精品一区二区免费开放| 老司机亚洲免费影院| 亚洲视频免费观看视频| 在线av久久热| 欧美亚洲日本最大视频资源| 性色av乱码一区二区三区2| av天堂久久9| 日韩制服丝袜自拍偷拍| a在线观看视频网站| 亚洲中文日韩欧美视频| 久久毛片免费看一区二区三区| 69精品国产乱码久久久| 欧美激情久久久久久爽电影 | 久久精品国产综合久久久| avwww免费| 69精品国产乱码久久久| 成年版毛片免费区| 国产av精品麻豆| 国产片内射在线| 久久人人97超碰香蕉20202| 久久人妻熟女aⅴ| 国产单亲对白刺激| 午夜福利影视在线免费观看| 飞空精品影院首页| 人妻 亚洲 视频| 少妇裸体淫交视频免费看高清 | 男人操女人黄网站| 久9热在线精品视频| 精品欧美一区二区三区在线| 9色porny在线观看| 老汉色∧v一级毛片| 欧美黑人欧美精品刺激| 国产熟女午夜一区二区三区| 日本精品一区二区三区蜜桃| 成年动漫av网址| 欧美精品人与动牲交sv欧美| 欧美人与性动交α欧美精品济南到| 国产一区二区 视频在线| 美女福利国产在线| 午夜精品国产一区二区电影| 久久av网站| 亚洲精华国产精华精| 欧美人与性动交α欧美精品济南到| 人人澡人人妻人| 天堂俺去俺来也www色官网| 欧美在线黄色| 精品免费久久久久久久清纯 | 国产男女超爽视频在线观看| 亚洲五月婷婷丁香| 久久影院123| 777久久人妻少妇嫩草av网站| 久久国产精品男人的天堂亚洲| 91av网站免费观看| 精品欧美一区二区三区在线| 美国免费a级毛片| 国产在线一区二区三区精| 精品视频人人做人人爽| 国产精品av久久久久免费| 国产精品美女特级片免费视频播放器 | 国产免费现黄频在线看| 日韩欧美三级三区| 久久久久国内视频| 久9热在线精品视频| 免费av中文字幕在线| 国产精品久久久久久精品电影小说| 精品国产一区二区三区四区第35| 少妇裸体淫交视频免费看高清 | 国产精品久久久人人做人人爽| 国产精品秋霞免费鲁丝片| 女同久久另类99精品国产91| 最近最新中文字幕大全免费视频| 亚洲欧美一区二区三区久久| 国产在线免费精品| 伦理电影免费视频| 午夜福利乱码中文字幕| 中文字幕色久视频| 新久久久久国产一级毛片| 亚洲五月婷婷丁香| 亚洲专区国产一区二区| 黄色视频,在线免费观看| 国产高清国产精品国产三级| 丝袜喷水一区| 99热国产这里只有精品6| 黄色视频在线播放观看不卡| 91av网站免费观看| netflix在线观看网站| 亚洲第一青青草原| 高清在线国产一区| 日韩精品免费视频一区二区三区| 久久这里只有精品19| 精品少妇一区二区三区视频日本电影| 欧美在线一区亚洲| 欧美人与性动交α欧美软件| 午夜福利在线免费观看网站| 窝窝影院91人妻| av视频免费观看在线观看| 亚洲成a人片在线一区二区| 国产无遮挡羞羞视频在线观看| 国产成人精品久久二区二区91| 精品少妇内射三级| 中文字幕人妻丝袜一区二区| 最近最新中文字幕大全电影3 | 亚洲成av片中文字幕在线观看| 国产亚洲欧美精品永久| 精品亚洲乱码少妇综合久久| 欧美日韩福利视频一区二区| 精品国产超薄肉色丝袜足j| 99国产精品一区二区蜜桃av | 91成人精品电影| 成年人黄色毛片网站| 亚洲午夜精品一区,二区,三区| 亚洲av电影在线进入| 国产在线精品亚洲第一网站| 国产精品国产av在线观看| 女人精品久久久久毛片| 精品人妻1区二区| 国产欧美日韩一区二区精品| 麻豆成人av在线观看| 狠狠狠狠99中文字幕| 国产成人免费无遮挡视频| 国产精品国产av在线观看| av又黄又爽大尺度在线免费看| www.自偷自拍.com| 99热国产这里只有精品6| 成人特级黄色片久久久久久久 | 夜夜骑夜夜射夜夜干| 美女高潮到喷水免费观看| 最新美女视频免费是黄的| 亚洲国产看品久久| 国产精品亚洲av一区麻豆| 中文字幕制服av| a在线观看视频网站| 美女国产高潮福利片在线看| 大型av网站在线播放| 91九色精品人成在线观看| 亚洲精品国产精品久久久不卡| 成人18禁高潮啪啪吃奶动态图| 久久国产精品男人的天堂亚洲| 日本vs欧美在线观看视频| 丝袜在线中文字幕| 最新的欧美精品一区二区| 欧美国产精品va在线观看不卡| 老司机影院毛片| 日日爽夜夜爽网站| 女人久久www免费人成看片| 国产日韩一区二区三区精品不卡| kizo精华| 精品第一国产精品| 久久精品国产亚洲av香蕉五月 | av不卡在线播放| 五月开心婷婷网| xxxhd国产人妻xxx| 欧美激情极品国产一区二区三区| 黄色怎么调成土黄色| 欧美老熟妇乱子伦牲交| 久久人妻av系列| 国产精品久久久人人做人人爽| 成人三级做爰电影| 日本欧美视频一区| 一级片'在线观看视频| 免费黄频网站在线观看国产| 久久精品国产99精品国产亚洲性色 | 日本精品一区二区三区蜜桃| 天堂8中文在线网| 欧美午夜高清在线| 在线播放国产精品三级| 蜜桃在线观看..| 多毛熟女@视频| 精品一区二区三卡| 法律面前人人平等表现在哪些方面| 亚洲第一av免费看| 美女扒开内裤让男人捅视频| 一区福利在线观看| 欧美国产精品一级二级三级| 久久毛片免费看一区二区三区| 91麻豆av在线| av在线播放免费不卡| 亚洲精品乱久久久久久| 美女午夜性视频免费| 日韩中文字幕视频在线看片| 精品一区二区三卡| 国产精品 欧美亚洲| 日韩精品免费视频一区二区三区| 国产精品 欧美亚洲| 午夜福利乱码中文字幕| 正在播放国产对白刺激| 少妇被粗大的猛进出69影院| 国产精品久久电影中文字幕 | 午夜福利,免费看| 搡老乐熟女国产| 亚洲视频免费观看视频| 青青草视频在线视频观看| 国产精品国产av在线观看| 色尼玛亚洲综合影院| 人妻 亚洲 视频| 777米奇影视久久| 在线观看一区二区三区激情| 少妇的丰满在线观看| 国产在线观看jvid| 欧美日韩亚洲国产一区二区在线观看 | 欧美性长视频在线观看| 国产成人精品久久二区二区91| 欧美精品啪啪一区二区三区| 51午夜福利影视在线观看| 欧美人与性动交α欧美软件| 亚洲av第一区精品v没综合| 不卡av一区二区三区| 国产有黄有色有爽视频| 成人永久免费在线观看视频 | 一二三四在线观看免费中文在| 国产精品 国内视频| 亚洲色图av天堂| 国产福利在线免费观看视频| 欧美精品啪啪一区二区三区| 亚洲精华国产精华精| 久热爱精品视频在线9| 久久久久国内视频| 一区二区三区激情视频| 久久久精品区二区三区| 麻豆国产av国片精品| 少妇粗大呻吟视频| 在线观看免费日韩欧美大片| 女性被躁到高潮视频| 天堂中文最新版在线下载| 一本久久精品| 亚洲精品粉嫩美女一区| 亚洲情色 制服丝袜| 91精品三级在线观看| 一级毛片女人18水好多| 亚洲人成77777在线视频| 9191精品国产免费久久| 天堂中文最新版在线下载| 五月天丁香电影| 69av精品久久久久久 | 69精品国产乱码久久久| 国产在线一区二区三区精| 国产成人免费观看mmmm| 男女无遮挡免费网站观看| 午夜成年电影在线免费观看| 亚洲男人天堂网一区| 精品国内亚洲2022精品成人 | 人成视频在线观看免费观看| 成年人午夜在线观看视频| av网站免费在线观看视频| 亚洲专区国产一区二区| 后天国语完整版免费观看| 免费在线观看影片大全网站| aaaaa片日本免费| 少妇精品久久久久久久| 欧美中文综合在线视频| 深夜精品福利| 91精品国产国语对白视频| 黄色视频,在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 国产不卡一卡二| www.精华液| 蜜桃在线观看..| 高潮久久久久久久久久久不卡| 在线天堂中文资源库| 热re99久久国产66热| 俄罗斯特黄特色一大片| 激情在线观看视频在线高清 | 最近最新中文字幕大全电影3 | 欧美在线一区亚洲| 天堂俺去俺来也www色官网| 别揉我奶头~嗯~啊~动态视频| 亚洲av成人不卡在线观看播放网| 成年女人毛片免费观看观看9 | 啦啦啦 在线观看视频| 色综合婷婷激情| 91av网站免费观看| 一二三四在线观看免费中文在| 色精品久久人妻99蜜桃| 麻豆成人av在线观看| 成年人黄色毛片网站| 91字幕亚洲| 日韩大片免费观看网站| 少妇被粗大的猛进出69影院| 国产99久久九九免费精品| videosex国产| 欧美激情久久久久久爽电影 | 午夜免费成人在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产免费av片在线观看野外av| 国产精品一区二区免费欧美| 亚洲中文av在线| 日韩免费高清中文字幕av| 久久亚洲真实| 婷婷成人精品国产| 日本vs欧美在线观看视频| 午夜福利欧美成人| 成人三级做爰电影| 国产又色又爽无遮挡免费看| 亚洲伊人久久精品综合| 久久久久久人人人人人| 国产精品久久久久久人妻精品电影 | 中文字幕精品免费在线观看视频| 久久久久精品人妻al黑| 国产成人av激情在线播放| 国产精品二区激情视频| 久久久久久久久免费视频了| 精品一品国产午夜福利视频| 免费久久久久久久精品成人欧美视频| 亚洲 国产 在线| 欧美日韩一级在线毛片| 每晚都被弄得嗷嗷叫到高潮| 久久免费观看电影| 女性被躁到高潮视频| 免费女性裸体啪啪无遮挡网站| 久久热在线av| 国产成人系列免费观看| 成人手机av| 欧美精品亚洲一区二区| 久久人妻熟女aⅴ| 亚洲专区字幕在线| 五月开心婷婷网| 在线观看免费高清a一片| 亚洲av美国av| 国产亚洲午夜精品一区二区久久| 国产日韩欧美在线精品| 亚洲成a人片在线一区二区| 18禁观看日本| 久9热在线精品视频| 午夜两性在线视频| 69av精品久久久久久 | 欧美日韩中文字幕国产精品一区二区三区 | 男女免费视频国产| 天堂8中文在线网| 亚洲精品久久成人aⅴ小说| 性少妇av在线| 男人舔女人的私密视频| 日本av免费视频播放| 一级黄色大片毛片| 少妇裸体淫交视频免费看高清 | 大香蕉久久成人网| 好男人电影高清在线观看| 精品久久蜜臀av无| 国产伦理片在线播放av一区| 99国产精品99久久久久| 亚洲国产精品一区二区三区在线| 动漫黄色视频在线观看| 免费不卡黄色视频| 久久久精品免费免费高清| 久久人人爽av亚洲精品天堂| 亚洲av电影在线进入| 亚洲午夜精品一区,二区,三区| 久久九九热精品免费| 自线自在国产av| 中文字幕色久视频| 欧美黄色片欧美黄色片| 色播在线永久视频| 母亲3免费完整高清在线观看| 久久婷婷成人综合色麻豆| 国产一区二区在线观看av| 亚洲国产看品久久| 久久国产精品人妻蜜桃| 免费高清在线观看日韩| 精品少妇一区二区三区视频日本电影| 国产极品粉嫩免费观看在线| 国产又色又爽无遮挡免费看| 精品少妇久久久久久888优播| 国内毛片毛片毛片毛片毛片| 国产视频一区二区在线看| 亚洲人成电影免费在线| 日韩 欧美 亚洲 中文字幕| 淫妇啪啪啪对白视频| 精品久久久久久久毛片微露脸| 精品福利观看| 日本一区二区免费在线视频| 精品视频人人做人人爽| 人人妻人人添人人爽欧美一区卜| 国产有黄有色有爽视频| 亚洲中文字幕日韩| 女人高潮潮喷娇喘18禁视频| 成人国语在线视频| 免费在线观看完整版高清| 午夜91福利影院| 一夜夜www| 精品一区二区三区视频在线观看免费 | 丝袜喷水一区| 搡老乐熟女国产| 香蕉丝袜av| 国产一区二区三区视频了| 大码成人一级视频| 中文字幕人妻丝袜一区二区| 色婷婷久久久亚洲欧美| 国产精品 国内视频| 精品国产亚洲在线| 欧美日韩视频精品一区| 久久国产精品人妻蜜桃| 国产日韩欧美在线精品| 久久人妻熟女aⅴ| 亚洲欧美精品综合一区二区三区| 大香蕉久久成人网| 国产高清视频在线播放一区| 别揉我奶头~嗯~啊~动态视频| 午夜福利视频在线观看免费| 午夜成年电影在线免费观看| 亚洲美女黄片视频| 高清av免费在线| 一本一本久久a久久精品综合妖精| 真人做人爱边吃奶动态| 狠狠狠狠99中文字幕| 国产日韩欧美视频二区| 99香蕉大伊视频| 在线观看人妻少妇| 久久久久久亚洲精品国产蜜桃av| 国产日韩欧美在线精品| 国产精品久久电影中文字幕 | 久久 成人 亚洲| 国产精品自产拍在线观看55亚洲 | av天堂久久9| 国产精品国产av在线观看| 看免费av毛片| 欧美日韩亚洲综合一区二区三区_| 欧美日韩精品网址| 男人操女人黄网站| 欧美日韩黄片免| 日韩一卡2卡3卡4卡2021年| 国产日韩欧美在线精品| 国产激情久久老熟女| 久久av网站| 久久国产亚洲av麻豆专区| 黄频高清免费视频| 黄色成人免费大全| 后天国语完整版免费观看| 欧美 日韩 精品 国产| 国产麻豆69| av在线播放免费不卡| 青青草视频在线视频观看| 一个人免费在线观看的高清视频| 色老头精品视频在线观看| 国产精品成人在线| 精品人妻1区二区| 熟女少妇亚洲综合色aaa.| 久久中文看片网| 9191精品国产免费久久| 男女边摸边吃奶| 国产欧美日韩一区二区三区在线| 欧美一级毛片孕妇| 亚洲久久久国产精品| 一本综合久久免费| 中国美女看黄片| 亚洲全国av大片| 老司机午夜十八禁免费视频| 欧美在线黄色| 制服人妻中文乱码| 女人久久www免费人成看片| 国产精品久久久久久精品电影小说| 日本av免费视频播放| 国产成人系列免费观看|