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

    熱風(fēng)供暖送風(fēng)方式對高大空間室內(nèi)污染物擴散的影響

    2015-08-24 08:53:12亢燕銘
    關(guān)鍵詞:溫差射流氣流

    葉 筱, 亢燕銘, 左 濱, 鐘 珂

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

    熱風(fēng)供暖送風(fēng)方式對高大空間室內(nèi)污染物擴散的影響

    葉筱, 亢燕銘, 左濱, 鐘珂

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

    利用數(shù)值模擬方法研究分析了混合通風(fēng)和碰撞射流通風(fēng)供暖時,高大空間內(nèi)的熱環(huán)境及污染物濃度分布特征.結(jié)果表明,兩種熱風(fēng)供暖方式在高大空間內(nèi)形成的氣流形態(tài)完全不同,導(dǎo)致污染物的空間分布也有很大差別.碰撞射流通風(fēng)時,不同類型污染物(CO2和2.5 μm顆粒物)在房間內(nèi)都能均勻分布,而混合通風(fēng)時,兩者在房間內(nèi)的分布極不均勻且表現(xiàn)出不同的分布特征,CO2集中分布于房間中部,2.5 μm顆粒物則集中分布在遠離外墻的近地面區(qū).碰撞射流通風(fēng)比混合通風(fēng)更有利于污染物的排除.此外,送風(fēng)溫差對混合通風(fēng)房間內(nèi)污染物的分布特征有明顯影響,但對碰撞射流通風(fēng)基本沒有影響.

    高大空間; 熱風(fēng)供暖; 混合通風(fēng); 碰撞射流通風(fēng); 主動式污染物; 被動式污染物

    通風(fēng)形式是影響房間內(nèi)污染物擴散的一個重要因素[1].傳統(tǒng)的冬季熱風(fēng)供暖房間,通常采用混合通風(fēng)(MV)系統(tǒng)中上送上回式的氣流組織形式.由于熱浮力的作用,送風(fēng)熱氣流將難以到達人員空間,導(dǎo)致供暖能量利用率偏低[2].碰撞射流通風(fēng)(IJV)作為一種新型的通風(fēng)方式,完全克服了混合通風(fēng)供暖時熱風(fēng)難以抵達人員空間的缺點.碰撞射流通風(fēng)系統(tǒng)中送風(fēng)氣流以較高的動量從垂直送風(fēng)口射出,氣流撞擊地板后向四周蔓延,冬季熱風(fēng)供暖時送風(fēng)氣流在慣性力的作用下與室內(nèi)空氣充分混合,供暖能量利用率明顯高于混合通風(fēng)[3-5].然而,由于以上兩種送風(fēng)方式在室內(nèi)形成的流場不同,室內(nèi)污染物擴散和分布特征也可能有很大區(qū)別.碰撞射流通風(fēng)能否提供優(yōu)于混合通風(fēng)房間的空氣品質(zhì),是其進一步推廣的關(guān)鍵.

    文獻[6-7]研究表明,房間內(nèi)污染物的運動主要受慣性、自然對流、溫度梯度及人員活動等因素影響.文獻[1]用數(shù)值模擬的方法研究了置換通風(fēng)和混合通風(fēng)房間內(nèi)氣溶膠粒子的沉積與濃度分布,結(jié)果表明,粒子的沉積及濃度分布主要受房間通風(fēng)條件的影響.文獻[8]指出,送、排風(fēng)口布置對室內(nèi)CO2的排出效果有重要影響.

    盡管關(guān)于碰撞射流通風(fēng)和混合通風(fēng)空間的溫度和流場分布的研究較多[2-5],但對高大空間在這兩種送風(fēng)方式下因人員呼吸產(chǎn)生的CO2(主動式污染物)和人員行走產(chǎn)生的顆粒物(被動式污染物)擴散特征的研究較少.本文在對混合通風(fēng)和碰撞射流通風(fēng)高大空間內(nèi)的熱環(huán)境及污染物濃度分布進行數(shù)值模擬的基礎(chǔ)上,分析比較兩種通風(fēng)方式下不同性質(zhì)污染物的濃度空間分布規(guī)律,為碰撞射流的推廣應(yīng)用提供理論依據(jù).

    1 數(shù)值計算模型

    1.1物理模型及邊界條件

    以一個層高為10 m的高大空間作為物理模型進行研究,圖1給出了房間的柱網(wǎng)布置情況,每個柱子的橫截面尺寸為0.6 m × 0.6 m,相鄰柱距間空調(diào)送、回風(fēng)口布置相同,根據(jù)對稱性理論,本文選取一個柱距單元作為研究對象,其平面分布見圖1灰色部分.采用MV時,送、回風(fēng)口均設(shè)置在外墻上部,送風(fēng)口尺寸為0.3 m × 0.3 m,回風(fēng)口尺寸為0.3 m × 0.2 m,送風(fēng)口高度為8.5 m,回風(fēng)口設(shè)在送風(fēng)口下方0.5 m處; 采用IJV時,送風(fēng)管緊貼柱子,送風(fēng)口尺寸為0.3 m × 0.2 m,距離地面的高度為0.3 m,回風(fēng)口尺寸為0.3 m × 0.3 m,位于屋頂中央.另外,房間內(nèi)設(shè)有3組站立的人群,每組人群尺寸為0.8 m × 0.8 m × 1.7 m(高).兩種送風(fēng)方式的風(fēng)口布置形式以及人員所在的空間位置如圖2所示.

    圖1 物理模型平面分布圖(單位:mm)Fig.1 Plan of the space researched(unit: mm)

    a—MV送風(fēng)口; b—MV回風(fēng)口; c,d—IJV送風(fēng)口; e—IJV回風(fēng)口; p—人員所在位置圖2 高大空間送、回風(fēng)口布置及典型面分布圖Fig.2 Locations of the inlet and outlet and typical surface in the large-height space

    研究空間采用非結(jié)構(gòu)化網(wǎng)格進行劃分,送、回風(fēng)口、近壁面處以及污染源附近的網(wǎng)格均進行加密處理,最小網(wǎng)格尺寸為0.03 m,總網(wǎng)格數(shù)為(310~400)萬.

    本文將外墻作為熱損失面.在模擬過程中,地板和屋頂均設(shè)為絕熱面.室內(nèi)熱源如人、燈和電腦等雖對房間供熱有利,但對氣流和溫度的局部分布影響很小[9],故在模擬中不考慮室內(nèi)熱源的影響.由于送風(fēng)空調(diào)房間通常保持正壓,所以在模擬計算時將不考慮冷風(fēng)滲透造成的熱損失.

    在不同熱負荷條件下,室內(nèi)的溫度和流場將出現(xiàn)不同的形態(tài).本文以室外溫度to為-5和4 ℃為例分別進行模擬計算,采用定風(fēng)量系統(tǒng),通過改變送風(fēng)溫差來應(yīng)對不同室外氣候條件下的負荷要求.根據(jù)民用建筑供暖通風(fēng)與空調(diào)設(shè)計規(guī)范[10],冬季室內(nèi)設(shè)計溫度取21 ℃.MV中熱風(fēng)以水平向下30°傾角射出,送風(fēng)速度vs為3.0 m/s.文獻[3]推薦IJV的送風(fēng)速度為1.4~3.0 m/s,因此,設(shè)熱風(fēng)以2.4 m/s的速度垂直向下射出.本文所計算工況的送風(fēng)參數(shù)和工作區(qū)(距離地面2 m以下空間)的平均溫度t列于表1中.

    表1 兩種送風(fēng)方式的送風(fēng)參數(shù)及工作區(qū)平均溫度Table 1 Supply parameters and temperature in occupied zones in the two air supply systems

    本文以人員呼出的CO2和人員行走產(chǎn)生的2.5 μm 顆粒物分別作為主動式污染物和被動式污染物的代表.模擬過程中,設(shè)室內(nèi)每組有4名人員,共3組,平均人員密度為0.3 人/m2.站立人員的CO2散發(fā)強度為0.022 m3/(h·人)[11],因此每組人群的CO2散發(fā)強度為0.088 m3/h.設(shè)人員行走引起的2.5 μm顆粒物全部來源于地面.文獻[12]的實測結(jié)果表明,每個人行走引起1~5 μm顆粒物的濃度增加1~2 μg/(m3·s),因此本文中設(shè)顆粒物散發(fā)強度為1.5 μg/(m3·s),密度設(shè)為1450 kg/m3.設(shè)CO2的送風(fēng)濃度為1.2 × 106μg/m3,顆粒物的送風(fēng)濃度為0 μg/m3.

    1.2數(shù)學(xué)模型

    本文所研究流體為三維連續(xù)不可壓縮流,采用Reynolds時均N-S方程(RANS)模擬計算室內(nèi)氣流的湍流流動,湍流模型選用RNGκ-ε模型.在研究過程中假設(shè)流體屬性不變,考慮到浮升力的影響,空氣的密度采用Boussinesq近似.數(shù)值計算時,采用SIMPLE算法控制壓力和速度的耦合.對于模型的離散化,除壓力項外其他項都選用二階迎風(fēng)格式,所有固體表面均設(shè)為無滲透和無滑移條件,送風(fēng)口邊界類型定義為velocity-inlet,回風(fēng)口定義為outflow,以保證進、出口質(zhì)量流量相等.

    目前主要有兩種方法處理多相流:歐拉-拉格朗日法和歐拉-歐拉法.歐拉-拉格朗日法用來跟蹤和分析粒子的運動軌跡,研究單個粒子的動態(tài)特性,而歐拉-歐拉法用來預(yù)測所有粒子的空間分布特征,以整個流場作為研究對象[13].本文的研究重點是污染物的空間分布規(guī)律,故采用歐拉-歐拉法中的混合模型進行模擬計算.模型中的粒子被處理成連續(xù)項,通過求解粒子與空氣之間的相對滑移速度方程確定重力對顆粒物分布的影響,本文通過引入用戶自定義函數(shù)(UDF)來確定粒子在不同方向壁面上的沉積量Jd,dA[13],如式(1)所示.

    Jd,dA=vd,dA·Cb,dA·dA

    (1)

    式中:Cb,dA為靠近壁面控制體中粒子相的體積分數(shù);dA為控制體壁面的面積;vd,dA為某一點粒子的沉積速度,其在不同壁面方向上的具體表達式見文獻[13].

    1.3數(shù)學(xué)模型的合理性驗證

    為了保證數(shù)值模擬方法的可靠性,需要驗證1.2 節(jié)中數(shù)學(xué)模型的合理性.首先對其流場及溫度場進行驗證.在東華大學(xué)環(huán)境科學(xué)與工程學(xué)院人工氣候?qū)嶒炇疫M行熱風(fēng)供暖實驗,對室內(nèi)的溫度場和速度場分布進行實測.人工氣候?qū)嶒炇夜膊贾?個測桿,每個測桿沿高度方向設(shè)置7個測點,測點布置情況如圖3所示.實驗時采用Humlog 20型溫度儀測量溫度,其分辨率為0.1 ℃,采用意大利Delta公司的萬向風(fēng)速儀來測量氣流速度,其測量范圍為0.05~5.00 m/s,精度為±3%,分辨率為0.01 m/s. 實驗過程中,為保證測量結(jié)果的穩(wěn)定,連續(xù)測量5~6 h,直到室內(nèi)空氣溫度和壁面溫度的波動范圍小于±5 ℃, 實驗數(shù)據(jù)均取自穩(wěn)定狀態(tài)下的實測值.

    圖3 人工氣候室測點布置情況(單位:mm)Fig.3 Arrangement diagram of monitoring point in the test chamber(unit: mm)

    針對上述實驗工況,本文采用數(shù)值計算的方法對其進行模擬,通過對比模擬結(jié)果和實測值,驗證該數(shù)學(xué)模型的合理性.圖4給出了1號和3號測桿數(shù)值模擬值與實測值的比較.由于篇幅限制,其他測桿的對比結(jié)果不再給出.

    (a) 1號測桿

    (b) 3號測桿

    由圖4可知,除某些位置數(shù)值模擬結(jié)果與實驗數(shù)據(jù)之間產(chǎn)生一些偏差外,整體上模擬值與實驗值能夠較好地吻合.這可能是由于門的開啟以及實驗過程中冷風(fēng)滲透等原因?qū)κ覂?nèi)的溫度場和速度場造成影響,進而導(dǎo)致模擬結(jié)果與實驗結(jié)果之間的偏差.可以認為本文所采用的計算流場和溫度場的數(shù)值方法是正確的.

    這里分別以文獻[14-15]的實驗參數(shù)(包括房間幾何參數(shù)、送風(fēng)參數(shù)以及污染物參數(shù)等),利用本文采用的污染物濃度分布的計算方法對二者的實驗及模擬結(jié)果進行了再現(xiàn).圖5是本文模擬結(jié)果與文獻[14]中給出的測點1和測點6處SF6濃度分布的比較. 圖5中,H*= 測點高度/房間高度,C*= (測點濃度-送風(fēng)濃度)/(排風(fēng)濃度-送風(fēng)濃度).由圖5可知,本文的模擬值比文獻[14]的模擬值更接近實測值,這充分說明本文所采用的計算主動式污染物濃度分布的方法是可行可靠的.

    圖6是文獻[15]給出的顆粒物濃度分布結(jié)果與本文模擬值的對比.圖6中,C= 測點濃度/排風(fēng)濃度.由圖6可以看出,除了某些位置與文獻中的結(jié)果有較小的偏差外,其他位置處的模擬值與文獻[15]給出的結(jié)果基本上能夠較好地吻合,可以認為本文所采用的計算顆粒物濃度分布的數(shù)值方法是合理的.

    (a) 測點1

    (b) 測點6

    (a) 測點3

    (b) 測點5

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

    對于定風(fēng)量系統(tǒng),房間供暖熱負荷越大,送風(fēng)溫差越大,送風(fēng)氣流受到的熱浮力作用越大,這可能影響到室內(nèi)的流場以及污染物分布特征.同時主動式污染物和被動式污染物的初始擴散速度不同造成二者的輸運特征有所區(qū)別.因此有必要對不同負荷條件下,研究分析MV和IJV兩種送風(fēng)方式在高大空間內(nèi)的溫度場和速度場,并對不同性質(zhì)污染物的濃度空間分布規(guī)律進行比較.

    2.1兩種送風(fēng)形式下室內(nèi)的溫度場與速度場分布

    圖7和8分別是兩種送風(fēng)方式在不同室外溫度條件下房間中心剖面也即過人體剖面(x=2.25 m) 上的氣流速度和流場分布.

    圖7 混合通風(fēng)房間x = 2.25 m平面的流場分布Fig.7 Air flow fields in the MV at the plane of x = 2.25 m

    圖8 碰撞射流通風(fēng)房間x = 2.25 m平面的流場分布Fig.8 Air flow fields in the IJV at the plane of x = 2.25 m

    由圖7可知,在MV中,熱風(fēng)射流以30°傾角向下射出后,在浮升力的作用下運動軌跡發(fā)生改變,在送風(fēng)氣流末端,送風(fēng)動量大幅度衰減,最終在房間近中心處形成明顯的上升氣流,送風(fēng)氣流上升至頂棚附近時,沿著頂棚向四周擴散,熱氣流難以到達下部空間,房間溫度出現(xiàn)上高下低的特征.比較圖7(a)和7(b)可以看出,送風(fēng)溫度越低,送風(fēng)熱氣流越容易到達房間下部,在房間上部形成的渦流范圍越大.這是因為送風(fēng)溫度較低時,熱浮力對熱氣流的抬升作用較弱,使得氣流向下輸運距離較遠,下部空間得到的熱量較多.

    由圖8可以看到,對于IJV,送風(fēng)熱氣流以較高動量從垂直送風(fēng)口向下射出并撞擊地面,隨后沿著地面水平擴散一定距離后在浮力的作用下脫離地面向上運動,形成了覆蓋房間大部分空間的大渦流,這使得房間內(nèi)的氣流在慣性力的主導(dǎo)作用下充分混合,從而使房間內(nèi)溫度均勻分布.對比圖8(a)和8(b)可知,送風(fēng)溫差較小時,上部空間渦流的水平影響范圍較寬.這是因為送風(fēng)溫差越小,熱風(fēng)受到的熱浮力作用越小,送風(fēng)熱氣流沿地面水平擴散的距離越大.

    對比圖7和8可以看到,MV中,近地面氣流以很低的速度曲折向上運動,而IJV中,近地面處氣流多以水平流動為主且速度較大.這種差別是由兩種送風(fēng)方式的送風(fēng)口位置決定的.此外,由圖7和8還可以看出,由于外墻內(nèi)表面溫度較低,兩種送風(fēng)方式下外墻附近的空氣均通過對流換熱被冷卻后下沉.但是,在相同室外溫度條件下,IJV冷空氣下沉速度明顯大于MV的情況,且在靠近風(fēng)口近地面處形成較大范圍的渦流區(qū).這是因為IJV房間內(nèi)空氣充分混合,外墻內(nèi)表面與周圍空氣的溫差較大而導(dǎo)致自然對流作用比MV強,同時,IJV熱風(fēng)在室內(nèi)形成的大范圍渦流(見圖8)在外墻附近的氣流方向與下沉冷氣流的運動方向一致,也增加了下沉冷空氣的速度.

    兩種熱風(fēng)供暖送風(fēng)方式的氣流形態(tài)差別,必然造成二者室內(nèi)溫度分布特征的不同.圖9給出了房間不同高度平面的平均溫度.由圖9可知,MV送風(fēng)方式時溫度在房間高度方向上產(chǎn)生明顯的分層現(xiàn)象,而IJV房間內(nèi)溫度分布比較均勻.這是由兩種送風(fēng)方式完全不同的氣流形態(tài)造成的(見圖7和8).由圖9還可以看出,MV供暖房間上下溫差隨著送風(fēng)溫差的減小而減小,兩種負荷時上下溫差分別為8和6 ℃.IJV供暖房間上下溫差受送風(fēng)溫差的影響很小,兩種負荷條件下均不足1 ℃.

    圖9 兩種送風(fēng)方式溫度沿高度方向的變化Fig.9 Vertical distribution of indoor temperature in the two air supply systems

    2.2兩種送風(fēng)方式下主動式污染物濃度空間分布規(guī)律

    圖10是兩種送風(fēng)方式在to= -5 ℃時x= 2.25 m 剖面上氣流形態(tài)與CO2質(zhì)量濃度分布圖.由圖10(a)可知,在MV中,CO2質(zhì)量濃度分布非常不均勻,上部的渦流阻止了CO2向上擴散,而下部空間的向上氣流使人體呼出的CO2難以向下運動,最終使CO2在房間中間高度處聚集.圖10(b)表明,由于IJV供暖房間內(nèi)空氣較充分地混合流動,CO2質(zhì)量濃度分布比較均勻,僅在渦流區(qū)出現(xiàn)小范圍的高濃度區(qū).

    圖10 x = 2.25 m剖面氣流形態(tài)與CO2質(zhì)量濃度分布情況Fig.10 Air flow patterns and CO2mass concentration distributions at x = 2.25 m plane

    圖11為不同高度平面上的CO2平均質(zhì)量濃度.由圖11可知,在MV送風(fēng)方式下CO2質(zhì)量濃度沿房間高度方向變化明顯,最高質(zhì)量濃度出現(xiàn)在高度2.5~ 3.0 m范圍內(nèi),這與圖10的結(jié)果是一致的.同時CO2的質(zhì)量濃度分布受送風(fēng)溫差的影響較大,送風(fēng)溫差越大,房間內(nèi)CO2質(zhì)量濃度越大.由圖11還可以看到,IJV送風(fēng)方式下CO2質(zhì)量濃度沿房間高度的變化雖然沒有混合通風(fēng)明顯,但在人員呼吸高度處即污染物散發(fā)位置處依然出現(xiàn)最大值.

    圖11 不同高度上CO2的平均質(zhì)量濃度Fig.11 Average mass concentrations of CO2in different heights

    通常用人員呼吸高度及工作區(qū)污染物的平均濃度分布情況來衡量工作區(qū)空氣品質(zhì)的優(yōu)劣,為此,圖12給出了不同送風(fēng)狀態(tài)下站立人體平均呼吸高度(z= 1.6 m)、2.0 m以下空間、整個空間內(nèi)CO2的平均質(zhì)量濃度以及排風(fēng)質(zhì)量濃度.圖12表明,在不同室外溫度條件下,IJV房間內(nèi)z= 1.6和2.0 m以下空間CO2平均質(zhì)量濃度基本相同,且與室外溫度為-5 ℃時MV房間的情況相近;而對于MV房間內(nèi)z= 1.6和2.0 m以下空間CO2質(zhì)量濃度隨送風(fēng)溫差的減小而減小.由圖12也可以看出,混合通風(fēng)中送風(fēng)溫差不僅影響CO2的空間分布特征,并且會改變2.0 m以下空間與整體空間CO2質(zhì)量濃度值的相對關(guān)系,即當送風(fēng)溫差較小時,整體空間的CO2質(zhì)量濃度大于2.0 m以下空間的值,當送風(fēng)溫差較大時,則出現(xiàn)相反的情況.這是因為送風(fēng)溫差越小,下部空間的空氣混合越均勻,越有利于CO2的擴散,CO2的平均質(zhì)量濃度越小.

    圖12 不同空間位置CO2的平均質(zhì)量濃度及排風(fēng)質(zhì)量濃度Fig.12 Average mass concentrations of CO2in different zones and the mass concentrations of CO2in exhaust air

    CO2排風(fēng)質(zhì)量濃度是影響污染物排除效率的一個重要參數(shù),由圖12可知,IJV中CO2的排風(fēng)質(zhì)量濃度與室內(nèi)平均質(zhì)量濃度接近,而MV的CO2排風(fēng)質(zhì)量濃度則遠小于室內(nèi)平均質(zhì)量濃度.

    2.3兩種送風(fēng)方式下被動式污染物濃度空間分布規(guī)律

    圖13是兩種送風(fēng)方式在to= -5 ℃時x= 2.25 m 平面上2.5 μm顆粒物的濃度及速度等值線圖.由圖13可知,IJV中2.5 μm顆粒物的濃度分布比較均勻,與CO2在房間中的分布情況基本一致,而MV中顆粒物主要集中在遠離外墻的房間下部區(qū)域.這是因為IJV房間中近地面處的氣流速度明顯高于MV,由地面產(chǎn)生的2.5 μm顆粒物在近地面處氣流的作用下更容易被攜帶至主流區(qū),在主流區(qū)內(nèi)慣性力對顆粒物起主要的誘導(dǎo)作用,從而使室內(nèi)的顆粒物均勻擴散,最終通過排風(fēng)口排出;而MV房間中近地面處氣流速度過小,約為0.006 m/s,不足以使2.5 μm顆粒物進入主流區(qū),另外靠近外墻附近下沉的冷空氣沿著地板表面水平擴散,在該氣流的作用下使顆粒物聚集在遠離外墻的近地面區(qū)域.

    圖13 x = 2.25 m剖面氣流形態(tài)與顆粒物的質(zhì)量濃度Cp分布圖Fig.13 Air flow pattern and particles mass concentration distributions at x = 2.25 m plane

    圖14是不同高度平面上2.5 μm顆粒物的平均質(zhì)量濃度.通過圖14可以看出,在不同熱負荷條件下,MV房間中2.5 μm顆粒物的平均質(zhì)量濃度沿高度的變化不同,當送風(fēng)溫差較小時,房間下部顆粒物的質(zhì)量濃度高于送風(fēng)溫差大的情況,在近地面處二者的顆粒物質(zhì)量濃度甚至相差4倍,而在房間上部則出現(xiàn)相反的結(jié)果.這是因為負荷小時送風(fēng)熱氣流更多地進入下部空間,將顆粒物攜帶至房間上部,并被及時排出.IJV房間中由于氣流充分混合且近地面處氣流速度較大,不同送風(fēng)溫差時2.5 μm顆粒物質(zhì)量濃度均隨高度變化不大.此外,由圖14可知,MV房間中在近地面處2.5 μm顆粒物的質(zhì)量濃度遠高于碰撞射流,這是因為IJV在近地面處的氣流速度遠大于MV.

    不同空間區(qū)域內(nèi)2.5 μm顆粒物的平均質(zhì)量濃度及排風(fēng)濃度大小關(guān)系如圖15所示.由圖15可知,對于MV,送風(fēng)溫差對顆粒物質(zhì)量濃度的空間分布有明顯影響,這與圖14的結(jié)果一致.在IJV中,不同送風(fēng)溫差下2.5 μm顆粒物的質(zhì)量濃度基本相等并且都遠小于混合通風(fēng)的情況.由此可以認為,IJV更有利于房間內(nèi)被動式污染物的擴散.

    圖15 不同空間顆粒物的平均質(zhì)量濃度Cp及排風(fēng)濃度Fig.15 Average mass concentrations of particles in different zones and the mass concentrations of particles in exhaust air

    3 兩種送風(fēng)方式下不同性質(zhì)污染物排污效率對比

    排污效率是衡量通風(fēng)系統(tǒng)性能優(yōu)劣的一個重要參數(shù),其大小反映了空調(diào)通風(fēng)系統(tǒng)排除污染物的能力.在穩(wěn)態(tài)通風(fēng)情況下,排污效率εp[16]為

    (2)

    式中:Ce、Cs和Ch分別為房間內(nèi)回風(fēng)口、送風(fēng)口和任一空間點污染物的質(zhì)量濃度,μg/m3.本文數(shù)值模擬中,CO2和顆粒物的送風(fēng)質(zhì)量濃度分別為1.2×106和0 μg/m3, 回風(fēng)口處以及空間任一點的質(zhì)量濃度可以由模擬計算結(jié)果得到.

    不同高度平面CO2和2.5 μm顆粒物的平均排污效率如圖16所示.由圖16(a)可以看出,IJV中CO2的平均排污效率高于MV,在房間下部近地面處排污效率甚至出現(xiàn)了大于1的情況,并且送風(fēng)溫差越小排污效率越大.由圖16(b)可知,兩種送風(fēng)方式對2.5 μm顆粒物的排除能力有明顯不同,MV中在不同高度平面上排污效率變化很大,在房間下部排污效率遠小于1,而在房間上部出現(xiàn)接近1的情況,IJV中不同高度上排污效率均接近于1.同時可以看出,兩種送風(fēng)方式下近地面處顆粒物的排污效率均低于房間上部,這是由于顆粒物自身重力作用的結(jié)果.總的來看,IJV比MV更有利于污染物的排除.

    (a) CO2排污效率

    (b) 顆粒物排污效率

    4 結(jié) 語

    本文利用數(shù)值模擬的方法研究分析了混合通風(fēng)供暖和碰撞射流通風(fēng)供暖時,高大空間內(nèi)的熱環(huán)境及污染物濃度分布特征,得到以下主要結(jié)論.

    (1) 混合通風(fēng)供暖時,熱浮力使得送風(fēng)氣流主要在高大空間上部循環(huán);碰撞射流通風(fēng)供暖時,直接進入人員空間的熱空氣在整個空間內(nèi)循環(huán)流動.兩種送風(fēng)方式不同的氣流形態(tài)導(dǎo)致室內(nèi)溫度沿高度分布特征完全不同,混合通風(fēng)呈明顯的下低上高特征,碰撞射流氣溫分布非常均勻.

    (2) CO2和2.5 μm顆粒物在碰撞射流和混合通風(fēng)中的分布有很大不同,碰撞射流通風(fēng)時兩者在房間內(nèi)均勻分布,而采用混合通風(fēng)時污染物質(zhì)量濃度分布很不均勻,CO2集中分布于房間中部,2.5 μm顆粒物聚集在遠離外墻的近地面區(qū)域.

    (3) 送風(fēng)溫差對碰撞射流通風(fēng)污染物的分布影響不大,而對混合通風(fēng)有明顯的影響,送風(fēng)溫差越小,混合通風(fēng)工作區(qū)內(nèi)CO2的質(zhì)量濃度越小,2.5 μm顆粒物的質(zhì)量濃度越大.

    (4) 對于主動式污染物,采用碰撞射流通風(fēng)時排污效率在近地面處甚至可能大于1.對于被動式污染物,碰撞射流通風(fēng)的排污效率在整個空間都基本上接近于1,而混合通風(fēng)在人員空間的排污效率小于0.3.由此表明碰撞射流通風(fēng)比混合通風(fēng)更有利于污染物的排除.

    [1] 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.

    [2] 亢燕銘,張艷霞,鐘珂.落地窗對熱風(fēng)供暖房間室內(nèi)熱環(huán)境的影響[J].暖通空調(diào),2014,44(4):104-108.

    [3] KARIMIPANAH T, AWBI H B. Theoretical and experimental investigation of impinging jet ventilation and comparison with wall displacement ventilation[J]. Building and Environment, 2002, 37(12):1329-1342.

    [4] CHEN H J, MOSHFEGH B, CEHLIN M.Investigation on the flow and thermal behavior of impinging jet ventilation systems in an office with different heat loads[J]. Building and Environment, 2013, 59:127-144.

    [5] VARODOMPUN J, NAVVAB M. HVAC ventilation strategies: The contribution for thermal comfort, energy efficiency, and indoor air quality[J]. Journal of Green Building, 2007, 2(2): 131-150.

    [6] MORAWSKA L, SALTHAMMER T. Indoor environment[J]. Weinhein: Wiley-VCH, 2003.

    [7] WILLIAM W N, GLEN R C. Mass-transport aspects of pollutant removal at indoor surfaces[J]. Environment International, 1989, 15(1/2/3/4/5/6):567-584.

    [8] CHUNG K C, HSU S P. Effect of ventilation pattern on room air and contaminant distribution[J]. Building and Environment, 2001, 36 (9):989-998.

    [9] COSTA J J, OLIVEIRA L A, BLAY D. Turbulent airflow in a room with a two-jet heating-ventilation system: A numerical parametric study[J]. Energy and Building, 2000, 32(3): 327-343.

    [10] GB 50736—2012 民用建筑供暖通風(fēng)與空調(diào)設(shè)計規(guī)范[S].北京:中國建筑工業(yè)出版社,2012.

    [11] 趙榮義,范存養(yǎng),薛殿華,等.空氣調(diào)節(jié)[M].4版.北京:中國建筑工業(yè)出版社,2009:116.

    [12] LUOMA M, BATTERMAN S A. Characterization of particulate emissions from occupant activities in offices[J]. Indoor Air, 2001, 11(1): 35-48.

    [13] ZHAO B, YANG C Q, YANG X D, et al. Particle dispersion and deposition in ventilated rooms: Testing and evaluation of

    different Eulerian and Lagrangian models[J]. Building and Environment, 2008, 43(4): 388-397.

    [14] YUAN X, CHEN Q, GLICKSMAN L R, et al. Measurements and computations of room airflow with displacement ventilation[J]. ASHRAE Transactions, 1999, 105(1): 340-352.

    [15] ZHANG Z, CHEN Q. Experimental measurements and numerical simulations of particle transport and distribution in ventilated rooms[J]. Atmospheric Environment, 2006, 40(18): 3396-3408.

    [16] 朱穎心.建筑環(huán)境學(xué)[M].2版.北京:中國建筑工業(yè)出版社,2005:201-202.

    Effects of the Warm-Air Heating Supply Modes on Pollutant Diffusion in Large-Height Spaces

    YEXiao,KANGYan-ming,ZUOBin,ZHONGKe

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

    The indoor thermal environment and contaminant concentration distribution are studied numerically in large-height spaces under mixing ventilation (MV) and impinging jet ventilation (IJV). The results show that the indoor air flow fields are different in the two ventilation strategies when they are used for room heating, and then results in a considerable difference in the spatial distribution of contaminants. The two types of pollutants (CO2and particles with size of 2.5 μm) are seemed to be uniformly distributed in the IJV, while the concentration field in the MV is different, i.e., the CO2is mainly distributed in the center of the room and the particles are distributed near the floor away from the external wall. The results also show that the IJV is more efficient for pollutant removal than the MV. In addition, the air supply temperature difference has a strong influence on the distribution characteristics of contaminations in the MV, and this effect can be neglected in the IJV.

    large-height space; warm-air heating; mixing ventilation; impinging jet ventilation; active pollutant; passive pollutant

    1671-0444(2015)06-0821-08

    2014-10-11

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

    葉筱(1990—),女,河南南陽人,博士研究生,研究方向為建筑環(huán)境與節(jié)能.E-mail: jhyexiao@126.com

    鐘珂(聯(lián)系人),女,教授,E-mail: zhongkeyx@dhu.edu.cn

    TU 831.3

    A

    猜你喜歡
    溫差射流氣流
    手表與溫差
    深海逃逸艙射流注水均壓過程仿真分析
    低壓天然氣泄漏射流擴散特性研究
    煤氣與熱力(2022年4期)2022-05-23 12:45:00
    氣流的威力
    北方冬季養(yǎng)豬當心“四差”
    溫差“催甜”等
    低級熱溫差發(fā)電器的開發(fā)應(yīng)用研究
    電子制作(2018年23期)2018-12-26 01:01:20
    固體運載火箭變軌發(fā)動機噴管氣流分離研究
    飛片下的空氣形成的“超強高速氣流刀”
    基于停留時間分布的氣流床氣化爐通用網(wǎng)絡(luò)模型
    性插视频无遮挡在线免费观看| 少妇熟女aⅴ在线视频| 12—13女人毛片做爰片一| 久久精品久久久久久噜噜老黄 | 观看美女的网站| 欧美乱妇无乱码| 亚洲美女搞黄在线观看 | 国产激情偷乱视频一区二区| 悠悠久久av| 97超视频在线观看视频| 性色avwww在线观看| 亚洲avbb在线观看| 黄色一级大片看看| 欧美bdsm另类| 日本 欧美在线| 一级av片app| 亚洲成av人片在线播放无| 999久久久精品免费观看国产| 黄色女人牲交| 天天躁日日操中文字幕| 黄色丝袜av网址大全| 激情在线观看视频在线高清| 国产高清激情床上av| 亚洲欧美日韩卡通动漫| 性欧美人与动物交配| 国产亚洲精品综合一区在线观看| 午夜福利在线在线| 国产精品女同一区二区软件 | 在线国产一区二区在线| 国产三级在线视频| 最近中文字幕高清免费大全6 | 国内久久婷婷六月综合欲色啪| 一个人看视频在线观看www免费| 99国产极品粉嫩在线观看| 麻豆国产97在线/欧美| 午夜精品在线福利| 高潮久久久久久久久久久不卡| 日本黄色视频三级网站网址| 亚洲18禁久久av| 俺也久久电影网| 精品人妻一区二区三区麻豆 | 尤物成人国产欧美一区二区三区| 国产精品1区2区在线观看.| 男人舔女人下体高潮全视频| 久久亚洲真实| 黄色丝袜av网址大全| av福利片在线观看| 亚州av有码| 免费观看的影片在线观看| 老女人水多毛片| 三级男女做爰猛烈吃奶摸视频| 日本 av在线| 欧美激情在线99| 成年女人毛片免费观看观看9| 日日摸夜夜添夜夜添av毛片 | 淫秽高清视频在线观看| 久久精品国产99精品国产亚洲性色| 国产亚洲精品av在线| 波多野结衣高清无吗| 九九热线精品视视频播放| 无人区码免费观看不卡| 在线观看舔阴道视频| 亚洲精品一区av在线观看| 亚洲无线在线观看| 99热这里只有是精品在线观看 | 久久午夜福利片| 欧美高清性xxxxhd video| 黄色一级大片看看| 美女黄网站色视频| 日韩亚洲欧美综合| 国产色婷婷99| 国产av在哪里看| 日韩欧美免费精品| 中文字幕免费在线视频6| 少妇高潮的动态图| 最新中文字幕久久久久| 老鸭窝网址在线观看| 特大巨黑吊av在线直播| 国产精品一区二区三区四区免费观看 | 一卡2卡三卡四卡精品乱码亚洲| 黄色配什么色好看| 久久中文看片网| 亚洲最大成人手机在线| 窝窝影院91人妻| 日韩免费av在线播放| 午夜福利免费观看在线| 91麻豆av在线| 精品不卡国产一区二区三区| 亚洲熟妇中文字幕五十中出| 国产亚洲欧美98| 国产成+人综合+亚洲专区| 无遮挡黄片免费观看| 听说在线观看完整版免费高清| 欧美色欧美亚洲另类二区| 国产成+人综合+亚洲专区| 欧美在线黄色| 一区二区三区免费毛片| 99热精品在线国产| 人人妻,人人澡人人爽秒播| 女生性感内裤真人,穿戴方法视频| 麻豆成人午夜福利视频| 欧美激情国产日韩精品一区| 国产私拍福利视频在线观看| 国内精品久久久久久久电影| 日韩精品中文字幕看吧| 欧美黑人欧美精品刺激| 很黄的视频免费| 中文资源天堂在线| 国产蜜桃级精品一区二区三区| 男插女下体视频免费在线播放| 精品一区二区三区av网在线观看| 男人舔奶头视频| 亚洲av.av天堂| 免费搜索国产男女视频| 熟女人妻精品中文字幕| 婷婷丁香在线五月| 观看免费一级毛片| 国产精品99久久久久久久久| 日韩免费av在线播放| 黄色视频,在线免费观看| 性欧美人与动物交配| 国产白丝娇喘喷水9色精品| 日韩欧美一区二区三区在线观看| 亚洲人成电影免费在线| 人妻制服诱惑在线中文字幕| 免费高清视频大片| aaaaa片日本免费| 国产在线精品亚洲第一网站| 日韩国内少妇激情av| 男人舔奶头视频| 亚洲av日韩精品久久久久久密| 国产成人欧美在线观看| 国产亚洲精品久久久久久毛片| 亚洲18禁久久av| 在线免费观看不下载黄p国产 | 亚洲精品亚洲一区二区| 老鸭窝网址在线观看| 国产成人a区在线观看| 91字幕亚洲| 日本免费a在线| 欧美日韩亚洲国产一区二区在线观看| 身体一侧抽搐| 18禁在线播放成人免费| 国产大屁股一区二区在线视频| 国内精品美女久久久久久| 90打野战视频偷拍视频| 日韩欧美在线二视频| 身体一侧抽搐| 黄色视频,在线免费观看| 精品99又大又爽又粗少妇毛片 | www.999成人在线观看| 日韩 亚洲 欧美在线| 内射极品少妇av片p| 精品一区二区免费观看| 精品一区二区三区视频在线| 色综合婷婷激情| 亚洲精品成人久久久久久| 成人无遮挡网站| av专区在线播放| 精华霜和精华液先用哪个| 欧美成人一区二区免费高清观看| 亚洲国产日韩欧美精品在线观看| 日韩欧美一区二区三区在线观看| 国产视频内射| 国产午夜福利久久久久久| 国产精品爽爽va在线观看网站| 亚洲av不卡在线观看| 99视频精品全部免费 在线| 成人特级av手机在线观看| 老司机午夜福利在线观看视频| 变态另类成人亚洲欧美熟女| 色哟哟·www| 国产成人欧美在线观看| 一区二区三区高清视频在线| 91久久精品国产一区二区成人| 真实男女啪啪啪动态图| 亚洲欧美日韩无卡精品| 最近在线观看免费完整版| 亚洲午夜理论影院| 欧美黄色片欧美黄色片| 亚洲av电影不卡..在线观看| 国产高潮美女av| 在线观看美女被高潮喷水网站 | 男女下面进入的视频免费午夜| 亚洲最大成人手机在线| 伦理电影大哥的女人| 一二三四社区在线视频社区8| 97超级碰碰碰精品色视频在线观看| 国语自产精品视频在线第100页| 十八禁网站免费在线| 色视频www国产| 国产老妇女一区| 色综合站精品国产| 真实男女啪啪啪动态图| 少妇丰满av| 99精品久久久久人妻精品| 欧美绝顶高潮抽搐喷水| 亚洲av日韩精品久久久久久密| 亚洲久久久久久中文字幕| 午夜福利高清视频| www.www免费av| 99久久99久久久精品蜜桃| 黄色女人牲交| 国产中年淑女户外野战色| 白带黄色成豆腐渣| 精品午夜福利在线看| 女人被狂操c到高潮| 在线天堂最新版资源| 黄色视频,在线免费观看| 国产在视频线在精品| 久久精品人妻少妇| 高清毛片免费观看视频网站| 亚洲国产日韩欧美精品在线观看| 九色成人免费人妻av| 校园春色视频在线观看| 婷婷精品国产亚洲av在线| 国产精品久久视频播放| 在线十欧美十亚洲十日本专区| 国产高清视频在线观看网站| 日本三级黄在线观看| 97碰自拍视频| a在线观看视频网站| 韩国av一区二区三区四区| 高清日韩中文字幕在线| 成年女人毛片免费观看观看9| 日本黄大片高清| 欧美zozozo另类| 日日摸夜夜添夜夜添小说| av欧美777| 中文字幕精品亚洲无线码一区| h日本视频在线播放| 久久久久久大精品| 人妻丰满熟妇av一区二区三区| www日本黄色视频网| 精品熟女少妇八av免费久了| 亚洲国产高清在线一区二区三| av在线天堂中文字幕| 国产又黄又爽又无遮挡在线| 午夜福利欧美成人| 两人在一起打扑克的视频| 欧美极品一区二区三区四区| 亚洲精品亚洲一区二区| 久9热在线精品视频| 欧美日韩瑟瑟在线播放| 深爱激情五月婷婷| 国产伦精品一区二区三区四那| 久久国产乱子免费精品| 麻豆成人午夜福利视频| 精品一区二区三区视频在线观看免费| 亚洲精华国产精华精| 亚洲午夜理论影院| 久久久久久国产a免费观看| 精品久久久久久久久久久久久| 日本与韩国留学比较| 亚洲最大成人手机在线| 99视频精品全部免费 在线| 日本免费a在线| 可以在线观看毛片的网站| 日本 av在线| 桃色一区二区三区在线观看| 日本免费一区二区三区高清不卡| 免费观看的影片在线观看| 欧美激情久久久久久爽电影| 美女xxoo啪啪120秒动态图 | 欧美在线一区亚洲| 内射极品少妇av片p| 国产一区二区激情短视频| 乱码一卡2卡4卡精品| 国产av在哪里看| 亚洲av熟女| 蜜桃久久精品国产亚洲av| 白带黄色成豆腐渣| 3wmmmm亚洲av在线观看| 亚洲精华国产精华精| 国产av一区在线观看免费| 国产老妇女一区| 精品久久久久久久久av| 亚洲人成伊人成综合网2020| 99国产精品一区二区三区| 午夜福利欧美成人| 在线观看午夜福利视频| 久久久精品大字幕| 午夜免费成人在线视频| 岛国在线免费视频观看| www.熟女人妻精品国产| 精品国产亚洲在线| 哪里可以看免费的av片| 亚洲欧美激情综合另类| 成人高潮视频无遮挡免费网站| 亚洲av.av天堂| 国产淫片久久久久久久久 | 国产黄色小视频在线观看| 赤兔流量卡办理| 国产精品一区二区三区四区久久| 又黄又爽又刺激的免费视频.| 日韩中字成人| 婷婷丁香在线五月| 亚洲精品成人久久久久久| 99热6这里只有精品| 国产午夜福利久久久久久| 精品久久久久久久末码| 伦理电影大哥的女人| 午夜亚洲福利在线播放| 国产精品亚洲美女久久久| 久久久久久久午夜电影| 天堂av国产一区二区熟女人妻| 美女高潮喷水抽搐中文字幕| 久久国产精品人妻蜜桃| 亚洲,欧美精品.| 国产久久久一区二区三区| 91av网一区二区| 老司机午夜十八禁免费视频| 悠悠久久av| 日本一二三区视频观看| 变态另类成人亚洲欧美熟女| 国产精品久久视频播放| 又紧又爽又黄一区二区| 日韩高清综合在线| 国产白丝娇喘喷水9色精品| 精品免费久久久久久久清纯| 午夜免费男女啪啪视频观看 | 91麻豆精品激情在线观看国产| 精品99又大又爽又粗少妇毛片 | 亚洲aⅴ乱码一区二区在线播放| 搡老熟女国产l中国老女人| 色综合站精品国产| 亚洲av中文字字幕乱码综合| 国产一区二区在线av高清观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩有码中文字幕| 成人毛片a级毛片在线播放| 午夜精品在线福利| 亚洲精品在线美女| 久久久色成人| 老熟妇仑乱视频hdxx| 禁无遮挡网站| 国产精品一区二区性色av| 久久性视频一级片| 身体一侧抽搐| 精品熟女少妇八av免费久了| 日韩免费av在线播放| 亚洲国产色片| 久久久久久九九精品二区国产| 国产亚洲精品久久久com| 亚洲专区国产一区二区| 久久久久国内视频| 精品国产亚洲在线| 久久久国产成人免费| 婷婷精品国产亚洲av| 亚洲国产精品成人综合色| 成人高潮视频无遮挡免费网站| 午夜亚洲福利在线播放| 国内精品一区二区在线观看| 免费人成视频x8x8入口观看| 亚洲av五月六月丁香网| 中亚洲国语对白在线视频| 成年版毛片免费区| 国模一区二区三区四区视频| 亚洲成人免费电影在线观看| 欧美成狂野欧美在线观看| 欧美区成人在线视频| 久久久久国内视频| 女同久久另类99精品国产91| 高清毛片免费观看视频网站| 国内精品美女久久久久久| 成年人黄色毛片网站| 亚洲成人免费电影在线观看| 亚洲欧美日韩高清专用| 中文字幕av成人在线电影| 欧美日韩福利视频一区二区| 午夜福利在线在线| 一级av片app| 丁香欧美五月| 蜜桃久久精品国产亚洲av| 热99在线观看视频| 久久人妻av系列| 国产精品一及| 精品国内亚洲2022精品成人| 成年女人永久免费观看视频| 国产精品不卡视频一区二区 | av天堂中文字幕网| 国产精华一区二区三区| 悠悠久久av| 熟女人妻精品中文字幕| 免费看a级黄色片| 欧美3d第一页| 美女免费视频网站| 精品午夜福利视频在线观看一区| 男女做爰动态图高潮gif福利片| 精品久久久久久成人av| 久久久久久久久大av| 男女那种视频在线观看| 国产av不卡久久| 免费一级毛片在线播放高清视频| 国产精品爽爽va在线观看网站| 9191精品国产免费久久| 欧美成人性av电影在线观看| 三级男女做爰猛烈吃奶摸视频| 最后的刺客免费高清国语| 国产老妇女一区| 国产亚洲av嫩草精品影院| 国内精品美女久久久久久| 别揉我奶头~嗯~啊~动态视频| 国内少妇人妻偷人精品xxx网站| 午夜激情欧美在线| 97人妻精品一区二区三区麻豆| 老司机福利观看| 国产不卡一卡二| 欧美高清性xxxxhd video| 亚洲综合色惰| 日日夜夜操网爽| 国产亚洲欧美在线一区二区| 国产精品99久久久久久久久| 亚洲av熟女| 变态另类成人亚洲欧美熟女| 欧美色视频一区免费| 中文字幕av成人在线电影| 国产视频内射| 一区福利在线观看| 国产大屁股一区二区在线视频| 波多野结衣高清作品| 天堂av国产一区二区熟女人妻| 最近视频中文字幕2019在线8| 精品一区二区免费观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av免费在线观看| 免费在线观看日本一区| 变态另类丝袜制服| 国产伦在线观看视频一区| 岛国在线免费视频观看| 亚洲成人免费电影在线观看| 国产精品美女特级片免费视频播放器| 久久国产乱子伦精品免费另类| 一卡2卡三卡四卡精品乱码亚洲| 国产免费av片在线观看野外av| 99热这里只有是精品50| a级毛片免费高清观看在线播放| 一本综合久久免费| 最近最新免费中文字幕在线| 麻豆久久精品国产亚洲av| 熟女人妻精品中文字幕| 亚洲国产精品999在线| 欧美高清性xxxxhd video| 黄色女人牲交| 人妻丰满熟妇av一区二区三区| 亚洲熟妇中文字幕五十中出| 亚洲avbb在线观看| 日本免费一区二区三区高清不卡| 久久久久国内视频| 成人高潮视频无遮挡免费网站| 九九热线精品视视频播放| 国产综合懂色| 搡女人真爽免费视频火全软件 | 国产极品精品免费视频能看的| 丝袜美腿在线中文| 特级一级黄色大片| 成人特级av手机在线观看| 国产成人a区在线观看| 久久这里只有精品中国| 99久国产av精品| 自拍偷自拍亚洲精品老妇| 无人区码免费观看不卡| 悠悠久久av| 国产精华一区二区三区| 精品不卡国产一区二区三区| 免费人成视频x8x8入口观看| 757午夜福利合集在线观看| 国产黄色小视频在线观看| 波野结衣二区三区在线| 国产欧美日韩一区二区精品| 午夜福利18| 亚洲久久久久久中文字幕| 亚洲av免费在线观看| 人人妻人人看人人澡| 欧美成人性av电影在线观看| 精品一区二区三区视频在线观看免费| 国内精品久久久久久久电影| 婷婷六月久久综合丁香| 成人欧美大片| 国产黄色小视频在线观看| 日本成人三级电影网站| 少妇人妻精品综合一区二区 | 久久99热6这里只有精品| 国产不卡一卡二| 亚洲av成人精品一区久久| 成年女人永久免费观看视频| 久久久久国产精品人妻aⅴ院| a在线观看视频网站| 欧美+亚洲+日韩+国产| 最新在线观看一区二区三区| 亚洲av.av天堂| 小蜜桃在线观看免费完整版高清| 亚洲av熟女| 欧美日韩黄片免| 最新在线观看一区二区三区| 国产精品一及| 国产亚洲精品av在线| 亚洲五月天丁香| 午夜福利高清视频| 校园春色视频在线观看| 亚洲熟妇熟女久久| 日韩国内少妇激情av| 嫁个100分男人电影在线观看| 国产精品综合久久久久久久免费| 国产视频内射| 免费电影在线观看免费观看| netflix在线观看网站| 久久久久久久精品吃奶| 日韩欧美 国产精品| 如何舔出高潮| 男女床上黄色一级片免费看| 91在线观看av| 高清毛片免费观看视频网站| 国产男靠女视频免费网站| 精品久久久久久久久久久久久| 丰满乱子伦码专区| av天堂中文字幕网| 亚洲人成伊人成综合网2020| 国产高清激情床上av| 国产av麻豆久久久久久久| 日日干狠狠操夜夜爽| 三级国产精品欧美在线观看| 成年免费大片在线观看| 97超视频在线观看视频| 国产亚洲精品综合一区在线观看| 在线观看午夜福利视频| av中文乱码字幕在线| 国产av一区在线观看免费| 综合色av麻豆| 欧美区成人在线视频| 97超级碰碰碰精品色视频在线观看| 免费av毛片视频| 国产精品不卡视频一区二区 | 国产精品久久久久久亚洲av鲁大| 欧美国产日韩亚洲一区| 一个人免费在线观看的高清视频| 麻豆成人av在线观看| 亚洲欧美激情综合另类| 亚洲av.av天堂| 搡老熟女国产l中国老女人| 午夜免费激情av| 村上凉子中文字幕在线| 成人美女网站在线观看视频| 91在线精品国自产拍蜜月| 两个人的视频大全免费| 亚洲男人的天堂狠狠| 有码 亚洲区| 亚洲欧美日韩东京热| av女优亚洲男人天堂| 国产真实伦视频高清在线观看 | 一进一出抽搐gif免费好疼| 欧美黄色淫秽网站| 亚洲熟妇中文字幕五十中出| 久久久久国产精品人妻aⅴ院| 国产人妻一区二区三区在| 一区福利在线观看| 少妇的逼水好多| 日韩欧美精品免费久久 | 久久久久精品国产欧美久久久| 久久久精品欧美日韩精品| 九九在线视频观看精品| 国产单亲对白刺激| 99在线人妻在线中文字幕| 国产白丝娇喘喷水9色精品| 在线免费观看不下载黄p国产 | 午夜免费激情av| 亚洲国产精品sss在线观看| 免费在线观看日本一区| 欧美黑人巨大hd| 此物有八面人人有两片| 九色国产91popny在线| 亚洲av第一区精品v没综合| 国产高清视频在线播放一区| 色尼玛亚洲综合影院| 国产精品自产拍在线观看55亚洲| 国产精品亚洲一级av第二区| 国产亚洲欧美98| 长腿黑丝高跟| 欧美不卡视频在线免费观看| x7x7x7水蜜桃| 国产精品一及| 少妇熟女aⅴ在线视频| 精品人妻熟女av久视频| 黄色女人牲交| 亚洲,欧美,日韩| 国产精品嫩草影院av在线观看 | 天堂√8在线中文| 国产精品98久久久久久宅男小说| 两人在一起打扑克的视频| 一二三四社区在线视频社区8| 老熟妇乱子伦视频在线观看| 久久久久性生活片| 欧美黄色片欧美黄色片| 国产精品一区二区三区四区免费观看 | 长腿黑丝高跟| 国产精品1区2区在线观看.| 国产亚洲欧美98| 少妇的逼好多水| 国产老妇女一区| 99久国产av精品| 90打野战视频偷拍视频| 亚洲精品在线美女| 免费黄网站久久成人精品 | 精品福利观看| 成年免费大片在线观看| 变态另类丝袜制服| 欧美区成人在线视频| 色综合欧美亚洲国产小说| 一二三四社区在线视频社区8| .国产精品久久| 色视频www国产| 成年女人永久免费观看视频| 日韩中文字幕欧美一区二区| 校园春色视频在线观看| av在线天堂中文字幕|