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

    微小水滴撞擊深水液池空腔運(yùn)動的數(shù)值模擬及機(jī)理研究?

    2018-12-18 05:58:20裴傳康魏炳乾
    物理學(xué)報(bào) 2018年22期
    關(guān)鍵詞:液池空腔液面

    裴傳康 魏炳乾

    (西安理工大學(xué),省部共建西北旱區(qū)生態(tài)水利國家重點(diǎn)實(shí)驗(yàn)室,西安 710048)

    (2018年7月28日收到;2018年9月20日收到修改稿)

    1 引 言

    液滴下落沖擊不同介質(zhì)的運(yùn)動過程在科技應(yīng)用及自然現(xiàn)象中廣泛存在.液滴沖擊液體表面的研究在犯罪取證、開發(fā)不可浸潤表面或完全可浸潤表面、高精度活化或表面污染物轉(zhuǎn)移等方面有諸多應(yīng)用[1],液-液界面的部分聚結(jié)過程也是許多復(fù)雜物理現(xiàn)象的組成部分,與地球物理學(xué)與土力學(xué)等有諸多關(guān)聯(lián)[2,3];對于海洋、湖泊等大面積水面,其廣泛的自然摻氣現(xiàn)象主要取決于水滴撞擊所引發(fā)的氣泡夾帶[4];在水利工程中,高水頭泄水建筑物的霧化、摻氣、消能等問題也與此息息相關(guān).因此,研究此類基本運(yùn)動對于理解自然界氣液流動的界面變形、改善液滴運(yùn)動在工程中的應(yīng)用具有重要的意義與作用.

    1963年,Worthington[5]通過實(shí)驗(yàn)首次系統(tǒng)地描述了液滴和固體小球撞擊液池的過程及其運(yùn)動規(guī)律.隨后的大量研究表明,當(dāng)液滴以較低的撞擊速度撞擊液面時,液滴與液面發(fā)生完全聚結(jié)現(xiàn)象,并在液面下方生成渦環(huán)[6,7],在這種情況下,自由液面并不會出現(xiàn)飛濺射流,而是呈現(xiàn)出平坦的狀態(tài).隨著液滴撞擊速度的增加,自由液面開始飛濺,液面出現(xiàn)較大變形,產(chǎn)生一個空腔,并在四周形成皇冠狀的射流[8],空腔塌陷后,液柱從撞擊中心升起,速度較大時,由于Plateau-Rayleigh不穩(wěn)定性的影響,會在液柱上方分離出一個或多個小液滴[9].一般情況下,液滴沖擊液池的運(yùn)動主要受重力、慣性力與表面張力的影響,因此對于弗勞德數(shù)(Fr)和韋伯?dāng)?shù)(We)十分敏感[10].

    傳統(tǒng)上對液滴運(yùn)動的探索多是以高速攝像機(jī)為主要工具進(jìn)行的實(shí)驗(yàn)研究,但由于液滴運(yùn)動的復(fù)雜性,實(shí)驗(yàn)攝影方式較難捕捉細(xì)致的自由面邊界,對諸如壓力和速度等物理量的測定也十分不便.隨著Hirt與Nicolas在1981年[11]提出利用流體體積法(VOF)分離兩種流體構(gòu)成的尖銳界面后,數(shù)值模擬在捕捉界面變形方面精度不斷提升[12,13],逐漸成為了研究液滴撞擊運(yùn)動的有力工具.Yue等[14]對液滴撞擊液面運(yùn)動進(jìn)行了基于相場法的數(shù)值模擬,并根據(jù)最大Oh數(shù)得到了部分融合發(fā)生的臨界標(biāo)準(zhǔn).Ray等[15]利用CLSVOF法對液滴撞擊液面進(jìn)行了更為具體的研究,結(jié)果表明慣性力和表面張力在很大程度上決定了液面的運(yùn)動過程.Castillo-Orozco等[16]基于VOF法和表面張力模型(CSF),討論了流體的物理化學(xué)性質(zhì)和液滴沖擊速度變化對二次液滴的斷裂以及冠狀射流的影響.戴劍鋒等[17]應(yīng)用CLSVOF法研究了液滴撞擊傾斜表面液膜后液膜的形態(tài)演化及飛濺過程,研究表明,空氣卷吸氣泡數(shù)量隨著沖擊角的增大而減小.黃虎等[18]則建立了液滴撞擊平面液膜的數(shù)學(xué)模型,采用格子Boltzmann方法探討了相對液膜厚度和表面張力等物理參數(shù)對界面運(yùn)動過程的影響.

    以往對液滴撞擊運(yùn)動的研究多針對毫米級直徑的液滴[19],驅(qū)動其空腔形成的要素主要為慣性力及重力[10],但在較小液滴尺度下,由于液滴質(zhì)量及體積的減小,毛細(xì)波運(yùn)動在驅(qū)動液-液界面變形及空腔形成過程中起到越發(fā)重要的作用.現(xiàn)有研究較少關(guān)注直徑在微米層級的液滴運(yùn)動及空腔形成過程以及決定微小液滴空腔運(yùn)動的主要驅(qū)動力,因此,本文主要關(guān)注高速運(yùn)動的微小水滴,運(yùn)用基于四叉樹結(jié)構(gòu)自適應(yīng)網(wǎng)格、VOF方法以及變密度不可壓縮Navier-Stocks方程的開源多相流程序Gerris[20?22],研究一定物理參數(shù)范圍下微小水滴撞擊深水液池的液面變形及毛細(xì)波運(yùn)動過程,并嘗試揭示水滴撞擊產(chǎn)生的空腔運(yùn)動及形成機(jī)理.

    2 計(jì)算區(qū)域及數(shù)值方法

    2.1 控制方程

    液滴撞擊液池的運(yùn)動過程可以使用帶有表面張力項(xiàng)的變密度、不可壓縮Navier-Stocks方程來描述,具體控制方程如下[20,21]:

    式中,ρ=ρ(x,t)為流體密度,u=(u,v,w)為流體速度,p為壓力,μ=μ(x,t)是流體的動力黏度,變形張量D定義為Dij=(?iuj+?jui)/2,σ為表面張力系數(shù),κ為界面曲率,狄拉克分布函數(shù)δs表示表面張力僅作用于兩相界面處,n為兩相界面的法向量.

    Gerris采用經(jīng)典的VOF方法追蹤相界面,對于兩相流動,引入計(jì)算網(wǎng)格中第一種流體的體積分?jǐn)?shù)c(x,t),并定義混合流體的密度和黏度為:

    式中,ρ1,ρ2,μ1,μ2分別是第一種流體和第二種流體的密度以及黏度;函數(shù)?c由體積分?jǐn)?shù)c平滑處理后得出,以便提高計(jì)算的穩(wěn)定性.

    密度對流方程可由等效的體積分?jǐn)?shù)對流方程替換

    2.2 數(shù)值方法

    本文數(shù)值模擬采用基于Linux的開源軟件Gerris進(jìn)行,該軟件使用基于四叉樹(二維)/八叉樹(三維)的自適應(yīng)空間離散方法,使用分步投影方法求解變密度不可壓縮的Navier-Stocks方程,使用VOF方法跟蹤相界面.高度函數(shù)和界面附近的自適應(yīng)網(wǎng)格細(xì)化可以精確表示表面張力作用,對流項(xiàng)使用Godunov格式求解,并行計(jì)算采用MPI庫進(jìn)行.

    如圖1所示,水滴撞擊深水液池的數(shù)值模擬在軸對稱坐標(biāo)系中進(jìn)行,Y軸為計(jì)算區(qū)域的對稱軸,D為初始水滴直徑,R=D/2,正方形計(jì)算區(qū)域的長度H=20D,水滴距離液面的距離H1=0.1D,液池深度H2=12D以消除底部對液滴撞擊運(yùn)動的影響,水滴在重力g和撞擊速度Vi的作用下撞擊液池.采用雷諾數(shù)、韋伯?dāng)?shù)和弗勞德數(shù)來描述液滴撞擊的運(yùn)動特征,三者分別表征液體慣性力與黏滯力間的關(guān)系、液體慣性力與表面張力間的關(guān)系以及液體慣性力與重力間的關(guān)系.三個無量綱參數(shù)的表達(dá)式如(7)式所示,主要物理參數(shù)如表1所列.

    表1 主要物理參數(shù)Table 1.Main physical parameters.

    圖1 計(jì)算區(qū)域簡圖Fig.1.Schematic diagram of the computational domain.

    2.3 自適應(yīng)網(wǎng)格技術(shù)

    采用數(shù)值方法對高速液滴的撞擊運(yùn)動進(jìn)行準(zhǔn)確的模擬極具挑戰(zhàn)性,因?yàn)檫\(yùn)動產(chǎn)生的微小界面變形、復(fù)雜的幾何形狀以及特征尺度的巨大差異需要足夠的網(wǎng)格分辨率來捕捉,從而大幅地增加了計(jì)算量與計(jì)算時間.目前針對該問題的一個有效解決方法是采用自適應(yīng)網(wǎng)格(adaptive mesh refinement,AMR)技術(shù).根據(jù)流動特征對網(wǎng)格進(jìn)行局部細(xì)化或粗化使得AMR技術(shù)可以將計(jì)算效率集中在最需要的區(qū)域,從而以最小的計(jì)算成本獲取精確的結(jié)果.

    本文采用Gerris進(jìn)行數(shù)值模擬,Gerris使用有限體積法(FVM)來求解控制方程,并根據(jù)四叉樹網(wǎng)格自適應(yīng)規(guī)則和條件將計(jì)算域離散為不同等級的計(jì)算網(wǎng)格.水滴撞擊深水液池?cái)?shù)值模擬的關(guān)鍵位置在于液-液界面的交接處以及相界面附近,本文依此設(shè)計(jì)如下網(wǎng)格自適應(yīng)規(guī)則,每一步更新一次計(jì)算網(wǎng)格,其中最大網(wǎng)格加密層數(shù)為11層,即在一個計(jì)算區(qū)域(box,Lbox=10)內(nèi)的最大網(wǎng)格數(shù)量為211.圖2為計(jì)算區(qū)域初始狀態(tài)自適應(yīng)網(wǎng)格的空間離散示意.

    1)計(jì)算初始加密水滴與液池接觸區(qū)域,即水滴與液池相界面處正負(fù)0.15內(nèi)的網(wǎng)格至11層.

    2)自動加密相界面附近體積分?jǐn)?shù)在0—1之間、梯度變化劇烈區(qū)域的網(wǎng)格,最大加密到11層,最小加密到6層,以最小化界面重建產(chǎn)生的誤差.

    3)自動加密渦量變化區(qū)域的網(wǎng)格,根據(jù)其變化劇烈程度最大加密到11層,最小加密到4層.

    4)根據(jù)U,V速度分量的變化自動加密網(wǎng)格,最大加密到11層.

    5)限制2)—4)條規(guī)則最小加密層數(shù)的加密區(qū)間為:Y向液滴中心上方2R至水面下方4R;x向?qū)ΨQ軸左右4R內(nèi)的矩形內(nèi),以提高計(jì)算效率.

    圖2 計(jì)算初始狀態(tài)的空間離散Fig.2.The spatial discretization of the initial simulation.

    2.4 模型驗(yàn)證與率定

    為了保證數(shù)值模擬結(jié)果的準(zhǔn)確性,本文選擇Morton的實(shí)驗(yàn)數(shù)據(jù)[23]對數(shù)值模型進(jìn)行驗(yàn)證,實(shí)驗(yàn)使用直徑為2.9 mm的液滴撞擊液池,弗勞德數(shù)及韋伯?dāng)?shù)分別為Fr=220,We=248.如圖3所示,照片為高速攝影機(jī)攝得的實(shí)驗(yàn)過程,白色線條表示相同時間節(jié)點(diǎn)下的數(shù)值模擬結(jié)果,t為物理時間乘以Vi/D后的無量綱時間.液滴下落后沖擊液池并產(chǎn)生了一個空腔,腔體在t=7.9時達(dá)到最大化.空腔塌陷后毛細(xì)波向中心處傳遞,并坍縮形成中心射流,使其高度不斷增大,在射流頂端斷裂生成二次液滴.由于實(shí)驗(yàn)環(huán)境的復(fù)雜性,模擬條件與實(shí)驗(yàn)條件無法完全一致,且本文采用軸對稱模型假定進(jìn)行模擬,無法捕捉非對稱運(yùn)動,因此模擬值與實(shí)驗(yàn)值存在一定差異.但數(shù)值模擬在界面變形、空腔的形成與成長、毛細(xì)波在空腔底部的傳播等方面與實(shí)驗(yàn)值取得了良好的一致性,且在空腔形成過程中給出了較實(shí)驗(yàn)更加詳盡的毛細(xì)波運(yùn)動細(xì)節(jié),中心射流最大高度以及空腔最大深度的誤差分別為1.7%,2.6%,表明數(shù)值模擬能夠較好地描述液滴撞擊液池的運(yùn)動.

    圖3 數(shù)值模擬與實(shí)驗(yàn)結(jié)果對比Fig.3.Comparison of the numerical and experimental results.

    3 計(jì)算結(jié)果與分析

    3.1 基本撞擊過程

    水滴撞擊深水液池后相界面的變化可以直觀反映出水滴撞擊的運(yùn)動特性、水滴與池水的混摻過程以及撞擊所夾帶的氣泡大小.因此,使用經(jīng)過驗(yàn)證的數(shù)值方法研究直徑為290μm的水滴以2.5—6.5 m/s的速度撞擊深水液池的運(yùn)動過程.

    圖4 不同工況下自由液面隨時間的運(yùn)動過程 (a)Fr=117.2,Re=725,We=25.2,Vi=2.5 m/s;(b)Fr=229.7,Re=1015,We=49.3,Vi=3.5 m/s;(c)Fr=379.7,Re=1305,We=81.6,Vi=4.5 m/s;(d)Fr=567.1,Re=1595,We=121.8,Vi=5.5 m/s;(e)Fr=792.1,Re=1885,We=170.2,Vi=6.5 m/sFig.4.Free surface profiles simulated at selected times:(a)Fr=117.2,Re=725,We=25.2,Vi=2.5 m/s;(b)Fr=229.7,Re=1015,We=49.3,Vi=3.5 m/s;(c)Fr=379.7,Re=1305,We=81.6,Vi=4.5 m/s;(d)Fr=567.1,Re=1595,We=121.8,Vi=5.5 m/s;(e)Fr=792.1,Re=1885,We=170.2,Vi=6.5 m/s.

    水滴以五種不同撞擊速度進(jìn)入水池后自由表面的模擬輪廓以及水滴分布如圖4所示.藍(lán)色為氣相,紅色為水池部分液相,天藍(lán)色為水滴部分液相,從工況(a)—(e),水滴的撞擊速度不斷增大.在圖4(a)中,Fr=117.2,We=25.2,Re=725,水滴以較低速度撞擊深水液池,發(fā)生完全聚結(jié)現(xiàn)象,并在水池中形成底部夾帶一個氣泡的空腔,水滴入水后基本附著在自由液面之下.隨著時間增加,空腔開始向內(nèi)坍縮,同時毛細(xì)波向空腔底部傳遞,使得水滴部分的水體向中心聚合并產(chǎn)生兩個對稱的渦,最終在空腔塌陷后對池水產(chǎn)生穿刺效應(yīng),池水在慣性力的作用下逐漸平復(fù).在圖4(b)—(d)中,撞擊產(chǎn)生的空腔隨著水滴撞擊時動能的增加越來越大,空腔底部夾帶的氣泡不斷變小,隨著撞擊速度的增大,在空腔塌陷后,一個短而粗的射流在撞擊中心產(chǎn)生,且射流最大高度不斷增加.由于毛細(xì)波的傳遞速度加快,水滴入水后附著自由液面的面積也逐漸變大,產(chǎn)生的渦距離撞擊中心越來越遠(yuǎn).在圖4(e)中,Fr=792.1,Re=1885,We=170.2,水滴撞擊后首先產(chǎn)生一個U形空腔,最大深度再次增加,隨后由于毛細(xì)波的快速傳播推動空腔底部部分側(cè)壁收縮閉合,截留形成一個較大的氣泡,并在閉合處產(chǎn)生極細(xì)的射流,射流由于豎向剪應(yīng)力較大,引起Plateau-Rayleigh不穩(wěn)定性,在尖端斷裂生成多個二次滴,自由液面變化更為劇烈.

    上述水滴撞擊水池的過程可以歸納為三個基本階段,第一階段為空腔的形成與膨脹;第二階段為毛細(xì)波傳播導(dǎo)致的空腔收縮;第三階段為空腔的回復(fù).液滴撞擊深水液池時的運(yùn)動分為完全聚結(jié)與部分聚結(jié)現(xiàn)象[24],在水滴與深水液池完全聚結(jié)時,水滴入水時由于水滴底部發(fā)生凹陷變形產(chǎn)生的氣泡夾帶隨著撞擊速度的增加而減小,而在部分聚結(jié)發(fā)生時氣泡的夾帶與截留行為則更為復(fù)雜.

    3.2 空腔運(yùn)動的基本特性

    為了深入探究微米級水滴撞擊深水液池后空腔運(yùn)動的基本特性,對本文五個模擬工況進(jìn)行定量研究.描述空腔幾何特征的基本參數(shù)如圖5所示,以初始靜水狀態(tài)下深水液池的液面高度為基準(zhǔn)線,hw為基準(zhǔn)線至波浪頂端的波浪隆起高度,h為基準(zhǔn)線至空腔最低點(diǎn)的空腔深度,r為空腔基準(zhǔn)線上軸對稱點(diǎn)處至空腔側(cè)壁的空腔寬度.下文涉及的所有長度單位均為實(shí)際長度除以初始水滴直徑D后的無量綱長度,時間t為實(shí)際時間乘以Vi/D后的無量綱時間.

    圖6為不同工況下空腔深度h隨時間的變化過程.由圖6可以看出,液滴撞擊深水液池后,由于動能、空腔側(cè)壁隆起部分重力勢能以及表面張力能的驅(qū)動,空腔深度先以較快速度增長,其后增勢逐漸放緩,在到達(dá)最大空腔深度后快速回彈,至接近原自由液面后回彈速度逐漸放緩,近乎趨于直線.弗勞德數(shù)的增加對空腔深度變化影響顯著,在工況(a)中,弗勞德數(shù)僅為117.2,空腔能量耗散時程較短,回復(fù)較為迅速,隨著弗勞德數(shù)的增大,撞擊動能增大,回復(fù)時程也逐漸增加.值得注意的是,在工況(e)中,由于空腔形狀由U形轉(zhuǎn)變?yōu)榻铺菪?空腔在t=1.5—2.6時先緩慢上升,隨后毛細(xì)波向空腔底部傳遞,促使其底部變?yōu)閳A柱形,且深度繼續(xù)增加,最后空腔底部側(cè)壁坍塌形成射流,因此在t=2.6后空腔回復(fù)速度較其他工況更為迅速.將數(shù)值結(jié)果運(yùn)用最小二乘法擬合,得到的擬合曲線表達(dá)式如(8)式所示,該式能夠在忽略毛細(xì)波運(yùn)動的前提下,在空腔深度為h=D至達(dá)到最大深度的范圍內(nèi)較好地描述了空腔深度隨時間的成長關(guān)系,Liow等[8]與Ray等[24]對液滴撞擊運(yùn)動中時間與空腔深度的關(guān)系進(jìn)行了建模,得出與本文相似的結(jié)論:

    圖5 空腔的幾何特征示意Fig.5.Geometric characteristics of the cavity.

    圖6 不同工況下空腔深度隨時間的變化Fig.6.Cavity depth with time under different conditions.

    圖7及圖8分別為不同工況下波浪隆起高度hw和空腔寬度r隨時間的變化過程.由圖7易得,水滴撞擊水池后波浪高度的變化過程也經(jīng)歷了從快速增長到逐漸回落,最終在慣性力的作用下自主運(yùn)動,趨于穩(wěn)定的過程.隨著弗勞德數(shù)的增大,更大的垂向速度在創(chuàng)造更深空腔的同時激發(fā)了更高的最大波浪隆起高度.從圖8可以看出,空腔寬度的增長與波浪的運(yùn)動息息相關(guān),在撞擊產(chǎn)生的動能與波浪自身重力勢能基本耗散后,空腔運(yùn)動主要由毛細(xì)波及其干擾驅(qū)動,最后在慣性下線性增長.毛細(xì)波現(xiàn)象在圖8工況(e)的t=0.8時非常顯著,其周期在圖6的t=1.5—2.5中也有體現(xiàn).

    圖7 不同工況下波浪隆起高度隨時間的變化Fig.7.Wave height with time under different conditions.

    圖8 不同工況下空腔寬度隨時間的變化Fig.8.Cavity width with time under different conditions.

    3.3 空腔形成以及毛細(xì)波傳播機(jī)理

    選取水滴完全聚結(jié)的工況(d)與水滴部分聚結(jié)的工況(e)對空腔形成以及毛細(xì)波傳播的機(jī)理進(jìn)行深入研究.圖9為不同時間節(jié)點(diǎn)下工況(d)及工況(e)水滴撞擊深水液池空腔運(yùn)動的等值線圖,其中左側(cè)為渦量場等值線圖,右側(cè)為壓力場等值線圖,黑色線條表示相界面.由圖9(a)可得,水滴撞擊液池后,空腔形狀由U形向V形轉(zhuǎn)變,最終腔底首先上升,形成沒有氣泡截留的射流.

    Berberovic等[25]的研究表明,毛細(xì)波的傳播路徑與低壓區(qū)的形成息息相關(guān).在圖9(a)壓力場等值線圖中可以看出,低壓區(qū)首先在空腔側(cè)壁與底部的交界處產(chǎn)生,隨后向空腔底部傳播,并形成一個相對尖銳的點(diǎn),在空腔底部上升后,毛細(xì)波開始向下方傳遞,并逐漸趨于平緩.

    在發(fā)生部分聚結(jié)現(xiàn)象的圖9(b)中,低壓區(qū)首先在波浪底部與側(cè)壁上交界處產(chǎn)生,并已經(jīng)在相界面上形成了一個尖銳的折點(diǎn),此時空腔形狀為半球形,隨后毛細(xì)波向空腔底部傳播,在底部中心空腔轉(zhuǎn)變?yōu)閳A柱狀,低壓區(qū)附著在圓柱下方與底面交點(diǎn),毛細(xì)波在空腔坍塌前并沒有到達(dá)空腔底部中心,而后空腔側(cè)壁坍塌形成氣泡截留,并在中心產(chǎn)生了快速射流液柱.

    渦量定義為流體速度矢量的旋度,描述液體流動的剪切運(yùn)動.在工況(d)中,如圖9(a)渦量場所示,流體渦量一直跟隨毛細(xì)波位置,當(dāng)空腔坍塌產(chǎn)生慢速射流時,渦量在靠近液面區(qū)域以及空腔底部靠近對稱軸的區(qū)域各產(chǎn)生一個較大的渦環(huán),而在工況(e)中,渦環(huán)的生成被抑制,流體僅在t=1.906時表現(xiàn)出一個靠近空腔底部的小渦環(huán),而后渦環(huán)迅速消失.由渦量場與壓力場對比可得,渦量較大的區(qū)域并不總是處于低壓區(qū)內(nèi),撞擊運(yùn)動初始自由液面的壓力差可能是渦量產(chǎn)生的誘因,但其后低壓區(qū)的運(yùn)動與渦量之間并無較強(qiáng)的相關(guān)性.

    4 結(jié) 論

    本文采用基于自適應(yīng)網(wǎng)格和VOF方法的開源程序Gerris對微小水滴撞擊深水液池后的流動過程以及空腔生長進(jìn)行了數(shù)值模擬,研究了不同F(xiàn)r數(shù)對撞擊后空腔毛細(xì)波運(yùn)動以及氣泡截留的影響,主要得到以下結(jié)論.

    1)在恰當(dāng)?shù)淖赃m應(yīng)條件下,Gerris程序能夠在節(jié)約計(jì)算資源的同時較為準(zhǔn)確地預(yù)測水滴撞擊深水液池的運(yùn)動,數(shù)值模擬所得的界面變形、空腔成長、毛細(xì)波的傳播以及中心射流過程與實(shí)驗(yàn)結(jié)果符合良好.

    圖9 不同時間節(jié)點(diǎn)下水滴撞擊深水液池空腔運(yùn)動的渦量場和壓力場等值線圖 (a)Fr=567.1,Re=1595,We=121.8;(b)Fr=792.1,Re=1885,We=170.2Fig.9.Vorticity and pressure contours at selected time:(a)Fr=567.1,Re=1595,We=121.8;(b)Fr=792.1,Re=1885,We=170.2.

    2)液滴下落撞擊深水液池后,液面擴(kuò)張形成一個空腔,其后隨著毛細(xì)波運(yùn)動逐漸回縮.液滴完全聚合時,空腔形狀往往由U形向V形轉(zhuǎn)變;在液滴部分聚合生成細(xì)長中心射流并產(chǎn)生氣泡截留時,空腔初始形狀則近似一個半球形,其后在底部變形為圓柱形.

    3)液滴撞擊深水液池后,空腔深度先以較快速度增長,在到達(dá)最大空腔深度后快速回彈,至接近原自由液面后速度逐漸放緩.在忽略毛細(xì)波作用、空腔深度為h=D至h=hmax范圍內(nèi)的前提下,空腔深度隨時間的成長關(guān)系可由(Vit)/D=0.15(h/D)5/2來描述,但最終空腔底部的形成是由毛細(xì)波運(yùn)動決定的;空腔寬度的增長主要由毛細(xì)波運(yùn)動及其干擾驅(qū)動,最后在慣性力作用下線性增長.

    4)毛細(xì)波運(yùn)動可由壓力場中低壓區(qū)的位置示蹤,在撞擊速度較低,液滴完全聚結(jié)時(Fr=567.1,Re=1595,We=121.8),低壓區(qū)首先在空腔側(cè)壁與底部交界處產(chǎn)生,隨后向空腔底部傳播,在靠近液面以及空腔底部靠近中心區(qū)域各產(chǎn)生一個較大的渦環(huán);在發(fā)生部分聚結(jié)現(xiàn)象,產(chǎn)生細(xì)長射流時(Fr=792.1,Re=1885,We=170.2),渦環(huán)的生成被抑制,低壓區(qū)首先在波浪底部與側(cè)壁上交界處產(chǎn)生,空腔底部變?yōu)閳A柱狀后,毛細(xì)波在空腔坍塌前并沒有到達(dá)空腔底部中心,導(dǎo)致空腔側(cè)壁首先坍塌形成氣泡截留.

    感謝西安理工大學(xué)水力學(xué)研究所的李林博、楊泓、薛博升和楊琰青在論文完成過程中的幫助與討論.

    猜你喜歡
    液池空腔液面
    原油管道泄漏擴(kuò)散影響因素模擬分析
    基于邊光滑有限元法的二維復(fù)合彈性空腔聲振特性分析
    吸管“喝”水的秘密
    基于DCS自動控制循環(huán)水液面的改造
    電子測試(2018年6期)2018-05-09 07:31:47
    水面LNG液池?cái)U(kuò)展模型的分析與對比研究*
    LNG船泄漏事故液池?cái)U(kuò)展計(jì)算及不確定性分析
    中國航海(2017年2期)2017-10-30 08:43:09
    空腔參數(shù)對重力壩穩(wěn)定的影響分析
    前置污水去油池
    前置污水去油池
    激光系統(tǒng)對液面信息的探測與研究
    河南科技(2014年14期)2014-02-27 14:11:57
    十分钟在线观看高清视频www | 国产成人91sexporn| 噜噜噜噜噜久久久久久91| av不卡在线播放| 水蜜桃什么品种好| 韩国高清视频一区二区三区| 日本猛色少妇xxxxx猛交久久| 国产有黄有色有爽视频| 日本-黄色视频高清免费观看| 亚洲电影在线观看av| 日韩一本色道免费dvd| 亚洲成人av在线免费| 国产欧美日韩精品一区二区| 秋霞伦理黄片| kizo精华| 九草在线视频观看| 高清午夜精品一区二区三区| 边亲边吃奶的免费视频| 久久99蜜桃精品久久| 国产伦精品一区二区三区视频9| 一本色道久久久久久精品综合| 国产精品人妻久久久影院| 91久久精品国产一区二区三区| 人人妻人人添人人爽欧美一区卜 | 乱系列少妇在线播放| 亚洲av在线观看美女高潮| 岛国毛片在线播放| 草草在线视频免费看| 国精品久久久久久国模美| 亚洲av综合色区一区| 日韩成人伦理影院| 国产 一区精品| 九色成人免费人妻av| 亚洲av免费高清在线观看| 熟妇人妻不卡中文字幕| h日本视频在线播放| 男女下面进入的视频免费午夜| 久久青草综合色| 热re99久久精品国产66热6| 亚洲精华国产精华液的使用体验| 日本与韩国留学比较| 亚洲av电影在线观看一区二区三区| 熟女人妻精品中文字幕| 久久av网站| 日本与韩国留学比较| 久久人人爽av亚洲精品天堂 | 亚洲欧洲国产日韩| 一级毛片我不卡| 极品少妇高潮喷水抽搐| 男女免费视频国产| 一级片'在线观看视频| 美女cb高潮喷水在线观看| 久久久a久久爽久久v久久| 久久久久久久久久人人人人人人| 丰满人妻一区二区三区视频av| 午夜老司机福利剧场| 日日啪夜夜爽| 三级国产精品欧美在线观看| 91久久精品国产一区二区三区| 51国产日韩欧美| 18+在线观看网站| 国产精品国产三级专区第一集| 狂野欧美激情性xxxx在线观看| 99热6这里只有精品| 久久久精品免费免费高清| 男的添女的下面高潮视频| 少妇的逼好多水| 国产精品久久久久久精品电影小说 | 91午夜精品亚洲一区二区三区| 啦啦啦在线观看免费高清www| 2021少妇久久久久久久久久久| 久久久a久久爽久久v久久| 免费观看av网站的网址| av一本久久久久| 麻豆乱淫一区二区| 亚洲性久久影院| 成人黄色视频免费在线看| 亚洲精品久久午夜乱码| 国产乱人偷精品视频| 最近中文字幕2019免费版| 亚洲av男天堂| 在线免费十八禁| 成人午夜精彩视频在线观看| 观看av在线不卡| 一二三四中文在线观看免费高清| 久久精品国产a三级三级三级| av福利片在线观看| 日本猛色少妇xxxxx猛交久久| 狂野欧美激情性bbbbbb| 日韩,欧美,国产一区二区三区| 精品久久久久久久久亚洲| 91精品国产国语对白视频| 99热这里只有精品一区| 高清午夜精品一区二区三区| 欧美另类一区| 久久久久久伊人网av| 18禁裸乳无遮挡免费网站照片| 少妇精品久久久久久久| 亚洲精品日韩在线中文字幕| av在线播放精品| 亚洲精品乱码久久久v下载方式| av天堂中文字幕网| 国产精品免费大片| 嘟嘟电影网在线观看| 国产男人的电影天堂91| 丰满少妇做爰视频| 国产色婷婷99| 日本av免费视频播放| 一级片'在线观看视频| 国产男女超爽视频在线观看| 一区二区三区精品91| av一本久久久久| 五月天丁香电影| 日本欧美视频一区| 亚洲国产高清在线一区二区三| 在线观看一区二区三区| 黄色怎么调成土黄色| 最近的中文字幕免费完整| 在线观看一区二区三区激情| 少妇高潮的动态图| 久久精品久久久久久噜噜老黄| 精品久久久久久久久av| 性色avwww在线观看| 亚洲精品456在线播放app| 精品一品国产午夜福利视频| 赤兔流量卡办理| 免费av不卡在线播放| 亚洲av国产av综合av卡| 亚洲综合色惰| 青青草视频在线视频观看| 在现免费观看毛片| 在线亚洲精品国产二区图片欧美 | 七月丁香在线播放| 久久97久久精品| 午夜激情福利司机影院| av在线老鸭窝| 国产黄片美女视频| 欧美另类一区| av线在线观看网站| av网站免费在线观看视频| 日韩强制内射视频| 精品人妻一区二区三区麻豆| 午夜福利在线在线| 一区二区三区乱码不卡18| 免费人成在线观看视频色| 久久精品人妻少妇| 少妇的逼好多水| 又大又黄又爽视频免费| 国产精品无大码| 自拍欧美九色日韩亚洲蝌蚪91 | 九色成人免费人妻av| 国语对白做爰xxxⅹ性视频网站| 干丝袜人妻中文字幕| 精品视频人人做人人爽| 亚洲经典国产精华液单| 久久99蜜桃精品久久| 国产成人精品婷婷| 成人无遮挡网站| 亚洲精品中文字幕在线视频 | 欧美成人午夜免费资源| 人人妻人人爽人人添夜夜欢视频 | 久久久久国产网址| 国产精品国产三级专区第一集| 国产亚洲欧美精品永久| 一区二区三区乱码不卡18| 亚洲精品乱码久久久v下载方式| 久久鲁丝午夜福利片| 中国三级夫妇交换| 99热这里只有是精品在线观看| 伦理电影大哥的女人| 久久久国产一区二区| 青春草国产在线视频| 三级经典国产精品| 国产精品爽爽va在线观看网站| 纯流量卡能插随身wifi吗| 男女免费视频国产| 啦啦啦中文免费视频观看日本| 精品亚洲成国产av| 国产真实伦视频高清在线观看| 大片免费播放器 马上看| 成人一区二区视频在线观看| 日本wwww免费看| 99久国产av精品国产电影| 蜜桃在线观看..| 26uuu在线亚洲综合色| 男人和女人高潮做爰伦理| 尾随美女入室| 一级毛片 在线播放| 国产中年淑女户外野战色| 亚洲av二区三区四区| 国产精品欧美亚洲77777| 久久久久久久国产电影| 国产黄色免费在线视频| 日本欧美国产在线视频| 亚洲av二区三区四区| 一级a做视频免费观看| 嘟嘟电影网在线观看| 国产免费福利视频在线观看| 国产精品爽爽va在线观看网站| 秋霞在线观看毛片| 欧美成人一区二区免费高清观看| 99久久综合免费| 久久精品久久久久久久性| 精品人妻一区二区三区麻豆| 女性生殖器流出的白浆| 亚洲欧美精品专区久久| 欧美国产精品一级二级三级 | 精品久久久精品久久久| 亚洲av在线观看美女高潮| 久久国产精品大桥未久av | 一本—道久久a久久精品蜜桃钙片| 天天躁日日操中文字幕| 欧美成人精品欧美一级黄| 伦精品一区二区三区| 国产精品成人在线| 亚洲成人中文字幕在线播放| 亚洲精品国产成人久久av| 秋霞在线观看毛片| 欧美激情国产日韩精品一区| 国产精品免费大片| 日本免费在线观看一区| 尾随美女入室| 国产成人一区二区在线| 黑丝袜美女国产一区| 久久久久久伊人网av| 少妇裸体淫交视频免费看高清| 王馨瑶露胸无遮挡在线观看| 亚洲四区av| 久久精品国产亚洲av天美| 亚洲欧美成人精品一区二区| 日本黄大片高清| 国产亚洲欧美精品永久| 国产片特级美女逼逼视频| 极品少妇高潮喷水抽搐| 国产精品久久久久久久电影| 国产精品一及| 18禁在线无遮挡免费观看视频| 精品99又大又爽又粗少妇毛片| 女人久久www免费人成看片| 国产 一区精品| 一二三四中文在线观看免费高清| 中文字幕免费在线视频6| 久久精品国产亚洲av涩爱| 蜜臀久久99精品久久宅男| 青春草亚洲视频在线观看| 久久久久网色| 精品亚洲成国产av| 成年美女黄网站色视频大全免费 | 美女cb高潮喷水在线观看| 国产精品一二三区在线看| 最近最新中文字幕免费大全7| 亚洲婷婷狠狠爱综合网| 18禁动态无遮挡网站| 成人午夜精彩视频在线观看| 天天躁夜夜躁狠狠久久av| 成人特级av手机在线观看| 国产亚洲一区二区精品| 亚洲最大成人中文| 在线观看一区二区三区| 大片免费播放器 马上看| 国产爽快片一区二区三区| 少妇高潮的动态图| 99久久精品热视频| 国语对白做爰xxxⅹ性视频网站| 最近最新中文字幕大全电影3| 免费观看在线日韩| 免费av中文字幕在线| 日本黄色片子视频| 女人久久www免费人成看片| 蜜臀久久99精品久久宅男| 亚洲av二区三区四区| 九九久久精品国产亚洲av麻豆| 亚洲欧美日韩卡通动漫| 干丝袜人妻中文字幕| 日韩一区二区视频免费看| 纵有疾风起免费观看全集完整版| 多毛熟女@视频| 日韩一区二区三区影片| 在线播放无遮挡| 在线观看免费日韩欧美大片 | 久久精品国产亚洲av涩爱| 国内揄拍国产精品人妻在线| av国产久精品久网站免费入址| 久久国产精品男人的天堂亚洲 | 噜噜噜噜噜久久久久久91| 在线观看一区二区三区激情| 18禁裸乳无遮挡动漫免费视频| 午夜精品国产一区二区电影| 亚洲天堂av无毛| 亚洲精品成人av观看孕妇| 亚洲电影在线观看av| 国产精品久久久久久av不卡| 黄色配什么色好看| 女人十人毛片免费观看3o分钟| 欧美日本视频| 国产精品一区二区在线不卡| 国产精品人妻久久久久久| 少妇的逼好多水| 国产高清三级在线| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久久久久免| 久久国产乱子免费精品| 嘟嘟电影网在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲国产色片| 边亲边吃奶的免费视频| 下体分泌物呈黄色| 久久热精品热| 97精品久久久久久久久久精品| 大话2 男鬼变身卡| 最近手机中文字幕大全| 久久久久久久久久人人人人人人| 一二三四中文在线观看免费高清| 女人久久www免费人成看片| 国产精品女同一区二区软件| a级一级毛片免费在线观看| 在线观看三级黄色| 国产色婷婷99| 日韩av免费高清视频| 亚洲av成人精品一区久久| www.色视频.com| 在线亚洲精品国产二区图片欧美 | 91午夜精品亚洲一区二区三区| 在线 av 中文字幕| 91久久精品电影网| 最近2019中文字幕mv第一页| 午夜日本视频在线| 秋霞在线观看毛片| 一区二区三区精品91| 亚洲欧美日韩另类电影网站 | 国产精品久久久久成人av| 久久久久久久久久人人人人人人| 97热精品久久久久久| 国产黄片美女视频| 丰满少妇做爰视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲高清免费不卡视频| 少妇高潮的动态图| 18禁在线无遮挡免费观看视频| 国产精品.久久久| 亚洲色图综合在线观看| 女人久久www免费人成看片| 国产亚洲午夜精品一区二区久久| 欧美97在线视频| 人体艺术视频欧美日本| 日日摸夜夜添夜夜爱| 国产精品99久久久久久久久| 日本爱情动作片www.在线观看| 日韩成人伦理影院| 天美传媒精品一区二区| 久久久久国产网址| 嫩草影院入口| 黄色日韩在线| 日韩免费高清中文字幕av| .国产精品久久| 嫩草影院入口| 在线天堂最新版资源| 欧美日韩国产mv在线观看视频 | av女优亚洲男人天堂| 国产免费一级a男人的天堂| 国产亚洲精品久久久com| 免费观看av网站的网址| 久久久久网色| 亚洲精品视频女| 午夜激情福利司机影院| 国产黄频视频在线观看| 亚洲精品中文字幕在线视频 | 97超视频在线观看视频| 亚洲三级黄色毛片| 精品久久国产蜜桃| 欧美老熟妇乱子伦牲交| 嫩草影院新地址| 亚洲在久久综合| 免费观看性生交大片5| 久久精品国产亚洲av涩爱| 日本av手机在线免费观看| 国产毛片在线视频| 女人十人毛片免费观看3o分钟| 九草在线视频观看| 欧美日韩国产mv在线观看视频 | 久久人人爽人人片av| 国产精品国产三级专区第一集| 国产乱来视频区| 国产黄频视频在线观看| 欧美日韩亚洲高清精品| 日韩精品有码人妻一区| 新久久久久国产一级毛片| 日韩大片免费观看网站| 男女无遮挡免费网站观看| 亚洲内射少妇av| 日韩一区二区视频免费看| 成人特级av手机在线观看| 亚洲精华国产精华液的使用体验| 久久久久久久精品精品| 人人妻人人澡人人爽人人夜夜| 亚洲婷婷狠狠爱综合网| 欧美最新免费一区二区三区| 18禁在线播放成人免费| 99热这里只有是精品50| 久久久久久人妻| 久久久久久九九精品二区国产| 卡戴珊不雅视频在线播放| 亚洲av成人精品一区久久| 精华霜和精华液先用哪个| 中文字幕亚洲精品专区| 国产精品一及| 久久婷婷青草| 久久综合国产亚洲精品| 另类亚洲欧美激情| 爱豆传媒免费全集在线观看| 午夜老司机福利剧场| 丝袜脚勾引网站| 中文资源天堂在线| 五月玫瑰六月丁香| 国产亚洲一区二区精品| 久久久久久九九精品二区国产| 99久久综合免费| 在线精品无人区一区二区三 | 久久精品夜色国产| 久久女婷五月综合色啪小说| 一级毛片电影观看| 国产大屁股一区二区在线视频| 国产精品福利在线免费观看| 日本一二三区视频观看| 99re6热这里在线精品视频| 久久人人爽人人爽人人片va| 欧美3d第一页| 两个人的视频大全免费| 王馨瑶露胸无遮挡在线观看| 国产精品99久久久久久久久| 有码 亚洲区| 乱码一卡2卡4卡精品| 少妇人妻一区二区三区视频| 国产精品久久久久久精品电影小说 | 亚洲欧美精品专区久久| 人妻少妇偷人精品九色| 伦理电影大哥的女人| 美女xxoo啪啪120秒动态图| 亚洲人成网站高清观看| 99热国产这里只有精品6| 男人爽女人下面视频在线观看| 大片免费播放器 马上看| 亚洲人成网站在线观看播放| 国产欧美另类精品又又久久亚洲欧美| 午夜福利在线观看免费完整高清在| 老女人水多毛片| 26uuu在线亚洲综合色| 成年免费大片在线观看| 亚洲激情五月婷婷啪啪| 国产精品久久久久久精品古装| 一本一本综合久久| 又粗又硬又长又爽又黄的视频| 中文字幕制服av| 午夜福利在线观看免费完整高清在| videossex国产| .国产精品久久| 久久久久久久精品精品| 免费看不卡的av| 国产精品三级大全| tube8黄色片| 亚洲国产精品专区欧美| 女的被弄到高潮叫床怎么办| 国产黄频视频在线观看| 久久久久国产精品人妻一区二区| 91久久精品电影网| 狂野欧美白嫩少妇大欣赏| 肉色欧美久久久久久久蜜桃| 免费观看性生交大片5| 免费看不卡的av| 一级a做视频免费观看| 日本欧美视频一区| 中国三级夫妇交换| 老师上课跳d突然被开到最大视频| 最黄视频免费看| 免费黄频网站在线观看国产| 噜噜噜噜噜久久久久久91| 欧美激情极品国产一区二区三区 | 伦精品一区二区三区| 嫩草影院入口| 国产爱豆传媒在线观看| 妹子高潮喷水视频| 十分钟在线观看高清视频www | 麻豆成人av视频| 大片免费播放器 马上看| 成人毛片60女人毛片免费| 国产69精品久久久久777片| 久久99热这里只有精品18| 在线观看一区二区三区激情| 欧美一区二区亚洲| 美女脱内裤让男人舔精品视频| 亚洲欧美一区二区三区国产| 欧美bdsm另类| 国产一区有黄有色的免费视频| 搡老乐熟女国产| 网址你懂的国产日韩在线| av女优亚洲男人天堂| 欧美日韩视频精品一区| 国产亚洲av片在线观看秒播厂| 亚洲欧洲日产国产| 国产精品.久久久| 亚洲国产色片| 欧美xxxx黑人xx丫x性爽| av在线蜜桃| 欧美3d第一页| 国产欧美另类精品又又久久亚洲欧美| 亚洲性久久影院| 午夜老司机福利剧场| 亚洲精品视频女| 人妻系列 视频| 91aial.com中文字幕在线观看| 国产亚洲午夜精品一区二区久久| 最近的中文字幕免费完整| 人妻夜夜爽99麻豆av| 男人爽女人下面视频在线观看| 中文天堂在线官网| 少妇 在线观看| 在线观看免费日韩欧美大片 | 黄色配什么色好看| av国产精品久久久久影院| 亚洲,欧美,日韩| 综合色丁香网| 看免费成人av毛片| 亚洲,一卡二卡三卡| 久久久久久久久久久免费av| 国产精品久久久久久av不卡| 黄色日韩在线| 亚洲国产最新在线播放| 亚洲国产高清在线一区二区三| av不卡在线播放| 亚洲av欧美aⅴ国产| 一级av片app| 日韩中文字幕视频在线看片 | 午夜福利影视在线免费观看| 自拍偷自拍亚洲精品老妇| 精品亚洲成国产av| 少妇裸体淫交视频免费看高清| 肉色欧美久久久久久久蜜桃| 国产精品伦人一区二区| 亚洲av福利一区| 久久久久国产网址| 黄色视频在线播放观看不卡| 亚洲欧美成人综合另类久久久| 80岁老熟妇乱子伦牲交| 久久青草综合色| 日韩欧美 国产精品| 成年av动漫网址| 97在线视频观看| 欧美最新免费一区二区三区| 人妻系列 视频| 日韩大片免费观看网站| 99国产精品免费福利视频| av不卡在线播放| 观看av在线不卡| 亚洲丝袜综合中文字幕| 99久久综合免费| 2022亚洲国产成人精品| 老女人水多毛片| 天美传媒精品一区二区| 亚洲精品aⅴ在线观看| 三级国产精品片| 婷婷色综合大香蕉| 国产乱人偷精品视频| 色综合色国产| 国产男女内射视频| 80岁老熟妇乱子伦牲交| 九色成人免费人妻av| 久久久精品免费免费高清| 多毛熟女@视频| 久久av网站| 国产成人免费无遮挡视频| 大香蕉久久网| 婷婷色综合www| av网站免费在线观看视频| 国产精品一二三区在线看| 精品人妻视频免费看| 国产高清三级在线| av在线播放精品| 久久精品人妻少妇| 精品亚洲成国产av| 免费观看性生交大片5| 精品人妻偷拍中文字幕| 中文在线观看免费www的网站| 午夜免费鲁丝| 日韩人妻高清精品专区| 国产大屁股一区二区在线视频| 精品少妇久久久久久888优播| 国产爱豆传媒在线观看| 中文字幕人妻熟人妻熟丝袜美| 校园人妻丝袜中文字幕| 精品人妻偷拍中文字幕| 国产亚洲av片在线观看秒播厂| 久久99热这里只频精品6学生| 久久久久久九九精品二区国产| 97在线人人人人妻| 观看美女的网站| 精品人妻偷拍中文字幕| 最近最新中文字幕大全电影3| 欧美少妇被猛烈插入视频| 欧美+日韩+精品| 精品久久久精品久久久| 日日摸夜夜添夜夜添av毛片| 26uuu在线亚洲综合色| 乱码一卡2卡4卡精品| 亚洲内射少妇av| 蜜桃在线观看..| 国产高清有码在线观看视频| 国产精品一区二区在线观看99| 老司机影院成人| 亚洲欧美精品专区久久| 久久久久久久久久久丰满| 一级片'在线观看视频| 国产毛片在线视频| 国产免费一区二区三区四区乱码| 成年美女黄网站色视频大全免费 | 久久久久人妻精品一区果冻|