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

    柱形料腔中超聲空化效應(yīng)的空間分布特性研究

    2020-05-09 08:34:34俞劍劉斌徐雪冬馮麗麗
    聲學(xué)技術(shù) 2020年2期
    關(guān)鍵詞:魯米諾柱形變幅

    俞劍,劉斌,徐雪冬,馮麗麗

    柱形料腔中超聲空化效應(yīng)的空間分布特性研究

    俞劍,劉斌,徐雪冬,馮麗麗

    (北京工商大學(xué)材料與機械工程學(xué)院,北京 100048)

    研究了20 kHz的圓柱形料腔中超聲空化效應(yīng)的形成及其空間分布特性。應(yīng)用柱貝塞爾函數(shù),推導(dǎo)獲取了柱形聲場內(nèi)超聲傳播的聲能密度的分布,并采用有限元方法進行仿真分析。針對頻率為20 kHz的功率超聲實驗,結(jié)合聲學(xué)測量方法和魯米諾聲致化學(xué)發(fā)光方法,對理論分析結(jié)果進行了驗證對照。結(jié)果表明:料腔半徑=50 mm,20 kHz諧振液位高度=90 mm時,若功放電流<40 mA,超聲空化效應(yīng)出現(xiàn)在變幅桿端部區(qū)域;若40 mA≤功放電流≤80 mA,空化效應(yīng)顯著增強,空化效應(yīng)的空間分布與場內(nèi)聲壓分布一致,空化效應(yīng)受聲模態(tài)影響,形成遠場空化效應(yīng)的分布特性;若功放電流>80 mA,受非線性因素影響,諧振液位時,空化效應(yīng)在聲流作用下呈柱形拖尾狀分布,并在底部壁面邊界形成平鋪狀分布;非諧振液位高度等于75 mm時,超聲空化效應(yīng)隨功率增加僅在變幅桿端部區(qū)域出現(xiàn),且呈現(xiàn)局域空化分布特性。

    柱形料腔;空化效應(yīng);功率超聲;諧振液位;魯米諾聲致化學(xué)發(fā)光法

    0 引言

    超聲波破碎是利用聲能對液體內(nèi)的固體顆?;蛏锝M織實施破碎的過程,與傳統(tǒng)破碎技術(shù)相比,超聲破碎可獲得更為均勻和細小的粒徑,并具有更高的能效[1-5]??栈?yīng)作為超聲破碎中的主要作用機制,理解超聲空化效應(yīng)的形成,如何增強空化效應(yīng)的作用范圍,提高超聲能效的利用已成為超聲研究的重要內(nèi)容[6-7]。

    超聲空化效應(yīng)是由料液聲場中的交變壓力所引發(fā),超聲波的輸入頻率[8-11]、振幅[8,11-12]、聲場條件[12-15]、料液空化閾值[11,16]是影響空化效應(yīng)的直接因素。當超聲波頻率、料液確定時,提高超聲波的輸入振幅,改變聲場條件是增強超聲空化的有效措施。提高超聲波輸入振幅主要有兩種途徑:(1) 增大變幅桿功率輸入[14];(2) 設(shè)計高效能的變幅桿結(jié)構(gòu)[17-19]。改變聲場條件則考慮料腔聲場的幾何形狀與邊界條件,以獲取更為高效的聲能密度分布,以利于空化效應(yīng)的增強[20-23]。

    已有研究表明,低頻功率超聲具有更強的超聲空化作用,現(xiàn)有的超聲破碎機頻率多采用20 kHz。本研究以頻率為20 kHz的浸入式圓柱形超聲破碎裝置為研究對象,從聲場能效增強利用的角度出發(fā),應(yīng)用柱貝塞爾函數(shù)推導(dǎo)獲取了柱形料液場的聲模態(tài),結(jié)合諧振強化,針對20 kHz超聲頻率設(shè)計了柱形容器的半徑和液位高度,在聲學(xué)有限元分析的基礎(chǔ)上,分析了該聲場條件下的聲能密度分布特性,并與非諧振液位進行了對比研究。最后采用直接聲場測量與間接魯米諾聲致化學(xué)發(fā)光方法相結(jié)合(Sono Chemi Luminescent, SCL),對理論分析結(jié)果進行了實驗驗證,獲取超聲空化效應(yīng)在柱形料腔中的空間分布特性。

    1 材料與方法

    1.1 材料與儀器

    1.1.1 材料與試劑

    魯米諾(luminol,C8H7N3O2,一種化學(xué)發(fā)光劑),純度98%,Ark Pharm公司;固體NaOH,純度96%,阿拉丁工業(yè)公司;反應(yīng)料液采用去離子水;半徑50 mm、高度120 mm的柱形透明容器。

    1.1.2 主要儀器設(shè)備

    超聲波處理單元構(gòu)成:超聲功率放大器,額定功率為200 W,變幅桿式超聲換能器,輸出端半徑為6.5 mm,輸出頻率為20 kHz,中國科學(xué)院聲學(xué)研究所定制;測振系統(tǒng):速度型激光測振儀,速度靈敏度為200 mm·s-1,B&K公司;多通道數(shù)字示波器:TDS 2024B數(shù)字彩色示波器,測試帶寬為200 MHz,采樣率為2 GHz,4通道,泰克公司。

    1.2 方法

    1.2.1 實驗方案及系統(tǒng)

    超聲波在水中傳播時,產(chǎn)生大量空化泡,空化泡潰滅時產(chǎn)生的高溫高壓會形成光輻射現(xiàn)象,即“聲致發(fā)光”。由于其光輻射能量較低,常溫狀態(tài)下肉眼難以觀察。為了強化其光輻射能量,可加入微量增強劑魯米諾(luminol)。其作用原理為超聲空化氣泡潰滅時產(chǎn)生的高溫高壓可水解產(chǎn)生羥基(OH自由基),并與魯米諾發(fā)生反應(yīng),羥基化的魯米醇顯示強烈的化學(xué)發(fā)光,這種現(xiàn)象被稱為聲致化學(xué)發(fā)光(SCL)[15]。發(fā)光強度由羥基數(shù)量決定,即超聲空化強度大小(反映了聲場能量的分布特性)決定。為了有效增加羥基數(shù)量,可在水溶液中加入適量NaOH或KOH,魯米諾濃度宜采用0.1 g·L-1,NaOH的添加量為1 g·L-1[12]。

    為觀察發(fā)光現(xiàn)象,實驗采用暗室條件,聲致化學(xué)發(fā)光可形成淡藍色的光束,肉眼可見,采用數(shù)碼相機記錄圖片,曝光時間>3 min。料腔工況條件為:底部懸置點支撐條件,使料腔底部與料腔徑向壁面的邊界條件統(tǒng)一為軟邊界條件。

    測試系統(tǒng)簡圖如圖1所示。圖中,為柱形料腔半徑,為變幅桿浸入水中深度,為液位高度。

    圖1 測試系統(tǒng)簡圖

    1.2.2 圓柱管腔內(nèi)的超聲傳播特性

    理想均質(zhì)流體媒介中小振幅聲波的波動方程可表示為

    超聲破碎料腔采用圓柱形管腔形式,根據(jù)水和空氣的聲阻抗特性及料腔壁厚,可將液體-空氣界面、液體-料腔壁-空氣界面均可簡化為軟壁面邊界條件,腔內(nèi)形成高次波的傳播形式,聲壓可表示為軸向波束和徑向波束的合成:

    由軟壁面邊界條件和對稱聲源特性,可以得到圓柱形管腔內(nèi)高次波貝塞爾函數(shù)為

    由此可得:

    研究表明,當料液腔高度取聲模態(tài)液位高度時,料液腔聲場中的聲能分布大大提高,由此可知,超聲破碎應(yīng)用中,當料腔中料液物性、邊界條件、聲源強度一定時,聲源頻率達到料腔聲場的聲模態(tài)頻率(諧振頻率),可增強聲傳播能量,提高空化能利用率并有效增加空化作用范圍。

    1.2.3 聲場有限元模型

    有限元分析采用COMSOL軟件。聲場幾何形狀為圓柱形腔體,聲場模型的網(wǎng)格劃分:根據(jù)軸對稱條件,可簡化為軸對稱二維平面,取一半進行四邊形網(wǎng)格劃分,網(wǎng)格最大邊長為4 mm,計算上限頻率約為62 kHz。傳聲媒介的特性取值:設(shè)定水的密度為1 000 kg·m-3,聲速取1 480 m·s-1。邊界條件:聲源邊界條件是變幅桿端部振動速度(實測值);壁面邊界條件是變幅桿側(cè)面-料液界面為剛性壁硬邊界全反射條件;料液-空氣界面為軟邊界全反射條件;料液-玻璃壁面-空氣可視為料液-空氣軟壁面邊界全反射條件,由于20 kHz超聲在中間層玻璃中的波長(縱波)約為280 mm,遠大于玻璃壁厚(2 mm),根據(jù)聲波通過中間層的傳播規(guī)律,玻璃層透射系數(shù)約為1,因此可忽略中間層玻璃的影響,簡化為料液-空氣邊界條件。

    2 結(jié)果與分析

    2.1 聲源邊界條件的確定

    由于聲源邊界條件為超聲變幅桿端部振動速度,因此需獲取超聲變幅桿的工況頻率和振速幅值作為仿真聲源邊界條件,以確保仿真的有效性。取功放電流為40 mA和70 mA,變幅桿端部振速幅值隨頻率的變化如圖2所示。

    圖2 聲源振速頻譜圖

    變幅桿出廠標稱頻率為20 kHz。由圖2可知,其工況頻率為20.008 kHz,頻率誤差為0.04%,對于20 kHz條件下模態(tài)液位計算的影響較小,對于平均聲能量密度幅值及分布的量化分析計算影響可忽略。仿真分析可采用工況頻率20 kHz。

    圖2表明,工況20 kHz頻率時,變幅桿端面的縱向振動速度達到最大,超聲功放電流為40 mA和70 mA時,振速分別為:0.516 m·s-1,0.899 m·s-1。

    2.2 諧振料腔聲場的聲能分布特性

    根據(jù)定容聲場條件,以半徑50 mm圓柱形料腔為例,取為1的柱貝塞爾函數(shù)解(高階解時,徑傳播的聲壓幅值衰減較大)。聲源工作頻率20 kHz,由式(5)中的諧振液位高度計算公式,可求得軸向簡正波數(shù)為=2時的液位高度=90 mm。另取一非諧振液位高度=75 mm進行對照??疾靸煞N液位高度下聲場中的聲能的分布特性。

    取功放電流為70 mA,變幅桿浸入料液腔深度為10mm作為聲源邊界條件。兩種液位高度的料腔聲場中的聲壓分布云圖如圖3所示。

    液位高度為75 mm時,該聲場條件下的一、二階諧振頻率與20 kHz偏差較大,在聲源頻率為20 kHz時不能形成諧振條件,場內(nèi)聲壓分布如圖3(a)所示。液位高度為90 mm時,該聲場條件下諧振頻率為20 kHz時,波束沿軸線方向形成2個簡正波數(shù),表現(xiàn)為二階聲模態(tài)特性,其聲壓分布呈雙橢圓上下對稱分布,兩聲壓極值位置分別位于料腔高度的1/3和2/3處,對稱面位于料腔高度的1/2處,且該對稱面處聲壓為0。聲壓幅值遠高于非諧振液位條件,且聲壓最大值分布表現(xiàn)為遠離聲源位置。

    圖3 兩種液位高度下的平均聲能密度

    2.3 不同功放電流條件下柱形液腔中的空化分布規(guī)律

    2.3.1 料腔聲場的聲能密度峰值特征

    為了驗證在諧振液位高度=75 mm和=90 mm下不同功率對超聲空化效應(yīng)空間分布特性的影響,功放電流取值范圍為40~80 mA,增量為10 mA(變幅桿浸入料液深度為10 mm,料腔半徑為50 mm,激勵頻率為20 kHz,液位高度為=75 mm和90 mm),兩種液位條件下料腔聲場中聲能密度峰值隨變幅桿輸入電流的變化特征如圖4所示。

    圖4 不同電流下的聲能密度峰值

    圖4中,在非諧振液位=75 mm,電流在40~80 mA的情況下,平均聲能密度的空間分布不變,聲能密度峰值整體較小,同時隨著電流增大,聲能密度峰值呈非線性增長。當柱形料腔處在諧振液位=90 mm時,電流在40~80 mA的情況下,平均聲能密度的空間分布不變,聲能輻射效果好,聲能的輻射便最大限度地得到了提高,并且是隨著電流的增大,聲能密度隨之增大。在諧振液位時,因柱形料腔中產(chǎn)生共振的現(xiàn)象明顯,共振效果好,導(dǎo)致空化效應(yīng)明顯,聲能密度較大;而非處于諧振液位時,共振效果不明顯,聲能密度峰值小,并且諧振液位是非諧振液位聲能密度峰值的10倍左右。

    為了更好地研究超聲破碎空化場的效果,將魯米諾發(fā)光實驗結(jié)果與仿真分析進行比對,以獲得超聲空化場的能量分布特性。

    魯米諾發(fā)光實驗中,變幅桿激勵頻率為20 kHz,功放電流范圍為40~200 mA,變幅桿浸入料液深度為10 mm,料腔半徑為50 mm。采用兩種液位條件:非諧振液位高度為75 mm;諧振液位高度為90 mm。

    2.3.2 非諧振液位料腔空化場

    柱形料腔中,在非諧振液位=75 mm時,隨功放電流變化,料腔場內(nèi)魯米諾發(fā)光圖像如圖5所示。當電流處在40~200 mA范圍時,空化區(qū)域不符合線性模態(tài)分布,此時空化區(qū)域呈現(xiàn)非線性分布。在此階段,隨著電流變化,柱形料腔中的空化效應(yīng)一直符合仿真的結(jié)果,聲能密度的大小隨之變化,料液中的物料在超聲處理下只在變幅桿端部附近能夠有效空化破碎,超聲能量的轉(zhuǎn)換效率比較低。非諧振液位時,聲壓最大值分布于變幅桿端部附近,超聲空化集中作用于變幅桿端部,一方面加劇了變幅桿端部表面的金屬腐蝕,另一方面使其附近料液的溫迅速明顯,大大降低了聲能的轉(zhuǎn)換效率和作用范圍,不利于實際應(yīng)用。實驗結(jié)果驗證了仿真得出的空化效應(yīng)的空間分布特性。

    圖5 非諧振液位H=75 mm時在不同電流下的SCL圖像

    2.3.3 諧振液位的空化場特征

    柱形料腔中,在諧振液位=90 mm時,隨功放電流變化,料腔場內(nèi)魯米諾發(fā)光圖像如圖6所示。當電流在40~80 mA范圍時,空化區(qū)域符合線性模態(tài)分布,在此階段,隨著電流變化,柱形料腔中的空化效應(yīng)一直符合仿真的結(jié)果,空化區(qū)域呈現(xiàn)模態(tài)分布,只有聲能密度的大小隨之變化,諧振液位時,聲壓最大值表現(xiàn)為遠離聲源位置,且超聲空化區(qū)域也明顯增大,聲能的轉(zhuǎn)換效率顯著提高,空化作用范圍有所擴大。當功放電流繼續(xù)增大,至80~200 mA時,空化區(qū)域的模態(tài)分布被破壞,呈現(xiàn)柱形束狀分布,隨著電流進一步增大,集束狀的發(fā)光區(qū)域擴展至料腔底部形成平鋪層狀。由此推斷,隨功放電流增大,聲源振幅的增大導(dǎo)致聲波的非線性傳播形成,出現(xiàn)有聲流現(xiàn)象,不再滿足小振幅聲模態(tài)規(guī)律,有待后續(xù)進一步研究分析,本文不再討論。

    2.3.4 空化域分析

    諧振液位=90 mm與非諧振液位=75 mm魯米諾發(fā)光面積占比百分數(shù)柱狀圖如圖7所示。

    圖6 諧振液位H=90 mm時在不同電流下的SCL圖像

    圖7 諧振與非諧振液位在不同電流下發(fā)光面積占比

    圖7中,柱形料腔在在非諧振液位=75 mm、電流在40~200 mA時,魯米諾發(fā)光面積占比百分數(shù)隨著電流的增大而增強,在諧振液位=90 mm,電流在40~200 mA時,隨著電流的增大,魯米諾發(fā)光面積占比百分數(shù)先增大后減小,此時的空化效應(yīng)效果明顯。并且圓柱形料腔聲場在諧振液位=90 mm時的發(fā)光面積百分比大于非諧振液位=75 mm時的發(fā)光面積百分比,即柱形料腔在諧振液位下,空化效應(yīng)遠強于非諧振液位,具有更高的能量轉(zhuǎn)化效率。

    3 結(jié)論

    本文基于有限元方法進行仿真,再應(yīng)用魯米諾聲致化學(xué)發(fā)光(SCL)的方法,對柱形料腔中的空化效應(yīng)進行半定量測量,研究了功率和液位高度對超聲空化效應(yīng)的空間分布特性的影響規(guī)律,得到以下結(jié)論:

    (1) 將實驗方法與有限元分析相結(jié)合,研究了定容料腔聲場的聲模態(tài)條件下的聲能分布特性。選取半徑為50mm的圓柱形料腔聲場,由實驗獲取的聲源頻率和振速作為聲源邊界條件,結(jié)合聲傳媒介的物料特性確定壁面邊界條件,應(yīng)用聲波導(dǎo)方程獲得該聲場條件下的諧振液位(二階)高度為90 mm,提供了定容料腔聲場內(nèi)超聲波聲能傳播及聲量轉(zhuǎn)化分析的理論基礎(chǔ)。

    (2) 對于定容料腔聲場而言,當容腔半徑一定時,不同液位高度中聲能量分布特性具有顯著差異。相同聲源條件下,非諧振液位聲場中平均聲能密度幅值有限,且聲能分布集中在聲源附近,作用區(qū)域有限,導(dǎo)致空化現(xiàn)象出現(xiàn)在變幅桿端部附近,且遠場空化效應(yīng)較弱;諧振液位條件下平均聲能密度幅值大大增強,且分布范圍擴大明顯,表現(xiàn)為空化效應(yīng)顯著增強,空化作用域的空間分布呈穩(wěn)定擴大趨勢,且具有明顯的遠場空化效應(yīng)。

    (3) 諧振液位條件下,功率大小(即功放電流大小)也是影響聲能密度幅值與空化效應(yīng)空間分布特性的重要因素。功率大小對聲能和空間分布特性影響顯著。在低功率條件下,空化效應(yīng)呈線性模態(tài)分布,超聲能量轉(zhuǎn)換更為高效;在高功率條件下,空化效應(yīng)受聲流影響呈非線性分布,空化區(qū)域呈柱形束狀分布,會導(dǎo)致液體溫升明顯,壁面邊界的空化腐蝕效果增強。

    本文研究了超聲空化效應(yīng)的空間分布特性,拓展了空化效應(yīng)理論的研究,為超聲破碎的實際應(yīng)用和化學(xué)反應(yīng)器的設(shè)計優(yōu)化提供了理論參考。

    [1] 何明霞, 夏大祥, 王康, 等. 低頻超聲激勵射流破碎過程模擬與實驗研究[J]. 天津大學(xué)學(xué)報(自然科學(xué)與工程技術(shù)版), 2013, 46(11): 1029-1033.

    HE Mingxia, XIA Daxiang, WANG Kang, et al. Jet breakup process simulation and experimental study of low-frequency ultrasonic excitation[J]. Journal of Tianjin University(Science and Technology), 2013, 46(11): 1029-1033.

    [2] 馮麗麗, 劉斌, 徐雪冬, 等. 液位對超聲破碎能量場的影響研究[J]. 農(nóng)產(chǎn)品加工, 2018(4): 15-17.

    FENG Lili, LIU Bin, XU Xuedong, et al. Research on the effect of liquid level on ultrasonic broken energy field[J]. Farm Products Processing, 2018(4): 15-17.

    [3] 嚴萬里, 陳曉明, 郭麗燕, 等. 超氧化物歧化酶活性測定的影響因素研究[J]. 生物學(xué)通報, 2011, 46(3): 50-53.

    YAN Wanli, CHEN Xiaoming, GUO Liyan, et al. Study on the influencing factors of superoxide dismutase activity determination [J]. Bulletin of Biology, 2011, 46(3): 50-53.

    [4] 郭璇, 楊艷玲, 李星, 等. 超聲空化及絮體破碎過程的模擬與試驗分析[J]. 中國環(huán)境科學(xué), 2015, 35(5): 1429-1435.

    GUO Xuan, YANG Yanling, LI Xing, et al. Simulation of ultrasound cavitation bubble and numerical/experimental analysis of flocs breakage process[J]. China Environmental Science, 2015, 35(5): 1429-1435.

    [5] 岑為, 夏先知, 劉月祥, 等. 基于超聲破碎測定聚丙烯球形催化劑強度的方法[J]. 石油煉制與化工, 2014, 45(8): 17-21.

    CEN Wei, XIA Xianzhi, LIU Yuexiang, et al. Determination of spherical catalyst strength for polypropylene by ultrasonication[J]. Petroleum Processing and Petrochemicals, 2014, 45(8): 17-21.

    [6] 張嬋, 鄭爽英. 超聲空化效應(yīng)及其應(yīng)用[J]. 水資源與水工程學(xué)報, 2009, 20(1): 136-138.

    ZHANG Chan, ZHENG Shuangying. Effect of ultrasonic cavitation and its application[J]. Journal of Water Resources and Water Engineering, 2009, 20(1): 136-138.

    [7] 應(yīng)崇福. 液體中的聲處理應(yīng)用和聲空化工程[J]. 應(yīng)用聲學(xué), 2006, 25(5): 261-264.

    YING Chongfu. Acoustic processing in liquids and acoustic cavitation engineering[J]. Applied Acoustics, 2006, 25(5): 261-264.

    [8] 崔方玲, 紀威. 超聲空化氣泡動力學(xué)仿真及其影響因素分析[J]. 農(nóng)業(yè)工程學(xué)報, 2013, 29(17): 24-29.

    CUI Fangling, JI Wei. Dynamic simulation of ultrasonic cavitation bubble and analysis of its influencing factors[J]. Transactions of the Chinese Society of Agricultural Engineering, 2013, 29(17): 24-29.

    [9] 劉麗艷, 劉芃宏, 楊洋, 等. 超聲頻率和液位對空化場的影響[J]. 天津大學(xué)學(xué)報(自然科學(xué)與工程技術(shù)版), 2016, 49(8): 802-808.

    LIU Liyan, LIU Penghong, YANG Yang, et al. Effect of ultrasonic frequency and liquid level on the cavitation field[J]. Journal of Tianjin University(Science and Technology), 2016, 49(8): 802-808.

    [10] 張艾萍, 楊釗, 夏榮濤, 等. 超聲頻率和管型對超聲空化效應(yīng)的影響研究[J]. 化工機械, 2016, 43(4): 517-521.

    ZANG Aiping, YANG Zhao, XIA Rongtao, et al. Influences of ultrasonic frequency and tube types on ultrasonic cavitation effect[J]. Chemical Engineering & Machinery, 2016, 43(4): 517-521.

    [11] 李爭彩, 林書玉. 超聲空化影響因素的數(shù)值模擬研究[J]. 陜西師范大學(xué)學(xué)報(自然科學(xué)版), 2008, 36(1): 38-42.

    LI Zhengcai, LIN Shuyu. Numerical simulation of the factors influencing ultrasonic cavitation[J]. Journal of Shaanxi Normal University(Natural Science Edition), 2008, 36(1): 38-42.

    [12] SON Y, LIM M, KHIM J, et al. Acoustic emission spectra and sonochemical activity in a 36 kHz sonoreactor[J]. Ultrasonics Sonochemistry, 2012, 19(1): 16-21.

    [13] KLíMA J, FRIAS-FERRER A, GONZáLEZ-GARCíA J, et al. Optimisation of 20 kHz sonoreactor geometry on the basis of numerical simulation of local ultrasonic intensity and qualitative comparison with experimental results[J]. Ultrasonics Sonochemistry, 2007, 14(1): 19-28.

    [14] ASHOKKUMAR M, LEE J, IIDA Y, et al. Spatial distribution of acoustic cavitation bubbles at different ultrasound frequencies[J]. Chemphyschem, 2010, 11(8): 1680-1684.

    [15] SON Y, LIM M, ASHOKKUMAR M, et al. Geometric optimization of sonoreactors for the enhancement of sonochemical activity[J]. Journal of Physical Chemistry C, 2011, 115(10): 4096-4103.

    [16] 楊日福, 丘泰球, 羅登林. 超臨界CO2流體中超聲空化閾值的研究[J]. 華南理工大學(xué)學(xué)報(自然科學(xué)版), 2005, 33(12): 100-104.

    YANG Rifu, QIU Taiqiu, LUO Denglin. Threshold investigation of ultrasonic cavitation in supercritical fluid CO2[J]. Journal of South China University of Technology(Natural Science Edition), 2005, 33(12): 100-104.

    [17] 林書玉, 鮮小軍. 功率超聲換能振動系統(tǒng)的優(yōu)化設(shè)計及其研究進展[J]. 陜西師范大學(xué)學(xué)報(自然科學(xué)版), 2014(6): 31-39.

    LIN Shuyu, XIAN Xiaojun. Progress and optimization design of high power piezoelectric ceramic ultrasonic vibrating system[J]. Journal of Shaanxi Normal University(Natural Science Edition), 2014(6): 31-39.

    [18] 劉浩東, 胡芳友, 李洪波, 等. 功率超聲技術(shù)的分類研究及應(yīng)用[J]. 電焊機, 2014, 44(12): 25-29.

    LIU Haodong, HU Fangyou, LI Hongbo, et al. Classification research and the application summarization of power ultrasonic technology[J]. Electric Welding Machine, 2014, 44(12): 25-29.

    [19] 潘巧生, 劉永斌, 賀良國, 等. 一種大振幅超聲變幅桿設(shè)計[J]. 振動與沖擊, 2014, 33(9): 1-5.

    PAN Qiaosheng, LIU Yongbin, HE Liangguo, et al. Design of an ultrasonic horn with high amplitude of longitudinal vibration[J]. Journal of Vibration and Shock, 2014, 33(9): 1-5.

    [20] 張穎, 劉曉宙, 龔秀芬. 兩個硬邊界附近的聲流的數(shù)值計算(英文)[J]. 南京大學(xué)學(xué)報(自然科學(xué)), 2012, 48(5): 565-571.

    ZHANG Ying, LIU Xiaozhou, GONG Xiufen. Numerical calculations of acoustic streaming near two rigid planes[J]. Journal of Nanjing University(Natural Sciences), 2012, 48(5): 565-571.

    [21] 劉丹曉, 周城光, 劉克. 諧振管形狀對熱聲發(fā)動機內(nèi)非線性聲場的影響[J]. 聲學(xué)學(xué)報, 2012(1): 36-44.

    LIU Danxiao, ZHOU Chengguang, LIU Ke. Influence of resonator shape on nonlinear acoustic field in a thermoacoustic engine[J]. Acta Acustica, 2012(1): 36-44.

    [22] 魏鑫, 吳勝舉, 許昊, 等. 槽壁邊界對超聲清洗聲場特性的影響[J]. 陜西師范大學(xué)學(xué)報(自然科學(xué)版), 2012, 40(3): 48-51.

    WEI Xin, WU Shengju, XU Hao, et al. Influence of tank wall boundary on characteristics of sound field in ultrasonic cleaning[J]. Journal of Shaanxi Normal University(Natural Science Edition), 2012, 40(3): 48-51.

    [23] 郭策, 祝錫晶, 王建青, 等. 超聲場下剛性界面附近潰滅空化氣泡的速度分析[J]. 物理學(xué)報, 2016, 65(4): 180-187.

    GUO Ce, ZHU Xijing, WANG Jianqing, et al. Velocity analysis for collapsing cavitation bubble near a rigid wall under an ultrasound field[J]. Acta Physica Sinica, 2016, 65(4): 180-187.

    Study on the spatial distribution of ultrasonic cavitation effect in cylindrical reactor

    YU Jian, LIU Bin, XU Xuedong, FENG Lili

    (College of Material and Mechanical Engineering, Beijing Technology and Business University, Beijing 100048, China)

    The formation and spatial distribution of ultrasonic cavitation in the cylindrical reactor of 20 kHz are studied. The distribution of the acoustic energy density in the cylindrical sound field is obtained by using Bessel function, and the finite element method is used for simulation analysis. By combining the acoustic measurement method and the sono-chemi luminescent (SCL) method, the theoretical analysis results are verified in the power ultrasonic experiment of 20 kHz. The results show that the ultrasonic cavitation effect appears in the end region of the horn when the cavity radius is=50 mm, the resonant liquid level of 20 kHz is=90mm and the current of Power Amplifier is <40 mA, in the case of 40 mA ≤ the current of Power Amplifier ≤ 80 mA, the cavitation effect is significantly enhanced; that the spatial distribution of the cavitation effect is consistent with the distribution of sound pressure in the field, the acoustic modes affect the cavitation effect to form the distribution characteristics of the far-field cavitation effect, and when the current of Power Amplifier is >80 mA, due to nonlinear factors, the cavitation effect is distributed in a cylindrical tail shape under the action of the acoustic flow at the resonant liquid level, and a tile-like distribution is formed at the boundary of the bottom wall; and that for the non-resonant liquid level of 75 mm, with the increase of power, the ultrasonic cavitation effect appears only in the end region of the horn and presents the local cavitation distribution characteristics.

    cylindrical rector; ultrasonic cavitation; power ultrasonic; resonance liquid level; sono chemi luminescent (SCL)

    O426

    A

    1000-3630(2020)-02-0127-07

    10.16300/j.cnki.1000-3630.2020.02.001

    2019-05-05;

    2019-07-12

    國家重點研發(fā)計劃資助(2016YFD0400305)、北京市自然科學(xué)基金-市教委聯(lián)合資助項目資助(KZ201810011017)

    俞劍(1995-), 男, 安徽蕪湖人, 碩士研究生, 研究方向為噪聲與振動控制。

    劉斌,E-mail: liubin@th.btbu.edu.cn

    猜你喜歡
    魯米諾柱形變幅
    聚(苯胺-魯米諾)/多壁碳納米管/金納米復(fù)合材料的合成及其發(fā)光性能研究
    基于超聲振動加工變幅桿的優(yōu)化設(shè)計
    N-烷基魯米諾的合成與光譜性質(zhì)
    河南化工(2020年11期)2020-12-10 05:39:06
    大端接圓柱桿的復(fù)合圓錐形變幅桿設(shè)計及應(yīng)用
    非柱形容器的壓力和壓強
    從“柱形凸透鏡成像”實驗感受體驗教學(xué)的魅力
    雙排斜槽縱-扭復(fù)合振動變幅桿仿真設(shè)計
    基于魯米諾-過氧化氫-CdTe的流動注射-化學(xué)發(fā)光體系及應(yīng)用于多巴胺的檢測
    超聲縱振動空心變幅桿的特性
    血催化魯米諾化學(xué)發(fā)光性能的研究
    精品熟女少妇八av免费久了| 亚洲av成人av| 日韩中字成人| av国产免费在线观看| 国产一级毛片七仙女欲春2| 极品教师在线视频| 欧美极品一区二区三区四区| 欧美色欧美亚洲另类二区| 一进一出抽搐动态| 一级黄片播放器| 一级av片app| 尤物成人国产欧美一区二区三区| 国产成人a区在线观看| 99久久精品一区二区三区| 国产伦一二天堂av在线观看| 又黄又爽又免费观看的视频| 亚洲精品一卡2卡三卡4卡5卡| 在线天堂最新版资源| 成人毛片a级毛片在线播放| 欧美日本亚洲视频在线播放| 欧美bdsm另类| 91在线精品国自产拍蜜月| 国产中年淑女户外野战色| 国产大屁股一区二区在线视频| 国产精品亚洲av一区麻豆| 久久久久久九九精品二区国产| 国产精品不卡视频一区二区 | 欧美bdsm另类| a级毛片免费高清观看在线播放| www.999成人在线观看| av在线天堂中文字幕| 日韩欧美精品免费久久 | 日韩人妻高清精品专区| 国产精品av视频在线免费观看| 欧美黑人欧美精品刺激| АⅤ资源中文在线天堂| 狠狠狠狠99中文字幕| 亚洲七黄色美女视频| 国产精品亚洲av一区麻豆| 麻豆成人午夜福利视频| 亚洲人与动物交配视频| 好男人在线观看高清免费视频| 一进一出抽搐动态| 一卡2卡三卡四卡精品乱码亚洲| 国产aⅴ精品一区二区三区波| 国产探花在线观看一区二区| 怎么达到女性高潮| 搡老熟女国产l中国老女人| 欧美激情久久久久久爽电影| 午夜两性在线视频| 国产爱豆传媒在线观看| 少妇人妻一区二区三区视频| 日韩欧美一区二区三区在线观看| www.999成人在线观看| 少妇裸体淫交视频免费看高清| 十八禁人妻一区二区| 亚洲人成伊人成综合网2020| 国产成人啪精品午夜网站| 国产成人欧美在线观看| 亚洲精品成人久久久久久| 丰满的人妻完整版| www.999成人在线观看| 一a级毛片在线观看| 国产高清激情床上av| 一进一出好大好爽视频| 久久久久久国产a免费观看| 99精品在免费线老司机午夜| 亚洲成av人片在线播放无| 伦理电影大哥的女人| 悠悠久久av| 精品久久久久久久久亚洲 | 午夜福利欧美成人| 欧美黑人欧美精品刺激| 男女之事视频高清在线观看| 午夜精品在线福利| 如何舔出高潮| 小蜜桃在线观看免费完整版高清| 亚洲国产日韩欧美精品在线观看| www日本黄色视频网| 91在线精品国自产拍蜜月| 成熟少妇高潮喷水视频| 一区二区三区激情视频| 欧美性猛交黑人性爽| 在线免费观看不下载黄p国产 | 免费看a级黄色片| 美女xxoo啪啪120秒动态图 | 久久精品影院6| 午夜精品久久久久久毛片777| 欧美一级a爱片免费观看看| 1000部很黄的大片| 99国产综合亚洲精品| 国产视频内射| 欧美激情国产日韩精品一区| 国产免费av片在线观看野外av| 在线观看午夜福利视频| 日韩大尺度精品在线看网址| 两个人的视频大全免费| 国产精品女同一区二区软件 | 国产探花极品一区二区| www.色视频.com| 亚洲精品久久国产高清桃花| 中文字幕久久专区| 国产男靠女视频免费网站| 在线观看美女被高潮喷水网站 | 国产精品影院久久| 中文字幕熟女人妻在线| 亚洲一区高清亚洲精品| 亚洲欧美日韩东京热| 欧美bdsm另类| 日韩亚洲欧美综合| 搡老熟女国产l中国老女人| 一级a爱片免费观看的视频| 美女xxoo啪啪120秒动态图 | 波多野结衣高清作品| 成人性生交大片免费视频hd| 特级一级黄色大片| 免费电影在线观看免费观看| 亚洲精品影视一区二区三区av| a在线观看视频网站| 老熟妇仑乱视频hdxx| ponron亚洲| 久久这里只有精品中国| 日日干狠狠操夜夜爽| 国产精品伦人一区二区| 免费一级毛片在线播放高清视频| 国产高潮美女av| 一二三四社区在线视频社区8| 欧美乱妇无乱码| 真实男女啪啪啪动态图| 亚洲专区国产一区二区| 不卡一级毛片| 久久精品久久久久久噜噜老黄 | 在线播放无遮挡| 在线天堂最新版资源| 美女免费视频网站| 国产精品日韩av在线免费观看| 亚洲五月婷婷丁香| 亚洲男人的天堂狠狠| 99久久成人亚洲精品观看| 亚洲av中文字字幕乱码综合| 久久精品人妻少妇| 国产精品精品国产色婷婷| 蜜桃亚洲精品一区二区三区| 国产精品1区2区在线观看.| 欧美黄色片欧美黄色片| 午夜a级毛片| 日韩高清综合在线| 日韩欧美在线乱码| 国产欧美日韩一区二区三| 麻豆成人午夜福利视频| 亚洲一区二区三区不卡视频| 久久精品国产清高在天天线| 美女 人体艺术 gogo| 国产精品99久久久久久久久| 日韩欧美三级三区| 欧美3d第一页| 日韩中文字幕欧美一区二区| 欧美日韩乱码在线| 男女那种视频在线观看| 人人妻人人澡欧美一区二区| 久久久久九九精品影院| 久久久久免费精品人妻一区二区| 久久中文看片网| 可以在线观看的亚洲视频| 久久久久久久久久成人| av欧美777| 久久久色成人| 在线播放无遮挡| 一区福利在线观看| 狂野欧美白嫩少妇大欣赏| 搡老熟女国产l中国老女人| 久久精品国产亚洲av香蕉五月| 国产亚洲av嫩草精品影院| 午夜日韩欧美国产| 亚洲国产精品999在线| 最近最新免费中文字幕在线| 91九色精品人成在线观看| 久久精品国产亚洲av香蕉五月| 99在线视频只有这里精品首页| 91在线观看av| 久久婷婷人人爽人人干人人爱| 国产人妻一区二区三区在| 国产伦一二天堂av在线观看| 99热只有精品国产| 成人亚洲精品av一区二区| 国产精品一区二区性色av| 日本免费a在线| 91麻豆av在线| 午夜精品一区二区三区免费看| 舔av片在线| 国产欧美日韩一区二区精品| 亚洲av不卡在线观看| 天堂av国产一区二区熟女人妻| 日本黄色视频三级网站网址| 一a级毛片在线观看| 国产一区二区亚洲精品在线观看| 亚洲av熟女| 中文字幕精品亚洲无线码一区| 国产 一区 欧美 日韩| 国产精品国产高清国产av| 精品一区二区三区人妻视频| 在线观看一区二区三区| а√天堂www在线а√下载| 一进一出抽搐动态| 高潮久久久久久久久久久不卡| av天堂中文字幕网| 亚洲精品成人久久久久久| 国产一区二区三区在线臀色熟女| 真人一进一出gif抽搐免费| 又黄又爽又刺激的免费视频.| 一个人看的www免费观看视频| www日本黄色视频网| 观看美女的网站| 国产在视频线在精品| 亚洲av成人不卡在线观看播放网| 嫁个100分男人电影在线观看| 在线国产一区二区在线| 美女黄网站色视频| 亚洲avbb在线观看| 日本五十路高清| 在线观看午夜福利视频| 日本熟妇午夜| 十八禁网站免费在线| 中出人妻视频一区二区| 波多野结衣高清无吗| 国模一区二区三区四区视频| 国产免费男女视频| 有码 亚洲区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人久久性| 91在线精品国自产拍蜜月| 亚洲在线观看片| 日韩成人在线观看一区二区三区| 国产美女午夜福利| 亚洲精品在线美女| 久久午夜福利片| 国产色婷婷99| 久久久成人免费电影| 99热6这里只有精品| 亚洲色图av天堂| 女生性感内裤真人,穿戴方法视频| 蜜桃久久精品国产亚洲av| 此物有八面人人有两片| 男人的好看免费观看在线视频| 好看av亚洲va欧美ⅴa在| 美女高潮喷水抽搐中文字幕| 久久草成人影院| 国产午夜福利久久久久久| 亚洲七黄色美女视频| 可以在线观看毛片的网站| 热99在线观看视频| 成年女人毛片免费观看观看9| 亚洲精品在线美女| 成人无遮挡网站| 日日摸夜夜添夜夜添av毛片 | 国产一区二区亚洲精品在线观看| 欧美一级a爱片免费观看看| 最新中文字幕久久久久| 成年女人毛片免费观看观看9| 一区二区三区高清视频在线| 给我免费播放毛片高清在线观看| 国产白丝娇喘喷水9色精品| 成人特级av手机在线观看| 国产久久久一区二区三区| 午夜影院日韩av| 久久午夜亚洲精品久久| 免费av不卡在线播放| 听说在线观看完整版免费高清| 国产高清三级在线| 欧美午夜高清在线| 国产成人啪精品午夜网站| 人人妻人人澡欧美一区二区| 三级毛片av免费| 熟女电影av网| 欧美成人一区二区免费高清观看| 国产午夜精品论理片| 亚洲18禁久久av| 亚洲不卡免费看| 自拍偷自拍亚洲精品老妇| 精品欧美国产一区二区三| 久久久久久久久久成人| 国产精品精品国产色婷婷| 亚洲三级黄色毛片| 99热只有精品国产| 乱码一卡2卡4卡精品| 午夜影院日韩av| 国产淫片久久久久久久久 | 深爱激情五月婷婷| 看片在线看免费视频| 亚洲av.av天堂| 久久精品国产99精品国产亚洲性色| 国内少妇人妻偷人精品xxx网站| 久久午夜福利片| 两性午夜刺激爽爽歪歪视频在线观看| 久久久精品欧美日韩精品| 成人美女网站在线观看视频| 宅男免费午夜| 午夜久久久久精精品| 3wmmmm亚洲av在线观看| 免费在线观看亚洲国产| 熟女电影av网| 性插视频无遮挡在线免费观看| 国产伦在线观看视频一区| 麻豆国产97在线/欧美| 男女那种视频在线观看| 香蕉av资源在线| 成人亚洲精品av一区二区| 最好的美女福利视频网| xxxwww97欧美| 午夜免费成人在线视频| aaaaa片日本免费| 一进一出抽搐gif免费好疼| 禁无遮挡网站| 亚洲成人免费电影在线观看| 最后的刺客免费高清国语| 国产精品国产高清国产av| 少妇裸体淫交视频免费看高清| 亚洲av中文字字幕乱码综合| 亚洲国产欧美人成| 亚洲人成网站高清观看| 亚洲自偷自拍三级| 久久久精品大字幕| 色吧在线观看| 国产又黄又爽又无遮挡在线| 亚洲av.av天堂| 亚洲avbb在线观看| 午夜a级毛片| 国产精品久久久久久久久免 | 欧美色视频一区免费| 精品不卡国产一区二区三区| 亚洲精品久久国产高清桃花| 国产老妇女一区| 亚洲精品影视一区二区三区av| 日韩精品中文字幕看吧| 在线观看美女被高潮喷水网站 | 床上黄色一级片| 国产精品一区二区免费欧美| 国产高清视频在线观看网站| 中文字幕久久专区| 国产精品三级大全| 真人做人爱边吃奶动态| 久久久久久久午夜电影| 麻豆一二三区av精品| 精品久久久久久久久av| 波多野结衣巨乳人妻| 精品久久久久久久久av| 波多野结衣巨乳人妻| 麻豆国产97在线/欧美| 看免费av毛片| 亚洲内射少妇av| 怎么达到女性高潮| 国产一区二区激情短视频| 免费观看精品视频网站| 少妇的逼好多水| 99热这里只有是精品50| АⅤ资源中文在线天堂| 真人做人爱边吃奶动态| 嫩草影院精品99| 国产亚洲精品综合一区在线观看| 色综合亚洲欧美另类图片| 在线看三级毛片| 国产精品一区二区三区四区免费观看 | 一区二区三区高清视频在线| 久久久久久久午夜电影| 欧美中文日本在线观看视频| 国产精品亚洲一级av第二区| 亚洲无线在线观看| 极品教师在线视频| 日本黄色视频三级网站网址| 乱人视频在线观看| 国内久久婷婷六月综合欲色啪| 国产极品精品免费视频能看的| 乱人视频在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲精品色激情综合| 又爽又黄无遮挡网站| 亚洲一区二区三区色噜噜| 99久久无色码亚洲精品果冻| 又紧又爽又黄一区二区| 久久九九热精品免费| 亚洲一区二区三区色噜噜| 国产精品嫩草影院av在线观看 | 久久国产乱子伦精品免费另类| 色视频www国产| 午夜免费成人在线视频| 老司机午夜福利在线观看视频| 丁香欧美五月| 欧美国产日韩亚洲一区| 色吧在线观看| 啦啦啦韩国在线观看视频| 国产精品三级大全| 中文资源天堂在线| 国产高清有码在线观看视频| 国语自产精品视频在线第100页| 夜夜爽天天搞| 亚洲成人免费电影在线观看| 麻豆成人午夜福利视频| 亚洲精品在线观看二区| 女人被狂操c到高潮| 美女免费视频网站| 国产精品永久免费网站| 免费电影在线观看免费观看| 99精品在免费线老司机午夜| 中文字幕熟女人妻在线| 日本熟妇午夜| 国模一区二区三区四区视频| 精品一区二区三区人妻视频| 中文字幕精品亚洲无线码一区| av在线观看视频网站免费| 国产单亲对白刺激| 黄色女人牲交| 亚洲在线观看片| 精品久久国产蜜桃| 国产伦精品一区二区三区四那| 亚洲av日韩精品久久久久久密| 国产精品不卡视频一区二区 | 亚洲 欧美 日韩 在线 免费| 可以在线观看的亚洲视频| 午夜福利在线在线| 亚洲精品456在线播放app | 日本 欧美在线| 欧美区成人在线视频| 成年人黄色毛片网站| 一区二区三区高清视频在线| 97超级碰碰碰精品色视频在线观看| 午夜福利18| 国产v大片淫在线免费观看| 精品国产三级普通话版| 99热这里只有是精品50| 综合色av麻豆| 久久久久国产精品人妻aⅴ院| 嫩草影院入口| 亚洲中文字幕日韩| 午夜福利在线观看吧| 免费看光身美女| 全区人妻精品视频| 日韩欧美精品v在线| 亚洲国产精品成人综合色| 亚洲av中文字字幕乱码综合| 亚洲熟妇熟女久久| 亚洲第一电影网av| av中文乱码字幕在线| 精品熟女少妇八av免费久了| 亚洲男人的天堂狠狠| 精品人妻偷拍中文字幕| 日韩亚洲欧美综合| 精品国产亚洲在线| 亚洲av免费高清在线观看| 不卡一级毛片| a在线观看视频网站| 精品福利观看| 欧美日韩乱码在线| 久久久久免费精品人妻一区二区| 97人妻精品一区二区三区麻豆| 桃红色精品国产亚洲av| 婷婷精品国产亚洲av| 成人午夜高清在线视频| 国产熟女xx| 最近视频中文字幕2019在线8| 啦啦啦韩国在线观看视频| 亚洲自偷自拍三级| 老女人水多毛片| 国产精品久久久久久精品电影| 99热这里只有精品一区| 欧美性感艳星| 伦理电影大哥的女人| 亚洲五月婷婷丁香| 国产三级中文精品| 人妻久久中文字幕网| 我的女老师完整版在线观看| 少妇的逼水好多| 俄罗斯特黄特色一大片| 免费av不卡在线播放| 久9热在线精品视频| 国产精品影院久久| 精品久久久久久成人av| 波多野结衣高清作品| 亚洲午夜理论影院| 亚洲精华国产精华精| 两个人的视频大全免费| 成年女人毛片免费观看观看9| 香蕉av资源在线| 国产精品野战在线观看| 亚洲欧美日韩卡通动漫| 午夜福利在线观看免费完整高清在 | 黄色丝袜av网址大全| 国产一区二区在线av高清观看| 韩国av一区二区三区四区| 91久久精品国产一区二区成人| 动漫黄色视频在线观看| 成年免费大片在线观看| 高清日韩中文字幕在线| 免费av观看视频| 啪啪无遮挡十八禁网站| 亚洲成人久久性| 性欧美人与动物交配| 91麻豆精品激情在线观看国产| 看片在线看免费视频| 欧美成人性av电影在线观看| 日本黄大片高清| 欧美丝袜亚洲另类 | 嫩草影院新地址| 亚洲成人免费电影在线观看| 欧美色欧美亚洲另类二区| 热99在线观看视频| a级一级毛片免费在线观看| 一区二区三区四区激情视频 | 在线免费观看不下载黄p国产 | 精品久久久久久久末码| 国产一级毛片七仙女欲春2| 美女高潮的动态| 色吧在线观看| 日日摸夜夜添夜夜添av毛片 | 老女人水多毛片| 噜噜噜噜噜久久久久久91| 伦理电影大哥的女人| av女优亚洲男人天堂| 少妇人妻一区二区三区视频| 成人美女网站在线观看视频| 在线观看美女被高潮喷水网站 | 久久亚洲真实| 亚洲 国产 在线| 欧洲精品卡2卡3卡4卡5卡区| 琪琪午夜伦伦电影理论片6080| 激情在线观看视频在线高清| 别揉我奶头~嗯~啊~动态视频| av中文乱码字幕在线| 级片在线观看| av在线天堂中文字幕| 国产亚洲精品久久久久久毛片| 在线看三级毛片| 欧美午夜高清在线| 中文资源天堂在线| 日本一二三区视频观看| 久久国产乱子伦精品免费另类| 又黄又爽又刺激的免费视频.| 一级毛片久久久久久久久女| 国产午夜福利久久久久久| 国产成人欧美在线观看| 亚洲国产欧洲综合997久久,| 99国产精品一区二区三区| 亚洲成av人片免费观看| 成人三级黄色视频| 高潮久久久久久久久久久不卡| 嫩草影院新地址| 国产免费一级a男人的天堂| 老司机午夜福利在线观看视频| 亚洲熟妇中文字幕五十中出| 亚洲国产精品合色在线| 国产真实乱freesex| 午夜a级毛片| 赤兔流量卡办理| 超碰av人人做人人爽久久| 久久久久国产精品人妻aⅴ院| 午夜福利免费观看在线| www.www免费av| 全区人妻精品视频| 久久久国产成人免费| 成年女人看的毛片在线观看| 久久亚洲真实| 全区人妻精品视频| 成人av在线播放网站| 午夜福利在线在线| 成年女人毛片免费观看观看9| 搡女人真爽免费视频火全软件 | 两性午夜刺激爽爽歪歪视频在线观看| 搡老妇女老女人老熟妇| 日韩精品中文字幕看吧| 色精品久久人妻99蜜桃| 精品一区二区三区视频在线观看免费| 成人性生交大片免费视频hd| 日本五十路高清| 婷婷丁香在线五月| 久久久成人免费电影| 夜夜躁狠狠躁天天躁| 丰满乱子伦码专区| 一个人看视频在线观看www免费| 99国产精品一区二区蜜桃av| 国产亚洲精品久久久久久毛片| av女优亚洲男人天堂| 婷婷精品国产亚洲av在线| 国产69精品久久久久777片| 久久精品久久久久久噜噜老黄 | 午夜精品久久久久久毛片777| 国产精品影院久久| 特大巨黑吊av在线直播| 丰满的人妻完整版| 可以在线观看毛片的网站| 黄色视频,在线免费观看| 国产毛片a区久久久久| 1024手机看黄色片| 日本 av在线| 9191精品国产免费久久| 亚洲av免费高清在线观看| 在线观看一区二区三区| 亚洲午夜理论影院| 亚洲成av人片免费观看| 最新中文字幕久久久久| 在线播放无遮挡| 精品欧美国产一区二区三| 热99在线观看视频| 国产午夜精品论理片| 99精品久久久久人妻精品| 一进一出好大好爽视频| 亚洲va日本ⅴa欧美va伊人久久| bbb黄色大片| 97超视频在线观看视频| 婷婷亚洲欧美| 免费在线观看亚洲国产| 亚洲国产欧洲综合997久久,| 国产精品98久久久久久宅男小说| 桃色一区二区三区在线观看| 波多野结衣高清作品| 免费av不卡在线播放| 久久精品国产99精品国产亚洲性色|