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

    過冷度影響海水結(jié)冰形狀與速度的相場模擬

    2021-06-01 06:21:16楊蘇杰
    上海交通大學(xué)學(xué)報 2021年5期
    關(guān)鍵詞:相場晶核枝晶

    白 旭, 楊蘇杰

    (江蘇科技大學(xué),船舶與海洋工程學(xué)院,江蘇 鎮(zhèn)江,212003)

    近年來,隨著全球氣候變暖,北極地區(qū)冰川消融,北極航道得到開發(fā).但是,由于北極地區(qū)特殊的地理環(huán)境,船舶在北極航行時極易受到結(jié)冰的影響[1].已有的針對船舶結(jié)冰的研究大多集中在宏觀層面,即不同氣候條件下,海洋結(jié)構(gòu)物各部位產(chǎn)生的結(jié)冰以及結(jié)冰造成的影響,結(jié)冰的原理也主要來自于經(jīng)驗公式和觀測結(jié)果總結(jié),很少有對其產(chǎn)生原因和微觀機理的分析[2-4].通過研究海水結(jié)冰的微觀機理,用實際環(huán)境參數(shù)來模擬海水結(jié)冰過程,可以得到更加真實的冰晶生長情況,為過冷條件下海水飛沫結(jié)冰提供微觀理論基礎(chǔ).針對海水的微觀結(jié)冰過程,現(xiàn)在國內(nèi)相關(guān)研究主要集中在單個晶核的凝固現(xiàn)象,而在實際情況下,外界環(huán)境所導(dǎo)致的結(jié)冰幾乎都是多晶核條件下的結(jié)冰.因此,本研究不僅會分析過冷度對海水結(jié)冰的影響,同時會涉及多晶核條件下的海水凝固情況.研究水滴凝固的微觀機理的重點在于追蹤冰晶凝固時固液界面的移動,而通過引入相場法,建立溫度、相場等多個微分方程,可以將真實環(huán)境條件引入模擬計算,模擬固液擴散界面,從而得到冰晶的生成演化過程[5-6].

    相場法最早是在20世紀(jì)80年代由Collins等提出[7].1993年,Kobayashi首先建立了關(guān)于相場和溫度的相場模型,在考慮各向異性的情況下,模擬了金屬凝固時的枝晶生長[8].Wheeler等[9]模擬了Ni的凝固過程,并提出了二元合金的相場方程,即WBM模型,并將相場模擬與實際物理量聯(lián)系起來,定量分析參數(shù)對凝固現(xiàn)象的影響.Fan 等[10]基于平均濃度場和熱、質(zhì)及動量輸運現(xiàn)象,模擬了蛋白質(zhì)溶液在冷凍過程中蛋白質(zhì)的分離與重分布.Yang等[11]使用冰模板,考慮石墨烯再分配和冰晶的各向異性,模擬了在水凍結(jié)時石墨烯氣凝膠的形成過程.

    國內(nèi)學(xué)者對于相場模型的研究也有一定進展.楊燕等[12]利用Wheeler模型模擬了過冷條件下不同過冷度對冰晶生長的影響.朱昌盛等[13]基于KKS模型模擬了在等溫和不等溫兩種條件下潛熱對二元合金凝固的影響.于艷梅等[14]對單相和二元合金分別進行相場模擬,分析空間尺度、外界噪聲及各向異性等參數(shù)對過冷熔體凝固過程的影響.韓端鋒等[15-16]分別使用Kobayashi模型和Wheeler模型,模擬了單晶核條件下純水和海水的凝固過程,定量分析了外界參數(shù)對冰晶凝固的影響.

    針對船舶結(jié)構(gòu)的宏觀結(jié)冰與微觀機理的研究方面,Ryerson等[17]利用船舶切割下的冰晶進行研究,將結(jié)冰的宏觀現(xiàn)象和微觀機理相聯(lián)系,得到了結(jié)構(gòu)結(jié)冰的物理性質(zhì)和結(jié)冰部位的微觀分析.Singh等[18]利用分子動力學(xué)定量研究了光滑粗糙表面上過冷水滴的行為,分析了在粗糙表面上凝固溫度和液滴黏附功與粗糙度之間的聯(lián)系.Sultana等[19]利用數(shù)值模擬分析了單個鹽水滴在撞擊極冷表面后的相變情況,探討了液相率、鹽度等物理性質(zhì)對液滴凍結(jié)的影響.

    本文基于Wheeler相場模型,運用有限差分法,對不同過冷度下海水結(jié)冰情況進行模擬,研究過冷度對海水凝固時冰晶生長形狀和速度的影響,為船舶結(jié)冰預(yù)報提供理論和技術(shù)支持.

    1 Wheeler相場模型理論

    1.1 相場方程

    相場法以金茲堡-道朗方程為研究基礎(chǔ)[20],通過引入相場參數(shù)φ來表示固液界面的移動.當(dāng)φ=1時表示固相,φ=0時表示液相,0<φ<1時表示固液界面,通過數(shù)值的變化來實現(xiàn)對固液界面的跟蹤.

    金茲堡-道朗自由能方程為

    (1)

    式中:f(φ,T)和β(T)分別為自由能密度函數(shù)和熱力學(xué)驅(qū)動力函數(shù),β(T)<1/2;W為單位體積能量;p(φ)為固相分數(shù),1-p(φ)為液相分數(shù),p(φ)=φ3(10-15φ+6φ2).

    相場方程為

    (2)

    其中:δ為特征長度;ω為長度參考尺寸;μ為界面遷移率;σ為表面能;TM為海水凝固點;K為熱擴散系數(shù);L為單位體積潛熱;cp為比熱容;ΔT為溶液過冷度,ΔT=TM-T0,T0為溶液初始溫度.

    溫度場方程為

    (3)

    式中:p(φ)的導(dǎo)數(shù)p′(φ)=30φ2(1-φ)2;無量綱溫度u=(T-TM)/(TM-T0).

    Wheeler相場將實際物理量與相場參數(shù)相結(jié)合,因此在計算中,應(yīng)該使用實際物理參數(shù)進行計算.但是為了將界面厚度控制在可計算范圍內(nèi),取ω為2.7×10-4cm,將其無量綱化,實際界面厚度等于無量綱界面厚度與長度參考尺寸的乘積,網(wǎng)格間距也按照同樣方法進行換算.

    1.2 溶質(zhì)場方程

    在模擬海水凝固時,將海水視為水和鹽的二元溶液,相比于一元溶液,二元溶液需要考慮溶質(zhì)的影響.溶質(zhì)對凝固的影響主要有兩種方式,一種是引入溶質(zhì)對過冷度的影響,另一種是在方程中添加濃度項[21].本文通過引入過冷度變化,來實現(xiàn)溶質(zhì)擴散對溶液凝固的影響.溶質(zhì)場方程為

    (4)

    (5)

    式中:C為海水中鹽的摩爾分數(shù),海水中鹽的質(zhì)量分數(shù)約為3.5%,則鹽的摩爾分數(shù)為0.011;k0為平衡分配系數(shù);D為無量綱擴散系數(shù);Ds為固相擴散率;Dl為液相擴散率.

    模擬過程中,計算域內(nèi)的過冷度受到溶質(zhì)濃度影響,實時過冷度為

    (6)

    式中:mL為液相線斜率;C0為溶液初始質(zhì)量分數(shù).

    1.3 參數(shù)取值

    本文所取物性參數(shù)均按照質(zhì)量分數(shù)為3.5%的海水溶液選取[15],具體取值如表1所示.

    表1 海水溶液物性參數(shù)Tab.1 Physical properties of seawater solution

    2 數(shù)值求解與可靠性驗證

    在得到相場模型后,需要對其進行數(shù)值求解.有限差分法原理簡單,易于編程,因此本文采用有限差分法來求解控制方程.有限差分法是將要求的函數(shù)在時間和空間上離散,用差分代替微分,并用離散的值來表示函數(shù)在計算域內(nèi)的變化,結(jié)合初始條件和邊界條件,求出函數(shù)在計算域內(nèi)各網(wǎng)格節(jié)點的值.

    2.1 有限差分法及離散方式

    求解相場控制方程時,空間離散采用中心差分方法,時間離散采用顯式差分方法.另外,對于方程中出現(xiàn)的拉普拉斯算子,采用九點差分方法進行離散:

    φi+1, j+1+φi-1, j-1+φi-1, j+1+φi+1, j-1]/

    [3(Δx)2]

    式中:Δx、Δy為網(wǎng)格間距;Δt為時間步長.對于溫度場和溶質(zhì)場方程,采用與相場方程同樣的差分格式.

    在計算域內(nèi),為了滿足計算精度,保證計算結(jié)果收斂,需要同時考慮相場、溫度場及溶質(zhì)場.計算域內(nèi)的時間步長和空間步長需滿足以下條件

    (7)

    (8)

    2.2 初始條件和邊界條件

    本文計算域設(shè)置為正方形,網(wǎng)格數(shù)設(shè)置為1 500×1 500,網(wǎng)格間距Δx和Δy在無量綱化后均為0.03,時間步長Δt=0.000 1 s.在計算域內(nèi),按照正方形設(shè)置4個初始晶核,晶核半徑為R.初始晶核距離計算域邊界均為50R,晶核之間距離為50R.晶核分布如圖1所示.

    圖1 初始晶核大小與位置Fig.1 Initial nuclei size and position

    晶核外相場參數(shù)φ設(shè)為1,無量綱溫度u為-Δ,質(zhì)量分數(shù)C為C0,具體晶核位置和條件設(shè)置如下式:

    (x-500)2+(y-500)2≤R2時,

    φ=0,u=-Δ,C=C0

    (9)

    (x-1 000)2+(y-500)2≤R2時,

    φ=0,u=-Δ,C=C0

    (10)

    (x-500)2+(y-1000)2≤R2時,

    φ=0,u=-Δ,C=C0

    (11)

    (x-1 000)2+(y-1 000)2≤R2時,

    φ=0,u=-Δ,C=C0

    (12)

    相場、溫度場和濃度場均采用Neumann邊界條件,即

    (13)

    2.3 模擬冰晶與自然冰晶的對比

    為了得到海水結(jié)冰的過程及規(guī)律,首先對冰晶凝固過程進行模擬,并將其和實際自然條件下生成的冰晶進行對比.使用上文設(shè)置參數(shù),網(wǎng)格劃分為500×500, 模擬時間為7 500Δt,得到結(jié)果如圖2所示.冰晶在凝固時,先由晶核發(fā)展為冰盤,隨后主要沿著6個方向發(fā)展,形成主枝,主枝之間的夾角約為60°,主枝隨后向兩側(cè)長出側(cè)枝和多次分枝,最終形成六角冰晶.通過比較模擬的冰晶形狀和現(xiàn)實的冰晶形狀[22],發(fā)現(xiàn)實驗結(jié)果吻合良好,可以進行進一步模擬計算.

    圖2 模擬冰晶與實際冰晶形狀對比Fig.2 Comparison of simulated ice crystal and actual ice crystal

    3 過冷度對晶核生長形狀及速度的影響

    3.1 過冷度對晶核生長形狀的影響

    根據(jù)凝固理論,過冷度是液體凝固的驅(qū)動力,在冰晶的形核和生長中起重要作用.本文通過分析冰晶形貌、生長速度以及枝晶間距來描述過冷度對冰晶生長的影響.

    選擇0.80、0.85、0.90、0.95及1.00這5個不同無量綱過冷度作為模擬溫度,各無量綱過冷度對應(yīng)實際過冷度分別為 -63.34、-67.30、-71.26、-75.22 及 -79.18 K.總模擬時間為 7 500Δt,模擬結(jié)果如圖3所示.

    從圖中可以看出,當(dāng)Δ=0.80時,枝晶生長緩慢,僅有少量二次分枝,枝晶直到模擬結(jié)束未發(fā)生交匯,仍然是晶核各自發(fā)展.隨著無量綱過冷度的增大,枝晶生長速度加快,主枝變粗,分枝數(shù)量增加,但枝晶之間仍未發(fā)生交匯.當(dāng)Δ=0.90時,生長速度進一步增加,出現(xiàn)二級以上多次分枝,枝晶在模擬過程中出現(xiàn)交匯.交匯后主枝停止生長,分枝未受影響,繼續(xù)向外生長.Δ=1.00時,粗大的主枝占據(jù)了極大的空間,使得分枝數(shù)量減少.當(dāng)主枝交匯后,主枝逐漸變粗,分枝繼續(xù)發(fā)展,逐漸填滿晶核間空隙,幾乎形成冰盤.

    圖4為各過冷度下分別模擬 10 000Δt得到結(jié)果.通過比較各過冷度凝固結(jié)果可以看出,在同時存在多個晶核時,過冷度越大,枝晶就越早接觸,枝晶內(nèi)部空隙就越小,以至于最終完全接觸,內(nèi)部不再存在空隙.過冷度較小時,橫向與縱向主枝有粗細差別,橫向主枝較粗,而過冷度較大時,橫向與縱向主枝粗細基本相同.

    圖3 冰晶凝固過程Fig.3 Solidification process of ice crystals

    圖4 10 000Δt各無量綱過冷度凝固結(jié)果Fig.4 Solidification results of each dimensionless undercooling at 10 000Δt

    3.2 過冷度對晶核生長速度的影響

    為了定量分析過冷度對海水凝固的影響,使用枝晶尖端生長速度vtip來描述冰晶生長情況,枝晶尖端生長速度與無量綱時間和無量綱過冷度之間關(guān)系如圖5所示.

    圖5(a)為無量綱過冷度為0.8時,尖端速度隨無量綱時間變化情況.根據(jù)計算結(jié)果,枝晶在生長過程中較為穩(wěn)定,尖端速度在開始時較大,隨后慢慢減小,近似呈正弦函數(shù)波動.枝晶尖端速度隨無量綱過冷度的變化如圖5(b)所示,在無量綱過冷度從0.8增至1.0的過程中,枝晶尖端平均速度分別為0.288、0.360、0.432、0.504、0.576 cm/s,枝晶生長速度隨著無量綱過冷度的增加而明顯增加.無量綱過冷度對冰晶凝結(jié)的主要影響在于,無量綱過冷度越大,界面厚度就越薄,凝固時散發(fā)的潛熱就越多,從而使冰晶凝固速度更快.因此,代入實際物理量后可知,過冷度越低,冰晶凝固越快,且凝固速度與過冷度基本為線性關(guān)系.

    圖5 枝晶尖端速度變化情況Fig.5 Changes of dendrite tip velocity

    圖6 枝晶尖端距離及枝晶距離斜率變化Fig.6 Changes of dendrite tip distance and dendrite distance slope

    3.3 過冷度對多晶核生長形狀的影響

    本文在計算域內(nèi)設(shè)置4個晶核,在晶核初始間距相同的情況下,通過計算不同過冷度下各晶核主枝之間的距離變化來分析過冷度對海水結(jié)冰的影響情況.

    不同過冷度下枝晶尖端距離Dtip隨無量綱時間的變化如圖6(a)所示.隨著無量綱時間的增加,枝晶尖端間距逐漸減小,且二者呈線性關(guān)系.當(dāng)無量綱過冷度達到0.90時,枝晶間距最終為0,代表不同晶核生長出的枝晶最終接觸.根據(jù)圖中趨勢,無量綱過冷度越大,接觸所需時間越短.圖6(b)為將枝晶尖端距離隨時間變化的曲線擬合為一次函數(shù)后的曲線斜率Ktip,其代表枝晶間距的變化快慢,即枝晶生長速度的大小.根據(jù)曲線所示,斜率與無量綱過冷度之間近似呈線性關(guān)系,說明了隨著無量綱過冷度的增大,枝晶生長速度的變化也會相應(yīng)地近似呈線性增大.

    4 結(jié)論

    基于Wheeler相場模型,使用有限差分法,將海水看作鹽-水二元溶液,把實際物理量轉(zhuǎn)換為方程中的無量綱數(shù),提取出枝晶生長速度作為定量分析海水結(jié)冰的指標(biāo),并通過不同過冷度下各晶核枝晶之間的間距變化來描述過冷度對海水結(jié)冰的影響.具體結(jié)論如下:

    (1) 無量綱過冷度為0.90以下時,枝晶較細,沒有或僅有二級分枝,枝晶生長較慢.無量綱過冷度為0.90和0.95時,枝晶較粗,枝晶生長速度較快,出現(xiàn)二級以上分枝.無量綱過冷度為1.00時,主枝與分枝明顯變粗,生長速度明顯加快,分枝之間幾乎無空隙.

    (2) 無量綱過冷度為0.80時,枝晶生長慢,形狀細,晶核內(nèi)部空隙大,最終形成各主枝組合成的枝狀冰晶.無量綱過冷度為0.90和0.95時,有多次分枝,分枝之間有一定空隙,且空隙隨著分枝發(fā)展而減少.無量綱過冷度1.00時,枝晶生長速度快,分枝之間空隙小,晶核內(nèi)部空隙被填滿,形成冰盤.

    (3) 在各過冷度下,枝晶生長速度較為穩(wěn)定,尖端最大速度隨無量綱時間呈正弦變化.無量綱過冷度從0.80增至1.00時,枝晶尖端生長速度從0.12 cm/s增至0.24 cm/s,生長速度與過冷度之間為線性關(guān)系.

    (4) 枝晶尖端間距隨著時間的增加而減小,且二者呈線性關(guān)系,間距曲線斜率與無量綱過冷度之間近似為線性關(guān)系,其反映了枝晶生長速度的與無量綱過冷度之間也近似為線性關(guān)系,與之前求得的枝晶尖端速度相吻合.

    猜你喜歡
    相場晶核枝晶
    蒸汽在含有不可溶核和可溶無機鹽的細顆粒物表面的核化特性
    基于子單元光滑有限元的混凝土相場損傷模型研究
    晶核劑對燒結(jié)法制備Li2O-A12O3-SiO2系微晶玻璃晶化過程的影響
    鑄件凝固微觀組織仿真程序開發(fā)
    基于相場理論的瀝青自愈合微觀進程與機理研究進展
    石油瀝青(2018年1期)2018-04-12 07:31:51
    不同晶核劑對聚甲醛MC90性能的影響研究
    基于COMSOL的相場模擬研究
    科技視界(2017年8期)2017-07-31 10:31:17
    A356合金近液相線半固態(tài)鑄造非枝晶組織與模鍛成形性能
    不同形狀橫向限制對枝晶間距影響作用的相場法模擬
    合金定向凝固一次枝晶間距模擬
    国产精品永久免费网站| 精品一区二区三区人妻视频| 一级a爱片免费观看的视频| 变态另类丝袜制服| 久久国内精品自在自线图片| 一本精品99久久精品77| 99视频精品全部免费 在线| 男女下面进入的视频免费午夜| 色哟哟·www| 深爱激情五月婷婷| 最近视频中文字幕2019在线8| 99久久精品热视频| 亚洲国产精品sss在线观看| 国产一区二区在线观看日韩| 欧美bdsm另类| 国产精品一及| 国产国拍精品亚洲av在线观看| 久久久精品94久久精品| 偷拍熟女少妇极品色| 男女之事视频高清在线观看| 国产成人a∨麻豆精品| 高清毛片免费观看视频网站| 人妻少妇偷人精品九色| 成年免费大片在线观看| 可以在线观看毛片的网站| 一个人看视频在线观看www免费| 亚洲中文字幕一区二区三区有码在线看| 久久精品国产清高在天天线| 美女cb高潮喷水在线观看| 国产精品永久免费网站| 亚洲人成网站在线播放欧美日韩| 亚洲自偷自拍三级| 国产淫片久久久久久久久| 亚洲在线观看片| 国产午夜福利久久久久久| 久久久久久久久大av| 搡老妇女老女人老熟妇| 蜜桃亚洲精品一区二区三区| 伦理电影大哥的女人| 成人三级黄色视频| 免费人成视频x8x8入口观看| 亚洲人成网站在线播| 少妇熟女aⅴ在线视频| 自拍偷自拍亚洲精品老妇| 综合色av麻豆| 女人被狂操c到高潮| 久久久久免费精品人妻一区二区| 日产精品乱码卡一卡2卡三| 精品一区二区免费观看| 嫩草影院精品99| 露出奶头的视频| 嫩草影院新地址| 俺也久久电影网| 国产成人aa在线观看| 女的被弄到高潮叫床怎么办| 国产aⅴ精品一区二区三区波| 亚洲三级黄色毛片| 国产精品人妻久久久久久| 国产毛片a区久久久久| 伦理电影大哥的女人| 亚洲欧美精品自产自拍| 人人妻,人人澡人人爽秒播| 少妇丰满av| 国产亚洲91精品色在线| 中国国产av一级| 成人午夜高清在线视频| 国产精华一区二区三区| 51国产日韩欧美| 一边摸一边抽搐一进一小说| 成人性生交大片免费视频hd| 欧美中文日本在线观看视频| 国产 一区 欧美 日韩| 亚洲精华国产精华液的使用体验 | 国产伦一二天堂av在线观看| 亚洲国产精品合色在线| 国产精品无大码| 日韩欧美免费精品| 亚洲国产精品久久男人天堂| 亚洲欧美精品自产自拍| a级毛片免费高清观看在线播放| 亚洲va在线va天堂va国产| 亚洲人成网站在线播放欧美日韩| 亚洲性久久影院| av在线亚洲专区| 插逼视频在线观看| 久久欧美精品欧美久久欧美| 在线观看一区二区三区| 国产国拍精品亚洲av在线观看| 国产精品女同一区二区软件| 热99在线观看视频| 成人特级黄色片久久久久久久| 久久久久久九九精品二区国产| 最好的美女福利视频网| 日本 av在线| 久久久久久久午夜电影| 亚洲av熟女| 亚洲成人久久爱视频| 久久国内精品自在自线图片| 国产精品嫩草影院av在线观看| 老司机福利观看| 国产乱人偷精品视频| 中文字幕久久专区| 国产成人freesex在线 | 99久久精品热视频| 中国国产av一级| 国产亚洲av嫩草精品影院| 内地一区二区视频在线| 国产伦在线观看视频一区| 久久久久久九九精品二区国产| 全区人妻精品视频| 人妻制服诱惑在线中文字幕| 高清午夜精品一区二区三区 | 欧美绝顶高潮抽搐喷水| 日韩精品中文字幕看吧| 搞女人的毛片| 伊人久久精品亚洲午夜| 不卡视频在线观看欧美| 日韩,欧美,国产一区二区三区 | 观看美女的网站| 91久久精品国产一区二区三区| 久久久久久久久久黄片| av免费在线看不卡| www.色视频.com| 亚洲一级一片aⅴ在线观看| 国产高清视频在线观看网站| 国产成年人精品一区二区| 免费不卡的大黄色大毛片视频在线观看 | 99热这里只有是精品50| 丝袜喷水一区| 国产成人a区在线观看| 亚洲人成网站高清观看| 国产av麻豆久久久久久久| 亚洲av中文字字幕乱码综合| 日韩成人伦理影院| 久久精品国产亚洲av涩爱 | 级片在线观看| 内地一区二区视频在线| 观看免费一级毛片| 成人一区二区视频在线观看| 小蜜桃在线观看免费完整版高清| 国产精品人妻久久久久久| 亚洲无线观看免费| 搡老熟女国产l中国老女人| 亚洲人成网站在线观看播放| 天堂av国产一区二区熟女人妻| 久久久精品大字幕| 午夜精品一区二区三区免费看| 国产成人a区在线观看| 国产探花在线观看一区二区| 午夜日韩欧美国产| 午夜免费男女啪啪视频观看 | 精华霜和精华液先用哪个| 男女做爰动态图高潮gif福利片| 欧美激情久久久久久爽电影| 少妇高潮的动态图| 插逼视频在线观看| 亚洲激情五月婷婷啪啪| 无遮挡黄片免费观看| 一区福利在线观看| 在线免费观看的www视频| 欧美一区二区国产精品久久精品| 搡老妇女老女人老熟妇| 午夜老司机福利剧场| 最近2019中文字幕mv第一页| 欧美丝袜亚洲另类| 人妻丰满熟妇av一区二区三区| 尾随美女入室| 国产精品亚洲一级av第二区| 亚洲真实伦在线观看| 久久中文看片网| 99九九线精品视频在线观看视频| 99久久精品国产国产毛片| 少妇人妻精品综合一区二区 | 我要看日韩黄色一级片| 亚洲五月天丁香| 秋霞在线观看毛片| 国内久久婷婷六月综合欲色啪| 免费观看精品视频网站| 日日摸夜夜添夜夜添小说| 天天一区二区日本电影三级| 欧美日本亚洲视频在线播放| 国国产精品蜜臀av免费| 国产淫片久久久久久久久| 天堂av国产一区二区熟女人妻| 久久久久久九九精品二区国产| 国产一区二区在线av高清观看| 最近2019中文字幕mv第一页| 啦啦啦韩国在线观看视频| 永久网站在线| 欧美精品国产亚洲| 性插视频无遮挡在线免费观看| 亚洲精品456在线播放app| 国产伦一二天堂av在线观看| 蜜臀久久99精品久久宅男| 五月玫瑰六月丁香| 99riav亚洲国产免费| 成人亚洲精品av一区二区| 97超碰精品成人国产| 久久热精品热| 高清日韩中文字幕在线| 夜夜爽天天搞| 成年av动漫网址| 精品不卡国产一区二区三区| 日本-黄色视频高清免费观看| 亚洲精品一区av在线观看| 亚洲国产欧美人成| 麻豆乱淫一区二区| 亚洲色图av天堂| 国内揄拍国产精品人妻在线| 亚洲婷婷狠狠爱综合网| 老女人水多毛片| 午夜精品一区二区三区免费看| 插逼视频在线观看| 亚洲精品国产成人久久av| 女人被狂操c到高潮| 黄色配什么色好看| 无遮挡黄片免费观看| 女人十人毛片免费观看3o分钟| 最新在线观看一区二区三区| 久久久久久久久久成人| 亚洲美女视频黄频| 俺也久久电影网| 亚洲天堂国产精品一区在线| 亚洲成人av在线免费| 在线免费观看的www视频| 午夜福利在线观看吧| 欧美极品一区二区三区四区| 国产精品亚洲美女久久久| 亚洲精品在线观看二区| 欧美一区二区国产精品久久精品| 最近手机中文字幕大全| 亚洲最大成人手机在线| a级毛色黄片| 国产在线精品亚洲第一网站| 亚洲中文字幕日韩| 免费电影在线观看免费观看| 中国美女看黄片| 高清毛片免费看| 亚洲性久久影院| 久久精品91蜜桃| 国产 一区精品| av视频在线观看入口| 91精品国产九色| 99久久九九国产精品国产免费| 国内揄拍国产精品人妻在线| 男女视频在线观看网站免费| 亚洲精品乱码久久久v下载方式| 人人妻人人澡人人爽人人夜夜 | 欧美色视频一区免费| 国产精品不卡视频一区二区| 少妇的逼好多水| 亚洲美女搞黄在线观看 | 久久久久久久久久久丰满| 亚洲精品亚洲一区二区| 毛片女人毛片| 国产精品亚洲一级av第二区| 校园人妻丝袜中文字幕| 中文资源天堂在线| 亚洲最大成人av| 久久久a久久爽久久v久久| 五月伊人婷婷丁香| 秋霞在线观看毛片| 男女啪啪激烈高潮av片| 成人av在线播放网站| 国产精品1区2区在线观看.| 免费大片18禁| 中文资源天堂在线| 少妇的逼好多水| 99久久成人亚洲精品观看| 国产精品无大码| 可以在线观看的亚洲视频| 一级黄色大片毛片| 嫩草影院入口| 国产白丝娇喘喷水9色精品| 俄罗斯特黄特色一大片| 小蜜桃在线观看免费完整版高清| 麻豆国产97在线/欧美| 亚洲精品乱码久久久v下载方式| 国产三级在线视频| 少妇人妻精品综合一区二区 | 人人妻,人人澡人人爽秒播| videossex国产| 国产精品电影一区二区三区| 99热只有精品国产| 久久久a久久爽久久v久久| 男人舔女人下体高潮全视频| 极品教师在线视频| 91av网一区二区| 非洲黑人性xxxx精品又粗又长| 永久网站在线| 高清毛片免费观看视频网站| 在线观看66精品国产| 日本免费a在线| 日韩一区二区视频免费看| 一进一出好大好爽视频| 在线播放国产精品三级| 亚洲欧美成人综合另类久久久 | 精品人妻熟女av久视频| 精品熟女少妇av免费看| 插阴视频在线观看视频| 一级黄色大片毛片| 97热精品久久久久久| 日韩精品中文字幕看吧| 欧美zozozo另类| 欧美人与善性xxx| 女生性感内裤真人,穿戴方法视频| 国产aⅴ精品一区二区三区波| 97人妻精品一区二区三区麻豆| 午夜激情欧美在线| 国产乱人视频| 乱系列少妇在线播放| 熟女电影av网| 男人的好看免费观看在线视频| 日韩精品中文字幕看吧| 日日摸夜夜添夜夜爱| 国产免费男女视频| 亚洲av一区综合| 国产精品人妻久久久久久| 精品久久久久久久久av| 啦啦啦韩国在线观看视频| 日韩三级伦理在线观看| 淫妇啪啪啪对白视频| 国产精品一及| 国产精品一区二区免费欧美| 亚洲精品色激情综合| 级片在线观看| 精品乱码久久久久久99久播| 亚洲人成网站高清观看| 亚洲国产欧洲综合997久久,| 免费观看精品视频网站| 麻豆国产97在线/欧美| 天堂网av新在线| 国产精品国产三级国产av玫瑰| 亚洲第一电影网av| 日韩制服骚丝袜av| 日韩欧美国产在线观看| 日韩欧美三级三区| 精品久久久噜噜| av免费在线看不卡| 99热这里只有是精品50| 亚洲熟妇中文字幕五十中出| 亚洲在线观看片| 亚洲av美国av| 亚洲在线自拍视频| 国产黄a三级三级三级人| 啦啦啦观看免费观看视频高清| 国产私拍福利视频在线观看| 国产大屁股一区二区在线视频| 丰满的人妻完整版| 十八禁国产超污无遮挡网站| 亚洲自偷自拍三级| 99热精品在线国产| 一级毛片电影观看 | 日本撒尿小便嘘嘘汇集6| 国产爱豆传媒在线观看| 日韩亚洲欧美综合| 精品人妻偷拍中文字幕| 中文字幕精品亚洲无线码一区| 国产乱人视频| 国产精品一二三区在线看| av中文乱码字幕在线| 狂野欧美白嫩少妇大欣赏| 久久热精品热| 99热这里只有是精品在线观看| 日本五十路高清| 国产日本99.免费观看| 99国产极品粉嫩在线观看| 免费在线观看成人毛片| 最好的美女福利视频网| 91精品国产九色| av国产免费在线观看| 国产单亲对白刺激| 婷婷色综合大香蕉| 精品99又大又爽又粗少妇毛片| 久久精品久久久久久噜噜老黄 | 一进一出抽搐动态| 此物有八面人人有两片| 成人鲁丝片一二三区免费| 深夜a级毛片| aaaaa片日本免费| .国产精品久久| 在线观看午夜福利视频| 免费观看在线日韩| 国产欧美日韩一区二区精品| 亚洲激情五月婷婷啪啪| 久久久国产成人精品二区| 色在线成人网| 成年免费大片在线观看| eeuss影院久久| 性欧美人与动物交配| 日本色播在线视频| 色尼玛亚洲综合影院| 国产私拍福利视频在线观看| 欧美最黄视频在线播放免费| a级毛片免费高清观看在线播放| 18+在线观看网站| 亚洲最大成人av| av在线亚洲专区| 搡老妇女老女人老熟妇| 国产精品1区2区在线观看.| 久久久久国内视频| 免费在线观看影片大全网站| 国产黄色视频一区二区在线观看 | 国产黄片美女视频| 午夜老司机福利剧场| 免费观看在线日韩| 狂野欧美白嫩少妇大欣赏| 日产精品乱码卡一卡2卡三| 女人被狂操c到高潮| 看非洲黑人一级黄片| 亚洲图色成人| av免费在线看不卡| 国产av麻豆久久久久久久| 国产精品爽爽va在线观看网站| 日韩在线高清观看一区二区三区| 两个人视频免费观看高清| 亚洲精品一卡2卡三卡4卡5卡| 久久久久九九精品影院| 国内少妇人妻偷人精品xxx网站| 日韩 亚洲 欧美在线| 久久欧美精品欧美久久欧美| 如何舔出高潮| 日韩精品有码人妻一区| 老师上课跳d突然被开到最大视频| 亚洲国产欧洲综合997久久,| 97超碰精品成人国产| 天美传媒精品一区二区| h日本视频在线播放| 在线观看午夜福利视频| 国产高清视频在线播放一区| 麻豆av噜噜一区二区三区| 少妇丰满av| 一本精品99久久精品77| 成年女人永久免费观看视频| 日日摸夜夜添夜夜添av毛片| 一个人观看的视频www高清免费观看| 亚洲精品一卡2卡三卡4卡5卡| 禁无遮挡网站| 中文字幕av成人在线电影| 麻豆av噜噜一区二区三区| 色噜噜av男人的天堂激情| 亚洲天堂国产精品一区在线| 成年女人永久免费观看视频| 91狼人影院| 成人午夜高清在线视频| av.在线天堂| 国产黄色视频一区二区在线观看 | 亚洲激情五月婷婷啪啪| 精品午夜福利视频在线观看一区| 色综合亚洲欧美另类图片| 中文字幕久久专区| 中国国产av一级| 中文字幕熟女人妻在线| 真实男女啪啪啪动态图| 一级av片app| 亚洲av中文av极速乱| 久久九九热精品免费| 亚洲成人av在线免费| av卡一久久| 国内精品一区二区在线观看| 免费观看人在逋| 尾随美女入室| 舔av片在线| 日韩一区二区视频免费看| 国产成年人精品一区二区| 天堂网av新在线| 我要看日韩黄色一级片| 人人妻人人澡人人爽人人夜夜 | 国产亚洲精品综合一区在线观看| av在线播放精品| 亚洲中文字幕日韩| 韩国av在线不卡| 婷婷色综合大香蕉| 亚洲在线观看片| av在线天堂中文字幕| 国语自产精品视频在线第100页| 毛片女人毛片| 国产午夜精品久久久久久一区二区三区 | 91狼人影院| 欧美xxxx黑人xx丫x性爽| 我的老师免费观看完整版| 久久精品国产亚洲网站| 搡女人真爽免费视频火全软件 | 久久这里只有精品中国| 美女被艹到高潮喷水动态| 久久综合国产亚洲精品| 亚洲人成网站在线播放欧美日韩| 一级黄片播放器| 看非洲黑人一级黄片| 99热这里只有是精品50| 免费看美女性在线毛片视频| 一级黄色大片毛片| 秋霞在线观看毛片| 日韩欧美一区二区三区在线观看| 亚洲中文日韩欧美视频| 六月丁香七月| 国产爱豆传媒在线观看| 成人亚洲精品av一区二区| 波多野结衣高清作品| 亚洲性久久影院| 日韩高清综合在线| 国产男靠女视频免费网站| 日本三级黄在线观看| 国产av不卡久久| 伊人久久精品亚洲午夜| 在现免费观看毛片| 美女免费视频网站| 亚洲欧美日韩高清在线视频| 日本爱情动作片www.在线观看 | 最近的中文字幕免费完整| 精华霜和精华液先用哪个| 亚洲人成网站高清观看| 日韩成人av中文字幕在线观看 | 麻豆av噜噜一区二区三区| 香蕉av资源在线| 久久久欧美国产精品| 亚洲五月天丁香| 日韩一本色道免费dvd| 免费黄网站久久成人精品| 久久人人爽人人爽人人片va| 国内精品一区二区在线观看| 亚洲国产精品成人久久小说 | 性插视频无遮挡在线免费观看| 国产爱豆传媒在线观看| 成熟少妇高潮喷水视频| 变态另类丝袜制服| 国产白丝娇喘喷水9色精品| 亚洲成人久久爱视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲三级黄色毛片| 91在线精品国自产拍蜜月| 人妻夜夜爽99麻豆av| 精品日产1卡2卡| 波多野结衣高清作品| 欧美+日韩+精品| 男插女下体视频免费在线播放| 亚洲美女黄片视频| 淫妇啪啪啪对白视频| 亚洲av第一区精品v没综合| 国产精品一区二区性色av| 成人二区视频| 国产av不卡久久| 黄色一级大片看看| 熟女人妻精品中文字幕| 日韩制服骚丝袜av| 免费一级毛片在线播放高清视频| 天堂√8在线中文| 亚洲丝袜综合中文字幕| 97碰自拍视频| 久久精品国产清高在天天线| 又爽又黄无遮挡网站| 日韩一本色道免费dvd| 51国产日韩欧美| 人人妻人人澡欧美一区二区| 午夜激情欧美在线| 午夜福利在线在线| 国产亚洲欧美98| 欧美日本亚洲视频在线播放| 熟女人妻精品中文字幕| 我的老师免费观看完整版| 18禁在线无遮挡免费观看视频 | 麻豆国产97在线/欧美| 日韩欧美在线乱码| а√天堂www在线а√下载| 精品少妇黑人巨大在线播放 | 尾随美女入室| 舔av片在线| 亚洲最大成人中文| 免费观看人在逋| 99热这里只有是精品50| 免费电影在线观看免费观看| 国产久久久一区二区三区| 免费av毛片视频| 国产高清三级在线| 亚洲自偷自拍三级| 欧美高清成人免费视频www| 国产黄色小视频在线观看| 亚洲色图av天堂| 成人无遮挡网站| 国产亚洲精品久久久久久毛片| 久久国内精品自在自线图片| 中文字幕av成人在线电影| 日本撒尿小便嘘嘘汇集6| 搡老岳熟女国产| 成熟少妇高潮喷水视频| 免费不卡的大黄色大毛片视频在线观看 | 99热6这里只有精品| 美女免费视频网站| 中文字幕人妻熟人妻熟丝袜美| 天堂√8在线中文| 国产亚洲精品av在线| 色播亚洲综合网| 亚洲中文字幕一区二区三区有码在线看| 久久精品夜色国产| 久久国产乱子免费精品| 国产在视频线在精品| 色综合站精品国产| 久久久久国产网址| 毛片女人毛片| 蜜臀久久99精品久久宅男| 精品99又大又爽又粗少妇毛片| 亚洲av美国av| 国产精品不卡视频一区二区| 免费搜索国产男女视频| 免费观看在线日韩| 久久久久性生活片| 久久久成人免费电影| av专区在线播放| 两个人视频免费观看高清| 国产精品久久久久久久电影| 热99re8久久精品国产| 亚洲精品日韩av片在线观看| 久久久久精品国产欧美久久久| 国产一区二区在线av高清观看| 两性午夜刺激爽爽歪歪视频在线观看|