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

    單孔射流流化床內顆?;旌咸匦缘臄抵的M

    2017-10-14 07:05:36張俊強紀律李斌陳海生朱瑞于文圣
    化工學報 2017年3期
    關鍵詞:床層流化床徑向

    張俊強,紀律,李斌,陳海生,朱瑞,于文圣

    ?

    單孔射流流化床內顆?;旌咸匦缘臄抵的M

    張俊強1,2,紀律2,李斌1,陳海生2,朱瑞1,2,于文圣1

    (1華北電力大學動力工程系,河北保定071003;2中國科學院工程熱物理研究所,北京100190)

    在歐拉-拉格朗日坐標系下,采用離散單元法對單孔射流流化床內顆?;旌咸匦赃M行了數值模擬。引入混合指數對床內軸向及徑向布置的顆?;旌腺|量進行定量分析,并研究了不同表觀氣速、不同彈性系數對顆?;旌咸匦缘挠绊?。模擬得到了顆粒軸向及徑向混合序列圖、氣體和顆粒速度分布、整床顆?;旌现笖捣植?、參量變化時整床顆?;旌现笖捣植肌=Y果表明:流化床床層內顆?;旌纤俣仁茴w粒內循環(huán)能力和顆粒擴散能力的綜合作用。單口射流噴動流化床顆粒軸向混合速度主要由顆粒內循環(huán)速度決定,顆粒徑向混合速度主要由顆粒擴散能力決定。表觀氣速增大時,顆粒內循環(huán)速度增加,從而加快了顆粒軸向混合進程,但對顆粒徑向混合影響微弱;彈性系數增大時,顆?;旌纤俣燃盎旌腺|量均下降,并且彈性系數增大對顆粒徑向混合進程影響小于顆粒軸向混合。

    流化床;顆粒;混合;混合指數;定量評價;數值模擬

    引 言

    氣固流化床以其燃料適應性廣、環(huán)保性能優(yōu)越、制備系統(tǒng)簡單等諸多優(yōu)點被廣泛應用于各工業(yè)領域。氣固流化床內的顆?;旌鲜欠从彻腆w顆粒在床內運動及傳熱、傳質過程的重要行為,對流化床反應器的研究、設計和操作具有重要作用[1]。流化床內稠密氣固兩相流體系中氣體與顆粒的作用機理復雜,直接影響了床內的動量、熱量和質量交換以及化學反應速率與反應效率,因此了解氣固兩相的動態(tài)混合行為對流化床反應器的設計與操作至關重要。

    流化床作為一種反應器,用實驗手段對噴動流化床氣固流動特性的研究較多,但由于受到實驗條件以及設備的限制,常規(guī)的技術難以捕捉到微觀層次上的顆粒運動信息,如顆粒速度、顆?;旌腺|量和基于顆粒尺度的平均碰撞頻率等[2-3],這是一個極為復雜的過程組合,其內在規(guī)律尚且不被人們完全了解。

    隨著計算機技術的發(fā)展,數值模擬在流化床內氣固兩相流的研究中發(fā)揮著越來越重要的作用。在眾多模擬計算方法中,離散單元法(discrete element method, DEM)以其獲取微觀信息豐富、模擬精確度高等優(yōu)點,被學者廣泛應用到氣固兩相流研究中[4-8]。目前,針對流化床內顆?;旌咸匦缘难芯渴橇骰驳难芯繜狳c之一[9-20],任立波等[14]對脈沖鼓泡床內鼓泡和顆?;旌咸匦赃M行了CFD-DEM數值模擬。Luo等[15]采用CFD-DEM方法對鼓泡流化床內固相的混合和分離行為進行了模擬研究。Olaofe等[16]采用CFD-DEM方法對流化床內顆粒的混合和偏析過程進行了數值模擬。Oke等[17]分別采用歐拉-歐拉和歐拉-拉格朗日模型對軸向布置顆粒的混合特性進行了研究。Renzo等[19]運用離散單元法研究了流化床內兩種不同密度顆粒的混合平衡過程。但通過數值模擬與實驗對比,且從微觀層次進行較為全面的定量研究還較少[21-27]。

    本文采用離散單元法對單孔射流流化床內顆?;旌咸匦赃M行了數值模擬,并與所搭實驗臺得到結果進行對比驗證。通過引入混合指數對床內軸向及徑向布置顆粒的混合質量進行定量評價,以歸納單孔射流流化床內主導顆粒軸向混合、徑向混合進程的作用機制,并從微觀層次揭示流化床內顆粒混合的機理。

    1 數學模型

    1.1 顆粒運動模型

    顆粒的運動主要受顆粒與顆粒(壁面)之間的碰撞力、氣體的曳力和自身重力作用[28-30],顆粒的運動控制方程為

    式中,F,j為流體對顆粒的曳力;p為顆粒間的碰撞力。

    對顆粒的模擬采用離散單元法軟球模型,顆粒與顆粒間的碰撞力包括彈性力和黏性阻尼力,法向力n,ij分解為法向彈性力cn,ij和法向黏性阻尼力dn,ij,切向力t,ij分解為彈性切向力ct,ij和切向黏性阻尼力dt,ij。同時顆粒還收到自身重力mg和流體對其的曳力f,i,由于本文在數值模擬時采用準三維模型,并選用直徑相同的球形顆粒,顆粒轉動對模擬結果的影響不太明顯,所以暫不考慮旋轉。具體受力分析如圖1所示。

    根據胡克定律可得出顆粒的碰撞力計算式。

    顆粒所受法向力

    cn,ij()=cn,ij(-D)-nDdn(2)

    n,ij=cn,ij+dn,ij(4)

    顆粒所受切向力

    ct,ij()=ct,ij(-D)-tDdt(5)

    t,ij=ct,ij+dt,ij(7)

    式中,為摩擦系數;n為法向彈性系數;t為切向彈性系數;n為法向黏性阻尼系數;t為切向黏性阻尼系數;n為兩顆粒間的法向相對位移;t為兩顆粒間的切向相對位移。

    1.2 氣相運動模型

    氣相模型采用考慮到氣固兩相耦合作用的Navier-Stocks方程,具體公式如下。

    連續(xù)性方程

    動量方程

    (10)

    式中,g為流體剪切黏度;t為湍流黏度,在旺盛湍流區(qū),剪切黏度可忽略不計;g為孔隙率;p,x、p,y為顆粒對流體的反作用力。

    氣相湍流模擬采用兩方程模型描述,其控制方程如下

    (13)

    其中,-方程各項構成及系數取值見文獻[4]。

    1.3 氣固兩相之間的耦合作用

    1.3.1 氣相對顆粒的曳力 采用Gidaspow的曳力模型[11],氣體對顆粒的曳力計算式如下。

    g>0.8時,利用Wen&Yu等方程得

    g≤0.8時,采用Ergun等式得

    (15)

    式中,g為氣體速度;p為顆粒速度;g為氣體剪切黏度;g為空隙率;g為氣體密度;p為顆粒體積;p為顆粒直徑;′d為有效曳力系數。

    (16)

    式中,為空隙率修正因子,一般的取值為4.65;d為單顆粒的曳力系數。

    當p≤1000時

    當p>1000時

    d=0.44 (18)

    顆粒Reynolds數的計算公式為

    空隙率的計算式為

    (20)

    式中,D為劃分區(qū)域的網格體積;p為顆粒體積;c為所劃分網格內顆粒的個數。

    1.3.2 顆粒對氣相的反作用力

    式中,F,i為氣體對顆粒的作用力。

    2 模擬對象及參數

    模擬對象為150 mm×4 mm×900 mm的單孔射流準三維矩形噴動床,床體底部中間位置設有一氣流噴口,氣流噴口寬度為10 mm,模擬對象尺寸如圖2所示。顆粒相采用2400個直徑為4 mm的球形顆粒,顆粒密度為2700 kg·m-3。氣相采用密度為1.205 kg·m-3的空氣,氣體黏度為1.8×10-5 kg·(m·s)-1,同時本模擬已進行網格無關性驗證,具體計算參數見表1。

    表1 數值模擬參數

    床體上部采用局部單向化處理,壁面處采用無滑移邊界條件。將顆粒均勻分為黑色和紅色2組,初始狀態(tài)采用軸向、徑向均分自然堆積。本節(jié)研究了單孔射流流化床內氣固兩相的流動特性及混合特性,從氣固兩相運動的微觀層面揭示顆粒的混合機制,并分析了參量變化對床內顆?;旌线M程的影響。

    3 模擬結果

    3.1 顆?;旌闲蛄袌D

    靜止床層時顆粒沿軸向和徑向均分成兩組,并標記為黑和紅兩種顏色。圖3和圖4為噴動氣速為30 m·s-1,模擬時間為2 s內床內顆粒軸向及徑向混合序列圖。由圖3顆粒軸向混合序列圖可得,噴動氣體進入后,噴射區(qū)底部的顆粒在氣流的曳力作用下向上運動,顆粒上行過程中與周圍的顆粒不斷相互碰撞,當向上運動的顆粒到達床層表面后,在噴泉區(qū)向兩側壁面進行擴散。壁面處氣速較小,擴散后的顆粒自身重力大于其所受曳力后開始沿兩側壁面下行,參與環(huán)隙區(qū)內顆?;旌线\動,最后再次進入噴射區(qū),形成了完整的顆粒內循環(huán)。循環(huán)周而復始,最終達到顆?;旌系膭討B(tài)平衡。同時由圖4顆粒徑向混合序列圖可知,顆粒的強烈混合首先發(fā)生于床層表面,位于床層表面的顆粒不斷在下部顆粒的作用下向兩側擴散,促進著顆粒的混合進程。綜上可知,床層顆?;旌纤俣仁茴w粒內循環(huán)能力和顆粒擴散能力的綜合作用,這兩項能力的改善將有效促進流化床內顆粒的混合進程。

    圖4 顆粒徑向混合序列圖

    3.2 顆粒及氣體速度分布

    3.2.1 顆粒速度分布 圖5為0.5~2 s內不同床層高度下顆粒時均速度分布。由圖5(a)顆粒水平方向速度分布可知,床層底部中心區(qū)域左側顆粒速度值為正,顆粒向右運動,底部中心區(qū)域右側顆粒速度值為負,顆粒向左運動。在接近床層表面處,中心區(qū)域左側顆粒速度為負,顆粒向左運動;中心區(qū)域右側顆粒速度為正,顆粒向右運動。由圖5(b)顆粒豎直方向速度分布可知,床層中心區(qū)域不同高度處顆粒速度值為正,顆粒向上運動,兩側壁面不同高度處顆粒速度值為負,顆粒向下運動。隨著床高增加,顆粒豎直方向速度呈先增大后減小的變化趨勢。由此分析,隨著中心射流噴入,床層底部噴射區(qū)內的顆粒受氣體卷吸作用向中間匯聚,再隨氣流向上運動,與上層顆粒不斷混合。顆粒到達床層表面噴泉區(qū)后向兩側壁面擴散,此時顆粒所受重力大于空氣對其的曳力作用,顆粒開始向下運動并進入環(huán)隙區(qū),與環(huán)隙區(qū)的顆粒混合后,緩慢向環(huán)隙區(qū)和噴射區(qū)交界面運動,進入下一個循環(huán),如此反復,混合程度不斷加深,最終達到一種動態(tài)平衡。

    3.2.2 氣體速度分布 圖6為0~2 s內不同床層高度下氣體時均速度分布。氣體速度分布規(guī)律與顆粒速度分布較為相似,水平方向上,床層底部氣體由兩層向中間噴口處匯聚,較高床層處氣體向兩側擴散,并隨著床高增加氣體速度值呈先減小后增大再減小的變化趨勢。豎直方向上,床層中心區(qū)域氣體向上運動,近壁面處氣體向下運動,并隨著床高增加氣體速度值亦逐漸減小。結合顆粒速度分布可知,顆粒的速度分布與氣體速度分布趨勢較為相似,說明氣體曳力作用下的顆粒運動是顆粒混合的主要原因。同時由氣體速度分布可以看出,氣體豎直速度數值在相同床層位置處明顯大于氣體水平速度,這也較好地解釋了軸向布置顆粒混合進程快于徑向混合的原因。

    圖7為=0.5 s顆粒速度矢量圖,顆粒矢量圖中顆粒的運動軌跡較明顯地形成了流化床內噴射區(qū)、噴泉區(qū)和環(huán)隙區(qū)結構。該圖得到的顆粒移動進程與顆粒混合序列圖及顆粒速度分布分析得出的顆粒運動進程相符,同時也證實了前文中對于氣固流動特性分析的正確性。

    為了驗證所用模型及數值模擬方法的準確性,本文進行了相關實驗研究,測量了相同條件下床內顆粒流場瞬時圖。圖8為顆粒流動模擬與實驗瞬時圖對比,可以得出模擬與實驗中床層軸向、徑向混合時,床內顆?;旌蠣顟B(tài)均具有較好的一致性。

    3.3 混合程度定量評價

    隨著顆?;旌线\動的不斷進行,單純在直觀上對混合程度分析已難以得出可靠的結論。本文引入Lacey提出的混合指數對顆粒的混合進程來深入研究。選定黑色顆粒為示蹤顆粒,同時將床體平均分成若干個取樣單元,對于有限取樣過程,示蹤顆粒濃度標準差的計算式為

    (23)

    式中,為取樣單元數;為所有取樣單元中示蹤顆粒平均濃度;c為取樣單元內的示蹤顆粒濃度。為更有效地比較不同場合內顆粒的混合進程,引入Lacey混合指數的計算式為

    為了更全面地認識流化床內顆?;旌咸匦约捌渲饕绊懸蛩?,本文將流化床分成5×30個取樣單元,數據結果輸出步長為0.01 s,研究了同一表觀氣速下的整床混合程度的變化進程,以及不同參量變化對顆粒混合程度的影響情況。

    圖9為進口氣速為30 m·s-1時示蹤顆粒軸向、徑向混合指數隨時間變化的過程曲線。由圖可知,軸向和徑向混合工況下顆粒運動起始階段,顆粒混合指數增長十分迅速,但軸向混合指數增速明顯快于徑向混合。后期顆?;旌现笖稻鏊僮兙?,混合指數值在0.9~1之間平穩(wěn)波動,此時可認為顆?;旌馅呌陔S機完全混合狀態(tài),而波動表明混合過程經歷著混合與偏析,由此達到動態(tài)平衡。在此對比工況下,軸向混合指數值在0.9以上規(guī)律波動的起始時間為0.85 s,而徑向混合的起始時間為1.45 s。同時在相同時刻下軸向混合指數均大于徑向。綜上可以得出,在相同氣速條件下,顆粒軸向布置的混合進程明顯快于徑向,顆?;旌纤俣认瓤旌舐詈筮_到顆?;旌系膭討B(tài)平衡。

    3.3.1 進口氣速對混合進程的影響 圖10為不同氣速條件下顆粒軸向、徑向混合指數隨時間的變化。由圖10(a)軸向混合指數分布與圖10(b)徑向混合指數分布對比可得,噴動氣速增加后顆粒軸向混合指數增速有所加快,但是顆粒徑向混合指數增速變化微弱。噴動氣速的提高可有效強化顆粒的軸向混合,但對徑向混合強化不明顯。分析原因可能為,顆粒內循環(huán)機制對上下布置的軸向均分顆粒混合進程影響較大,而對左右布置的徑向均分顆粒影響較小。噴動氣速提高引起顆粒內循環(huán)速度加快,有助于促進床內顆粒混合,同時噴動氣速提高時,軸向混合強于徑向混合。這一結論與文獻[18]研究結果一致。

    3.3.2 彈性系數對混合進程的影響 圖11為不同彈性系數下顆粒軸向、徑向混合指數隨時間的變化。由圖11(a)顆粒軸向混合指數分布可知,隨著顆粒彈性系數不斷變大,導致顆粒軸向混合指數增速變緩,混合指數值達到隨機完全混合范圍的時間變長,顆?;旌线M程隨之延長。而由圖11(b)顆粒徑向混合指數分布可得,雖然顆粒彈性指數增大亦導致徑向混合速度變慢,但其所受影響明顯弱于軸向混合,同時相同彈性系數下顆粒混合質量優(yōu)于軸向混合。分析原因可能為,當顆粒彈性系數較大時,位于床層底部的顆粒會形成一定厚度的“互鎖”床層,且互鎖床層厚度隨彈性系數增大而增大,導致射流很難擊穿床層,互鎖床層在射流的左右有一定的整體上移趨勢,這影響了噴射區(qū)和環(huán)隙區(qū)內顆粒的混合進程,從而造成顆粒軸向、徑向混合速度均隨彈性系數增大而減慢。

    當顆粒彈性系數增大時,床層底部顆粒間存在“互鎖”作用更明顯,床層底部的不動區(qū)厚度增加,顆粒軸向混合的進程變慢。氣流通過顆粒間空隙向上運動,上部顆粒由于壓力較小,因此在噴泉區(qū)內顆?;旌线\動較為劇烈,顆粒徑向混合進程較快,這很好地解釋了彈性系數增大對顆粒徑向混合速度影響小于顆粒軸向混合的原因。同時顆粒整體混合質量會隨彈性系數的增大而下降。

    4 結 論

    (1)噴動流化床床層內顆粒混合速度受顆粒內循環(huán)能力和顆粒擴散能力的綜合作用。

    (2)單口射流噴動流化床顆粒軸向混合速度主要由顆粒內循環(huán)速度決定,內循環(huán)速度增加,軸向混合速度加快。而顆粒徑向混合速度主要由顆粒擴散能力決定,床層擴散能力提高,徑向混合速度加快。

    (3)在單孔射流流化床內,相同氣速條件下顆粒軸向混合速度快于徑向混合。增大表觀氣速提高了顆粒內循環(huán)速度,從而加快了顆粒軸向混合進程,但對顆粒徑向混合影響微弱。

    (4)在單孔射流流化床內,顆粒彈性系數增大時,床層底部顆粒間存在“互鎖”作用明顯,顆粒混合質量會隨彈性系數的增大而下降,同時彈性系數增大對顆粒徑向混合速度影響小于顆粒軸向混合。

    (5)提高進口噴動氣速或顆粒彈性系數,不僅影響著顆粒混合速度,還會造成全床顆?;旌腺|量下降。

    References

    [1] 周月桂, 彭軍, 朱賢, 等. 循環(huán)流化床煙氣脫硫氣固兩相流動特性的試驗研究[J]. 動力工程學報, 2009, 29(6): 560-564. ZHOU Y G, PENG J, ZHU X,. Experimental research on gas-solid two-phase flow characteristic for flue gas desulfurization in circulating fluidized bed[J]. Journal of Power Engineering, 2009, 29(6): 560-564.

    [2] WANG X J, JIN B S, WAN Y Y,. Three-dimensional multi-phase simulation of the mixing and segregation of binary particle mixtures in a two-jet spout fluidized bed[J]. Particuology, 2015, (22): 185-193.

    [3] 張瑞卿, 楊海瑞, 呂俊復. 應用于循環(huán)流化床鍋爐氣固流動和燃燒的CPFD數值模擬[J]. 中國電機工程學報, 2013, 33(23): 75-83. ZHANG R Q, YANG H R, Lü J F. Application of CPFD approach on gas-solid flow and combustion in industrial CFB boilers[J]. Proceedings of the CSEE, 2013, 33(23): 75-83.

    [4] 紀律. 循環(huán)流化床流動及磨損特性的 DEM 數值模擬[D]. 保定: 華北電力大學, 2011. JI L. Numerical simulation of hydrodynamics and abrasion characteristic in circulating fluidized bed by discrete element method[D]. Baoding: North China Electric Power University, 2011.

    [5] HERNáNDEZ-JIMéNEZ F, LI T, CANO-PLEITE E,. Characterization of the particle–wall frictional forces in pseudo-2D fluidized bedsusing DEM[J]. Powder Technology, 2014, 116: 136-143.

    [6] SAIDI M, TABRIZI H B, GRACE J R,. Hydrodynamic investigation of gas-solid flow in rectangular spout-fluid bed using CFD-DEM modeling[J]. Powder Technology, 2015, 284: 355-364.

    [7] ALOBAID F, BARAKI N, EPPLE B. Investigation into improving the efficiency and accuracy of CFD/DEM simulations[J]. Particuology, 2014, (16): 41-53.

    [8] 張勇, 金保昇, 鐘文琪, 等. 噴動流化床顆?;旌咸匦缘娜S直接數值模擬[J]. 中國電機工程學報, 2008, 28(2): 33-38. ZHANG Y, JIN B S, ZHONG W Q,. Three-dimensional DEM simulation on particle mixing characteristics of spout-fluid bed[J]. Proceedings of the CSEE, 2008, 28(2): 33-38.

    [9] 李斌, 紀律. 流化床爐內顆?;旌系碾x散單元法數值模擬[J]. 中國電機工程學報, 2012, 32(20): 42-48. LI B, JI L. Numerical simulation of particle mixing in circulating fluidized bed with discrete element method[J]. Proceedings of the CSEE, 2012, 32(20): 42-48.

    [10] 朱潤孺, 朱衛(wèi)兵, 邢力超, 等. 矩形噴動床混合特性的三維數值研究[J]. 中國電機工程學報, 2010, 30(17): 12-16. ZHU R R, ZHU W B, XING L C,.A three-dimensional numerical investigation on particle mixing characteristics in rectangular spouted beds[J].Proceedings of the CSEE, 2010, 30(17): 12-16.

    [11] 李斌, 宋小龍. 循環(huán)流化床內顆?;旌咸匦缘臄抵的M[J]. 動力工程學報, 2013, 33(10): 759-764. Li B, SONG X L. Numerical simulation on mixing characteristics of particles in circulating fluidized bed[J]. Journal of Chinese Society of Power Engineering, 2013, 33(10):759-764.

    [12] 劉偉偉, 盧春喜, 范怡平, 等. 氣固流化床中雙組分混合顆粒的流態(tài)化特性[J]. 化工學報, 2008, 59(8): 1971-1978. LIU W W, LU C X, FAN Y P,. Flow behavior of binary mixture particles in gas-solid fluidized beds[J]. Journal of Chemical Industry and Engineering (China), 2008, 59(8): 1971-1978.

    [13] 楊旭, 劉雅寧, 余劍, 等. 微型流化床內混合特性的數值模擬[J]. 化工學報, 2014, 65(9): 3323-3330. YANG X, LIU Y N, YU J,. Numerical simulation of mixing characteristics of trace sample and bed material in micro fluidized bed reaction analyzer[J]. CIESC Journal, 2014, 65(9): 3323-3330.

    [14] 任立波, 尚立寶, 閆日雄, 等. 脈沖鼓泡床內鼓泡和顆?;旌咸匦缘腃FD-DEM數值模擬[J]. 山東大學學報, 2015, 45(2): 62-66. REN L B, SHANG L B, YAN R X,. CFD-DEM simulation of bubbling and particle mixing properties in pulsed jet fluidized bed[J]. Journal of Shandong University, 2015, 45(2): 62-66.

    [15] LUO K, WU F, YANG S L,. CFD-DEM study of mixing and dispersion behaviors of solid phase in a bubbling fluidized bed[J]. Powder Technology, 2015, 274: 482-493.

    [16] OLAOFE O, PATIL A, DEEN N, VAN DER HOEF M, KUIPERS J. Simulation of particle mixing and segregation in bidisperse gas fluidized beds[J]. Chemical Engineering Science, 2014, 108: 258-269.

    [17] REN B, SHAO Y, ZHONG W,. Investigation of mixing behaviors in a spouted bed with different density particles using discrete element method[J]. Powder Technology, 2012, 22:85-94.

    [18] 張勇, 金保昇, 鐘文琪. 噴動氣固流化床顆粒混合規(guī)律的實驗研究[J]. 中國電機工程學報, 2008, 28(20): 8-14. ZHANG Y, JIN B S, ZHONG W Q.Experimental investigation on particle mixing in spout-fluid bed[J]. Proceedings of the CSEE, 2008, 28(20): 8-14.

    [19] RENZO A D, MAIO F P D, GIRIMONTE R,. DEM simulation of the mixing equilibrium in fluidized beds of two solids differing in density[J]. Powder Technology, 2008, 184: 214-223.

    [20] REN B, SHAO Y, ZHONG W,. Investigation of mixing behaviors in a spouted bed with different density particles using discrete element method[J]. Powder Technology, 2012, 22:85-94.

    [21] 張科. 復雜稠密氣固兩相流動的CFD-DEM模擬研究[D]. 杭州: 浙江大學, 2012. ZHANG K, CFD-DEM simulation of complex dense gas-solid flow[D]. Hangzhou: Zhejiang University, 2012.

    [22] 彭麗, 吳迎亞, 李佳瑤, 等. 基于DEM模擬氣固鼓泡床中顆粒碰撞參數對流場間歇性的影響[J]. 化工學報, 2015, 66(6): 2041-2048. PENG L, WU Y Y, LI J Y,. Effect of granular collision parameters on DEM simulation of flow field intermittency in gas-solids bubbling fluidized bed[J]. CIESC Journal, 2015, 66(6): 2041-2048.

    [23] 任立波, 韓吉田. 基于CFD-DEM耦合并行算法的錐形噴動床內離散顆粒數值模擬[J]. 東南大學學報, 2014, 44(5): 993-998.REN L B, HAN J T. Numerical simulation of discrete parcicles in conical-base spouted bed based on parallel coupled CFD-DEM model[J]. Journal of Southeast University, 2014, 44(5): 993-998.

    [24] FENG Y Q, YU A B. Microdynamic modelling and analysis of the mixing and segregation of binary mixtures of particles in gas fluidization[J]. Chemical Engineering Science, 2007, 62(1): 256-268.

    [25] 陶賀, 鐘文琪, 張勇, 等. 不同顆粒物性對異徑混合非球形顆粒分離特性影響的數值模擬[J]. 動力工程學報, 2015, 35(8): 652-658. TAO H, ZHONG W Q, ZHANG Y,. Effects of particle properties on separation behavior of differently-sized non-spherical particle mixtures[J]. Journal of Chinese Society of Power Engineering, 2015, 35(8): 652-658.

    [26] LI T, ZHANG Y, HERNáNDEZ-JIMéNEZD F. Investigation of particle-wall interaction in a pseudo-2D ?uidized bedusing CFD-DEM simulations[J]. Particuology, 2016, 25: 10-22.

    [27] 卜昌盛, 陳曉平, 劉道銀, 等. 基于顆粒尺度的離散顆粒傳熱模型[J]. 化工學報, 2012, 63(3): 698-704. BU C S, CHEN X P, LIU D Y,. Heat transfer model for particles with discrete element method[J]. CIESC Journal, 2012, 63(3): 698-704.

    [28] MüLLER C R, SCOTT S A, HOLLAND D J,. Validation of a discrete element model using magnetic resonance measurements[J]. Particuology, 2009, 7(4): 297-306.

    [29] SAIDI M, TABRIZI H B, GRACE R,. Hydrodynamic investigation of gas-solid ?ow in rectangular spout-?uid bed using CFD-DEM Modeling[J]. Powder Technology, 2015, 284: 355-364.

    [30] YANG W J, ZHOU Z Y, YU A B. Discrete particle simulation of solid flow in a three-dimensional blast furnace sector model[J]. Chemical Engineering Journal, 2015, 278: 339-352.

    Numerical simulation of particle mixing in single jet fluidized bed

    ZHANG Junqiang1,2, JI Lü2, LI Bin1, CHEN Haisheng2, ZHU Rui1,2, YU Wensheng1

    (1Department of Energy and Power Engineering, North China Electric Power Uninversity, Baoding 071003, Heibei, China;2Institude of Engineering Themophsics, Chinese Academy of Sciences, Beijing 100190, China)

    In the Euler-Lagrange coordinate system, mixing characteristics of particles in a single jet fluidized bed were studied with discrete element method. A mixing index was introduced to quantitatively analyze axial and radial mixing quality in the bed, as well as to investigate influnce of superficial gas velocity and spring constant on mixing charcteristics. The simulation provided sequence diagrams of axial and radial mixing, distributions of gas and particle velocities, and mixing index distributions of particles across the bed at various parameter settings. The results showed that particle mixing process in the bed was determined by capacity of particle circulation and diffusion. Axial mixing rate was mainly controlled by internal circulation speed of particles, whereas radial mixing rate was mainly controlled by particle diffusion capacity. With the increase of superficial gas velocity, the internal circulation speed of particles was enhanced and thus axial mixing process was accelerated but radial mixing was weakly influenced. When spring constant was increased, both mixing speed and quality of particles were decreased as well as radial mixing process was less affected than that axial mixing by spring constant.

    fluidized bed; particle; mixing; mixing index; quantitative evaluation; numerical simulation

    10.11949/j.issn.0438-1157.20161321

    TQ 051.1

    A

    0438—1157(2017)03—0879—10

    國家自然科學基金優(yōu)秀青年基金項目(51522605);國家國際科技合作專項項目(2014DFA60600);中國科學院可再生能源高效利用創(chuàng)新交叉團隊項目。

    2016-09-21收到初稿,2016-12-07收到修改稿。

    聯(lián)系人:李斌。第一作者:張俊強(1991—),男,碩士研究生。

    2016-09-21.

    LI Bin, binli_871@163.com

    supported by the National Natural Science Foundation of China (51522605), the International S&T Cooperation Program of China (2014DFA60600) and the CAS Interdisciplinary Innovation Team.

    猜你喜歡
    床層流化床徑向
    木屑顆粒堆積參數對液體滲流特性的影響規(guī)律
    徑向分層對隨機堆積固定床內氣體流動特性的影響研究?
    淺探徑向連接體的圓周運動
    RN上一類Kirchhoff型方程徑向對稱正解的存在性
    基于PID+前饋的3MN徑向鍛造機控制系統(tǒng)的研究
    重型機械(2020年3期)2020-08-24 08:31:40
    SiCl4冷氫化反應器床層密度的計算
    四川化工(2020年2期)2020-05-20 10:46:42
    一類無窮下級整函數的Julia集的徑向分布
    微纖維- 活性炭雙床層對苯蒸汽吸附動力學研究
    兵工學報(2019年6期)2019-08-06 01:06:46
    流化床丙烷脫氫反應段的模擬及優(yōu)化
    關于循環(huán)流化床鍋爐集控運行研究
    91精品国产九色| 国产黄片视频在线免费观看| 蜜臀久久99精品久久宅男| 亚洲中文av在线| 国产午夜精品一二区理论片| 老熟女久久久| 国产 精品1| 日韩欧美精品免费久久| 亚洲av电影在线观看一区二区三区| 在线观看免费日韩欧美大片 | 插逼视频在线观看| 午夜免费男女啪啪视频观看| 蜜桃久久精品国产亚洲av| 亚洲人与动物交配视频| 女性被躁到高潮视频| 美女福利国产在线 | 欧美成人a在线观看| 狂野欧美激情性xxxx在线观看| 九九久久精品国产亚洲av麻豆| 亚洲美女黄色视频免费看| 亚洲精品aⅴ在线观看| 少妇高潮的动态图| 久久毛片免费看一区二区三区| 在线观看一区二区三区激情| 99热这里只有是精品在线观看| 免费黄色在线免费观看| 亚洲色图综合在线观看| 99九九线精品视频在线观看视频| 精品久久久久久久久av| 亚洲欧洲国产日韩| 免费人妻精品一区二区三区视频| 国产av精品麻豆| 黄色日韩在线| 不卡视频在线观看欧美| 小蜜桃在线观看免费完整版高清| 亚洲国产毛片av蜜桃av| 国产中年淑女户外野战色| 五月天丁香电影| 卡戴珊不雅视频在线播放| 2018国产大陆天天弄谢| 日韩大片免费观看网站| 久久99蜜桃精品久久| 亚洲国产高清在线一区二区三| 国产精品一区二区在线观看99| 亚洲久久久国产精品| a级一级毛片免费在线观看| 精品久久久噜噜| 午夜精品国产一区二区电影| 久久99热6这里只有精品| 欧美高清性xxxxhd video| 人妻少妇偷人精品九色| 亚洲婷婷狠狠爱综合网| 观看美女的网站| 91久久精品国产一区二区三区| 一级毛片 在线播放| 国产极品天堂在线| 91久久精品电影网| 免费观看在线日韩| 色网站视频免费| 国产欧美亚洲国产| 日本黄色日本黄色录像| 高清视频免费观看一区二区| 日韩欧美一区视频在线观看 | 免费看av在线观看网站| 亚洲精品成人av观看孕妇| 免费黄色在线免费观看| 最近2019中文字幕mv第一页| 亚洲欧美清纯卡通| 一区二区三区乱码不卡18| 最后的刺客免费高清国语| 毛片一级片免费看久久久久| 久久青草综合色| 多毛熟女@视频| 欧美国产精品一级二级三级 | xxx大片免费视频| 久久精品久久久久久噜噜老黄| 在线观看一区二区三区| 五月伊人婷婷丁香| 亚洲成人中文字幕在线播放| 国产午夜精品久久久久久一区二区三区| 免费人妻精品一区二区三区视频| 国产精品无大码| av国产免费在线观看| 亚洲天堂av无毛| 国产精品秋霞免费鲁丝片| 国产精品秋霞免费鲁丝片| 九九久久精品国产亚洲av麻豆| xxx大片免费视频| 一边亲一边摸免费视频| 在线观看国产h片| 欧美日韩视频高清一区二区三区二| 欧美一区二区亚洲| av卡一久久| 国产在视频线精品| 女性被躁到高潮视频| 久久久久国产精品人妻一区二区| 水蜜桃什么品种好| 在线 av 中文字幕| 熟女人妻精品中文字幕| 啦啦啦啦在线视频资源| 日韩一区二区三区影片| 欧美成人精品欧美一级黄| 国产精品蜜桃在线观看| av黄色大香蕉| 国产精品麻豆人妻色哟哟久久| 国产在线男女| 亚洲国产精品999| 伊人久久精品亚洲午夜| av天堂中文字幕网| 亚洲精品乱码久久久久久按摩| 日本与韩国留学比较| 草草在线视频免费看| 日韩制服骚丝袜av| 最近最新中文字幕免费大全7| 97超视频在线观看视频| 精品一区在线观看国产| 一级二级三级毛片免费看| 久久精品久久精品一区二区三区| 91久久精品国产一区二区成人| 能在线免费看毛片的网站| 一级毛片我不卡| 亚洲va在线va天堂va国产| 久久久久网色| 国产人妻一区二区三区在| 在线免费观看不下载黄p国产| 国产视频内射| 国产一级毛片在线| 老熟女久久久| 成人一区二区视频在线观看| 人妻一区二区av| 一个人看视频在线观看www免费| 尾随美女入室| 亚洲经典国产精华液单| 美女脱内裤让男人舔精品视频| 成人毛片a级毛片在线播放| 久久99蜜桃精品久久| 久久久久精品久久久久真实原创| 一区二区av电影网| 精品久久久久久久久av| 毛片一级片免费看久久久久| a 毛片基地| 建设人人有责人人尽责人人享有的 | 自拍偷自拍亚洲精品老妇| 日本午夜av视频| 美女国产视频在线观看| 五月玫瑰六月丁香| 噜噜噜噜噜久久久久久91| 亚洲高清免费不卡视频| 亚洲精品中文字幕在线视频 | 国产高清国产精品国产三级 | 十八禁网站网址无遮挡 | 波野结衣二区三区在线| 老司机影院毛片| 大陆偷拍与自拍| 在线观看免费日韩欧美大片 | 精品亚洲成国产av| 国产伦在线观看视频一区| 国产一区二区三区综合在线观看 | 亚洲美女黄色视频免费看| 爱豆传媒免费全集在线观看| 天美传媒精品一区二区| 99久国产av精品国产电影| 亚洲欧美一区二区三区国产| 中文字幕人妻熟人妻熟丝袜美| kizo精华| 伦理电影大哥的女人| 亚洲一区二区三区欧美精品| 国产精品秋霞免费鲁丝片| 国产高清国产精品国产三级 | 九九在线视频观看精品| 久久99精品国语久久久| 亚洲伊人久久精品综合| 视频中文字幕在线观看| 2021少妇久久久久久久久久久| 欧美成人精品欧美一级黄| kizo精华| 肉色欧美久久久久久久蜜桃| 国产精品福利在线免费观看| 精品午夜福利在线看| 天天躁日日操中文字幕| 秋霞在线观看毛片| 成人综合一区亚洲| 国产在线免费精品| 我的老师免费观看完整版| 亚洲美女视频黄频| 亚洲av电影在线观看一区二区三区| 99久国产av精品国产电影| 久久精品国产亚洲av天美| 亚洲av中文字字幕乱码综合| h日本视频在线播放| 亚洲精品日韩在线中文字幕| 精品午夜福利在线看| 毛片一级片免费看久久久久| 狂野欧美白嫩少妇大欣赏| 高清毛片免费看| 国产精品人妻久久久久久| 精品熟女少妇av免费看| 国产黄频视频在线观看| 国产精品国产三级国产av玫瑰| 91久久精品国产一区二区三区| 国产精品久久久久成人av| av线在线观看网站| 日韩一区二区视频免费看| av福利片在线观看| 久久精品国产自在天天线| 在线观看国产h片| 美女cb高潮喷水在线观看| 亚洲av不卡在线观看| 99久久精品一区二区三区| 国产高清不卡午夜福利| 国产精品成人在线| 一本—道久久a久久精品蜜桃钙片| 大片电影免费在线观看免费| 久久精品久久精品一区二区三区| 少妇熟女欧美另类| 一级爰片在线观看| 亚洲最大成人中文| 亚洲激情五月婷婷啪啪| 2018国产大陆天天弄谢| 国产黄片视频在线免费观看| av国产久精品久网站免费入址| 亚洲成人一二三区av| 成人午夜精彩视频在线观看| 99视频精品全部免费 在线| 亚洲熟女精品中文字幕| 亚洲精品乱久久久久久| 少妇裸体淫交视频免费看高清| 91精品伊人久久大香线蕉| a 毛片基地| 亚洲精品国产av蜜桃| 国产在线男女| 国产av码专区亚洲av| 日本午夜av视频| 久久久久人妻精品一区果冻| 亚洲国产色片| 国产深夜福利视频在线观看| 国产成人午夜福利电影在线观看| 久久久久久久久久人人人人人人| 亚洲真实伦在线观看| 亚洲av.av天堂| 成人漫画全彩无遮挡| av女优亚洲男人天堂| 国产91av在线免费观看| 日日摸夜夜添夜夜添av毛片| 久久久午夜欧美精品| 色哟哟·www| 国产精品熟女久久久久浪| 精品亚洲成国产av| 精品国产一区二区三区久久久樱花 | 人体艺术视频欧美日本| 国产永久视频网站| 欧美精品国产亚洲| 亚洲四区av| 国产淫语在线视频| 久久久久久久久久人人人人人人| 蜜臀久久99精品久久宅男| 国产国拍精品亚洲av在线观看| 高清欧美精品videossex| 精品久久久噜噜| 天堂中文最新版在线下载| 午夜激情久久久久久久| 亚洲四区av| 成年美女黄网站色视频大全免费 | 99热这里只有精品一区| 伦精品一区二区三区| 久久久亚洲精品成人影院| 麻豆国产97在线/欧美| 久久久久精品性色| av视频免费观看在线观看| 中文字幕av成人在线电影| 最近最新中文字幕免费大全7| 国产免费视频播放在线视频| 欧美极品一区二区三区四区| 欧美激情极品国产一区二区三区 | 九色成人免费人妻av| 在线免费十八禁| 午夜免费男女啪啪视频观看| 91午夜精品亚洲一区二区三区| 男女国产视频网站| 下体分泌物呈黄色| 亚洲国产毛片av蜜桃av| 一级毛片aaaaaa免费看小| .国产精品久久| 中文字幕av成人在线电影| 精品午夜福利在线看| 免费观看的影片在线观看| 老师上课跳d突然被开到最大视频| 国产黄片美女视频| 久久久久久伊人网av| 国产免费一区二区三区四区乱码| 特大巨黑吊av在线直播| 国产精品不卡视频一区二区| 欧美精品一区二区大全| 国产亚洲欧美精品永久| 毛片女人毛片| 丝瓜视频免费看黄片| 日日啪夜夜爽| 在线看a的网站| 久久久精品94久久精品| 久久久久国产精品人妻一区二区| av专区在线播放| 男女边摸边吃奶| 国产精品精品国产色婷婷| 亚洲精品中文字幕在线视频 | 国产无遮挡羞羞视频在线观看| 久久久久久久久久成人| 亚洲真实伦在线观看| 国产高清有码在线观看视频| 国产69精品久久久久777片| 日韩在线高清观看一区二区三区| 这个男人来自地球电影免费观看 | 99久久精品国产国产毛片| 成年av动漫网址| 色网站视频免费| 免费不卡的大黄色大毛片视频在线观看| 我要看黄色一级片免费的| 免费人成在线观看视频色| 久久热精品热| 免费久久久久久久精品成人欧美视频 | 久热久热在线精品观看| 91狼人影院| 观看av在线不卡| 久久久久久伊人网av| 99久久综合免费| 中文在线观看免费www的网站| 久久亚洲国产成人精品v| 在线观看免费日韩欧美大片 | 亚洲精品一区蜜桃| 亚洲欧美日韩东京热| 亚洲精品自拍成人| 久久久精品94久久精品| 插阴视频在线观看视频| 免费黄网站久久成人精品| h视频一区二区三区| 永久免费av网站大全| av国产免费在线观看| 高清不卡的av网站| 极品少妇高潮喷水抽搐| www.av在线官网国产| 精品一区二区免费观看| 亚洲,一卡二卡三卡| 欧美精品亚洲一区二区| 亚洲人成网站高清观看| 伊人久久精品亚洲午夜| 汤姆久久久久久久影院中文字幕| 亚洲美女视频黄频| 久久青草综合色| 免费黄网站久久成人精品| 人人妻人人爽人人添夜夜欢视频 | 卡戴珊不雅视频在线播放| 国产精品一二三区在线看| 亚洲第一av免费看| 欧美+日韩+精品| 蜜臀久久99精品久久宅男| 国国产精品蜜臀av免费| 成年人午夜在线观看视频| 国产午夜精品一二区理论片| 亚洲第一区二区三区不卡| 亚洲欧美成人综合另类久久久| 亚洲色图综合在线观看| 少妇人妻精品综合一区二区| 男女无遮挡免费网站观看| 精品熟女少妇av免费看| 大码成人一级视频| 我要看黄色一级片免费的| 免费av中文字幕在线| 麻豆乱淫一区二区| 麻豆成人av视频| 国产精品一及| 午夜老司机福利剧场| 午夜免费观看性视频| 高清黄色对白视频在线免费看 | 亚洲经典国产精华液单| 久久99精品国语久久久| 天堂8中文在线网| 少妇熟女欧美另类| 97超碰精品成人国产| 永久免费av网站大全| 免费不卡的大黄色大毛片视频在线观看| 91久久精品国产一区二区成人| 丝瓜视频免费看黄片| 久久女婷五月综合色啪小说| 亚洲欧洲国产日韩| 毛片一级片免费看久久久久| 午夜精品国产一区二区电影| 亚洲国产日韩一区二区| 少妇被粗大猛烈的视频| 中文精品一卡2卡3卡4更新| 我要看黄色一级片免费的| 免费黄频网站在线观看国产| 肉色欧美久久久久久久蜜桃| 少妇的逼好多水| 春色校园在线视频观看| 国产一级毛片在线| 国产在线免费精品| 一区二区av电影网| 少妇猛男粗大的猛烈进出视频| 国产免费又黄又爽又色| 岛国毛片在线播放| 亚洲av成人精品一二三区| 一本久久精品| 久久久久久久大尺度免费视频| 久久99蜜桃精品久久| 免费av中文字幕在线| 亚洲欧美中文字幕日韩二区| 国产日韩欧美亚洲二区| 亚洲av日韩在线播放| 日日撸夜夜添| 99久久精品国产国产毛片| 精品亚洲成a人片在线观看 | 亚洲精品视频女| 国产高清有码在线观看视频| 天堂俺去俺来也www色官网| 亚洲精华国产精华液的使用体验| 汤姆久久久久久久影院中文字幕| 搡老乐熟女国产| 少妇熟女欧美另类| 熟女av电影| 99视频精品全部免费 在线| 丰满乱子伦码专区| 少妇人妻精品综合一区二区| kizo精华| 街头女战士在线观看网站| 精品一区在线观看国产| 99热全是精品| .国产精品久久| 日韩中文字幕视频在线看片 | 国产高清三级在线| 女性生殖器流出的白浆| 日韩精品有码人妻一区| 亚洲av免费高清在线观看| 免费人成在线观看视频色| 男的添女的下面高潮视频| 国产成人精品婷婷| 亚洲av电影在线观看一区二区三区| 肉色欧美久久久久久久蜜桃| 人人妻人人爽人人添夜夜欢视频 | 国产伦理片在线播放av一区| 久久精品熟女亚洲av麻豆精品| 亚洲三级黄色毛片| 欧美人与善性xxx| 亚洲精品日韩av片在线观看| 一本色道久久久久久精品综合| 夫妻午夜视频| 免费黄色在线免费观看| 秋霞伦理黄片| 亚洲一区二区三区欧美精品| 一级二级三级毛片免费看| 在线亚洲精品国产二区图片欧美 | 亚洲精品自拍成人| 欧美性感艳星| 国产v大片淫在线免费观看| 黄色一级大片看看| 爱豆传媒免费全集在线观看| 日韩亚洲欧美综合| 日韩视频在线欧美| 精品人妻偷拍中文字幕| 卡戴珊不雅视频在线播放| 国产乱人偷精品视频| 精品国产露脸久久av麻豆| 一本—道久久a久久精品蜜桃钙片| 亚洲va在线va天堂va国产| 一级毛片aaaaaa免费看小| 国产高潮美女av| 久久国产亚洲av麻豆专区| 成年av动漫网址| 三级国产精品欧美在线观看| 夫妻午夜视频| 性色av一级| 成人亚洲欧美一区二区av| 激情五月婷婷亚洲| 晚上一个人看的免费电影| 亚洲电影在线观看av| 欧美精品国产亚洲| 美女cb高潮喷水在线观看| 在现免费观看毛片| 亚洲av成人精品一区久久| 日日啪夜夜撸| 久久综合国产亚洲精品| 女的被弄到高潮叫床怎么办| 日韩人妻高清精品专区| 综合色丁香网| 欧美变态另类bdsm刘玥| 日韩欧美 国产精品| 美女福利国产在线 | 国产精品国产av在线观看| 国产精品无大码| 男人和女人高潮做爰伦理| 久久婷婷青草| 久久99热这里只有精品18| 久久久久网色| 少妇高潮的动态图| 免费不卡的大黄色大毛片视频在线观看| 欧美成人精品欧美一级黄| 久久久久久久久久人人人人人人| 久久综合国产亚洲精品| 三级国产精品片| 大陆偷拍与自拍| 免费av中文字幕在线| 国产精品一区www在线观看| 亚洲最大成人中文| 黄片无遮挡物在线观看| 天天躁日日操中文字幕| 久久精品熟女亚洲av麻豆精品| 两个人的视频大全免费| 国产成人a∨麻豆精品| 久久久精品免费免费高清| 成人毛片a级毛片在线播放| 一二三四中文在线观看免费高清| 丝袜喷水一区| 午夜日本视频在线| 在现免费观看毛片| 人人妻人人澡人人爽人人夜夜| 亚洲美女黄色视频免费看| 王馨瑶露胸无遮挡在线观看| 男女无遮挡免费网站观看| 亚洲精品久久久久久婷婷小说| 成年人午夜在线观看视频| 尾随美女入室| 国产精品欧美亚洲77777| 久久久久久久大尺度免费视频| 国产男女超爽视频在线观看| 亚洲精品国产成人久久av| 国产精品国产三级国产专区5o| 女性生殖器流出的白浆| 我要看日韩黄色一级片| 国产爽快片一区二区三区| 久久女婷五月综合色啪小说| 日韩强制内射视频| 亚洲人成网站在线播| 人体艺术视频欧美日本| 国产精品一区二区在线不卡| 欧美高清性xxxxhd video| av播播在线观看一区| 欧美日韩国产mv在线观看视频 | 精品久久久久久电影网| 2021少妇久久久久久久久久久| 日韩欧美一区视频在线观看 | 精品一区二区三卡| 久久99精品国语久久久| 久久 成人 亚洲| 秋霞伦理黄片| 亚洲美女黄色视频免费看| 一本久久精品| 色5月婷婷丁香| 久久国产乱子免费精品| 少妇猛男粗大的猛烈进出视频| 国产日韩欧美在线精品| 国产黄色免费在线视频| 亚洲精品国产av成人精品| 亚洲在久久综合| 麻豆国产97在线/欧美| 国产成人精品一,二区| 少妇丰满av| 久久久久国产精品人妻一区二区| 人妻夜夜爽99麻豆av| 欧美少妇被猛烈插入视频| 男女无遮挡免费网站观看| 3wmmmm亚洲av在线观看| 美女高潮的动态| 欧美xxⅹ黑人| 91在线精品国自产拍蜜月| 麻豆成人av视频| 一边亲一边摸免费视频| 国产精品麻豆人妻色哟哟久久| 亚洲激情五月婷婷啪啪| 女性被躁到高潮视频| 日本黄色日本黄色录像| 91久久精品国产一区二区成人| 日韩电影二区| 成人黄色视频免费在线看| 超碰97精品在线观看| 全区人妻精品视频| 最近中文字幕高清免费大全6| 另类亚洲欧美激情| 精品久久久久久电影网| 九九爱精品视频在线观看| 婷婷色麻豆天堂久久| 一级黄片播放器| 伦理电影大哥的女人| 亚洲国产欧美人成| 老女人水多毛片| 日产精品乱码卡一卡2卡三| 国内精品宾馆在线| 99久久精品一区二区三区| 免费av中文字幕在线| 中文在线观看免费www的网站| 久久热精品热| 亚洲va在线va天堂va国产| 在线精品无人区一区二区三 | 简卡轻食公司| 永久网站在线| 一二三四中文在线观看免费高清| 伦精品一区二区三区| 国产亚洲91精品色在线| 蜜桃久久精品国产亚洲av| 欧美亚洲 丝袜 人妻 在线| 我要看黄色一级片免费的| 黄色一级大片看看| 少妇猛男粗大的猛烈进出视频| 久久久久久久亚洲中文字幕| 久久久成人免费电影| 狂野欧美白嫩少妇大欣赏| 日韩人妻高清精品专区| 国模一区二区三区四区视频| 国产白丝娇喘喷水9色精品| 91在线精品国自产拍蜜月| 亚洲经典国产精华液单| 多毛熟女@视频| 久久久精品免费免费高清| 久久久久视频综合| 中文字幕亚洲精品专区| 男人舔奶头视频| 成人二区视频| 亚洲欧美精品专区久久|