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

    黃原膠水溶液管道流動減阻特性的試驗

    2015-12-21 07:43:24李昌烽禹燕飛趙文斌李明義侯金亮

    李昌烽,禹燕飛,趙文斌,李明義,侯金亮

    黃原膠水溶液管道流動減阻特性的試驗

    李昌烽,禹燕飛,趙文斌,李明義,侯金亮

    (江蘇大學(xué)能源與動力工程學(xué)院,江蘇 鎮(zhèn)江 212013)

    對連續(xù)循環(huán)光滑管道(直徑分別為5,10,20 mm)中黃原膠水溶液流動的減阻特性進(jìn)行了試驗,分析了黃原膠減阻的濃度效應(yīng)、管徑效應(yīng)以及抗剪切特性,得到了減阻率與黃原膠水溶液濃度的關(guān)系曲線、雷諾數(shù)(Re)與減阻率的影響曲線以及減阻率與剪切時間的變化曲線.結(jié)果表明:黃原膠是很好的減阻劑,在較低Re流動時黃原膠在相對較小直徑(5,10 mm)管道中表現(xiàn)為B型減阻特性,而在較大管徑(20 mm)中則為A型減阻;在高Re時黃原膠水溶液在3種管徑管道流動中皆表現(xiàn)為B型減阻,由于其主要具有B型減阻特性,使得黃原膠具有較好的抗剪切特性.

    黃原膠溶液;管道流動;湍流減阻;A型減阻;B型減阻

    添加劑湍流減阻技術(shù)是指在管道中的液體湍流流動中加入微量添加劑,使湍流摩擦阻力顯著降低的流動輸送技術(shù).1948年B.A.Toms在第一屆國際流變學(xué)會議上發(fā)表了關(guān)于高分子聚合物減阻的論文[1].研究發(fā)現(xiàn):在湍流內(nèi)流的氯苯中,溶解少量的聚甲基丙烯酸甲酯(PMMA),可大幅度降低流動阻力.因此高分子減阻又稱為Toms效應(yīng).自從發(fā)現(xiàn)高分子減阻效應(yīng)以來,P.S.Virk等[1]對減阻效應(yīng)進(jìn)行了廣泛研究,在此期間發(fā)現(xiàn)了最大減阻漸近線(maximum drag reduction asymptote,MDRA)的存在,找到了影響高分子減阻效應(yīng)的許多重要參數(shù),包括高聚物的分子量、分子長度、柔性和伸展性等.20世紀(jì)90年代以來,隨著測試與顯示技術(shù)的發(fā)展,各國學(xué)者對減阻流體進(jìn)行了大量試驗研究.F.T.Pinho等[2-3]用激光多普勒測速儀(laser doppler velocimetry,LDV)測量了管道和槽道流動的湍流統(tǒng)計量,發(fā)現(xiàn)高分子并不是抑制所有方向上的湍流運(yùn)動,而是減弱了法向方向的猝發(fā)和湍流強(qiáng)度,但流向方向的湍流強(qiáng)度有所加強(qiáng).M.D.Warholic等[4-6]用粒子圖像測速技術(shù)(particle image velocimetry,PIV)對高分子減阻流體的流場進(jìn)行測量.M.D.Warholic等[4]采用PIV技術(shù)測量了2維湍流流場,證實了在最大減阻漸近線附近(MDRA)下存在零雷諾切應(yīng)力的現(xiàn)象.O.Cadot等[7]通過有肋條和無肋條的旋轉(zhuǎn)圓盤試驗,研究了邊界層對高分子流體減阻的影響,表明只有邊界層的存在才能發(fā)生減阻.張根廣等[8]研究了細(xì)菌多糖黃原膠的減阻特性,分別在15,22℃下對黃原膠多種濃度進(jìn)行了減阻測試,得到了范寧阻力系數(shù)和雷諾數(shù)的關(guān)系曲線.G.C.Garwood等[9]通過試驗發(fā)現(xiàn)聚合物在高剪切力作用下出現(xiàn)裂變,高分子長鏈斷裂,使減阻消失且不可恢復(fù).A. Jaafar等[10]對天然高聚物(殼聚糖、黃原膠)溶液的減阻特性進(jìn)行了研究,并用LDA技術(shù)測量了減阻溶液的雷諾應(yīng)力.姚世勇等[11]用2維高分辯率TRPIV對高分子聚合物溶液減阻性能進(jìn)行了試驗研究,發(fā)現(xiàn)高分子聚合物溶液能明顯地抑制壁面湍流相關(guān)結(jié)構(gòu)的噴射和掃掠的強(qiáng)度.

    此外,P.S.Virk等[12-13]詳細(xì)討論過 2種減阻類型,發(fā)現(xiàn)聚合物溶液中NaCl質(zhì)量分?jǐn)?shù)從10-8到10-11的變化可能會導(dǎo)致減阻類型從A型變化到B型.這種差異與聚合物分子的構(gòu)象有關(guān),這反過來又似乎取決于溶液中離子的濃度(例如聚丙烯酰胺溶液中的NaCl的量)等的不同而呈現(xiàn)很大差異.A型減阻的聚合物分子是典型的卷曲狀(在高含鹽量情況下),這需要在減阻起始點(diǎn)(onset)前存在一定程度的壁面切應(yīng)力.由于流速進(jìn)一步增大減阻增強(qiáng),直到后來降解變得明顯,減阻隨之下降.另一方面,B型減阻的聚合物分子是伸展?fàn)顟B(tài)的(例如具有非常低的含鹽量的聚丙烯酰胺溶液),流動狀態(tài)從層流過渡到湍流之后立即表現(xiàn)出漸近減阻.隨著流速增大,漸近摩擦系數(shù)保持不變,直到達(dá)到一個減阻再起始點(diǎn)(retro-onset),在這之后對于給定管徑,隨著速度的進(jìn)一步增加減阻水平保持大致恒定.

    筆者選擇細(xì)菌多糖黃原膠,在不同直徑的光滑有機(jī)玻璃管道中,進(jìn)行減阻特性試驗,主要研究黃原膠水溶液的濃度效應(yīng)、管徑效應(yīng)以及流速對減阻率DR的影響,為研究天然高分子聚合物的減阻機(jī)理提供參考依據(jù).

    1 試驗裝置和方法

    1.1 試驗裝置

    試驗在設(shè)計的裝置[14](如圖1所示)上進(jìn)行,裝置和回路主要由儲液箱、離心泵、變頻器、電磁流量計、恒壓水箱、壓力表、差壓變送器、測量槽以及各種閥門管路組成.其中測試段AE包括1組3個不同直徑的管道(5,10,20 mm),在壓力測量點(diǎn)B之前,為充分達(dá)到穩(wěn)定狀態(tài),需要足夠長的發(fā)展段AB,根據(jù)F.M.White[15]提出的經(jīng)驗公式,AB長度應(yīng)滿足不小于138倍的管徑,因此5,10,20 mm這3種直徑的管道發(fā)展段分別為0.7,1.4,2.8 m,3種管道的測試段BC長度均為2.0 m.本試驗采用差壓變送器測量測試段BC兩端的壓差,其測量精度可達(dá)0.25%;采用體積法與電磁流量計2種方式測量流量,相對誤差小于0.1%,由變頻器調(diào)節(jié)離心泵的轉(zhuǎn)速,并由閥門控制管道流速.

    圖1 試驗裝置示意圖

    1.2 試驗方法

    對細(xì)菌黃原膠的減阻特性進(jìn)行試驗研究,試驗中的控制變量為添加劑質(zhì)量分?jǐn)?shù)和流體速度.試驗過程中需測量體積流量Q以及B,C點(diǎn)的壓力差Δp,試驗的具體操作步驟:①保持溶液濃度不變,調(diào)節(jié)流速,分別測量3種直徑管道中在同一濃度下不同流速對應(yīng)的壓力差Δp;②保持流速恒定,改變?nèi)芤簼舛?,分別測量3種直徑管道中不同濃度時同一流速對應(yīng)的壓力差Δp;③在溶液保持最佳減阻濃度時,對減阻溶液進(jìn)行連續(xù)循環(huán)剪切,每隔1 h記錄1次Q和Δp.

    1.3 減阻液的配制

    試驗所用黃原膠是由淄博中軒生物有限公司提供的商業(yè)食品級產(chǎn)品.配制質(zhì)量分?jǐn)?shù)分別為(50~550)×10-6(每次增加50×10-6)的黃原膠水溶液.試驗所需的黃原膠水溶液在試驗前3天,用自來水將固體黃原膠配制成質(zhì)量分?jǐn)?shù)為1 000×10-6母液,靜置3 d,待試驗開始前1天,用自來水將其稀釋成所需濃度溶液,靜置1 d,保證試驗時黃原膠水溶液充分溶解并混合均勻.

    1.4 試驗數(shù)據(jù)處理方法

    減阻率的具體計算需要用到以下公式:圓管沿程阻力系數(shù)為

    式中:ρ為流體的密度;l為BC段長度;D為管徑.范寧摩擦阻力系數(shù)為

    Re的計算公式為

    式中:u為管道斷面的平均速度;μ為溶液的動力黏度;υ為溶液的運(yùn)動黏度.試驗中雷諾數(shù)的變化范圍為2 300~50 000.

    同一Re下減阻率DR的表達(dá)式為

    式中:f為高聚物溶液的范寧摩擦阻力系數(shù):fw為清水牛頓流體的范寧摩擦阻力系數(shù).

    1.5 管徑的率定

    圓管流壓降為

    體積流量為

    將式(5),(6)聯(lián)立,并對等號兩邊取對數(shù),然后再微分可得

    式(7)等號右邊的各項可看作各物理量單項測量的相對誤差.從式(7)可以看出:管徑和流量的測量誤差傳播系數(shù)分別為5和2,即如果試驗管徑相差1%,則由此項帶來的誤差就達(dá)5%.因此,精確率定試驗管徑非常重要.

    在進(jìn)行添加劑減阻試驗之前,先進(jìn)行清水試驗,根據(jù)已知的牛頓流體相關(guān)理論來率定管徑.從理論上講,在層流區(qū),清水的阻力系數(shù)符合Hagen-Poiseuille定律:

    式中:Re包含管徑,因此可以精確計算出試驗管道的直徑.然后用清水進(jìn)行湍流試驗,根據(jù)層流試驗計算得到的管徑,以及湍流試驗求出的沿程阻力系數(shù)λ,驗證其是否滿足Prandtl-Karman定律:

    基于上述管徑率定原則,對20℃的自來水進(jìn)行試驗,精確率定后管徑分別為5.28,11.30,20.00 mm.率定結(jié)果如圖2所示,用率定后管徑處理的試驗結(jié)果層流區(qū)基本符合Hagen-Poiseuille定律,湍流區(qū)也較好地符合Prandtl-Karman定律,說明試驗管徑的精確度可以滿足本試驗的要求,為進(jìn)一步進(jìn)行減阻試驗提供了可靠的依據(jù).

    圖2 清水試驗管徑率定結(jié)果

    2 試驗結(jié)果與討論

    2.1 濃度(質(zhì)量分?jǐn)?shù))效應(yīng)

    高聚物的減阻效果與其在水溶液中的質(zhì)量分?jǐn)?shù)有很大關(guān)系,為了評價黃原膠水溶液的減阻效果,在雷諾數(shù)為25 000時,測量質(zhì)量分?jǐn)?shù)w范圍為(50~550)×10-6的黃原膠水溶液的減阻性能,試驗結(jié)果如圖3所示.

    圖3 不同質(zhì)量分?jǐn)?shù)黃原膠水溶液的減阻率

    由圖3可以看出:黃原膠水溶液的濃度效應(yīng)相當(dāng)顯著,在低質(zhì)量分?jǐn)?shù)時,減阻效果較低,隨著質(zhì)量分?jǐn)?shù)不斷升高,減阻率明顯增大,高質(zhì)量分?jǐn)?shù)時,減阻率增長趨勢變緩,直到黃原膠質(zhì)量分?jǐn)?shù)達(dá)到500×10-6,減阻率基本保持不變,不再隨著質(zhì)量分?jǐn)?shù)的增大而增加,5,10,20 mm管徑黃原膠水溶液流動對應(yīng)的最大減阻率分別為55.1%,47.6%,46.4%.這也從另一個方面表現(xiàn)出管徑效應(yīng),即管徑越小,減阻效果越好.

    2.2 減阻率DR與雷諾數(shù)Re的關(guān)系

    為了探索黃原膠水溶液的減阻率DR與雷諾數(shù)的關(guān)系,測量了不同質(zhì)量分?jǐn)?shù)溶液在不同雷諾數(shù)下的減阻率,結(jié)果如圖4所示.圖4a,b中,在同一質(zhì)量分?jǐn)?shù)下,較小管徑(5,10 mm)流動中的減阻率基本不隨雷諾數(shù)的增加而變化.由圖4c可以看出:低雷諾數(shù)時(約15 000以下)隨著雷諾數(shù)的增加,減阻率不斷增大,當(dāng)雷諾數(shù)到一定程度(15 000~20 000)后,減阻率基本不再受雷諾數(shù)影響.

    圖4 不同質(zhì)量分?jǐn)?shù)黃原膠水溶液在不同雷諾數(shù)下的減阻率曲線

    2.3 管徑效應(yīng)

    P.S.Virk等[12-13]依據(jù)不同的減阻特性把高分子減阻劑分為A型和B型減阻劑,A型和B型減阻中的阻力系數(shù)變化特性如圖5所示.

    為方便說明和對比,把3種不同直徑(5,10,20 mm)管道中減阻試驗數(shù)據(jù)也整理成與圖5一樣的Prandtl-Karman曲線,如圖6所示.

    對比圖5可發(fā)現(xiàn):在圖6a,b中,小管徑(5和10 mm)流動中摩擦阻力系數(shù)曲線表現(xiàn)出明顯的B型減阻特性,在Ref1/2較低時,阻力系數(shù)曲線接近平行于L線(牛頓流體層流摩擦阻力系數(shù)曲線);在Ref1/2較高時,阻力系數(shù)曲線基本平行于N線(牛頓流體湍流摩擦阻力系數(shù)曲線).而在20 mm大直徑的管流中,在低速時,黃原膠的阻力系數(shù)曲線卻表現(xiàn)為A型減阻,它與N線有1個交點(diǎn).在此后其摩擦阻力系數(shù)曲線呈扇形分布,曲線斜率隨著黃原膠水溶液濃度的增加而增大,減阻水平隨著Ref1/2的增大而提高.而在較高流速下,其阻力系數(shù)曲線又基本平行于N線,表現(xiàn)為B型減阻.這與K.Gasljevic等[16]的試驗觀測結(jié)果一致,給定的流體在給定的較大管道流動中可以從低Re下的A型減阻轉(zhuǎn)變?yōu)楦逺e下的B型減阻.

    圖5 A型和B型減阻特性示意圖

    圖6 黃原膠水溶液在不同直徑管道中的摩擦阻力系數(shù)曲線

    黃原膠屬于生物高聚物,不同于合成高聚物,由于糖鏈上含有羧基和硫酸酯從而具有電解性,通常是強(qiáng)陰離子.它們所帶的強(qiáng)負(fù)電荷使得高分子聚合物受到的電荷之間的排斥力強(qiáng)于其他分子對它的親和力,從而使得聚合物分子在流體中維持著一定的平衡拉伸狀態(tài).在小管徑(5,10 mm)流動中,在給定Re下流動強(qiáng)度(剪切)可能已經(jīng)足夠大來拉伸黃原膠分子,使得層流-湍流轉(zhuǎn)捩區(qū)后就立即產(chǎn)生B型漸近減阻.但對大管徑(20 mm)來說,相比小管徑情形,在給定Re下剪切強(qiáng)度是小的,湍流絕對尺度是大的,在層流-湍流轉(zhuǎn)捩區(qū)后仍然需要剪切來激發(fā)黃原膠溶液減阻效應(yīng),所以這時表現(xiàn)為A型減阻;隨著雷諾數(shù)的提高,大管徑中的流動剪切強(qiáng)度增加,足以來拉伸黃原膠分子,減阻水平隨著分子進(jìn)一步伸展而增大.試驗結(jié)果表明:在較高Re下黃原膠水溶液在3種大小管徑管道流動中都表現(xiàn)為B型減阻,摩擦阻力與管徑大小無關(guān),減阻程度也近似與Re無關(guān).所以在大管道(20 mm)流動中,在低、高Re下黃原膠水溶液從A型轉(zhuǎn)變?yōu)锽型減阻特性,應(yīng)該理解為黃原膠分子與流動之間的相互作用程度的不同,而不是減阻本質(zhì)不同的現(xiàn)象.

    2.4 抗剪切性能

    高聚物水溶液在連續(xù)循環(huán)系統(tǒng)中,主要受泵的機(jī)械剪切作用,大分子鏈斷裂,發(fā)生降解,不斷失去原有的減阻效果,這是影響高聚物在工程上廣泛應(yīng)用的障礙之一.因此,抗剪切性能是評價減阻劑性能的重要指標(biāo)之一.在雷諾數(shù)為 40 000左右,對質(zhì)量分?jǐn)?shù)為500×10-6的黃原膠水溶液進(jìn)行連續(xù)循環(huán)剪切.第1天連續(xù)剪切6 h,靜置18 h;第2天連續(xù)剪切11 h,靜置18 h;第3天連續(xù)剪切10 h,靜置18 h;第4天連續(xù)剪切9 h.在連續(xù)剪切過程中,每隔1 h記錄1次數(shù)據(jù),試驗結(jié)果如圖7所示.

    圖7 質(zhì)量分?jǐn)?shù)為500×10-6黃原膠水溶液在機(jī)械剪切后的減阻率

    由圖7可以看出:在10 mm直徑的管道中,黃原膠水溶液降解較慢,直到第3天才表現(xiàn)出顯著降解,剪切45h之后,減阻率從最初的47.6%降到38.2%;而在20mm直徑的管道中,在第2天黃原膠水溶液就表現(xiàn)出明顯降解,到第3天降解趨勢變緩,剪切45h之后,減阻率從46.4%降到35.5%;在每次長時間的靜置之后,減阻率都有個小幅度的回升,但是在連續(xù)剪切作用力下,回升很快消失,溶液繼續(xù)發(fā)生降解,說明黃原膠水溶液的機(jī)械降解是不完全可逆的.綜上所述,黃原膠具有較好的抗降解性能,且機(jī)械剪切對較大管徑流動減阻率的影響稍大一些.

    3 結(jié) 論

    黃原膠水溶液具有明顯的濃度效應(yīng),即濃度越高減阻性能越好,在試驗范圍內(nèi)得到質(zhì)量分?jǐn)?shù)為500×10-6為它的最大飽和減阻濃度;黃原膠水溶液的減阻性能與管徑有關(guān),管徑越小,減阻性能越好;黃原膠水溶液在小直徑管道(5,10 mm)中表現(xiàn)為 B型減阻;在大直徑管道(20 mm)中,低雷諾數(shù)時為A型減阻,高雷諾數(shù)時為B型減阻;黃原膠主要具有B型減阻特性,有較好的抗剪切特性,是一種很有潛力的天然高聚物減阻劑.

    [1] Virk P S,Merrill EW,Mickley H S,et al.The Toms phenomenon:turbulent pipe flow of dilute polymer solutions[J].Journal of Fluid Mechanics,1967,30:305-328.

    [2] Pinho F T,Whitelaw JH.Flow of non-Newtonian fluids in a pipe[J].Journal of Non-Newtonian Fluid Mechanics,1990,34:129-144.

    [3] Harder K J,Tiederman W G.Drag reduction and turbulent structure in two-dimensional channel flows[J]. Phil Trans Math Phys Eng Sci,1991,336:19-34.

    [4] Warholic M D,Heist D K,Katcher M,et al.A study with particle-image velocimetry of the influence of drag reducing polymers on the structure of turbulence[J]. Experiments in Fluids,2001,31:474-483.

    [5] White C M,Somandepalli V S R,Mungal M G.The turbulence structure of drag reduced boundary layer flow[J].Experiments in Fluids,2004,36:62-69.

    [6] Hou Y X,Somandepalli V SR,Mungal M G.A technique to determine total shear stress and polymer stress profiles in drag reduced boundary layer flows[J].Experiments in Fluids,2006,40:589-600.

    [7] Cadot O,Bonn D,Douady S.Turbulent drag reduction in a closed flow system:boundary layer versus bulk effects[J].Physics of Fluids,1998,10(2):426-436.

    [8] 張根廣,張鳴遠(yuǎn),楊萬鷹,等.黃原膠溶液減阻特性試驗[J].西安交通大學(xué)學(xué)報,2007,41(9):1092-1095. Zhang Genguang,Zhang Mingyuan,Yang Wanying,et al.Drag reduction in turbulent pipe flows of aqueous xanthan gum solutions[J].Journal of Xi'an Jiaotong U-niversity,2007,41(9):1092-1095.(in Chinese)

    [9] Garwood G C,Winkel E S,Vanapalli S,et al.Drag reduction by a homogenous polymer solution in large diameter high shear pipe flow[C]∥2nd International Symposium on Seawater Drag Reduction,2005:1-7.

    [10] Jaafar A,Poole R J.Drag reduction of biopolymer flows[J].Journal of Applied Sciences,2011,11(9):1544-1551.

    [11] 姚世勇,管新蕾,姜 楠.高分子聚合物溶液減阻特性的TR-PIV實驗研究[C]∥第7屆全國流體力學(xué)學(xué)術(shù)會議論文摘要集.北京:中國力學(xué)學(xué)會流體力學(xué)專業(yè)委員會,2012:125-126.

    [12] Virk P S,Wagger D L.Aspects of mechanisms in type B drag reduction[C]∥Structure of Turbulence and Drag Reduction.Berlin:Springer,1989:201-203.

    [13] Virk P S,Wagger D L.The effect of salinity on turbulent drag reduction by polyelectrolyte additives at high Reynolds numbers[C]∥Addendum to the Proceedings of the International Symposium on Seawater Drag Reduction,1998:128-132.

    [14] 禹燕飛,李昌烽,趙文斌,等.一種高分子聚合物溶液全流態(tài)減阻特性實驗裝置:中國,201320028552.4[P].2013-04.

    [15] White F M.Fluid Mechanics[M].New York:McGraw-Hill Series,2010.

    [16] Gasljevic K,Aguilar G,Matthys E F.On two distinct types of drag-reducing fluids,diameter scaling,and turbulent profiles[J].J Non-Newtonian Fluid Mech,2001,96:405-425.

    (責(zé)任編輯 賈國方)

    Experiment on drag reduction characteristics of xanthan gum solution in pipe flow

    Li Changfeng,Yu Yanfei,Zhao Wenbin,LiMingyi,Hou Jinliang
    (School of Energy and Power Engineering,Jiangsu University,Zhenjiang,Jiangsu 212013,China)

    The systematic experiments on drag reduction characteristics of xanthan gum solution in continuous cycle smooth pipes with diameters of 5,10 and 20 mm were respectively conducted.The concentration and pipe diameter effects and the withstanding mechanical degradation of xanthan gum solution on drag reduction characteristics were analyzed to obtain relationship between drag reductions and concentration of xanthan gum solution,influence of Reynolds number on drag reductions and percent drag reduction changing with shearing time.The experimental results show that xanthan gum solution is good drag reduction agent.At low Reynolds number,xanthan gum solution performs type-B drag reduction characteristics in small pipes with diameter of 5 mm and 10 mm,but shows type-A drag reductions characteristics in large pipe with diameter of 20 mm.At high Reynolds number,the xanthan gum solution performs type-B drag reduction characteristics in all pipe flows with three diameters.Xanthan gum solution exhibits high resistance to mechanical degradation because of the essential type-B drag reduction characteristics.

    xanthan gum solution;pipe flow;turbulent drag reduction;type-A drag reduction;type-B drag reduction

    O130.25

    A

    1671-7775(2015)01-0030-06

    李昌烽,禹燕飛,趙文斌,等.黃原膠水溶液管道流動減阻特性的試驗[J].江蘇大學(xué)學(xué)報:自然科學(xué)版,2015,36(1):030-035.

    10.3969/j.issn.1671-7775.2015.01.006

    2014-02-27

    國家自然科學(xué)基金資助項目(10672069,11072091);教育部科學(xué)技術(shù)研究重點(diǎn)項目(210078)

    李昌烽(1969-),男,江蘇徐州人,教授(cfli@ujs.edu.cn),主要從事復(fù)雜流體動力學(xué)、湍流減阻與流動控制研究.禹燕飛(1988-),女,甘肅渭源人,碩士研究生(yu.yanfei1988@163.com),主要從事湍流減阻與流動控制研究.

    精品一区二区免费观看| 3wmmmm亚洲av在线观看| 黄色配什么色好看| 中文资源天堂在线| 老熟妇乱子伦视频在线观看| 不卡一级毛片| 久久热精品热| 日本与韩国留学比较| 久久国产乱子免费精品| av又黄又爽大尺度在线免费看 | 亚洲欧美日韩高清在线视频| 国产久久久一区二区三区| 女的被弄到高潮叫床怎么办| 99热这里只有精品一区| 国产日韩欧美在线精品| 国产日本99.免费观看| 一本久久精品| 亚洲最大成人手机在线| 中文字幕精品亚洲无线码一区| 亚洲一区高清亚洲精品| 国产精品永久免费网站| 天堂网av新在线| 亚洲精华国产精华液的使用体验 | 青春草亚洲视频在线观看| а√天堂www在线а√下载| 日本免费a在线| 三级毛片av免费| 免费观看精品视频网站| 日本-黄色视频高清免费观看| 日韩精品有码人妻一区| 中国美女看黄片| 亚洲人成网站在线播| 亚洲在久久综合| 午夜福利在线在线| 舔av片在线| 国产麻豆成人av免费视频| 成人综合一区亚洲| 久久精品91蜜桃| 亚洲综合色惰| 一本久久精品| av专区在线播放| 国产成人精品久久久久久| 日日撸夜夜添| 好男人在线观看高清免费视频| 老女人水多毛片| 一级黄色大片毛片| 国产高清视频在线观看网站| 久久综合国产亚洲精品| 在线免费观看的www视频| 日本一本二区三区精品| 久久精品夜色国产| 桃色一区二区三区在线观看| 亚洲国产精品国产精品| 亚洲三级黄色毛片| 亚洲,欧美,日韩| 欧美区成人在线视频| 国内少妇人妻偷人精品xxx网站| 日本色播在线视频| 深夜精品福利| 国产精品女同一区二区软件| 亚洲va在线va天堂va国产| 少妇人妻精品综合一区二区 | 欧美+日韩+精品| 男的添女的下面高潮视频| 亚洲欧美成人综合另类久久久 | videossex国产| 九九爱精品视频在线观看| 欧美不卡视频在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费观看精品视频网站| 人人妻人人澡欧美一区二区| 日日干狠狠操夜夜爽| 国产精品蜜桃在线观看 | 嫩草影院精品99| 久久久久九九精品影院| 97在线视频观看| 久久精品国产亚洲av天美| 女人被狂操c到高潮| 在线免费观看不下载黄p国产| 热99re8久久精品国产| 成年av动漫网址| 日本一本二区三区精品| 搞女人的毛片| 我的老师免费观看完整版| 免费观看a级毛片全部| 国产亚洲91精品色在线| 狠狠狠狠99中文字幕| 国产片特级美女逼逼视频| 九九爱精品视频在线观看| 在线播放无遮挡| 12—13女人毛片做爰片一| 国产伦精品一区二区三区视频9| 免费av毛片视频| 中国美女看黄片| 久久久a久久爽久久v久久| 国产高清不卡午夜福利| 久久久久久久午夜电影| 精品日产1卡2卡| 91狼人影院| 国产 一区精品| av在线播放精品| 国产91av在线免费观看| 日本黄大片高清| 亚洲国产精品国产精品| 欧美+亚洲+日韩+国产| 天堂√8在线中文| 国产人妻一区二区三区在| 看片在线看免费视频| 亚洲精品久久久久久婷婷小说 | 精品国产三级普通话版| 国产老妇伦熟女老妇高清| 久久久精品大字幕| 国产极品天堂在线| av在线天堂中文字幕| 久久精品久久久久久久性| 熟妇人妻久久中文字幕3abv| 我的女老师完整版在线观看| 热99re8久久精品国产| 国产伦一二天堂av在线观看| 少妇丰满av| 午夜老司机福利剧场| 成人综合一区亚洲| 日本色播在线视频| 欧美xxxx黑人xx丫x性爽| 国产日韩欧美在线精品| 国产真实伦视频高清在线观看| 女人十人毛片免费观看3o分钟| 欧洲精品卡2卡3卡4卡5卡区| 日本免费一区二区三区高清不卡| 久久6这里有精品| 久久精品国产鲁丝片午夜精品| 国产一级毛片七仙女欲春2| 国产精品三级大全| 麻豆精品久久久久久蜜桃| 国产免费男女视频| 97超视频在线观看视频| 亚洲成a人片在线一区二区| 亚洲欧美日韩东京热| 日韩欧美一区二区三区在线观看| 老司机影院成人| 赤兔流量卡办理| 国产午夜精品久久久久久一区二区三区| 69av精品久久久久久| 国产乱人偷精品视频| 不卡一级毛片| 亚洲欧美日韩卡通动漫| 日本一本二区三区精品| 日本色播在线视频| 在线免费十八禁| 亚洲欧美成人精品一区二区| 精品一区二区免费观看| 99热这里只有是精品在线观看| 日韩,欧美,国产一区二区三区 | 午夜福利视频1000在线观看| 日日摸夜夜添夜夜添av毛片| 久久精品久久久久久噜噜老黄 | 欧美一级a爱片免费观看看| 中文字幕av成人在线电影| 综合色丁香网| 高清在线视频一区二区三区 | 97人妻精品一区二区三区麻豆| av卡一久久| 午夜精品一区二区三区免费看| 欧美+亚洲+日韩+国产| 观看美女的网站| 午夜福利成人在线免费观看| 中文字幕久久专区| 精品久久久久久久久av| 久久99热这里只有精品18| 九色成人免费人妻av| 一级毛片我不卡| 亚洲av不卡在线观看| 亚洲性久久影院| 日本与韩国留学比较| 在线观看美女被高潮喷水网站| 少妇人妻精品综合一区二区 | 成人亚洲欧美一区二区av| 久久精品夜夜夜夜夜久久蜜豆| 搡老妇女老女人老熟妇| 麻豆av噜噜一区二区三区| 中文字幕av在线有码专区| 国产爱豆传媒在线观看| 五月玫瑰六月丁香| 亚洲最大成人手机在线| 成人美女网站在线观看视频| 日本-黄色视频高清免费观看| АⅤ资源中文在线天堂| 亚洲婷婷狠狠爱综合网| 别揉我奶头 嗯啊视频| 你懂的网址亚洲精品在线观看 | 国产亚洲欧美98| 国产精品一区二区在线观看99 | 波多野结衣高清无吗| 欧美在线一区亚洲| 少妇的逼水好多| 亚洲欧美成人综合另类久久久 | 亚洲人成网站在线观看播放| 久久久久久久久久黄片| 日本黄色片子视频| 日本黄大片高清| 欧美日韩乱码在线| 免费黄网站久久成人精品| 天天一区二区日本电影三级| 欧美日韩在线观看h| 欧美激情在线99| 一级av片app| 欧美区成人在线视频| 国产精品三级大全| 一进一出抽搐动态| 亚洲欧美日韩高清在线视频| 插逼视频在线观看| 亚洲自偷自拍三级| 国产精品美女特级片免费视频播放器| 99热这里只有是精品50| 欧美精品一区二区大全| 91久久精品国产一区二区三区| 欧美丝袜亚洲另类| 亚洲欧美日韩东京热| 日本色播在线视频| 免费观看人在逋| 蜜桃亚洲精品一区二区三区| 欧美+日韩+精品| 亚洲真实伦在线观看| 一级毛片电影观看 | 国产高清有码在线观看视频| 亚洲不卡免费看| 伦理电影大哥的女人| 久久精品久久久久久噜噜老黄 | 在线播放国产精品三级| 国产黄片美女视频| 亚洲欧美成人综合另类久久久 | 国产精品一二三区在线看| 亚洲中文字幕日韩| 色播亚洲综合网| a级一级毛片免费在线观看| 男女做爰动态图高潮gif福利片| 三级毛片av免费| 美女脱内裤让男人舔精品视频 | av专区在线播放| 国产精品av视频在线免费观看| 热99re8久久精品国产| 美女脱内裤让男人舔精品视频 | 少妇的逼水好多| 久久久久久伊人网av| 欧美丝袜亚洲另类| 精品人妻偷拍中文字幕| 日韩高清综合在线| 国产成人一区二区在线| 欧美一区二区国产精品久久精品| 日韩视频在线欧美| 色吧在线观看| 高清毛片免费看| 国产极品天堂在线| 国产美女午夜福利| 色哟哟哟哟哟哟| 男女视频在线观看网站免费| 特级一级黄色大片| 乱人视频在线观看| 麻豆一二三区av精品| 国产精品久久久久久精品电影| 久久99热6这里只有精品| 亚洲国产高清在线一区二区三| av在线播放精品| 波多野结衣巨乳人妻| 国产精品,欧美在线| 干丝袜人妻中文字幕| 一区二区三区高清视频在线| 亚洲最大成人手机在线| 亚洲中文字幕日韩| 一级毛片我不卡| 欧美人与善性xxx| 亚洲熟妇中文字幕五十中出| 十八禁国产超污无遮挡网站| 国产精品av视频在线免费观看| 色视频www国产| 日韩欧美 国产精品| 成人无遮挡网站| 国产精品乱码一区二三区的特点| 国产精品精品国产色婷婷| 精品久久久噜噜| 精品国产三级普通话版| 中文字幕制服av| 波多野结衣高清无吗| 亚洲欧美日韩高清专用| 在线观看66精品国产| 观看美女的网站| 欧美变态另类bdsm刘玥| 亚洲内射少妇av| 国产又黄又爽又无遮挡在线| 久久久久久久久大av| 亚洲av电影不卡..在线观看| 欧美日韩精品成人综合77777| 51国产日韩欧美| 你懂的网址亚洲精品在线观看 | 99久久人妻综合| 国产 一区精品| 99热全是精品| 国产三级中文精品| 可以在线观看的亚洲视频| 免费观看在线日韩| 国产亚洲av嫩草精品影院| 久久久久久久久久久免费av| av视频在线观看入口| 日韩欧美在线乱码| 综合色av麻豆| 99久久精品一区二区三区| av免费观看日本| 国内精品宾馆在线| 国产欧美日韩精品一区二区| 身体一侧抽搐| 91久久精品国产一区二区三区| 性欧美人与动物交配| 亚洲欧美日韩东京热| 亚洲综合色惰| 亚洲欧美清纯卡通| 亚洲五月天丁香| 亚洲精品国产av成人精品| 在线国产一区二区在线| 亚洲欧美精品专区久久| 中文字幕熟女人妻在线| 麻豆av噜噜一区二区三区| 九九在线视频观看精品| av在线天堂中文字幕| 真实男女啪啪啪动态图| av免费观看日本| 亚洲国产色片| 国产免费男女视频| 欧美bdsm另类| 国产人妻一区二区三区在| 三级男女做爰猛烈吃奶摸视频| 色哟哟·www| 国产精品永久免费网站| 亚洲激情五月婷婷啪啪| 精品免费久久久久久久清纯| 最近最新中文字幕大全电影3| 亚洲真实伦在线观看| 亚洲精品久久国产高清桃花| 国产综合懂色| 波多野结衣高清无吗| 亚洲高清免费不卡视频| 真实男女啪啪啪动态图| 男女做爰动态图高潮gif福利片| 免费看日本二区| 在线观看免费视频日本深夜| 一个人免费在线观看电影| 久久久精品大字幕| 欧美日本亚洲视频在线播放| 91久久精品国产一区二区成人| 天堂网av新在线| 黄片无遮挡物在线观看| 三级男女做爰猛烈吃奶摸视频| 最新中文字幕久久久久| 亚洲五月天丁香| 成人综合一区亚洲| 午夜精品国产一区二区电影 | 免费看av在线观看网站| 亚洲中文字幕一区二区三区有码在线看| 欧美性感艳星| 亚洲丝袜综合中文字幕| 免费av毛片视频| 亚洲欧洲日产国产| 色播亚洲综合网| 又黄又爽又刺激的免费视频.| 色噜噜av男人的天堂激情| 国产真实乱freesex| 两性午夜刺激爽爽歪歪视频在线观看| 成人特级av手机在线观看| 欧美极品一区二区三区四区| 小说图片视频综合网站| 国产一区二区三区在线臀色熟女| 精品久久久噜噜| avwww免费| 91久久精品国产一区二区三区| 99久久成人亚洲精品观看| 夫妻性生交免费视频一级片| 一区福利在线观看| 成人性生交大片免费视频hd| 久久久久久久亚洲中文字幕| 一进一出抽搐gif免费好疼| 一个人免费在线观看电影| 91精品国产九色| 亚洲成人久久性| 欧美成人a在线观看| 赤兔流量卡办理| 亚洲av男天堂| 国产单亲对白刺激| 身体一侧抽搐| av卡一久久| 国产亚洲精品av在线| 高清日韩中文字幕在线| 国产单亲对白刺激| 欧美精品一区二区大全| 久久韩国三级中文字幕| 免费一级毛片在线播放高清视频| 久久精品国产亚洲av香蕉五月| 日韩一区二区三区影片| 亚洲国产高清在线一区二区三| 免费看美女性在线毛片视频| 久久精品国产鲁丝片午夜精品| 国产精品一区二区三区四区免费观看| 黄片wwwwww| 狂野欧美白嫩少妇大欣赏| 成人高潮视频无遮挡免费网站| 男女边吃奶边做爰视频| 啦啦啦观看免费观看视频高清| 搡女人真爽免费视频火全软件| 亚洲欧美日韩高清专用| 婷婷亚洲欧美| 亚洲av免费在线观看| 婷婷色av中文字幕| 在线a可以看的网站| 一级二级三级毛片免费看| 十八禁国产超污无遮挡网站| www日本黄色视频网| 国产精品一区www在线观看| 国产精品1区2区在线观看.| 黄色视频,在线免费观看| 午夜精品一区二区三区免费看| 少妇丰满av| 日韩一区二区三区影片| 成人毛片60女人毛片免费| 亚洲婷婷狠狠爱综合网| 特级一级黄色大片| 国产一区二区激情短视频| 久久久a久久爽久久v久久| 日韩精品有码人妻一区| 国产精品一二三区在线看| 国产精品综合久久久久久久免费| 99riav亚洲国产免费| 内射极品少妇av片p| 久久久久久久午夜电影| 亚洲第一电影网av| 日本色播在线视频| 美女cb高潮喷水在线观看| 丰满的人妻完整版| 成人亚洲精品av一区二区| 只有这里有精品99| 中国美女看黄片| 国产成人一区二区在线| 综合色av麻豆| 亚洲丝袜综合中文字幕| 亚洲国产欧洲综合997久久,| 国产在视频线在精品| 欧美日韩在线观看h| 亚洲欧美中文字幕日韩二区| 美女脱内裤让男人舔精品视频 | 一级毛片久久久久久久久女| 成人综合一区亚洲| 99热网站在线观看| 国产美女午夜福利| 日本成人三级电影网站| 狠狠狠狠99中文字幕| 亚洲人成网站在线观看播放| 国产成人精品婷婷| 丝袜喷水一区| 草草在线视频免费看| 中文字幕制服av| 美女大奶头视频| 国产视频内射| 亚洲av中文av极速乱| 在线免费十八禁| 国产成人精品久久久久久| 18禁裸乳无遮挡免费网站照片| 亚洲经典国产精华液单| av视频在线观看入口| 亚洲在久久综合| 最近视频中文字幕2019在线8| 午夜视频国产福利| 国产亚洲91精品色在线| 观看免费一级毛片| 久久国产乱子免费精品| 麻豆乱淫一区二区| 午夜福利在线观看免费完整高清在 | 欧美3d第一页| 国产精品伦人一区二区| 国产精品一二三区在线看| 成人毛片60女人毛片免费| 午夜精品在线福利| 免费在线观看成人毛片| 国产美女午夜福利| 亚洲精品乱码久久久久久按摩| 韩国av在线不卡| 久久人人精品亚洲av| 又爽又黄无遮挡网站| 亚洲成人精品中文字幕电影| 亚洲欧美精品综合久久99| 少妇被粗大猛烈的视频| 国产91av在线免费观看| 搡女人真爽免费视频火全软件| 免费黄网站久久成人精品| 久久精品久久久久久噜噜老黄 | 色视频www国产| 欧美日韩在线观看h| 女的被弄到高潮叫床怎么办| 国产伦精品一区二区三区视频9| 老师上课跳d突然被开到最大视频| 美女黄网站色视频| 国产精品电影一区二区三区| 欧美性猛交黑人性爽| 国内精品久久久久精免费| 亚洲综合色惰| 久久精品久久久久久久性| .国产精品久久| 国内揄拍国产精品人妻在线| 久久久久网色| 在线a可以看的网站| 欧美xxxx黑人xx丫x性爽| 亚洲中文字幕日韩| 能在线免费看毛片的网站| 国产免费一级a男人的天堂| 精品无人区乱码1区二区| 在线观看av片永久免费下载| 日本黄色片子视频| 久久精品国产自在天天线| 国产精品永久免费网站| 男女做爰动态图高潮gif福利片| 又爽又黄a免费视频| 我的女老师完整版在线观看| 精品久久久久久久久久久久久| 久久这里只有精品中国| 日韩欧美国产在线观看| 99久久中文字幕三级久久日本| 国产在线精品亚洲第一网站| 国产国拍精品亚洲av在线观看| 国产精品一区二区三区四区免费观看| 最近视频中文字幕2019在线8| 日韩成人伦理影院| 搡老妇女老女人老熟妇| 国产精品一区二区三区四区免费观看| 亚洲高清免费不卡视频| 免费av不卡在线播放| 国产高清有码在线观看视频| 网址你懂的国产日韩在线| 性色avwww在线观看| 人体艺术视频欧美日本| 午夜免费激情av| 国语自产精品视频在线第100页| 插阴视频在线观看视频| 成人特级黄色片久久久久久久| 久久韩国三级中文字幕| av天堂中文字幕网| 久久久久九九精品影院| 亚洲中文字幕日韩| 禁无遮挡网站| 亚洲成人av在线免费| 久久久久九九精品影院| 人妻夜夜爽99麻豆av| 可以在线观看的亚洲视频| 2022亚洲国产成人精品| 国产精品久久久久久久久免| 国产高清激情床上av| 久久人妻av系列| 国产一区二区激情短视频| 91狼人影院| 99久久无色码亚洲精品果冻| 永久网站在线| 精品一区二区三区视频在线| 国产三级中文精品| 干丝袜人妻中文字幕| 看非洲黑人一级黄片| 亚洲av不卡在线观看| 国产不卡一卡二| 国产v大片淫在线免费观看| avwww免费| 最近最新中文字幕大全电影3| 又粗又硬又长又爽又黄的视频 | 国产免费一级a男人的天堂| 乱系列少妇在线播放| 2021天堂中文幕一二区在线观| 嘟嘟电影网在线观看| 亚洲欧洲国产日韩| 一级毛片aaaaaa免费看小| 少妇猛男粗大的猛烈进出视频 | 亚洲av不卡在线观看| 男女边吃奶边做爰视频| 亚洲成人av在线免费| 啦啦啦啦在线视频资源| 哪里可以看免费的av片| 亚洲av二区三区四区| 欧美色欧美亚洲另类二区| 久久中文看片网| 一个人看视频在线观看www免费| 国产精品久久久久久精品电影小说 | 丰满的人妻完整版| 国产伦一二天堂av在线观看| 精品一区二区三区视频在线| 亚洲美女视频黄频| 尾随美女入室| 国产淫片久久久久久久久| 国产精品av视频在线免费观看| 中文资源天堂在线| 欧美成人a在线观看| 大型黄色视频在线免费观看| 中出人妻视频一区二区| 91麻豆精品激情在线观看国产| 国产探花在线观看一区二区| а√天堂www在线а√下载| 亚洲第一区二区三区不卡| 欧美丝袜亚洲另类| 日日撸夜夜添| 亚洲国产日韩欧美精品在线观看| 一级毛片我不卡| 午夜福利高清视频| 国产精品久久电影中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产激情偷乱视频一区二区| 亚洲国产精品成人综合色| 22中文网久久字幕| 成年av动漫网址| 国产一区二区激情短视频| 成年版毛片免费区| 日韩 亚洲 欧美在线| 久久婷婷人人爽人人干人人爱| 在线国产一区二区在线| 日韩欧美精品免费久久|