• <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
    欧美日本亚洲视频在线播放| 99精品在免费线老司机午夜| 美女黄网站色视频| 婷婷亚洲欧美| 啦啦啦观看免费观看视频高清| 麻豆精品久久久久久蜜桃| 久久99热这里只有精品18| 亚洲最大成人中文| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品成人综合色| 丰满人妻一区二区三区视频av| 97碰自拍视频| 国产精品野战在线观看| 性色avwww在线观看| 精品乱码久久久久久99久播| 搡老岳熟女国产| 99热这里只有是精品50| 一本久久中文字幕| 国产一区二区亚洲精品在线观看| 大型黄色视频在线免费观看| 伦精品一区二区三区| 亚洲av不卡在线观看| 日本一二三区视频观看| 国产一区二区三区在线臀色熟女| 色噜噜av男人的天堂激情| 午夜福利在线观看免费完整高清在 | 99久国产av精品国产电影| 十八禁网站免费在线| .国产精品久久| 国产色爽女视频免费观看| 中文字幕精品亚洲无线码一区| 免费搜索国产男女视频| 我的女老师完整版在线观看| 3wmmmm亚洲av在线观看| 91av网一区二区| 日本五十路高清| 午夜福利在线在线| 久久鲁丝午夜福利片| 亚洲美女搞黄在线观看 | 啦啦啦观看免费观看视频高清| 麻豆精品久久久久久蜜桃| a级毛色黄片| 男人和女人高潮做爰伦理| 欧美zozozo另类| 欧美日韩乱码在线| 69av精品久久久久久| 久久久久久久久久黄片| 日日摸夜夜添夜夜添小说| 性欧美人与动物交配| 麻豆久久精品国产亚洲av| 熟女电影av网| 黄色一级大片看看| 免费看光身美女| 看免费成人av毛片| 伦精品一区二区三区| 午夜免费男女啪啪视频观看 | 少妇熟女aⅴ在线视频| a级一级毛片免费在线观看| av专区在线播放| 18禁裸乳无遮挡免费网站照片| 一a级毛片在线观看| 99久久九九国产精品国产免费| 老女人水多毛片| 国产精品一区www在线观看| 国产亚洲91精品色在线| 国产一级毛片七仙女欲春2| 国产高清三级在线| 在线免费观看不下载黄p国产| 舔av片在线| 国产男靠女视频免费网站| 国产伦一二天堂av在线观看| 99久国产av精品| 日韩高清综合在线| 我的女老师完整版在线观看| 插逼视频在线观看| 亚洲欧美日韩高清在线视频| 国产日本99.免费观看| 日日摸夜夜添夜夜爱| 成人国产麻豆网| av福利片在线观看| 亚洲国产精品合色在线| 成人一区二区视频在线观看| 久久久国产成人精品二区| 日韩av在线大香蕉| 热99在线观看视频| 国产亚洲精品久久久久久毛片| 免费一级毛片在线播放高清视频| 免费看av在线观看网站| 校园人妻丝袜中文字幕| 男人舔女人下体高潮全视频| 色综合色国产| 一级毛片aaaaaa免费看小| 简卡轻食公司| 国产一区二区激情短视频| 日本五十路高清| 夜夜看夜夜爽夜夜摸| 亚洲丝袜综合中文字幕| 少妇人妻精品综合一区二区 | 免费观看在线日韩| 国产精品99久久久久久久久| 日本黄色片子视频| 亚洲av熟女| 女同久久另类99精品国产91| 免费黄网站久久成人精品| 偷拍熟女少妇极品色| 在线观看一区二区三区| 亚洲国产高清在线一区二区三| 亚洲精华国产精华液的使用体验 | 欧美潮喷喷水| 少妇熟女欧美另类| 校园人妻丝袜中文字幕| 亚洲内射少妇av| 国内揄拍国产精品人妻在线| 99视频精品全部免费 在线| 大又大粗又爽又黄少妇毛片口| 日本在线视频免费播放| 麻豆成人午夜福利视频| a级毛片a级免费在线| 黄片wwwwww| 日本五十路高清| av天堂中文字幕网| 亚洲欧美日韩高清在线视频| 欧美丝袜亚洲另类| 性色avwww在线观看| 国产亚洲91精品色在线| 人人妻人人澡欧美一区二区| 免费观看在线日韩| 日日啪夜夜撸| 神马国产精品三级电影在线观看| 男插女下体视频免费在线播放| 18禁在线播放成人免费| 国产亚洲精品久久久久久毛片| av在线老鸭窝| 麻豆精品久久久久久蜜桃| 国产麻豆成人av免费视频| 性插视频无遮挡在线免费观看| 亚洲18禁久久av| 日韩成人av中文字幕在线观看 | 精品久久久久久久末码| 免费av毛片视频| 色哟哟哟哟哟哟| 久久久久久久久久成人| 国产激情偷乱视频一区二区| 欧美3d第一页| 给我免费播放毛片高清在线观看| 久久天躁狠狠躁夜夜2o2o| 噜噜噜噜噜久久久久久91| 在线观看一区二区三区| 日韩一本色道免费dvd| 日本-黄色视频高清免费观看| 亚洲精品一卡2卡三卡4卡5卡| 婷婷精品国产亚洲av| 三级男女做爰猛烈吃奶摸视频| 精品不卡国产一区二区三区| 22中文网久久字幕| 男女下面进入的视频免费午夜| 天堂动漫精品| 观看美女的网站| 午夜福利在线观看吧| av女优亚洲男人天堂| 成人二区视频| 黄色欧美视频在线观看| 免费av观看视频| 人人妻人人澡人人爽人人夜夜 | 日本撒尿小便嘘嘘汇集6| 韩国av在线不卡| 日本黄色片子视频| 国产亚洲欧美98| 亚洲人成网站在线播| 99久久成人亚洲精品观看| 亚洲国产精品久久男人天堂| 特大巨黑吊av在线直播| 成人鲁丝片一二三区免费| 欧美高清性xxxxhd video| 国产乱人偷精品视频| 国产色爽女视频免费观看| 美女被艹到高潮喷水动态| 精品久久久久久久人妻蜜臀av| 免费人成在线观看视频色| 国产伦在线观看视频一区| 麻豆久久精品国产亚洲av| 高清日韩中文字幕在线| 成人午夜高清在线视频| 亚洲高清免费不卡视频| 99久国产av精品国产电影| 毛片女人毛片| 久久国内精品自在自线图片| 在线免费十八禁| 一进一出抽搐动态| 欧美色欧美亚洲另类二区| 最新在线观看一区二区三区| 人妻久久中文字幕网| 在线播放国产精品三级| 男女那种视频在线观看| 国产探花在线观看一区二区| 特级一级黄色大片| 国产白丝娇喘喷水9色精品| 国产熟女欧美一区二区| 女同久久另类99精品国产91| 亚洲激情五月婷婷啪啪| 一级毛片久久久久久久久女| 国产伦精品一区二区三区视频9| 少妇丰满av| 又粗又爽又猛毛片免费看| 国产不卡一卡二| 九九久久精品国产亚洲av麻豆| 免费高清视频大片| 卡戴珊不雅视频在线播放| 午夜爱爱视频在线播放| 午夜福利成人在线免费观看| 人人妻人人看人人澡| 黄色一级大片看看| 老熟妇乱子伦视频在线观看| 亚洲性久久影院| 亚洲电影在线观看av| 日韩av不卡免费在线播放| 一级a爱片免费观看的视频| 久久精品国产清高在天天线| 国产精品免费一区二区三区在线| 床上黄色一级片| 午夜亚洲福利在线播放| 久久午夜亚洲精品久久| 亚洲欧美精品综合久久99| 老女人水多毛片| 欧美日本视频| 波多野结衣高清无吗| 日韩欧美三级三区| 不卡视频在线观看欧美| 日韩欧美 国产精品| 亚洲美女黄片视频| 国产女主播在线喷水免费视频网站 | 99热只有精品国产| 男人狂女人下面高潮的视频| 国产成人影院久久av| 久久精品夜色国产| 嫩草影院新地址| 99九九线精品视频在线观看视频| 亚洲av第一区精品v没综合| a级毛片a级免费在线| 欧美zozozo另类| 亚洲av熟女| 免费搜索国产男女视频| 亚洲国产色片| 成人性生交大片免费视频hd| 精品99又大又爽又粗少妇毛片| 99久国产av精品| 亚洲国产精品国产精品| 国产综合懂色| 国产色爽女视频免费观看| 国产精品国产高清国产av| 久久鲁丝午夜福利片| 欧美+日韩+精品| 亚洲av五月六月丁香网| 婷婷亚洲欧美| 一区二区三区高清视频在线| 不卡一级毛片| 国产爱豆传媒在线观看| 欧美日韩在线观看h| 亚洲精品日韩av片在线观看| 色吧在线观看| 久久中文看片网| 哪里可以看免费的av片| 亚洲av一区综合| 观看美女的网站| 免费人成视频x8x8入口观看| 六月丁香七月| 亚洲va在线va天堂va国产| 成年免费大片在线观看| 国产精品伦人一区二区| 中文字幕免费在线视频6| 国内久久婷婷六月综合欲色啪| 国产在视频线在精品| 午夜免费男女啪啪视频观看 | 热99在线观看视频| 日韩制服骚丝袜av| 亚洲av免费在线观看| 少妇裸体淫交视频免费看高清| 日日干狠狠操夜夜爽| 日本欧美国产在线视频| 精品人妻熟女av久视频| 夜夜夜夜夜久久久久| 一本精品99久久精品77| 亚洲三级黄色毛片| 欧美日韩国产亚洲二区| 人妻制服诱惑在线中文字幕| 免费av不卡在线播放| 22中文网久久字幕| 2021天堂中文幕一二区在线观| 我的女老师完整版在线观看| 免费高清视频大片| 美女xxoo啪啪120秒动态图| 毛片一级片免费看久久久久| 最后的刺客免费高清国语| 中文字幕免费在线视频6| 国产高清三级在线| 亚洲欧美成人精品一区二区| 成人高潮视频无遮挡免费网站| .国产精品久久| 午夜福利18| 在线a可以看的网站| 伦理电影大哥的女人| 日本成人三级电影网站| 国产高潮美女av| 国产精品av视频在线免费观看| 亚洲av.av天堂| av视频在线观看入口| 高清毛片免费观看视频网站| 热99re8久久精品国产| 亚洲欧美精品综合久久99| 久久人妻av系列| 又黄又爽又免费观看的视频| 欧美日韩综合久久久久久| 97碰自拍视频| 一区福利在线观看| 精品久久国产蜜桃| 亚洲aⅴ乱码一区二区在线播放| 日本一二三区视频观看| 高清毛片免费看| 欧美最黄视频在线播放免费| 婷婷精品国产亚洲av| 欧美一区二区精品小视频在线| 国产一区二区在线观看日韩| 黑人高潮一二区| 偷拍熟女少妇极品色| 精品一区二区三区视频在线| 欧美另类亚洲清纯唯美| 久久精品91蜜桃| 俄罗斯特黄特色一大片| 久久精品国产亚洲av香蕉五月| 有码 亚洲区| 国产麻豆成人av免费视频| 少妇的逼好多水| 亚洲精品粉嫩美女一区| 男人的好看免费观看在线视频| 三级毛片av免费| 国产综合懂色| 搡老熟女国产l中国老女人| av.在线天堂| 国产午夜精品久久久久久一区二区三区 | 亚洲精品在线观看二区| 天堂动漫精品| 一级毛片久久久久久久久女| 麻豆国产av国片精品| 老司机午夜福利在线观看视频| 99热这里只有是精品50| 日韩,欧美,国产一区二区三区 | 嫩草影院新地址| 女人被狂操c到高潮| 精品一区二区三区av网在线观看| 看十八女毛片水多多多| 成人毛片a级毛片在线播放| 国产一区二区激情短视频| 久99久视频精品免费| 中文字幕久久专区| 成人av一区二区三区在线看| av在线天堂中文字幕| 精品日产1卡2卡| 亚洲精品一卡2卡三卡4卡5卡| 高清日韩中文字幕在线| 禁无遮挡网站| 国产成人a区在线观看| 亚洲欧美精品自产自拍| 有码 亚洲区| 少妇人妻精品综合一区二区 | 好男人在线观看高清免费视频| 欧美极品一区二区三区四区| 欧美最新免费一区二区三区| 欧美bdsm另类| 午夜老司机福利剧场| 国产成人a区在线观看| 国产女主播在线喷水免费视频网站 | 色综合色国产| 久久精品久久久久久噜噜老黄 | 午夜精品一区二区三区免费看| 成人高潮视频无遮挡免费网站| 天天躁夜夜躁狠狠久久av| 少妇的逼好多水| 可以在线观看的亚洲视频| 成人综合一区亚洲| 99久久久亚洲精品蜜臀av| 亚洲内射少妇av| 国产精品美女特级片免费视频播放器| 精品国内亚洲2022精品成人| 免费在线观看影片大全网站| 亚洲内射少妇av| 欧美xxxx性猛交bbbb| 亚洲成av人片在线播放无| 国产精品无大码| АⅤ资源中文在线天堂| 别揉我奶头~嗯~啊~动态视频| 亚洲在线自拍视频| 国产高潮美女av| or卡值多少钱| 蜜桃亚洲精品一区二区三区| 免费无遮挡裸体视频| 精品乱码久久久久久99久播| videossex国产| 亚洲人成网站在线播| 免费av不卡在线播放| 婷婷色综合大香蕉| 国模一区二区三区四区视频| 成人综合一区亚洲| 亚洲欧美日韩无卡精品| 夜夜爽天天搞| 欧美一级a爱片免费观看看| 亚洲精品一区av在线观看| 亚洲18禁久久av| 成人无遮挡网站| 午夜福利成人在线免费观看| 99热精品在线国产| 波多野结衣高清无吗| 日韩成人av中文字幕在线观看 | 12—13女人毛片做爰片一| 国产精品爽爽va在线观看网站| 综合色丁香网| 麻豆av噜噜一区二区三区| 色综合色国产| 男女边吃奶边做爰视频| 成人二区视频| 欧美区成人在线视频| 亚洲国产欧美人成| 国产v大片淫在线免费观看| 精品不卡国产一区二区三区| 综合色av麻豆| 干丝袜人妻中文字幕| 午夜福利在线观看免费完整高清在 | 亚洲色图av天堂| 欧美性猛交黑人性爽| 精品一区二区三区视频在线| 国产色婷婷99| 色哟哟哟哟哟哟| 长腿黑丝高跟| 亚洲欧美中文字幕日韩二区| 99热精品在线国产| 两个人的视频大全免费| 麻豆国产97在线/欧美| 日韩 亚洲 欧美在线| 麻豆一二三区av精品| 村上凉子中文字幕在线| 少妇裸体淫交视频免费看高清| 在线播放国产精品三级| 婷婷精品国产亚洲av| 亚洲18禁久久av| 人妻少妇偷人精品九色| 国产亚洲精品久久久久久毛片| 露出奶头的视频| 别揉我奶头~嗯~啊~动态视频| av女优亚洲男人天堂| 色5月婷婷丁香| 最好的美女福利视频网| 精品乱码久久久久久99久播| 九色成人免费人妻av| 最近在线观看免费完整版| 日韩在线高清观看一区二区三区| av中文乱码字幕在线| 国产黄色小视频在线观看| 亚洲av不卡在线观看| 国语自产精品视频在线第100页| 久久婷婷人人爽人人干人人爱| 在线免费十八禁| av.在线天堂| 99riav亚洲国产免费| 内射极品少妇av片p| 亚洲综合色惰| av视频在线观看入口| 夜夜夜夜夜久久久久| 又黄又爽又免费观看的视频| 国产精品一及| 简卡轻食公司| 婷婷亚洲欧美| 国产精品国产三级国产av玫瑰| 日韩 亚洲 欧美在线| 日日摸夜夜添夜夜添av毛片| 欧美bdsm另类| 久久九九热精品免费| av专区在线播放| eeuss影院久久| 春色校园在线视频观看| 国产免费男女视频| 久久99热这里只有精品18| 久久九九热精品免费| 在线国产一区二区在线| av中文乱码字幕在线| 午夜福利在线观看吧| 精品久久久噜噜| 黄片wwwwww| 给我免费播放毛片高清在线观看| 午夜视频国产福利| 久久精品久久久久久噜噜老黄 | 国产又黄又爽又无遮挡在线| 日韩av不卡免费在线播放| 免费观看人在逋| 亚洲人成网站高清观看| 此物有八面人人有两片| 天天躁日日操中文字幕| 大又大粗又爽又黄少妇毛片口| 国产三级在线视频| 97热精品久久久久久| a级一级毛片免费在线观看| 99热网站在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美绝顶高潮抽搐喷水| av免费在线看不卡| 久久久久性生活片| 国产视频一区二区在线看| 女生性感内裤真人,穿戴方法视频| 欧美极品一区二区三区四区| 成年女人看的毛片在线观看| 成人漫画全彩无遮挡| 国产 一区 欧美 日韩| 男女啪啪激烈高潮av片| 一级毛片我不卡| 在线国产一区二区在线| 国产蜜桃级精品一区二区三区| 国产人妻一区二区三区在| 亚洲国产精品合色在线| 中国国产av一级| 国产又黄又爽又无遮挡在线| 欧美日本亚洲视频在线播放| 久久久久九九精品影院| av在线亚洲专区| 国产一区二区激情短视频| 欧美另类亚洲清纯唯美| 久久热精品热| 一边摸一边抽搐一进一小说| 亚洲中文日韩欧美视频| 久久中文看片网| 亚洲一区高清亚洲精品| 国内精品美女久久久久久| 插阴视频在线观看视频| 日韩成人伦理影院| 国产成年人精品一区二区| 午夜福利成人在线免费观看| 国产精品综合久久久久久久免费| 欧美xxxx黑人xx丫x性爽| 欧美日本视频| 国产精品一区二区三区四区久久| 亚洲不卡免费看| 最好的美女福利视频网| 亚洲最大成人手机在线| 久久久久国产精品人妻aⅴ院| 色哟哟·www| 亚洲欧美日韩高清在线视频| 日韩中字成人| 国产欧美日韩精品亚洲av| 欧美成人免费av一区二区三区| 老熟妇仑乱视频hdxx| 一本久久中文字幕| 不卡视频在线观看欧美| 久久精品综合一区二区三区| 亚洲在线自拍视频| 日本 av在线| 国产视频内射| 日本黄色视频三级网站网址| 69人妻影院| 搡老熟女国产l中国老女人| 久久精品影院6| 尤物成人国产欧美一区二区三区| 成人国产麻豆网| 美女高潮的动态| 22中文网久久字幕| 能在线免费观看的黄片| 九九在线视频观看精品| 日日摸夜夜添夜夜爱| 特级一级黄色大片| 97在线视频观看| 欧美一级a爱片免费观看看| 91麻豆精品激情在线观看国产| 极品教师在线视频| 自拍偷自拍亚洲精品老妇| 可以在线观看的亚洲视频| 免费观看在线日韩| 国产精品久久视频播放| 黄色配什么色好看| 免费看美女性在线毛片视频| a级一级毛片免费在线观看| 国产精品人妻久久久影院| 亚洲第一电影网av| 菩萨蛮人人尽说江南好唐韦庄 | 免费观看精品视频网站| 国产一区二区在线av高清观看| 赤兔流量卡办理| 中国国产av一级| 久久鲁丝午夜福利片| 婷婷色综合大香蕉| 国内精品宾馆在线| 22中文网久久字幕| 亚洲欧美中文字幕日韩二区| 亚洲激情五月婷婷啪啪| a级毛色黄片| 亚洲av五月六月丁香网| 少妇被粗大猛烈的视频| 一级av片app| 五月伊人婷婷丁香| 桃色一区二区三区在线观看| 精品人妻视频免费看| 91精品国产九色| 一a级毛片在线观看| 亚洲国产精品久久男人天堂| 久久久久久久久大av| 国产aⅴ精品一区二区三区波| 国产私拍福利视频在线观看| 日本一本二区三区精品| 亚洲内射少妇av| 国产伦在线观看视频一区| 日韩强制内射视频| 一个人看的www免费观看视频| 性色avwww在线观看| 老熟妇仑乱视频hdxx| 午夜激情福利司机影院| 国产一区二区三区av在线 | 高清毛片免费观看视频网站| 在线天堂最新版资源| 国产人妻一区二区三区在| 最好的美女福利视频网|