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

    雙模態(tài)顆粒聲凝并微觀(guān)行為的數(shù)值模擬

    2023-06-30 05:38:10楊娜娜凡鳳仙胡曉紅趙豪
    關(guān)鍵詞:數(shù)值模擬

    楊娜娜 凡鳳仙 胡曉紅 趙豪

    摘要:考慮顆粒間多種相互作用機(jī)理及顆粒碰撞后的凝并與反彈機(jī)制,建立聲凝并微觀(guān)動(dòng)力學(xué)模型,對(duì)大粒徑外加顆粒及兩個(gè)鄰近的亞微米細(xì)顆粒之間的聲凝并微觀(guān)行為進(jìn)行數(shù)值模擬研究。結(jié)果表明,雙模態(tài)顆粒聲凝并行為主要表現(xiàn)為外加顆粒僅與一個(gè)細(xì)顆粒發(fā)生凝并,以及外加顆粒先后與兩個(gè)細(xì)顆粒發(fā)生凝并,聲尾流效應(yīng)、重力沉降作用是雙模態(tài)顆粒聲凝并的重要機(jī)理。隨著外加顆粒初始位置與波節(jié)間距的增大,外加顆粒與較近細(xì)顆粒的平均凝并時(shí)間變化很小,兩者形成的團(tuán)聚體與較遠(yuǎn)細(xì)顆粒的平均凝并時(shí)間顯著縮短;外加顆粒距離波節(jié)很近時(shí),其仍能與較近的細(xì)顆粒發(fā)生凝并。對(duì)于外加顆粒與較近細(xì)顆粒以及兩者形成的團(tuán)聚體與較遠(yuǎn)細(xì)顆粒的凝并,隨著外加顆粒直徑的增加,平均凝并時(shí)間顯著縮短,凝并概率先增加而后趨于恒定。

    關(guān)鍵詞:細(xì)顆粒;外加顆粒;聲凝并;微觀(guān)行為;數(shù)值模擬

    中圖分類(lèi)號(hào):O 359??????????? 文獻(xiàn)標(biāo)志碼:A

    Numerical simulation on microscopic behavior of acoustic agglomeration of bimodal particles

    YANG Nana, FAN Fengxian, HU Xiaohong, ZHAO Hao

    (School ofEnergy and Power Engineering, University of Shanghaifor Science and Technology, Shanghai 200093, China)

    Abstract: A microscopic dynamic model for particle agglomeration was developed with simultaneous consideration of multiple particle interaction mechanisms as well as mechanisms of agglomeration or rebound upon inter-particle collision. The microscopic behaviors of acoustic agglomeration of a large- sizedadditionalparticleanditstwoneighboringsubmicron-sizedfineparticleswerenumerically investigated. The results show that the major acoustic agglomeration behaviors of the bimodal particles are characterized by the agglomeration between the additional particle and either fine particle and the successiveagglomerate between theaddition particleand thecloser andfartherfine particles. The gravitational effect and acoustic wake effect play an important role in the acoustic agglomeration of bimodal particles. When the distance between the initial position of the additional particle and wave node increases, small variation in the average agglomeration time between the additional particle and the closerfine particleisobserved, whereas theaverageagglomeration time between theagglomerateformed? by? the? said? particles? and? the? farther? fine? particle? decreases? obviously.? When? theadditionalparticle is initially close to the wave node, it is possible for the particle to agglomerate with its closerfine particle. Concerning the agglomeration between the additional particle and its closer fine particle aswell? as? the? agglomerate? formed? by? the? said? particles? and? the? farther? fine? particle,? the? averageagglomeration? time? increases? while? the? agglomeration? probability? increases? first? and? then? tends? to? aconstant value with the increasing additional particle diameter.

    Keywords: fine particle; additional particle; acoustic agglomeration; microscopic behavior;numerical simulation

    細(xì)顆粒( PM2.5)對(duì)人體健康和大氣環(huán)境帶來(lái)嚴(yán)重危害,高溫燃燒過(guò)程的排放是細(xì)顆粒的主要源頭。隨著環(huán)保標(biāo)準(zhǔn)的日趨嚴(yán)格,細(xì)顆粒排放控制日益成為能源與環(huán)境領(lǐng)域的研究熱點(diǎn)。鑒于細(xì)顆粒粒徑微小,常規(guī)除塵裝置對(duì)其脫除難以奏效,一些學(xué)者提出了聲凝并[1-2]、蒸汽相變凝結(jié)[3-4]等細(xì)顆粒預(yù)處理技術(shù),以促使其粒徑增大,從而提高常規(guī)除塵裝置的效率。其中,聲凝并利用外加聲場(chǎng)促進(jìn)顆粒凝并在一起,該技術(shù)因裝置簡(jiǎn)單、適用性強(qiáng)而備受研究者關(guān)注[1-2,5-9]。

    聲凝并效果與顆粒粒徑分布密切相關(guān),已有研究表明,對(duì)于粒徑集中在亞微米尺度的細(xì)顆粒,聲凝并效果仍不夠理想[6-9]。為了有效提高亞微米細(xì)顆粒的聲凝并效果,發(fā)展出了雙模態(tài)顆粒聲凝并技術(shù),即向煙氣中添加的微米尺度外加顆粒作為一個(gè)模態(tài)、煙氣中的細(xì)顆粒作為另一個(gè)模態(tài),利用兩者在聲場(chǎng)中動(dòng)力學(xué)行為的差異,從而強(qiáng)化聲凝并效果[6-9]。然而,目前對(duì)雙模態(tài)顆粒聲凝并微觀(guān)行為及機(jī)理的掌握仍很欠缺,難以對(duì)其效果優(yōu)化提供有效指導(dǎo)。

    為揭示雙模態(tài)顆粒聲凝并中涉及的顆粒間相互作用直至碰撞、凝并的微觀(guān)行為及機(jī)理,本文綜合考慮同向相互作用、重力沉降作用、聲尾流效應(yīng)、互散射效應(yīng)這4種顆粒間相互作用機(jī)理,以及顆粒碰撞后凝并與反彈的物理機(jī)制,建立雙模態(tài)顆粒凝并過(guò)程的微觀(guān)動(dòng)力學(xué)模型;利用數(shù)值模擬方法,研究顆粒凝并微觀(guān)行為隨外加顆粒和細(xì)顆粒中心連線(xiàn)與聲波波動(dòng)方向的初始夾角的變化特性,探討外加顆粒的初始位置和直徑對(duì)平均凝并時(shí)間和凝并概率的影響規(guī)律。本文研究可為細(xì)顆粒低成本、高效脫除提供科學(xué)依據(jù),對(duì)細(xì)顆粒超低排放具有重要意義。

    1模型與方法

    1.1物理模型

    針對(duì)水平駐波聲場(chǎng)中外加顆粒與鄰近的2個(gè)細(xì)顆粒開(kāi)展研究,顆粒的相對(duì)位置如圖1所示。圖中,顆粒1為外加顆粒,顆粒2和顆粒3為細(xì)顆粒; x 向?yàn)槁暡ú▌?dòng)的正方向, y 向?yàn)橹亓Ψ较颍?di 為顆粒i( i=1,2,3)的直徑;rij為顆粒i與 j 的中心距, rji=rij;θij為顆粒i指向顆粒 j 的向量與 x 向的夾角,θji=π+θij。

    1.2數(shù)學(xué)模型

    為建模的方便,作出如下簡(jiǎn)化假設(shè):

    a.認(rèn)為顆粒發(fā)生二維運(yùn)動(dòng),顆粒碰撞后僅能發(fā)生凝并和反彈;

    b.將顆粒和顆粒凝并后生成的團(tuán)聚體視為球形,認(rèn)為顆粒不帶電荷;

    c.為著重探討聲場(chǎng)的作用效果,不考慮布朗力的影響,認(rèn)為范德華力僅在顆粒碰撞接觸過(guò)程中發(fā)揮作用。

    1.2.1顆粒運(yùn)動(dòng)模型

    考慮顆粒所受曳力、重力與浮力,顆粒i的運(yùn)動(dòng)方程可寫(xiě)為[1, 10]

    1.3數(shù)值計(jì)算方法

    數(shù)值模擬流程如圖2所示。其中,顆粒間發(fā)生首次凝并后,繼續(xù)跟蹤團(tuán)聚體和未參與碰撞顆粒的運(yùn)動(dòng),直至再次凝并(即模擬顆粒總數(shù)Nt=1)或達(dá)到計(jì)算終止時(shí)間te=5 s 時(shí),計(jì)算停止。數(shù)值模擬時(shí),采用的時(shí)間步長(zhǎng)為?t=10?7 s ,利用單顆粒運(yùn)動(dòng)速度的解析解設(shè)置顆粒的初始速度[10]。

    2結(jié)果與討論

    目前尚缺少雙模態(tài)顆粒聲凝并微觀(guān)行為的實(shí)驗(yàn)結(jié)果,給模型的直接驗(yàn)證帶來(lái)困難。文獻(xiàn)[1,10]已對(duì)水平駐波聲場(chǎng)中2個(gè)等粒徑顆粒相互作用直至碰撞的微觀(guān)行為進(jìn)行了充分驗(yàn)證。本文將在3個(gè)顆粒中心共線(xiàn)的初始條件下,以1°為間隔在?180°~180°范圍改變顆粒中心連線(xiàn)與聲波波動(dòng)方向的初始夾角(簡(jiǎn)稱(chēng)初始夾角),進(jìn)行數(shù)值模擬。考慮到雙模態(tài)顆粒聲凝并的實(shí)際應(yīng)用中,外加顆粒采用CaO顆粒有望實(shí)現(xiàn)煙氣中細(xì)顆粒和 SO2的協(xié)同脫除,在數(shù)值模擬中,外加顆粒和細(xì)顆粒的物性參數(shù)和相互作用參數(shù),如密度、哈默克數(shù)、最小接觸距離、極限接觸應(yīng)力等依據(jù)CaO顆粒以及燃煤飛灰顆粒進(jìn)行設(shè)置[15]。表1給出了數(shù)值模擬采用的參數(shù)。需要說(shuō)明的是:由于細(xì)顆粒直徑在亞微米尺度,屬于零慣性顆粒,其對(duì)流場(chǎng)的擾動(dòng)可以忽略,且受外加顆粒影響發(fā)生運(yùn)動(dòng)的軌跡幾乎一致,并因其黏附性強(qiáng),與外加顆粒發(fā)生碰撞后,能夠最終與外加顆粒凝并在一起。因此,雙模態(tài)顆粒聲凝并行為對(duì)細(xì)顆粒直徑的變化不敏感。基于此,數(shù)值模擬中保持細(xì)顆粒直徑不變。

    2.1微觀(guān)行為分析

    圖3給出了不同初始夾角下的凝并時(shí)間。其中,凝并時(shí)間 ta 為發(fā)生凝并所需要的聲波作用時(shí)間。數(shù)值模擬中采用的外加顆粒直徑 d1=10μm ,初始位置 x0,1=0.1λ,λ=c/f=0.17 m,λ為聲波波長(zhǎng)。圖3的內(nèi)插圖給出了初始時(shí)刻顆粒的相對(duì)位置關(guān)系,初始時(shí)刻3個(gè)顆粒的中心共線(xiàn),顆粒2和顆粒3分別為初始時(shí)刻距離外加顆粒(顆粒1)較近和較遠(yuǎn)的細(xì)顆粒。需要說(shuō)明的是:在本文模擬中,大多數(shù)情況下顆粒首次碰撞后,經(jīng)歷數(shù)次“反彈?遠(yuǎn)離?靠近?碰撞”過(guò)程而發(fā)生凝并,少數(shù)情況下顆粒首次碰撞即引發(fā)凝并。前者發(fā)生顆粒首次碰撞速度較大時(shí),此時(shí)顆粒碰撞后的動(dòng)能能夠克服顆粒間黏附能,導(dǎo)致顆粒反彈。反彈后由于聲尾流與互散射作用,顆粒將再次靠近。由于碰撞過(guò)程中能量的耗散,碰撞速度降低,碰撞后反彈速度隨之降低。如此反復(fù),最終顆粒碰撞后的動(dòng)能小于顆粒間的黏附能,顆粒發(fā)生凝并。后者發(fā)生在首次碰撞速度較小時(shí),此時(shí)顆粒首次碰撞后的動(dòng)能無(wú)法克服顆粒間的黏附能,引起顆粒凝并在一起。

    通過(guò)分析可知,外加顆粒難以被聲波夾帶,與氣相之間存在較大的相對(duì)速度,從而引起強(qiáng)烈的聲尾流和互散射效應(yīng);細(xì)顆粒由于粒徑小,對(duì)流場(chǎng)的擾動(dòng)可以忽略??紤]到聲尾流與互散射效應(yīng)的作用距離,施加聲場(chǎng)后,外加顆粒的聲尾流效應(yīng)開(kāi)始發(fā)揮作用,而當(dāng)細(xì)顆粒接近外加顆粒表面時(shí),互散射效應(yīng)才發(fā)揮作用。聲尾流效應(yīng)表現(xiàn)為:顆粒中心連線(xiàn)與聲波方向接近平行時(shí),顆粒間發(fā)生吸引,且在平行時(shí)吸引作用最強(qiáng);顆粒中心連線(xiàn)與聲波方向接近垂直時(shí),顆粒間因發(fā)生排斥,且在垂直時(shí)排斥作用最強(qiáng)。

    圖3中的結(jié)果顯示,當(dāng)?12°≤θ0,12≤0°時(shí),顆粒不發(fā)生凝并。其原因是由于外加顆粒的重力沉降作用,細(xì)顆粒在其聲尾流效應(yīng)吸引區(qū)的停留時(shí)間較短,不足以使顆粒間發(fā)生碰撞。當(dāng)?25°≤θ0,12

    圖3中的結(jié)果顯示,當(dāng)?90°≤θ0,12

    圖3中的結(jié)果還表明,當(dāng)?180°≤θ0,12

    2.2初始位置的影響

    圖4給出了外加顆粒直徑 d1=10μm 時(shí)初始位置 x0,1對(duì)平均凝并時(shí)間 a 和凝并概率 p 的影響。數(shù)值模擬中,初始時(shí)刻顆粒的相對(duì)位置關(guān)系與圖3相同。 T 為聲波周期, T=1/f=0.5 ms;平均凝并時(shí)間定義為能夠發(fā)生凝并的初始夾角下顆粒凝并時(shí)間的平均值;凝并概率定義為能夠發(fā)生凝并的初始夾角個(gè)數(shù)與總的初始夾角個(gè)數(shù)的比值。需要說(shuō)明的是:本文針對(duì)顆粒微觀(guān)凝并行為定義了凝并概率,其不同于顆粒凝并宏觀(guān)效果建模中通常采用的凝并率,凝并率的定義是單位時(shí)間、單位體積內(nèi)顆粒凝并次數(shù)[16]。為確保計(jì)算結(jié)果獨(dú)立于夾角個(gè)數(shù)的選擇,在不同的夾角個(gè)數(shù)下進(jìn)行計(jì)算,發(fā)現(xiàn)當(dāng)夾角個(gè)數(shù)達(dá)到360時(shí),繼續(xù)增加夾角個(gè)數(shù),凝并時(shí)間與凝并概率的計(jì)算結(jié)果幾乎不發(fā)生改變,這表明選擇360個(gè)夾角進(jìn)行計(jì)算是合理的。

    由圖4可見(jiàn),外加顆粒與初始距離較近的細(xì)顆粒的平均凝并時(shí)間對(duì)初始位置的變化不敏感。其原因是:當(dāng)初始位置靠近波節(jié)點(diǎn)(聲波波動(dòng)速度始終為0的點(diǎn),如 x =0位置)時(shí),聲波引起的氣體介質(zhì)振動(dòng)較弱,因而外加顆粒的聲尾流效應(yīng)較弱,細(xì)顆粒處于外加顆粒聲尾流效應(yīng)的吸引區(qū)時(shí),向著外加顆粒運(yùn)動(dòng)的速度較低。初始位置遠(yuǎn)離波節(jié)點(diǎn)時(shí),雖然外加顆粒的聲尾流效應(yīng)更強(qiáng),但此時(shí)細(xì)顆粒的位移振幅增大,以致細(xì)顆粒受聲波夾帶往復(fù)運(yùn)動(dòng)過(guò)程中交替出現(xiàn)在外加顆粒尾流效應(yīng)的吸引區(qū)與排斥區(qū),降低了細(xì)顆粒向外加顆??拷乃俣取S蓤D4還可以看出,隨著初始位置遠(yuǎn)離波節(jié)點(diǎn),外加顆粒和初始距離較近的細(xì)顆粒形成的團(tuán)聚體與初始距離較遠(yuǎn)的細(xì)顆粒間的平均凝并時(shí)間縮短,這是由于初始位置越靠近波腹點(diǎn),初始距離較遠(yuǎn)的細(xì)顆粒受團(tuán)聚體聲尾流效應(yīng)的吸引作用越強(qiáng)的緣故。

    圖4中的結(jié)果還表明,外加顆粒與初始距離較近的細(xì)顆粒的凝并概率最大,并且在初始位置較靠近波節(jié)點(diǎn)時(shí)也能發(fā)生凝并,形成團(tuán)聚體;該團(tuán)聚體與初始距離較遠(yuǎn)的細(xì)顆粒的凝并概率居中;外加顆粒與初始距離較遠(yuǎn)的細(xì)顆粒的凝并概率最小,且初始位置距離波節(jié)點(diǎn)較近時(shí),無(wú)法發(fā)生凝并。其原因可解釋為:外加顆粒的聲尾流效應(yīng)對(duì)初始距離較近的細(xì)顆粒的吸引更強(qiáng),因而兩者發(fā)生凝并的概率最大;兩者形成團(tuán)聚體后,團(tuán)聚體的聲尾流效應(yīng)對(duì)初始距離較遠(yuǎn)的細(xì)顆粒發(fā)揮作用,相應(yīng)的凝并概率取決于該細(xì)顆粒在團(tuán)聚體的聲尾流效應(yīng)吸引區(qū)的停留時(shí)間、聲尾流效應(yīng)的強(qiáng)度,以及團(tuán)聚體與細(xì)顆粒接近時(shí)互散射效應(yīng)的作用效果,因此,凝并概率明顯低于外加顆粒與初始距離較近的細(xì)顆粒之間的凝并概率。外加顆粒與初始距離較遠(yuǎn)的細(xì)顆粒發(fā)生碰撞和凝并是建立在初始距離較近的細(xì)顆粒受外加顆粒聲尾流效應(yīng)的吸引而接近外加顆粒表面時(shí),互散射效應(yīng)使得細(xì)顆粒未能到達(dá)外加顆粒表面的基礎(chǔ)之上,而互散射效應(yīng)的作用距離短,作用效果受到顆粒中心連線(xiàn)與聲波波動(dòng)方向偏離程度的影響,相比于聲尾流效應(yīng)在顆粒碰撞中發(fā)揮的主導(dǎo)作用,互散射效應(yīng)起到輔助作用,因而外加顆粒僅與初始距離較遠(yuǎn)的細(xì)顆粒的凝并概率最低。此外,初始位置距離波節(jié)點(diǎn)較近時(shí),聲尾流的吸引效應(yīng)較弱,初始距離較遠(yuǎn)的細(xì)顆粒無(wú)法到達(dá)外加顆粒表面。凝并時(shí)間短、凝并概率大,意味著單位時(shí)間內(nèi)有更多細(xì)顆粒被外加顆粒捕集,細(xì)顆粒凝并的宏觀(guān)效果好?;诖耍梢圆捎梅蔷鶆蚣恿?,在凝并效果差的位置多添加外加顆粒,從而強(qiáng)化細(xì)顆粒凝并長(zhǎng)大效果。

    2.3外加顆粒直徑的影響

    圖5給出了外加顆粒直徑對(duì)平均凝并時(shí)間和凝并概率的影響。數(shù)值模擬采用的外加顆粒初始位置 x0,1=0.1λ。數(shù)值模擬中,初始時(shí)刻顆粒的相對(duì)位置關(guān)系與圖3相同。由圖5可以看出,顆粒凝并行為除呈現(xiàn)出圖4中的4種情況外,還存在外加顆粒先與初始距離較遠(yuǎn)的細(xì)顆粒凝并,形成團(tuán)聚體,之后團(tuán)聚體與初始距離較近的細(xì)顆粒發(fā)生凝并的現(xiàn)象。數(shù)值模擬顯示,這種現(xiàn)象發(fā)生在?180°≤θ0,12

    圖5中的結(jié)果還表明,隨著外加顆粒直徑的增大,外加顆粒與初始距離較近的細(xì)顆粒以及兩者形成的團(tuán)聚體與初始距離較遠(yuǎn)的細(xì)顆粒的凝并概率呈現(xiàn)先增加而后趨于恒定的特性。這是由于外加顆粒直徑越大,其聲尾流效應(yīng)的吸引作用越強(qiáng),外加顆?;驁F(tuán)聚體凝并細(xì)顆粒的能力增強(qiáng);然而,外加顆粒粒徑增大到一定程度,由于聲尾流效應(yīng)的吸引作用已足夠引起細(xì)顆粒處于吸引區(qū)時(shí)發(fā)生凝并,繼續(xù)增大外加顆粒直徑,凝并概率變化很小。由于外加顆粒與初始距離較遠(yuǎn)的細(xì)顆粒以及兩者形成的團(tuán)聚體與初始距離較近的細(xì)顆粒的凝并行為與互散射效應(yīng)有關(guān),互散射效應(yīng)在顆粒凝并中發(fā)揮輔助作用,因此,凝并概率很低。實(shí)際應(yīng)用中,可適當(dāng)增加外加顆粒直徑,以縮短凝并時(shí)間,提高凝并效率[17-18],從而提高細(xì)顆粒凝并的宏觀(guān)效果。此外,聲場(chǎng)條件也會(huì)影響凝并時(shí)間和凝并效率,后續(xù)的研究中可以對(duì)顆粒凝并動(dòng)力學(xué)模型進(jìn)行無(wú)量綱化,從而給出微觀(guān)凝并行為受無(wú)量綱參數(shù)的影響,以深層次揭示雙模態(tài)顆粒聲凝并的行為規(guī)律。

    3結(jié)論

    同時(shí)考慮多種顆粒間相互作用機(jī)理,以及顆粒碰撞后凝并與反彈的機(jī)制,建立聲凝并微觀(guān)動(dòng)力學(xué)模型,對(duì)雙模態(tài)顆粒聲凝并微觀(guān)行為受初始夾角、初始位置和直徑的影響進(jìn)行數(shù)值模擬研究,得到以下結(jié)論:

    a.雙模態(tài)顆粒凝并微觀(guān)行為主要表現(xiàn)為外加顆粒僅與初始距離較近的細(xì)顆粒凝并、外加顆粒僅與初始距離較遠(yuǎn)的細(xì)顆粒凝并、外加顆粒先后與初始距離較近和較遠(yuǎn)的細(xì)顆粒凝并,聲尾流效應(yīng)、重力沉降作用在雙模態(tài)顆粒聲凝并中發(fā)揮重要作用。

    b.隨著外加顆粒初始位置距波節(jié)點(diǎn)的距離增大,其與初始距離較近的細(xì)顆粒的平均凝并時(shí)間的變化很小,兩者形成的團(tuán)聚體與初始距離較遠(yuǎn)的細(xì)顆粒的平均凝并時(shí)間顯著縮短;外加顆粒與初始距離較近的細(xì)顆粒的凝并概率最大,且外加顆粒初始位置距離波節(jié)點(diǎn)很近時(shí),仍能發(fā)生凝并。

    c.對(duì)于在雙模態(tài)顆粒聲凝并中占主要地位的外加顆粒與初始位置較近的細(xì)顆粒凝并以及兩者形成的團(tuán)聚體與初始位置較遠(yuǎn)的細(xì)顆粒的凝并,隨著外加顆粒直徑的增加,平均凝并時(shí)間顯著縮短,凝并概率呈現(xiàn)先增加而后趨于恒定的特性。

    參考文獻(xiàn):

    [1] FAN F X, XU X, ZHANG S H, et al. Modeling of particle interactiondynamicsinstandingwaveacousticfield[J]. AerosolScienceandTechnology, 2019, 53(10):1204–1216.

    [2] WU Z H, FAN F X, YAN J P, et al. An adaptable direct simulationMonteCarlomethodforsimulatingacoustic agglomeration of solid particles[J]. Chemical Engineering Science, 2022, 249:117298.

    [3] FANF X, ZHANGSH,PENG ZB,etal. Numerical investigation of heterogeneous nucleation of water vapour on PM10 for particulate abatement[J]. The Canadian Journal of Chemical Engineering, 2019, 97(4):930–939.

    [4] FAN F X, ZHANG S H, WANG W Y, et al. Numerical investigation of PM2.5 size enlargement by heterogeneous condensation for particulate abatement[J]. ProcessSafety and Environmental Protection, 2019, 125:197–206.

    [5] KILIKEVI?IEN?K,KA?IANAUSKASR,KILIKEVI- ?IUSA,etal. Experimentalinvestigationofacoustic agglomeration of diesel engine exhaust particles using new createdacousticchamber[J]. PowderTechnology, 2020,360:421–429.

    [6] ZHANG G X, ZHOU T T, ZHANG L L, et al. Improving acousticagglomerationefficiencyofcoal-firedfly-ash particlesbyadditionofliquidbinders[J]. Chemical Engineering Journal, 2018, 334:891–899.

    [7] YAN J P, LIN Q, ZHAO S H, et al. Effect of seed nucleicombined with acoustic field on fine particles removal[J]. Powder Technology, 2018, 340:8–16.

    [8] ZHANG G X, ZHANG L L, WANG J, et al. Improving acoustic agglomeration efficiency by addition of sprayed liquiddroplets[J]. PowderTechnology, 2017, 317:181–188.

    [9] YAN J P, CHEN L Q, LI Z. Removal of fine particles from coalcombustioninthecombinedeffectofacoustic agglomerationandseeddropletswithwettingagent[J]. Fuel, 2016, 165:316–323.

    [10]屈廣寧, 凡鳳仙, 張斯宏, 等.駐波聲場(chǎng)中單分散細(xì)顆粒的相互作用特性[J].物理學(xué)報(bào), 2020, 69(6):064704.

    [11] ZHANGGX,LIUJZ,WANGJ,etal. Numericalsimulationofacousticwakeeffectinacoustic agglomerationunderOseenflowcondition[J]. Chinese Science Bulletin, 2012, 57(19):2404–2412.

    [12] SONG L M. Modeling of acoustic agglomeration of fine aerosolparticles[D]. StateCollege: PennsylvaniaState University, 1990.

    [13] SGROTTOL,SOMMERFELDM. Influenceof inter-particlecollisionsandagglomerationoncyclone performanceandcollectionefficiency[J]. TheCanadian Journal of Chemical Engineering, 2019, 97(2):511–522.

    [14] CROWE C T, SCHWARZKOPF J D, SOMMERFELD M, et al. Multiphase flows with droplets and particles[M].2nd ed. New York: CRC Press, 2011.

    [15] MARSHALLJS,LISQ. Adhesiveparticleflow: a discrete-elementapproach[M]. NewYork: Cambridge University Press, 2014.

    [16] SHENGCD,SHENXL. Modellingofacoustic agglomeration processes using the direct simulation Monte Carlo method[J]. Journal of Aerosol Science, 2006, 37(1):16–36.

    [17]蔣瑜, 賈楠, 蘇明旭.基于改進(jìn)差分進(jìn)化算法的超聲衰減譜反演計(jì)算[J].上海理工大學(xué)學(xué)報(bào) , 2020, 42(4):333–338.

    [18]凡鳳仙, 張斯宏, 白鵬博, 等.基于聲凝并的 PM2.5脫除技術(shù)研究進(jìn)展(Ⅱ):聲凝并與其他機(jī)制聯(lián)合作用[J].能源研究與信息, 2017, 33(4):205–210.

    (編輯:石瑛)

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場(chǎng)信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測(cè)
    鉆孔灌注樁樁底沉渣對(duì)樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測(cè)試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無(wú)支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開(kāi)切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    小说图片视频综合网站| 国产精品1区2区在线观看.| 国产高清不卡午夜福利| 久久久久久久久久久免费av| 亚洲欧美日韩无卡精品| 欧美精品国产亚洲| 亚洲国产色片| 午夜精品一区二区三区免费看| 国产高清不卡午夜福利| 九九久久精品国产亚洲av麻豆| 青青草视频在线视频观看| 午夜爱爱视频在线播放| 啦啦啦韩国在线观看视频| 精品酒店卫生间| 精品人妻视频免费看| 中文字幕熟女人妻在线| 国产麻豆成人av免费视频| 中文字幕人妻熟人妻熟丝袜美| 色综合色国产| 看非洲黑人一级黄片| 91久久精品国产一区二区三区| 男插女下体视频免费在线播放| 成年av动漫网址| 美女被艹到高潮喷水动态| 最近中文字幕2019免费版| 成人国产麻豆网| 欧美性猛交黑人性爽| 亚洲国产色片| 一级爰片在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av中文av极速乱| 亚洲av.av天堂| 亚洲av熟女| 少妇人妻精品综合一区二区| 尾随美女入室| 在线观看av片永久免费下载| av免费在线看不卡| 国产精华一区二区三区| 欧美另类亚洲清纯唯美| 黄色配什么色好看| 91av网一区二区| 日本免费a在线| 亚洲成人久久爱视频| 高清在线视频一区二区三区 | 午夜日本视频在线| 国产黄a三级三级三级人| 亚洲国产最新在线播放| 欧美高清性xxxxhd video| 少妇高潮的动态图| 国产不卡一卡二| 亚洲最大成人av| 99热6这里只有精品| 国产精品麻豆人妻色哟哟久久 | av在线观看视频网站免费| 美女黄网站色视频| 免费观看的影片在线观看| 午夜福利成人在线免费观看| 国产伦一二天堂av在线观看| 免费播放大片免费观看视频在线观看 | 久久6这里有精品| 国产美女午夜福利| 精品国内亚洲2022精品成人| 国产亚洲5aaaaa淫片| 又粗又硬又长又爽又黄的视频| 一个人观看的视频www高清免费观看| 最近2019中文字幕mv第一页| 精品久久久久久久久久久久久| 尾随美女入室| 国内揄拍国产精品人妻在线| 国产成人a∨麻豆精品| 欧美成人午夜免费资源| 亚洲成色77777| 少妇人妻一区二区三区视频| 乱人视频在线观看| 欧美成人免费av一区二区三区| 亚洲最大成人中文| 七月丁香在线播放| 亚洲精品影视一区二区三区av| 精品99又大又爽又粗少妇毛片| 国产激情偷乱视频一区二区| 夜夜看夜夜爽夜夜摸| 日韩三级伦理在线观看| 亚洲高清免费不卡视频| 日本一本二区三区精品| 两个人的视频大全免费| 免费看美女性在线毛片视频| 国产黄片视频在线免费观看| 国模一区二区三区四区视频| videossex国产| 久久午夜福利片| 亚洲va在线va天堂va国产| 久久久久九九精品影院| 非洲黑人性xxxx精品又粗又长| 国产乱来视频区| 内地一区二区视频在线| 国产精品野战在线观看| 亚洲成av人片在线播放无| 亚洲,欧美,日韩| 内射极品少妇av片p| 国产精品1区2区在线观看.| 免费观看在线日韩| 日韩精品青青久久久久久| 青春草亚洲视频在线观看| 久久久国产成人免费| 乱码一卡2卡4卡精品| 国产免费福利视频在线观看| 变态另类丝袜制服| 国产麻豆成人av免费视频| 国产亚洲精品av在线| 婷婷色麻豆天堂久久 | a级一级毛片免费在线观看| 亚洲精品乱码久久久久久按摩| 国产成人福利小说| 狂野欧美白嫩少妇大欣赏| 久久亚洲国产成人精品v| 久久欧美精品欧美久久欧美| 岛国在线免费视频观看| 国产精品永久免费网站| 1000部很黄的大片| 国产亚洲最大av| 熟女人妻精品中文字幕| www.av在线官网国产| 少妇人妻一区二区三区视频| 国产男人的电影天堂91| 成人漫画全彩无遮挡| 我要搜黄色片| 国产乱来视频区| 黄色配什么色好看| 亚洲av.av天堂| 麻豆成人午夜福利视频| 国产成人午夜福利电影在线观看| 99久久精品热视频| 日韩亚洲欧美综合| 国产亚洲精品av在线| 国产亚洲5aaaaa淫片| 欧美激情在线99| 如何舔出高潮| 精品国产露脸久久av麻豆 | 成人毛片a级毛片在线播放| 久久久精品94久久精品| 亚洲国产精品成人久久小说| 内射极品少妇av片p| 国产真实伦视频高清在线观看| 七月丁香在线播放| 欧美3d第一页| 99久久九九国产精品国产免费| 我的女老师完整版在线观看| 国产精品久久久久久精品电影小说 | 亚洲成人av在线免费| 看非洲黑人一级黄片| 男人的好看免费观看在线视频| 午夜精品国产一区二区电影 | 嫩草影院新地址| 亚洲国产精品久久男人天堂| 三级经典国产精品| 亚洲国产精品专区欧美| 国内精品宾馆在线| 99热精品在线国产| 久久精品久久久久久噜噜老黄 | 欧美不卡视频在线免费观看| 看非洲黑人一级黄片| 国产精华一区二区三区| 亚洲av中文字字幕乱码综合| eeuss影院久久| 亚洲aⅴ乱码一区二区在线播放| 丝袜喷水一区| 国产乱人视频| 狠狠狠狠99中文字幕| 特大巨黑吊av在线直播| av在线蜜桃| 亚洲内射少妇av| 99热这里只有是精品在线观看| 尾随美女入室| 日本wwww免费看| 欧美色视频一区免费| 国产欧美日韩精品一区二区| 成人国产麻豆网| 欧美xxxx黑人xx丫x性爽| 精品国内亚洲2022精品成人| 欧美日韩在线观看h| 丰满少妇做爰视频| 秋霞伦理黄片| 热99在线观看视频| 久久韩国三级中文字幕| av专区在线播放| 国产三级中文精品| 亚洲五月天丁香| 久久精品国产鲁丝片午夜精品| 中文天堂在线官网| 视频中文字幕在线观看| 成人美女网站在线观看视频| av在线观看视频网站免费| videossex国产| 亚洲在线自拍视频| 91aial.com中文字幕在线观看| 一区二区三区四区激情视频| 午夜视频国产福利| 国产 一区精品| 亚洲综合精品二区| 国产色婷婷99| 成人国产麻豆网| 国产在视频线精品| 99久久九九国产精品国产免费| 久久这里只有精品中国| 人妻系列 视频| 熟妇人妻久久中文字幕3abv| 精品久久久久久久久av| 成人二区视频| 久久久亚洲精品成人影院| 国产黄a三级三级三级人| 超碰97精品在线观看| 欧美激情国产日韩精品一区| 国产精品三级大全| 精品人妻熟女av久视频| 国产精品不卡视频一区二区| 婷婷色综合大香蕉| 综合色av麻豆| 长腿黑丝高跟| 欧美激情在线99| av女优亚洲男人天堂| 久久精品国产99精品国产亚洲性色| 真实男女啪啪啪动态图| 国产精品.久久久| 久久久久精品久久久久真实原创| 男女国产视频网站| 观看美女的网站| 能在线免费看毛片的网站| 级片在线观看| 99久久九九国产精品国产免费| 国产三级中文精品| 91久久精品国产一区二区三区| 日韩亚洲欧美综合| 啦啦啦观看免费观看视频高清| 一个人看的www免费观看视频| 亚洲乱码一区二区免费版| 视频中文字幕在线观看| 99热全是精品| 一个人看的www免费观看视频| 99在线人妻在线中文字幕| 嘟嘟电影网在线观看| 成人亚洲欧美一区二区av| 97在线视频观看| 又爽又黄无遮挡网站| 身体一侧抽搐| 一区二区三区四区激情视频| 亚洲欧美成人精品一区二区| 色综合亚洲欧美另类图片| 国产中年淑女户外野战色| 伊人久久精品亚洲午夜| 国产精品日韩av在线免费观看| 干丝袜人妻中文字幕| 亚洲欧美精品自产自拍| 中国国产av一级| 亚洲天堂国产精品一区在线| 自拍偷自拍亚洲精品老妇| 亚洲精品国产av成人精品| a级一级毛片免费在线观看| 久久久久久久久大av| 亚洲欧美中文字幕日韩二区| 国产精品1区2区在线观看.| 久久精品影院6| 亚洲欧洲国产日韩| 熟女人妻精品中文字幕| 三级男女做爰猛烈吃奶摸视频| 最近视频中文字幕2019在线8| 亚洲av电影不卡..在线观看| 男女啪啪激烈高潮av片| 欧美xxxx黑人xx丫x性爽| 熟妇人妻久久中文字幕3abv| 天堂√8在线中文| 日本三级黄在线观看| 国产女主播在线喷水免费视频网站 | 亚洲熟妇中文字幕五十中出| 亚洲一区高清亚洲精品| 日韩人妻高清精品专区| 91在线精品国自产拍蜜月| 国产乱人偷精品视频| 日韩人妻高清精品专区| 精品一区二区三区视频在线| 久久鲁丝午夜福利片| av免费在线看不卡| 免费看日本二区| 久久欧美精品欧美久久欧美| 99久久精品国产国产毛片| 白带黄色成豆腐渣| 男女国产视频网站| 国产伦精品一区二区三区四那| 亚洲精品成人久久久久久| 热99re8久久精品国产| 精品酒店卫生间| 免费看美女性在线毛片视频| 亚洲三级黄色毛片| 久久久成人免费电影| 午夜激情欧美在线| 亚洲激情五月婷婷啪啪| 毛片女人毛片| 波野结衣二区三区在线| 视频中文字幕在线观看| 中文字幕av在线有码专区| 午夜a级毛片| 男人狂女人下面高潮的视频| 精品人妻偷拍中文字幕| 一本一本综合久久| 人人妻人人澡人人爽人人夜夜 | 国产高清国产精品国产三级 | 91午夜精品亚洲一区二区三区| 两个人视频免费观看高清| 国产午夜精品论理片| 国内揄拍国产精品人妻在线| 少妇熟女欧美另类| 寂寞人妻少妇视频99o| 国产精品久久视频播放| 日韩国内少妇激情av| 国产高清视频在线观看网站| 亚洲av不卡在线观看| 国产爱豆传媒在线观看| 亚洲最大成人手机在线| 伦理电影大哥的女人| 18禁在线播放成人免费| 看片在线看免费视频| 亚洲av.av天堂| 欧美丝袜亚洲另类| 成人毛片60女人毛片免费| 久久6这里有精品| 色尼玛亚洲综合影院| 久久国产乱子免费精品| 水蜜桃什么品种好| 久久久亚洲精品成人影院| 日韩高清综合在线| 欧美最新免费一区二区三区| 69人妻影院| 亚洲性久久影院| 黑人高潮一二区| 观看免费一级毛片| 在线观看一区二区三区| 99久国产av精品| 成人毛片a级毛片在线播放| 久久久色成人| av福利片在线观看| 亚洲国产精品成人综合色| 国产欧美日韩精品一区二区| 国产av不卡久久| 精品99又大又爽又粗少妇毛片| ponron亚洲| 亚洲国产色片| 22中文网久久字幕| 欧美一区二区国产精品久久精品| 久久久久久久亚洲中文字幕| 看非洲黑人一级黄片| 免费黄色在线免费观看| 欧美三级亚洲精品| 99在线视频只有这里精品首页| 精品国产一区二区三区久久久樱花 | 国产国拍精品亚洲av在线观看| 高清在线视频一区二区三区 | 麻豆成人午夜福利视频| 尾随美女入室| 亚洲国产精品专区欧美| 成人无遮挡网站| 99热这里只有精品一区| 国产精品.久久久| 亚洲国产精品专区欧美| 麻豆成人午夜福利视频| 小蜜桃在线观看免费完整版高清| 日韩制服骚丝袜av| 国产一级毛片在线| 欧美日韩精品成人综合77777| 非洲黑人性xxxx精品又粗又长| 又粗又爽又猛毛片免费看| 一区二区三区乱码不卡18| 国语自产精品视频在线第100页| 超碰av人人做人人爽久久| 日韩欧美国产在线观看| 亚洲18禁久久av| 成人亚洲精品av一区二区| 少妇丰满av| 国产亚洲一区二区精品| www日本黄色视频网| 国产一区二区亚洲精品在线观看| 日韩欧美在线乱码| 国产亚洲91精品色在线| 国产欧美日韩精品一区二区| 亚洲av福利一区| 国产大屁股一区二区在线视频| 免费电影在线观看免费观看| 国产精品永久免费网站| 五月玫瑰六月丁香| 亚洲图色成人| 色噜噜av男人的天堂激情| 一卡2卡三卡四卡精品乱码亚洲| 免费看a级黄色片| 国模一区二区三区四区视频| 特级一级黄色大片| 91久久精品国产一区二区三区| 啦啦啦观看免费观看视频高清| 国产国拍精品亚洲av在线观看| 国产综合懂色| av免费观看日本| 国产亚洲一区二区精品| 哪个播放器可以免费观看大片| 少妇人妻精品综合一区二区| 国产真实伦视频高清在线观看| 免费黄色在线免费观看| 成人三级黄色视频| 欧美另类亚洲清纯唯美| 精品久久久久久久久亚洲| 亚洲精品乱码久久久v下载方式| 人人妻人人澡欧美一区二区| 干丝袜人妻中文字幕| 午夜福利高清视频| 国产乱来视频区| 亚洲欧美日韩无卡精品| 国产单亲对白刺激| 男人舔女人下体高潮全视频| 亚洲欧美精品综合久久99| 少妇丰满av| 99久国产av精品| 国产精品av视频在线免费观看| .国产精品久久| 国产黄片视频在线免费观看| 三级男女做爰猛烈吃奶摸视频| 国产黄片视频在线免费观看| 女人久久www免费人成看片 | 亚洲av二区三区四区| 18禁裸乳无遮挡免费网站照片| 黄色配什么色好看| videossex国产| 国产黄片美女视频| 精品欧美国产一区二区三| 一级黄色大片毛片| 成人av在线播放网站| 久久久亚洲精品成人影院| 国产人妻一区二区三区在| 久久精品人妻少妇| 水蜜桃什么品种好| 亚洲精品乱码久久久久久按摩| 乱系列少妇在线播放| 日日啪夜夜撸| 少妇熟女欧美另类| 男女下面进入的视频免费午夜| 又爽又黄a免费视频| 亚洲伊人久久精品综合 | 男女啪啪激烈高潮av片| 干丝袜人妻中文字幕| 2022亚洲国产成人精品| 男人狂女人下面高潮的视频| 亚洲欧美精品专区久久| 久久草成人影院| 日韩人妻高清精品专区| 一级黄片播放器| 五月伊人婷婷丁香| 亚洲av电影不卡..在线观看| 亚洲av免费在线观看| 国产色爽女视频免费观看| 国产黄a三级三级三级人| 久久久久免费精品人妻一区二区| 我要搜黄色片| 人体艺术视频欧美日本| 国产精品爽爽va在线观看网站| 熟妇人妻久久中文字幕3abv| 美女内射精品一级片tv| 在线a可以看的网站| 在线免费十八禁| 欧美日本亚洲视频在线播放| 成人鲁丝片一二三区免费| 男人舔女人下体高潮全视频| 人妻夜夜爽99麻豆av| 禁无遮挡网站| 亚洲精品亚洲一区二区| 成年免费大片在线观看| 国产成人91sexporn| 99久久人妻综合| 久久这里只有精品中国| 国产亚洲av片在线观看秒播厂 | 久久久精品94久久精品| 看片在线看免费视频| 搡老妇女老女人老熟妇| 伦精品一区二区三区| 欧美xxxx黑人xx丫x性爽| 一区二区三区高清视频在线| 最近的中文字幕免费完整| 哪个播放器可以免费观看大片| 又爽又黄a免费视频| av免费在线看不卡| 亚洲精华国产精华液的使用体验| 嫩草影院精品99| h日本视频在线播放| 精品一区二区三区视频在线| 男女啪啪激烈高潮av片| 蜜桃久久精品国产亚洲av| 又粗又硬又长又爽又黄的视频| 免费人成在线观看视频色| 亚洲av一区综合| 国内精品一区二区在线观看| 好男人在线观看高清免费视频| 日本与韩国留学比较| 国产v大片淫在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 九九热线精品视视频播放| 看黄色毛片网站| 色播亚洲综合网| 国产高清不卡午夜福利| 国产精品美女特级片免费视频播放器| 亚洲精品自拍成人| 亚洲人成网站在线播| 亚洲伊人久久精品综合 | 国产av在哪里看| 麻豆精品久久久久久蜜桃| 国产伦精品一区二区三区视频9| 成人一区二区视频在线观看| 我要看日韩黄色一级片| 精品不卡国产一区二区三区| 久久国内精品自在自线图片| av在线蜜桃| 国产美女午夜福利| 国产精品福利在线免费观看| 欧美成人免费av一区二区三区| 久久99热这里只有精品18| 桃色一区二区三区在线观看| 久久精品国产自在天天线| 夜夜爽夜夜爽视频| 久久99热这里只频精品6学生 | 精品久久久久久久久av| 岛国毛片在线播放| 免费观看性生交大片5| 成人亚洲欧美一区二区av| 三级毛片av免费| 国产真实伦视频高清在线观看| 国产高清三级在线| 有码 亚洲区| 高清视频免费观看一区二区 | 综合色av麻豆| 小蜜桃在线观看免费完整版高清| 国产亚洲最大av| 91精品一卡2卡3卡4卡| 我要看日韩黄色一级片| 国产午夜精品论理片| 免费观看性生交大片5| 日韩国内少妇激情av| 国产成人a∨麻豆精品| 国产老妇女一区| 日韩亚洲欧美综合| 日韩人妻高清精品专区| 免费在线观看成人毛片| 欧美日韩在线观看h| 国产69精品久久久久777片| 国产一区二区在线观看日韩| 日本三级黄在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 哪个播放器可以免费观看大片| 免费黄网站久久成人精品| 国产成人免费观看mmmm| 精品久久久噜噜| 国产一区二区在线av高清观看| 欧美变态另类bdsm刘玥| 嫩草影院精品99| 亚洲最大成人中文| 久久午夜福利片| 国产一级毛片七仙女欲春2| 水蜜桃什么品种好| 国产三级在线视频| 欧美成人午夜免费资源| 国产男人的电影天堂91| av在线亚洲专区| 亚洲国产最新在线播放| 91久久精品电影网| 久久久国产成人精品二区| 亚洲国产欧洲综合997久久,| 大香蕉久久网| 非洲黑人性xxxx精品又粗又长| 国产成人免费观看mmmm| 精品少妇黑人巨大在线播放 | 免费观看a级毛片全部| 夜夜爽夜夜爽视频| 男女视频在线观看网站免费| 日本免费在线观看一区| 麻豆精品久久久久久蜜桃| 国产欧美另类精品又又久久亚洲欧美| 成人高潮视频无遮挡免费网站| 午夜免费激情av| 少妇人妻一区二区三区视频| 天堂影院成人在线观看| 在线观看美女被高潮喷水网站| 国产精品一区二区在线观看99 | 亚洲婷婷狠狠爱综合网| 伦精品一区二区三区| 国产淫片久久久久久久久| 久久欧美精品欧美久久欧美| 国产高清三级在线| 欧美性猛交黑人性爽| 高清午夜精品一区二区三区| 国产成人aa在线观看| 2021天堂中文幕一二区在线观| 99热这里只有精品一区| 国产黄a三级三级三级人| 三级男女做爰猛烈吃奶摸视频| 欧美3d第一页| 国产一区二区在线av高清观看| 看免费成人av毛片| 男人和女人高潮做爰伦理| 精品国产三级普通话版| 真实男女啪啪啪动态图| 日韩强制内射视频| 一级毛片aaaaaa免费看小| 久久综合国产亚洲精品| 国产白丝娇喘喷水9色精品| 日韩一区二区三区影片| 国产亚洲91精品色在线| 熟妇人妻久久中文字幕3abv| 国产成人aa在线观看| av播播在线观看一区| 久久久精品大字幕| 天天躁日日操中文字幕| 男女啪啪激烈高潮av片| 国产一区有黄有色的免费视频 |