• <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)流化床鍋爐集控運行研究
    69av精品久久久久久| 麻豆精品久久久久久蜜桃| 欧美成人一区二区免费高清观看| 欧美日韩一区二区视频在线观看视频在线 | 国产爱豆传媒在线观看| www.av在线官网国产| 2022亚洲国产成人精品| 国产视频首页在线观看| 久99久视频精品免费| 国产 一区 欧美 日韩| 一个人看视频在线观看www免费| 久久精品人妻少妇| 色哟哟·www| 久久国内精品自在自线图片| 又黄又爽又刺激的免费视频.| 九草在线视频观看| 久久精品夜夜夜夜夜久久蜜豆| 一级毛片久久久久久久久女| 国产69精品久久久久777片| 一级毛片aaaaaa免费看小| 国产老妇女一区| 男女下面进入的视频免费午夜| 少妇裸体淫交视频免费看高清| 亚洲精品成人av观看孕妇| 国产色婷婷99| 国产精品综合久久久久久久免费| 精品久久久久久久久av| 亚洲精品第二区| 一级爰片在线观看| 黄色日韩在线| 成人美女网站在线观看视频| 纵有疾风起免费观看全集完整版 | 国产精品伦人一区二区| 黄色配什么色好看| 亚洲欧美成人综合另类久久久| 一夜夜www| 51国产日韩欧美| 欧美高清性xxxxhd video| 美女高潮的动态| 国产不卡一卡二| 少妇丰满av| 国产成人福利小说| 国产爱豆传媒在线观看| 日本wwww免费看| 精品一区在线观看国产| videossex国产| 韩国av在线不卡| 特大巨黑吊av在线直播| 最近最新中文字幕大全电影3| 80岁老熟妇乱子伦牲交| 看免费成人av毛片| 成年女人看的毛片在线观看| 久久99蜜桃精品久久| 亚洲国产av新网站| 久久久久久久午夜电影| 69人妻影院| 亚洲精品日本国产第一区| 成人二区视频| 免费黄频网站在线观看国产| 免费av观看视频| 国产白丝娇喘喷水9色精品| av在线蜜桃| 男人舔女人下体高潮全视频| av一本久久久久| 亚洲国产精品成人综合色| 一本一本综合久久| 麻豆av噜噜一区二区三区| 亚洲国产精品国产精品| 日日摸夜夜添夜夜添av毛片| 三级男女做爰猛烈吃奶摸视频| 婷婷色麻豆天堂久久| 亚洲欧洲国产日韩| 日韩强制内射视频| 日韩欧美精品免费久久| 精品人妻一区二区三区麻豆| 我的老师免费观看完整版| 777米奇影视久久| 久久久色成人| 成年女人在线观看亚洲视频 | 日韩中字成人| 国语对白做爰xxxⅹ性视频网站| 伊人久久国产一区二区| 男的添女的下面高潮视频| 床上黄色一级片| 午夜激情欧美在线| 欧美不卡视频在线免费观看| 国产 亚洲一区二区三区 | 亚洲高清免费不卡视频| 国产亚洲av片在线观看秒播厂 | 久久精品久久精品一区二区三区| 国产综合懂色| 欧美三级亚洲精品| 狂野欧美激情性xxxx在线观看| 午夜激情福利司机影院| 大香蕉久久网| 亚洲自偷自拍三级| 在线免费十八禁| 午夜免费观看性视频| 中文字幕av成人在线电影| 综合色av麻豆| 免费看av在线观看网站| 国产高清不卡午夜福利| 国产黄a三级三级三级人| 熟女人妻精品中文字幕| 免费看光身美女| 婷婷色综合大香蕉| 高清欧美精品videossex| 免费少妇av软件| 两个人的视频大全免费| 国产综合精华液| 99久国产av精品| 看十八女毛片水多多多| 久久人人爽人人爽人人片va| 久热久热在线精品观看| 全区人妻精品视频| 97在线视频观看| 日韩一区二区三区影片| 国产精品一区二区在线观看99 | 免费黄网站久久成人精品| 欧美最新免费一区二区三区| 亚洲18禁久久av| 久99久视频精品免费| 日日摸夜夜添夜夜添av毛片| 色综合亚洲欧美另类图片| 亚洲最大成人手机在线| 久久久久免费精品人妻一区二区| 国产成人91sexporn| 午夜福利视频1000在线观看| 午夜视频国产福利| 久久精品久久久久久噜噜老黄| 国产成人a∨麻豆精品| av线在线观看网站| 尾随美女入室| 国产老妇女一区| 亚洲av日韩在线播放| 只有这里有精品99| 在线免费观看的www视频| 欧美xxxx黑人xx丫x性爽| 国产色爽女视频免费观看| 纵有疾风起免费观看全集完整版 | 国产在线一区二区三区精| 老女人水多毛片| 偷拍熟女少妇极品色| 性插视频无遮挡在线免费观看| 成年人午夜在线观看视频 | 国产精品久久久久久久电影| 国产亚洲最大av| 女人被狂操c到高潮| 欧美不卡视频在线免费观看| 国产亚洲午夜精品一区二区久久 | 黄色配什么色好看| 高清av免费在线| 国产片特级美女逼逼视频| 看黄色毛片网站| 国产一区二区在线观看日韩| 中文乱码字字幕精品一区二区三区 | 日韩精品有码人妻一区| 日韩 亚洲 欧美在线| 国产成人一区二区在线| 免费播放大片免费观看视频在线观看| 亚洲va在线va天堂va国产| 成年女人在线观看亚洲视频 | 国精品久久久久久国模美| 精品人妻一区二区三区麻豆| 日产精品乱码卡一卡2卡三| 91久久精品国产一区二区成人| 国产一区有黄有色的免费视频 | 国产精品嫩草影院av在线观看| 看黄色毛片网站| 久久精品久久久久久噜噜老黄| 日本欧美国产在线视频| 日韩国内少妇激情av| 免费电影在线观看免费观看| 国产亚洲一区二区精品| 成年女人看的毛片在线观看| 免费观看性生交大片5| 亚洲天堂国产精品一区在线| 日日啪夜夜撸| 亚洲最大成人手机在线| 欧美激情久久久久久爽电影| 男女边吃奶边做爰视频| 亚洲不卡免费看| 日本黄大片高清| av在线天堂中文字幕| 成年av动漫网址| 国产高潮美女av| 日本爱情动作片www.在线观看| 最近最新中文字幕大全电影3| 国产极品天堂在线| 国产成人a∨麻豆精品| 美女被艹到高潮喷水动态| 日韩欧美精品v在线| 成人无遮挡网站| 99久久精品一区二区三区| 看非洲黑人一级黄片| 久久国内精品自在自线图片| 成人高潮视频无遮挡免费网站| 亚洲自偷自拍三级| 国产91av在线免费观看| 日产精品乱码卡一卡2卡三| 午夜福利在线观看吧| 亚洲丝袜综合中文字幕| 亚洲av男天堂| 免费看日本二区| av.在线天堂| 18禁在线播放成人免费| 免费少妇av软件| 美女大奶头视频| 高清欧美精品videossex| 国产伦一二天堂av在线观看| 一级片'在线观看视频| 婷婷色av中文字幕| 国产女主播在线喷水免费视频网站 | 99久久精品一区二区三区| 男人舔女人下体高潮全视频| 建设人人有责人人尽责人人享有的 | 人妻少妇偷人精品九色| 九九在线视频观看精品| 观看美女的网站| 最近最新中文字幕大全电影3| 国产黄频视频在线观看| 亚洲最大成人手机在线| 一个人看的www免费观看视频| 成人特级av手机在线观看| 内地一区二区视频在线| 2022亚洲国产成人精品| 国产精品久久久久久精品电影小说 | 中文字幕亚洲精品专区| 亚洲真实伦在线观看| 能在线免费看毛片的网站| 尾随美女入室| 一级毛片电影观看| 国产成人精品一,二区| 国产真实伦视频高清在线观看| 久久久久九九精品影院| 黄色一级大片看看| 视频中文字幕在线观看| 日韩欧美一区视频在线观看 | 午夜免费观看性视频| 久久这里只有精品中国| 极品教师在线视频| 亚洲精品久久久久久婷婷小说| av国产久精品久网站免费入址| 人人妻人人看人人澡| 九九久久精品国产亚洲av麻豆| 最近的中文字幕免费完整| 波野结衣二区三区在线| 日本av手机在线免费观看| 夜夜爽夜夜爽视频| 丰满人妻一区二区三区视频av| 麻豆av噜噜一区二区三区| 寂寞人妻少妇视频99o| 国产免费视频播放在线视频 | 18禁在线无遮挡免费观看视频| 午夜视频国产福利| 久久久国产一区二区| 亚洲av电影不卡..在线观看| 中文字幕av成人在线电影| 欧美日韩国产mv在线观看视频 | 欧美变态另类bdsm刘玥| 国产成年人精品一区二区| 国产v大片淫在线免费观看| 国内精品美女久久久久久| 成人特级av手机在线观看| 七月丁香在线播放| 国产精品av视频在线免费观看| 综合色丁香网| 国产一级毛片在线| 特级一级黄色大片| 97人妻精品一区二区三区麻豆| 一区二区三区四区激情视频| 午夜精品在线福利| 亚洲欧美日韩东京热| 五月伊人婷婷丁香| 欧美xxxx黑人xx丫x性爽| 午夜精品一区二区三区免费看| 全区人妻精品视频| 热99在线观看视频| 女人十人毛片免费观看3o分钟| 国产精品久久久久久久久免| 欧美激情国产日韩精品一区| 久久99蜜桃精品久久| 99九九线精品视频在线观看视频| 秋霞在线观看毛片| 在线播放无遮挡| 欧美日韩在线观看h| 久久6这里有精品| 久久久久久久久久久丰满| 欧美日韩综合久久久久久| 最近最新中文字幕免费大全7| 亚洲精品国产av成人精品| 国产免费又黄又爽又色| 国产极品天堂在线| 18禁在线无遮挡免费观看视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美xxxx性猛交bbbb| 国产精品久久久久久久电影| 最后的刺客免费高清国语| 国产毛片a区久久久久| 99九九线精品视频在线观看视频| 国产亚洲精品av在线| 精品久久国产蜜桃| 亚洲丝袜综合中文字幕| 亚洲内射少妇av| 又黄又爽又刺激的免费视频.| 亚洲av成人精品一二三区| 亚洲精品中文字幕在线视频 | 欧美xxxx黑人xx丫x性爽| 少妇熟女aⅴ在线视频| 亚洲av电影在线观看一区二区三区 | 91狼人影院| 少妇的逼好多水| 亚洲天堂国产精品一区在线| 七月丁香在线播放| 老师上课跳d突然被开到最大视频| 26uuu在线亚洲综合色| 女人被狂操c到高潮| 亚洲av中文av极速乱| 成人av在线播放网站| 国产有黄有色有爽视频| 能在线免费看毛片的网站| 成年版毛片免费区| 一夜夜www| 国产伦理片在线播放av一区| 深爱激情五月婷婷| 久久久久免费精品人妻一区二区| av在线天堂中文字幕| 看黄色毛片网站| 国产黄a三级三级三级人| 亚洲av成人精品一区久久| 亚洲国产av新网站| ponron亚洲| 国产在视频线精品| 精品国产露脸久久av麻豆 | 天天一区二区日本电影三级| 国产高清有码在线观看视频| 国产午夜精品论理片| 噜噜噜噜噜久久久久久91| 国产淫语在线视频| 少妇的逼好多水| 天堂网av新在线| 国产成人免费观看mmmm| 天堂影院成人在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人一区二区视频在线观看| 亚洲欧洲日产国产| 看黄色毛片网站| 男女国产视频网站| 亚洲国产日韩欧美精品在线观看| 国产综合精华液| 国产成人一区二区在线| 免费大片18禁| 99re6热这里在线精品视频| 国产高清国产精品国产三级 | 国产 一区 欧美 日韩| 免费av毛片视频| 成年av动漫网址| 男人舔女人下体高潮全视频| 色尼玛亚洲综合影院| 久久精品久久久久久噜噜老黄| 亚洲国产成人一精品久久久| 亚洲美女视频黄频| 午夜久久久久精精品| 国产成人免费观看mmmm| 国产成人a区在线观看| 日日撸夜夜添| 国产高清三级在线| 2018国产大陆天天弄谢| 欧美潮喷喷水| 嫩草影院精品99| 成人午夜精彩视频在线观看| a级毛片免费高清观看在线播放| 可以在线观看毛片的网站| 亚洲av成人av| 99久久精品一区二区三区| 国产精品日韩av在线免费观看| 日韩电影二区| or卡值多少钱| 搞女人的毛片| 男女视频在线观看网站免费| 九九在线视频观看精品| 在线 av 中文字幕| 日韩精品青青久久久久久| 非洲黑人性xxxx精品又粗又长| 国产成人精品久久久久久| 日日摸夜夜添夜夜爱| 午夜精品一区二区三区免费看| 国产一级毛片在线| 精品国产一区二区三区久久久樱花 | 免费av观看视频| 日韩精品青青久久久久久| 亚洲国产精品成人综合色| 青春草视频在线免费观看| 免费观看在线日韩| 久久综合国产亚洲精品| 亚洲国产精品国产精品| 免费高清在线观看视频在线观看| 亚洲精品影视一区二区三区av| 99久久精品一区二区三区| 又大又黄又爽视频免费| 狂野欧美白嫩少妇大欣赏| 黄色日韩在线| 偷拍熟女少妇极品色| 精品久久久噜噜| 欧美日韩精品成人综合77777| 男人爽女人下面视频在线观看| 国产伦理片在线播放av一区| 国产精品一区二区三区四区久久| 亚洲久久久久久中文字幕| 国产伦精品一区二区三区视频9| 国产久久久一区二区三区| 亚洲怡红院男人天堂| 乱码一卡2卡4卡精品| 久久久久久久久久黄片| 老师上课跳d突然被开到最大视频| 国产精品人妻久久久影院| 成人综合一区亚洲| a级毛片免费高清观看在线播放| 久久亚洲国产成人精品v| 深爱激情五月婷婷| 成年版毛片免费区| 毛片女人毛片| 亚洲av国产av综合av卡| 国产精品久久久久久久久免| 日日撸夜夜添| 成人午夜精彩视频在线观看| 国产黄色小视频在线观看| 国产白丝娇喘喷水9色精品| 性色avwww在线观看| 精品一区二区三区人妻视频| 国产精品国产三级国产专区5o| 日本色播在线视频| 亚洲图色成人| 2018国产大陆天天弄谢| 人妻一区二区av| 久久99热这里只频精品6学生| 少妇熟女aⅴ在线视频| 身体一侧抽搐| 大香蕉久久网| 欧美bdsm另类| 老司机影院成人| 少妇高潮的动态图| 91久久精品国产一区二区三区| 黄色欧美视频在线观看| 最新中文字幕久久久久| 蜜桃亚洲精品一区二区三区| 自拍偷自拍亚洲精品老妇| 久久久久久伊人网av| 毛片女人毛片| 亚洲欧美精品自产自拍| 插阴视频在线观看视频| 26uuu在线亚洲综合色| 国产黄色小视频在线观看| 人人妻人人看人人澡| 欧美另类一区| 精品久久久久久久人妻蜜臀av| 天堂√8在线中文| 午夜精品一区二区三区免费看| 爱豆传媒免费全集在线观看| 亚洲成人久久爱视频| 亚洲欧美一区二区三区黑人 | 97热精品久久久久久| 大香蕉久久网| 亚洲aⅴ乱码一区二区在线播放| 久久鲁丝午夜福利片| 视频中文字幕在线观看| 色5月婷婷丁香| 久久久精品欧美日韩精品| 亚洲欧美成人精品一区二区| 亚洲精品国产成人久久av| 国产欧美日韩精品一区二区| 又粗又硬又长又爽又黄的视频| 亚洲天堂国产精品一区在线| 久久国产乱子免费精品| 欧美3d第一页| 国产极品天堂在线| 99热全是精品| 国产精品国产三级专区第一集| 99热这里只有是精品在线观看| 亚洲真实伦在线观看| 久久久久性生活片| 日日啪夜夜撸| 熟女人妻精品中文字幕| 22中文网久久字幕| 五月天丁香电影| 欧美精品一区二区大全| 成年人午夜在线观看视频 | 我要看日韩黄色一级片| 久久精品熟女亚洲av麻豆精品 | 国产精品一二三区在线看| 中文字幕av成人在线电影| 午夜激情福利司机影院| 免费高清在线观看视频在线观看| 久久99热6这里只有精品| 国产精品福利在线免费观看| 麻豆av噜噜一区二区三区| 亚洲乱码一区二区免费版| 三级国产精品片| 日本与韩国留学比较| 美女高潮的动态| 97人妻精品一区二区三区麻豆| 国产v大片淫在线免费观看| 国产不卡一卡二| 亚洲精品乱码久久久v下载方式| 欧美日韩精品成人综合77777| 亚洲av电影不卡..在线观看| 又大又黄又爽视频免费| 国产在线一区二区三区精| 国产精品一及| 久久久a久久爽久久v久久| 七月丁香在线播放| 十八禁国产超污无遮挡网站| 日韩三级伦理在线观看| 免费观看精品视频网站| 男女边吃奶边做爰视频| 97超碰精品成人国产| 国产综合精华液| 欧美激情久久久久久爽电影| 在线观看免费高清a一片| 夫妻性生交免费视频一级片| 色播亚洲综合网| 成年版毛片免费区| 六月丁香七月| 美女脱内裤让男人舔精品视频| 亚洲va在线va天堂va国产| 亚洲av日韩在线播放| 水蜜桃什么品种好| 亚洲精品乱码久久久久久按摩| 非洲黑人性xxxx精品又粗又长| 久久99蜜桃精品久久| 欧美3d第一页| 国产不卡一卡二| 亚洲精品一区蜜桃| 美女黄网站色视频| 联通29元200g的流量卡| 国产黄a三级三级三级人| 国产在线男女| 丰满乱子伦码专区| 在线播放无遮挡| 日本免费a在线| 亚洲欧美日韩东京热| 国产精品久久久久久av不卡| 美女脱内裤让男人舔精品视频| 亚洲四区av| 欧美三级亚洲精品| av女优亚洲男人天堂| 日日啪夜夜撸| 蜜桃亚洲精品一区二区三区| 成人国产麻豆网| 真实男女啪啪啪动态图| 高清av免费在线| 国产日韩欧美在线精品| 高清日韩中文字幕在线| av在线老鸭窝| 国产精品一区二区三区四区久久| 国产欧美日韩精品一区二区| 最近视频中文字幕2019在线8| 精品少妇黑人巨大在线播放| 91久久精品国产一区二区成人| 大陆偷拍与自拍| 亚洲最大成人手机在线| 22中文网久久字幕| 久久久久精品性色| 亚洲av国产av综合av卡| 狂野欧美白嫩少妇大欣赏| 自拍偷自拍亚洲精品老妇| 国产精品蜜桃在线观看| 婷婷六月久久综合丁香| 国产日韩欧美在线精品| 久久精品综合一区二区三区| 国产午夜精品论理片| 一二三四中文在线观看免费高清| 精品国产露脸久久av麻豆 | 日韩成人伦理影院| 国产黄片视频在线免费观看| 免费av不卡在线播放| 日韩亚洲欧美综合| 22中文网久久字幕| 女人被狂操c到高潮| 最新中文字幕久久久久| 老女人水多毛片| 国产精品日韩av在线免费观看| 91狼人影院| 国产精品久久久久久av不卡| 日本-黄色视频高清免费观看| 97精品久久久久久久久久精品| 2021天堂中文幕一二区在线观| 精品99又大又爽又粗少妇毛片| 亚洲av男天堂| 免费少妇av软件| 亚洲最大成人av| 免费大片黄手机在线观看| 亚洲综合精品二区| ponron亚洲| 婷婷色av中文字幕| www.色视频.com| 午夜精品一区二区三区免费看| 亚洲精品国产成人久久av| 成人高潮视频无遮挡免费网站| 一区二区三区乱码不卡18| 亚洲精品aⅴ在线观看| 韩国高清视频一区二区三区| 免费播放大片免费观看视频在线观看| 观看免费一级毛片| a级毛片免费高清观看在线播放| 久久久久国产网址| 最近手机中文字幕大全| ponron亚洲| 晚上一个人看的免费电影| 免费不卡的大黄色大毛片视频在线观看 | 欧美xxⅹ黑人| 非洲黑人性xxxx精品又粗又长| 国产一区二区在线观看日韩| 国产成人a区在线观看| 美女xxoo啪啪120秒动态图|