• <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ò)模型
    精品国产乱码久久久久久男人| 免费av中文字幕在线| 最近中文字幕2019免费版| h视频一区二区三区| 国产黄色视频一区二区在线观看| 美女xxoo啪啪120秒动态图| 日韩精品有码人妻一区| 亚洲欧美一区二区三区黑人 | 成年女人毛片免费观看观看9 | 高清黄色对白视频在线免费看| 免费在线观看完整版高清| 十八禁高潮呻吟视频| 曰老女人黄片| 久久久久人妻精品一区果冻| 欧美精品一区二区免费开放| 大陆偷拍与自拍| 日韩熟女老妇一区二区性免费视频| 青春草视频在线免费观看| 999精品在线视频| 在线观看www视频免费| 国产片特级美女逼逼视频| 色婷婷久久久亚洲欧美| 我要看黄色一级片免费的| 亚洲成人手机| av国产久精品久网站免费入址| 日本wwww免费看| 欧美老熟妇乱子伦牲交| 国产白丝娇喘喷水9色精品| 国产精品成人在线| 秋霞伦理黄片| av在线app专区| 国产精品嫩草影院av在线观看| 男女边吃奶边做爰视频| 看非洲黑人一级黄片| 永久免费av网站大全| 亚洲久久久国产精品| 欧美精品高潮呻吟av久久| 国产乱来视频区| 性色avwww在线观看| 成人漫画全彩无遮挡| 色吧在线观看| 91aial.com中文字幕在线观看| 你懂的网址亚洲精品在线观看| 1024香蕉在线观看| 久久精品国产a三级三级三级| 一级a爱视频在线免费观看| 91精品伊人久久大香线蕉| 亚洲av福利一区| 欧美人与性动交α欧美精品济南到 | av在线播放精品| 啦啦啦在线观看免费高清www| 日本色播在线视频| 亚洲欧美一区二区三区黑人 | 国产高清国产精品国产三级| 99热网站在线观看| 亚洲精品,欧美精品| 日本免费在线观看一区| 国产又色又爽无遮挡免| 极品少妇高潮喷水抽搐| 韩国高清视频一区二区三区| 最近手机中文字幕大全| 久久久久久久精品精品| 国产精品免费大片| 亚洲欧美成人综合另类久久久| 亚洲精品,欧美精品| 高清黄色对白视频在线免费看| 伦理电影大哥的女人| 亚洲视频免费观看视频| 精品亚洲成国产av| 久久久精品区二区三区| 国产成人精品福利久久| 超碰成人久久| 久久97久久精品| 一级毛片电影观看| 99久久中文字幕三级久久日本| 国产一区二区三区av在线| 男女下面插进去视频免费观看| 成人亚洲精品一区在线观看| 亚洲精华国产精华液的使用体验| 91成人精品电影| 中文字幕亚洲精品专区| 欧美成人午夜精品| 99久久精品国产国产毛片| 一二三四中文在线观看免费高清| 97在线视频观看| 亚洲精品乱久久久久久| 国产极品粉嫩免费观看在线| 1024视频免费在线观看| 在线天堂最新版资源| 午夜福利,免费看| 日韩一区二区视频免费看| 精品亚洲成国产av| 成年av动漫网址| 国产精品久久久久久精品古装| 97人妻天天添夜夜摸| 亚洲精品国产av成人精品| 啦啦啦在线免费观看视频4| 激情五月婷婷亚洲| 国产一区亚洲一区在线观看| 精品一区二区三区四区五区乱码 | 亚洲色图综合在线观看| 精品一区在线观看国产| 中文字幕色久视频| 久久精品久久久久久噜噜老黄| 免费人妻精品一区二区三区视频| av网站免费在线观看视频| 黄频高清免费视频| 欧美日韩综合久久久久久| 黄色配什么色好看| 亚洲第一青青草原| 少妇被粗大的猛进出69影院| 国产又爽黄色视频| 国产野战对白在线观看| 天天躁夜夜躁狠狠躁躁| 婷婷色av中文字幕| 久久精品久久久久久久性| 麻豆乱淫一区二区| 99久久精品国产国产毛片| 国产免费视频播放在线视频| 午夜免费鲁丝| 国产亚洲午夜精品一区二区久久| av片东京热男人的天堂| 日韩精品有码人妻一区| 成人黄色视频免费在线看| 亚洲精品久久久久久婷婷小说| 天天躁夜夜躁狠狠久久av| 最近最新中文字幕大全免费视频 | 日韩免费高清中文字幕av| 91精品伊人久久大香线蕉| 黄色毛片三级朝国网站| 国精品久久久久久国模美| 91久久精品国产一区二区三区| 2018国产大陆天天弄谢| av又黄又爽大尺度在线免费看| 日韩大片免费观看网站| 91在线精品国自产拍蜜月| 校园人妻丝袜中文字幕| 黄片播放在线免费| 国产一区二区 视频在线| 青春草亚洲视频在线观看| 精品一区在线观看国产| 男女国产视频网站| 18禁裸乳无遮挡动漫免费视频| 性色avwww在线观看| 免费看av在线观看网站| 在现免费观看毛片| 亚洲图色成人| 成人亚洲精品一区在线观看| 26uuu在线亚洲综合色| 国产精品久久久久久av不卡| 亚洲成人一二三区av| 欧美97在线视频| videosex国产| 80岁老熟妇乱子伦牲交| 久久久欧美国产精品| 欧美97在线视频| 亚洲精品久久成人aⅴ小说| 国产不卡av网站在线观看| 十八禁网站网址无遮挡| 日韩av在线免费看完整版不卡| 高清av免费在线| 观看美女的网站| 丝袜脚勾引网站| 久久精品夜色国产| 成人二区视频| 国产精品免费视频内射| 亚洲国产欧美网| 亚洲欧洲精品一区二区精品久久久 | 亚洲美女搞黄在线观看| 中文字幕制服av| 欧美 日韩 精品 国产| 999久久久国产精品视频| 国产精品蜜桃在线观看| 日韩制服丝袜自拍偷拍| 久久精品久久久久久久性| 久久人人爽人人片av| 性色av一级| 蜜桃国产av成人99| 精品国产超薄肉色丝袜足j| 少妇精品久久久久久久| 亚洲国产色片| 欧美xxⅹ黑人| 天天躁夜夜躁狠狠久久av| 国产极品天堂在线| 日本午夜av视频| 看非洲黑人一级黄片| 亚洲av男天堂| 国产97色在线日韩免费| 久久久国产一区二区| av一本久久久久| 国产精品免费视频内射| 久久av网站| 国产精品国产av在线观看| 自线自在国产av| 亚洲精品在线美女| 亚洲精品在线美女| 久久人人爽人人片av| 国产精品免费视频内射| 欧美精品亚洲一区二区| 一级a爱视频在线免费观看| 亚洲国产精品成人久久小说| 在线观看免费高清a一片| 国产日韩一区二区三区精品不卡| 国产免费一区二区三区四区乱码| 久久久久精品性色| 一边摸一边做爽爽视频免费| 另类精品久久| 精品亚洲乱码少妇综合久久| 国产日韩欧美在线精品| 精品视频人人做人人爽| 亚洲第一青青草原| 国产一区二区激情短视频 | 十八禁网站网址无遮挡| 色94色欧美一区二区| 免费女性裸体啪啪无遮挡网站| 天天操日日干夜夜撸| 亚洲av中文av极速乱| 国产精品一国产av| 亚洲综合色惰| 亚洲欧洲日产国产| 9色porny在线观看| av卡一久久| www日本在线高清视频| 麻豆乱淫一区二区| 伊人久久大香线蕉亚洲五| 在线观看人妻少妇| 国产av一区二区精品久久| 国产在线视频一区二区| 美女脱内裤让男人舔精品视频| 纵有疾风起免费观看全集完整版| 欧美亚洲 丝袜 人妻 在线| 日本午夜av视频| 十分钟在线观看高清视频www| 夫妻午夜视频| av国产精品久久久久影院| 亚洲欧美清纯卡通| 久久久久国产精品人妻一区二区| 999精品在线视频| 免费黄频网站在线观看国产| 亚洲精品久久成人aⅴ小说| 尾随美女入室| 免费高清在线观看日韩| 伊人亚洲综合成人网| 波多野结衣一区麻豆| av免费在线看不卡| 伦理电影免费视频| 青春草亚洲视频在线观看| 国产成人a∨麻豆精品| 国产精品嫩草影院av在线观看| 国产亚洲一区二区精品| 亚洲图色成人| 国产av一区二区精品久久| 黄色一级大片看看| 成人毛片60女人毛片免费| 欧美精品一区二区免费开放| 十分钟在线观看高清视频www| 夫妻午夜视频| 午夜福利视频精品| 97人妻天天添夜夜摸| 午夜精品国产一区二区电影| 欧美少妇被猛烈插入视频| 黄色怎么调成土黄色| 99久久人妻综合| 视频在线观看一区二区三区| 啦啦啦中文免费视频观看日本| 国产xxxxx性猛交| 狠狠婷婷综合久久久久久88av| 久久久久久久久久久久大奶| 欧美日韩亚洲国产一区二区在线观看 | 色视频在线一区二区三区| 久久精品国产综合久久久| 好男人视频免费观看在线| 久久97久久精品| 黄色视频在线播放观看不卡| 丝袜美足系列| 精品少妇久久久久久888优播| 国产成人a∨麻豆精品| 在现免费观看毛片| 国产一级毛片在线| 国产成人一区二区在线| 黄频高清免费视频| 人人妻人人澡人人爽人人夜夜| 国产精品无大码| 精品国产一区二区久久| 人人妻人人澡人人看| 在线观看一区二区三区激情| 亚洲经典国产精华液单| 两性夫妻黄色片| 久久99一区二区三区| 赤兔流量卡办理| 你懂的网址亚洲精品在线观看| 久久99一区二区三区| 国产一区二区在线观看av| 精品国产一区二区久久| 美女国产视频在线观看| 在线观看美女被高潮喷水网站| 久热久热在线精品观看| 女的被弄到高潮叫床怎么办| 永久网站在线| 亚洲成色77777| 国产国语露脸激情在线看| 精品一区二区三区四区五区乱码 | 国产男女超爽视频在线观看| 久久久久久伊人网av| 亚洲精品视频女| 久久精品久久久久久噜噜老黄| 午夜激情久久久久久久| 中文字幕亚洲精品专区| 大话2 男鬼变身卡| 两个人看的免费小视频| 黑人巨大精品欧美一区二区蜜桃| 尾随美女入室| freevideosex欧美| 亚洲成国产人片在线观看| 国产成人免费无遮挡视频| 亚洲精品中文字幕在线视频| 欧美变态另类bdsm刘玥| 精品国产露脸久久av麻豆| 99香蕉大伊视频| 色哟哟·www| 又黄又粗又硬又大视频| 国产成人欧美| 考比视频在线观看| 寂寞人妻少妇视频99o| 久久久精品国产亚洲av高清涩受| 国产一区二区 视频在线| 曰老女人黄片| 涩涩av久久男人的天堂| 久久鲁丝午夜福利片| 18在线观看网站| 大片电影免费在线观看免费| 国产精品熟女久久久久浪| 老熟女久久久| 亚洲国产日韩一区二区| 亚洲av国产av综合av卡| 国产麻豆69| 久久97久久精品| 在线亚洲精品国产二区图片欧美| 精品少妇黑人巨大在线播放| 青草久久国产| 国产精品三级大全| 亚洲av免费高清在线观看| 十八禁高潮呻吟视频| 伦理电影大哥的女人| 最黄视频免费看| 热re99久久国产66热| 亚洲美女搞黄在线观看| 日韩,欧美,国产一区二区三区| 亚洲精品,欧美精品| xxx大片免费视频| 夜夜骑夜夜射夜夜干| 亚洲国产精品国产精品| 性高湖久久久久久久久免费观看| 一边亲一边摸免费视频| 五月天丁香电影| 美女午夜性视频免费| 最近最新中文字幕大全免费视频 | 久久精品国产亚洲av天美| videosex国产| 久久久国产一区二区| 亚洲欧美中文字幕日韩二区| 丝袜喷水一区| 丝袜美腿诱惑在线| 亚洲国产av影院在线观看| 婷婷色综合www| av视频免费观看在线观看| 日韩人妻精品一区2区三区| 高清av免费在线| 9热在线视频观看99| 欧美变态另类bdsm刘玥| 国产一区二区三区av在线| 久久青草综合色| 桃花免费在线播放| 在现免费观看毛片| 日本91视频免费播放| 国产av码专区亚洲av| 亚洲av福利一区| 日韩一区二区三区影片| 热99久久久久精品小说推荐| 视频在线观看一区二区三区| 亚洲av电影在线观看一区二区三区| 亚洲成人一二三区av| 久久久久久久久久久免费av| 亚洲欧美成人综合另类久久久| 91精品伊人久久大香线蕉| 高清在线视频一区二区三区| 中文天堂在线官网| 黄色毛片三级朝国网站| 午夜福利一区二区在线看| 国产精品亚洲av一区麻豆 | 好男人视频免费观看在线| 女人被躁到高潮嗷嗷叫费观| 天堂中文最新版在线下载| 国产高清国产精品国产三级| 天堂俺去俺来也www色官网| 欧美变态另类bdsm刘玥| 1024视频免费在线观看| 国产免费一区二区三区四区乱码| 国产乱人偷精品视频| 精品国产乱码久久久久久男人| 中国三级夫妇交换| 亚洲五月色婷婷综合| 国产精品二区激情视频| 69精品国产乱码久久久| 精品国产国语对白av| 亚洲欧美色中文字幕在线| 91午夜精品亚洲一区二区三区| 最近最新中文字幕免费大全7| 2021少妇久久久久久久久久久| 国产色婷婷99| 美女主播在线视频| 久久午夜福利片| 欧美av亚洲av综合av国产av | 成人免费观看视频高清| 久久久久久伊人网av| 欧美 亚洲 国产 日韩一| 熟女av电影| 国产熟女午夜一区二区三区| 美女福利国产在线| 午夜福利一区二区在线看| 亚洲精品,欧美精品| 国产av国产精品国产| 黄片小视频在线播放| 高清av免费在线| 亚洲国产精品国产精品| 一区二区三区乱码不卡18| 少妇精品久久久久久久| 精品国产一区二区久久| 一区二区三区激情视频| 十八禁网站网址无遮挡| 一边亲一边摸免费视频| 丝袜喷水一区| 久久久精品国产亚洲av高清涩受| 久久久久久久大尺度免费视频| 亚洲熟女精品中文字幕| 久久精品久久久久久久性| 久久99精品国语久久久| 亚洲欧美日韩另类电影网站| 99国产综合亚洲精品| 久久久久久久久免费视频了| 亚洲精品一二三| 国产男女内射视频| 欧美在线黄色| 女人精品久久久久毛片| 老熟女久久久| 极品人妻少妇av视频| 伊人久久大香线蕉亚洲五| 日本色播在线视频| 青春草国产在线视频| 男的添女的下面高潮视频| 久久国内精品自在自线图片| 精品国产露脸久久av麻豆| 久久国产亚洲av麻豆专区| 国产精品二区激情视频| 十分钟在线观看高清视频www| 亚洲国产精品国产精品| 黑人猛操日本美女一级片| 一级a爱视频在线免费观看| 亚洲久久久国产精品| 七月丁香在线播放| 啦啦啦啦在线视频资源| 制服丝袜香蕉在线| 天堂8中文在线网| 国产色婷婷99| 日韩中字成人| 五月天丁香电影| 99精国产麻豆久久婷婷| 日日摸夜夜添夜夜爱| 亚洲,欧美,日韩| 最新中文字幕久久久久| 妹子高潮喷水视频| 免费观看在线日韩| xxxhd国产人妻xxx| 男人添女人高潮全过程视频| 精品久久久久久电影网| 日本欧美国产在线视频| 国产淫语在线视频| 两性夫妻黄色片| 一本久久精品| 亚洲,欧美精品.| 国产精品欧美亚洲77777| av线在线观看网站| 亚洲精品久久久久久婷婷小说| 午夜老司机福利剧场| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 满18在线观看网站| 日日摸夜夜添夜夜爱| 欧美变态另类bdsm刘玥| 日日撸夜夜添| 热re99久久精品国产66热6| 亚洲欧美一区二区三区黑人 | 考比视频在线观看| av免费在线看不卡| 黄色 视频免费看| 中文字幕人妻丝袜制服| 啦啦啦中文免费视频观看日本| 久久这里有精品视频免费| 久久久久久伊人网av| 婷婷色av中文字幕| 国产精品 国内视频| 一区二区三区乱码不卡18| 老汉色av国产亚洲站长工具| 国产精品欧美亚洲77777| 国产乱人偷精品视频| 中文欧美无线码| 97人妻天天添夜夜摸| 熟女电影av网| 一区二区三区乱码不卡18| 伦精品一区二区三区| 中国国产av一级| 国产av一区二区精品久久| 777米奇影视久久| 亚洲欧美色中文字幕在线| 成人国语在线视频| 丁香六月天网| 久久久久久伊人网av| 女人高潮潮喷娇喘18禁视频| 少妇人妻 视频| 黑人巨大精品欧美一区二区蜜桃| 91aial.com中文字幕在线观看| 麻豆精品久久久久久蜜桃| 男女边摸边吃奶| 人妻少妇偷人精品九色| 亚洲av日韩在线播放| av不卡在线播放| 你懂的网址亚洲精品在线观看| 好男人视频免费观看在线| 成人免费观看视频高清| 国产1区2区3区精品| 美女高潮到喷水免费观看| 久久综合国产亚洲精品| 欧美变态另类bdsm刘玥| √禁漫天堂资源中文www| 80岁老熟妇乱子伦牲交| 又粗又硬又长又爽又黄的视频| 国产亚洲精品第一综合不卡| 精品国产一区二区久久| 亚洲伊人色综图| 亚洲综合色惰| 丁香六月天网| 国产一区亚洲一区在线观看| 国产免费又黄又爽又色| xxx大片免费视频| 男女无遮挡免费网站观看| 国产精品三级大全| 赤兔流量卡办理| 亚洲精品国产色婷婷电影| 久热这里只有精品99| 爱豆传媒免费全集在线观看| 丝袜喷水一区| 又粗又硬又长又爽又黄的视频| 午夜影院在线不卡| 欧美日韩亚洲国产一区二区在线观看 | 777久久人妻少妇嫩草av网站| 你懂的网址亚洲精品在线观看| 免费在线观看视频国产中文字幕亚洲 | 国产高清不卡午夜福利| 一区二区三区乱码不卡18| 久久人妻熟女aⅴ| 男人爽女人下面视频在线观看| 精品人妻偷拍中文字幕| 免费黄色在线免费观看| 亚洲人成电影观看| 蜜桃国产av成人99| av福利片在线| 伊人久久国产一区二区| 99久久中文字幕三级久久日本| 国产精品香港三级国产av潘金莲 | 波野结衣二区三区在线| 亚洲色图综合在线观看| 人人妻人人澡人人爽人人夜夜| 精品久久久精品久久久| 日本-黄色视频高清免费观看| 免费日韩欧美在线观看| 最新的欧美精品一区二区| 三上悠亚av全集在线观看| 一区二区av电影网| 免费黄频网站在线观看国产| 有码 亚洲区| 91精品三级在线观看| videossex国产| 亚洲第一区二区三区不卡| 少妇人妻精品综合一区二区| 亚洲精品,欧美精品| 国产成人欧美| 亚洲国产av影院在线观看| 自线自在国产av| av免费在线看不卡| 久久久久久久亚洲中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 亚洲国产最新在线播放| 中文字幕另类日韩欧美亚洲嫩草| 免费在线观看视频国产中文字幕亚洲 | 不卡视频在线观看欧美| av在线播放精品| 亚洲成人av在线免费| 日本欧美国产在线视频| 天天操日日干夜夜撸| 久久人人爽人人片av| 国产熟女午夜一区二区三区| 久久精品国产亚洲av天美| 久久这里只有精品19| 在线观看一区二区三区激情| 国产精品三级大全| 18禁国产床啪视频网站| 观看美女的网站| 亚洲中文av在线| 久久国内精品自在自线图片| 午夜福利网站1000一区二区三区| 欧美精品亚洲一区二区| 中文欧美无线码| 久久久久久免费高清国产稀缺| 国产极品粉嫩免费观看在线| 久久久a久久爽久久v久久| 日韩av免费高清视频|