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

    轉(zhuǎn)盤表面黏性薄液膜穩(wěn)態(tài)流動特性數(shù)值模擬

    2017-06-05 01:21:40王東祥凌祥彭浩楊新俊崔政偉
    化工學報 2017年6期
    關(guān)鍵詞:水躍量綱徑向速度

    王東祥,凌祥,彭浩,楊新俊,崔政偉

    ?

    轉(zhuǎn)盤表面黏性薄液膜穩(wěn)態(tài)流動特性數(shù)值模擬

    王東祥1,2,凌祥2,彭浩2,楊新俊1,2,崔政偉1

    (1江南大學機械工程學院,江蘇省食品先進制造裝備技術(shù)重點實驗室,江蘇無錫214122;2南京工業(yè)大學江蘇省過程強化與新能源裝備技術(shù)重點實驗室,江蘇南京211800)

    轉(zhuǎn)盤表面薄膜流廣泛存在于轉(zhuǎn)盤反應器、造粒和分子蒸餾等化工領(lǐng)域,其流動特性對造粒、反應以及熱質(zhì)傳遞具有重要影響。建立了穩(wěn)態(tài)薄膜流動特性理論模型,并對比實驗數(shù)據(jù)與數(shù)值模擬結(jié)果,研究了等效Froude數(shù)、量綱1特征厚度與長度比值和澆注尺寸r對薄膜厚度分布的影響,導出了水躍和同步區(qū)半徑模型并通過實驗驗證。結(jié)果表明:等效Froude數(shù)不會對薄膜厚度分布產(chǎn)生明顯影響;量綱1和澆注尺寸r是水躍發(fā)生與否的決定條件,增大或縮小r有助于水躍出現(xiàn),水躍半徑穩(wěn)定于=0.85;平均徑向速度呈典型三分區(qū)特征,縮小澆注尺寸r將引起澆注區(qū)范圍加大而加速區(qū)范圍顯著縮小,導致薄膜無明顯加速現(xiàn)象而直接進入同步區(qū),但同步區(qū)半徑維持在=1.53。研究結(jié)果為轉(zhuǎn)盤反應器以及離心?;鞯鹊脑O(shè)計與優(yōu)化提供了可借鑒的理論與應用基礎(chǔ)。

    流體力學;模型;實驗驗證;流域;轉(zhuǎn)盤;穩(wěn)態(tài)流

    引 言

    液體在轉(zhuǎn)盤表面因離心作用以薄層形式向外流動形成薄膜流,具有高熱質(zhì)傳遞系數(shù)等獨特優(yōu)點,廣泛存在于分子蒸餾[1]、轉(zhuǎn)盤反應器[2-8]、離心造粒[9-11]以及高溫熔融渣料余熱回收[12-14]等領(lǐng)域。薄膜流體自由液面形狀或者液膜厚度對造粒和熱質(zhì)傳遞具有重要影響[15-18]。如可調(diào)控液膜厚度以避免傳熱表面發(fā)生干斑或干區(qū)而導致熱敏性物料變味變質(zhì)、非熱敏性物料結(jié)焦,以至覆蓋傳熱界面引起傳熱能力或者反應速率的顯著下降[19-20];在轉(zhuǎn)盤粒化法回收熔融渣料余熱工藝中,液膜厚度決定了流體在轉(zhuǎn)盤邊緣的破碎形態(tài)[21-23],進而作用于熔融渣顆粒的凝固換熱和物相轉(zhuǎn)變[24],是影響余熱回收效率和渣粒資源化再利用的重要因素。

    Emslie等[25]最早探索轉(zhuǎn)盤表面的薄膜流動行為,研究發(fā)現(xiàn)若忽略薄膜與轉(zhuǎn)盤間的速度滑移(即忽略科式力),薄膜流理論上由黏性力與離心力共同作用,較大Ekman數(shù)(>1)條件下的液膜厚度與轉(zhuǎn)盤徑向距離間滿足∝-2/3,提出了離心模型。在一些特殊條件下,某一轉(zhuǎn)盤區(qū)域的液膜厚度會顯著上升,形成水躍[26-29],但實現(xiàn)水躍半徑的理論預測還鮮有涉及。Wang等[30]在離心模型基礎(chǔ)上考慮了速度滑移的影響,研究發(fā)現(xiàn)離心模型只適用于水躍之后的薄膜流,科式力會顯著提高水躍區(qū)域的薄膜厚度。Woods[31]分析了強慣性力(<1)薄膜流,提出了Pigford模型。Burns等[32]采用電阻法對薄膜流的平均徑向速度進行測量,發(fā)現(xiàn)徑向速度分布可分為澆注區(qū)、加速區(qū)以及同步區(qū),并提出黏性因子,對Pigford模型進行了修正。Sisoev等[33]研究發(fā)現(xiàn)空氣誘導剪切力引起薄膜表面可能伴隨波動,但平均液膜厚度不會顯著改變,轉(zhuǎn)盤軸心和邊緣區(qū)域的液膜近似呈穩(wěn)態(tài)流動[34]。

    現(xiàn)有研究工作主要針對轉(zhuǎn)盤軸心或邊緣區(qū)域的薄膜流動行為。本文針對整個轉(zhuǎn)盤表面薄膜流體的穩(wěn)態(tài)流動特性,采用量綱1分析與積分平均法簡化Navier-Stokes(N-S)方程,建立穩(wěn)態(tài)薄膜流動特性理論模型,研究薄膜自由液面形狀和液膜厚度分布,同時探討水躍和同步區(qū)演化的一般規(guī)律,為轉(zhuǎn)盤反應器以及離心?;鞯鹊脑O(shè)計與優(yōu)化提供可借鑒的理論與應用基礎(chǔ)。

    1 穩(wěn)態(tài)薄膜流動模型

    1.1 控制方程量綱1化

    工質(zhì)以恒定流量沿轉(zhuǎn)盤軸心澆注,液膜呈旋轉(zhuǎn)周期性流動??紤]一無限大轉(zhuǎn)盤,水平放置,以恒定角速度沿軸旋轉(zhuǎn)。工質(zhì)澆注口半徑為in,流量為,牛頓流體,不可壓縮,完全潤濕轉(zhuǎn)盤表面。液膜厚度分布如圖1所示,r為水躍半徑。

    參考系固定于轉(zhuǎn)盤,微團所受體積力為

    式中,為工質(zhì)密度,為轉(zhuǎn)盤轉(zhuǎn)速,為液膜與轉(zhuǎn)盤之間相對速度,為重力加速度。

    將式(1)代入Navier-Stokes(N-S)方程并進行量綱1化處理化簡方程,引入薄膜特征厚度0以及特征長度0[29,34]

    式中,為工質(zhì)運動黏度。基于相似定律,速度尺度和壓力尺度分別為

    0=0=l0,0=0,0=220(3)

    結(jié)合式(1)~式(3)處理連續(xù)性方程和N-S方程可得式(4)~式(7)。、、、分別為量綱1速度和壓力,=0/0,=20/。

    (5)

    (6)

    1.2 薄膜近似

    澆注口尺寸一般滿足in?0,薄膜特征厚度0相比其特征長度0為微量,可采用薄膜近似[30],即=0/0?1,式(4)~式(7)中的2和4項均可忽略;采用積分平均法將式(4)~式(7)中的物理量按照處理。變換后,式(4)~式(7)進一步簡化為式(8)。為便于理解,式(8)中已略去所有物理量的上劃線。

    式中,=s+,為等效Froude數(shù),s為大氣壓強。

    1.3 邊界條件

    薄膜流表面需滿足運動學以及應力平衡邊界條件。Kim等[29]認為r處的切向速度可采用式(9)計算,但并未給出r的計算方法。而工質(zhì)澆注條件(澆注尺寸或澆注速度)均有可能影響薄膜流特性,可定義r=in/0。若忽略工質(zhì)沖擊轉(zhuǎn)盤引起的微弱能量損失,動能由軸向完全轉(zhuǎn)變?yōu)閺较颍?i>r處平均徑向速度可以近似為式(10),=r處的液膜厚度應滿足連續(xù)性方程式(11)。綜上,式(8)需滿足邊界條件式(9)~式(11)。

    (10)

    (11)

    1.4 模型驗證

    薄膜流厚度一般在數(shù)百微米,測量方法一般為探針法[28]和電阻法[26,32]。圖2、圖3為式(8)預測值分別與Leshev等[28]的探針法和Miyasaka[26]的電阻法測量值對比,工質(zhì)均為甘油水溶液。Miyasaka[26]的實驗條件為=1200 kg·m-3,=(3~8)×10-3Pa·s,進口in為200,=0.1~0.3??梢钥闯?,薄膜的厚度分布與實驗測量吻合得很好。Miyasaka在實驗中觀察到了水躍現(xiàn)象,從圖3可以看出,式(8)準確地預測到了水躍的發(fā)生半徑,并且薄膜厚度與測量值也完全吻合。

    Burns等[32]采用電阻法對薄膜平均徑向速度進行了測量,發(fā)現(xiàn)徑向速度分布可分為澆注區(qū)、加速區(qū)以及同步區(qū)。由圖4可以看出,式(8)對薄膜流徑向速度的預測值與測量值吻合得很好,完全可以體現(xiàn)速度的三分區(qū)特征。

    2 結(jié)果與討論

    2.1 液膜厚度分布

    由式(8)~式(11)可以看出,薄膜流的流動特性與等效Froude數(shù)、澆注口半徑r以及有關(guān)。圖5為=0.01與r=0.1時不同的薄膜厚度(或自由液面形狀)沿徑向的分布??梢钥闯?,隨著由0增至0.04,薄膜厚度自澆注口后先呈微弱下降,而后在水躍半徑處匯聚,在水躍之后的薄膜流同樣先微弱下降而后完全重合。水躍半徑r則不受控制始終維持在0.85左右,另外對于>1.5的轉(zhuǎn)盤區(qū)域,薄膜流自由液面完全重合。因此只對水躍附近的薄膜流產(chǎn)生了微弱影響,而不會改變水躍半徑以及遠離轉(zhuǎn)盤軸心的薄膜流特性,忽略后在=0.48處引入的最大誤差小于3%。從式(8)可以看出,等效Froude數(shù)對薄膜流的影響主要通過薄膜流體壓強來體現(xiàn),但由于薄膜流厚度一般在數(shù)百微米,其內(nèi)部壓強近似等于自由液面的大氣壓強而沒有顯著改變,即使Froude數(shù)在0.25~∝范圍時也不會對薄膜流自由液面形狀產(chǎn)生顯著影響。

    圖6為=0與r=0.1時不同值的薄膜厚度徑向分布。隨著由0.001增至0.2,薄膜流從無水躍現(xiàn)象逐漸過渡到明顯的發(fā)生水躍,水躍半徑最終穩(wěn)定于=0.85的轉(zhuǎn)盤位置;當<1.5時,薄膜厚度呈顯著下降直至>0.05后趨于穩(wěn)定,而>1.5轉(zhuǎn)盤區(qū)域的薄膜流則不受影響,完全重合。由于相同澆注條件下主要體現(xiàn)了轉(zhuǎn)速對薄膜流的影響,因此水躍現(xiàn)象是否發(fā)生與轉(zhuǎn)盤轉(zhuǎn)速直接相關(guān)。圖7所示為=0與=0.01時不同澆注口尺寸r的薄膜厚度分布。其變化規(guī)律與圖6類似,澆注口半徑r<0.1時水躍現(xiàn)象開始于=0.85處出現(xiàn),當r降至0.02以下時,水躍穩(wěn)定且自由液面開始趨于一致。Wang等[30]采用CFD方法對不同澆注口尺寸in時的薄膜流特性進行了數(shù)值分析,其中=45 ml·s-1,=2590 kg·m-3,=0.7 Pa·s,=0.478 N·m-1,=209.4 r·s-1。從圖8同樣可以看出,相同工況條件下澆注口尺寸顯著影響水躍是否發(fā)生和水躍前的薄膜厚度,而對水躍后區(qū)域則無影響。

    結(jié)合圖6和圖7可以看出,相同工況下,工質(zhì)進口條件(澆注口尺寸、澆注流量)以及轉(zhuǎn)盤轉(zhuǎn)速直接決定了水躍現(xiàn)象是否發(fā)生,并且影響<1.5轉(zhuǎn)盤區(qū)域的薄膜流特性。但水躍的發(fā)生半徑始終維持在=0.85左右,>1.5轉(zhuǎn)盤區(qū)域的薄膜流特性也不受影響。原因主要在于以=1.5為界,>1.5薄膜流徑向流動特性主要受離心力和黏性力控制,而<1.5的區(qū)域作用力主要為慣性力和黏性力。對于給定的澆注流量,澆注口尺寸縮小引起工質(zhì)沖擊轉(zhuǎn)盤軸心速度的顯著上升,在黏性力作用下澆注區(qū)速度急劇縮小,工質(zhì)由急流轉(zhuǎn)變?yōu)榫徚?,部分動能轉(zhuǎn)化為位能而導致水躍現(xiàn)象的發(fā)生。而增大轉(zhuǎn)盤轉(zhuǎn)速時,慣性力影響區(qū)域減小,離心作用區(qū)域加大而引起薄膜快速進入加速區(qū)而無水躍現(xiàn)象發(fā)生。

    2.2 同步區(qū)半徑

    Burns等[32]采用電阻法對薄膜流的平均徑向速度進行測量,發(fā)現(xiàn)徑向速度分布可分為澆注區(qū)、加速區(qū)以及同步區(qū)。以圖9中r=0.15為例,工質(zhì)澆注至轉(zhuǎn)盤軸心后由于黏性力作用引起速度下降,在水躍附近速度降至最低,而后離心力引起薄膜切向速度逐漸與轉(zhuǎn)盤一致,徑向速度也同時上升直至=1.53處達到最大,最終進入同步區(qū)。但從圖9也可以看出,隨著澆注口尺寸的縮小或薄膜流慣性力的增大,澆注區(qū)范圍逐漸加大,工質(zhì)速度因黏性減速效果也越明顯,而加速區(qū)范圍顯著縮小,薄膜無明顯加速現(xiàn)象,但不影響同步區(qū)尺寸。圖10所示為根據(jù)文獻測量數(shù)據(jù)得到的同步區(qū)半徑回歸方程??梢钥闯觯烤V1同步區(qū)半徑=1.47,與本文的理論分析結(jié)果基本一致,這也驗證了本文的研究結(jié)論。

    3 結(jié) 論

    本文以轉(zhuǎn)盤表面薄膜流穩(wěn)態(tài)流動特性為研究對象,采用量綱1分析與積分平均法化簡N-S方程,建立了量綱1薄膜流理論模型,研究了等效Froude數(shù)、量綱1和澆注口尺寸r對量綱1薄膜厚度分布的影響,探討了薄膜流水躍半徑以及同步區(qū)半徑與操作參數(shù)的一般規(guī)律。得到以下結(jié)論。

    (1)等效Froude數(shù)在0~0.04范圍內(nèi)時,對薄膜流的影響主要通過薄膜流體壓強來體現(xiàn),但由于薄膜流厚度一般在數(shù)百微米,其內(nèi)部壓強近似等于自由液面的大氣壓強而沒有顯著改變,不會對量綱1薄膜厚度產(chǎn)生明顯影響,進而水躍半徑與無關(guān)。

    (2)量綱1以及澆注口尺寸r直接決定了水躍現(xiàn)象是否發(fā)生。澆注口尺寸縮小引起工質(zhì)沖擊轉(zhuǎn)盤軸心速度顯著上升,黏性力引起澆注區(qū)速度由急流急劇向緩流轉(zhuǎn)變,流體位能增大而有利于水躍發(fā)生;而慣性力影響區(qū)域隨轉(zhuǎn)盤轉(zhuǎn)速增大而減小,離心作用區(qū)域加大而引起薄膜快速進入加速區(qū)而無水躍現(xiàn)象發(fā)生。增大或者縮小r均有利于水躍的發(fā)生,量綱1水躍半徑始終維持在=0.85。

    (3)>1.53區(qū)域的薄膜流徑向流動特性主要受離心力和黏性力控制,而<1.53的區(qū)域作用力主要為慣性力和黏性力,由此導致薄膜的平均徑向速度呈現(xiàn)典型三分區(qū)特征,縮小澆注口尺寸r將引起澆注區(qū)范圍加大而加速區(qū)范圍顯著縮小,導致薄膜無明顯加速現(xiàn)象而直接進入同步區(qū),但同步區(qū)半徑始終維持在=1.53,與實驗測量結(jié)果=1.47基本一致。

    符 號 說 明

    F——力,N g——重力加速度,m·s-2 h——薄液膜厚度,m l0——薄膜特征長度,m p——量綱1薄膜流壓強 ps——大氣壓強,Pa Q——工質(zhì)流量,m3·s-1 r——量綱1轉(zhuǎn)盤半徑 ri——量綱1進口半徑 rin,rj,rs——分別為工質(zhì)進口半徑、水躍半徑與同步區(qū)半徑,m u,v,w——分別為薄膜徑向、切向與軸向速度,m·s-1 u0,v0,w0——分別為量綱1薄膜徑向、切向與軸向速度 d0——薄膜特征厚度,m e——薄膜特征厚度與長度比值 μ——工質(zhì)動力黏度,Pa·s ρ——工質(zhì)密度,kg·m-3 υ——工質(zhì)運動黏度,m2·s-1 w——轉(zhuǎn)盤轉(zhuǎn)速,rad·s-1 下角標 i——澆注進口 j——水躍半徑 0——特征參數(shù)

    References

    [1] BATISTELLA C B, MORAES E B, MACIEL W M R. Comparing centrifugal and falling film molecular stills using reflux and cascade for fine chemical separations[J]. Computers & Chemical Engineering, 1999, 23: S767-S770.

    [2] BOIARKINA I, NORRIS S, PATTERSON D A. Investigation into the effect of flow structure on the photocatalytic degradation of methylene blue and dehydroabietic acid in a spinning disc reactor[J]. Chemical Engineering Journal, 2013, 222(4): 159-171.

    [3] 韓進, 朱彤, 今井剛, 等. 基于高速轉(zhuǎn)盤法的剩余污泥可溶化處理[J]. 化工學報, 2008, 59(2): 478-483. HAN J, ZHU T, JIN J G,. Solubilization of excess sludge by high speed rotary disk[J]. Journal of Chemical Industry and Engineering (China), 2008, 59(2): 478-483.

    [4] 李丹, 張大皓, 譚天偉. 轉(zhuǎn)盤反應器固定根霉重復批次發(fā)酵生產(chǎn)脂肪酶[J]. 化工學報, 2007, 58(9): 2347-2351. LI D, ZHANG D H, TAN T W. Product ion of lipase by repeated batch fermentation with immobilizedusing rotating biological contactor[J]. Journal of Chemical Industry and Engineering (China), 2007, 58(9): 2347-2351.

    [5] 任南琪, 周顯嬌, 郭婉茜, 等. 染料廢水處理技術(shù)研究進展[J]. 化工學報, 2013, 64(1): 84-94. REN N Q, ZHOU X J, GUO W Q,. A review on treatment methods of dye wastewater[J]. CIESC Journal, 2013, 64(1): 84-94.

    [6] 孫志成, 韓進, 張明楊, 等. 機械轉(zhuǎn)盤聯(lián)合超聲破解剩余污泥[J]. 化工學報, 2016, 67(12): 5229-5236. SUN Z C, HAN J, ZHANG M Y,. Disintegration of excess sludge by mechanical rotary disk combination with ultrasonication[J]. CIESC Journal, 2016, 67(12): 5229-5236.

    [7] BOIARKINA I, NORRIS S, PATTERSON D A. Investigation into the effect of flow structure on the photocatalytic degradation of methylene blue and dehydroabietic acid in a spinning disc reactor[J]. Chemical Engineering Journal, 2013, 222 (4): 159-171.

    [8] 魏冰, 張巧玲, 劉有智, 等. 旋轉(zhuǎn)盤反應器光催化降解苯酚廢水動力學[J]. 環(huán)境工程學報, 2016, 10(3): 1305-1309. WEI B, ZHANG Q L, LIU Y Z,. Kinetics of phenol by photocatalysis degradation in spinning disk reactor[J]. Chinese Journal of Environmental Engineering, 2016, 10(3): 1305-1309.

    [9] WANG D X, LING X, PENG H,. Experimental investigation of ligament formation dynamics of thin viscous liquid film at spinning disk edge[J]. Industrial & Engineering Chemistry Research, 2016, 55(34): 9267-9275.

    [10] 吳君軍, 王宏, 朱恂, 等. 轉(zhuǎn)盤離心?;薪z狀成粒特性[J]. 化工學報, 2015, 66(7): 2474-2480. WU J J, WANG H, ZHU X,. Characteristic of ligament in centrifugal granulation by spinning disc[J]. CIESC Journal, 2015, 66(7): 2474-2480.

    [11] AHMED M, YOUSSEF M S. Influence of spinning cup and disk atomizer configurations on droplet size and velocity characteristics[J]. Chemical Engineering Science, 2014, 107(14): 149-157.

    [12] ZHANG H, WANG H, ZHU X,. A review of waste heat recovery technologies towards molten slag in steel industry[J]. Applied Energy, 2013, 112(4): 956-966.

    [13] BARATI M, ESFAHANI S, UTIGARD T A. Energy recovery from high temperature slags[J]. Energy, 2011, 36(9): 5440-5449.

    [14] WANG D X, LING X, PENG H. Simulation of ligament mode breakup of molten slag by spinning disk in the dry granulation process[J]. Applied Thermal Engineering, 2015, 84: 437-447.

    [15] JACOBSEN N C, HINRICHSEN O. Micromixing efficiency of a spinning disk reactor[J]. Industrial & Engineering Chemistry Research, 2012, 51(36): 11643-11652.

    [16] AHMED M, YOUSSEF M S. Characteristics of mean droplet size produced by spinning disk atomizers[J]. Journal of Fluid Engineering, 2012, 134(7): 1-9.

    [17] BOIARKINA I, NORRIS S, PATTERSON D A. The case for the photocatalytic spinning disc reactor as a process intensification technology: comparison to an annular reactor for the degradation of methylene blue[J]. Chemical Engineering Journal, 2013, 225(3): 752-765.

    [18] CAPRARIIS B D, RITA D M, STOLLER M,. Reaction-precipitation by a spinning disc reactor: influence of hydrodynamics on nanoparticles production[J]. Chemical Engineering Science, 2012, 76(28): 73-80.

    [19] 閻維平, 葉學民, 李洪濤. 液體薄膜流的流動和傳熱特性[J]. 華北電力大學學報, 2005, 32(1): 59-65. YAN W P, YE X M, LI H T. Hydrodynamics and heat transfer of liquid thin films[J]. Journal of North China Electric Power University, 2005, 32(1): 59-65.

    [20] CHIANG C, CHANG M, LIU H,. Process intensification in the production of photocatalysts for solar hydrogen generation[J]. Industrial & Engineering Chemistry Research, 2012, 51(14): 5207-5215.

    [21] PENG H, WANG N, WANG D X,. Experimental study on critical characteristics of liquid atomization by spinning disk[J]. Industrial & Engineering Chemistry Research, 2016, 55(21): 6175-6185.

    [22] 王東祥, 凌祥, 彭浩. 轉(zhuǎn)盤離心粒化熔渣液膜流動特性數(shù)值模擬研究[J]. 南京工業(yè)大學學報(自然科學版), 2015, 37(3): 67-73. WANG D X, LING X, PENG H. Numerical simulation of film flow characteristics of molten slag on spinning disk in centrifugal atomization process[J]. Journal of Nanjing Tech University (Natural Science Edition), 2015, 37(3): 67-73.

    [23] WANG D X, LING X, PENG H,. Experimental investigation of ligament formation dynamics of thin viscous liquid film at spinning disk edge[J]. Industrial & Engineering Chemistry Research, 2016, 55(34): 9267-9275.

    [24] 何先琰, 王宏, 朱恂, 等. 鉛錫合金熔融顆粒風冷相變換熱特性實驗研究[J]. 工程熱物理學報, 2015, 36(8): 1748-1751.HE X Y, WANG H, ZHU X,. Experiment study on air-cooled phase change heat transfer characteristics of Sn-Pb alloy droplets[J]. Journal of Engineering Thermophysics, 2015, 36(8): 1748-1751.

    [25] EMSLIE A G, BONNER F T, PECK L G. Flow of a viscous liquid on a rotating disk[J]. Journal of Applied Physics, 1958, 29(5): 858-862.

    [26] MIYASAKA Y. On the flow of a viscous free boundary jet on a rotating disk[J]. Bulletin of the JSME, 1974, 40(331): 1461-1468.

    [27] THOMAS S, FAGHRI A, HANKEY W. Experimental analysis and flow visualization of a thin liquid film on a stationary and rotating disk[J]. Journal of Fluids Engineering, 1991, 113(1): 73-80.

    [28] LESHEV I, PEEV G. Film flow on a horizontal rotating disk[J]. Chemical Engineering and Processing, 2003, 42 (11): 925-929.

    [29] KIM T S, KIM M U. The flow and hydrodynamic stability of a liquid film on a rotating disc[J]. Fluid Dynamics Research, 2009, 41(3): 475-479.

    [30] WANG D X, LING X, PENG H. Theoretical analysis of free-surface film flow on the rotary granulating disk in waste heat recovery process of molten slag[J]. Applied Thermal Engineering, 2014, 63(1): 387-395.

    [31] WOODS W P. The hydrodynamics of thin liquid films flowing over a rotating disk[D]. United Kingdom: University of Newcastle, 1995.

    [32] BURNS J R, RAMSHAW C, JACHUCK R J. Measurement of liquid film thickness and the determination of spin-up radius on a rotating disc using an electrical resistance technique[J]. Chemical Engineering Science, 2003, 58 (11): 2245-2253.

    [33] SISOEV G M, MATER O K, LAWRENCE C J. Axisymmetric wave regimes in viscous liquid film flow over a spinning disk[J]. Journal of Fluid Mechanics, 2003, 495(495): 385-411.

    [34] PRIELING D, STEINER H. Unsteady thin film flow on spinning disks at large ekman numbers using an integral boundary layer method[J]. International Journal of Heat and Mass Transfer, 2013, 65(7): 10-22.

    Numerical simulation of stable flow dynamics of viscous film on spinning disk surface

    WANG Dongxiang1,2, LING Xiang2, PENG Hao2, YANG Xinjun1,2, CUI Zhengwei1

    (1Jiangsu Key Laboratory of Advanced Food Manufacturing Equipment & Technology, School of Mechanical Engineering, Jiangnan University, Wuxi 214122, Jiangsu, China;2Jiangsu Key Laboratory of Process Enhancement and New Energy Equipment Technology, Nanjing Tech University, Nanjing 211800, Jiangsu, China)

    Film flow on spinning disk surface exists extensively in chemical engineering operations, such as centrifugal graining, molecular distillation and spinning disk reactors. Flow dynamics of the film flow has a major impact on graining, reaction, heat and mass transfer rate. Theoretical model of stable film flow dynamics was established and verified by comparison of experimental results and numerical simulations. The effect of equivalent Froude number, dimensionless ratio of characteristic thickness over length,and casting sizeron dimensionless film thickness distribution was studied, which a model of dimensionless hydraulic jump and synchronized zone radius was derived and verified by experiments. The results show that equivalent Froude number has little effect on film thickness distribution and occurrence of hydraulic jump phenomenon mainly depends uponandr. Increasingor decreasingrleads to appearance of hydraulic jump with a dimensionless hydraulic jump radius always at=0.85. Mean radial velocity exhibited features of typical three zone distribution. Reducing casting size would extend injection zone and shrink acceleration zone, such that film flow goes directly to synchronized zone at dimensionless synchronized radius of 1.53 without evident acceleration. The results will provide theoretical reference for design and optimization of spinning disk reactors and centrifugal pelletizers.

    fluid mechanics; model; experimental validation; flow regimes; spinning disk; steady flow

    10.11949/j.issn.0438-1157.20170040

    TQ 021.1

    A

    0438—1157(2017)06—2321—07

    凌祥。

    王東祥(1985—),男,博士,講師。

    國家自然科學基金項目(51606086,51406078)。

    2017-01-10收到初稿,2017-02-16收到修改稿。

    2017-01-10.

    Prof. LING Xiang, xling@njtech.edu.cn

    supported by the National Natural Science Foundation of China (51606086, 51406078).

    猜你喜歡
    水躍量綱徑向速度
    那吉航運樞紐大壩消力池損毀成因和搶修方案分析
    量綱分析在熱力學統(tǒng)計物理中的教學應用*
    物理通報(2022年6期)2022-05-27 09:43:14
    淺談量綱法推導物理公式的優(yōu)勢
    ——以勻加速直線運動公式為例
    突然擴散水躍方程的改進與比較
    漸擴式消力池水躍長度的分析與計算
    科技論文編輯加工中的量綱問題
    非圓形光纖研究進展
    航空兵器(2017年6期)2018-01-24 15:00:10
    臺風威馬遜造成云南文山州強降水天氣雷達回波分析
    距離頻率ML方法無模糊估計動目標徑向速度
    矩形平底明渠水躍長度公式的分析與應用
    亚洲人成网站在线播放欧美日韩| 久久人人97超碰香蕉20202| 日韩 欧美 亚洲 中文字幕| 免费在线观看完整版高清| 欧美成人性av电影在线观看| 在线av久久热| 中文字幕高清在线视频| 91字幕亚洲| 亚洲国产精品sss在线观看| 美女高潮到喷水免费观看| 亚洲五月婷婷丁香| 一卡2卡三卡四卡精品乱码亚洲| 色综合婷婷激情| 亚洲一区二区三区色噜噜| 麻豆久久精品国产亚洲av| av欧美777| 久久精品国产亚洲av香蕉五月| 91成年电影在线观看| 国产欧美日韩一区二区精品| 可以在线观看毛片的网站| 动漫黄色视频在线观看| 日韩欧美三级三区| 久久性视频一级片| 人人澡人人妻人| 国产乱人伦免费视频| 无遮挡黄片免费观看| av有码第一页| 97超级碰碰碰精品色视频在线观看| 日日爽夜夜爽网站| 亚洲av成人不卡在线观看播放网| 久久中文字幕一级| 99久久国产精品久久久| 咕卡用的链子| 搡老岳熟女国产| 亚洲性夜色夜夜综合| 人人妻人人澡人人看| 久久青草综合色| 少妇 在线观看| 美女大奶头视频| 亚洲av电影在线进入| 午夜精品久久久久久毛片777| 精品第一国产精品| 亚洲成人免费电影在线观看| 国产一卡二卡三卡精品| 久久久久九九精品影院| 成人18禁在线播放| av免费在线观看网站| 波多野结衣av一区二区av| 色综合站精品国产| 97超级碰碰碰精品色视频在线观看| 欧美一级a爱片免费观看看 | 亚洲中文av在线| 亚洲精华国产精华精| 久久久久久亚洲精品国产蜜桃av| 国产成+人综合+亚洲专区| 日韩免费av在线播放| 国产成年人精品一区二区| 母亲3免费完整高清在线观看| 热99re8久久精品国产| 午夜免费鲁丝| 亚洲熟女毛片儿| 日本 av在线| 两个人视频免费观看高清| 99久久精品国产亚洲精品| 亚洲人成77777在线视频| 亚洲色图 男人天堂 中文字幕| 午夜成年电影在线免费观看| 亚洲专区国产一区二区| 免费少妇av软件| 黑丝袜美女国产一区| 国产极品粉嫩免费观看在线| 久久精品国产综合久久久| 女人被躁到高潮嗷嗷叫费观| 一区福利在线观看| 免费人成视频x8x8入口观看| 亚洲天堂国产精品一区在线| 国产午夜福利久久久久久| 国产精品久久久人人做人人爽| 国内毛片毛片毛片毛片毛片| 欧美色欧美亚洲另类二区 | 老汉色av国产亚洲站长工具| 久99久视频精品免费| 国产精品一区二区免费欧美| 亚洲欧美日韩另类电影网站| 亚洲国产看品久久| 亚洲欧美一区二区三区黑人| 亚洲一区高清亚洲精品| 老司机深夜福利视频在线观看| 妹子高潮喷水视频| 88av欧美| 1024视频免费在线观看| 夜夜躁狠狠躁天天躁| 成年版毛片免费区| 在线观看免费视频网站a站| 国产一区二区三区在线臀色熟女| 日韩有码中文字幕| 久久精品国产综合久久久| 国产精品久久久久久亚洲av鲁大| 制服丝袜大香蕉在线| 亚洲免费av在线视频| 咕卡用的链子| avwww免费| 男女下面插进去视频免费观看| 成熟少妇高潮喷水视频| 欧美精品亚洲一区二区| 午夜免费观看网址| 99香蕉大伊视频| 99精品在免费线老司机午夜| 午夜视频精品福利| 女生性感内裤真人,穿戴方法视频| 亚洲人成电影免费在线| 一级毛片女人18水好多| av天堂在线播放| 黑丝袜美女国产一区| 日韩一卡2卡3卡4卡2021年| 成人精品一区二区免费| 午夜影院日韩av| 狠狠狠狠99中文字幕| 老司机午夜福利在线观看视频| 熟妇人妻久久中文字幕3abv| 国产主播在线观看一区二区| 一区二区三区精品91| 制服人妻中文乱码| 身体一侧抽搐| 黄色视频,在线免费观看| 精品久久蜜臀av无| 亚洲欧美精品综合一区二区三区| 国产精品一区二区三区四区久久 | 国产精品久久久久久人妻精品电影| 免费久久久久久久精品成人欧美视频| 一区二区三区激情视频| 成人免费观看视频高清| 亚洲avbb在线观看| 亚洲中文字幕一区二区三区有码在线看 | 国产av又大| 中文字幕久久专区| 伊人久久大香线蕉亚洲五| 精品国产乱码久久久久久男人| 精品久久久久久成人av| 久久精品亚洲精品国产色婷小说| 美女免费视频网站| 亚洲一码二码三码区别大吗| 亚洲免费av在线视频| 中文字幕精品免费在线观看视频| 丰满的人妻完整版| 亚洲人成电影观看| 精品久久久精品久久久| 午夜福利,免费看| 热99re8久久精品国产| 天天一区二区日本电影三级 | 黄片播放在线免费| 自线自在国产av| 精品免费久久久久久久清纯| 免费在线观看日本一区| 中文字幕人成人乱码亚洲影| 18禁黄网站禁片午夜丰满| 搞女人的毛片| 男女床上黄色一级片免费看| 亚洲欧美激情综合另类| 夜夜看夜夜爽夜夜摸| 成人免费观看视频高清| 午夜成年电影在线免费观看| 亚洲专区中文字幕在线| 99国产精品一区二区蜜桃av| 国产精品一区二区免费欧美| 亚洲精品一区av在线观看| 亚洲精品中文字幕一二三四区| 欧美大码av| 日韩免费av在线播放| 亚洲欧美一区二区三区黑人| 可以在线观看毛片的网站| 精品久久久久久,| 在线观看午夜福利视频| 男男h啪啪无遮挡| 97人妻精品一区二区三区麻豆 | 久久九九热精品免费| 丝袜美足系列| 免费在线观看完整版高清| 很黄的视频免费| 午夜激情av网站| 亚洲熟妇中文字幕五十中出| 法律面前人人平等表现在哪些方面| 身体一侧抽搐| 亚洲国产欧美一区二区综合| 色播在线永久视频| 宅男免费午夜| 中出人妻视频一区二区| 亚洲电影在线观看av| 久久中文字幕人妻熟女| 国产精品久久电影中文字幕| 亚洲精品中文字幕在线视频| 亚洲av熟女| 人人澡人人妻人| 亚洲自拍偷在线| 国产成人影院久久av| 看片在线看免费视频| 久久狼人影院| x7x7x7水蜜桃| 久久热在线av| 91成人精品电影| 十八禁人妻一区二区| 国产成人精品久久二区二区免费| 18禁国产床啪视频网站| 国产成人欧美在线观看| 国产精品综合久久久久久久免费 | 男男h啪啪无遮挡| av欧美777| 国产精品美女特级片免费视频播放器 | 国产蜜桃级精品一区二区三区| 波多野结衣一区麻豆| 无限看片的www在线观看| 欧美黄色片欧美黄色片| 国产成人啪精品午夜网站| 悠悠久久av| 99国产综合亚洲精品| 国产aⅴ精品一区二区三区波| 国产精品综合久久久久久久免费 | 一区二区三区国产精品乱码| 亚洲,欧美精品.| 久久久精品国产亚洲av高清涩受| 日韩精品中文字幕看吧| 久久精品人人爽人人爽视色| 欧美久久黑人一区二区| 国产亚洲精品第一综合不卡| 91国产中文字幕| 成人永久免费在线观看视频| 欧美黄色片欧美黄色片| 午夜成年电影在线免费观看| 久久久国产欧美日韩av| 国产精品美女特级片免费视频播放器 | 欧美激情 高清一区二区三区| 欧美色欧美亚洲另类二区 | 日本一区二区免费在线视频| 侵犯人妻中文字幕一二三四区| videosex国产| 两个人看的免费小视频| 精品国产一区二区三区四区第35| 给我免费播放毛片高清在线观看| 国产精品99久久99久久久不卡| 女性生殖器流出的白浆| 亚洲天堂国产精品一区在线| 午夜福利欧美成人| 亚洲成a人片在线一区二区| 午夜两性在线视频| 亚洲九九香蕉| 欧美一级a爱片免费观看看 | 亚洲国产精品sss在线观看| 两性夫妻黄色片| 超碰成人久久| 免费高清在线观看日韩| 黄色毛片三级朝国网站| 宅男免费午夜| 亚洲第一电影网av| 国产精品久久久av美女十八| 97超级碰碰碰精品色视频在线观看| 成人特级黄色片久久久久久久| 免费女性裸体啪啪无遮挡网站| 久久青草综合色| 国产精品香港三级国产av潘金莲| 99久久精品国产亚洲精品| 不卡一级毛片| www.www免费av| 国产成+人综合+亚洲专区| 99国产精品免费福利视频| 大码成人一级视频| 村上凉子中文字幕在线| 国产麻豆成人av免费视频| 亚洲专区国产一区二区| 免费人成视频x8x8入口观看| 又黄又粗又硬又大视频| 女性生殖器流出的白浆| 中文字幕久久专区| 美女高潮喷水抽搐中文字幕| 中出人妻视频一区二区| 涩涩av久久男人的天堂| 亚洲人成电影免费在线| 麻豆久久精品国产亚洲av| 欧美在线一区亚洲| 亚洲一区二区三区不卡视频| 亚洲精品国产精品久久久不卡| 中文字幕最新亚洲高清| 一a级毛片在线观看| 国产av又大| 97人妻天天添夜夜摸| 久久人妻福利社区极品人妻图片| 国产成人免费无遮挡视频| 丰满的人妻完整版| www.精华液| 免费看a级黄色片| 好男人在线观看高清免费视频 | 亚洲精品久久成人aⅴ小说| 欧美色欧美亚洲另类二区 | 久久久久久亚洲精品国产蜜桃av| www.自偷自拍.com| 久久午夜亚洲精品久久| 国产精品1区2区在线观看.| 精品久久久久久久毛片微露脸| 国产在线精品亚洲第一网站| 日韩精品免费视频一区二区三区| 免费在线观看亚洲国产| 纯流量卡能插随身wifi吗| 国产成人影院久久av| 757午夜福利合集在线观看| 成人亚洲精品一区在线观看| 人人澡人人妻人| 国产免费av片在线观看野外av| 久久精品国产99精品国产亚洲性色 | 亚洲激情在线av| 国产单亲对白刺激| 午夜福利18| 国产麻豆69| 国产成人欧美| 午夜精品在线福利| 亚洲一区二区三区色噜噜| 久久天躁狠狠躁夜夜2o2o| 欧美 亚洲 国产 日韩一| 久久久久久大精品| 亚洲色图综合在线观看| 亚洲精品粉嫩美女一区| 欧美一级a爱片免费观看看 | 久久精品aⅴ一区二区三区四区| 欧美黄色淫秽网站| 成人国产一区最新在线观看| 午夜福利在线观看吧| 亚洲 欧美一区二区三区| 午夜免费成人在线视频| 一本久久中文字幕| 可以在线观看毛片的网站| 国产精品一区二区在线不卡| 在线观看舔阴道视频| 国产又爽黄色视频| 夜夜爽天天搞| 国产精品电影一区二区三区| 神马国产精品三级电影在线观看 | 午夜福利影视在线免费观看| 久久精品国产清高在天天线| 欧美日韩福利视频一区二区| 老司机深夜福利视频在线观看| 免费在线观看完整版高清| 亚洲七黄色美女视频| 成人手机av| 伦理电影免费视频| 日本免费一区二区三区高清不卡 | 午夜精品在线福利| 亚洲自偷自拍图片 自拍| 亚洲精品一区av在线观看| 欧美日韩亚洲综合一区二区三区_| 自线自在国产av| av天堂在线播放| 男女下面插进去视频免费观看| www.熟女人妻精品国产| 成人av一区二区三区在线看| 黄色视频不卡| 精品日产1卡2卡| 亚洲成人免费电影在线观看| 亚洲av日韩精品久久久久久密| 久久草成人影院| 亚洲精品中文字幕在线视频| 国内毛片毛片毛片毛片毛片| 久久影院123| 自线自在国产av| 色播亚洲综合网| 亚洲色图综合在线观看| 国产xxxxx性猛交| 给我免费播放毛片高清在线观看| 精品国产国语对白av| 在线永久观看黄色视频| 久99久视频精品免费| 香蕉久久夜色| 国产av精品麻豆| 亚洲av成人一区二区三| 99国产综合亚洲精品| 欧美乱码精品一区二区三区| 国产精品精品国产色婷婷| 精品高清国产在线一区| 天天添夜夜摸| 熟妇人妻久久中文字幕3abv| 国内毛片毛片毛片毛片毛片| 欧美绝顶高潮抽搐喷水| 国产成人精品久久二区二区91| 无限看片的www在线观看| 国产伦人伦偷精品视频| 亚洲一区二区三区不卡视频| 免费女性裸体啪啪无遮挡网站| 亚洲一区二区三区不卡视频| 亚洲av熟女| 亚洲国产精品合色在线| 1024视频免费在线观看| 亚洲欧美激情在线| 如日韩欧美国产精品一区二区三区| 亚洲中文字幕日韩| 国产色视频综合| 国产一区二区三区视频了| 久久九九热精品免费| 亚洲精品中文字幕在线视频| 在线观看免费午夜福利视频| 伊人久久大香线蕉亚洲五| 极品人妻少妇av视频| 激情视频va一区二区三区| 亚洲一码二码三码区别大吗| 可以在线观看毛片的网站| 亚洲国产看品久久| bbb黄色大片| 97人妻天天添夜夜摸| 91成年电影在线观看| 制服人妻中文乱码| 97人妻精品一区二区三区麻豆 | 国产午夜精品久久久久久| 色av中文字幕| 自线自在国产av| a在线观看视频网站| 欧美绝顶高潮抽搐喷水| 亚洲精品美女久久av网站| 亚洲精品国产一区二区精华液| 精品电影一区二区在线| 色播亚洲综合网| 黑人操中国人逼视频| 国产精品一区二区三区四区久久 | 大码成人一级视频| 欧美黑人精品巨大| 91大片在线观看| 深夜精品福利| 久久青草综合色| 视频区欧美日本亚洲| 波多野结衣巨乳人妻| 久久久久久国产a免费观看| 一级毛片高清免费大全| 制服丝袜大香蕉在线| www日本在线高清视频| 欧美成人性av电影在线观看| 一级a爱视频在线免费观看| 91国产中文字幕| 丝袜在线中文字幕| 精品无人区乱码1区二区| 精品人妻1区二区| 成年女人毛片免费观看观看9| 美女大奶头视频| 精品久久久久久成人av| 亚洲男人天堂网一区| 国产成人av激情在线播放| 啦啦啦 在线观看视频| 中文字幕高清在线视频| 狠狠狠狠99中文字幕| 国产精品99久久99久久久不卡| 国产成人影院久久av| 一区二区三区高清视频在线| 极品人妻少妇av视频| 亚洲国产精品999在线| √禁漫天堂资源中文www| 国产精品久久久久久精品电影 | 亚洲精品一卡2卡三卡4卡5卡| 日韩高清综合在线| 国产精品综合久久久久久久免费 | 不卡一级毛片| 波多野结衣一区麻豆| 亚洲一区二区三区色噜噜| 操出白浆在线播放| 久久久国产欧美日韩av| 国产av一区在线观看免费| 久久精品国产清高在天天线| 亚洲色图综合在线观看| 叶爱在线成人免费视频播放| 老司机在亚洲福利影院| 两性夫妻黄色片| 国产精品久久久久久亚洲av鲁大| 久久久国产成人精品二区| av欧美777| 成人18禁高潮啪啪吃奶动态图| 日韩三级视频一区二区三区| 99久久99久久久精品蜜桃| 性色av乱码一区二区三区2| 国产精品1区2区在线观看.| 成人手机av| 神马国产精品三级电影在线观看 | 波多野结衣一区麻豆| 婷婷精品国产亚洲av在线| 国产亚洲欧美精品永久| 色综合亚洲欧美另类图片| 欧美久久黑人一区二区| 90打野战视频偷拍视频| 欧美日韩瑟瑟在线播放| 久久天躁狠狠躁夜夜2o2o| 999久久久精品免费观看国产| 欧美日韩福利视频一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 高清在线国产一区| 亚洲精品粉嫩美女一区| 欧美精品啪啪一区二区三区| 一本久久中文字幕| www.自偷自拍.com| 一卡2卡三卡四卡精品乱码亚洲| 18禁国产床啪视频网站| 看片在线看免费视频| 精品免费久久久久久久清纯| 国产精品1区2区在线观看.| 后天国语完整版免费观看| 99国产极品粉嫩在线观看| 大码成人一级视频| 亚洲第一青青草原| 久久国产亚洲av麻豆专区| 日韩欧美免费精品| 午夜激情av网站| 99精品久久久久人妻精品| 男女午夜视频在线观看| 在线播放国产精品三级| av有码第一页| 岛国在线观看网站| 午夜福利视频1000在线观看 | e午夜精品久久久久久久| av视频免费观看在线观看| 亚洲 欧美一区二区三区| 免费女性裸体啪啪无遮挡网站| 侵犯人妻中文字幕一二三四区| 日本一区二区免费在线视频| 日韩 欧美 亚洲 中文字幕| 国内精品久久久久精免费| 日韩视频一区二区在线观看| av免费在线观看网站| 一级作爱视频免费观看| 亚洲国产看品久久| 精品国产乱子伦一区二区三区| 欧美日韩乱码在线| 欧美成人午夜精品| 精品一区二区三区四区五区乱码| 男女做爰动态图高潮gif福利片 | 国产91精品成人一区二区三区| 国产精品野战在线观看| 日本a在线网址| 国产亚洲精品久久久久久毛片| 欧美绝顶高潮抽搐喷水| 国产成人影院久久av| 亚洲精品在线观看二区| 色尼玛亚洲综合影院| 国产精品九九99| 国产精品国产高清国产av| 男人操女人黄网站| 97人妻精品一区二区三区麻豆 | 国产1区2区3区精品| 亚洲av第一区精品v没综合| 激情视频va一区二区三区| 操美女的视频在线观看| 国产蜜桃级精品一区二区三区| 久久 成人 亚洲| 亚洲七黄色美女视频| 每晚都被弄得嗷嗷叫到高潮| 国产亚洲欧美98| 欧美午夜高清在线| 国产亚洲精品一区二区www| 18禁裸乳无遮挡免费网站照片 | 日韩 欧美 亚洲 中文字幕| 91老司机精品| 免费女性裸体啪啪无遮挡网站| 亚洲国产高清在线一区二区三 | 日本三级黄在线观看| 国产午夜精品久久久久久| 欧美激情久久久久久爽电影 | 亚洲自拍偷在线| 久久国产精品人妻蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 亚洲 欧美 日韩 在线 免费| e午夜精品久久久久久久| 免费在线观看完整版高清| 亚洲一区二区三区不卡视频| 欧美成狂野欧美在线观看| 亚洲美女黄片视频| xxx96com| 两个人视频免费观看高清| 亚洲国产欧美日韩在线播放| 人妻久久中文字幕网| АⅤ资源中文在线天堂| 97人妻天天添夜夜摸| 三级毛片av免费| 国产91精品成人一区二区三区| 国产熟女xx| 久久性视频一级片| 国产成人系列免费观看| 可以在线观看毛片的网站| 亚洲五月天丁香| 大型黄色视频在线免费观看| 午夜久久久久精精品| 丝袜美腿诱惑在线| 精品国产超薄肉色丝袜足j| 可以免费在线观看a视频的电影网站| 国产成年人精品一区二区| 亚洲va日本ⅴa欧美va伊人久久| 国产熟女午夜一区二区三区| 中文字幕人成人乱码亚洲影| 一级毛片女人18水好多| a级毛片在线看网站| 嫩草影视91久久| 欧美大码av| 久久精品国产亚洲av高清一级| 一进一出抽搐gif免费好疼| 亚洲国产精品久久男人天堂| or卡值多少钱| 亚洲国产高清在线一区二区三 | www.精华液| 一个人免费在线观看的高清视频| 国产精品99久久99久久久不卡| 热99re8久久精品国产| 丰满的人妻完整版| 亚洲自偷自拍图片 自拍| 日韩有码中文字幕| 国产精品久久久人人做人人爽| 成人18禁高潮啪啪吃奶动态图| 亚洲狠狠婷婷综合久久图片| 久久精品91蜜桃| 美女扒开内裤让男人捅视频| 国产亚洲欧美在线一区二区| 欧美成人午夜精品| 中文字幕色久视频| 久久香蕉国产精品| 亚洲五月色婷婷综合| 狂野欧美激情性xxxx| 久久香蕉国产精品|