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

    攪拌對冰漿存儲冰晶粒徑演化的影響

    2017-10-13 12:13:40唐藝芳劉志強王肖肖康威張靜雅
    中南大學學報(自然科學版) 2017年5期
    關鍵詞:槽內冰晶槳葉

    唐藝芳,劉志強, 2,王肖肖,康威,張靜雅

    ?

    攪拌對冰漿存儲冰晶粒徑演化的影響

    唐藝芳1,劉志強1, 2,王肖肖1,康威1,張靜雅1

    (1. 中南大學能源科學與工程學院,湖南長沙,410083;2. 湖南省建筑節(jié)能與環(huán)境控制關鍵技術協(xié)同創(chuàng)新中心,湖南株洲,412007)

    為研究攪拌作用下冰漿存儲過程中冰晶粒徑演化的規(guī)律,基于冰晶動力學事件核模型,利用Fluent軟件、耦合多相流模型與數(shù)群平衡模型對冰漿存儲過程中冰晶粒徑演化進行模擬,將平均粒徑模擬值與實驗值進行對比分析,并對模型有效性進行驗證?;谛拚谋恿W事件核模型對攪拌作用下冰漿存儲過程中的冰晶粒徑演化進行模擬,分析單螺旋槳槳葉直徑及位置、雙螺旋槳槳葉間距對冰漿存儲過程中冰晶體積分數(shù)分布、冰晶粒徑數(shù)量密度分布及冰晶平均粒徑演化的影響。研究結果表明:槳葉直徑對冰晶體積分數(shù)分布、冰晶粒徑數(shù)量密度分布及冰晶平均粒徑均有影響;槳葉間距對冰晶體積分數(shù)分布影響顯著;雙螺旋槳能夠使蓄冰槽內冰晶分布更加均勻,更能抑制冰晶粒徑增長。

    儲能;冰漿;多相流模型;數(shù)群平衡模型

    冰漿是由載流溶液和冰晶組成的固液兩相流,具有流動與傳熱特性良好、蓄能密度大、無污染等優(yōu)點,現(xiàn)已應用于區(qū)域供冷、蓄冷空調系統(tǒng)及工業(yè)冷卻等領域,且在化工、醫(yī)藥、環(huán)保等領域具有廣闊的應用前景[1]。冰漿在存儲過程中,在成核、生長、破碎、團聚等作用下,冰晶粒徑及分布會隨著時間發(fā)生變化。冰晶粒徑及分布對冰漿流動與換熱特性有重要影響,如冰晶的富集容易導致冰漿釋冷速度慢且不能完全融化、泵送困難等問題[2?4],人們對冰晶粒徑控制等方面進行了研究,目前研究主要集中在冰晶粒徑生長抑制添加劑種類及濃度方面,冰漿存儲過程特性研究則主要集中在冰漿的富集與分層等宏觀方面[5?8]。DELAHAYE等[9]采取向蓄冰槽內通入空氣的措施抑制冰晶富集,但該方法容易造成設備腐蝕。加入添加劑溶液能抑制冰晶聚集增長,但會明顯降低溶液冷凝溫度。而在蓄冰槽中增加擾動能夠有效地抑制冰晶聚集,同時又不會引起腐蝕設備與降低溶液冷凝溫度等問題[3],因此,螺旋槳攪拌得到了較廣泛應用。目前,基于微觀層面對冰晶粒徑演化已有少量研究。PRONK等[10?11]運用結晶動力學理論對冰漿存儲過程冰晶粒徑演化機理進行理論分析和實驗研究,理論上認為冰晶粒徑分布主要是Ostwald熟化作用的結果[3]。趙騰磊等[12?13]運用數(shù)群平衡模型研究了破碎與團聚對冰漿存儲過程冰晶粒徑演化的影響,但均未考慮存儲方式及空間等因素的影響。目前,基于微觀層面的冰晶粒徑演化的研究較少,且已有的冰晶聚集核模型及破碎核模型的理論計算結果與實驗結果之間仍存在明顯差異。數(shù)群平衡模型(population balance model, PBM)是模擬離散系統(tǒng)中離散單元動力學演化過程的有效方法[14],本文綜合考慮冰晶的成核、生長、團聚及破碎等微觀動力學行為,基于已有的冰晶動力學事件核模型建立冰晶數(shù)群平衡模型,耦合多相流模型與數(shù)群平衡模型對冰漿存儲過程中冰晶粒徑演化進行數(shù)值模擬,并對模型的有效性進行驗證,最后著重分析單螺旋槳槳葉直徑及位置、雙螺旋槳槳葉間距對冰漿存儲過程中冰晶體積分數(shù)分布、冰晶粒徑數(shù)量密度分布及冰晶平均粒徑演化的影響。

    1 模型

    1.1 多相流模型

    冰漿在存儲過程中,冰晶粒子懸浮在溶液中,不會發(fā)生兩相完全分層的現(xiàn)象,故選用Euler模型描述冰漿流動。

    1.2 冰晶動力學事件核模型

    1.2.1 冰晶成核模型

    MICHAEL等[3,15]對冰晶初次成核進行了理論分析與實驗研究,得出冰晶初次成核率的計算式為

    同質成核時,的表達式為

    (2)

    異質成核時,?的表達式為

    式中:L為單位體積內分子數(shù);為波爾茲曼常數(shù),J/K;為溶液熱力學溫度,K;為普朗克常數(shù),J?s;?A為水分子活化能,J/kg;?為成核能,J/kg;iw為冰水之間表面能,mJ/m2;eq為溶液相變平衡熱力學溫度,K;ice為冰晶的密度,kg/m3;f為凝固潛熱,J/kg;?為過冷度,K;()為晶核在異質基體表面的接觸角、異質基體粒度與臨界晶核粒度之比的函數(shù)。

    1.2.2 冰晶生長核模型

    目前應用較多的是PRONK[16]提出的計算冰晶生長率模型,該模型能較準確地計算冰晶的生長率:

    式中:V和A分別為體積常數(shù)和面積常數(shù);b為溶液中水的質量分數(shù),%;d為質量傳遞常數(shù),m/s;lid為溶液密度,kg/m3;m為熱量傳遞函數(shù),W/(m2?K)。

    1.2.3 冰晶聚集核模型

    用布朗聚集核模型表示體積為的顆粒團聚為體積為的速率表達式為

    式中:為冰漿中載流溶液的動力黏度,Pa?s。

    1.2.4 冰晶破碎核模型

    粒子的破碎受到機械攪動、溶液物性等因素影響,可用下式計算冰晶的破碎率[17]:

    (7)

    式中:為粒子的形狀因子,本文假設粒子為球形,即=1。

    2 模型驗證

    基于PRONK的冰漿存儲實驗[8]對模型進行驗證。依據其中實驗3,建立無攪拌蓄冰槽數(shù)值模型,直徑為640 mm、高為860 mm的圓柱形蓄冰槽內為非均勻存儲的冰漿溶液。依據其中實驗13,建立有攪拌蓄冰槽數(shù)值模型,直徑為125 mm、高為130 mm的圓柱形蓄冰槽內為均勻存儲的冰漿溶液,蓄冰槽內配備有螺旋槳攪拌器。分別選用顆粒動力學理論模型、標準?混合湍流模型描述冰晶粒子間的相互作用及冰漿湍流作用,建立流體力學模型。

    2.1 模擬條件

    根據實驗條件對模型物性參數(shù)進行設置,冰漿的物性參數(shù)如表1所示。

    表1 冰漿的物性參數(shù)

    2.1.1 邊界條件

    1) 假設冰漿存儲過程中與外界沒有質量和能量交換,蓄冰槽壁面可以認為是理想絕熱壁面。2) 假設冰晶成核發(fā)生在最小粒徑處,添加劑為乙醇、氯化鈉時,冰晶粒徑計算區(qū)域分別取0.01~900 μm和0.01~ 1 200 μm。不計計算區(qū)域以外的冰晶顆粒,即假設在計算區(qū)域以外的冰晶粒徑數(shù)量密度為0 個/m4。

    2.1.2 初始條件

    1) 假設初始時刻每個空間網格節(jié)點處冰晶粒徑分布相同,將實驗測定的冰晶粒徑分布進行擬合得到初始時刻冰晶粒徑數(shù)量密度的對數(shù)正態(tài)分布的數(shù)學表達式如下:

    式中:為冰晶粒徑,μm;0為初始時刻冰晶平均粒徑,μm。0根據實驗結果給定,實驗3、實驗13中初始時刻冰晶平均粒徑分別為97 μm和207 μm。

    2) 假定冰晶在粒徑最小處的成核率為1 m-3/s。冰晶的生長率與冰晶粒徑相關,但在同一種添加劑溶液、過冷度一定的條件下其變化不大,故假設不同粒徑的冰晶生長率相同,為10?10 m/s。

    2.2 模型驗證

    冰晶在布朗擴散作用下會發(fā)生碰撞導致破碎,粒徑越大,其破碎作用越明顯,故應考慮布朗擴散的破碎作用,引入布朗擴散耗損率0。布朗擴散作用下破碎率公式為

    式中:0為布朗擴散耗損率,其值為0.000 1 W/kg。

    通過選用合適的參數(shù),得到無攪拌作用下冰漿存儲2 h后冰晶平均粒徑實驗擬合曲線與數(shù)值模擬曲線,如圖1所示。從圖1可見:冰晶平均粒徑實驗擬合曲線與數(shù)值模擬曲線基本吻合,驗證了模型的有效性。

    圖1 無攪拌作用下冰漿存儲2 h后冰晶平均粒徑演化的實驗擬合結果與數(shù)值模擬結果

    在攪拌作用下,冰晶的聚集受布朗作用與湍流作用影響,故應考慮湍流聚集率,引入湍流聚集參數(shù)3。

    通過選用合適的破碎率參數(shù)及湍流聚集參數(shù),得到攪拌作用下冰漿存儲2 h后冰晶平均粒徑演化實驗擬合曲線與數(shù)值模擬曲線如圖2所示,可見兩者吻合較好,進一步驗證了模型的有效性。

    聚集核模型和破碎核模型的選用參數(shù)如下:破碎參數(shù)1為4.81×10?6;破碎參數(shù)2為0.080 0;湍流聚集參數(shù)3為0.030 9。

    圖2 攪拌作用下冰漿存儲2 h后冰晶平均粒徑演化的實驗擬合結果與數(shù)值模擬結果

    3 結果與分析

    3.1 單螺旋槳對冰漿存儲的影響

    按螺旋槳槳葉直徑、位置不同,設置CaseⅠ~ Case Ⅳ共4種槳葉參數(shù),見表2。

    表2 單螺旋槳參數(shù)

    3.1.1 槳葉直徑對冰漿存儲的影響

    在CaseⅠ,Case Ⅱ和Case Ⅲ下,當存儲時間分別為1.5 h和3.0 h時蓄冰槽內冰晶體積分數(shù)分布見 圖3。

    從圖3可見:現(xiàn)槳葉直徑越大,蓄冰槽上部和槳葉下邊緣冰晶體積分數(shù)越小,蓄冰槽底部冰晶體積分數(shù)越大,頂部和底部冰晶體積分數(shù)之差越小。由于槳葉直徑越大,其旋轉造成周圍流體的流動越劇烈且擾動范圍越大,冰晶向上的流動速度隨之增大,冰晶粒子的破碎作用顯著增強,因此,蓄冰槽上部冰晶粒子數(shù)量會隨槳葉直徑增大而增大,但冰晶總體積反而減小,冰晶富集現(xiàn)象減弱。

    在Case Ⅰ,Case Ⅱ和Case Ⅲ下,冰晶粒徑數(shù)量密度分布曲線見圖4。從圖4可知:槳葉直徑不同,但冰晶數(shù)量密度分布曲線均類似于對數(shù)正態(tài)分布曲線,且隨時間逐漸向右推移,出現(xiàn)由“高窄”到“矮寬”的變化,只是推移的程度不同。對比存儲時間= 3 h可發(fā)現(xiàn):槳葉直徑越大,同時刻小粒徑冰晶數(shù)量密度越大,大粒徑冰晶數(shù)量密度越?。槐Я綌?shù)量密度的峰值越大,其數(shù)量密度峰值對應的粒徑越小;當槳葉直徑增大時,冰晶粒徑數(shù)量密度分布曲線向右推移會出現(xiàn)遲緩現(xiàn)象,且槳葉直徑越大,遲緩現(xiàn)象越嚴重。這是由于槳葉直徑越大其擾動作用越強,冰晶粒子之間、粒子與蓄冰槽壁面的碰撞概率及能量越大,冰晶碰撞后破碎概率也就越大,因此,螺旋槳槳葉直徑越大,冰晶粒子破碎作用越強。但聚集作用始終大于破碎作用,只是兩者之差會逐漸減小。由此可見,增大槳葉直徑能夠增強對冰晶粒徑增長的抑制作用。

    (a) Case Ⅰ,1.5 h;(b) Case Ⅰ,3.0 h;(c) Case Ⅱ,1.5 h;(d) Case Ⅱ,3.0 h;(e) Case Ⅲ,1.5 h;(f) Case Ⅲ,3.0 h

    1—3 h, d1; 2—3 h, 1.1d1; 3—3 h, 1.2d1, 4—0 h。

    在Case Ⅰ,Case Ⅱ和Case Ⅲ下,冰晶平均粒徑隨存儲時間的變化見圖5。從圖5可知:冰晶平均粒徑隨時間逐漸增大,但增長速率隨時間逐漸減小。這主要是因為在冰漿儲存初期,蓄冰槽內冰晶粒徑較小,此時,冰晶粒子的聚集速率較大,破碎速率較小,所以,冰晶的粒徑迅速增大;而存儲時間越長,冰晶粒徑在逐漸增大的同時其破碎速率迅速增大,導致冰晶粒子聚集作用與破碎作用之差逐漸減小,因此,冰晶粒徑會增大,但增大速率逐漸降低。經對比可知:槳葉直徑越大,冰晶的平均粒徑越小,驗證了增大槳葉直徑能夠抑制冰晶粒徑增大,且槳葉直徑越大,其抑制作用越強。

    1—d1; 2—1.1d1; 3—1.2d1。

    3.1.2 槳葉位置對冰漿存儲的的影響

    在Case Ⅳ下,當蓄冰槽內存儲時間為1.5 h和 3.0 h時冰晶體積分數(shù)分布見圖6。對比圖6與圖3可知:槳葉距蓄冰槽底部高度越大,蓄冰槽上部冰晶粒子體積分數(shù)越小,底部冰晶體積分數(shù)越大,冰晶富集現(xiàn)象越輕微。這是由于擾動區(qū)域越接近冰晶堆積區(qū)域,更能夠削弱富冰層的形成,且由于冰晶流動在槳葉上部均會形成旋渦,越大,旋渦越大,旋渦處冰晶流速越大,冰晶越不容易向上流動,故螺旋槳周圍冰晶粒子體積分數(shù)均呈現(xiàn)上少下多的分布情況;隨增大,其槳葉上下濃度差值變小且冰晶富集分布區(qū)域變小。

    在CaseⅠ和Case Ⅳ下,冰晶粒徑數(shù)量密度分布見圖7。從圖7可知:冰晶數(shù)量密度分布曲線均隨時間向右推移,出現(xiàn)由“高窄”變?yōu)椤鞍帧钡默F(xiàn)象,以=3 h為例,冰晶粒徑數(shù)量密度分布曲線幾乎重合。這是由于當變化、槳葉直徑不變只會引起擾動區(qū)域變化,而擾動強度不變,故冰晶動力學事件的發(fā)生速率基本沒有變化??梢姡簶~位置的改變主要引起冰晶粒徑在空間分布上不同,而對冰晶數(shù)量密度分布影響很小。

    在CaseⅠ和Case Ⅳ下,冰晶平均粒徑隨存儲時間的變化見圖8。從圖8可見:槳葉高度不同,冰晶平均粒徑均隨時間而逐漸增大,但平均粒徑隨時間的變化曲線幾乎重疊,故的變化對冰晶平均粒徑的影響很小。

    (a) Case Ⅳ 1.5 h;(b) Case Ⅳ 3 h

    1—0 h; 2—3 h, 槳葉高度為35 mm; 3—3 h, 槳葉高度為65 mm。

    槳葉高度/mm:1—35;2—65。

    3.1.3 雙螺旋槳對冰漿存儲影響分析

    按螺旋槳槳葉間距不同,設置CaseⅤ,CaseⅥ和CaseⅦ這3種槳葉參數(shù),見表3。

    表3 雙螺旋槳參數(shù)

    3.1.4 槳葉間距對冰漿存儲的影響

    在Case Ⅴ,Case Ⅵ和Case Ⅶ下,存儲時間分別為1.5 h和3.0 h時蓄冰槽內冰晶體積分數(shù)分布見圖9。由圖9可知:在雙螺旋漿攪拌下,冰晶體積分數(shù)在蓄冰槽頂部和槳葉下邊緣比較大,在蓄冰槽底部和槳葉上邊緣比較小;下槳葉下邊緣更容易產生冰晶堆積,且下槳葉上、下邊緣冰晶體積分數(shù)差值較上槳葉大。由此可以看出:螺旋槳旋轉在一定程度上抑制了冰晶向上部堆積,下槳葉會阻止一部分冰晶在上槳葉下邊緣聚集;在同一時刻,槳葉間距1不同,蓄冰槽內冰晶體積分數(shù)分布大體一致,蓄冰槽頂部和底部冰晶體積分數(shù)之差先隨槳葉間距增大而減小,隨著槳葉間距繼續(xù)增大,其值反而增大。由于槳葉間距越大,對蓄冰槽上、下部冰漿的攪拌作用越強,但對蓄冰槽中部攪拌作用減弱,導致蓄冰槽內冰晶體積分數(shù)分布均勻性先隨槳葉間距增大而越來越好,當槳葉間距達到某一值后冰晶體積分數(shù)分布均勻性反而越來越差。因此,存在1個最佳的槳葉間距1使蓄冰槽內冰晶分布 均勻。

    (a) Case Ⅴ,1.5 h;(b) Case Ⅴ,3.0 h;(c) Case Ⅵ,1.5 h;(d) Case Ⅵ,3.0 h;(e) Case Ⅶ,1.5 h;(f) Case Ⅶ,3.0 h

    Case Ⅴ,Case Ⅵ和Case Ⅶ下冰晶粒徑數(shù)量密度分布見圖10。從圖10可知:冰晶數(shù)量密度分布曲線均類似于對數(shù)正態(tài)分布,隨時間向右推移,出現(xiàn)由“高窄”變?yōu)椤鞍珜挕钡默F(xiàn)象。與=3 h時對比發(fā)現(xiàn)不同槳葉間距下冰晶粒徑數(shù)量密度分布曲線幾乎重疊。這是由于槳葉間距1的變化主要引起攪拌區(qū)域發(fā)生變化,影響冰晶在蓄冰槽內的空間分布,而對擾動強度、冰晶粒子總的聚集作用和破碎作用幾乎沒有影響。

    在Case Ⅴ,Case Ⅵ和 Case Ⅶ下冰晶平均粒徑隨存儲時間的變化見圖11。從圖11可見:冰晶平均粒徑均隨時間逐漸增大,其增長速率隨時間逐漸減??;槳葉間距1不同,冰晶平均粒徑隨時間變化曲線幾乎重合,表明1對于冰晶平均粒徑演化幾乎沒有影響。

    1—3 h, 槳葉高度為30 mm; 2—3 h, 槳葉高度為50 mm; 3—3 h, 槳葉高度為70 mm; 4—0 h。

    槳葉高度/mm:1—30;2—50;3—70。

    3.2 單、雙螺旋槳對冰漿存儲影響的對比分析

    選取槳葉直徑相同、冰晶分布均勻性最好的 Case Ⅳ與Case Ⅵ,對比圖6中 Case Ⅳ與圖9中 Case Ⅵ的冰晶體積分數(shù)分布發(fā)現(xiàn):由于雙螺旋槳增大了擾動強度和擾動區(qū)域,對冰晶粒子的破碎作用增強,使得雙螺旋槳蓄冰槽內冰晶數(shù)量更大,但每個冰晶粒子的體積會比較小,雙螺旋槳蓄冰槽內冰晶總體積比較小,有效地抑制冰晶因受到升力等作用向上運動形成富冰層的趨勢,因此,雙槳葉螺旋能夠使蓄冰槽內冰晶分布更加均勻。

    在Case Ⅳ和Case Ⅵ下,冰晶粒徑數(shù)量密度分布見圖12。從圖12可知:冰晶數(shù)量密度分布曲線均隨時間向右推移,經過相同時間其推移量不同。與單螺旋槳相比,由于雙螺旋槳有效地增大了擾動強度和擾動區(qū)域,雙槳葉螺旋槳蓄冰槽內小粒徑冰晶數(shù)量密度隨時間的減小量和大粒徑冰晶粒徑數(shù)量密度增加量均較小,同一時刻冰晶粒徑數(shù)量密度峰值較大,峰值對應的冰晶粒徑比較小。因此,增加槳葉個數(shù)能夠起到抑制冰晶生長的作用。

    1—0 h;2—3 h,單螺旋槳;3—3 h,雙螺旋槳。

    單、雙螺旋槳蓄冰槽內冰晶平均粒徑隨存儲時間的變化見圖13。從圖13可見:單、雙螺旋槳蓄冰槽內冰晶平均粒徑均隨時間逐漸增大;與單螺旋槳相比,雙螺旋槳蓄冰槽中冰晶平均粒徑及其增長速率均較小。其原因是雙槳葉螺旋槳對冰漿的擾動強度及區(qū)域均比單槳葉螺旋槳的大,因此,對冰晶粒徑增長的抑制作用越強,冰晶越不容易長大。

    1—單螺旋槳;2—雙螺旋槳。

    4 結論

    1) 單螺旋槳槳葉直徑對蓄冰槽內冰晶體積分數(shù)、冰晶粒徑數(shù)量密度分布、冰晶平均粒徑均有影響。槳葉直徑越大,其擾動作用越強,對冰晶的破碎作用越明顯,蓄冰槽頂部越不容易形成富冰層。

    2) 單螺旋槳槳葉位置對冰晶流動、冰晶體積分數(shù)分布影響較大,槳葉距蓄冰槽底部的高度越大,冰晶越不容易在蓄冰槽上部聚集形成富冰層。但槳葉位置對冰晶數(shù)量密度分布、冰晶平均粒徑影響很小。

    3) 雙螺旋槳槳葉間距影響蓄冰槽內冰晶分布均勻性,并非槳葉間距越大均勻性越好,而是存在1個最佳值。但槳葉間距對于冰晶平均粒徑影響很小。

    4) 單、雙螺旋槳槳葉直徑相同,雙螺旋槳能更好地抑制冰晶粒徑增長,對冰晶生長的抑制作用更明顯。

    [1] EGOLF P W, KAUFFELD M. From physical properties of ice slurries to industrial ice slurry applications[J]. International Journal of Refrigeration, 2005, 28(1): 4?12.

    [2] KUMANO H, HIRATA T, HAGIWARA Y, et al. Effects of storage on flow and heat transfer characteristics of ice slurry[J]. International Journal of Refrigeration, 2012, 35(1): 122?129.

    [3] MICHAEL K, MASAHIRO K. Handbook on ice slurries: fundamentals and engineering[M]. Paris: International Institute of Refrigeration Press, 2005: 11?30.

    [4] MIKA L. Ice slurry flow in a poppet-type flow control value[J]. Experimental Thermal and Fluid Science, 2013(45): 128?135.

    [5] MASAYUKI T, KOZAWA Y, HIJIKATA K, et al. Prediction of ice storage process in dynamic-type ice storage system[J]. Transactions of the Japan Society of Refrigerating and Air Conditioning Engineers, 2011, 14(1): 87?95.

    [6] 青春耀, 肖睿, 宋文吉, 等. 冰漿在蓄冰槽內的蓄冰特性及其均勻度研究[J]. 低溫與超導, 2009, 37(5): 41?46. QING Chunyao, XIAO Rui, SONG Wenji, et al. Study on ice slurry storing characteristics and its uniformity in ice storage tank[J]. Cryogenics & Superconductivity, 2009, 37(5): 41?46.

    [7] KANETOSHI H, KASZA K. A method for measuring ice slurry particle agglomeration in storage tanks[J]. Dallas in ASHRAE Transactions, 2000(106): 117?123.

    [8] PRONK P, HANSEN T M, INFANTE FERREIRA C A, et al. Time-dependent behavior of different ice slurries during storage[J]. International Journal of Refrigeration, 2005, 28(1): 27?36.

    [9] DELAHAYE A, FOURNAISON L, GUILPART J. Characterisation of ice and THF hydrate slurry crystal size distribution by microscopic observation method[J]. International Journal of Refrigeration, 2010, 33(8): 1639?1647.

    [10] PRONK P, FERREIRA C A I, WITKAMP G J. Effect of long-term ice slurry storage on crystal size distribution[C]// 5thWorkshop on Ice Slurries of the IIR. Stockholm, Sweden, 2002: 30?31.

    [11] HANSEN T M, RADOSEVIC M, KAUFFELD M, et al. Investigation of ice crystal growth and geometrical characteristics in ice slurry[J]. International Journal of HVAC and R Research, 2003, 9(1): 19?32.

    [12] 趙騰磊. 冰漿存儲過程冰晶動力學演化研究[D]. 長沙: 中南大學能源科學與工程學院, 2013: 21?30.ZHAO Tenglei. Research on the ice crystal dynamics evolution in ice slurry storage process[D]. Changsha: Central South University. School of Energy Science and Engineering, 2013: 21?30.

    [13] XU Aixiang, LIU Zhiqiang, ZHAO Tenglei, et al. Population balance model of ice crystals size distribution during ice slurry storage[J]. International Journal of Air-Conditioning and Refrigeration, 2013, 22(2): 1440001-1?1440001-6.

    [14] 趙海波, 鄭楚光. 離散系統(tǒng)動力學演變過程的顆粒群平衡模擬[M]. 北京: 科學出版社, 2008: 13?18. ZHAO Haibo, ZHENG Chuguang. Population balance modeling for the process of the dynamic evolution in dispersed systems[M]. Beijing: Science Press, 2008: 13?18.

    [15] JONES G A. Crystallization process systems[M]. London: University College London, 2002: 34?45.

    [16] PRONK P. Fluidized bed heat exchangers to prevent fouling in ice slurry systems and industrial crystallizers[D]. Delft: Delft University of Technology, 2006: 22?25.

    [17] COULALOGLOU C A, TAVLARIDES L. Description of interaction processes in agitated liquid-liquid dispersions[J]. Chemical Engineering Science, 1977, 32(11): 1289?1297.

    (編輯 陳燦華)

    Effects of stirring action on ice crystal size evolution during storage

    TANG Yifang1, LIU Zhiqiang1, 2, WANG Xiaoxiao1, KANG Wei1, ZHANG Jingya1

    (1. School of Energy Science and Engineering, Central South University, Changsha 410083, China;2. Collaborative Innovation Center of Building Energy Conservation & Environmental Control, Zhuzhou, 412007, China)

    In order to study the evolution of ice crystal size with stirring action during the process of ice slurry storage, the multiphase flow model and the population balance model (PBM) were coupled to simulate the evolution of ice crystal size based on the ice crystal event nucleus models. The ice crystal event nuclear model parameters were modified based on the comparison of numerical results and experiment results. The effects of the stirring action on the evolution of ice crystal size were studied by analyzing the effect of the blade size and position of single screw, the blade spacing of twin screw on the volume fraction distribution of ice crystal, number density of ice crystal size and average particle size. The results show that the blade size has influence on volume fraction distribution of ice crystal, number density of ice crystal size and average particle size, and the blade position and spacing have significant influence on volume fraction distribution. Twin screw can make the distribution of ice crystal more uniform in the storage tank, and can be more effective to inhibit the growth of ice crystal size.

    energy storage; ice slurry; multiphase flow model; population balance model

    10.11817/j.issn.1672-7207.2017.05.033

    TB65

    A

    1672?7207(2017)05?1376?08

    2016?06?10;

    2016?08?21

    國家自然科學基金資助項目(51376198);湖南省自然科學基金資助項目(11JJ22029) (Project(51376198) supported by the National Natural Science Foundation of China; Project(11JJ22029) supported by the Natural Science Foundation of Hunan Province)

    劉志強,博士,教授,從事儲能系統(tǒng)的熱動力學研究;E-mail: liuzq@csu.edu.cn

    猜你喜歡
    槽內冰晶槳葉
    冰晶世界,凝凍微觀的美麗
    探究奇偶旋翼對雷達回波的影響
    雪為什么只在冬天下
    渠首閘后防沖槽的拋石粒徑與消能試驗研究
    立式捏合機槳葉結構與槳葉變形量的CFD仿真*
    小泥人冰晶畫 蘊藏大商機
    迷你甜甜圈制作機
    KD405:多功能折疊筷子
    冰晶奇域
    琴童(2016年5期)2016-05-14 15:43:58
    直升機槳葉/吸振器系統(tǒng)的組合共振研究
    深夜精品福利| av女优亚洲男人天堂| 国产欧美日韩一区二区三区在线| 不卡视频在线观看欧美| 欧美亚洲 丝袜 人妻 在线| 日韩三级伦理在线观看| 老司机亚洲免费影院| 亚洲欧美成人精品一区二区| 亚洲精品日韩在线中文字幕| 男女高潮啪啪啪动态图| 午夜激情av网站| 国产日韩欧美在线精品| av.在线天堂| 日韩制服骚丝袜av| 久久 成人 亚洲| 天堂俺去俺来也www色官网| 日韩欧美一区视频在线观看| 欧美精品一区二区大全| √禁漫天堂资源中文www| 久久国内精品自在自线图片| 一级爰片在线观看| 美女脱内裤让男人舔精品视频| 26uuu在线亚洲综合色| 在线观看美女被高潮喷水网站| 亚洲精品成人av观看孕妇| 成年女人在线观看亚洲视频| 男女边摸边吃奶| 国产国拍精品亚洲av在线观看| 性色avwww在线观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲国产精品一区三区| 侵犯人妻中文字幕一二三四区| 午夜福利视频在线观看免费| 99国产综合亚洲精品| 精品国产一区二区久久| 91精品国产国语对白视频| 亚洲成人一二三区av| 欧美激情 高清一区二区三区| 日韩一本色道免费dvd| 一级片'在线观看视频| 国产亚洲精品久久久com| av黄色大香蕉| 亚洲国产精品一区二区三区在线| 国产亚洲午夜精品一区二区久久| 国产黄色视频一区二区在线观看| kizo精华| 婷婷色综合大香蕉| 大片免费播放器 马上看| 在线亚洲精品国产二区图片欧美| 久久久久久久久久人人人人人人| 日日撸夜夜添| 一本大道久久a久久精品| 十八禁高潮呻吟视频| 国产成人精品在线电影| 爱豆传媒免费全集在线观看| 国产精品一区www在线观看| 欧美97在线视频| 国产免费福利视频在线观看| 一级毛片电影观看| 狂野欧美激情性bbbbbb| 欧美激情 高清一区二区三区| 汤姆久久久久久久影院中文字幕| 最近中文字幕高清免费大全6| 999精品在线视频| 99re6热这里在线精品视频| 在线观看免费日韩欧美大片| 美女xxoo啪啪120秒动态图| 欧美成人午夜精品| 日本黄色日本黄色录像| 久久精品久久久久久噜噜老黄| 视频中文字幕在线观看| 天天躁夜夜躁狠狠躁躁| 九九在线视频观看精品| 亚洲欧洲日产国产| 最近2019中文字幕mv第一页| 精品少妇黑人巨大在线播放| 午夜av观看不卡| 青春草亚洲视频在线观看| 韩国av在线不卡| 精品人妻一区二区三区麻豆| 捣出白浆h1v1| 欧美日韩视频精品一区| 国产精品人妻久久久影院| 国内精品宾馆在线| 午夜免费男女啪啪视频观看| 国产熟女午夜一区二区三区| 久久精品国产综合久久久 | 国产亚洲欧美精品永久| a级毛色黄片| 免费在线观看完整版高清| 两性夫妻黄色片 | 秋霞伦理黄片| 黄片无遮挡物在线观看| 中文字幕另类日韩欧美亚洲嫩草| 夫妻午夜视频| 黄色 视频免费看| 精品酒店卫生间| 久久久久精品久久久久真实原创| 毛片一级片免费看久久久久| 青春草视频在线免费观看| 精品人妻偷拍中文字幕| 午夜免费鲁丝| 97在线人人人人妻| av网站免费在线观看视频| 欧美少妇被猛烈插入视频| 一边摸一边做爽爽视频免费| 国产欧美亚洲国产| 久久久久久伊人网av| 国产在视频线精品| 国产成人精品一,二区| 纵有疾风起免费观看全集完整版| 国产精品女同一区二区软件| 亚洲av在线观看美女高潮| 精品国产一区二区三区久久久樱花| 成人国产av品久久久| 十分钟在线观看高清视频www| 亚洲欧美清纯卡通| 国产不卡av网站在线观看| 国产精品成人在线| 好男人视频免费观看在线| 亚洲少妇的诱惑av| 国产女主播在线喷水免费视频网站| 97在线视频观看| 青春草国产在线视频| 99re6热这里在线精品视频| 国产国语露脸激情在线看| 久久韩国三级中文字幕| 久久精品熟女亚洲av麻豆精品| 亚洲激情五月婷婷啪啪| 欧美性感艳星| 精品99又大又爽又粗少妇毛片| 欧美最新免费一区二区三区| 国产精品欧美亚洲77777| 国产乱来视频区| 99国产精品免费福利视频| 成人国产麻豆网| 中文字幕精品免费在线观看视频 | 日本av手机在线免费观看| 亚洲国产精品一区三区| 日日啪夜夜爽| 22中文网久久字幕| 久久99热这里只频精品6学生| 高清毛片免费看| 全区人妻精品视频| 亚洲国产精品一区三区| 亚洲美女搞黄在线观看| 一级片免费观看大全| 成人18禁高潮啪啪吃奶动态图| 国产成人欧美| 在线看a的网站| 搡老乐熟女国产| 亚洲精品国产av蜜桃| 国产av精品麻豆| 2018国产大陆天天弄谢| 国产乱来视频区| 女性生殖器流出的白浆| 男女边吃奶边做爰视频| 久久精品国产a三级三级三级| 成年美女黄网站色视频大全免费| 综合色丁香网| 精品亚洲成a人片在线观看| 美女大奶头黄色视频| 在线观看美女被高潮喷水网站| 丝瓜视频免费看黄片| 天堂中文最新版在线下载| 久久人人爽人人爽人人片va| 国产精品蜜桃在线观看| 永久免费av网站大全| 成人漫画全彩无遮挡| 男人舔女人的私密视频| 国产永久视频网站| av在线播放精品| 亚洲国产毛片av蜜桃av| 丝袜人妻中文字幕| 一本久久精品| 91精品三级在线观看| 日韩不卡一区二区三区视频在线| 水蜜桃什么品种好| 色网站视频免费| 欧美少妇被猛烈插入视频| 久久久久久久大尺度免费视频| 九九在线视频观看精品| 精品少妇久久久久久888优播| 日本91视频免费播放| 2022亚洲国产成人精品| 国内精品宾馆在线| 国产毛片在线视频| 一级片'在线观看视频| 亚洲av综合色区一区| 精品国产一区二区三区久久久樱花| 欧美另类一区| av在线app专区| 91午夜精品亚洲一区二区三区| 久久青草综合色| 国产一区有黄有色的免费视频| 国产黄色免费在线视频| 又黄又粗又硬又大视频| 国产精品人妻久久久久久| 成人国产麻豆网| 下体分泌物呈黄色| 夫妻性生交免费视频一级片| 免费大片18禁| 精品视频人人做人人爽| 国产一区二区在线观看av| 9热在线视频观看99| 在线天堂中文资源库| 爱豆传媒免费全集在线观看| 精品99又大又爽又粗少妇毛片| 亚洲精品aⅴ在线观看| 久久久久久久精品精品| 日韩欧美一区视频在线观看| 中国三级夫妇交换| 另类精品久久| 永久免费av网站大全| 国产av国产精品国产| 成人综合一区亚洲| 亚洲国产av新网站| 婷婷成人精品国产| 天天躁夜夜躁狠狠久久av| 91精品伊人久久大香线蕉| 在线观看免费视频网站a站| 看非洲黑人一级黄片| 中文乱码字字幕精品一区二区三区| 亚洲在久久综合| 爱豆传媒免费全集在线观看| 日韩制服丝袜自拍偷拍| av播播在线观看一区| 国产精品久久久久成人av| 亚洲久久久国产精品| 黑人高潮一二区| 另类精品久久| 日本黄大片高清| 久久青草综合色| 美女中出高潮动态图| 69精品国产乱码久久久| 亚洲欧美日韩另类电影网站| 九九在线视频观看精品| 黄片播放在线免费| 秋霞伦理黄片| 精品久久蜜臀av无| 亚洲av在线观看美女高潮| 18在线观看网站| 国语对白做爰xxxⅹ性视频网站| 色5月婷婷丁香| 亚洲精品色激情综合| 欧美亚洲日本最大视频资源| 亚洲av中文av极速乱| 久久精品国产综合久久久 | 精品国产国语对白av| 国产又色又爽无遮挡免| 日韩av免费高清视频| 在线亚洲精品国产二区图片欧美| 亚洲精品,欧美精品| av不卡在线播放| 国产一区二区三区综合在线观看 | 看免费成人av毛片| 中国国产av一级| 免费黄频网站在线观看国产| 亚洲欧洲国产日韩| 国产成人欧美| 国产黄频视频在线观看| 亚洲四区av| 国产精品.久久久| 青春草亚洲视频在线观看| 亚洲国产毛片av蜜桃av| 美女国产视频在线观看| 亚洲美女搞黄在线观看| 美女主播在线视频| 三级国产精品片| 欧美日韩综合久久久久久| 国产成人av激情在线播放| 夜夜爽夜夜爽视频| 国产成人免费观看mmmm| av线在线观看网站| 久久99精品国语久久久| 免费黄色在线免费观看| 久久久久久久久久久久大奶| av不卡在线播放| 日韩熟女老妇一区二区性免费视频| 搡老乐熟女国产| 国产色婷婷99| 精品国产一区二区久久| 国产一级毛片在线| 日本与韩国留学比较| 亚洲国产色片| 免费高清在线观看日韩| 91精品伊人久久大香线蕉| 亚洲性久久影院| 男女边吃奶边做爰视频| 久久精品熟女亚洲av麻豆精品| 亚洲综合色惰| 亚洲美女搞黄在线观看| 国产欧美日韩综合在线一区二区| 国产成人欧美| 欧美亚洲 丝袜 人妻 在线| 久久久久网色| 制服丝袜香蕉在线| 又黄又粗又硬又大视频| 全区人妻精品视频| 午夜福利影视在线免费观看| 成人亚洲欧美一区二区av| 啦啦啦中文免费视频观看日本| 久久毛片免费看一区二区三区| 精品久久蜜臀av无| 国产高清三级在线| 欧美日韩亚洲高清精品| 香蕉国产在线看| 国产黄色视频一区二区在线观看| 欧美精品人与动牲交sv欧美| av女优亚洲男人天堂| 99视频精品全部免费 在线| 久久青草综合色| 尾随美女入室| 精品一区在线观看国产| 考比视频在线观看| 97在线视频观看| 国产免费又黄又爽又色| 欧美精品国产亚洲| 成人综合一区亚洲| 一级片'在线观看视频| 一个人免费看片子| 亚洲av.av天堂| 亚洲精品国产av蜜桃| 国产黄频视频在线观看| 国产亚洲一区二区精品| 国产黄频视频在线观看| 97精品久久久久久久久久精品| 男女免费视频国产| 欧美 亚洲 国产 日韩一| 黄色一级大片看看| 高清不卡的av网站| 欧美+日韩+精品| 成人影院久久| 人妻 亚洲 视频| 大码成人一级视频| 国产xxxxx性猛交| 18禁裸乳无遮挡动漫免费视频| 国产成人欧美| 中文字幕人妻熟女乱码| 久热久热在线精品观看| 一区在线观看完整版| 激情视频va一区二区三区| 国产精品 国内视频| av不卡在线播放| 国产av码专区亚洲av| 少妇猛男粗大的猛烈进出视频| 国产一区二区三区综合在线观看 | 99视频精品全部免费 在线| 丝袜美足系列| 亚洲av电影在线进入| 久久国产亚洲av麻豆专区| 久久久久国产网址| 视频中文字幕在线观看| 国产一区二区三区综合在线观看 | 国产熟女欧美一区二区| 99久久精品国产国产毛片| 国产熟女欧美一区二区| 99久国产av精品国产电影| 国产熟女欧美一区二区| 国内精品宾馆在线| 精品福利永久在线观看| a级毛片黄视频| 在线看a的网站| 在现免费观看毛片| 亚洲av欧美aⅴ国产| 久久久久网色| 国产永久视频网站| 亚洲精品成人av观看孕妇| 亚洲精品视频女| 国产成人aa在线观看| 午夜福利在线观看免费完整高清在| 人妻人人澡人人爽人人| 亚洲欧美一区二区三区黑人 | 国产xxxxx性猛交| 亚洲国产欧美在线一区| 久久婷婷青草| 久久ye,这里只有精品| 久久毛片免费看一区二区三区| 免费av中文字幕在线| 80岁老熟妇乱子伦牲交| 激情五月婷婷亚洲| 国产毛片在线视频| 久久97久久精品| 大香蕉久久网| 最后的刺客免费高清国语| 精品国产露脸久久av麻豆| av国产久精品久网站免费入址| 在线天堂中文资源库| 高清不卡的av网站| 国产高清不卡午夜福利| 亚洲欧美清纯卡通| 天天操日日干夜夜撸| 97人妻天天添夜夜摸| 久久免费观看电影| 激情视频va一区二区三区| 香蕉国产在线看| 黄色 视频免费看| 精品国产乱码久久久久久小说| 熟妇人妻不卡中文字幕| 亚洲国产精品专区欧美| 巨乳人妻的诱惑在线观看| 国产国拍精品亚洲av在线观看| 男女无遮挡免费网站观看| 男女边吃奶边做爰视频| 久久av网站| 视频中文字幕在线观看| 久久久久久久久久久免费av| 日韩免费高清中文字幕av| 欧美xxⅹ黑人| 国产黄色视频一区二区在线观看| 91国产中文字幕| 亚洲av综合色区一区| 欧美日韩精品成人综合77777| 女的被弄到高潮叫床怎么办| av片东京热男人的天堂| 久久久久久久亚洲中文字幕| 日韩av免费高清视频| av视频免费观看在线观看| 欧美精品av麻豆av| 捣出白浆h1v1| 丰满少妇做爰视频| 亚洲国产精品一区二区三区在线| 少妇人妻精品综合一区二区| 国产一区二区三区av在线| 三上悠亚av全集在线观看| 精品久久久久久电影网| 日本欧美视频一区| 国产高清国产精品国产三级| 久久鲁丝午夜福利片| 久久 成人 亚洲| 人体艺术视频欧美日本| av女优亚洲男人天堂| 亚洲av国产av综合av卡| 啦啦啦啦在线视频资源| 亚洲欧洲国产日韩| 天天影视国产精品| 成人影院久久| 国产成人精品福利久久| 春色校园在线视频观看| 久久久亚洲精品成人影院| 国产 精品1| 免费高清在线观看视频在线观看| 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品专区欧美| 精品一品国产午夜福利视频| 久久综合国产亚洲精品| 亚洲av电影在线进入| 免费人妻精品一区二区三区视频| 国产乱来视频区| 婷婷色综合大香蕉| 狠狠精品人妻久久久久久综合| 亚洲国产av新网站| 精品一区二区三卡| 成人18禁高潮啪啪吃奶动态图| 亚洲国产日韩一区二区| 美国免费a级毛片| 人成视频在线观看免费观看| 麻豆精品久久久久久蜜桃| 免费观看a级毛片全部| 亚洲图色成人| 高清av免费在线| 尾随美女入室| 18禁在线无遮挡免费观看视频| 欧美xxⅹ黑人| 久久99热6这里只有精品| 久久久久精品久久久久真实原创| 又大又黄又爽视频免费| 精品人妻一区二区三区麻豆| 一级a做视频免费观看| 免费在线观看完整版高清| 国产成人免费无遮挡视频| 你懂的网址亚洲精品在线观看| 久久久久精品性色| 国产精品成人在线| 国产乱人偷精品视频| 啦啦啦在线观看免费高清www| 国产午夜精品一二区理论片| 看免费av毛片| 日本免费在线观看一区| 一边摸一边做爽爽视频免费| 少妇人妻 视频| 午夜影院在线不卡| 欧美性感艳星| 色网站视频免费| 国产成人aa在线观看| 国产熟女午夜一区二区三区| 亚洲国产毛片av蜜桃av| 亚洲国产精品一区三区| 夫妻性生交免费视频一级片| 亚洲欧美一区二区三区国产| 精品国产乱码久久久久久小说| 少妇精品久久久久久久| 欧美精品av麻豆av| 成人亚洲精品一区在线观看| 亚洲国产精品成人久久小说| 一区在线观看完整版| 午夜福利视频在线观看免费| 国产精品一区二区在线不卡| 国产免费现黄频在线看| 曰老女人黄片| 99香蕉大伊视频| 国产一区有黄有色的免费视频| 国产高清国产精品国产三级| 人人妻人人澡人人爽人人夜夜| 色婷婷久久久亚洲欧美| 亚洲熟女精品中文字幕| 精品亚洲成国产av| 亚洲天堂av无毛| 日日啪夜夜爽| 成人影院久久| 黑丝袜美女国产一区| 国产成人aa在线观看| 久久精品国产自在天天线| 国产免费又黄又爽又色| xxxhd国产人妻xxx| 日本-黄色视频高清免费观看| 国产精品一二三区在线看| www.色视频.com| 国产午夜精品一二区理论片| 毛片一级片免费看久久久久| a级毛片在线看网站| 女人被躁到高潮嗷嗷叫费观| 啦啦啦在线观看免费高清www| 亚洲成av片中文字幕在线观看 | 亚洲五月色婷婷综合| 免费在线观看黄色视频的| 精品一区二区三卡| 91在线精品国自产拍蜜月| 大话2 男鬼变身卡| 飞空精品影院首页| 久久人人爽人人爽人人片va| 九草在线视频观看| 欧美日韩av久久| 国产熟女午夜一区二区三区| 日韩视频在线欧美| 亚洲色图综合在线观看| 国产免费一级a男人的天堂| 亚洲国产欧美日韩在线播放| 777米奇影视久久| 久久久久精品人妻al黑| 中文乱码字字幕精品一区二区三区| 一区在线观看完整版| 久久国产精品男人的天堂亚洲 | 精品熟女少妇av免费看| 久久久精品区二区三区| 色吧在线观看| 色5月婷婷丁香| 高清欧美精品videossex| 中文欧美无线码| 乱码一卡2卡4卡精品| 9色porny在线观看| 精品熟女少妇av免费看| 欧美亚洲日本最大视频资源| 又粗又硬又长又爽又黄的视频| 亚洲精品美女久久av网站| a级毛片黄视频| 久久久久精品久久久久真实原创| 人妻 亚洲 视频| 精品卡一卡二卡四卡免费| 这个男人来自地球电影免费观看 | 久久久久久久大尺度免费视频| 高清av免费在线| 女性生殖器流出的白浆| 少妇 在线观看| av片东京热男人的天堂| 深夜精品福利| 久久久久网色| 人人妻人人澡人人看| 久久久欧美国产精品| 最后的刺客免费高清国语| 街头女战士在线观看网站| 又黄又爽又刺激的免费视频.| 精品少妇久久久久久888优播| 亚洲精品aⅴ在线观看| 亚洲美女搞黄在线观看| 久久午夜综合久久蜜桃| 免费观看在线日韩| 欧美日韩亚洲高清精品| 欧美另类一区| 在线天堂最新版资源| 欧美亚洲日本最大视频资源| 日本猛色少妇xxxxx猛交久久| 久久久国产精品麻豆| 美女国产高潮福利片在线看| 亚洲欧美成人综合另类久久久| 国产午夜精品一二区理论片| av电影中文网址| 国产 一区精品| 男人添女人高潮全过程视频| 精品少妇内射三级| 一区二区三区精品91| 精品人妻一区二区三区麻豆| 人人妻人人添人人爽欧美一区卜| 亚洲av日韩在线播放| 亚洲美女搞黄在线观看| 青春草亚洲视频在线观看| 大话2 男鬼变身卡| 欧美成人午夜免费资源| 欧美性感艳星| 永久免费av网站大全| 久久精品国产a三级三级三级| 亚洲av电影在线观看一区二区三区| 制服丝袜香蕉在线| 自线自在国产av| 精品第一国产精品| 男女下面插进去视频免费观看 | 黑人猛操日本美女一级片| 大陆偷拍与自拍| 久热久热在线精品观看| 久久久精品94久久精品| 你懂的网址亚洲精品在线观看| 女的被弄到高潮叫床怎么办| 欧美精品av麻豆av| 寂寞人妻少妇视频99o| 九九爱精品视频在线观看|