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

    激光驅(qū)動(dòng)液滴遷移的機(jī)理研究1)

    2022-03-20 15:52:26葉致君
    力學(xué)學(xué)報(bào) 2022年2期
    關(guān)鍵詞:溫度梯度母液毛細(xì)

    葉致君 段 俐 康 琦

    (中國科學(xué)院力學(xué)研究所微重力實(shí)驗(yàn)室,北京 100190)

    (中國科學(xué)院大學(xué)工程科學(xué)學(xué)院,北京 100049)

    引言

    液滴是自然界中的一類典型和普遍存在的物質(zhì)形態(tài).液滴的平衡、遷移和運(yùn)動(dòng)、以及液滴之間的相互作用等一直都是流體力學(xué)的研究熱點(diǎn).在微重力環(huán)境下,液滴的行為具有新的特性.由于有效重力作用極小,液滴受外部壓力場均勻,可以保持較好的球形,并且液滴可以在空中懸浮,為研究液滴的行為提供了極好的條件[1].在微重力環(huán)境中,重力場不再是驅(qū)動(dòng)液滴運(yùn)動(dòng)的主要作用場,溫度場引發(fā)的液滴界面張力不均勻,成為驅(qū)動(dòng)液滴運(yùn)動(dòng)的主要作用場,使液滴產(chǎn)生熱毛細(xì)遷移運(yùn)動(dòng).液滴的運(yùn)動(dòng)規(guī)律與眾多空間應(yīng)用技術(shù)密切相關(guān),在微重力條件下,對液滴運(yùn)動(dòng)規(guī)律與行為的研究在空間合金制備[2]、兩相分離[3]、晶體生長[4-5]、空間玻璃制備[6]、空間液滴燃燒與消防安全[7]、空間流體管理[1,8]等方面有著巨大的應(yīng)用價(jià)值,其在空間應(yīng)用中的重要性將得到進(jìn)一步的提升.對液滴運(yùn)動(dòng)規(guī)律的探索,實(shí)現(xiàn)液滴的主動(dòng)驅(qū)動(dòng)控制,是近幾年提出的新理論和新思想,不僅有重要的學(xué)術(shù)意義,同時(shí)具有巨大的工業(yè)價(jià)值.

    1959 年Young 等[9]針對氣泡或液滴熱毛細(xì)遷移問題提出YGB (Young-Goldstein-Block)理論,忽略了慣性項(xiàng)與對流項(xiàng)等非線性項(xiàng)的影響.介于此問題Subramanian[10]提出了準(zhǔn)穩(wěn)態(tài)假設(shè),并使用漸進(jìn)展開的方法得到液滴的遷移速度結(jié)果[11].Balasubramaniam 和Chai[12]考慮慣性項(xiàng)和對流項(xiàng)影響,得到液滴周圍的速度場,并得到液滴遷移速度以及液滴內(nèi)外的溫度場分布,發(fā)現(xiàn)在Re 和Ma 均較大時(shí),液滴周圍存在邊界層.1990 年Barton 和Subramanian[13]分析得到了液滴的熱毛細(xì)遷移速度、液滴周圍的速度場與溫度場.對于雙液滴熱毛細(xì)遷移研究,Keh 和Chen[14]分析了軸對稱雙液滴的熱毛細(xì)遷移過程,得到母液中的速度場和液滴遷移速度.

    對于單液滴穩(wěn)態(tài)遷移的定常數(shù)值模擬研究,隨著運(yùn)動(dòng)界面追蹤數(shù)值計(jì)算方法的發(fā)展而進(jìn)展迅速.Haj-Hariri 等[15]對三維可變形黏性液滴的熱毛細(xì)遷移進(jìn)行數(shù)值模擬,探討了慣性項(xiàng)、對流項(xiàng)、液滴變形和流體介質(zhì)物理屬性對液滴遷移速度和液滴表面溫度分布的影響.對于單液滴熱毛細(xì)遷移的非定常數(shù)值模擬也開展了很多相應(yīng)的研究.Yin 等[16]研究了不變形液滴的熱毛細(xì)遷移過程,討論了無量綱參數(shù)及初始條件對遷移的影響.Brady 等[17]分析了三維可變形液滴的熱毛細(xì)遷移過程,并且發(fā)現(xiàn)液池壁面形狀最終影響液滴遷移過程.對于雙液滴熱毛細(xì)遷移定常數(shù)值模擬研究,Loewenberg 和Davis[18]計(jì)算得到了軸對稱排列的雙液滴不同半徑比條件下的熱毛細(xì)遷移速度.Lavrenteva 和Nir[19]討論了重力與熱毛細(xì)力共同作用下的軸對稱雙液滴遷移及相互作用,存在熱邊界層和熱尾跡.對于雙液滴熱毛細(xì)遷移非定常數(shù)值模擬研究.Berejnov 等[20]給出了兩個(gè)液滴中心間距隨時(shí)間的演變,討論了不同半徑比的兩個(gè)液滴熱毛細(xì)遷移過程及相互作用.Yin 和Li[21]研究了具有不同半徑比雙液滴熱毛細(xì)遷移及相互作用,分析了后繼大液滴與前導(dǎo)小液滴的遷移速度.

    與理論研究和數(shù)值模擬相比,液滴遷移的相關(guān)實(shí)驗(yàn)較少.對于單液滴熱毛細(xì)遷移,Barton 和Subramanian[22]實(shí)驗(yàn)發(fā)現(xiàn),遷移速度與液滴半徑成正比,與YGB 理論預(yù)測的趨勢相同.對于雙液滴熱毛細(xì)遷移,H?hnel等[23]觀察到大液滴超越小液滴的過程.Kang 等[24]實(shí)驗(yàn)研究了雙液滴熱毛細(xì)遷移及其相互作用,發(fā)現(xiàn)前導(dǎo)小液滴出現(xiàn)傾斜“8”字的遷移軌跡;近年利用數(shù)字全息干涉方法實(shí)時(shí)觀測了單液滴和雙液滴遷移過程[25-26],從全新的視角研究溫度場對液滴遷移速度的影響.Vincent 和Delville[27]實(shí)驗(yàn)研究了微流道中通過局部激光限制生長液滴前沿前進(jìn)的熱毛細(xì)遷移,發(fā)現(xiàn)激光能夠控制微流道中液滴的流動(dòng),起到“軟門”(soft gate)的作用.

    本文工作利用COMSOL 軟件對激光驅(qū)動(dòng)單液滴遷移進(jìn)行的仿真計(jì)算,研究了影響單液滴遷移的參數(shù)及單液滴持續(xù)遷移的條件,為地面實(shí)驗(yàn)及空間實(shí)驗(yàn)提供了依據(jù),基于仿真計(jì)算在實(shí)驗(yàn)中利用激光驅(qū)動(dòng)液滴遷移,初步驗(yàn)證激光驅(qū)動(dòng)液滴遷移的可行性,為實(shí)現(xiàn)對液滴的有效驅(qū)動(dòng)和控制打下基礎(chǔ).

    1 驅(qū)動(dòng)液滴遷移原理和仿真計(jì)算模型

    1.1 激光驅(qū)動(dòng)液滴遷移原理

    在微重力環(huán)境中,溫度場引發(fā)的液滴界面張力不均勻,成為驅(qū)動(dòng)液滴運(yùn)動(dòng)的主要作用場,由于界面張力一般會隨著溫度升高而減小,這種界面溫度的不均勻?qū)е陆缑鎻埩Φ牟痪鶆?使得液滴界面上受力不均勻,從而驅(qū)動(dòng)界面兩側(cè)的流體發(fā)生流動(dòng),使液滴產(chǎn)生熱毛細(xì)遷移運(yùn)動(dòng).不考慮重力的情況下激光沿液滴直徑豎直向下照射時(shí),母液與液滴均會吸收激光能量,激光與液滴界面接觸上下兩點(diǎn)的溫度分別為T1和T2,母液溫度為T0,T1>T2>T0,由于激光照射在液滴界面附近產(chǎn)生了溫度差,導(dǎo)致液滴界面附近產(chǎn)生表面張力梯度,同時(shí)引起母液流動(dòng),在界面熱毛細(xì)力與母液流動(dòng)的共同作用下驅(qū)動(dòng)液滴遷移.原理如圖1 所示.

    圖1 激光驅(qū)動(dòng)液滴遷移原理圖Fig.1 Schematic of laser-driven droplet migration

    1.2 仿真計(jì)算模型

    利用COMSOL 軟件進(jìn)行仿真計(jì)算,激光沿液滴直徑豎直向下照射,二維模型示意圖見圖2.激光能量由顏色表示,紅色表能量最強(qiáng),藍(lán)色表能量最弱.光強(qiáng)在光束橫截面上為高斯分布,符合實(shí)際激光器激光光束光強(qiáng)分布情況.由于液滴與母液不混溶,二者之間界面的存在導(dǎo)致了表面張力梯度的產(chǎn)生,使得液滴在母液中遷移,則選用了動(dòng)網(wǎng)格為兩種不混溶流體的層流建模,可以跟蹤流體-流體界面的位置.進(jìn)行仿真模擬時(shí)液滴界面附近的網(wǎng)格相比其他區(qū)域需要進(jìn)行細(xì)化,如圖3 所示.液滴與母液的界面處設(shè)置了隨溫度線性減小的表面張力,符合實(shí)際情況;母液四邊設(shè)置無滑移條件,模擬母液與液池壁面接觸的邊界條件.計(jì)算時(shí)采用了液滴不變形假設(shè).為了研究母液尺寸、液滴直徑、母液和液滴對激光的吸收系數(shù)[28-29]等參數(shù)對液滴遷移的影響,本文選擇了母液尺寸為寬50 mm 高50 mm 及寬50 mm 高20 mm,液滴直徑為3,4,4.5,5 和6 mm,吸收系數(shù)0.06 m-1和47 m-1進(jìn)行了詳細(xì)的仿真計(jì)算和分析.

    圖2 激光驅(qū)動(dòng)液滴遷移模型Fig.2 Laser-driven droplet migration model

    圖3 網(wǎng)格劃分Fig.3 Mesh division

    ρi表示液體密度,νi表示液體運(yùn)動(dòng)學(xué)黏性系數(shù),μi表示液體動(dòng)力學(xué)黏性系數(shù),кi表示液體熱擴(kuò)散系數(shù),Λi表示液體熱傳導(dǎo)系數(shù),Г表示母液中的溫度梯度,R 表示液滴半徑,σT表示界面張力隨溫度變化系數(shù),Ar表示液滴直徑與母液寬度之比,下標(biāo)1 和2 分別代表母液及液滴.

    Re 表示慣性力與黏性力之比,定義如下

    馬朗戈尼數(shù)(Ma)表示熱對流傳輸與熱擴(kuò)散之比,定義如下

    利用密度匹配法,母液選擇30 cSt 硅油,液滴選擇水和乙醇的混合液,減小重力引起的浮力效應(yīng),故模擬計(jì)算時(shí)選用了相應(yīng)的物性參數(shù),見表1.

    表1 物性參數(shù)Table 1 Physical property parameters

    2 模擬結(jié)果與分析

    2.1 時(shí)間步長無關(guān)性驗(yàn)證

    對于瞬態(tài)數(shù)值模擬,需要確定計(jì)算中時(shí)間步長與計(jì)算結(jié)果之間無關(guān)聯(lián)性,即要進(jìn)行時(shí)間步長無關(guān)性驗(yàn)證.對同一工況選取時(shí)間步長Δt 分別為0.3 s,0.5 s 和1.0 s 的計(jì)算.得到不同時(shí)間步長下同一工況液滴界面附近母液流速及液滴速度圖,圖4 中藍(lán)線均表示時(shí)間步長1.0 s 時(shí)的速度,綠色中空方塊線均表示時(shí)間步長0.5 s 時(shí)的速度,紅色中空圓形線均表示時(shí)間步長0.3 s 時(shí)的速度.

    由圖4 可看出,時(shí)間步長選用0.3 s 和0.5 s,速度曲線重合度極高.時(shí)間步長選用0.5 s 和1.0 s 時(shí)液滴上界面附近母液流動(dòng)最大速度與有一定差別,但其他兩種速度沒有太大差異.由于步長取得太小會大大增加計(jì)算周期,取得太大又會影響計(jì)算精度,綜合考慮,在本文計(jì)算中取計(jì)算的時(shí)間步長大部分為1.0 s,當(dāng)選用時(shí)間步長1.0 s 時(shí)速度曲線平滑度較差的情況下選用時(shí)間步長為0.5 s.

    2.2 母液流動(dòng)對液滴遷移的影響

    當(dāng)母液與液滴吸收系數(shù)很小時(shí),激光穿過母液射到液滴界面上時(shí),溫度變化很小,引起非常小的界面張力變化,此時(shí)液滴行為主要受母液流動(dòng)的影響.選取母液和液滴的吸收系數(shù)均為0.06 m-1,選擇液滴直徑與母液寬度之比Ar均為50/300,分別研究液滴處于不同初始位置時(shí)的母液流動(dòng)狀態(tài)對液滴遷移行為的影響,探討母液流動(dòng)的作用機(jī)理.

    液滴初始位置在母液上部,液滴直徑R=5 mm,Γmax=1/30(K·mm-1).計(jì)算發(fā)現(xiàn)(見圖5),液滴界面附近母液流動(dòng)較快,且液滴下方母液流速小于液滴上方母液流速;黑色箭頭表示流場施加于液滴的壓力,由伯努利原理可知,此時(shí)液滴所受壓力的合力向上,用紅色箭頭表示,因此液滴向上遷移但遷移速度較慢.圖6 給出液滴初始位置在母液上部時(shí)液滴界面附近母液流速及液滴速度,其中綠線表示液滴上邊界附近母液的最大流速,紅線表示液滴下邊界附近母液的最大流速,藍(lán)線表示液滴速度.圖7 給出液滴界面附近的壓力圖,藍(lán)線表示上邊界附近最大壓力,綠線表示下邊界附近最大壓力.由于吸收系數(shù)小,溫度梯度引起的熱毛細(xì)力較小,母液流動(dòng)導(dǎo)致液滴上邊界的最大壓力小于液滴下邊界的最大壓力,合力方向向上,液滴向上遷移.

    圖5 液滴初始位置在母液上部時(shí)的流線圖Fig.5 Streamline when the initial droplet position droplet is in the upper part of the mother liquor

    圖6 液滴初始位置在母液上部液滴界面附近母液流速及液滴速度圖Fig.6 Velocity of droplet and mother liquor around droplet interface when initial droplet position is in the upper part of the mother liquor

    圖7 液滴初始位置在母液上部液滴界面附近壓力圖Fig.7 Pressure variation around droplet interface when initial droplet position is in the upper part of the mother liquor

    當(dāng)液滴初始位置在母液正中,Γmax=2.61/30(K·mm-1).計(jì)算發(fā)現(xiàn)(見圖8),液滴周圍形成對稱的流場,下方母液流速等于液滴上方母液流速,流場施加于液滴的合力為0,因此液滴始終處于初始位置.圖9 給出液滴初始位置在母液正中時(shí)液滴界面附近母液流速及液滴速度,其中綠線表示液滴上邊界附近母液的最大流速,紅線表示液滴下邊界附近母液的最大流速.圖10 給出液滴界面附近的壓力圖,藍(lán)線表示上邊界附近最大壓力,綠線表示下邊界附近最大壓力.液滴上邊界的最大壓力等于液滴下邊界的最大壓力,液滴受到的壓力合力為零,液滴始終處于初始位置.

    圖8 液滴初始位置在母液正中時(shí)的流線圖Fig.8 Streamline when the initial droplet position droplet is in the middle of the mother liquor

    圖9 液滴初始位置在母液正中液滴界面附近母液流速圖Fig.9 Velocity of droplet and mother liquor around droplet interface when initial droplet position is in the middle of the mother liquor

    圖10 液滴初始位置在母液正中液滴界面附近壓力圖Fig.10 Pressure variation around droplet interface when initial droplet position is in the middle of the mother liquor

    當(dāng)液滴初始位置在母液下部,Γmax=3.64/30(K·mm-1).計(jì)算發(fā)現(xiàn)(見圖11),液滴上方母液流速小于液滴下方母液流速;黑色箭頭表示流場施加于液滴的壓力,由伯努利原理可知,此時(shí)液滴所受壓力的合力向下,用紅色箭頭表示,因此液滴向下微微遷移.圖12 給出液滴初始位置在母液下部時(shí)液滴界面附近母液流速及液滴速度,其中綠線表示液滴上邊界附近母液的最大流速,紅線表示液滴下邊界附近母液的最大流速,藍(lán)線表示液滴速度.圖13 給出液滴界面附近的壓力圖,藍(lán)線表示上邊界附近最大壓力,綠線表示下邊界附近最大壓力.在吸收系數(shù)較小的情況下,溫度梯度引起的熱毛細(xì)力較小,母液流動(dòng)對液滴遷移的影響占主導(dǎo)地位.液滴上邊界的最大壓力大于液滴下邊界的最大壓力,液滴受到的壓力合力方向向下,液滴向下遷移.

    圖11 液滴初始位置在母液下部時(shí)的流線圖Fig.11 Streamline when the initial droplet position droplet is in the lower part of the mother liquor

    圖12 液滴初始位置在母液下部液滴界面附近母液 流速及液滴速度圖Fig.12 Velocity of droplet and mother liquor around droplet interface when initial droplet position is in the lower part of the mother liquor

    圖13 液滴初始位置在母液下部液滴界面附近壓力Fig.13 Pressure variation around droplet interface when initial droplet position is in the lower part of the mother liquor

    在母液、液滴吸收系數(shù)均較小的情況下,界面張力梯度很小,熱毛細(xì)驅(qū)動(dòng)力很小,此時(shí)母液流動(dòng)對液滴在豎直方向上運(yùn)動(dòng)的影響占主導(dǎo)地位.這種情況下,液滴向上或向下運(yùn)動(dòng)取決于液滴的初始位置.

    2.3 界面張力對液滴遷移的影響

    當(dāng)母液吸收系數(shù)小而液滴吸收系數(shù)大時(shí),液滴界面上形成較大的溫度梯度,影響液滴遷移行為的主要因素是由該溫度梯度引起的界面張力梯度.選取母液吸收系數(shù)為0.06 m-1,液滴吸收系數(shù)為47 m-1進(jìn)行計(jì)算,當(dāng)液滴初始位置在母液下部,Ar=1/10,液滴周圍 Γmax=50/2.5(K·mm-1) .發(fā)現(xiàn)液滴向上遷移,并于6 s 時(shí)接觸母液頂部,見圖14.初始時(shí)刻母液中有4 個(gè)渦存在,與母液和液滴吸收系數(shù)均很小時(shí)相比,液滴下方母液中產(chǎn)生的渦較小,且在液滴遷移過程中母液中演化成沿y 方向?qū)ΨQ的兩個(gè)渦.圖15 給出液滴初始位置在母液下部時(shí)液滴界面附近母液流速及液滴速度,其中綠線表示液滴上邊界附近母液的最大流速,紅線表示液滴下邊界附近母液的最大流速,藍(lán)線表示液滴速度.由速度圖可以看出,液滴上邊界母液流速大于液滴下邊界母液流速,當(dāng)液滴即將接觸母液頂部時(shí),由于邊界的限制,液滴速度及液滴上下邊界母液流速均會驟減.

    圖14 液滴初始位置在母液下部時(shí)的流線圖Fig.14 Streamline when the initial droplet position droplet is in the lower part of the mother liquor

    圖15 液滴初始位置在母液下部液滴界面附近母液流速及液滴速度圖Fig.15 Velocity of droplet and mother liquor around droplet interface when initial droplet position is in the lower part of the mother liquor

    當(dāng)母液吸收系數(shù)小而液滴吸收系數(shù)大時(shí),即使液滴初始位置在母液下部,液滴仍將向上遷移.因此,實(shí)現(xiàn)激光驅(qū)動(dòng)液滴遷移,需要選取母液吸收系數(shù)小而液滴吸收系數(shù)大的液滴遷移模型.

    2.4 液滴直徑與母液寬度比Ar 對液滴遷移的影響

    液滴直徑與母液寬度比Ar也是液滴遷移的重要影響因素,因?yàn)槟敢毫鲃?dòng)和界面張力的耦合作用控制液滴遷移行為.選取與2.3 小節(jié)相同的吸收系數(shù),即母液吸收系數(shù)為0.06 m-1以及液滴吸收系數(shù)為47 m-1,通過選取不同的Ar計(jì)算流場狀態(tài),考察Ar對液滴遷移的影響,進(jìn)而探索母液流動(dòng)和界面張力的耦合作用機(jī)理.

    當(dāng)母液吸收系數(shù)遠(yuǎn)小于液滴吸收系數(shù)時(shí),在液滴直徑與母液寬度比Ar<1/10 的情況下,液滴位于母液下部時(shí)仍會向上遷移,但此時(shí)液滴會出現(xiàn)往復(fù)運(yùn)動(dòng).通過速度圖發(fā)現(xiàn)(見圖16 和圖17),液滴上邊界處母液流動(dòng)速度遠(yuǎn)大于液滴遷移速度,母液對流影響了液滴持續(xù)向上遷移的趨勢,造成液滴失穩(wěn),產(chǎn)生往復(fù)運(yùn)動(dòng).但從液滴球心距母液底端距離的趨勢中發(fā)現(xiàn)液滴球心距離母液頂端的距離在減小,推測計(jì)算時(shí)間足夠長的情況下,液滴最終會到達(dá)母液頂端.

    圖16 不同Ar 時(shí)液滴球心距母液底端距離Fig.16 Distance variation from the center of the droplet sphere to the bottom of the mother liquor with different Ar

    圖17 液滴初始位置在母液下部液滴界面附近母液流速及液滴速度圖Fig.17 Velocity of droplet and mother liquor around droplet interface when initial droplet position is in the lower part of the mother liquor

    當(dāng)母液吸收系數(shù)遠(yuǎn)小于液滴吸收系數(shù)時(shí),在1/10≤Ar<1 的情況下,當(dāng)液滴初始位置處于母液下部時(shí),在初始階段,由于此時(shí)液滴界面周圍溫度梯度較小,液滴會首先向下遷移,由于液滴吸收系數(shù)遠(yuǎn)大于母液吸收系數(shù),較短時(shí)間后熱毛細(xì)驅(qū)動(dòng)力起主導(dǎo)作用,液滴持續(xù)向上遷移并很快到達(dá)母液頂端,見圖18和圖19.

    圖18 不同Ar 時(shí)液滴球心距母液底端距離Fig.18 Distance variation from the center of the droplet sphere to the bottom of the mother liquor with different Ar

    圖19 液滴初始位置在母液下部液滴界面附近母液流速及液滴速度圖Fig.19 Velocity of droplet and mother liquor around droplet interface when initial droplet position is in the lower part of the mother liquor

    仿真結(jié)果與原理分析一致,由于溫度梯度與液滴界面的存在,母液中最初會形成4 個(gè)渦.在母液吸收系數(shù)小、液滴吸收系數(shù)大的情況下,液滴界面周圍的溫度梯度較大.當(dāng)母液吸收系數(shù)小而液滴吸收系數(shù)大時(shí),激光照射液滴使液滴界面上形成較大的溫度梯度而引起較大的界面張力梯度,進(jìn)而驅(qū)動(dòng)液滴遷移.在液滴遷移過程中,母液中的渦逐漸由4 個(gè)發(fā)展成2 個(gè).液滴的初始位置對液滴遷移的方向影響較小,熱毛細(xì)力對液滴遷移方向的影響更大.液滴界面附近溫度梯度越大,液滴的運(yùn)動(dòng)速度越大,液滴周圍母液流速越大.

    3 YGB 理論與仿真結(jié)果的對比

    YGB 理論忽略了流體運(yùn)動(dòng)的黏性項(xiàng)和流體傳熱的對流項(xiàng)等非線性項(xiàng),也不包含時(shí)間項(xiàng),解決的是定常線性問題.但實(shí)際上液滴在遷移過程中,由于液體中各種的物性參數(shù)隨溫度變化,周圍條件的不斷改變使得液滴遷移運(yùn)動(dòng)最終不一定能達(dá)到穩(wěn)定狀態(tài).Subramanian 介于此問題提出了一個(gè)準(zhǔn)穩(wěn)態(tài)假設(shè),如果氣泡遷移足夠慢,使得相對于由氣泡運(yùn)動(dòng)而產(chǎn)生的周圍母液物性參數(shù)變化所花費(fèi)的時(shí)間,速度場和溫度場重新到其穩(wěn)定狀態(tài)所花費(fèi)的時(shí)間非常小,那么可以認(rèn)為此狀態(tài)為準(zhǔn)穩(wěn)態(tài).液滴遷移系統(tǒng)也同樣適用準(zhǔn)穩(wěn)態(tài)假設(shè),如果液滴遷移系統(tǒng)符合準(zhǔn)穩(wěn)態(tài)假設(shè),則仍然可以認(rèn)為液滴最終可以達(dá)到一個(gè)穩(wěn)定狀態(tài),關(guān)于液滴遷移的研究大都是建立在這個(gè)假設(shè)的基礎(chǔ)上的,在分析過程中所有的物性參數(shù)都以常數(shù)來表示,不隨溫度而變化.

    通過YGB 理論計(jì)算得到的氣泡遷移速度一般也可以應(yīng)用在液滴熱毛細(xì)遷移中,理論遷移速度為

    仿真計(jì)算得到的液滴遷移速度Vmax與YGB 理論計(jì)算速度VYGB見表2.

    表2 不同直徑不同溫度梯度液滴遷移仿真計(jì)算速度與YGB 理論遷移速度及Re,Ma 對照表Table 2 Size of droplets,simulation calculation migration velocity,YGB model velocity,Re and Ma

    在液滴直徑和母液尺寸相同的情況下,溫度梯度越大,液滴遷移速度越大,與YGB 理論預(yù)測一致.YGB 理論考慮的是小Re 和小Ma 情況,仿真計(jì)算考慮了黏性項(xiàng),對流項(xiàng),Ma 并沒有遠(yuǎn)小于1,線性假設(shè)不能成立,能量損耗會使液滴遷移速度降低.當(dāng)液滴吸收系數(shù)遠(yuǎn)大于母液吸收系數(shù)時(shí),液滴界面附近有局部強(qiáng)溫度梯度,但整體溫度梯度遠(yuǎn)小于最大溫度梯度,導(dǎo)致理論遷移速度大于仿真計(jì)算速度.

    4 結(jié)論

    利用COMSOL 軟件計(jì)算微重力環(huán)境下激光驅(qū)動(dòng)液滴遷移過程,發(fā)現(xiàn)當(dāng)母液和液滴對激光的吸收系數(shù)都很小時(shí),母液流動(dòng)對液滴在豎直方向上運(yùn)動(dòng)的影響占主導(dǎo)地位.這種情況下,液滴向上或向下運(yùn)動(dòng)取決于液滴的初始位置.當(dāng)母液吸收系數(shù)遠(yuǎn)小于液滴吸收系數(shù)時(shí),液滴的最終趨勢均為向上遷移,但液滴遷移行為受液滴直徑與母液寬度之比Ar的影響,當(dāng)Ar<1/10 時(shí),液滴會有往復(fù)運(yùn)動(dòng)的行為,當(dāng)1/10≤Ar<1 時(shí),液滴不會出現(xiàn)往復(fù)運(yùn)動(dòng)的現(xiàn)象.

    將仿真結(jié)果與理論結(jié)果對比,發(fā)現(xiàn)在液滴直徑和母液尺寸相同的情況下,溫度梯度越大,液滴遷移速度越大,與YGB 理論一致.由于YGB 理論未考慮非線性項(xiàng),且溫度梯度隨時(shí)間變化,仿真計(jì)算得到的液滴遷移速度小于理論計(jì)算速度,但在母液尺寸相同的情況下,液滴直徑、溫度梯度越大,液滴遷移速度越大,與YGB 理論計(jì)算的趨勢一致.仿真計(jì)算為激光驅(qū)動(dòng)液滴遷移的地面實(shí)驗(yàn)及空間實(shí)驗(yàn)提供了研究基礎(chǔ),基于仿真計(jì)算結(jié)果進(jìn)行了激光器選型,選擇了近紅外光波段激光器,搭建實(shí)驗(yàn)臺進(jìn)行了初步實(shí)驗(yàn),實(shí)現(xiàn)激光驅(qū)動(dòng)液滴遷移,為解決工業(yè)生產(chǎn)問題提供了思路.

    猜你喜歡
    溫度梯度母液毛細(xì)
    “毛細(xì)”政務(wù)號關(guān)停背后
    廉政瞭望(2024年5期)2024-05-26 13:21:07
    溫度對中小跨徑混凝土梁橋應(yīng)力及變形的影響研究
    多孔建筑材料毛細(xì)吸水系數(shù)連續(xù)測量的方法
    能源工程(2020年6期)2021-01-26 00:55:16
    嚴(yán)寒地區(qū)混凝土箱梁實(shí)測溫度梯度分析
    山西建筑(2019年21期)2019-12-02 02:35:08
    溫度梯度場對聲表面波器件影響研究
    電子制作(2018年23期)2018-12-26 01:01:20
    出現(xiàn)憋喘 可能是毛細(xì)支氣管炎!
    媽媽寶寶(2017年3期)2017-02-21 01:22:16
    高滲鹽水霧化吸入治療毛細(xì)支氣管炎的療效觀察
    高速鐵路CRTSⅢ型板式無砟軌道溫度梯度試驗(yàn)研究
    鐵道建筑(2015年3期)2015-12-26 09:41:44
    母液吸氨器安裝高度的探討
    氯乙酸母液再利用技術(shù)最新進(jìn)展
    中國氯堿(2014年11期)2014-02-28 01:05:08
    日韩亚洲欧美综合| 亚洲欧美日韩东京热| ponron亚洲| 久久久久久伊人网av| 男女之事视频高清在线观看| xxxwww97欧美| 色5月婷婷丁香| 99久久精品国产国产毛片| 日韩欧美一区二区三区在线观看| 欧美性猛交╳xxx乱大交人| 日韩亚洲欧美综合| 亚洲在线观看片| 三级男女做爰猛烈吃奶摸视频| 少妇丰满av| 国产精品自产拍在线观看55亚洲| 国产精品精品国产色婷婷| 亚州av有码| 观看免费一级毛片| 亚洲第一电影网av| 人妻久久中文字幕网| 两人在一起打扑克的视频| 亚洲性夜色夜夜综合| 国产一区二区激情短视频| 嫩草影视91久久| 一级黄色大片毛片| 男女啪啪激烈高潮av片| 免费电影在线观看免费观看| 色精品久久人妻99蜜桃| 亚洲欧美日韩高清在线视频| 在线观看av片永久免费下载| 噜噜噜噜噜久久久久久91| 午夜福利视频1000在线观看| 天堂网av新在线| 国产极品精品免费视频能看的| 免费av毛片视频| 亚洲黑人精品在线| 又爽又黄a免费视频| 亚洲av五月六月丁香网| 久久精品国产清高在天天线| 最近最新免费中文字幕在线| 99久久成人亚洲精品观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品sss在线观看| 黄色视频,在线免费观看| 热99re8久久精品国产| 国产精品电影一区二区三区| 久久天躁狠狠躁夜夜2o2o| 久久久久久大精品| 两人在一起打扑克的视频| 亚洲av电影不卡..在线观看| 日韩,欧美,国产一区二区三区 | 色综合站精品国产| 欧美xxxx性猛交bbbb| 啦啦啦韩国在线观看视频| 国产一区二区三区视频了| 丰满人妻一区二区三区视频av| 亚洲五月天丁香| 日韩欧美在线二视频| 欧美最新免费一区二区三区| 他把我摸到了高潮在线观看| 久久久久久久精品吃奶| av黄色大香蕉| 三级国产精品欧美在线观看| 亚洲内射少妇av| 麻豆国产av国片精品| 午夜精品一区二区三区免费看| 成年女人看的毛片在线观看| 国产精品国产三级国产av玫瑰| 在线观看午夜福利视频| 啦啦啦观看免费观看视频高清| 不卡一级毛片| 十八禁网站免费在线| 日本欧美国产在线视频| 搡老熟女国产l中国老女人| 自拍偷自拍亚洲精品老妇| 嫩草影院精品99| xxxwww97欧美| 免费av不卡在线播放| 搡女人真爽免费视频火全软件 | 禁无遮挡网站| 精品久久久久久久人妻蜜臀av| 中文字幕av成人在线电影| 综合色av麻豆| 精品午夜福利视频在线观看一区| 亚洲成人免费电影在线观看| 国产成人福利小说| 最新在线观看一区二区三区| 欧美激情国产日韩精品一区| 久久国产乱子免费精品| 亚洲精品久久国产高清桃花| 老师上课跳d突然被开到最大视频| 在线观看免费视频日本深夜| 麻豆国产97在线/欧美| bbb黄色大片| 天堂av国产一区二区熟女人妻| 成人国产一区最新在线观看| 日日干狠狠操夜夜爽| 九九爱精品视频在线观看| 国产精品亚洲一级av第二区| 国产一区二区激情短视频| 亚洲欧美激情综合另类| 黄色日韩在线| 亚洲av.av天堂| 高清毛片免费观看视频网站| 精品久久久久久久人妻蜜臀av| 亚洲精品一区av在线观看| 亚洲av一区综合| 老司机午夜福利在线观看视频| 国产精品乱码一区二三区的特点| 久久中文看片网| 国产精品嫩草影院av在线观看 | 丰满人妻一区二区三区视频av| 一本一本综合久久| 一个人看视频在线观看www免费| 婷婷色综合大香蕉| 亚洲最大成人手机在线| 欧美日本视频| 高清日韩中文字幕在线| 99热6这里只有精品| 国产成人一区二区在线| 国产精品免费一区二区三区在线| 美女 人体艺术 gogo| 中文资源天堂在线| 国产私拍福利视频在线观看| 久久精品国产亚洲网站| 伊人久久精品亚洲午夜| 99热只有精品国产| 能在线免费观看的黄片| 老司机深夜福利视频在线观看| 3wmmmm亚洲av在线观看| 久久久久九九精品影院| 国产精品伦人一区二区| 日本三级黄在线观看| 真人做人爱边吃奶动态| 12—13女人毛片做爰片一| 日本一本二区三区精品| 变态另类成人亚洲欧美熟女| 波多野结衣高清无吗| 噜噜噜噜噜久久久久久91| av女优亚洲男人天堂| 女的被弄到高潮叫床怎么办 | 日本五十路高清| 我的女老师完整版在线观看| 俄罗斯特黄特色一大片| 午夜老司机福利剧场| 国产中年淑女户外野战色| 精品99又大又爽又粗少妇毛片 | 免费高清视频大片| 亚洲七黄色美女视频| 别揉我奶头~嗯~啊~动态视频| 干丝袜人妻中文字幕| 十八禁国产超污无遮挡网站| 国产精品野战在线观看| 欧美激情久久久久久爽电影| 999久久久精品免费观看国产| 日韩精品有码人妻一区| 色吧在线观看| 免费看美女性在线毛片视频| 好男人在线观看高清免费视频| 中文字幕人妻熟人妻熟丝袜美| 国产免费男女视频| 日日啪夜夜撸| 国产91精品成人一区二区三区| 免费高清视频大片| 久久精品国产清高在天天线| 嫩草影院新地址| 高清毛片免费观看视频网站| 久久久久性生活片| 中文字幕人妻熟人妻熟丝袜美| 国产激情偷乱视频一区二区| 国产白丝娇喘喷水9色精品| 国产午夜福利久久久久久| 午夜福利欧美成人| 在线免费十八禁| 色5月婷婷丁香| 亚洲综合色惰| 成年女人毛片免费观看观看9| 三级男女做爰猛烈吃奶摸视频| 久久久久久国产a免费观看| 女的被弄到高潮叫床怎么办 | xxxwww97欧美| 999久久久精品免费观看国产| 国产精品精品国产色婷婷| 中文字幕免费在线视频6| 亚洲电影在线观看av| 亚洲精品国产成人久久av| 99热这里只有精品一区| 999久久久精品免费观看国产| 欧美性感艳星| 99久久成人亚洲精品观看| 亚洲专区中文字幕在线| 亚洲欧美清纯卡通| 亚洲欧美日韩东京热| 欧美黑人欧美精品刺激| 色视频www国产| 亚洲精品一卡2卡三卡4卡5卡| 十八禁网站免费在线| 久久草成人影院| 中文亚洲av片在线观看爽| 听说在线观看完整版免费高清| 欧美高清性xxxxhd video| 国产精品综合久久久久久久免费| 亚洲精华国产精华精| 成人欧美大片| 国产伦精品一区二区三区视频9| 成人性生交大片免费视频hd| 欧美精品啪啪一区二区三区| 亚洲最大成人av| 日本三级黄在线观看| 欧美+亚洲+日韩+国产| 日本 欧美在线| 国产激情偷乱视频一区二区| 男女做爰动态图高潮gif福利片| 午夜精品久久久久久毛片777| 亚洲最大成人中文| 成人特级av手机在线观看| 欧美高清性xxxxhd video| 大型黄色视频在线免费观看| 亚洲黑人精品在线| 久久久久久久久久成人| 国产黄a三级三级三级人| 99久国产av精品| 欧美日韩黄片免| 岛国在线免费视频观看| 色综合婷婷激情| 久久久久久伊人网av| 99热这里只有是精品在线观看| 三级男女做爰猛烈吃奶摸视频| or卡值多少钱| 国产在线精品亚洲第一网站| 国内揄拍国产精品人妻在线| 欧美色欧美亚洲另类二区| 国产精品亚洲一级av第二区| 日韩av在线大香蕉| 国产午夜精品论理片| 我的女老师完整版在线观看| 99热这里只有精品一区| 少妇人妻一区二区三区视频| 很黄的视频免费| 动漫黄色视频在线观看| 欧美成人a在线观看| 赤兔流量卡办理| 国产伦精品一区二区三区视频9| 精品无人区乱码1区二区| 欧美激情久久久久久爽电影| 91麻豆精品激情在线观看国产| 欧美成人一区二区免费高清观看| 午夜福利成人在线免费观看| 国产探花极品一区二区| 精品人妻视频免费看| 久久人人爽人人爽人人片va| 毛片一级片免费看久久久久 | 日韩,欧美,国产一区二区三区 | 淫秽高清视频在线观看| 在线看三级毛片| a级一级毛片免费在线观看| 日本黄色视频三级网站网址| 久久99热这里只有精品18| 一夜夜www| 特大巨黑吊av在线直播| 中国美白少妇内射xxxbb| 九九久久精品国产亚洲av麻豆| 美女 人体艺术 gogo| 一区二区三区激情视频| 97热精品久久久久久| 欧美人与善性xxx| 天堂动漫精品| 国产高清三级在线| 草草在线视频免费看| 非洲黑人性xxxx精品又粗又长| 久久精品91蜜桃| 最近中文字幕高清免费大全6 | 99久久中文字幕三级久久日本| 国产精品伦人一区二区| 2021天堂中文幕一二区在线观| 嫁个100分男人电影在线观看| 成人高潮视频无遮挡免费网站| 麻豆国产97在线/欧美| 日韩欧美国产在线观看| 黄色女人牲交| 中出人妻视频一区二区| 精品人妻1区二区| 欧美成人性av电影在线观看| netflix在线观看网站| 国产男人的电影天堂91| 精品久久久久久久久久免费视频| 日韩欧美一区二区三区在线观看| 国产久久久一区二区三区| av在线观看视频网站免费| 午夜精品久久久久久毛片777| 国产视频一区二区在线看| 少妇裸体淫交视频免费看高清| 看片在线看免费视频| 亚洲精品色激情综合| 校园春色视频在线观看| 少妇被粗大猛烈的视频| 又爽又黄无遮挡网站| 欧美+亚洲+日韩+国产| 国产精品嫩草影院av在线观看 | 嫩草影院新地址| 免费av观看视频| 床上黄色一级片| 尤物成人国产欧美一区二区三区| 国产蜜桃级精品一区二区三区| 又爽又黄无遮挡网站| 神马国产精品三级电影在线观看| 国产精品美女特级片免费视频播放器| 亚洲一区二区三区色噜噜| 国产精品久久久久久av不卡| 日韩精品中文字幕看吧| 一个人免费在线观看电影| 看片在线看免费视频| 国产亚洲91精品色在线| 观看免费一级毛片| 亚洲av免费在线观看| 黄色丝袜av网址大全| avwww免费| 啦啦啦啦在线视频资源| 国产真实伦视频高清在线观看 | 久久精品夜夜夜夜夜久久蜜豆| aaaaa片日本免费| 国产免费一级a男人的天堂| 午夜福利在线观看吧| 婷婷亚洲欧美| 色播亚洲综合网| 少妇熟女aⅴ在线视频| 舔av片在线| 欧美潮喷喷水| 久久久久久久午夜电影| 天堂av国产一区二区熟女人妻| 精品久久久久久久久av| 亚洲三级黄色毛片| 有码 亚洲区| 日韩大尺度精品在线看网址| 国产精品一区二区三区四区久久| 黄色丝袜av网址大全| 婷婷六月久久综合丁香| 精品久久久久久久久久久久久| 我的老师免费观看完整版| 乱人视频在线观看| 伦精品一区二区三区| 我的女老师完整版在线观看| 免费看av在线观看网站| 亚洲专区中文字幕在线| x7x7x7水蜜桃| 久久热精品热| 国产麻豆成人av免费视频| 精品一区二区三区人妻视频| 少妇猛男粗大的猛烈进出视频 | 亚洲无线观看免费| 黄色视频,在线免费观看| 久久久久久久久中文| 精品久久久久久久久av| 成人一区二区视频在线观看| 床上黄色一级片| 免费人成在线观看视频色| 长腿黑丝高跟| 窝窝影院91人妻| 亚州av有码| 久久久久久久亚洲中文字幕| 国产一区二区三区av在线 | 国产免费一级a男人的天堂| 国产男人的电影天堂91| 成年免费大片在线观看| 1000部很黄的大片| 51国产日韩欧美| 日韩欧美 国产精品| 国产精品久久视频播放| 一边摸一边抽搐一进一小说| 国产黄a三级三级三级人| 成人午夜高清在线视频| 美女高潮的动态| av天堂中文字幕网| 色在线成人网| 男女之事视频高清在线观看| av在线老鸭窝| 日本精品一区二区三区蜜桃| 日韩一区二区视频免费看| 免费高清视频大片| 赤兔流量卡办理| 亚洲最大成人av| 亚洲国产日韩欧美精品在线观看| 精品免费久久久久久久清纯| 国语自产精品视频在线第100页| 69av精品久久久久久| 看十八女毛片水多多多| 亚洲va在线va天堂va国产| 亚洲成人免费电影在线观看| 欧美色视频一区免费| 日韩欧美 国产精品| 人妻少妇偷人精品九色| 亚洲精品一区av在线观看| 啦啦啦啦在线视频资源| 一a级毛片在线观看| 校园春色视频在线观看| 高清日韩中文字幕在线| 亚洲人成网站高清观看| 琪琪午夜伦伦电影理论片6080| 日韩欧美在线乱码| 日本爱情动作片www.在线观看 | av在线蜜桃| 国产人妻一区二区三区在| 一区二区三区高清视频在线| aaaaa片日本免费| 精品久久久久久,| 女人被狂操c到高潮| 国产真实伦视频高清在线观看 | 十八禁网站免费在线| 美女被艹到高潮喷水动态| 色哟哟·www| 中文字幕高清在线视频| 搡老妇女老女人老熟妇| 一级黄色大片毛片| 国产欧美日韩一区二区精品| 国产乱人视频| 老司机午夜福利在线观看视频| 少妇丰满av| 国产av一区在线观看免费| 久久婷婷人人爽人人干人人爱| 国产又黄又爽又无遮挡在线| 精品久久国产蜜桃| 欧美激情国产日韩精品一区| 嫩草影院入口| 日本熟妇午夜| 无遮挡黄片免费观看| 国产爱豆传媒在线观看| 啪啪无遮挡十八禁网站| 亚洲人成网站高清观看| 日本与韩国留学比较| 免费看美女性在线毛片视频| 亚洲国产精品合色在线| 毛片一级片免费看久久久久 | 日日撸夜夜添| 亚洲欧美日韩高清在线视频| 舔av片在线| 国产精品人妻久久久久久| av在线天堂中文字幕| 日韩一区二区视频免费看| 国产91精品成人一区二区三区| 一级毛片久久久久久久久女| 亚洲av第一区精品v没综合| 国产欧美日韩精品一区二区| 国产蜜桃级精品一区二区三区| 久久久国产成人精品二区| 国产高清视频在线播放一区| 日韩高清综合在线| 不卡视频在线观看欧美| 九九热线精品视视频播放| 午夜精品一区二区三区免费看| 久久精品国产清高在天天线| 91久久精品国产一区二区成人| 少妇熟女aⅴ在线视频| 国产精品久久视频播放| 国产精品亚洲一级av第二区| 国产人妻一区二区三区在| 男女啪啪激烈高潮av片| 亚洲av不卡在线观看| 91麻豆精品激情在线观看国产| 欧美最黄视频在线播放免费| 亚洲狠狠婷婷综合久久图片| 老熟妇乱子伦视频在线观看| 免费观看的影片在线观看| 成人毛片a级毛片在线播放| 国产亚洲欧美98| 热99在线观看视频| 黄色丝袜av网址大全| 99在线人妻在线中文字幕| 中国美白少妇内射xxxbb| 日韩av在线大香蕉| 国产精品一区二区性色av| 97超级碰碰碰精品色视频在线观看| 一进一出抽搐gif免费好疼| 国产精华一区二区三区| 人人妻,人人澡人人爽秒播| 成年女人毛片免费观看观看9| 他把我摸到了高潮在线观看| 国产精品电影一区二区三区| 黄色视频,在线免费观看| 国产乱人伦免费视频| 日韩亚洲欧美综合| 国产精品无大码| 国产精品久久久久久久久免| 深夜a级毛片| 免费大片18禁| 国产淫片久久久久久久久| 在线观看免费视频日本深夜| 一区二区三区四区激情视频 | 亚洲18禁久久av| 精品人妻一区二区三区麻豆 | 日韩人妻高清精品专区| 给我免费播放毛片高清在线观看| 亚洲成人久久性| 很黄的视频免费| 能在线免费观看的黄片| 美女大奶头视频| 国产av一区在线观看免费| 国产一区二区激情短视频| 天天一区二区日本电影三级| 俺也久久电影网| 国产精品久久久久久亚洲av鲁大| 亚洲avbb在线观看| 有码 亚洲区| 欧美精品啪啪一区二区三区| 真人一进一出gif抽搐免费| 村上凉子中文字幕在线| 观看免费一级毛片| 深爱激情五月婷婷| 99久久九九国产精品国产免费| 国产精品嫩草影院av在线观看 | 国产麻豆成人av免费视频| 一a级毛片在线观看| 啦啦啦观看免费观看视频高清| 日韩欧美免费精品| 动漫黄色视频在线观看| 欧美中文日本在线观看视频| 国产主播在线观看一区二区| 日韩欧美在线二视频| 亚洲av.av天堂| 乱码一卡2卡4卡精品| 88av欧美| 淫妇啪啪啪对白视频| 色吧在线观看| 国产大屁股一区二区在线视频| 久久久久久久久大av| 黄色丝袜av网址大全| 国产私拍福利视频在线观看| 乱码一卡2卡4卡精品| 亚洲精品色激情综合| 日韩大尺度精品在线看网址| 精品久久久久久久久av| 99热精品在线国产| 日本免费a在线| 永久网站在线| 99久国产av精品| 十八禁网站免费在线| 亚洲va日本ⅴa欧美va伊人久久| 成年人黄色毛片网站| 成人三级黄色视频| 亚洲三级黄色毛片| 久久久久久久精品吃奶| 精品人妻偷拍中文字幕| 日本三级黄在线观看| 久久久久久大精品| 在线观看免费视频日本深夜| 成年免费大片在线观看| 亚洲自偷自拍三级| 特级一级黄色大片| 性色avwww在线观看| 两个人的视频大全免费| 国产淫片久久久久久久久| 日韩一本色道免费dvd| 99热6这里只有精品| 国产午夜精品久久久久久一区二区三区 | 国产不卡一卡二| 免费观看精品视频网站| 色播亚洲综合网| 亚洲精品亚洲一区二区| 久久久久精品国产欧美久久久| 听说在线观看完整版免费高清| 国产精品久久久久久久电影| 国产亚洲欧美98| 深夜精品福利| 日韩精品中文字幕看吧| 三级毛片av免费| 国产精品一及| 国产探花极品一区二区| 国产人妻一区二区三区在| 国产高清激情床上av| 亚洲成av人片在线播放无| 一本久久中文字幕| 午夜福利18| 大又大粗又爽又黄少妇毛片口| 国产伦人伦偷精品视频| 国产精品综合久久久久久久免费| 亚洲avbb在线观看| 狂野欧美激情性xxxx在线观看| 变态另类成人亚洲欧美熟女| 成年免费大片在线观看| 国产亚洲91精品色在线| 免费无遮挡裸体视频| 无人区码免费观看不卡| 日韩欧美三级三区| 亚洲国产精品合色在线| 午夜激情欧美在线| 久久久久久久久久黄片| 免费看光身美女| 日韩 亚洲 欧美在线| 在线观看一区二区三区| 看片在线看免费视频| 亚洲va在线va天堂va国产| 一级a爱片免费观看的视频| 免费在线观看日本一区| 91久久精品国产一区二区成人| 国产精品久久久久久久电影| 国产精品一区www在线观看 | 亚洲一区高清亚洲精品| 国产免费av片在线观看野外av| 亚洲人成网站在线播| 精品不卡国产一区二区三区| 亚洲成人久久爱视频| 久久精品人妻少妇| 亚洲av不卡在线观看| 天堂√8在线中文| 国产精品不卡视频一区二区| 亚洲国产精品成人综合色| 两个人的视频大全免费| 国产精品久久久久久久久免| 亚洲国产欧洲综合997久久,| 搡老熟女国产l中国老女人| 午夜老司机福利剧场| 日韩精品青青久久久久久| 欧美国产日韩亚洲一区| АⅤ资源中文在线天堂| 欧美绝顶高潮抽搐喷水| 亚洲av免费高清在线观看|