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

    二維平行平板內(nèi)分散氣泡對液相表觀黏度影響的數(shù)值研究

    2016-10-22 07:14:23徐磊龐明軍
    化工進展 2016年10期
    關(guān)鍵詞:毛細表觀氣泡

    徐磊,龐明軍

    (常州大學機械工程學院,江蘇 常州 213164)

    二維平行平板內(nèi)分散氣泡對液相表觀黏度影響的數(shù)值研究

    徐磊,龐明軍

    (常州大學機械工程學院,江蘇 常州 213164)

    泡狀液體廣泛存在于各種工業(yè)領(lǐng)域,因此深入了解泡狀液體的物理特性對產(chǎn)品質(zhì)量改進和過程強化具有重要意義。本文為了理解氣泡對液相表觀黏度的影響,揭示液相表觀黏度變化的原因,利用流體體積函數(shù)(VOF)與動網(wǎng)格結(jié)合的方法對二維平行板模型內(nèi)氣泡對液相表觀黏度的影響進行了深入的數(shù)值研究;分別研究了相同體積分數(shù)、不同毛細數(shù)工況下,氣泡加入對液相表觀黏度的影響。結(jié)果表明:毛細數(shù)和氣泡體積分數(shù)對液相的表觀黏度具有巨大的影響作用;相同的體積分數(shù)下,當毛細數(shù)較大(Ca>1)時,氣泡的加入會使液相的相對黏度減小,體積分數(shù)越大,液相的相對黏度減小得越嚴重;另外,相同的體積分數(shù)下,當毛細數(shù)較小(Ca<1)時,氣泡的加入會使液相的相對黏度增大,體積分數(shù)越高,液相的相對黏度增大得越明顯。

    泡狀液體;分散氣泡;相對黏度;流體體積函數(shù);動網(wǎng)格;二維平行平板

    泡狀液體廣泛存在于化工、食品、流動輸送、生物醫(yī)學等工業(yè)領(lǐng)域。已有的研究表明,氣泡的加入會極大地影響工業(yè)過程的質(zhì)量、熱量和動量傳遞效率,甚至會對泡狀設施的運行效率和產(chǎn)品質(zhì)量帶來巨大的影響[1-11]。為了有效地設計泡狀設施和高效地控制泡狀設施的運行過程,深入理解泡狀液體的水動力學特性是非常必要的。特別是理解氣泡加入對液相流變特性的影響,對于理解氣泡湍流減阻機理和改進食品工業(yè)的產(chǎn)品質(zhì)量具有重要的意義。

    目前,國外研究者針對泡狀液體的流變特性已開展了一定的實驗和數(shù)值研究。RUST和 MICHAEL[1]基于旋轉(zhuǎn)流變儀原理測量了氣泡對玉米糖漿黏度的影響。YUICHI和HIROSHI[2]利用落球黏度計原理詳細測試了氣泡瞬態(tài)變形對硅油黏度的影響。LLEWELLIN等[3]學者在實驗研究的基礎上,通過理論分析給出一組泡狀流流變特性的本構(gòu)模型。文獻[4-8]報道了氣泡對巖漿流變特性影響的研究。文獻[9-10]給出了低毛細數(shù)下,泡狀流相對黏度的本構(gòu)方程。文獻[7,11]運用數(shù)值方法對泡狀流的流變特性開展了研究。

    綜上所述,發(fā)現(xiàn)盡管針對泡狀液體的流變特性開展了一定的研究,但限于問題的復雜性,仍未獲得 一個很好的機理解釋,不同的學者給出不同的結(jié)果。為了進一步理解氣泡存在時是如何影響液相的表 觀黏度,并希望獲得一些液相表觀黏度變化的原因,本文采用二維平行板模型,對分散氣泡對液相表觀黏度的影響進行了數(shù)值模擬。

    1 幾何模型和計算條件

    1.1 幾何模型、物性參數(shù)與邊界條件

    鑒于二維平行平板流動(即庫埃流動模型)簡單,在低雷諾數(shù)下,能給出理論解,且結(jié)果便于分析和對比,為此選用二維平行平板作為計算的幾何模型,如圖1所示。其中平板長度為L=1m,板間距為h=0.5m;計算時,將氣泡置于二維平行平板中心,氣泡半徑為r=0.1m、0.13m、0.15m。為了減小計算工作量,減小計算區(qū)域,在左右邊界設置為周期邊界條件;讓下平板靜止不動,而上平板以恒定的速度u=0.1m/s向右作勻速運動。為了便于計算收斂,將氣液相的密度均設置為ρ=0.1kg/m3,液相的動力黏度設置為μl=1000Pa·s,氣相的動力黏度設置為μb=100Pa·s。

    1.2 網(wǎng)格劃分

    對于流體體積函數(shù)(VOF)方法而言,選用非結(jié)構(gòu)化網(wǎng)格,計算精度較高、便于界面重構(gòu);且非結(jié)構(gòu)化網(wǎng)格由于網(wǎng)格單元形狀的多樣性以及單元之間連接關(guān)系的自由性、適應強,為此本文選用三角形網(wǎng)格進行網(wǎng)格劃分。其中,網(wǎng)格節(jié)點數(shù)為2145,網(wǎng)格單元數(shù)有4096,網(wǎng)格具體劃分情況見圖2。

    圖1 計算幾何模型

    圖2 網(wǎng)格劃分

    1.3 計算工況

    目前的研究表明,氣泡影響液相表觀黏度的因素可概括為毛細數(shù)(Ca)和體積分數(shù)(Φ)兩個量綱歸一化參數(shù)。為此,根據(jù)目前的計算模型和計算能力,依據(jù)這兩個參數(shù)設計了如下計算工況,詳見表 1。體積分數(shù)的改變是通過調(diào)整氣泡半徑的大小實現(xiàn)的,毛細數(shù)的大小,本文是通過改變表面張力系數(shù)實現(xiàn)的。計算時,為了便于計算收斂(即便于界面重構(gòu)),氣泡半徑選取得較大。

    其中毛細數(shù)是衡量氣泡變形大小的一個參數(shù)[1-5],具體指的是氣泡所受剪切力與表面張力的相對大小。計算表達式詳見式(1)。

    式中,μ0為純液體的動力黏度,Pa·s;為剪切率,s-1;r為沒有變形時氣泡的半徑,m;表示的是液相的表面張力系數(shù),N/m。

    表1 計算工況(含對應工況下的上平板所受剪應力)

    2 控制方程和計算方法

    2.1 控制方程

    因目前氣液兩相的運動速度較小,所以計算時均不考慮其可壓縮性。氣液兩相混合流體的控制方程表述如式(2)、式(3)。

    式中,ρ是氣液混合相的密度,kg/m3;p為壓強,Pa;t是時間,s;u為速度矢量,m/s;μ是黏度,Pa·s;F表示的是表面張力,N。

    2.2 VOF方法

    精確描述氣液兩相的運動界面是研究氣泡行為的關(guān)鍵。由于 VOF 方法追蹤的是網(wǎng)格中的流體體積,不是追蹤流體質(zhì)點的運動,因而具有計算速度較快和界面追蹤精度高等優(yōu)點[12-16]。因此,本文首先考慮采用VOF方法追蹤氣泡-流體界面。VOF方法的核心是相函數(shù)的輸運方程,所有的界面構(gòu)造和流體輸運都基于此方法進行計算,其表達式為

    式(4)。

    式中,u為速度,m/s;f為定義在整個計算區(qū)域內(nèi)的相函數(shù)。

    公式(4)可以看作f的守恒輸運公式,以二維流場為例,此時若將流場離散為網(wǎng)格單元,將f值在單個網(wǎng)格面積上求積分,該積分值與當前網(wǎng)格單元面積之比即為流體的局部體積分數(shù)(ψ),即式(6)。

    式中,i是網(wǎng)格單元編號;Ai是編號為i的網(wǎng)格單元面積。該公式描述了網(wǎng)格局部體積分數(shù)ψ與相函數(shù)f之間的關(guān)系,其中f是流體計算的對象,而ψ是界面追蹤的直觀指標以及界面構(gòu)造方案的關(guān)鍵。

    當計算單元中是液體和氣體時,方程中的流體物性為相應液體或氣體的物性。當計算單元內(nèi)包含兩相界面時,流體物性按照兩相體積分數(shù)的加權(quán)平均進行計算,即

    2.3 表面張力

    氣泡對液相的影響是通過表面張力實現(xiàn)的,要求解動量方法,必須計算出動量方程(3)中的表面張力項。

    本文采用連續(xù)表面張力模型(即CSF模型)來計算作用在氣液兩相界面的上力,其計算表達式為式(8)。

    式中,k為流體自由面的局部曲率;δ為狄拉克分布函數(shù);n為自由面的局部單位法向量。

    2.4 界面重構(gòu)

    界面重構(gòu)是VOF模型最關(guān)鍵的部分,這里應用分段線性界面計算方法(PLIC)[13],對氣液兩相界面進行重構(gòu)。該界面重構(gòu)方法的特點就是在單個網(wǎng)格內(nèi)用直線段近似氣液兩相的界面。首先確定界面的法向量n。界面的法向?qū)嶋H上就是相函數(shù)f的梯度方向,可以通過求解網(wǎng)格單元中氣泡局部體積分數(shù)(ψ)的梯度獲得[15]。界面法向單位法向量的計算公式如式(9)。

    然后根據(jù)法向和當前網(wǎng)格中的氣泡相的局部體積分數(shù)計算出當前氣液界面線的確切位置。最后將所有網(wǎng)格的氣泡相的局部體積分數(shù)加上流入的流體再減去流出的流體就是每個網(wǎng)格中下一時刻新的氣泡相的局部體積分數(shù),依據(jù)新的氣泡相的局部體積分數(shù)就可以重新對整個流場的自由面進行界面重構(gòu)。

    2.5 數(shù)值求解過程

    對于目前的計算,因上平板作勻速運動,所以采用動網(wǎng)格進模型進行計算。選用層流模型進行計算,求解時速度與壓力場的耦合采用SIMPLE算法;單元中心的變量梯度采用Least-Squares cell-Based,壓力的離散采用 PESTO格式,動量方程采用二階迎風格式進行離散,體積分數(shù)采用Geo-Reconstruct格式進行離散;時間步長Δt=0.001,時間的推進采用一階隱式格式。

    3 數(shù)值模擬結(jié)果分析

    為了驗證計算結(jié)果的準確性,先對單相流動進行了模擬,將模擬結(jié)果與理論結(jié)果進行了對比,詳見圖 4??梢钥闯觯壳暗挠嬎憬Y(jié)果與理論結(jié)果吻合的很好,表明計算結(jié)果是可信的。為了觀察計算結(jié)果是否取得統(tǒng)計穩(wěn)態(tài),計算兩相流動時,考察了上平板上的剪切力隨計算時間的變化情況,如圖5所示??梢钥闯觯?000次以后,上平板上的剪切力基本趨于穩(wěn)定,表明計算結(jié)果可以進行統(tǒng)計處理。本文選取迭代5000次以后的值進行平均,從而得到上平板的剪切力。然后再利用牛頓流體剪應力的計算公式反推得到液相的表觀黏度。

    圖4 單相流的速度計算曲線

    圖5 監(jiān)控體積分數(shù)為Φ=0.063與毛細數(shù)Ca=10的上平板上的剪切力

    3.1 氣泡加入對液相表觀黏度的影響

    為了研究氣泡加入對液相表觀黏度的影響,圖6給出不同工況下液相表觀黏度的變化情況。其中橫坐標為毛細數(shù),縱坐標為液相的相對黏度(為加入氣泡后液相的表觀黏度除以純液體對應狀態(tài)下的黏度)。為了對比,圖6也給出了文獻[8]中的實驗結(jié)果,可以看出目前的計算結(jié)果與文獻[8]的實驗結(jié)果定性一致,定量的差別或許是因目前計算時給出的假設引起的。

    圖6 相對黏度變化情況

    從圖6可以看出,當毛細數(shù)小于0.5時(Ca<0.5),隨著毛細數(shù)的減小,氣泡加入使液相的相對黏度增大,當毛細數(shù)小到一定程度時(對于目前的計算Ca≈0.2),液相相對黏度趨于穩(wěn)定、不再隨著毛細數(shù)的減小而再增大;然而,當毛細數(shù)大于 0.5時(Ca>0.5),隨著毛細數(shù)的增大,氣泡加入使液相的相對黏度減小,當毛細數(shù)大到一定程度時(對于目前的計算Ca≈10),液相相對黏度趨于穩(wěn)定、不再隨著毛細數(shù)的增大而再減小。而且在相同的毛細數(shù)下,氣泡的體積分數(shù)越大,氣泡對液相相對黏度的影響就越大。

    由此看來,氣泡對液相表觀黏度的影響與氣泡的體積分數(shù)和毛細數(shù)有直接的關(guān)系。而毛細數(shù)的大小又直接決定了氣泡形變的大小,所以氣泡形狀的改變又會對液相相對黏度的變化產(chǎn)生直接的影響。為此,選取了最大體積分數(shù)(Φ=0.14)下氣泡形狀隨毛細數(shù)的變化情況,詳見圖7。圖7為計算時間為5s時的氣液相相分布云圖,其中藍色區(qū)域表示的是氣泡,紅色為液相所占區(qū)域??梢钥闯觯S著毛細數(shù)的增大,氣泡由圓形逐漸變?yōu)榧氶L的橢圓形,且橢圓形氣泡的長軸與平板法向的夾角與最大剪切率與法向的夾角基本一致。

    當毛細數(shù)較小時(對于目前的計算 Ca<0.2),氣泡所受的剪切力小于其所受的表面張力,氣泡保持為圓形,這時液相的相對黏度由于氣泡的加入而增大。當毛細數(shù)進一步減小時,氣泡的形狀不會發(fā)生變化、一直保持為球形,所以對應的液相的相對黏度也不會再增大、而趨于一個定值。然而,當毛細數(shù)較大時(對于目前的計算Ca>10),氣泡所受的剪切力大于表面張力,由于強的剪切作用,氣泡不能再保持球形而變?yōu)榧氶L的橢圓形,這時液相的表觀黏度由于氣泡的加入而減小。因這時氣泡長軸與平板法向的夾角已近似等于最大剪切率與法向的夾角,當毛細數(shù)再增大,橢圓形氣泡長軸與法向的夾角也不會變化,所以液相的相對黏度也趨于一個定值、不會隨著毛細數(shù)的增大而減小。

    圖7 氣泡體積分數(shù)Φ=0.14在5s時氣液相分布云圖

    3.2 氣泡加入對液相流場的影響

    為了理解毛細數(shù)(即氣泡形狀)對液相表觀黏度的影響,下面分析了氣泡加入對液相流場的影響情況。圖8給出最大體積分數(shù)為Φ=0.14下毛細管為0.1和1時的壓力場分布,圖9給出速度矢量場和渦量場的分布情況。可以看出,當毛細數(shù)為 0.1時,由于較大的表面張力,氣泡幾乎保持為球形,氣泡周圍的壓力小于氣泡內(nèi)部的壓力,在氣泡周圍最大剪切方向上產(chǎn)生大面積的負壓區(qū),對應的在負壓區(qū)域出現(xiàn)大面積的渦量區(qū)(如圖9所示),導致氣泡與液體在此區(qū)域的黏性耗散增加,所以宏觀上表現(xiàn)為液相相對黏度的增加。

    圖8 氣泡體積分數(shù)Φ=0.14的壓力場

    當毛細數(shù)為1時,由于剪切力大于表面張力,氣泡發(fā)生嚴重的變形,在流體流動方向[如圖 8(b)中箭頭所示]上,氣泡的左右兩側(cè)出現(xiàn)了小面積的負壓區(qū),對應的在此區(qū)域也出現(xiàn)了一定的渦量區(qū);負壓區(qū)沒有出現(xiàn)在最大剪切方向上,所以在該方向上沒有渦運動(如圖9所示),即氣泡形變提供了自由滑移面,所以氣泡加入使黏性耗散增加的不明顯。反而,由于氣泡形變后沿著最大剪切方向呈現(xiàn)為細長橢圓形,而氣體的黏度又小于液體的黏度,即低黏度流體占據(jù)了大面積的高剪切率區(qū)域,所以形變的氣泡導致黏性耗散大幅度減小??傊?,對于形變的氣泡而言,其黏性耗散減小的程度大于黏性耗散增加的程度,最終宏觀表現(xiàn)為液相的相對黏度減小。

    為了直觀得觀察不同毛細數(shù)下,氣泡對液相流場的影響情況,圖10給出了不同毛細數(shù)下流場流線分布情況??梢钥闯?,與單液體相相比,毛細數(shù)較小時[如圖10(b)所示],氣泡形變較?。ㄒ簿褪钦f氣泡幾乎保持為圓形),氣泡加入,使流線彎曲的越厲害,導致黏性耗散增加,表現(xiàn)為液相相對黏度增大;反之,毛細數(shù)較大時[如圖10(d)所示],氣泡形變較大、為流體的流動提供了滑移面,氣泡加入使流線的彎曲不明顯,導致總體黏性耗散減小,表現(xiàn)為液相相對黏度減小。這與文獻[1,5]給出的解釋相符。

    圖9 氣泡體積分數(shù)Φ=0.14的速度場和渦量場

    圖10 不同毛細數(shù)下流線分布

    4 結(jié) 論

    本文通過對二維平行板模型內(nèi)分散氣泡對液相表觀黏度的影響進行了詳細的數(shù)值模擬,可以得到如下幾點結(jié)論。

    (1)當毛細數(shù)小于0.5時,氣泡所受的剪切力小于表面張力,氣泡幾乎保持為圓形,氣泡的加入使液相的相對黏度增大。當毛細數(shù)進一步減小時,由于氣泡的形狀不會發(fā)生變化、一直保持為球形,所以對應液相的相對黏度也不會再增大,而趨于一個定值。而且在相同的毛細數(shù)下,氣泡的體積分數(shù)越大,氣泡對液相相對黏度增加的越大。

    (2)當毛細數(shù)大于0.5,氣泡所受的剪切力大于表面張力,在剪切作用下,氣泡變?yōu)榧氶L的橢圓形,氣泡的加入導致液相表觀黏度的減小。這時氣泡長軸與平板法向的夾角已近似等于最大剪切率與法向的夾角,當毛細數(shù)再增大,橢圓形氣泡長軸與法向的夾角也不會變化,所以液相的相對黏度也趨于一個定值,不會隨著毛細數(shù)的增大而減小。然而隨著體積分數(shù)的增加,液相表觀黏度的減小幅度增大。

    (3)在低毛細數(shù)下,由于較大的表面張力,氣泡幾乎保持為球形,氣泡周圍的壓力小于氣泡內(nèi)部的壓力,在氣泡周圍最大剪切方向上產(chǎn)生大面積的負壓區(qū),對應的在負壓區(qū)域出現(xiàn)大面積的渦量區(qū),并且氣泡的加入使流線彎曲嚴重,導致黏性耗散增加,所以表現(xiàn)為液相相對黏度增大。

    (4)在高毛細數(shù)下,由于氣泡所受的剪切力大于表面張力,氣泡發(fā)生嚴重的變形,在流體流動方向上,氣泡的左右兩側(cè)出現(xiàn)了小面積的負壓區(qū),對應的在此區(qū)域也出現(xiàn)了一定的渦量區(qū);負壓區(qū)沒有出現(xiàn)在最大剪切方向上,所以在該方向上沒有渦運動,即氣泡形變提供了自由滑移面。所以形變的氣泡導致黏性耗散大幅度減小。表現(xiàn)為液相相對黏度減小。

    [1] RUST A C,MICHAEL M. Effects of bubble deformation on the viscosity of dilute suspensions[J]. J. Non-Newtonian Fluid Mech., 2002,104:53-63.

    [2] YUICHI M,HIROSHI O. Increase of effective viscosity in bubbly liquids from transient bubble deformation[J]. Fluid Dynamics Research,2008,40:565-575.

    [3] LLEWELLIN E W,MADER H M,WILSON S D R. The rheology of a bubbly liquid[J]. Proc. R. Soc. A,2002,458:987-1016.

    [4] MANGA M,LOEWENBERG M. Viscosity of magmas containing highly deformable bubbles[J]. Journal of Volcanology and Geothermal Research,2001,105:19-24.

    [5] MADER H M,LLEWELLIN E W,MUELLER S P. The rheology of two-phase magmas:a review and analysis[J]. Journal of Volcanology and Geothermal Research,2013,257:135-158.

    [6] LEJEUNE A M,BOTTINGA Y,TRULL T W,et al. Rheology of bubble-baring magmas:reply to a critical comment by Spera and Stein[J]. Earth and Planetary Science Letters,2000,175:333-334.

    [7] MANGA M,CASTRO J,CASHMAN K V,et al. Rheology of bubble-bearing magmas[J]. J. Volcanol. Geotherm. Res.,1998,87:15-28.

    [8] SUNG W J,SEUNG H L,JAE R Y. Rheological behavior of polydispersed bubble suspensions in shear flows[J]. Polymer Engineering and Science,2010,50(1):128-137

    [9] PAL R. rheological constitutive equation for bubbly suspensions[J]. Ind. Eng. Chem. Res.,2004,43:5372-5379.

    [10] FRANKEL N A,ACRIVOS A. The constitutive equation for a dilute emulsion[J]. J. Fluid Mech.,1970,44:65-78.

    [11] YUN M L,DONG J S,JAE R Y. Rheological behavior of dilute bubble suspensions in polyol[J]. Korea-Australia Rheology Journal,2004,16:47-54.

    [12] HIRT C W,NICHOLS B D. Volume of fluid(VOF)method for the dynamics of free boundaries[J]. Journal of Computational Physics, 1981,39:201-205.

    [13] RIDER W J,KOTHE D B. Reconstructing volume tracking[J]. Journal of Computational Physics,1997,141:112-152.

    [14] 馬斗. 單孔氣泡動力學行為的 VOF 數(shù)值模擬[D]. 天津:天津大學化工學院,2009.

    [15] 黃萌. 電場作用下油水乳液中水滴的運動和形變特性研究[D]. 西安:西安交通大學,2015.

    [16] 李宜樂,劉應中. 繆國平. 二階精度的VOF自由面追蹤方法及其應用[J]. 船舶力學,1999,3(1):44-52.

    Numerical study on influence of dispersed bubbles on liquid-phase apparent viscosity in two-dimensional parallel plate

    XU Lei,PANG Mingjun
    (School of Mechanical Engineering,Changzhou University,Changzhou 213164,Jiangsu,China)

    Bubble liquid exists widely in various industrial fields,so the detailed understanding of physical properties of bubbly liquid is of great significance for improving the quality of products and for strengthening industrial processes. In this paper,in order to understand the influence of bubbles on the liquid-phase apparent viscosity and to explain the changing reason of the liquid-phase apparent viscosity,based on two-dimensional parallel plate model,the effect of large bubbles on the liquid-phase apparent viscosity was deeply investigated with the volume of fluid (VOF) method combined with the dynamic mesh. The effect of bubbles on the liquid-phase apparent viscosity was numerically studied in detail under the conditions with the same void fraction but different capillary numbers. The present study showed that the capillary number and the bubble volume fraction has the great effect on liquid-phase apparent viscosity. For the same volume fraction,when the capillary number is relatively large (Ca>1),the bubble injection causes the decrease of the relative viscosity of the liquid phase,and the larger the volume fraction is. The relative viscosity of the liquid phase decreases more sharper. In addition,under the same volume fraction,when the capillary number is relatively small (Ca<1),the relative viscosity of the liquid phase increases due to the addition of bubbles,and the higher the volume fraction is. The relative viscosity of the liquid phase increases more obviously.

    bubbly liquid; dispersed bubble;relative viscosity;volume of fluid(VOF);moving grid; two-dimensional parallel plate

    O 37

    A

    1000-6613(2016)10-3065-07

    10.16085/j.issn.1000-6613.2016.10.008

    2016-03-03;修改稿日期:2016-03-31。

    國家自然科學基金(51376026)及江蘇省青藍工程項目。

    徐磊 (1990—),碩士研究生,研究方向為泡狀流流變特性。聯(lián)系人:龐明軍,博士,副教授,研究方向為氣液兩相流動與湍流減阻。E-mail pangmj@cczu.edu.cn。

    猜你喜歡
    毛細表觀氣泡
    “毛細”政務號關(guān)停背后
    廉政瞭望(2024年5期)2024-05-26 13:21:07
    檸檬氣泡水
    欣漾(2024年2期)2024-04-27 15:19:49
    綠盲蝽為害與赤霞珠葡萄防御互作中的表觀響應
    河北果樹(2021年4期)2021-12-02 01:14:50
    SIAU詩杭便攜式氣泡水杯
    新潮電子(2021年7期)2021-08-14 15:53:12
    浮法玻璃氣泡的預防和控制對策
    多孔建筑材料毛細吸水系數(shù)連續(xù)測量的方法
    能源工程(2020年6期)2021-01-26 00:55:16
    鋼結(jié)構(gòu)表觀裂紋監(jiān)測技術(shù)對比與展望
    上海公路(2019年3期)2019-11-25 07:39:28
    例析對高中表觀遺傳學的認識
    冰凍氣泡
    出現(xiàn)憋喘 可能是毛細支氣管炎!
    媽媽寶寶(2017年3期)2017-02-21 01:22:16
    精品亚洲乱码少妇综合久久| 少妇被粗大猛烈的视频| 高清视频免费观看一区二区| kizo精华| 久热这里只有精品99| 精品亚洲成a人片在线观看| 中文字幕久久专区| 岛国毛片在线播放| 免费高清在线观看视频在线观看| 成人毛片60女人毛片免费| 国产白丝娇喘喷水9色精品| 国产精品久久久久久久久免| 精品卡一卡二卡四卡免费| a级毛片免费高清观看在线播放| 国产精品 国内视频| a级片在线免费高清观看视频| 新久久久久国产一级毛片| 国产精品久久久久久久久免| 水蜜桃什么品种好| 99久久综合免费| 青春草亚洲视频在线观看| 18禁裸乳无遮挡动漫免费视频| 五月伊人婷婷丁香| 国产极品天堂在线| 大话2 男鬼变身卡| 日本色播在线视频| 欧美国产精品一级二级三级| 亚洲欧美日韩卡通动漫| 在线免费观看不下载黄p国产| 国产欧美日韩一区二区三区在线 | 永久网站在线| 少妇丰满av| 男女边吃奶边做爰视频| 性色avwww在线观看| 狂野欧美激情性bbbbbb| 九九爱精品视频在线观看| 精品一区二区三卡| 一个人看视频在线观看www免费| 全区人妻精品视频| 日韩精品有码人妻一区| 国产精品欧美亚洲77777| 国产视频首页在线观看| 热99国产精品久久久久久7| 性色av一级| 亚洲av在线观看美女高潮| 男女无遮挡免费网站观看| 在线观看免费视频网站a站| 国产成人freesex在线| 肉色欧美久久久久久久蜜桃| 观看av在线不卡| 在线观看三级黄色| 精品卡一卡二卡四卡免费| 亚洲精品第二区| 99热这里只有精品一区| 赤兔流量卡办理| 夫妻性生交免费视频一级片| 2022亚洲国产成人精品| 午夜久久久在线观看| 人人澡人人妻人| videossex国产| 男人爽女人下面视频在线观看| 亚洲国产欧美日韩在线播放| 少妇被粗大的猛进出69影院 | 在线观看三级黄色| 一级毛片我不卡| 久久99热6这里只有精品| 精品亚洲乱码少妇综合久久| 在线精品无人区一区二区三| 成人毛片a级毛片在线播放| 婷婷色综合大香蕉| 国产亚洲最大av| 春色校园在线视频观看| 亚洲精品乱码久久久v下载方式| 精品一区二区三卡| 老司机影院毛片| 汤姆久久久久久久影院中文字幕| 91aial.com中文字幕在线观看| 亚州av有码| xxxhd国产人妻xxx| 久久久久国产精品人妻一区二区| 大又大粗又爽又黄少妇毛片口| 岛国毛片在线播放| 国产精品一国产av| 99九九在线精品视频| 3wmmmm亚洲av在线观看| 欧美日韩国产mv在线观看视频| 九色成人免费人妻av| 亚洲精品日本国产第一区| 成年av动漫网址| 性高湖久久久久久久久免费观看| 七月丁香在线播放| 欧美精品亚洲一区二区| 十八禁网站网址无遮挡| 久久久久久伊人网av| 日韩免费高清中文字幕av| 国产精品麻豆人妻色哟哟久久| 三级国产精品欧美在线观看| 国产精品欧美亚洲77777| 亚洲五月色婷婷综合| 国产精品不卡视频一区二区| 亚洲国产av新网站| 99九九在线精品视频| 亚洲av成人精品一区久久| 我的老师免费观看完整版| 有码 亚洲区| 久久人人爽人人片av| 男人添女人高潮全过程视频| 国产精品一区二区在线观看99| 国产一区二区在线观看av| 国产国拍精品亚洲av在线观看| 一本久久精品| 久久久欧美国产精品| 91国产中文字幕| 91在线精品国自产拍蜜月| 久久精品国产亚洲av天美| 女性生殖器流出的白浆| 国产永久视频网站| 最新的欧美精品一区二区| 国产伦精品一区二区三区视频9| 亚洲欧美成人精品一区二区| 欧美激情国产日韩精品一区| 亚洲人成77777在线视频| 午夜福利网站1000一区二区三区| 免费播放大片免费观看视频在线观看| 青春草亚洲视频在线观看| 热99国产精品久久久久久7| 99热这里只有精品一区| 不卡视频在线观看欧美| 精品久久久噜噜| 成人毛片60女人毛片免费| 丰满乱子伦码专区| 国产爽快片一区二区三区| 波野结衣二区三区在线| 欧美另类一区| 免费黄网站久久成人精品| 一级毛片黄色毛片免费观看视频| 搡老乐熟女国产| 如何舔出高潮| 亚洲av日韩在线播放| 免费av中文字幕在线| 天堂俺去俺来也www色官网| 国产欧美亚洲国产| 国国产精品蜜臀av免费| 精品亚洲乱码少妇综合久久| 亚洲精品乱码久久久久久按摩| 国产精品女同一区二区软件| 老熟女久久久| 亚洲激情五月婷婷啪啪| 成人无遮挡网站| 成年女人在线观看亚洲视频| 中文欧美无线码| 午夜福利视频精品| 免费少妇av软件| 亚洲av国产av综合av卡| 国产免费视频播放在线视频| h视频一区二区三区| 久久99精品国语久久久| 精品久久久精品久久久| 99国产精品免费福利视频| 亚洲欧洲国产日韩| 最近的中文字幕免费完整| 免费大片黄手机在线观看| 久久久久久久久大av| 男人爽女人下面视频在线观看| 日日摸夜夜添夜夜添av毛片| 26uuu在线亚洲综合色| 婷婷色综合www| 久久久精品免费免费高清| 午夜av观看不卡| 亚洲成人手机| 大又大粗又爽又黄少妇毛片口| 成人手机av| 丝袜美足系列| 18+在线观看网站| 3wmmmm亚洲av在线观看| av不卡在线播放| 国产伦理片在线播放av一区| 日韩中文字幕视频在线看片| 嘟嘟电影网在线观看| 九九在线视频观看精品| 男女高潮啪啪啪动态图| 免费黄频网站在线观看国产| av在线观看视频网站免费| 一二三四中文在线观看免费高清| 亚洲丝袜综合中文字幕| 日本vs欧美在线观看视频| 亚洲国产色片| 国产高清不卡午夜福利| 成人18禁高潮啪啪吃奶动态图 | kizo精华| 精品久久蜜臀av无| 成年av动漫网址| 人成视频在线观看免费观看| 卡戴珊不雅视频在线播放| 欧美一级a爱片免费观看看| 欧美日韩视频精品一区| 九色亚洲精品在线播放| 欧美激情极品国产一区二区三区 | 国产成人免费观看mmmm| 欧美亚洲日本最大视频资源| 久久久久久久久久久久大奶| 国产免费一区二区三区四区乱码| 亚洲国产精品999| 少妇人妻久久综合中文| 欧美丝袜亚洲另类| 两个人的视频大全免费| 日韩欧美精品免费久久| 天天影视国产精品| 国产女主播在线喷水免费视频网站| 国产亚洲最大av| 亚洲av不卡在线观看| 夜夜看夜夜爽夜夜摸| 午夜激情久久久久久久| 午夜激情福利司机影院| 国产探花极品一区二区| 国产国语露脸激情在线看| 大片免费播放器 马上看| av在线老鸭窝| 国产永久视频网站| 久久青草综合色| av黄色大香蕉| 国产永久视频网站| 亚洲内射少妇av| 麻豆乱淫一区二区| 成人免费观看视频高清| 亚洲一级一片aⅴ在线观看| 成年人午夜在线观看视频| 赤兔流量卡办理| 久久99蜜桃精品久久| av天堂久久9| 亚洲国产成人一精品久久久| 免费看光身美女| 成人国产麻豆网| 婷婷色综合www| 欧美变态另类bdsm刘玥| 久久99精品国语久久久| 在线观看免费高清a一片| 内地一区二区视频在线| 精品酒店卫生间| 欧美日韩视频精品一区| 亚洲欧美日韩卡通动漫| 九九爱精品视频在线观看| 国产高清不卡午夜福利| 精品一区二区免费观看| 91精品国产国语对白视频| 国产在线一区二区三区精| 熟女电影av网| 在线观看三级黄色| 下体分泌物呈黄色| 热99久久久久精品小说推荐| 大香蕉久久成人网| 久久人人爽av亚洲精品天堂| 特大巨黑吊av在线直播| 两个人的视频大全免费| 国产精品久久久久成人av| 欧美激情国产日韩精品一区| 久久午夜综合久久蜜桃| 成人午夜精彩视频在线观看| 又粗又硬又长又爽又黄的视频| 美女视频免费永久观看网站| 一级,二级,三级黄色视频| 一本久久精品| 国产高清不卡午夜福利| 毛片一级片免费看久久久久| 99国产精品免费福利视频| 有码 亚洲区| 精品亚洲成a人片在线观看| 男人操女人黄网站| 日韩一区二区三区影片| 亚洲精品视频女| 国产在视频线精品| 色94色欧美一区二区| 国产精品女同一区二区软件| 国产精品国产av在线观看| 成人国产av品久久久| 国产成人精品在线电影| 国产免费视频播放在线视频| 人妻制服诱惑在线中文字幕| 日本爱情动作片www.在线观看| 少妇的逼水好多| 精品午夜福利在线看| 亚洲图色成人| 精品久久久久久久久亚洲| 99热全是精品| 久久久亚洲精品成人影院| 精品久久久精品久久久| 欧美日韩综合久久久久久| 久久97久久精品| 久久亚洲国产成人精品v| 精品午夜福利在线看| 久久久久久久久久久丰满| 69精品国产乱码久久久| 亚洲精品中文字幕在线视频| 日韩av不卡免费在线播放| 99久国产av精品国产电影| 免费高清在线观看日韩| 美女内射精品一级片tv| 少妇高潮的动态图| xxxhd国产人妻xxx| 免费看光身美女| 国产精品一区二区在线不卡| 精品久久久噜噜| 七月丁香在线播放| 国产亚洲av片在线观看秒播厂| 一区在线观看完整版| 最近的中文字幕免费完整| 国产免费福利视频在线观看| 春色校园在线视频观看| 中文字幕人妻丝袜制服| 王馨瑶露胸无遮挡在线观看| 国产高清不卡午夜福利| 99九九在线精品视频| 女人久久www免费人成看片| 观看美女的网站| 蜜桃在线观看..| 欧美精品国产亚洲| 美女中出高潮动态图| 久久久久久久久大av| 亚洲美女搞黄在线观看| 一区二区av电影网| 熟妇人妻不卡中文字幕| 日本wwww免费看| 亚洲精品国产av蜜桃| 一区二区三区精品91| 精品人妻熟女av久视频| 伊人亚洲综合成人网| 美女cb高潮喷水在线观看| 黑丝袜美女国产一区| 毛片一级片免费看久久久久| 免费av不卡在线播放| av在线播放精品| 精品少妇黑人巨大在线播放| 久久ye,这里只有精品| 久久久精品免费免费高清| 国产国拍精品亚洲av在线观看| 久久久久久久亚洲中文字幕| 欧美精品国产亚洲| 秋霞伦理黄片| 国产亚洲午夜精品一区二区久久| 夜夜爽夜夜爽视频| 国产亚洲欧美精品永久| 丝袜喷水一区| 一本—道久久a久久精品蜜桃钙片| 内地一区二区视频在线| 26uuu在线亚洲综合色| 丝袜在线中文字幕| 国产欧美日韩综合在线一区二区| 亚洲精品国产av蜜桃| 精品久久久久久久久亚洲| 999精品在线视频| 免费高清在线观看视频在线观看| 国产精品一二三区在线看| 日韩av在线免费看完整版不卡| av女优亚洲男人天堂| 我的老师免费观看完整版| 国产精品嫩草影院av在线观看| 熟女电影av网| 亚洲精品日本国产第一区| 日本黄大片高清| 免费观看的影片在线观看| 99国产综合亚洲精品| 日日摸夜夜添夜夜爱| 国产在线免费精品| 亚洲综合色惰| 国产在线一区二区三区精| 日韩制服骚丝袜av| 免费黄频网站在线观看国产| 国产男人的电影天堂91| 亚洲国产毛片av蜜桃av| 亚洲精品乱久久久久久| 亚洲欧美色中文字幕在线| 18禁裸乳无遮挡动漫免费视频| 国产亚洲精品第一综合不卡 | 亚洲国产精品一区三区| 99国产精品免费福利视频| 青青草视频在线视频观看| 欧美人与善性xxx| 国产精品女同一区二区软件| 国产成人一区二区在线| 天堂8中文在线网| 午夜福利,免费看| 青春草亚洲视频在线观看| 亚洲无线观看免费| 亚洲欧美日韩卡通动漫| 国产男女超爽视频在线观看| 亚洲四区av| 成年美女黄网站色视频大全免费 | 欧美日韩成人在线一区二区| 久久99热6这里只有精品| 中文欧美无线码| 美女视频免费永久观看网站| 精品国产一区二区久久| 又黄又爽又刺激的免费视频.| 中国三级夫妇交换| 极品少妇高潮喷水抽搐| 18禁裸乳无遮挡动漫免费视频| 少妇精品久久久久久久| 国产亚洲午夜精品一区二区久久| 日韩成人伦理影院| 校园人妻丝袜中文字幕| 波野结衣二区三区在线| 热99国产精品久久久久久7| 国产成人精品无人区| 中文乱码字字幕精品一区二区三区| 中文天堂在线官网| 天堂俺去俺来也www色官网| 欧美日韩视频精品一区| 乱人伦中国视频| 中国三级夫妇交换| 久久这里有精品视频免费| 五月玫瑰六月丁香| 日日啪夜夜爽| 亚洲欧洲国产日韩| 午夜影院在线不卡| 97超视频在线观看视频| 久久精品国产鲁丝片午夜精品| 视频区图区小说| 中文字幕人妻熟人妻熟丝袜美| 国产一区亚洲一区在线观看| 日韩一区二区视频免费看| 久久久久人妻精品一区果冻| 精品少妇黑人巨大在线播放| 亚洲婷婷狠狠爱综合网| 人妻夜夜爽99麻豆av| 少妇人妻精品综合一区二区| 国产熟女午夜一区二区三区 | 午夜激情久久久久久久| 欧美老熟妇乱子伦牲交| 最近中文字幕高清免费大全6| 大香蕉久久网| 春色校园在线视频观看| 亚洲激情五月婷婷啪啪| 如何舔出高潮| 最近最新中文字幕免费大全7| 国产精品麻豆人妻色哟哟久久| 九九久久精品国产亚洲av麻豆| 国产无遮挡羞羞视频在线观看| av播播在线观看一区| 少妇人妻精品综合一区二区| 寂寞人妻少妇视频99o| 能在线免费看毛片的网站| 免费看光身美女| 国产精品女同一区二区软件| 一区二区日韩欧美中文字幕 | 久久热精品热| 国产在视频线精品| 色94色欧美一区二区| 91在线精品国自产拍蜜月| 国产精品不卡视频一区二区| √禁漫天堂资源中文www| 精品一区二区三卡| 各种免费的搞黄视频| 内地一区二区视频在线| 哪个播放器可以免费观看大片| 少妇熟女欧美另类| 久久热精品热| 曰老女人黄片| 夜夜看夜夜爽夜夜摸| 久久午夜福利片| 免费大片黄手机在线观看| 国产成人av激情在线播放 | 成人国产麻豆网| 欧美成人午夜免费资源| 久久久国产一区二区| 免费人成在线观看视频色| 中文字幕最新亚洲高清| 国产白丝娇喘喷水9色精品| 大香蕉97超碰在线| 亚洲国产成人一精品久久久| 国产欧美日韩综合在线一区二区| av在线老鸭窝| 成人18禁高潮啪啪吃奶动态图 | √禁漫天堂资源中文www| 久久毛片免费看一区二区三区| 亚洲情色 制服丝袜| 精品少妇黑人巨大在线播放| 最近手机中文字幕大全| 国产极品天堂在线| 久久久久久久久久久丰满| 国产男女内射视频| 精品少妇黑人巨大在线播放| 夜夜骑夜夜射夜夜干| 亚洲av中文av极速乱| 久久久久久久国产电影| 一级片'在线观看视频| 日日摸夜夜添夜夜添av毛片| 欧美日韩国产mv在线观看视频| 人妻夜夜爽99麻豆av| 日本黄色片子视频| 精品一区二区免费观看| 亚洲精品久久成人aⅴ小说 | kizo精华| 好男人视频免费观看在线| av网站免费在线观看视频| 超碰97精品在线观看| 少妇人妻久久综合中文| 另类亚洲欧美激情| 99久久综合免费| 中文天堂在线官网| 午夜精品国产一区二区电影| 2018国产大陆天天弄谢| 看免费成人av毛片| 2021少妇久久久久久久久久久| av在线观看视频网站免费| 乱人伦中国视频| 夫妻午夜视频| 国产白丝娇喘喷水9色精品| 亚洲欧美中文字幕日韩二区| 国产 一区精品| 自线自在国产av| 精品久久国产蜜桃| 女人久久www免费人成看片| 亚洲av免费高清在线观看| 蜜臀久久99精品久久宅男| 久久久久久久久久久丰满| 丰满乱子伦码专区| 夫妻性生交免费视频一级片| 亚洲精品第二区| 亚洲怡红院男人天堂| 午夜视频国产福利| 久久97久久精品| 两个人的视频大全免费| 国产欧美日韩一区二区三区在线 | 免费看光身美女| 99久国产av精品国产电影| 国产一区二区三区av在线| 国产高清三级在线| 欧美性感艳星| 你懂的网址亚洲精品在线观看| 久久99热6这里只有精品| 色吧在线观看| 人妻一区二区av| 高清欧美精品videossex| 国产午夜精品久久久久久一区二区三区| 涩涩av久久男人的天堂| 黑人欧美特级aaaaaa片| 日本-黄色视频高清免费观看| 久久精品国产亚洲网站| 十分钟在线观看高清视频www| 国产伦理片在线播放av一区| 免费观看的影片在线观看| 日本av免费视频播放| tube8黄色片| 99九九线精品视频在线观看视频| 免费观看在线日韩| 免费看av在线观看网站| 亚洲精品色激情综合| 伦精品一区二区三区| 热99国产精品久久久久久7| 日本午夜av视频| 日韩精品免费视频一区二区三区 | 免费观看的影片在线观看| 如日韩欧美国产精品一区二区三区 | 蜜桃国产av成人99| 国产亚洲欧美精品永久| 少妇高潮的动态图| 51国产日韩欧美| 色婷婷av一区二区三区视频| 18+在线观看网站| 国产在视频线精品| 一区在线观看完整版| 国产淫语在线视频| 中文字幕av电影在线播放| 精品人妻偷拍中文字幕| 国产乱人偷精品视频| 国产亚洲最大av| 高清午夜精品一区二区三区| 亚洲精品456在线播放app| 2018国产大陆天天弄谢| 色94色欧美一区二区| 国产日韩欧美视频二区| 日韩中文字幕视频在线看片| 黄色配什么色好看| 国产精品一区www在线观看| 九九久久精品国产亚洲av麻豆| 日韩成人伦理影院| 免费日韩欧美在线观看| 免费播放大片免费观看视频在线观看| 国产又色又爽无遮挡免| 一本大道久久a久久精品| 69精品国产乱码久久久| 两个人免费观看高清视频| 丰满少妇做爰视频| 成人黄色视频免费在线看| 亚洲精品久久成人aⅴ小说 | 国产成人精品福利久久| 亚洲精品成人av观看孕妇| 国产一区二区在线观看av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美日韩视频精品一区| 亚洲丝袜综合中文字幕| 亚洲国产最新在线播放| 日本色播在线视频| 少妇的逼水好多| 免费观看的影片在线观看| 精品国产乱码久久久久久小说| 热99久久久久精品小说推荐| 我要看黄色一级片免费的| 日韩伦理黄色片| av视频免费观看在线观看| 中文字幕精品免费在线观看视频 | 嫩草影院入口| 黄片无遮挡物在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲婷婷狠狠爱综合网| 国产精品一国产av| 80岁老熟妇乱子伦牲交| 我要看黄色一级片免费的| 久久精品久久久久久噜噜老黄| 欧美老熟妇乱子伦牲交| 美女国产高潮福利片在线看| 亚洲综合色惰| 欧美激情 高清一区二区三区| 在线 av 中文字幕| 日本vs欧美在线观看视频|