• <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)鑄造非枝晶組織與模鍛成形性能
    不同形狀橫向限制對枝晶間距影響作用的相場法模擬
    合金定向凝固一次枝晶間距模擬
    男女边摸边吃奶| 又爽又黄无遮挡网站| 免费电影在线观看免费观看| 婷婷六月久久综合丁香| 国产男人的电影天堂91| 国产黄a三级三级三级人| 亚洲av电影在线观看一区二区三区 | 国产亚洲av片在线观看秒播厂 | 纵有疾风起免费观看全集完整版 | 欧美成人一区二区免费高清观看| 男插女下体视频免费在线播放| 亚洲国产精品成人久久小说| 欧美 日韩 精品 国产| 国产日韩欧美在线精品| 亚洲人成网站在线观看播放| 国产精品久久久久久精品电影小说 | 街头女战士在线观看网站| 边亲边吃奶的免费视频| 国产免费一级a男人的天堂| 国产女主播在线喷水免费视频网站 | 精品午夜福利在线看| 亚洲av福利一区| 国产精品日韩av在线免费观看| 国产精品三级大全| 免费看美女性在线毛片视频| 在线观看免费高清a一片| 亚洲精品成人av观看孕妇| 欧美一级a爱片免费观看看| 最近中文字幕2019免费版| 日韩中字成人| 欧美成人一区二区免费高清观看| 日本一本二区三区精品| 麻豆久久精品国产亚洲av| 久久久久久久久久黄片| 性插视频无遮挡在线免费观看| 久久99热6这里只有精品| 18+在线观看网站| 国内少妇人妻偷人精品xxx网站| 日韩欧美三级三区| 毛片女人毛片| 国产欧美日韩精品一区二区| 91久久精品电影网| av一本久久久久| 精品人妻熟女av久视频| 亚洲精华国产精华液的使用体验| 亚洲精品,欧美精品| 淫秽高清视频在线观看| 国产视频内射| 伦精品一区二区三区| 日本与韩国留学比较| 嫩草影院精品99| 免费在线观看成人毛片| 一个人看视频在线观看www免费| 久久国产乱子免费精品| 69av精品久久久久久| 中文字幕亚洲精品专区| 黄色日韩在线| 国产成年人精品一区二区| 美女cb高潮喷水在线观看| av黄色大香蕉| 日本-黄色视频高清免费观看| 国产精品国产三级国产av玫瑰| 欧美日韩国产mv在线观看视频 | 久久精品久久久久久久性| 国产又色又爽无遮挡免| 国产真实伦视频高清在线观看| 欧美性猛交╳xxx乱大交人| 性色avwww在线观看| 亚洲高清免费不卡视频| 午夜日本视频在线| 欧美精品一区二区大全| 国产一区二区三区av在线| 少妇被粗大猛烈的视频| 国产亚洲91精品色在线| 日本免费a在线| 大香蕉97超碰在线| 干丝袜人妻中文字幕| 国产久久久一区二区三区| 美女被艹到高潮喷水动态| 日韩国内少妇激情av| 免费观看在线日韩| 男人舔奶头视频| 日本wwww免费看| 日韩强制内射视频| 一本久久精品| a级毛色黄片| 99热这里只有是精品50| 国产精品1区2区在线观看.| 两个人的视频大全免费| 午夜激情久久久久久久| 欧美丝袜亚洲另类| 内射极品少妇av片p| 亚洲精品成人久久久久久| 国产高清有码在线观看视频| 久久精品国产鲁丝片午夜精品| 国产高清不卡午夜福利| 久久久久久久久久黄片| 日本黄色片子视频| 人人妻人人澡人人爽人人夜夜 | 亚洲精品日本国产第一区| 亚洲av成人精品一区久久| 国产成人精品久久久久久| 欧美精品国产亚洲| 午夜激情欧美在线| 亚洲精品一二三| 久久99精品国语久久久| 观看美女的网站| 女的被弄到高潮叫床怎么办| av天堂中文字幕网| 久久久精品94久久精品| 18禁动态无遮挡网站| 国产男女超爽视频在线观看| 日韩欧美三级三区| 成年免费大片在线观看| 大香蕉久久网| 777米奇影视久久| 午夜激情欧美在线| 国产精品1区2区在线观看.| 爱豆传媒免费全集在线观看| 麻豆乱淫一区二区| 高清视频免费观看一区二区 | 国产日韩欧美在线精品| 日韩av不卡免费在线播放| 青春草亚洲视频在线观看| 国产v大片淫在线免费观看| 三级经典国产精品| 精品久久久久久久久av| 亚洲精品第二区| 国产成人freesex在线| 日韩精品有码人妻一区| 国产亚洲精品久久久com| 国产淫语在线视频| 2022亚洲国产成人精品| 最后的刺客免费高清国语| 国产日韩欧美在线精品| 永久免费av网站大全| 国产成人福利小说| 国产视频内射| 亚洲电影在线观看av| 国产成人精品婷婷| 国产精品一区www在线观看| 国产亚洲精品av在线| 日韩伦理黄色片| 午夜老司机福利剧场| 欧美97在线视频| 国产av在哪里看| 一区二区三区四区激情视频| 午夜精品一区二区三区免费看| 91午夜精品亚洲一区二区三区| 国产精品女同一区二区软件| 真实男女啪啪啪动态图| 日本与韩国留学比较| 18禁在线无遮挡免费观看视频| 国产 一区 欧美 日韩| 国产免费一级a男人的天堂| 搡老乐熟女国产| 亚洲精品国产av蜜桃| 欧美区成人在线视频| 日本猛色少妇xxxxx猛交久久| 国产精品国产三级国产专区5o| 天天一区二区日本电影三级| 在线免费观看的www视频| 成人av在线播放网站| 久久久久久久久久成人| 国产精品一区二区三区四区久久| 国产精品久久久久久久久免| 国产高潮美女av| 深爱激情五月婷婷| 国产精品麻豆人妻色哟哟久久 | 熟妇人妻不卡中文字幕| 日韩av在线免费看完整版不卡| 18禁在线播放成人免费| 男女那种视频在线观看| 亚洲国产精品sss在线观看| 午夜福利在线观看吧| 久久这里有精品视频免费| 亚洲精品成人久久久久久| 99久国产av精品国产电影| 人人妻人人看人人澡| 3wmmmm亚洲av在线观看| 日韩亚洲欧美综合| 99久久九九国产精品国产免费| 亚洲精品色激情综合| 少妇人妻精品综合一区二区| 最近中文字幕2019免费版| 久久久精品94久久精品| 久久久a久久爽久久v久久| 日韩av在线免费看完整版不卡| 一二三四中文在线观看免费高清| 国产探花在线观看一区二区| 亚洲欧美成人综合另类久久久| 国产片特级美女逼逼视频| 国产伦在线观看视频一区| 亚洲欧洲国产日韩| 久久99热6这里只有精品| 精品国产露脸久久av麻豆 | 午夜精品国产一区二区电影 | 国产毛片a区久久久久| 国产精品久久久久久久电影| 日日撸夜夜添| 亚洲国产最新在线播放| 又爽又黄a免费视频| 日韩国内少妇激情av| 91av网一区二区| 91久久精品国产一区二区成人| 极品教师在线视频| 亚洲,欧美,日韩| 91aial.com中文字幕在线观看| 日本爱情动作片www.在线观看| 国产淫片久久久久久久久| 真实男女啪啪啪动态图| 欧美最新免费一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品久久午夜乱码| 麻豆av噜噜一区二区三区| 亚洲av中文字字幕乱码综合| 日本免费在线观看一区| 欧美bdsm另类| 日韩av在线大香蕉| 精品久久久久久成人av| 看免费成人av毛片| 国产精品国产三级专区第一集| 精品午夜福利在线看| 美女cb高潮喷水在线观看| 国产 一区 欧美 日韩| 一级毛片黄色毛片免费观看视频| 久久99热这里只有精品18| 卡戴珊不雅视频在线播放| 亚洲在线自拍视频| 日产精品乱码卡一卡2卡三| 久久99热这里只频精品6学生| 最近最新中文字幕大全电影3| 亚洲欧美日韩无卡精品| 乱系列少妇在线播放| 伦理电影大哥的女人| 麻豆国产97在线/欧美| 欧美高清性xxxxhd video| 日韩精品青青久久久久久| 插阴视频在线观看视频| 欧美激情久久久久久爽电影| 亚洲av电影在线观看一区二区三区 | 极品少妇高潮喷水抽搐| 深夜a级毛片| 在线免费十八禁| www.av在线官网国产| 精品久久久久久久久av| 天堂网av新在线| 免费播放大片免费观看视频在线观看| 三级毛片av免费| 插逼视频在线观看| 日韩av在线大香蕉| 天天一区二区日本电影三级| 国产精品女同一区二区软件| 国产在线男女| 国产男人的电影天堂91| av在线老鸭窝| 一级毛片 在线播放| 蜜桃亚洲精品一区二区三区| av在线观看视频网站免费| 一本一本综合久久| 高清午夜精品一区二区三区| 国产成人a∨麻豆精品| 午夜免费男女啪啪视频观看| 麻豆国产97在线/欧美| 日韩欧美国产在线观看| 3wmmmm亚洲av在线观看| 国产精品伦人一区二区| 边亲边吃奶的免费视频| 97热精品久久久久久| 日本熟妇午夜| 亚洲国产欧美人成| 精品久久久噜噜| 久久精品久久久久久噜噜老黄| 国产麻豆成人av免费视频| 欧美日韩一区二区视频在线观看视频在线 | 日本猛色少妇xxxxx猛交久久| 尤物成人国产欧美一区二区三区| 久久精品国产自在天天线| a级毛色黄片| 国产女主播在线喷水免费视频网站 | 日韩成人伦理影院| 亚洲最大成人av| 成人综合一区亚洲| 亚洲国产欧美在线一区| 国产高清三级在线| 亚洲一级一片aⅴ在线观看| 中文字幕制服av| 免费看a级黄色片| 成人av在线播放网站| 九九在线视频观看精品| av卡一久久| 亚洲性久久影院| 九色成人免费人妻av| 在线观看人妻少妇| 亚洲欧洲国产日韩| 九九爱精品视频在线观看| 国产黄频视频在线观看| 欧美xxxx性猛交bbbb| 永久网站在线| 2021少妇久久久久久久久久久| 亚洲av二区三区四区| 偷拍熟女少妇极品色| 天堂av国产一区二区熟女人妻| 亚洲精品色激情综合| 日日啪夜夜爽| 国产伦精品一区二区三区视频9| 国产精品久久久久久久久免| 91av网一区二区| 韩国高清视频一区二区三区| 神马国产精品三级电影在线观看| 精品人妻熟女av久视频| 在线观看av片永久免费下载| 亚洲四区av| av天堂中文字幕网| 欧美日韩国产mv在线观看视频 | 伦精品一区二区三区| 精品一区在线观看国产| 久久久精品免费免费高清| 国产精品一区二区三区四区久久| 夜夜爽夜夜爽视频| 午夜激情久久久久久久| 成年女人在线观看亚洲视频 | 99久国产av精品国产电影| 精品不卡国产一区二区三区| 一区二区三区乱码不卡18| 国产毛片a区久久久久| 日韩成人av中文字幕在线观看| 亚洲精华国产精华液的使用体验| av一本久久久久| 美女被艹到高潮喷水动态| 五月玫瑰六月丁香| 大又大粗又爽又黄少妇毛片口| 美女国产视频在线观看| 久久精品国产鲁丝片午夜精品| 久久精品久久久久久噜噜老黄| 一级av片app| 免费高清在线观看视频在线观看| 秋霞在线观看毛片| 久久久久久九九精品二区国产| 禁无遮挡网站| 我的老师免费观看完整版| 久久午夜福利片| 国产精品.久久久| 两个人视频免费观看高清| 人妻一区二区av| 国产视频首页在线观看| 尾随美女入室| 免费大片18禁| 欧美高清成人免费视频www| 午夜激情久久久久久久| 天堂av国产一区二区熟女人妻| 九色成人免费人妻av| 中文字幕久久专区| 一级a做视频免费观看| 国产成人a区在线观看| 久久久久久久久大av| 久久精品国产亚洲av天美| 国产一区二区在线观看日韩| 丝瓜视频免费看黄片| 亚洲欧洲国产日韩| videos熟女内射| 我的女老师完整版在线观看| 最近最新中文字幕大全电影3| 亚洲av.av天堂| 午夜久久久久精精品| 国产不卡一卡二| 国产一级毛片在线| 久久精品久久久久久噜噜老黄| 国产91av在线免费观看| 国产精品久久久久久精品电影| 高清毛片免费看| 久久久久精品久久久久真实原创| 精品亚洲乱码少妇综合久久| 一级毛片 在线播放| 最近最新中文字幕大全电影3| 国产 亚洲一区二区三区 | 国产91av在线免费观看| 精品一区二区免费观看| 国产在线一区二区三区精| 亚洲综合色惰| 欧美一级a爱片免费观看看| 亚洲av福利一区| 免费人成在线观看视频色| 一级a做视频免费观看| 80岁老熟妇乱子伦牲交| 免费看美女性在线毛片视频| 国内精品一区二区在线观看| 久久韩国三级中文字幕| 亚洲国产日韩欧美精品在线观看| 国产 亚洲一区二区三区 | 女人十人毛片免费观看3o分钟| 婷婷色综合www| 色5月婷婷丁香| 亚洲av.av天堂| 免费看av在线观看网站| 寂寞人妻少妇视频99o| 午夜亚洲福利在线播放| 人人妻人人看人人澡| 国产成人精品婷婷| 真实男女啪啪啪动态图| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美精品一区二区大全| 一夜夜www| 国产午夜精品论理片| 丰满乱子伦码专区| 一级爰片在线观看| 国产精品福利在线免费观看| av一本久久久久| 午夜视频国产福利| 丰满人妻一区二区三区视频av| 九九爱精品视频在线观看| 岛国毛片在线播放| 欧美一区二区亚洲| 秋霞伦理黄片| 日本午夜av视频| 能在线免费观看的黄片| 国产熟女欧美一区二区| 美女内射精品一级片tv| 特级一级黄色大片| 国产视频首页在线观看| 国产免费又黄又爽又色| 国产成人91sexporn| 狠狠精品人妻久久久久久综合| 日产精品乱码卡一卡2卡三| 2021少妇久久久久久久久久久| 亚洲成人av在线免费| av卡一久久| 三级经典国产精品| 国精品久久久久久国模美| 色5月婷婷丁香| 深爱激情五月婷婷| 联通29元200g的流量卡| 亚洲精品乱久久久久久| 午夜精品国产一区二区电影 | 国产成人精品久久久久久| 一区二区三区四区激情视频| 欧美97在线视频| 亚洲av成人精品一二三区| 亚洲精品中文字幕在线视频 | 亚洲av国产av综合av卡| 亚洲精品国产av成人精品| 国产高清不卡午夜福利| 91久久精品国产一区二区成人| 日日摸夜夜添夜夜添av毛片| 国产亚洲最大av| 精品国内亚洲2022精品成人| 丰满人妻一区二区三区视频av| 亚洲av在线观看美女高潮| 美女大奶头视频| 亚洲欧美中文字幕日韩二区| 99久久九九国产精品国产免费| 非洲黑人性xxxx精品又粗又长| 国产精品伦人一区二区| 亚洲欧美日韩无卡精品| 成人性生交大片免费视频hd| 亚洲av二区三区四区| 欧美一区二区亚洲| 波多野结衣巨乳人妻| 亚洲自偷自拍三级| 亚洲精品一区蜜桃| 日本色播在线视频| 亚洲四区av| 日日摸夜夜添夜夜添av毛片| 精品午夜福利在线看| 久久久久性生活片| 免费观看av网站的网址| 青春草亚洲视频在线观看| 男人狂女人下面高潮的视频| 国产极品天堂在线| 能在线免费看毛片的网站| 亚洲在久久综合| av免费观看日本| 国产亚洲av嫩草精品影院| 国产成人一区二区在线| 国产精品日韩av在线免费观看| 亚洲图色成人| 汤姆久久久久久久影院中文字幕 | 国产乱人偷精品视频| 亚洲av中文av极速乱| 久久久色成人| 国产精品一区二区在线观看99 | 国产单亲对白刺激| 国产一区二区在线观看日韩| 99热这里只有是精品50| 深爱激情五月婷婷| 精品国产露脸久久av麻豆 | 天堂中文最新版在线下载 | 久久精品国产鲁丝片午夜精品| 国产国拍精品亚洲av在线观看| 国产伦在线观看视频一区| 欧美xxxx黑人xx丫x性爽| 22中文网久久字幕| 精品99又大又爽又粗少妇毛片| 亚洲国产欧美在线一区| 内射极品少妇av片p| 亚洲va在线va天堂va国产| 亚洲av国产av综合av卡| 国产久久久一区二区三区| 国产成人福利小说| 国产精品一区二区三区四区免费观看| 国内精品宾馆在线| 亚洲国产精品sss在线观看| 最近最新中文字幕大全电影3| 国产精品爽爽va在线观看网站| 看十八女毛片水多多多| 中文乱码字字幕精品一区二区三区 | 男插女下体视频免费在线播放| 国产黄色免费在线视频| 国产乱人偷精品视频| 亚洲精品乱久久久久久| 日韩视频在线欧美| 久久久久久久午夜电影| 亚洲经典国产精华液单| 午夜精品在线福利| 国产爱豆传媒在线观看| 日韩一区二区三区影片| 久久这里只有精品中国| 亚洲人成网站在线观看播放| 亚洲婷婷狠狠爱综合网| 国产黄色免费在线视频| 又爽又黄a免费视频| 国产视频首页在线观看| 亚洲成人av在线免费| 99热6这里只有精品| 日本色播在线视频| 欧美成人午夜免费资源| 成人国产麻豆网| 看非洲黑人一级黄片| 亚洲欧美清纯卡通| 99视频精品全部免费 在线| 69av精品久久久久久| 久久久久国产网址| 日日干狠狠操夜夜爽| 国内精品一区二区在线观看| 亚洲内射少妇av| 久久这里有精品视频免费| 国国产精品蜜臀av免费| 精品久久久久久久人妻蜜臀av| 国产国拍精品亚洲av在线观看| 免费观看精品视频网站| 麻豆国产97在线/欧美| 日韩欧美一区视频在线观看 | 天天躁日日操中文字幕| 夫妻性生交免费视频一级片| 欧美成人一区二区免费高清观看| 国产有黄有色有爽视频| 免费观看的影片在线观看| 欧美区成人在线视频| 免费观看的影片在线观看| 亚洲电影在线观看av| 乱系列少妇在线播放| 国产亚洲5aaaaa淫片| 国产69精品久久久久777片| 精品久久久久久电影网| 夜夜爽夜夜爽视频| 18禁在线播放成人免费| 人妻制服诱惑在线中文字幕| 亚洲高清免费不卡视频| 午夜日本视频在线| 丝袜喷水一区| 亚洲精品国产av成人精品| 欧美日韩在线观看h| 精品一区二区三卡| 久久韩国三级中文字幕| 久久久午夜欧美精品| 观看免费一级毛片| 国产综合精华液| 国产精品一区二区性色av| 晚上一个人看的免费电影| 国产亚洲午夜精品一区二区久久 | 国产成人免费观看mmmm| 欧美+日韩+精品| 男人爽女人下面视频在线观看| 欧美日韩精品成人综合77777| 国产激情偷乱视频一区二区| 亚洲精品影视一区二区三区av| 日韩av在线免费看完整版不卡| 乱码一卡2卡4卡精品| 午夜亚洲福利在线播放| 成人毛片a级毛片在线播放| 中文字幕免费在线视频6| 国产一区二区亚洲精品在线观看| 亚洲成人av在线免费| 免费看av在线观看网站| 国产黄a三级三级三级人| av天堂中文字幕网| 成人鲁丝片一二三区免费| 少妇的逼好多水| 老司机影院毛片| 婷婷色av中文字幕| 亚洲国产欧美人成| 高清午夜精品一区二区三区| 特大巨黑吊av在线直播| 一级a做视频免费观看| 国产老妇女一区| 99久国产av精品| 男女国产视频网站| 永久网站在线| 亚洲图色成人| 观看美女的网站| 久久精品人妻少妇| 十八禁国产超污无遮挡网站| 日日摸夜夜添夜夜爱| 国产免费福利视频在线观看| 建设人人有责人人尽责人人享有的 | a级一级毛片免费在线观看| 少妇裸体淫交视频免费看高清| 九九在线视频观看精品| 人妻系列 视频| 亚洲综合色惰| 亚洲国产日韩欧美精品在线观看| 午夜福利在线观看吧| 日日啪夜夜撸| 精品亚洲乱码少妇综合久久| 一级a做视频免费观看|