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

    超聲波作用下液滴的冷卻凍結(jié)規(guī)律

    2017-11-22 10:20:08高蓬輝張夢杜玉吉程博張東海周國慶
    化工學(xué)報 2017年11期
    關(guān)鍵詞:傳質(zhì)液滴氣泡

    高蓬輝,張夢,杜玉吉,2,程博,張東海,周國慶

    (1中國礦業(yè)大學(xué)力學(xué)與土木工程學(xué)院,江蘇 徐州 221116;2中節(jié)能城市節(jié)能研究院有限公司,江蘇 常州 213001)

    超聲波作用下液滴的冷卻凍結(jié)規(guī)律

    高蓬輝1,張夢1,杜玉吉1,2,程博1,張東海1,周國慶1

    (1中國礦業(yè)大學(xué)力學(xué)與土木工程學(xué)院,江蘇 徐州 221116;2中節(jié)能城市節(jié)能研究院有限公司,江蘇 常州 213001)

    超聲波液體輔助凍結(jié)在食品、海水淡化及冰蓄冷等方面的應(yīng)用得到了廣泛關(guān)注。結(jié)合聲場理論,在對液滴界面處熱質(zhì)傳遞分析的基礎(chǔ)上,建立了超聲波作用下的液滴冷卻凍結(jié)數(shù)學(xué)模型,討論了不同超聲波頻率、強(qiáng)度及作用時間對液滴界面處熱質(zhì)傳遞的影響,得到了超聲波作用下液滴的冷卻凍結(jié)規(guī)律。結(jié)果表明:超聲波作用下界面處的傳質(zhì)系數(shù)隨著超聲波強(qiáng)度的增加而增大,但隨超聲波頻率的增加而減小;在液滴冷卻凍結(jié)過程中,質(zhì)量傳遞與熱量傳遞的方向相同,液滴的冷卻凍結(jié)在超聲波作用下得到了強(qiáng)化。在超聲波頻率為20000 Hz(強(qiáng)度為400 W?m-2),經(jīng)過相同的時間(60 s),超聲波作用下的液滴溫度比無超聲作用的溫度低2.0~2.5℃。研究將有助于深入理解超聲波輔助冷卻凍結(jié)機(jī)理,并為其工程應(yīng)用提供了指導(dǎo)和參考。

    超聲波;凍結(jié);傳質(zhì);相變;溶液

    引 言

    超聲波輔助液體凍結(jié)具有促進(jìn)液相凍結(jié)、控制液相晶體生長速度的作用[1-2],并逐漸在食品工業(yè)、制冰及醫(yī)學(xué)領(lǐng)域中得到了廣泛應(yīng)用[3-5]。超聲波分為兩種,一種為高頻低能的診斷超聲波(100×104~3000×104Hz);另一種為低頻高能的功率超聲波(2×104~10×104Hz),用于強(qiáng)化溶液結(jié)晶。超聲波強(qiáng)化液相結(jié)晶主要是通過超聲波在液相中引起的機(jī)械空化效應(yīng)產(chǎn)生作用,并影響冰晶的成核和發(fā)展。

    近年來對超聲波作用下液體凍結(jié)的研究受到了廣泛的關(guān)注。余德洋等[6-8]通過超聲結(jié)晶實驗分析了超聲波作用下的液體結(jié)晶成核,明確超聲空化效應(yīng)是影響溶液結(jié)晶的主要因素。文獻(xiàn)[9-12]中利用低溫電鏡等手段分析了超聲波作用參數(shù)對冰晶尺寸、凍干速率以及溶液強(qiáng)化凍結(jié)的影響,分析表明超聲波能夠有效強(qiáng)化晶核的生成,對溶液的凍結(jié)過程起到促進(jìn)作用。王全海[13]、Saclier等[14]和Jordens等[15]從水的結(jié)晶理論出發(fā)研究了超聲波對結(jié)晶成核及過程的影響,結(jié)果表明超聲波具有細(xì)化晶粒的作用和促進(jìn)晶核生成的作用。Inada等[1]、Gondrexon等[16]和Wohlgemuth等[17]進(jìn)一步通過實驗觀察研究了超聲波作用下水凍結(jié)的過程,表明超聲振動是影響凍結(jié)相變的主要因素,但沒有進(jìn)行相關(guān)的理論分析。Fen等[18]對除氣水和未除氣水在超聲波作用下的冷凍過程進(jìn)行了研究,結(jié)果表明未除氣水在超聲波作用下冰晶更易生成。同時,文獻(xiàn)[19-22]中分析研究了超聲波作用下冰晶成核溫度對冰晶的一次成核、二次成核以及凍結(jié)的影響,以實現(xiàn)對凍結(jié)狀態(tài)的控制。對于凍結(jié)過程,如何通過加載超聲波實現(xiàn)對冰體生成速度和結(jié)構(gòu)的控制,成為制冰、冷凍法海水淡化及食品工業(yè)中所關(guān)注的問題之一。上述對超聲波作用下溶液凍結(jié)的研究,明確了空化效應(yīng)對溶液凍結(jié)的效果和作用,但未揭示超聲波對溶液凍結(jié)的內(nèi)部作用機(jī)制,未從凍結(jié)自身所必須滿足的能量角度對凍結(jié)過程進(jìn)行分析討論。

    本研究將在聲場理論的基礎(chǔ)上,結(jié)合超聲波在溶液內(nèi)的空化效應(yīng),通過對液滴凍結(jié)過程中熱質(zhì)傳遞的分析,在能量平衡的基礎(chǔ)上建立超聲波作用下液滴冷卻凍結(jié)的數(shù)學(xué)模型,并由此討論不同超聲波頻率、強(qiáng)度及作用時間對液滴冷卻凍結(jié)的影響,為深入理解超聲波輔助冷卻凍結(jié)機(jī)理及工程應(yīng)用提供理論基礎(chǔ)和指導(dǎo)。

    1 超聲波作用下的液滴凍結(jié)過程及數(shù)學(xué)模型

    1.1 聲場作用下的傳質(zhì)過程

    在超聲波作用下,液滴內(nèi)部由于空化效應(yīng)將產(chǎn)生大量的氣泡,而氣泡在液滴表面處溢出將有助于界面處的熱質(zhì)傳遞過程,從而影響液滴的凍結(jié)過程。在液滴凍結(jié)過程中,熱量傳遞的方向為從液滴內(nèi)部向液滴外部;傳質(zhì)的方向為由內(nèi)到外,兩者方向一致,傳質(zhì)強(qiáng)化了傳熱,將有助于液滴的凍結(jié)。超聲波作用下液滴內(nèi)氣泡的溢出見圖1。

    圖1 液滴內(nèi)氣泡的溢出Fig.1 Bubbles escaping from liquid surface

    空化效應(yīng)是由于聲波在液相中傳播時產(chǎn)生正負(fù)壓交替變化,而介質(zhì)中的裂隙由此產(chǎn)生,當(dāng)壓力為正時,迅速合并湮沒;壓力為負(fù)時,迅速膨脹[23]。在分析液滴凍結(jié)過程中,認(rèn)為氣泡為球形且分布均勻。

    超聲波作用在液相上的壓力變化為

    超聲波角頻率 ω=2πf=2π/T。

    超聲波強(qiáng)度為

    超聲波作用下液體內(nèi)在空化效應(yīng)作用下產(chǎn)生氣泡的最大和最小半徑分別為Rmax與Rmin[14],Rmax可由式(3)獲得

    Rmin可由式(4)獲得

    超聲波空化效應(yīng)產(chǎn)生的氣泡數(shù)為[24]

    式中,A是常數(shù)(10-10)。

    考慮到所有氣泡溢出是不可能的,令φ為界面處氣泡的溢出率,因此在液體界面處溢出的氣泡數(shù)為

    當(dāng)一個氣泡離開液體表面時,所引起的表面更新面積為

    由式(6)、式(7),可得總的表面更新面積為

    根據(jù)Fick定律,液滴界面處的質(zhì)量傳遞滿足

    式中,Cw為界面處的水蒸氣濃度;CAf為計算位置處的水蒸氣濃度。

    根據(jù)表面更新理論及文獻(xiàn)[25-26]中對傳質(zhì)系數(shù)的計算,可得液滴表面對流傳質(zhì)系數(shù)為

    在式(9)、式(10)中,超聲波作用下可以認(rèn)為τc=1/f。

    球狀液滴表面的對流傳質(zhì)系數(shù)為[27]

    因此,液滴表面的對流傳質(zhì)系數(shù)可以認(rèn)為由兩部分組成:一部分為氣泡溢出引起的部分,可以采用式(10)進(jìn)行計算;其他部分采用式(11)計算??偟囊旱伪砻娴膶α鱾髻|(zhì)系數(shù)為

    式中,φ為液滴表面更新比例,φ=SN/SL,SL為液滴外表面,

    1.2 聲場作用下的液滴凍結(jié)模型

    液滴與周圍環(huán)境的對流傳熱系數(shù)為[27-28]

    超聲波在介質(zhì)中傳播時存在熱效應(yīng),可表示為[29]

    式中,α為吸聲系數(shù),α=0.26f1.1。

    圖2 超聲波作用下液滴凍結(jié)物理模型Fig.2 Physical model of droplet freezing

    2 超聲波對液滴凍結(jié)模擬結(jié)果及分析

    為了揭示超聲波對液滴凍結(jié)過程的影響,先后分別對不同頻率、強(qiáng)度的超聲波作用下,液滴內(nèi)氣泡數(shù)、氣泡半徑、液滴表面更新比、液滴凍結(jié)半徑等的變化規(guī)律進(jìn)行了分析研究。在分析計算過程中主要計算參數(shù)見表1。

    表1 計算中的主要參數(shù)Table 1 Main parameters in calculation

    超聲波強(qiáng)度和頻率對液滴內(nèi)氣泡最大半徑的影響見圖3。圖3(a)、(b)表明在一定的超聲波頻率下,液滴內(nèi)在超聲波作用下產(chǎn)生的氣泡最大半徑隨著超聲波強(qiáng)度的增大而增大;氣泡最大半徑在相同的超聲波強(qiáng)度下,隨著超聲波頻率的增強(qiáng)而減小。原因為超聲波強(qiáng)度有利于氣泡的生成發(fā)展,而超聲波頻率的增大將導(dǎo)致作用在液滴內(nèi)的壓力發(fā)生較快的變化,導(dǎo)致氣泡更加劇烈地相互碰撞湮沒,從而使氣泡半徑減小。因此在超聲波輔助凍結(jié)中,為了強(qiáng)化界面處的傳質(zhì)進(jìn)而促進(jìn)凍結(jié),應(yīng)選擇低頻高強(qiáng)度的超聲波。

    圖3 不同超聲波強(qiáng)度和頻率下液滴內(nèi)氣泡的最大半徑Fig.3 Maximal bubbles radius in different ultrasound intensity and ultrasound frequency

    不同超聲波強(qiáng)度和頻率對液滴內(nèi)氣泡數(shù)的影響見圖4。從圖4(a)可以看出在一定的超聲波頻率下,隨著超聲波強(qiáng)度的增加,液滴內(nèi)產(chǎn)生的氣泡數(shù)在下降;而在一定的超聲波強(qiáng)度下,隨著超聲波頻率的增加,氣泡數(shù)在增加。原因為在一定的超聲波頻率下,高強(qiáng)度超聲波強(qiáng)化了氣泡的碰撞湮沒,從而使氣泡數(shù)下降;而超聲波強(qiáng)度一定時,高頻超聲波導(dǎo)致壓力發(fā)生較快的變化,從而使液滴內(nèi)氣泡數(shù)增加。

    圖5給出了不同超聲波頻率和強(qiáng)度下,液滴表面更新率的變化規(guī)律。從圖5可以看出,在一定的超聲波頻率下,隨著超聲波強(qiáng)度的增大,液滴表面更新率逐漸增加;當(dāng)超聲波頻率增大時,表面更新率下降??傮w來說,液滴表面更新率與液滴內(nèi)氣泡數(shù)的多少相關(guān),氣泡數(shù)越多,液滴表面更新率越大;氣泡數(shù)越少,液滴表面更新率越小。

    圖4 不同超聲波強(qiáng)度和頻率對液滴內(nèi)氣泡數(shù)的影響Fig.4 Bubbles numbers in different ultrasound intensity and ultrasound frequency

    圖5 不同超聲波強(qiáng)度對液滴表面更新率的影響Fig.5 Variation of surface renewal ratio in different ultrasonic intensity

    圖6給出了超聲波對液滴表面?zhèn)髻|(zhì)系數(shù)的影響。從圖中可以發(fā)現(xiàn),超聲波作用下液滴表面的傳質(zhì)得到了強(qiáng)化,當(dāng)超聲波頻率為20000 Hz、強(qiáng)度為200 W?m-2時,液滴表面對流傳質(zhì)系數(shù)為 0.0122 m?s-1,比無超聲波作用時液滴表面的對流傳質(zhì)系數(shù)要大。原因在于超聲波能夠強(qiáng)化液滴表面氣泡的脫離,從而促進(jìn)液滴表面的傳質(zhì),使得液滴表面的傳質(zhì)系數(shù)增加。

    圖6 超聲波對液滴表面?zhèn)髻|(zhì)系數(shù)的影響Fig.6 Variation of convective mass transfer coefficient in different ultrasonic intensity

    為了反映出超聲波導(dǎo)致的傳質(zhì)引起的傳熱占全部傳熱量的份額,引入ψ=qm/(qh+qm)表示。

    圖7給出了由超聲波引起的傳質(zhì)而導(dǎo)致的傳熱量占總傳熱量的比例,從圖中可以看出,當(dāng)超聲波頻率為 20000 Hz、超聲波強(qiáng)度為 200 W?m-2時,ψ=0.1695;當(dāng)超聲波頻率為80000 Hz、超聲波強(qiáng)度為200 W?m-2時,ψ=0.1645。表明超聲波引起的質(zhì)擴(kuò)散導(dǎo)致的傳熱量占總傳熱量的16%左右,說明在液滴凍結(jié)過程中加載超聲波將有助于液滴的凍結(jié)。

    圖7 超聲波引起的傳質(zhì)導(dǎo)致的傳熱量占全部傳熱量的份額Fig.7 Proportion between heat quantity by mass transfer and total heat

    從圖8可以看出超聲波作用下液滴的溫度比無超聲波作用下的液滴溫度下降得快,說明超聲波作用強(qiáng)化了液滴表面的傳質(zhì),由于傳質(zhì)與傳熱同向,傳質(zhì)強(qiáng)化了傳熱,因此超聲波作用下液滴的溫度下降速度要比無超聲波時快。

    圖8 超聲波與無超聲波作用下液滴的溫度變化規(guī)律Fig.8 Variation of droplet temperature in effect of ultrasonic and un-ultrasonic

    當(dāng)液滴的半徑為0.004 m、液滴初始溫度為10℃、超聲波頻率為20000 Hz時,圖8給出超聲波和無超聲波作用下液滴的溫度變化規(guī)律;當(dāng)超聲波(超聲波強(qiáng)度為400 W?m-2)作用時間為30 s時液滴的溫度比無超聲波作用的液滴溫度低約1.0℃,當(dāng)超聲波(超聲波強(qiáng)度為 400 W?m-2)作用時間為60 s時液滴的溫度比無超聲波作用低1.9℃,當(dāng)超聲波作用時間為90 s時液滴的溫度比無超聲波作用低2.5℃。

    圖9給出了超聲波頻率為20000 Hz、超聲波強(qiáng)度為200 W?m-2時,不同大小液滴在超聲波作用下的溫度變化規(guī)律。從圖9可以看出,液滴半徑為0.01 m經(jīng)過200 s后,溫度從10℃下降至0.4℃;液滴半徑為0.015 m時,溫度從10℃下降至5.1℃;液滴半徑為0.02 m時,溫度從10℃下降至7.2℃。說明在相同的超聲波作用下,液滴半徑越小,越有利于液滴的快速降溫凍結(jié)。

    圖9 超聲波作用下液滴溫度的變化規(guī)律Fig.9 Variation of droplet temperature with time

    超聲波作用下液滴的凍結(jié)過程是從液滴外表逐漸向里發(fā)展,直至完全凍結(jié)。圖10給出了當(dāng)超聲波頻率為20000 Hz、超聲波強(qiáng)度為100 W?m-2、液滴外氣流速度為 0.08 m?s-1時,在液滴凍結(jié)過程中液滴內(nèi)部液固界面的漸變規(guī)律。從圖10可以看出,隨著液滴半徑的減小,凍結(jié)時間也減小。液滴半徑為0.01 m,液固界面經(jīng)過50 s可以發(fā)展至液滴中心位置;液滴半徑為0.015 m,液固界面經(jīng)過 220 s可以發(fā)展至液滴中心位置;液滴半徑為0.02 m,需要經(jīng)過630 s液固界面可以發(fā)展至液滴中心位置。

    圖11給出了當(dāng)超聲波頻率為20000 Hz、超聲波強(qiáng)度100 W?m-2、液滴半徑為0.02 m時,液滴凍結(jié)過程中液滴內(nèi)液固界面的漸變規(guī)律。從圖11中可以發(fā)現(xiàn),外部氣流速度為0.8 m?s-1時,液固界面經(jīng)過200 s可以發(fā)展至液滴中心位置;外部氣流速度為0.08 m?s-1時,液固界面需經(jīng)過600 s才可以發(fā)展至液滴中心位置。說明增加外部氣流速度可以強(qiáng)化液滴的凍結(jié)過程。

    在液滴半徑為0.02 m,當(dāng)超聲波頻率為20000 Hz、液滴外部氣流速度為0.08 m?s-1時,不同超聲波強(qiáng)度下,液滴內(nèi)部液固界面的漸變規(guī)律見圖12。從圖可以看出,不同超聲波強(qiáng)度下,液滴內(nèi)液固界面的漸變曲線非常接近,表明在超聲波輔助凍結(jié)工程應(yīng)用中,應(yīng)采取低強(qiáng)度超聲波。

    圖11 不同外部氣流速度對液滴凍結(jié)過程中液固界面的影響Fig.11 Solid-liquid interface position in droplet freezing process for different gas velocity

    圖12 不同超聲波強(qiáng)度對液滴凍結(jié)過程中液固界面的影響Fig.12 Solid-liquid interface position in droplet freezing process for different ultrasonic intensity

    3 超聲波作用下的液滴凍結(jié)實驗

    超聲波用下液滴的凍結(jié)實驗系統(tǒng),見圖13。實驗系統(tǒng)中:恒溫恒濕箱(1)尺寸為 300 mm×300 mm×400 mm,液滴被置于熱電偶(4)上,通過熱電偶(4)可以獲得液滴的溫度變化。液滴的狀態(tài)可以通過高度攝像機(jī)(7)和電子顯微鏡(8)進(jìn)行觀察。溫度傳感器(2)用來測量恒溫恒濕箱(1)中的溫度,不同頻率和強(qiáng)度的超聲波由超聲波發(fā)生器(11)來提供。

    圖13 超聲波作用下的液滴凍結(jié)實驗系統(tǒng)Fig.13 Droplet freezing system by ultrasound

    表2 主要測試儀器參數(shù)Table 2 Technical specification of instrumentations &devices used in experiment system

    實驗過程中,液滴半徑大小為0.01~0.02 m,恒溫恒濕機(jī)中環(huán)境溫度為-10.0℃,考慮到液滴溫度的測量是在恒溫恒濕箱中進(jìn)行,實驗過程中由于液滴初始溫度為10.0℃,液滴與恒溫恒濕箱體環(huán)境之間存在溫差傳熱,進(jìn)而引起周圍空氣的微運(yùn)動,因此針對實驗工況的理論分析中將液滴周圍的氣流速度選取為0.005 m?s-1。實驗過程中,液滴周圍超聲波強(qiáng)度通過TD-YP0511C進(jìn)行測量,實驗過程中的主要測試儀器參數(shù)見表2。

    對超聲波和無超聲波作用下液體中的氣泡進(jìn)行了觀測,見圖14。從圖14可以看出,在超聲波和無超聲波作用下液體的狀態(tài)是不同的,超聲波作用下液體中產(chǎn)生的氣泡數(shù)明顯多于無超聲波作用的。原因在于超聲波空化效應(yīng)有助于液體中氣泡的生成。

    Fig.14 超聲波和無超聲波作用下液體中的氣泡Fig.14 Bubbles in liquid without ultrasound and with ultrasound

    為了理解超聲波作用下的凍結(jié)機(jī)理,利用放大倍數(shù)300倍的顯微鏡對超聲波頻率為20000 Hz和無超聲波作用時液滴的凍結(jié)進(jìn)行了對比觀測,見圖15、圖16。從圖15、圖16的對比來看,超聲波作用下能夠明顯觀測到液體中的氣泡,而無超聲波作用的液滴內(nèi)氣泡較少(較為平靜)。在此過程中,液滴外表面處的傳質(zhì)與液滴的傳熱方向一致(均為由內(nèi)到外),因此傳質(zhì)強(qiáng)化了傳熱,有利于液滴的冷卻凍結(jié)。

    為了驗證超聲波凍結(jié)理論模型的正確性,選取液滴溫度的變化來反映超聲波作用下液滴的凍結(jié)過程,對超聲波作用下液滴(液滴半徑分別為0.01和0.015 m)凍結(jié)過程中的溫度變化進(jìn)行了理論和實驗對比,并給出了無超聲波作用時液滴的溫度變化見圖17。從圖17可以看出,理論曲線與實驗測量值吻合,當(dāng)液滴半徑為0.01 m時,理論值和實驗值的平均誤差為 18.3%;與無超聲波作用的液滴溫度變化相比,超聲波能夠明顯加快液滴的冷卻冷凍凍結(jié),在此過程中當(dāng)時間為160 s時,超聲波作用下的液滴溫度(實驗值)比無超聲波作用下的液滴溫度(實驗值)低3.13℃;當(dāng)液滴半徑為0.015 m時,理論值和實驗值的平均誤差為 10.2%;在此過程中當(dāng)時間為160 s時,超聲波作用下的液滴溫度(實驗值)比無超聲波作用下的液滴溫度(實驗值)低1.65℃。這一對比表明了超聲波作用下的熱質(zhì)傳遞及凍結(jié)模型的正確性,可以利用該模型來對超聲波作用下的液體凍結(jié)相關(guān)工程應(yīng)用進(jìn)行分析與指導(dǎo)。

    圖15 超聲波作用下液滴的凍結(jié)狀態(tài)Fig.15 State variation in freezing process with ultrasound

    圖16 無超聲波作用下液滴的凍結(jié)狀態(tài)Fig.16 State variation in freezing process without ultrasound

    圖17 超聲波作用下液滴溫度的變化Fig.17 Variation of droplet temperature with time in ultrasound intensity 100 W?m-2

    4 結(jié) 論

    結(jié)合聲場理論,在對超聲波作用下液滴冷卻凍結(jié)過程中界面處熱質(zhì)傳遞分析的基礎(chǔ)上,建立了超聲波作用液滴的冷卻凍結(jié)模型,通過分析得到了超聲波作用下液滴的冷卻凍結(jié)規(guī)律,結(jié)論如下。

    (1)超聲波空化效應(yīng)有助于溶液界面處的傳質(zhì),由于在冷卻凍結(jié)過程中與傳熱方向一致,強(qiáng)化了傳熱從而有利于液滴的冷卻凍結(jié);

    (2)超聲波作用下液滴界面處的傳質(zhì)系數(shù)隨著超聲波強(qiáng)度的增加而增大,但隨超聲波頻率的增加而減??;

    (3)由超聲波導(dǎo)致的傳質(zhì)引起的傳熱占總傳熱量份額為15%~18%。當(dāng)超聲波頻率為20000 Hz、超聲波強(qiáng)度為200 W?m-2時,超聲波引起的強(qiáng)化傳熱效果為ψ=0.1695;

    (4)當(dāng)超聲波頻率為20000 Hz、超聲波強(qiáng)度為400 W?m-2時,經(jīng)過相同的時間(60 s后),超聲波作用下液滴的溫度比無超聲作用液滴的溫度低約2.0~2.5℃;

    (5)相同頻率超聲波下,超聲波強(qiáng)度對液滴冷卻凍結(jié)的影響不明顯;與高頻超聲波相比,低頻超聲波更有利于溶液液滴的冷卻凍結(jié)。

    符 號 說 明

    CAf——計算位置處的水蒸氣濃度,kg?m-3

    Cw——界面處的水蒸氣濃度,kg?m-3

    c ——聲波速率,m?s-1

    D ——分子擴(kuò)散系數(shù),m2?s

    f ——超聲波頻率,Hz

    I ——超聲波強(qiáng)度,W?m-2

    Le——蒸發(fā)潛熱,J?kg-1

    Lf——凝固潛熱,J?kg-1

    n ——氣泡數(shù)

    P ——超聲波功率,W

    pA——超聲波壓力(振幅),Pa

    pa——超聲波作用在液相上的壓力,Pa

    ph——流體靜壓,Pa

    R ——半徑,m

    S' ——超聲波作用面積,m2

    SL——液滴外表面,m2

    T ——溫度,K

    Ta——環(huán)境溫度,K

    Td——液滴的溫度,K

    t ——超聲波作用時間,s

    λ ——熱導(dǎo)率,W?m-1·K-1

    ν ——運(yùn)動黏度,N?s?m-2

    ρ ——密度,kg?m-3

    ρva——周圍環(huán)境水蒸氣濃度,kg?m-3

    ρvas——液滴表面水蒸氣濃度,kg?m-3

    σ ——表面張力,N?m-1

    τc——暴露時間,s

    φ ——界面處氣泡的溢出率

    φ ——液滴表面更新比例

    ω ——超聲波角頻率,rad?s-1

    下角標(biāo)

    L ——液滴

    max ——最大

    min ——最小

    [1] INADA T, ZHANG X, YABE A, et al. Active control of phase change from supercooled water to ice by ultrasonic vibration (Ⅰ):Control of freezing temperature[J]. International Journal of Heat and Mass Transfer, 2001, 44(23): 4523-4531.

    [2] KIANI H, ZHANG Z H, DELGADO A, et al. Ultrasound assisted nucleation of some liquid and solid model foods during freezing[J].Food Research International, 2011, 44(9): 2915-2921.

    [3] OLMO A, BAENA R, RISCO R. Use of a droplet nucleation analyzer in the study of water freezing kinetics under the influence of ultrasound waves[J]. International Journal of Refrigeration, 2008,31(2): 262-269.

    [4] KIANI H, ZHANG Z H, SUN D W. Effect of ultrasound irradiation on ice crystal size distribution in frozen agargel samples[J].Innovative Food Science & Emerging Technologies, 2013, 18(2):126-131.

    [5] AMIRA J H, ROMAN P, PIERRE L. Ultrasonically triggered freezing of aqueous solutions: influence of initial oxygen content on ice crystals’ size distribution[J]. Journal of Crystal Growth, 2014,402(9): 78-82.

    [6] YU D Y, LIU B L, WANG B C. The effect of ultrasonic waves on the nucleation of pure water and degassed water[J]. Ultrason Sonochem,2012, 19(3): 459-463.

    [7] 余德洋, 劉寶林, 王伯春. 超聲場中聲壓與空化對冰晶分裂的影響[J]. 制冷學(xué)報, 2011, 32(6): 30-34.YU D Y, LIU B L, WANG B C. The effect of acoustic pressure and cavitation on the secondary nucleation of ice[J]. Journal of Refrigeration, 2011, 32(6): 30-34.

    [8] 余德洋. 超聲波強(qiáng)化溶液凍結(jié)的機(jī)理研究[D]. 上海: 上海理工大學(xué), 2012.YU D Y. The study on the mechanism of enhancement of solution crystallization by power ultrasound[D]. Shanghai: University of Shanghai for Science and Technology, 2012.

    [9] 周新麗, 騰蕓, 戴澄. 超聲波平板冷凍提高胡蘿卜凍干速率[J]. 農(nóng)業(yè)工程學(xué)報, 2017, 33(1): 256-261.ZHOU X L, TENG Y, DAI C. Contact ultrasound freezing improving freeze drying rate of carrot[J]. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(1): 256-261.

    [10] XU B G, ZHANG M, BHANDARI B, et al. Influence of power ultrasound on ice nucleation of radish cylinders during ultrasound-assisted immersion freezing[J]. International Journal of Refrigeration, 2014, 46(1): 1-8.

    [11] GIELEN B, KUSTERS P, JORDENS J, et al. Energy efficient crystallization of paracetamol using pulsed ultrasound[J]. Chemical Engineering and Processing: Process Intensification, 2017, 114(1):55-66.

    [12] LOUHI-KULTANEN M, KARJALAINEN M, RANTANEN J, et al.Crystallization of glycine with ultrasound[J]. International Journal of Pharmaceutics, 2006, 320(1): 23-29.

    [13] 王全海. 過冷水動態(tài)結(jié)晶的超聲機(jī)理研究[D]. 洛陽: 河南科技大學(xué), 2014.WANG Q H. The study on the mechanism of enhancement of supercooled water crystallization by power ultrasound[D]. Luoyang:Henan University of Science and Technology, 2014.

    [14] SACLIER M, PECZALSKI R, ANDRIEU J. A theoretical model for ice primary nucleation induced by acoustic cavitation[J]. Ultrasonics Sonochemistry, 2010, 17(17): 98-105.

    [15] JORDENS J, GIELEN B, BRAEKEN L, et al. Determination of the effect of the ultrasonic frequency on the cooling crystallization of paracetamol[J]. Chemical Engineering and Processing: Process Intensification, 2014, 84(84): 38-44.

    [16] GONDREXON N, CHEZE L, JIN Y, et al. Intensification of heat and mass transfer by ultrasound: application to heat exchangers and membarane separation processes[J]. Ultrasonic Sonochemistry, 2015,25(1): 40-50.

    [17] WOHLGEMUTH K, KORDYLLA A, RUETHER F, et al.Experimental study of the effect of bubbles on nucleation during batch cooling crystallization[J]. Chemical Engineering Science, 2009,64(19): 4155-4163.

    [18] FEN H, SUN D W, GAO W H, et al. Effects of pre-existing bubbles on ice nucleation and crystallization during ultrasound-assisted freezing of water and sucrose solution[J]. Innovative Food Science and Emerging Technologies, 2013, 20(4): 161-166.

    [19] MORRIS G J, ACTON E. Controlled ice nucleation in cryopreservation - a review[J]. Cryobiology, 2013, 66(2): 85-92.

    [20] GEIDOBLER R, WINTER G. Controlled ice nucleation in the field of freeze drying: fundamentals and technology review[J]. Eur. J. Pharm.Biopharm., 2013, 85(2): 214-222.

    [21] PETZOLD G, AGUILERA J M. Ice morphology: fundamentals and technological applications in foods[J]. Food Biophysics, 2009, 4(4):378-396.

    [22] SEARLES J A, CARPENTER J F, RANDOLPH T W. The ice nucleation temperature determines the primary drying rate of lyophilization for samples frozen on a temperature-controlled shelf[J].Pharmaceutical Sciences, 2001, 90(7): 860-871.

    [23] ZHENG L Y, SUN D W. Innovative applications of power ultrasound during food freezing process-a review[J]. Trends in Food Science and Technology, 2006, 17(1): 16-23.

    [24] WALTON A J, REYNOLDS G T. Sonoluminescence[J]. Advance in Phycics, 1984, 33(6): 595-660.

    [25] 高蓬輝, 紀(jì)邵斌, 衡文佳, 等. 濕度差驅(qū)動下溶液蒸發(fā)冷凍過程中冰體的發(fā)展規(guī)律[J]. 化工學(xué)報, 2013, 64(8): 2820-2826.GAO P H, JI S B, HENG W J, et al. Icing evolution in process of liquor evaporation-refrigeration driven by humidity difference[J].CIESC Journal, 2013, 64(8): 2820-2826.

    [26] 高蓬輝, 衡文佳, 周興業(yè), 等. 臨界條件(0℃)下溶液蒸發(fā)冷凍過程中的傳質(zhì)規(guī)律[J]. 化工學(xué)報, 2013, 64(9): 3206-3212.GAO P H, HENG W J, ZHOU X Y, et al. Mass transfer of liquor in evaporation-refrigeration process under critical condition[J]. CIESC Journal, 2013, 64(9): 3206-3212.

    [27] SPENGLER J D, GOKHALE N R. Freezing of freely suspended,supercooled water drops in a large vertical wind tunnel[J]. Journal of Applied Meteorology, 1972, 11(11): 1101-1107.

    [28] HINDMARSH J P, RUSSELL A B, CHEN X D. Experimental and numerical analysis of the temperature transition of a suspended freezing water droplet[J]. International Journal of Heat and Mass Transfer, 2003, 46(7): 1199-1213.

    [29] MASON T J, LORIMER P J. Sonochemistry: Theory, Application and Uses of Ultrasonic in Chemistry[M]. New York: Ellis Horwood Press, 1988: 96-132.

    date:2017-04-10.

    GAO Penghui, gaopenghui2004@126.com

    supported by the National Natural Science Foundation of China (51106176) and the Subject Front Research of China University of Mining and Technology (2015XKQY16).

    Cooling and freezing law for liquid drop in ultrasound wave

    GAO Penghui1, ZHANG Meng1, DU Yuji1,2, CHENG Bo1, ZHANG Donghai1, ZHOU Guoqing1
    (1School of Architecture and Civil Engineering,China University of Mining and Technology,Xuzhou221116,Jiangsu,China;2China Energy Conservation and Environmental Protection City Energy Conservation Company Limited,Changzhou213001,Jiangsu,China)

    The application of ultrasound to liquid freezing got much attention over the last few years and its potential seems very promising. In order to make clear droplet freezing assisted by ultrasound, the heat and mass transfer characteristic was studied based on ultrasound theory, penetration theory of mass transfer and energy conservation. The ultrasound frequency, ultrasound intensity and operation time were studied in the process of the liquid drop freezing. The results showed that ultrasound could accelerate mass transfer and make droplet rapidly cooling. In the effect of ultrasound, the bubble size in the droplet was decreased with ultrasound frequency, and the bubble number in the droplet was increased with ultrasound frequency. Mass transfer coefficient of droplet was increased with ultrasound intensity and reduced with ultrasound frequency. For the mass transfer and heat transfer,the direction were same in the droplet freezing process, the heat transfer was strengthened by mass transfer in the droplet freezing process. Comparing with no ultrasound, droplet temperature with ultrasound (ultrasound frequency 20000 Hz and ultrasound intensity 400 W?m-2) was lower 2.0—2.5℃ after the same time (60 s). Hence the ultrasound helps to cool and freeze droplet. This study is favor to understanding the freezing by ultrasound and its application.

    ultrasound; freezing; mass transfer; phase transformation; liquid

    TB 66

    A

    0438—1157(2017)11—4095—10

    10.11949/j.issn.0438-1157.20170381

    2017-04-10收到初稿,2017-07-06收到修改稿。

    聯(lián)系人及第一作者:高蓬輝(1979—),男,博士。

    國家自然科學(xué)基金項目(51106176);中國礦業(yè)大學(xué)學(xué)科前沿專項項目(2015XKQY16)。

    猜你喜歡
    傳質(zhì)液滴氣泡
    檸檬氣泡水
    欣漾(2024年2期)2024-04-27 15:19:49
    SIAU詩杭便攜式氣泡水杯
    新潮電子(2021年7期)2021-08-14 15:53:12
    浮法玻璃氣泡的預(yù)防和控制對策
    液滴間相互碰撞融合與破碎的實驗研究
    噴淋液滴在空氣環(huán)境下的運(yùn)動特性
    冰凍氣泡
    氨基酸鹽吸收二氧化碳過程的傳質(zhì)特性
    PTFE膜吸收CO2的工藝參數(shù)對傳質(zhì)性能的影響
    清潔轉(zhuǎn)向酸H+表面?zhèn)髻|(zhì)行為實驗研究
    煤顆粒熱解的傳熱傳質(zhì)分析
    計算物理(2014年1期)2014-03-11 17:00:36
    熟女av电影| 亚洲三区欧美一区| 国产 精品1| 自线自在国产av| 成人午夜精彩视频在线观看| 中文字幕人妻丝袜制服| xxx大片免费视频| 日本av手机在线免费观看| 久久久久国产精品人妻一区二区| 高清在线视频一区二区三区| 十八禁网站网址无遮挡| 91国产中文字幕| 欧美中文综合在线视频| 日韩免费高清中文字幕av| 国产黄色视频一区二区在线观看| 欧美av亚洲av综合av国产av | 一二三四在线观看免费中文在| 两个人免费观看高清视频| 亚洲三区欧美一区| 叶爱在线成人免费视频播放| 美女福利国产在线| 欧美+日韩+精品| 国产在视频线精品| 看非洲黑人一级黄片| 在线亚洲精品国产二区图片欧美| 久久精品国产亚洲av高清一级| 丰满迷人的少妇在线观看| 国产精品不卡视频一区二区| 亚洲三级黄色毛片| 99re6热这里在线精品视频| 久久久久精品人妻al黑| 女人久久www免费人成看片| 母亲3免费完整高清在线观看 | 亚洲情色 制服丝袜| 日本爱情动作片www.在线观看| 亚洲精品国产av成人精品| 最近2019中文字幕mv第一页| 国产精品香港三级国产av潘金莲 | 一本大道久久a久久精品| 日本欧美国产在线视频| 亚洲欧美中文字幕日韩二区| 青春草视频在线免费观看| 91精品三级在线观看| 波野结衣二区三区在线| 免费播放大片免费观看视频在线观看| 亚洲欧洲日产国产| 在线观看免费高清a一片| 天堂俺去俺来也www色官网| 国产乱人偷精品视频| 日韩欧美一区视频在线观看| 亚洲中文av在线| 亚洲第一区二区三区不卡| 汤姆久久久久久久影院中文字幕| 亚洲一级一片aⅴ在线观看| 91精品国产国语对白视频| 一本大道久久a久久精品| 国产乱来视频区| 亚洲欧美中文字幕日韩二区| 国产日韩欧美视频二区| 亚洲国产最新在线播放| 亚洲中文av在线| a 毛片基地| 国产97色在线日韩免费| 日本黄色日本黄色录像| 美女国产高潮福利片在线看| 亚洲国产精品一区二区三区在线| 极品少妇高潮喷水抽搐| 亚洲精品av麻豆狂野| 亚洲精品久久久久久婷婷小说| 久久久久视频综合| 三上悠亚av全集在线观看| 午夜福利视频在线观看免费| 五月开心婷婷网| 母亲3免费完整高清在线观看 | 伦理电影大哥的女人| 韩国高清视频一区二区三区| 看免费av毛片| 日本色播在线视频| 国产白丝娇喘喷水9色精品| 18禁动态无遮挡网站| 亚洲一码二码三码区别大吗| 日韩在线高清观看一区二区三区| 日韩免费高清中文字幕av| av免费观看日本| 免费看不卡的av| 一区二区三区乱码不卡18| 久久亚洲国产成人精品v| 黑人猛操日本美女一级片| 国产 精品1| 亚洲美女黄色视频免费看| 十八禁高潮呻吟视频| 9191精品国产免费久久| 国产福利在线免费观看视频| av在线播放精品| 亚洲av欧美aⅴ国产| 亚洲国产精品国产精品| √禁漫天堂资源中文www| 免费观看av网站的网址| av线在线观看网站| 亚洲综合色网址| 日本免费在线观看一区| 777久久人妻少妇嫩草av网站| 久久99热这里只频精品6学生| 精品午夜福利在线看| 少妇人妻 视频| 免费观看a级毛片全部| av片东京热男人的天堂| 黄色毛片三级朝国网站| 亚洲国产精品成人久久小说| 国产午夜精品一二区理论片| 纵有疾风起免费观看全集完整版| 亚洲美女视频黄频| 观看美女的网站| 午夜福利视频精品| 啦啦啦啦在线视频资源| 日韩在线高清观看一区二区三区| 电影成人av| 国产精品免费大片| 亚洲国产色片| 黑人巨大精品欧美一区二区蜜桃| 黄色一级大片看看| 成人黄色视频免费在线看| 欧美少妇被猛烈插入视频| 韩国av在线不卡| 久久婷婷青草| 精品福利永久在线观看| 久久青草综合色| 久久人人爽人人片av| 国产日韩欧美在线精品| 免费人妻精品一区二区三区视频| 黄片播放在线免费| 婷婷成人精品国产| 国产有黄有色有爽视频| 亚洲精品乱久久久久久| 亚洲av免费高清在线观看| 男女无遮挡免费网站观看| 色网站视频免费| 久久精品久久久久久噜噜老黄| videossex国产| 18禁观看日本| 日韩熟女老妇一区二区性免费视频| 七月丁香在线播放| 精品久久久精品久久久| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av中文av极速乱| 国产男人的电影天堂91| www日本在线高清视频| 午夜福利视频在线观看免费| 午夜福利乱码中文字幕| 亚洲少妇的诱惑av| 亚洲图色成人| 亚洲精品一区蜜桃| 亚洲图色成人| 久久亚洲国产成人精品v| 九草在线视频观看| 久久久久久久久久人人人人人人| 免费观看在线日韩| 男女下面插进去视频免费观看| 一级片免费观看大全| 在线天堂最新版资源| 久久99蜜桃精品久久| 免费在线观看黄色视频的| 18禁裸乳无遮挡动漫免费视频| 99久久综合免费| 热re99久久精品国产66热6| 高清视频免费观看一区二区| 成人亚洲精品一区在线观看| 99久久综合免费| 熟女少妇亚洲综合色aaa.| 免费观看无遮挡的男女| 欧美日韩精品网址| 免费在线观看视频国产中文字幕亚洲 | 欧美少妇被猛烈插入视频| 天美传媒精品一区二区| 大片免费播放器 马上看| 麻豆乱淫一区二区| videos熟女内射| 免费黄频网站在线观看国产| 巨乳人妻的诱惑在线观看| 日韩在线高清观看一区二区三区| 最近中文字幕2019免费版| 欧美bdsm另类| 777米奇影视久久| 最新中文字幕久久久久| 大码成人一级视频| 最近中文字幕高清免费大全6| 免费在线观看视频国产中文字幕亚洲 | 成年女人在线观看亚洲视频| 午夜福利,免费看| 在线天堂最新版资源| 成人黄色视频免费在线看| 精品午夜福利在线看| 亚洲情色 制服丝袜| 麻豆乱淫一区二区| 午夜影院在线不卡| 午夜福利乱码中文字幕| 精品一品国产午夜福利视频| 女人高潮潮喷娇喘18禁视频| 亚洲成人一二三区av| tube8黄色片| 一区二区日韩欧美中文字幕| 只有这里有精品99| 精品人妻偷拍中文字幕| 精品酒店卫生间| 另类精品久久| 黄色毛片三级朝国网站| 亚洲av.av天堂| 26uuu在线亚洲综合色| 午夜av观看不卡| 国产精品 欧美亚洲| 精品久久久精品久久久| 巨乳人妻的诱惑在线观看| 人人妻人人澡人人爽人人夜夜| 精品人妻偷拍中文字幕| 2018国产大陆天天弄谢| 久久久久精品性色| 乱人伦中国视频| 精品亚洲成a人片在线观看| 人人妻人人爽人人添夜夜欢视频| 大陆偷拍与自拍| 国产免费福利视频在线观看| 男女啪啪激烈高潮av片| 日产精品乱码卡一卡2卡三| 七月丁香在线播放| 国产免费一区二区三区四区乱码| 国产精品免费视频内射| 国产乱人偷精品视频| 十八禁网站网址无遮挡| 国产精品偷伦视频观看了| 亚洲国产欧美在线一区| 国产97色在线日韩免费| 色视频在线一区二区三区| 天堂8中文在线网| 最近手机中文字幕大全| 欧美 日韩 精品 国产| 精品人妻在线不人妻| 最近中文字幕2019免费版| 美女主播在线视频| 春色校园在线视频观看| 18+在线观看网站| 美女主播在线视频| 亚洲国产欧美网| 日日啪夜夜爽| 久久免费观看电影| 久久午夜综合久久蜜桃| 自线自在国产av| 人人妻人人澡人人爽人人夜夜| 可以免费在线观看a视频的电影网站 | 有码 亚洲区| 亚洲美女搞黄在线观看| 日本-黄色视频高清免费观看| 人人妻人人澡人人看| av国产久精品久网站免费入址| 国产黄色免费在线视频| 永久免费av网站大全| 巨乳人妻的诱惑在线观看| 亚洲av免费高清在线观看| 晚上一个人看的免费电影| 涩涩av久久男人的天堂| 亚洲精品美女久久久久99蜜臀 | 午夜福利网站1000一区二区三区| 日韩大片免费观看网站| 热re99久久国产66热| 国产成人aa在线观看| 性色av一级| 少妇被粗大的猛进出69影院| 国产熟女欧美一区二区| 777米奇影视久久| 天天躁日日躁夜夜躁夜夜| www日本在线高清视频| 亚洲精品一区蜜桃| 五月开心婷婷网| 少妇人妻 视频| 国产综合精华液| 国产精品一国产av| 大香蕉久久成人网| 极品人妻少妇av视频| 日韩大片免费观看网站| 国产精品免费大片| av国产久精品久网站免费入址| 国产老妇伦熟女老妇高清| 在线观看一区二区三区激情| 99国产综合亚洲精品| 2022亚洲国产成人精品| 日韩中字成人| 欧美国产精品一级二级三级| 国产在视频线精品| 成年av动漫网址| 少妇的逼水好多| 永久网站在线| 菩萨蛮人人尽说江南好唐韦庄| 久久人人爽av亚洲精品天堂| 麻豆乱淫一区二区| 日韩欧美一区视频在线观看| 街头女战士在线观看网站| 亚洲av福利一区| 咕卡用的链子| 亚洲国产欧美日韩在线播放| 午夜福利在线免费观看网站| 日韩伦理黄色片| 国产精品久久久久成人av| 久久精品国产亚洲av涩爱| 日日撸夜夜添| 国产精品一二三区在线看| 边亲边吃奶的免费视频| 丰满少妇做爰视频| 亚洲第一青青草原| 成人18禁高潮啪啪吃奶动态图| 精品亚洲乱码少妇综合久久| 亚洲精品视频女| 一本大道久久a久久精品| 国产精品国产av在线观看| 啦啦啦在线观看免费高清www| 建设人人有责人人尽责人人享有的| 亚洲av电影在线进入| 综合色丁香网| 男人添女人高潮全过程视频| 欧美bdsm另类| 看免费av毛片| 丝袜在线中文字幕| 亚洲精品视频女| 亚洲国产欧美在线一区| 妹子高潮喷水视频| 国产一区二区激情短视频 | 男女边摸边吃奶| 久久精品国产亚洲av高清一级| 亚洲欧洲日产国产| 性色av一级| 波野结衣二区三区在线| 亚洲精品久久久久久婷婷小说| xxx大片免费视频| 久久久久网色| kizo精华| 亚洲一区二区三区欧美精品| 国产精品.久久久| 久久久久久伊人网av| 国精品久久久久久国模美| 99热全是精品| 人体艺术视频欧美日本| 一区在线观看完整版| 午夜免费鲁丝| 亚洲国产欧美日韩在线播放| 久久精品人人爽人人爽视色| 成人影院久久| 亚洲精品中文字幕在线视频| 亚洲,欧美,日韩| 超色免费av| 亚洲欧美一区二区三区久久| 久久久久久久亚洲中文字幕| 制服丝袜香蕉在线| 男女国产视频网站| 建设人人有责人人尽责人人享有的| 国产欧美亚洲国产| 日本午夜av视频| 国产精品一区二区在线不卡| 99热全是精品| 91精品三级在线观看| 最新中文字幕久久久久| 一区二区三区四区激情视频| 欧美成人精品欧美一级黄| 尾随美女入室| 亚洲欧美一区二区三区黑人 | 中文乱码字字幕精品一区二区三区| 亚洲美女搞黄在线观看| 欧美日本中文国产一区发布| 亚洲人成电影观看| 丰满乱子伦码专区| 久久国内精品自在自线图片| 丝瓜视频免费看黄片| 欧美黄色片欧美黄色片| 街头女战士在线观看网站| 国产综合精华液| 婷婷色av中文字幕| 亚洲精品久久成人aⅴ小说| 日本-黄色视频高清免费观看| 欧美日韩精品成人综合77777| 免费在线观看视频国产中文字幕亚洲 | 色网站视频免费| 久久久久视频综合| 国产成人aa在线观看| 国产成人一区二区在线| 一级毛片 在线播放| 一区二区三区乱码不卡18| 欧美97在线视频| 在线看a的网站| 少妇人妻 视频| 另类精品久久| 亚洲人成电影观看| 蜜桃国产av成人99| 只有这里有精品99| 青春草国产在线视频| 久久人妻熟女aⅴ| www.熟女人妻精品国产| 欧美精品高潮呻吟av久久| 亚洲av欧美aⅴ国产| 免费观看无遮挡的男女| 天天影视国产精品| 亚洲人成电影观看| 精品国产超薄肉色丝袜足j| 国产欧美亚洲国产| 久久精品国产综合久久久| 亚洲,欧美,日韩| 少妇的丰满在线观看| 亚洲欧美色中文字幕在线| 一区在线观看完整版| 一区二区三区四区激情视频| 国产av一区二区精品久久| 18禁裸乳无遮挡动漫免费视频| 男人舔女人的私密视频| 亚洲成人av在线免费| 国产一区有黄有色的免费视频| 午夜福利在线免费观看网站| 成人毛片60女人毛片免费| 国产精品亚洲av一区麻豆 | 欧美日韩视频高清一区二区三区二| av网站在线播放免费| 啦啦啦在线观看免费高清www| 热re99久久国产66热| 午夜福利影视在线免费观看| 亚洲欧美清纯卡通| 一二三四在线观看免费中文在| 亚洲三区欧美一区| 美女视频免费永久观看网站| 亚洲欧美成人精品一区二区| 欧美激情高清一区二区三区 | 综合色丁香网| 精品国产国语对白av| 午夜免费鲁丝| 国产精品久久久久久久久免| 一级毛片黄色毛片免费观看视频| 少妇的丰满在线观看| 亚洲欧美中文字幕日韩二区| 久热久热在线精品观看| 日本黄色日本黄色录像| 伊人亚洲综合成人网| 精品久久蜜臀av无| 好男人视频免费观看在线| 丝袜喷水一区| 国产成人精品在线电影| 国产精品无大码| 久久人人爽av亚洲精品天堂| 亚洲欧美中文字幕日韩二区| 天天躁狠狠躁夜夜躁狠狠躁| 黄色一级大片看看| 国产黄色免费在线视频| 中文字幕色久视频| 一区二区三区四区激情视频| 少妇 在线观看| 欧美97在线视频| 午夜福利在线观看免费完整高清在| 国产精品无大码| 国产精品三级大全| www.精华液| 成人18禁高潮啪啪吃奶动态图| 久久精品国产综合久久久| 国产精品99久久99久久久不卡 | 人妻 亚洲 视频| 亚洲成色77777| 日韩精品有码人妻一区| 丰满迷人的少妇在线观看| www.自偷自拍.com| 大话2 男鬼变身卡| 国产男女内射视频| 激情视频va一区二区三区| 国产伦理片在线播放av一区| 国产av码专区亚洲av| 中文天堂在线官网| 精品午夜福利在线看| 啦啦啦在线观看免费高清www| 午夜激情久久久久久久| 午夜日韩欧美国产| 男女午夜视频在线观看| 亚洲欧美一区二区三区黑人 | 精品国产超薄肉色丝袜足j| 交换朋友夫妻互换小说| 777久久人妻少妇嫩草av网站| 日产精品乱码卡一卡2卡三| 亚洲欧美一区二区三区国产| 男女边摸边吃奶| 热re99久久精品国产66热6| 91国产中文字幕| 中文乱码字字幕精品一区二区三区| 亚洲一区二区三区欧美精品| 最近的中文字幕免费完整| 男女免费视频国产| 韩国av在线不卡| 日韩一本色道免费dvd| 国产片特级美女逼逼视频| 一区二区三区精品91| 99热全是精品| 视频在线观看一区二区三区| 中文字幕亚洲精品专区| 最黄视频免费看| 激情视频va一区二区三区| 精品国产一区二区久久| 老熟女久久久| 高清av免费在线| 亚洲三级黄色毛片| 熟女少妇亚洲综合色aaa.| 黄片小视频在线播放| 成人18禁高潮啪啪吃奶动态图| 午夜免费男女啪啪视频观看| 啦啦啦在线免费观看视频4| 日本黄色日本黄色录像| 日日摸夜夜添夜夜爱| 亚洲人成77777在线视频| 大片免费播放器 马上看| 捣出白浆h1v1| 日本爱情动作片www.在线观看| 国产精品二区激情视频| 成人18禁高潮啪啪吃奶动态图| 久久久国产精品麻豆| 国产精品久久久久久av不卡| 高清黄色对白视频在线免费看| 曰老女人黄片| 日韩中字成人| 午夜福利视频在线观看免费| 各种免费的搞黄视频| 久久久精品免费免费高清| www.熟女人妻精品国产| 黄色一级大片看看| 观看美女的网站| 午夜福利在线免费观看网站| 九色亚洲精品在线播放| 18在线观看网站| 男人舔女人的私密视频| www.精华液| 岛国毛片在线播放| 曰老女人黄片| 26uuu在线亚洲综合色| 国产成人精品一,二区| 亚洲国产精品一区三区| 亚洲经典国产精华液单| 水蜜桃什么品种好| 日韩精品有码人妻一区| 一区二区日韩欧美中文字幕| 一区二区三区四区激情视频| 永久免费av网站大全| 国产精品久久久久久精品古装| 成人国语在线视频| 亚洲国产精品999| 欧美日韩国产mv在线观看视频| 欧美xxⅹ黑人| 女性被躁到高潮视频| 街头女战士在线观看网站| 国产在线一区二区三区精| 男女国产视频网站| 国产精品一区二区在线观看99| 亚洲欧美中文字幕日韩二区| 日韩伦理黄色片| 精品一品国产午夜福利视频| 欧美精品av麻豆av| 欧美激情极品国产一区二区三区| 夜夜骑夜夜射夜夜干| 久久国产精品大桥未久av| 韩国av在线不卡| 美女主播在线视频| 新久久久久国产一级毛片| 久久精品久久久久久噜噜老黄| 久久人人爽人人片av| 色94色欧美一区二区| 咕卡用的链子| 欧美激情极品国产一区二区三区| 日本vs欧美在线观看视频| 亚洲精品中文字幕在线视频| 亚洲激情五月婷婷啪啪| 国产精品熟女久久久久浪| 人妻一区二区av| 久久精品久久精品一区二区三区| 人妻 亚洲 视频| 国产男人的电影天堂91| 韩国高清视频一区二区三区| 91aial.com中文字幕在线观看| 久久精品夜色国产| 成年美女黄网站色视频大全免费| 91精品国产国语对白视频| 久久精品人人爽人人爽视色| 成人影院久久| 欧美精品高潮呻吟av久久| 欧美亚洲 丝袜 人妻 在线| 日本欧美国产在线视频| 国产一区二区激情短视频 | 狂野欧美激情性bbbbbb| 九草在线视频观看| 两个人看的免费小视频| 伊人久久国产一区二区| 亚洲欧洲国产日韩| 久久精品久久久久久久性| 日本-黄色视频高清免费观看| 我的亚洲天堂| 亚洲,欧美精品.| 国产高清国产精品国产三级| 人妻一区二区av| 亚洲情色 制服丝袜| 国产成人午夜福利电影在线观看| 免费观看a级毛片全部| 999久久久国产精品视频| 欧美日韩成人在线一区二区| www日本在线高清视频| 久久久国产精品麻豆| 青春草亚洲视频在线观看| 亚洲久久久国产精品| 91久久精品国产一区二区三区| 国产乱人偷精品视频| 亚洲美女视频黄频| 亚洲三级黄色毛片| 日本猛色少妇xxxxx猛交久久| 亚洲精品久久成人aⅴ小说| 午夜激情久久久久久久| 精品亚洲成a人片在线观看| 欧美人与性动交α欧美精品济南到 | 亚洲美女搞黄在线观看| 久久精品久久久久久噜噜老黄| 人妻系列 视频| 69精品国产乱码久久久|