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

    農(nóng)藥霧滴空間運(yùn)行中的變形特征分析

    2020-10-20 02:12:02馬學(xué)虎薛士東孫桐奚溪宋小沫趙峻逸蘭忠郝婷婷
    化工進(jìn)展 2020年10期
    關(guān)鍵詞:曳力橢球表面張力

    馬學(xué)虎,薛士東,孫桐,奚溪,宋小沫,趙峻逸,蘭忠,郝婷婷

    (遼寧省化工資源清潔利用重點(diǎn)實(shí)驗(yàn)室,大連理工大學(xué)化學(xué)工程研究所,遼寧大連116024)

    化學(xué)農(nóng)藥在重大病蟲草害防控、保證國家糧食安全方面發(fā)揮著不可替代的作用,合理使用農(nóng)藥有助于提高農(nóng)業(yè)生產(chǎn)力和保障農(nóng)作物質(zhì)量,使用農(nóng)藥可減少全球高達(dá)45%的糧食損失[1]。據(jù)國家統(tǒng)計(jì)局?jǐn)?shù)據(jù)顯示,2014 年中國農(nóng)藥的使用量為180.69×104t,但是由于噴霧技術(shù)和噴霧管理的限制,農(nóng)藥的利用率僅為20%~30%[2]。

    農(nóng)藥的使用過程一般分為以下幾個(gè)步驟[3-4]:①藥液通過噴頭完成霧化過程;②霧滴在空氣中運(yùn)動(dòng)擴(kuò)散,沉積到靶標(biāo)冠層或飄失到靶標(biāo)外;③霧滴在冠層表面撞擊、鋪展;④霧滴被冠層表面吸收,傳導(dǎo)到為害部位并發(fā)揮藥效。根據(jù)ISO 22856 標(biāo)準(zhǔn)所定義,霧滴飄移是指“在施藥過程中,農(nóng)藥霧滴被氣流攜帶脫離靶標(biāo)冠層、流失到地面或者懸浮在大氣中的行為”,被認(rèn)為是導(dǎo)致農(nóng)藥污染的“罪魁禍?zhǔn)住?,?duì)周圍河流中的水生動(dòng)植物、大氣環(huán)境、人類身體健康造成不可扭轉(zhuǎn)的危害[5-10]。霧滴飄移是諸多因素共同作用的結(jié)果,包括藥液的理化性質(zhì)、霧滴的粒徑及分布、作物結(jié)構(gòu)(排與排的間距、冠層結(jié)構(gòu)、冠層穿透性)、操作條件(噴霧流量、操作壓力、農(nóng)機(jī)行駛速度)、氣候條件(風(fēng)速、風(fēng)向、溫度、濕度)等[11-18]。

    研究霧滴飄移的手段主要包括風(fēng)洞試驗(yàn)、田間試驗(yàn)、數(shù)值模擬等。在風(fēng)速可調(diào)、溫濕度可控的低速風(fēng)洞中,研究噴頭類型、操作壓力、助劑種類等對(duì)霧滴群飄移比例的影響[19-21]。在田間試驗(yàn)中,通過檢測靶標(biāo)冠層內(nèi)示蹤劑的含量來評(píng)價(jià)霧滴群的飄移沉積比例[22-24]。在數(shù)值模擬中,采用歐拉-拉格朗日離散相模型追蹤霧滴群的運(yùn)動(dòng)軌跡,從而獲得特定條件下霧滴群的飄移率[25-27]。

    以往的數(shù)值模擬研究中多將農(nóng)藥霧滴視為不會(huì)發(fā)生形變的剛性顆粒,忽略了霧滴的形變過程。在霧滴撞擊靶標(biāo)和界面蒸發(fā)的研究中也多以完美球狀作為霧滴的初始狀態(tài)[28-30]。然而大量研究表明,霧滴形狀對(duì)霧滴的聚并與破碎、靶標(biāo)撞擊、所受的曳力大小等都有著重要作用,會(huì)直接影響霧滴的飄移損失與界面流失等。QUAN 等[31]利用耦合了移動(dòng)網(wǎng)格相界面追蹤的三維有限體積交錯(cuò)網(wǎng)格法模擬了液滴在加速氣流中的運(yùn)動(dòng)情況,分析了液滴的形變狀態(tài)及附近氣流的速度場變化。結(jié)果表明,變形過程中液滴的非穩(wěn)態(tài)曳力總是比同等雷諾數(shù)下液滴的穩(wěn)態(tài)曳力大,同時(shí)降低表面張力和黏度會(huì)增大液滴的形變量,導(dǎo)致曳力增加,而密度比對(duì)液滴動(dòng)力學(xué)影響不顯著。QU 等[32]用ANSYS Fluent 軟件對(duì)非對(duì)稱流中的減速液滴進(jìn)行模擬,結(jié)果表明隨著We 的增大,液滴的形變量及曳力會(huì)急劇增大,振蕩周期延長;隨著Oh 的增大,形變量與曳力略有增大而振蕩周期幾乎不變。

    目前針對(duì)單個(gè)農(nóng)藥霧滴在風(fēng)場中的運(yùn)動(dòng)研究鮮有報(bào)道,尤其是對(duì)微米或者亞毫米級(jí)別農(nóng)藥霧滴的空間變形特征更是缺乏研究。單霧滴是構(gòu)成霧滴群的基礎(chǔ)元素,明晰單霧滴的運(yùn)行飄移規(guī)律,對(duì)于深入研究整個(gè)霧滴群的運(yùn)動(dòng)規(guī)律具有重要的指導(dǎo)意義。因此本文以農(nóng)藥噴霧中的單霧滴為研究對(duì)象,采用數(shù)值模擬的手段研究了霧滴在橫向風(fēng)場中的空間形變特性,揭示了不同因素作用下單霧滴的形態(tài)演化規(guī)律,并建立了多因素協(xié)同下霧滴形變量與形變周期的預(yù)測模型。

    1 數(shù)值模擬計(jì)算模型與方法

    1.1 幾何模型

    基于COMSOL Multiphysics 5.2建立數(shù)值模擬過程,選擇帶相初始化的層流-水平集物理場,對(duì)氣液相界面進(jìn)行追蹤,獲得霧滴的空間形變特性。如圖1 所示,針對(duì)直徑為100μm 的霧滴,選取1200μm×300μm二維矩形計(jì)算域,霧滴的初始位置坐標(biāo)為(100μm,1100μm)。計(jì)算域的尺寸隨著霧滴的大小進(jìn)行相應(yīng)調(diào)整,以保證計(jì)算域的合理利用并盡可能降低計(jì)算成本。

    圖1 幾何模型

    1.2 控制方程及邊界

    對(duì)含有分界面的不可壓縮、不混溶的氣液兩相流動(dòng)系統(tǒng)進(jìn)行模擬時(shí),關(guān)鍵問題在于如何將界面離散化并且準(zhǔn)確地追蹤界面的演化過程。Level Set方法通過設(shè)置一個(gè)距離函數(shù)隱式追蹤自由面,并且將界面條件融入控制方程中,可在整個(gè)域內(nèi)用統(tǒng)一的差分格式進(jìn)行求解。與移動(dòng)網(wǎng)格方法和相場方法相比,可以避免顯式地追蹤界面,不需要進(jìn)行復(fù)雜的界面重構(gòu),可方便地追蹤界面的拓?fù)浣Y(jié)構(gòu)改變,兩相控制方程組如式(1)、式(2)所示。

    連續(xù)性方程

    Navier-Stokes方程

    式中,F(xiàn)σ為單位體積流體所受的表面張力;Fd為單位體積流體所受的曳力。由于氣相密度遠(yuǎn)小于液相密度,因此忽略霧滴所受浮力的作用。

    霧滴在橫向風(fēng)場中所受的曳力大小如式(3)所示。

    式中,A為霧滴在風(fēng)向上的投影面積,在二維計(jì)算中降維為投影長度;ud為霧滴在空間某一位置處的速度矢量;ua為相應(yīng)位置處空氣的速度矢量;Cd為曳力系數(shù),計(jì)算依據(jù)選取文獻(xiàn)[33]中提供的農(nóng)藥噴灑過程中霧滴雷諾數(shù)在0.25<Re<104范圍內(nèi)的經(jīng)驗(yàn)公式,如式(4)所示。所受曳力方向與霧滴的相對(duì)運(yùn)動(dòng)方向相反。

    表面張力可以在整個(gè)計(jì)算域中連續(xù)求解,表達(dá)為一個(gè)光滑函數(shù)的形式,見式(5)。

    其中,κ為相界面曲率,其值見式(6)。

    n為氣液相界面法向向量,其值見式(7)。

    φ為經(jīng)Heaviside函數(shù)平滑后的水平集函數(shù),代表網(wǎng)格內(nèi)空氣相的體積分?jǐn)?shù),其在界面一側(cè)的值恒為1,另一側(cè)恒為0,在相界面很小的寬度內(nèi)由1過渡到0。分別求解各相內(nèi)的流體物性,在相界面處借助水平集函數(shù)對(duì)物性參數(shù)進(jìn)行定義,使物性在界面處得以連續(xù)。計(jì)算過程見式(8)、式(9)。

    δ表示規(guī)整后的delta函數(shù),其值見式(10)。

    δ為一維Dirac函數(shù),通過δ可將邊界上的積分改寫為體積上的積分。δ只在相界面附近有值,從而使表面張力的作用限定在界面處。

    1.3 邊界條件

    二維矩形計(jì)算域的左側(cè)設(shè)置為速度入口,如式(11)。

    右側(cè)邊界為壓力出口,與大氣連通,如式(12)。

    上下邊界定義為周期性流動(dòng)條件,以減小計(jì)算域的面積,縮短計(jì)算時(shí)間,減少占用內(nèi)存。域內(nèi)空氣初速度大小等于其入口速度ud,0,霧滴在初始化位置具有豎直向下的初速度,如式(13)。

    本文中所有算例霧滴的初速度設(shè)定為10m/s[34]。

    1.4 網(wǎng)格劃分

    為提高網(wǎng)格規(guī)整度,防止液滴初始化導(dǎo)致的網(wǎng)格不規(guī)則分布,使用網(wǎng)格控制邊對(duì)計(jì)算域分區(qū),在各子域內(nèi)劃分網(wǎng)格,如圖1中虛線所示。霧滴初始化位置四周使用自由三角形網(wǎng)格,余下區(qū)域應(yīng)用較規(guī)整的四邊形單元map網(wǎng)格,以提高計(jì)算精度并減少所需計(jì)算資源。本模擬中通過設(shè)置網(wǎng)格最大單元尺寸和最小單元尺寸,在COMSOL 中自動(dòng)進(jìn)行網(wǎng)格的生成,計(jì)算前進(jìn)行網(wǎng)格無關(guān)性檢驗(yàn),以優(yōu)化網(wǎng)格的劃分。

    1.5 后處理

    通過霧滴在空間運(yùn)行中的軌跡變化來判斷網(wǎng)格選擇的合理性,需要得到不同時(shí)刻霧滴的位置。如圖2所示,將計(jì)算得到的體積分?jǐn)?shù)圖像用MATLAB進(jìn)行編程處理,將圖像轉(zhuǎn)化為二值圖像,提取相界面的輪廓。通過計(jì)算輪廓各點(diǎn)坐標(biāo)的均值可獲得各個(gè)時(shí)刻霧滴的質(zhì)心坐標(biāo),從而連成霧滴的運(yùn)動(dòng)軌跡。

    霧滴在空間運(yùn)行中會(huì)由球狀變成橢球狀,為了描述霧滴在橫向風(fēng)場中下落時(shí)的形態(tài)演化規(guī)律,定義霧滴的形變比E=b/a,其中a 代表橢圓長軸的長度,b代表橢圓短軸的長度。定義霧滴的形變量為e=1-E,當(dāng)霧滴呈現(xiàn)球狀時(shí)e=0,橢球狀時(shí)0<e<1。

    圖2 后處理過程

    2 結(jié)果與討論

    2.1 網(wǎng)格無關(guān)性驗(yàn)證

    以直徑為100μm 的霧滴為例,設(shè)定最小及最大單元尺寸自動(dòng)劃分網(wǎng)格。霧滴下落是一個(gè)瞬態(tài)過程,以霧滴在0.1ms內(nèi)的運(yùn)動(dòng)軌跡及速度變化作為評(píng)價(jià)指標(biāo)來驗(yàn)證網(wǎng)格無關(guān)性,結(jié)果如圖3所示。

    當(dāng)網(wǎng)格數(shù)達(dá)到45569時(shí),再提高網(wǎng)格精度對(duì)霧滴的軌跡及速度影響不大,可認(rèn)為該網(wǎng)格精度足夠高,具有較好的計(jì)算性價(jià)比。此時(shí)網(wǎng)格最大單元尺寸為3μm,最小為0.01μm。所有算例的域內(nèi)最大速度基本一致,對(duì)不同尺寸的霧滴選取相同的網(wǎng)格劃分精度,保證了所有霧滴都具有相同的界面厚度。

    2.2 霧滴空間形變過程

    具有向下初速度的霧滴在左側(cè)來流橫向風(fēng)中運(yùn)動(dòng)的形態(tài)演化過程,如圖4所示。霧滴在一個(gè)振蕩周期Δtd內(nèi)經(jīng)歷了“最初的球→橢球→球→橢球→球”的形變過程,橢球長軸方向與豎直方向有一定的夾角。霧滴由起始時(shí)刻的球形變?yōu)闄E球狀再回到球形的這一過程簡稱為“初次變形”,所經(jīng)歷的時(shí)間記為Δt1;繼而再由球形變?yōu)闄E球到振蕩周期結(jié)束這一過程簡稱為“二次變形”,所經(jīng)歷的時(shí)間記為Δt2,Δtd=Δt1+Δt2。霧滴受到左側(cè)來流橫向風(fēng)的作用,同時(shí)自身具有向下的初動(dòng)能,在霧滴下緣產(chǎn)生一個(gè)低壓區(qū),從左側(cè)來流的橫向風(fēng)更容易進(jìn)入霧滴下側(cè),氣流對(duì)霧滴產(chǎn)生剪切作用,使霧滴形變?yōu)殚L軸偏向上風(fēng)向的橢球形。隨著霧滴不斷下降,霧滴的形狀逐漸變?yōu)榫哂休^大曲率的橢球形,導(dǎo)致霧滴的表面張力增大,逐漸增大的表面張力與曳力相平衡使得形變程度慢慢減弱,在一段時(shí)間后幾乎又恢復(fù)到球狀,但由于霧滴的形變?nèi)栽诶^續(xù),二次形變中霧滴的拉伸與壓縮方向與第一次相反,曳力的阻礙作用導(dǎo)致形變程度相對(duì)于第一次發(fā)生了衰減。隨著霧滴繼續(xù)下落,再一次恢復(fù)到近似球狀,至此一個(gè)振蕩周期結(jié)束。

    圖3 網(wǎng)格無關(guān)性驗(yàn)證(ua=3m/s;μd=1.0mPa·s;σ=72.7mN/m)

    圖4 150μm的霧滴在橫向風(fēng)中運(yùn)動(dòng)的形變過程

    霧滴在橫向風(fēng)場中運(yùn)動(dòng)的形變過程是多種因素共同作用下的結(jié)果。如圖5所示,為不同條件下霧滴達(dá)到最大形變量時(shí)的形態(tài)特征。下面依次討論霧滴直徑、橫向風(fēng)速、藥液的理化性質(zhì)等對(duì)霧滴形變特性的影響。

    圖5 不同條件下霧滴達(dá)到最大形變量的計(jì)算結(jié)果圖

    2.2.1 霧滴直徑對(duì)霧滴形變的影響

    根據(jù)美國農(nóng)業(yè)工程師學(xué)會(huì)(ASAE)對(duì)農(nóng)藥霧滴尺寸的分類標(biāo)準(zhǔn),取直徑為100μm、150μm、200μm、250μm、300μm 的細(xì)霧滴、中等霧滴、粗霧滴,將它們置于速度為3m/s的橫向風(fēng)場中,液相黏度為1.0mPa·s,表面張力為72.7mN/m??疾祆F滴直徑對(duì)空間形變特性的影響,結(jié)果如圖6所示。

    圖6 不同直徑霧滴在橫向風(fēng)場中形變量e隨時(shí)間的變化

    結(jié)果表明,霧滴的振蕩周期Δtd和形變量e 均隨著直徑的增大而增大。當(dāng)直徑增大時(shí),霧滴的重力和曳力均會(huì)增加,在與黏性力及表面張力的競爭中取得優(yōu)勢(shì),保持原有運(yùn)動(dòng)狀態(tài)的趨勢(shì)增強(qiáng),霧滴形變的周期會(huì)延長,同時(shí)霧滴自身重力的提高會(huì)增大形變程度。在一個(gè)振蕩周期Δtd內(nèi),形變量的第二個(gè)峰值小于第一個(gè)峰值emax1>emax2,同時(shí)有Δt1>Δt2,表明振蕩形變發(fā)生了衰減。一方面是由于霧滴在下落過程中速度減慢,在一定程度上抑制了形變過程;另一方面是由于氣流方向向右,空氣對(duì)霧滴的拖拽作用總是沿著左上右下這一方向,使得霧滴在此方向上被拉長得更多。另外可以發(fā)現(xiàn),霧滴從橢球變?yōu)榍虻臅r(shí)間略短于從球變?yōu)闄E球的時(shí)間,這是因?yàn)閷?duì)于橢球而言,由于它發(fā)生了形變,在風(fēng)的投影方向長度拉長,空氣給予橢球霧滴的曳力要大于球形霧滴,所以在霧滴呈現(xiàn)橢球形時(shí)速度衰減更快,形變也就更快。

    2.2.2 橫向風(fēng)速對(duì)霧滴形變的影響

    取直徑為100μm、液相黏度為1.0mPa·s、表面張力為72.7mN/m 的霧滴,分別將其置于橫向風(fēng)速為1m/s、3m/s、5m/s、7m/s的穩(wěn)定氣流場中,考察橫向風(fēng)速對(duì)霧滴形變特性的影響。

    如圖7 所示,為霧滴形變量e 隨時(shí)間的變化情況。隨著橫向風(fēng)速的增加,霧滴的形變程度也呈現(xiàn)增大的趨勢(shì)。霧滴受到空氣的曳力作用而發(fā)生形變,曳力大小與氣液相對(duì)速度呈正相關(guān),風(fēng)速越大意味著空氣拖拽霧滴的能力越強(qiáng),形變?cè)綇?qiáng)烈。尤其當(dāng)風(fēng)速達(dá)到7m/s 時(shí),形變量有明顯提升,第一次形變幅值emax1達(dá)到了風(fēng)速為5m/s時(shí)的1.6倍左右,第二次形變幅值emax2達(dá)到了5m/s 時(shí)的1.2 倍左右。此外,分析霧滴的振蕩周期可以發(fā)現(xiàn),隨著風(fēng)速的增大,Δt1、Δt2和Δtd幾乎不發(fā)生改變,即4 種情況下霧滴的振蕩周期都為1.32×10-4s,達(dá)到第一個(gè)形變峰值的時(shí)間都在0.048ms左右。結(jié)果表明,霧滴的振蕩周期受橫向風(fēng)速的影響不明顯,振蕩周期為一個(gè)本征量,由霧滴自身特性決定。

    圖7 直徑為100μm的霧滴在不同橫向風(fēng)場速度下形變量e隨時(shí)間的變化

    2.2.3 液相黏度對(duì)霧滴形變的影響

    液體的黏度反映出流體內(nèi)部的黏性效應(yīng),黏度越大的液體表明其內(nèi)部的凝聚力越大,因此黏度的大小勢(shì)必影響著霧滴的形變程度??疾祆F滴初始直徑為150μm,橫向風(fēng)速3m/s,表面張力為72.7mN/m,液相黏度分別為0.3mPa·s、0.5mPa·s、1.0mPa·s、2.0mPa·s時(shí),霧滴的空間形變特性。

    如圖8所示,隨著液相黏度增大,霧滴的形變量減小,但振蕩周期Δtd沒有明顯的變化,說明霧滴的振蕩周期與液相黏度幾乎無關(guān)。當(dāng)液相黏度為0.3mPa?s 和0.5mPa?s 時(shí),Δt1與Δt2幾乎相等,各占一個(gè)完整周期的一半,同時(shí)霧滴的振蕩幅度也沒有衰減,兩次形變量最大值幾近相同。然而當(dāng)液相黏度為1.0mPa?s 和2.0mPa?s 時(shí),emax1>emax2,Δt1>Δt2,形變程度發(fā)生了衰減。分析原因認(rèn)為,當(dāng)液相黏度較大時(shí),霧滴的形變主要由兩相間的剪切作用決定,故前半周期比后半周期長;而當(dāng)液相黏度較小時(shí),內(nèi)部凝聚力小,為影響形變的主導(dǎo)作用,形變的衰減現(xiàn)象不再顯著。此外,黏度較小的霧滴在第一個(gè)振蕩周期結(jié)束時(shí)難以恢復(fù)到近球形,黏度越小這一趨勢(shì)越明顯,這是由于霧滴自身的黏度較小,內(nèi)聚力較弱,在運(yùn)動(dòng)過程中難以與曳力相抗衡,很難恢復(fù)到最初的形狀,只能維持一定的形變量在空氣中繼續(xù)運(yùn)動(dòng)。

    圖8 不同黏度的150μm霧滴在橫向風(fēng)場中形變量e隨時(shí)間的變化

    2.2.4 表面張力對(duì)霧滴形變的影響

    根據(jù)文獻(xiàn)[16,35-36]報(bào)道,按照一定的配比添加不同種類的助劑可以調(diào)控霧滴的表面張力在20~73mN/m 之間,在田間實(shí)際操作的范圍內(nèi)。考察霧滴直徑為150μm,橫向風(fēng)速為3m/s,液相黏度為1.0mPa·s,表面張力分別為28.5mN/m、37.4mN/m、56.9mN/m、72.7mN/m時(shí)空間形變特性。

    根據(jù)能量最低原理,在氣液相界面處的液相表面分子會(huì)自發(fā)趨向于液相內(nèi)部而非表面,這使得表面張力對(duì)塑造霧滴形態(tài)起著十分關(guān)鍵的作用[37]。如圖9所示,不同表面張力下霧滴的形變規(guī)律與前述算例相同,形變幅度逐漸衰減(emax1>emax2),形變時(shí)間逐漸縮短(Δt1>Δt2)。當(dāng)液相表面張力較小時(shí),在與曳力、重力的競爭中處于弱勢(shì)地位,霧滴恢復(fù)到表面能較低的球形狀態(tài)更加困難,從而延長了霧滴振蕩形變周期,增加了形變幅度。無論是前半個(gè)周期結(jié)束時(shí)還是整個(gè)振蕩周期結(jié)束時(shí),表面張力越小的霧滴越難以恢復(fù)到初始的球形,如表面張力為28.5mN/m的霧滴在0.392ms時(shí)完成整個(gè)振蕩周期,此時(shí)形變量仍有0.111,呈現(xiàn)出長軸與水平方向的橢球狀,這一點(diǎn)與低黏度的影響結(jié)果較為相似。

    圖9 不同表面張力的150μm霧滴在橫向風(fēng)場中形變量e隨時(shí)間的變化

    2.3 霧滴形變預(yù)測關(guān)聯(lián)式

    本文中影響霧滴空間變形特征的因素有霧滴直徑、環(huán)境風(fēng)速、霧滴黏度和表面張力,并將霧滴的空間形變歸結(jié)于重力、曳力、黏性力和表面張力的共同作用。而Re 表征慣性力與黏性力之比,其中慣性力反映出霧滴速度對(duì)空間形變過程的影響,同時(shí)霧滴與空氣的相對(duì)速度又影響著曳力,因此曳力對(duì)形變過程的影響可通過慣性力進(jìn)行描述;We 表征慣性力與表面張力之比;Oh 表征黏性力與表面張力之比。結(jié)合本文中影響霧滴空間形變的物理參數(shù)以及上述量綱為1 參數(shù)的意義,選擇了韋伯?dāng)?shù)We、奧內(nèi)佐格數(shù)Oh和雷諾數(shù)Re的關(guān)聯(lián)式對(duì)形變量和形變周期作出預(yù)測。計(jì)算公式分別見式(14)~式(16)。

    其中,當(dāng)量直徑De選取為霧滴的初始直徑,ud和ua分別為霧滴的初始速度和橫向空氣速度。經(jīng)過計(jì)算后,擬合所涉及的量綱為1 數(shù)范圍依次為:138.874<We<573.466, 66.730<Re<207.950,2.873×10-3<Oh<1.916×10-2。

    采用多元非線性回歸方法對(duì)計(jì)算結(jié)果進(jìn)行擬合,結(jié)果如式(17)所示。

    其中確定系數(shù)R2= 0.9497,如圖10 所示對(duì)比分析了回歸模型的結(jié)果與實(shí)際模擬的結(jié)果,霧滴形變量預(yù)測模型的平均相對(duì)偏差為8.40%。霧滴的最大形變量與液相黏度和表面張力呈負(fù)相關(guān),說明液相黏性力和表面張力能夠阻礙霧滴形變,使其維持原有形態(tài),而最大形變比與相對(duì)速度和初始直徑呈正相關(guān),反映出曳力與重力對(duì)形變的增強(qiáng)作用。擬合分析結(jié)果與上文討論結(jié)果基本一致。

    圖10 霧滴形變量模擬結(jié)果與回歸模型的偏差比較

    同時(shí),對(duì)初次形變時(shí)間占一次振蕩形變周期的比例(Δt1/Δtd) 進(jìn)行擬合,得到如下結(jié)果,見式(18)。

    其中確定系數(shù)R2=0.9150,如圖11 所示,對(duì)比分析了回歸模型的結(jié)果與實(shí)際模擬的結(jié)果,該預(yù)測模型的平均相對(duì)偏差為5.83%。當(dāng)霧滴特征參數(shù)及環(huán)境風(fēng)速已知時(shí),可以通過式(17)和式(18)分別預(yù)測出霧滴的最大形變量和形變周期,以此明晰霧滴的空間形態(tài)演化規(guī)律,為提出霧滴群的減飄手段奠定理論基礎(chǔ)。

    圖11 霧滴初次形變時(shí)間占比模擬結(jié)果與回歸模型偏差比較

    3 結(jié)論

    通過數(shù)值模擬手段對(duì)多因素耦合下單霧滴的空間形態(tài)演化過程進(jìn)行研究,結(jié)論如下。

    (1)霧滴在下落過程中的形變存在振蕩周期,且形變會(huì)發(fā)生衰減。霧滴的尺寸嚴(yán)重影響著振蕩周期和幅度,直徑越大,霧滴的重力和曳力越大,形變周期越長,形變幅度越大;橫向風(fēng)速幾乎不改變形變周期,但風(fēng)速增大會(huì)加劇霧滴形變程度;霧滴黏度主要影響形變程度,黏度越小,形變程度越大,同時(shí)在低黏度下霧滴完成一個(gè)形變周期后難以恢復(fù)到初始的球形,只能保持橢球狀繼續(xù)運(yùn)動(dòng);表面張力值增大導(dǎo)致形變周期縮短,形變程度減弱,當(dāng)表面張力很小時(shí)霧滴也難以恢復(fù)球形。

    (2)霧滴的形變量受霧滴直徑、橫向風(fēng)速、藥液黏度、表面張力等因素的共同影響,即形變量是一個(gè)由重力、曳力、黏性力及表面張力共同作用的結(jié)果;曳力和重力促使霧滴發(fā)生形變,而黏性力和表面張力阻礙霧滴發(fā)生形變。

    (3)通過回歸分析得到的霧滴形變量與形變周期的量綱為1預(yù)測模型,可以計(jì)算出霧滴的空間形變量及形變周期,可明晰霧滴的空間形態(tài)演化規(guī)律,為提出農(nóng)藥的減飄手段奠定理論基礎(chǔ)。

    符號(hào)說明

    A—— 霧滴在風(fēng)向的投影長度,m

    a—— 橢圓長軸,m

    b—— 橢圓短軸,m

    Cd—— 曳力系數(shù)

    De—— 霧滴的當(dāng)量直徑,m

    E—— 形變比

    e—— 形變量

    emax1—— 第一個(gè)振蕩周期內(nèi),初次變形過程的最大形變量

    Emax2—— 第一個(gè)振蕩周期內(nèi),二次變形過程的最大形變量

    Fd—— 單位體積霧滴所受曳力,N/m3

    Fσ—— 單位體積流體的表面張力,N/m3

    g—— 重力加速度,m/s2

    n—— 相界面法向向量

    Oh—— 奧內(nèi)佐格數(shù)

    p—— 單位體積壓力,N/m3

    patm—— 大氣壓,Pa

    pout—— 出口壓力,Pa

    R2—— 擬合相關(guān)系數(shù)

    Re—— 雷諾數(shù)

    u—— 流體的速度,m/s

    ua—— 空氣的速度,m/s

    ua,0—— 空氣的初始速度,m/s

    ud—— 霧滴的速度,m/s

    ud,0—— 霧滴的初始速度,m/s

    uin—— 入口速度,m/s

    We—— 韋伯?dāng)?shù)

    x—— 到相界面距離,m

    Δt1—— 第一個(gè)振蕩周期內(nèi),初次變形時(shí)間,s

    Δt2—— 第一個(gè)振蕩周期內(nèi),二次變形時(shí)間,s

    Δtd—— 第一個(gè)振蕩周期,s

    δ—— Delta函數(shù)

    κ—— 相界面曲率

    μ—— 流體黏度,Pa?s

    μa—— 空氣的黏度,Pa?s

    μd—— 霧滴的黏度,Pa?s

    ρ—— 流體的密度,kg/m3

    ρa(bǔ)—— 空氣的密度,kg/m3

    ρd—— 霧滴的密度,kg/m3

    σ—— 表面張力,N/m

    φ—— 水平集函數(shù)

    下角標(biāo)

    a—— 空氣相

    d—— 霧滴相

    max—— 最大值

    猜你喜歡
    曳力橢球表面張力
    預(yù)測天然氣斜井臨界攜液流量新方法
    獨(dú)立坐標(biāo)系橢球變換與坐標(biāo)換算
    循環(huán)流化床鍋爐爐膛流動(dòng)特性數(shù)值模擬進(jìn)展
    橢球槽宏程序編制及其Vericut仿真
    智能制造(2021年4期)2021-11-04 08:54:44
    橢球精加工軌跡及程序設(shè)計(jì)
    基于外定界橢球集員估計(jì)的純方位目標(biāo)跟蹤
    神奇的表面張力
    小布老虎(2016年4期)2016-12-01 05:46:08
    MgO-B2O3-SiO2三元體系熔渣表面張力計(jì)算
    上海金屬(2016年2期)2016-11-23 05:34:45
    基于EMMS模型的攪拌釜內(nèi)氣液兩相流數(shù)值模擬
    CaF2-CaO-Al2O3-MgO-SiO2渣系表面張力計(jì)算模型
    上海金屬(2014年3期)2014-12-19 13:09:06
    精品国产乱子伦一区二区三区 | 亚洲国产欧美一区二区综合| 亚洲,欧美精品.| 人妻一区二区av| 在线天堂中文资源库| 日韩三级视频一区二区三区| 国产有黄有色有爽视频| 大型av网站在线播放| 黑人欧美特级aaaaaa片| 国产片内射在线| 国产成人欧美在线观看 | 老熟妇仑乱视频hdxx| 在线观看免费高清a一片| 窝窝影院91人妻| 国产免费一区二区三区四区乱码| 日韩免费高清中文字幕av| 中文字幕精品免费在线观看视频| 高清欧美精品videossex| 亚洲精品国产av成人精品| www日本在线高清视频| 色婷婷久久久亚洲欧美| 国产精品一区二区精品视频观看| 国产视频一区二区在线看| 另类精品久久| 午夜免费成人在线视频| 国产又爽黄色视频| 精品乱码久久久久久99久播| 国产av又大| 亚洲精品国产色婷婷电影| 久久久久网色| avwww免费| 国内毛片毛片毛片毛片毛片| 久久久国产一区二区| 无限看片的www在线观看| 成人国产一区最新在线观看| 午夜福利免费观看在线| 一进一出抽搐动态| 在线观看舔阴道视频| 韩国精品一区二区三区| 久久中文字幕一级| 精品欧美一区二区三区在线| 欧美人与性动交α欧美软件| 啦啦啦免费观看视频1| 天堂中文最新版在线下载| 色婷婷av一区二区三区视频| 国产精品久久久久成人av| 久久久久久久国产电影| 色94色欧美一区二区| 操出白浆在线播放| 日韩免费高清中文字幕av| 精品亚洲成国产av| 19禁男女啪啪无遮挡网站| 亚洲一码二码三码区别大吗| 国产精品亚洲av一区麻豆| 美国免费a级毛片| 精品人妻1区二区| 99精国产麻豆久久婷婷| 免费不卡黄色视频| 色婷婷久久久亚洲欧美| 亚洲国产看品久久| 日本黄色日本黄色录像| 国产成人av激情在线播放| 操美女的视频在线观看| 狂野欧美激情性xxxx| 成人三级做爰电影| 正在播放国产对白刺激| 制服人妻中文乱码| 啦啦啦在线免费观看视频4| 国产精品久久久久成人av| 秋霞在线观看毛片| 亚洲av成人不卡在线观看播放网 | 91国产中文字幕| 18禁黄网站禁片午夜丰满| av又黄又爽大尺度在线免费看| 亚洲第一欧美日韩一区二区三区 | av网站免费在线观看视频| 高清av免费在线| 一边摸一边抽搐一进一出视频| 国产又爽黄色视频| 一区二区三区精品91| 精品第一国产精品| 999久久久精品免费观看国产| 中文字幕人妻丝袜一区二区| 日韩制服骚丝袜av| 一区二区三区乱码不卡18| 午夜福利视频在线观看免费| 这个男人来自地球电影免费观看| 欧美精品av麻豆av| 老司机影院成人| 嫩草影视91久久| 久热这里只有精品99| 99re6热这里在线精品视频| 90打野战视频偷拍视频| 免费观看av网站的网址| 老鸭窝网址在线观看| 麻豆国产av国片精品| 一本—道久久a久久精品蜜桃钙片| 国产高清国产精品国产三级| 午夜91福利影院| 国产日韩欧美视频二区| 人人妻,人人澡人人爽秒播| 国产欧美日韩一区二区三 | 亚洲国产av新网站| www.熟女人妻精品国产| 国产亚洲精品一区二区www | 法律面前人人平等表现在哪些方面 | 我要看黄色一级片免费的| 国产精品麻豆人妻色哟哟久久| cao死你这个sao货| 丰满迷人的少妇在线观看| 超碰97精品在线观看| 亚洲精品国产精品久久久不卡| 亚洲天堂av无毛| 国产亚洲午夜精品一区二区久久| 深夜精品福利| 国产av精品麻豆| 99久久精品国产亚洲精品| 自拍欧美九色日韩亚洲蝌蚪91| 50天的宝宝边吃奶边哭怎么回事| 丰满少妇做爰视频| 少妇 在线观看| 老熟女久久久| 成年女人毛片免费观看观看9 | 午夜免费成人在线视频| 免费在线观看完整版高清| 免费观看av网站的网址| 精品熟女少妇八av免费久了| 久久ye,这里只有精品| 国产在视频线精品| 久久人人爽人人片av| 人妻 亚洲 视频| 中亚洲国语对白在线视频| 亚洲av电影在线观看一区二区三区| 不卡一级毛片| 国产免费福利视频在线观看| 久久99热这里只频精品6学生| www.自偷自拍.com| 一级,二级,三级黄色视频| 久久久久久久久免费视频了| 一个人免费看片子| 亚洲熟女精品中文字幕| 成年人午夜在线观看视频| 久久天堂一区二区三区四区| 极品人妻少妇av视频| 99国产精品一区二区三区| 人妻 亚洲 视频| 亚洲天堂av无毛| 精品久久久精品久久久| 久久毛片免费看一区二区三区| 国产av又大| 欧美乱码精品一区二区三区| 在线永久观看黄色视频| 欧美日韩福利视频一区二区| 手机成人av网站| www.999成人在线观看| 国产成人av教育| 亚洲自偷自拍图片 自拍| 人妻 亚洲 视频| 免费观看a级毛片全部| 精品国产国语对白av| 久久亚洲精品不卡| 午夜老司机福利片| 亚洲国产精品成人久久小说| 欧美日韩中文字幕国产精品一区二区三区 | 嫁个100分男人电影在线观看| 丁香六月欧美| 搡老岳熟女国产| 一区二区三区精品91| 亚洲av成人一区二区三| 国产不卡av网站在线观看| av福利片在线| 老司机深夜福利视频在线观看 | 精品国产乱子伦一区二区三区 | 欧美大码av| 成人av一区二区三区在线看 | 水蜜桃什么品种好| 久久久久精品人妻al黑| 黑人猛操日本美女一级片| 欧美亚洲日本最大视频资源| 最近最新中文字幕大全免费视频| 99香蕉大伊视频| 精品国产乱码久久久久久小说| 天堂中文最新版在线下载| 国产精品秋霞免费鲁丝片| 18禁黄网站禁片午夜丰满| 黑人欧美特级aaaaaa片| 男人爽女人下面视频在线观看| 亚洲全国av大片| 久久久精品免费免费高清| 国产一区二区三区在线臀色熟女 | 午夜老司机福利片| 亚洲情色 制服丝袜| 日韩 亚洲 欧美在线| 久久亚洲精品不卡| 最近最新中文字幕大全免费视频| 99久久国产精品久久久| 9色porny在线观看| 麻豆乱淫一区二区| 黄片小视频在线播放| 国产国语露脸激情在线看| 岛国在线观看网站| 亚洲性夜色夜夜综合| 免费av中文字幕在线| 悠悠久久av| 精品免费久久久久久久清纯 | 亚洲精品国产一区二区精华液| 久久久久久免费高清国产稀缺| 香蕉丝袜av| 国产片内射在线| 丝瓜视频免费看黄片| 美女大奶头黄色视频| 久久99热这里只频精品6学生| 亚洲国产毛片av蜜桃av| 免费在线观看黄色视频的| 欧美激情久久久久久爽电影 | 日本一区二区免费在线视频| 久热爱精品视频在线9| 日韩熟女老妇一区二区性免费视频| 午夜激情久久久久久久| 日本五十路高清| 亚洲人成电影观看| 午夜91福利影院| 久久中文看片网| 99久久精品国产亚洲精品| 国产日韩欧美视频二区| 性色av乱码一区二区三区2| 高潮久久久久久久久久久不卡| videos熟女内射| 亚洲人成电影观看| 中文字幕精品免费在线观看视频| 大香蕉久久网| 成人18禁高潮啪啪吃奶动态图| 人人妻,人人澡人人爽秒播| 最新的欧美精品一区二区| 亚洲欧美一区二区三区久久| 日本黄色日本黄色录像| 久久午夜综合久久蜜桃| 丁香六月天网| 一区二区三区精品91| 高清欧美精品videossex| 91精品三级在线观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲人成电影免费在线| 久热爱精品视频在线9| 大陆偷拍与自拍| 丰满迷人的少妇在线观看| tube8黄色片| 国产一区二区三区av在线| 欧美国产精品va在线观看不卡| 亚洲精品美女久久av网站| 多毛熟女@视频| 成人国产一区最新在线观看| 色婷婷久久久亚洲欧美| 日韩三级视频一区二区三区| 青青草视频在线视频观看| 亚洲欧美色中文字幕在线| 欧美成狂野欧美在线观看| 免费在线观看视频国产中文字幕亚洲 | 国产男人的电影天堂91| 女人高潮潮喷娇喘18禁视频| 麻豆国产av国片精品| 老司机在亚洲福利影院| 亚洲av日韩在线播放| 啦啦啦中文免费视频观看日本| 青春草视频在线免费观看| 黄色片一级片一级黄色片| 一边摸一边做爽爽视频免费| 亚洲精华国产精华精| 亚洲国产欧美网| 精品乱码久久久久久99久播| 国产男人的电影天堂91| 国产一区有黄有色的免费视频| 亚洲一卡2卡3卡4卡5卡精品中文| 建设人人有责人人尽责人人享有的| 国产伦人伦偷精品视频| 亚洲av片天天在线观看| 视频区欧美日本亚洲| 中国国产av一级| 男人添女人高潮全过程视频| 成年av动漫网址| 精品久久久久久久毛片微露脸 | 亚洲 国产 在线| 十八禁网站免费在线| 国产精品秋霞免费鲁丝片| 热re99久久国产66热| 国产亚洲精品第一综合不卡| 久久人妻福利社区极品人妻图片| 搡老乐熟女国产| 欧美 日韩 精品 国产| 美女大奶头黄色视频| 人人妻人人添人人爽欧美一区卜| 日本欧美视频一区| 国产黄色免费在线视频| 午夜福利视频精品| 日本五十路高清| 国产精品av久久久久免费| 首页视频小说图片口味搜索| 欧美精品亚洲一区二区| 国产主播在线观看一区二区| 亚洲专区中文字幕在线| 国产精品久久久av美女十八| 国产福利在线免费观看视频| 午夜激情久久久久久久| 精品少妇一区二区三区视频日本电影| 国产主播在线观看一区二区| 一区二区三区精品91| 黄色怎么调成土黄色| 日韩欧美国产一区二区入口| 久久人妻福利社区极品人妻图片| 丰满人妻熟妇乱又伦精品不卡| 热re99久久精品国产66热6| 亚洲精品粉嫩美女一区| 十分钟在线观看高清视频www| a级毛片黄视频| 亚洲,欧美精品.| 三级毛片av免费| 亚洲精品国产av成人精品| 成年动漫av网址| 免费人妻精品一区二区三区视频| 男男h啪啪无遮挡| 亚洲精品一二三| 亚洲精品国产一区二区精华液| 亚洲成人免费av在线播放| 国产国语露脸激情在线看| 国产一级毛片在线| 99精品久久久久人妻精品| 激情视频va一区二区三区| 交换朋友夫妻互换小说| 午夜激情av网站| 亚洲国产精品一区二区三区在线| 色老头精品视频在线观看| 91av网站免费观看| 亚洲专区字幕在线| 免费女性裸体啪啪无遮挡网站| 99久久99久久久精品蜜桃| 啦啦啦 在线观看视频| 午夜激情久久久久久久| 岛国在线观看网站| xxxhd国产人妻xxx| 亚洲精品在线美女| 91麻豆精品激情在线观看国产 | 国产麻豆69| 色老头精品视频在线观看| 国产91精品成人一区二区三区 | 亚洲天堂av无毛| 无限看片的www在线观看| 久久久欧美国产精品| 国产熟女午夜一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 久久亚洲精品不卡| 美女主播在线视频| 亚洲国产欧美网| 久久亚洲精品不卡| 日本wwww免费看| 曰老女人黄片| 亚洲精品美女久久久久99蜜臀| 亚洲国产欧美在线一区| 建设人人有责人人尽责人人享有的| 欧美精品人与动牲交sv欧美| 水蜜桃什么品种好| 亚洲人成电影观看| 妹子高潮喷水视频| 免费在线观看黄色视频的| 欧美日韩亚洲高清精品| 三级毛片av免费| 在线亚洲精品国产二区图片欧美| 99香蕉大伊视频| 久久久久久久大尺度免费视频| 久9热在线精品视频| 亚洲国产欧美日韩在线播放| 日本一区二区免费在线视频| 午夜免费成人在线视频| 首页视频小说图片口味搜索| 青春草视频在线免费观看| 老司机福利观看| 国产淫语在线视频| 亚洲欧洲日产国产| 宅男免费午夜| 人人澡人人妻人| 午夜福利,免费看| 黑人欧美特级aaaaaa片| 欧美日本中文国产一区发布| 午夜福利视频精品| 久久中文字幕一级| 曰老女人黄片| 又大又爽又粗| 国产成人欧美| 精品国产一区二区久久| 亚洲精品一区蜜桃| 精品久久蜜臀av无| 国产精品久久久av美女十八| 亚洲国产中文字幕在线视频| 宅男免费午夜| 啦啦啦中文免费视频观看日本| 亚洲精品国产精品久久久不卡| 欧美久久黑人一区二区| 天天躁狠狠躁夜夜躁狠狠躁| av有码第一页| 99精国产麻豆久久婷婷| 亚洲国产精品一区三区| 久久国产精品男人的天堂亚洲| 国产黄色免费在线视频| 美国免费a级毛片| 欧美黑人欧美精品刺激| 久久久久久亚洲精品国产蜜桃av| 国产精品1区2区在线观看. | 精品久久蜜臀av无| 在线十欧美十亚洲十日本专区| 日本一区二区免费在线视频| 精品一区在线观看国产| 亚洲精品一区蜜桃| 99国产精品一区二区蜜桃av | 美女中出高潮动态图| 国产极品粉嫩免费观看在线| 久久久久国产精品人妻一区二区| 日韩视频一区二区在线观看| 正在播放国产对白刺激| av在线老鸭窝| 欧美日韩一级在线毛片| 午夜激情av网站| 免费高清在线观看日韩| 国产精品一区二区在线不卡| 久久国产精品男人的天堂亚洲| 操出白浆在线播放| 一二三四在线观看免费中文在| 91av网站免费观看| 丁香六月欧美| 久久亚洲国产成人精品v| 国产精品一区二区免费欧美 | 色婷婷久久久亚洲欧美| 这个男人来自地球电影免费观看| 久久久久久久大尺度免费视频| 热99re8久久精品国产| 国产亚洲av片在线观看秒播厂| 中文精品一卡2卡3卡4更新| 少妇被粗大的猛进出69影院| 黄色怎么调成土黄色| 90打野战视频偷拍视频| 我要看黄色一级片免费的| 母亲3免费完整高清在线观看| 97精品久久久久久久久久精品| 日本av免费视频播放| 久久久久久免费高清国产稀缺| 99精品久久久久人妻精品| 97在线人人人人妻| 国产精品免费大片| 99re6热这里在线精品视频| 国产成人a∨麻豆精品| 亚洲九九香蕉| 一二三四在线观看免费中文在| 国产精品偷伦视频观看了| 18禁观看日本| 国产精品亚洲av一区麻豆| 午夜福利乱码中文字幕| 人人澡人人妻人| 午夜福利免费观看在线| 成年人午夜在线观看视频| 久久久久国内视频| 在线观看www视频免费| 日本一区二区免费在线视频| 国产真人三级小视频在线观看| 伦理电影免费视频| 在线天堂中文资源库| 91九色精品人成在线观看| 国产精品久久久久久人妻精品电影 | 久久久久国产精品人妻一区二区| a在线观看视频网站| 久久久久久久国产电影| 久久久久国产一级毛片高清牌| 午夜福利一区二区在线看| 两个人看的免费小视频| 亚洲人成77777在线视频| a级毛片黄视频| 一区二区三区精品91| 嫩草影视91久久| 欧美一级毛片孕妇| 欧美黑人欧美精品刺激| 99九九在线精品视频| 美女大奶头黄色视频| 欧美日韩黄片免| 少妇人妻久久综合中文| 久久精品熟女亚洲av麻豆精品| 精品一区二区三区av网在线观看 | 亚洲精华国产精华精| 欧美亚洲日本最大视频资源| 国产成人欧美在线观看 | 成年动漫av网址| 人人澡人人妻人| 国产精品一区二区在线观看99| 手机成人av网站| 日本av免费视频播放| 丁香六月天网| 搡老乐熟女国产| 色94色欧美一区二区| 美女国产高潮福利片在线看| 久久香蕉激情| 巨乳人妻的诱惑在线观看| 久久久久久久久免费视频了| 99国产综合亚洲精品| 水蜜桃什么品种好| 精品国产乱码久久久久久小说| 亚洲国产毛片av蜜桃av| 两个人看的免费小视频| 美女脱内裤让男人舔精品视频| 下体分泌物呈黄色| 69精品国产乱码久久久| 国产xxxxx性猛交| 50天的宝宝边吃奶边哭怎么回事| 香蕉丝袜av| 日韩中文字幕欧美一区二区| 日本vs欧美在线观看视频| 欧美变态另类bdsm刘玥| 窝窝影院91人妻| 高清黄色对白视频在线免费看| 精品人妻在线不人妻| 啦啦啦中文免费视频观看日本| 日本a在线网址| 在线观看www视频免费| 国内毛片毛片毛片毛片毛片| 丝袜脚勾引网站| 亚洲人成电影观看| 国产日韩欧美亚洲二区| 国产日韩一区二区三区精品不卡| 久久人妻福利社区极品人妻图片| 亚洲中文字幕日韩| 一本—道久久a久久精品蜜桃钙片| 蜜桃在线观看..| 9191精品国产免费久久| 中国美女看黄片| 亚洲精品一区蜜桃| 久久影院123| 久久中文字幕一级| 亚洲精品在线美女| 精品久久蜜臀av无| 日本黄色日本黄色录像| 美女国产高潮福利片在线看| 国产欧美日韩综合在线一区二区| 久久精品aⅴ一区二区三区四区| 免费一级毛片在线播放高清视频 | 婷婷丁香在线五月| 老汉色av国产亚洲站长工具| 欧美 亚洲 国产 日韩一| 欧美精品一区二区免费开放| 免费少妇av软件| 我的亚洲天堂| 免费观看av网站的网址| 老司机靠b影院| 老司机福利观看| 啦啦啦 在线观看视频| 免费在线观看视频国产中文字幕亚洲 | 精品久久久精品久久久| 亚洲七黄色美女视频| 欧美日韩av久久| 国产一区二区三区av在线| 日本欧美视频一区| av天堂久久9| 亚洲精品美女久久av网站| 大片免费播放器 马上看| 欧美精品人与动牲交sv欧美| 久久精品国产a三级三级三级| 精品一区在线观看国产| 国内毛片毛片毛片毛片毛片| 热re99久久精品国产66热6| 777米奇影视久久| 啦啦啦 在线观看视频| 亚洲av成人一区二区三| 亚洲七黄色美女视频| 久热爱精品视频在线9| 成人手机av| 黄色 视频免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品中文字幕在线视频| tocl精华| 午夜精品国产一区二区电影| 操美女的视频在线观看| 国产精品久久久av美女十八| 美女扒开内裤让男人捅视频| 国产精品久久久人人做人人爽| 国产片内射在线| 亚洲成av片中文字幕在线观看| 一区二区三区四区激情视频| 19禁男女啪啪无遮挡网站| 亚洲av成人一区二区三| 国产伦人伦偷精品视频| 午夜激情av网站| 免费观看a级毛片全部| 欧美+亚洲+日韩+国产| av福利片在线| 狠狠精品人妻久久久久久综合| 国产成人av教育| 国产亚洲精品第一综合不卡| 欧美一级毛片孕妇| 久久久精品免费免费高清| 天堂8中文在线网| 国产三级黄色录像| 亚洲欧美清纯卡通| 亚洲欧美精品自产自拍| 亚洲专区国产一区二区| 别揉我奶头~嗯~啊~动态视频 | 99热国产这里只有精品6| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧洲日产国产| 日本黄色日本黄色录像| 欧美日韩国产mv在线观看视频| 纵有疾风起免费观看全集完整版| 久久国产精品影院| 99香蕉大伊视频| av在线老鸭窝| 伦理电影免费视频| 欧美激情极品国产一区二区三区| 一本大道久久a久久精品| 麻豆av在线久日| 久久久久久久国产电影| 伦理电影免费视频| 精品免费久久久久久久清纯 |