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

    計(jì)算機(jī)輔助催化劑效率因子的模擬計(jì)算

    2023-10-12 05:46:50董松濤董天諾焦博陽
    石油煉制與化工 2023年10期
    關(guān)鍵詞:等距模擬計(jì)算長距離

    董松濤,董天諾,焦博陽

    (1.中石化石油化工科學(xué)研究院有限公司,北京 100083;2.南京大學(xué)軟件學(xué)院;3.大連理工大學(xué)船舶工程學(xué)院 工業(yè)裝備結(jié)構(gòu)分析國家重點(diǎn)實(shí)驗(yàn)室)

    催化劑效率因子對(duì)催化劑的應(yīng)用效果具有重要影響,對(duì)此,之前有不少有價(jià)值的研究和探索[1-4],劉奕等[5]用效率因子來衡量本征速率和表觀速率之間的差別。

    (-RA)宏=η(-RA)本

    (1)

    式中:(-RA)宏為宏觀反應(yīng)速率;(-RA)本為本征反應(yīng)速率;η為效率因子。

    效率因子表示催化劑受內(nèi)擴(kuò)散影響的程度。效率因子越接近1,表明催化劑內(nèi)外濃度差越小,內(nèi)擴(kuò)散影響越小。對(duì)于等溫顆粒而言,效率因子可以表示為西勒模數(shù)的函數(shù),關(guān)系為:

    (2)

    式中,φ為西勒模數(shù)。

    西勒模數(shù)是表征內(nèi)擴(kuò)散過程對(duì)化學(xué)反應(yīng)影響的重要參數(shù),由式(2)可知,可以通過減小西勒模數(shù)來增大反應(yīng)效率因子。西勒模數(shù)是擴(kuò)散距離、反應(yīng)速率常數(shù)和有效擴(kuò)散系數(shù)的函數(shù)。

    (3)

    式中:km為反應(yīng)速率常數(shù);L為擴(kuò)散路徑長度;DA,eff為有效擴(kuò)散系數(shù)。

    對(duì)于一個(gè)給定條件下的反應(yīng),反應(yīng)速率常數(shù)km為定值,通過式(3)可知,提高催化劑孔道的有效擴(kuò)散系數(shù)和縮短擴(kuò)散路徑均可減小西勒模數(shù)??赏ㄟ^優(yōu)化催化劑孔道結(jié)構(gòu)及工藝條件來增加有效擴(kuò)散系數(shù),或者縮短擴(kuò)散路徑。

    為了求解催化劑的效率因子,進(jìn)行了很多有益的探索和研究[6-8]。這些工作都為計(jì)算反應(yīng)效率提供了很好的思路和方法,但這些方法都是基于規(guī)范形狀,如果用于一些非規(guī)范形狀的催化劑,或者實(shí)際生產(chǎn)的存在外形變形,甚至內(nèi)部包含空腔或裂縫的催化劑,就存在一定的局限性。

    進(jìn)一步考慮式(3),根號(hào)內(nèi)有兩項(xiàng),有效擴(kuò)散系數(shù)和反應(yīng)速率常數(shù),如果原材料和制備方法相同,催化劑的孔道結(jié)構(gòu)也相同,則有效擴(kuò)散系數(shù)是一個(gè)定值。對(duì)于一個(gè)特定的反應(yīng),在給定的催化劑和反應(yīng)條件下,反應(yīng)速率常數(shù)也是一個(gè)定值,因此,式(3)可簡化為:

    φ=C×L

    (4)

    式中,C為包含了反應(yīng)速率常數(shù)和有效擴(kuò)散系數(shù)的常數(shù)。

    也就是說,對(duì)于給定的催化劑和反應(yīng)條件,西勒模數(shù)只與擴(kuò)散路徑的長度有關(guān),更確切地說,僅與距表面的距離有關(guān)。因此,對(duì)催化劑外形設(shè)計(jì)而言,并非一定要采用某個(gè)特定的形狀,而是要在保證物理性能的基礎(chǔ)上,使催化劑的活性組分盡可能地貼近表面。活性中心貼近表面的程度可用其距邊界距離與等距點(diǎn)的數(shù)量或計(jì)數(shù)關(guān)系曲線來描述。

    距邊界距離與等距點(diǎn)計(jì)數(shù)之間的關(guān)系可通過2種方法獲得:①對(duì)于一些簡單、規(guī)范的形狀,可以基于幾何原理來計(jì)算;②對(duì)于復(fù)雜的形狀,可以采用計(jì)算機(jī)輔助模擬方法得到,比如地圖處理軟件ArcGIS,圖形處理軟件python+openCV、FIJI等。基本方法為:先確定邊界,然后計(jì)算點(diǎn)到邊界的距離,統(tǒng)計(jì)距邊界等距點(diǎn)的計(jì)數(shù)。為了提升計(jì)算機(jī)輔助模擬方法的可靠性,可用規(guī)范形狀的計(jì)算結(jié)果來校驗(yàn)?zāi)M結(jié)果,對(duì)相關(guān)參數(shù)進(jìn)行優(yōu)化,使得模型實(shí)現(xiàn)盡可能高的實(shí)用性。

    為了直觀表示距離分布情況,用不同的顏色深度來表示點(diǎn)與邊界的距離。為了改善催化劑性能,需要找到對(duì)催化劑效率影響最大的點(diǎn)。用人工方法可以繪制一些簡單形狀的距離與等距點(diǎn)數(shù)量分布圖,如果是更復(fù)雜的形狀,人工繪制難度極大,也少見文獻(xiàn)報(bào)道,這就需要采用計(jì)算機(jī)輔助技術(shù)。計(jì)算機(jī)輔助技術(shù)可以得到幾乎任意圖形的距邊界距離與等距點(diǎn)計(jì)數(shù)的幾何分布圖?;谶@樣的分布圖,能夠直觀地找到效率低點(diǎn)的位置,方便下一步進(jìn)行優(yōu)化。

    1 實(shí)驗(yàn)方法

    對(duì)于一些相對(duì)簡單的外形,可基于幾何原理測定點(diǎn)與邊界的距離。假定二維圖形的周長或三維圖形的外表面積大小與活性點(diǎn)的數(shù)量或計(jì)數(shù)有對(duì)應(yīng)關(guān)系,按照從外向里“剝洋蔥”的方式,每次減去足夠小的“厚度”,則可以持續(xù)得到距邊界距離和該距離下點(diǎn)的計(jì)數(shù)關(guān)系曲線。

    計(jì)算機(jī)輔助模擬方法基本過程為:①取一個(gè)平面圖形(自行繪制或照相);②以像素點(diǎn)為基本單位,找出邊界點(diǎn),獲得所有邊界點(diǎn)坐標(biāo)值;③選擇圖中某像素點(diǎn),并得到其坐標(biāo)值,計(jì)算該像素點(diǎn)與所有邊界點(diǎn)的距離,取這些距離的最小值作為該點(diǎn)距邊界的距離值,遍歷圖形內(nèi)所有像素點(diǎn);④將所有點(diǎn)距邊界的距離值進(jìn)行統(tǒng)計(jì)匯總,得到距邊界距離與等距點(diǎn)計(jì)數(shù)數(shù)組;⑤將計(jì)算得到的距離值進(jìn)行歸一化顯示,用顏色的深淺代表點(diǎn)與邊界距離的長短。具體工具方面,很多軟件都可以實(shí)現(xiàn)。如:python+openCV組合,主要基于openCV中的pointPolygonTest這個(gè)API功能。也可以采用上文提到的圖形或地理軟件,或者采用其他具有類似功能的軟件來實(shí)現(xiàn)。

    整體效率因子根據(jù)文獻(xiàn)的反應(yīng)速率常數(shù)和擴(kuò)散系數(shù)計(jì)算,或者直接設(shè)定式(4)的C值,計(jì)算出距邊界距離下的催化劑效率因子。對(duì)于某種形狀的催化劑,通過計(jì)算機(jī)輔助模擬方法,得到距邊界距離及等距點(diǎn)計(jì)數(shù)曲線,將所有距邊界距離的計(jì)數(shù)和效率因子相乘,取加權(quán)平均值,得到整體效率因子。

    (5)

    式中:η總為總有效率;Nd為距邊界距離為d單位時(shí)的等距點(diǎn)計(jì)數(shù);ηd為距邊界距離為d單位時(shí)的效率因子。

    將所有距邊界的距離按照大小進(jìn)行歸并,為方便與八位灰度圖的繪圖要求相匹配,將數(shù)據(jù)歸并到最大為28=256個(gè)“簇”中,實(shí)際計(jì)算時(shí),僅需對(duì)這些“簇”進(jìn)行計(jì)算,這樣就可以采用常用的辦公軟件實(shí)現(xiàn)整體效率因子的模擬和計(jì)算。由于采用了模擬計(jì)算方法,簡化并歸并到“簇”,因此得到的效率因子是一個(gè)近似結(jié)果。

    距邊界距離幾何分布圖按照灰度的方式表示點(diǎn)距邊界的距離,其中顏色越深,距離邊界越遠(yuǎn),反之亦然。

    2 結(jié)果與討論

    2.1 簡單形狀

    催化劑均具有三維形狀,直接進(jìn)行模擬計(jì)算的難度較大。大部分催化劑是條柱狀,如果假設(shè)催化劑長度為無限長,那么距邊界距離與等距點(diǎn)計(jì)數(shù)之間的關(guān)系僅僅與截面的形狀有關(guān),這樣就實(shí)現(xiàn)了模擬過程的簡化處理。

    有幾種形狀相對(duì)簡單,如圓形、三角形和正方形。

    圓柱形的截面是圓形,也是目前常用的催化劑形狀之一?;趫A的幾何特性,一個(gè)半徑為R的圓形,圓內(nèi)任一點(diǎn)距圓周的距離有兩個(gè)極限值:①圓周上,距離邊界的距離為0,計(jì)數(shù)數(shù)量為周長,即2πR;②圓心,為1個(gè)點(diǎn),距離邊界的距離為R;如果圓中其他點(diǎn)距邊界的距離為R1,則具有相同距離點(diǎn)的數(shù)量為2π(R-R1),從半徑和圓的周長公式可知,距邊界距離和等距點(diǎn)的數(shù)量為線性關(guān)系,如圖1所示。其中,由兩個(gè)坐標(biāo)軸和直線組成的圖形的面積和原圖面積相等,其他形狀也具有同樣特征。

    圖1 圓形等距點(diǎn)幾何分布

    采用計(jì)算機(jī)輔助模擬計(jì)算的分布情況如圖2所示。

    圖2 計(jì)算機(jī)輔助模擬得到的圓形等距點(diǎn)幾何分布

    比較圖1和圖2可見,距邊界距離與等距點(diǎn)計(jì)數(shù)的關(guān)系圖,采用幾何關(guān)系時(shí)預(yù)計(jì)是一條直線(圖1),采用計(jì)算機(jī)輔助模擬得到的曲線(圖2)整體與圖1類似,但是在距離邊界較近距離區(qū)域計(jì)數(shù)點(diǎn)較分散,但仍保持了很好的趨勢。由此表明,計(jì)算機(jī)輔助模擬得到的結(jié)果與采用幾何特性計(jì)算的結(jié)果具有良好的一致性。

    在實(shí)際模擬中發(fā)現(xiàn),由于采用圖形的像素大小問題,圓形在放大后,邊界變成了鋸齒形,這可能是距邊界較短距離區(qū)計(jì)數(shù)點(diǎn)較分散的原因之一。三角形和正方形的情況與圓形類似,分別如圖3和圖4所示。

    圖3 計(jì)算機(jī)輔助模擬得到的三角形等距點(diǎn)幾何分布

    圖4 計(jì)算機(jī)輔助模擬得到的正方形等距點(diǎn)幾何分布

    由圖3和圖4可見,三角形等距點(diǎn)的分布曲線中在距邊界較近距離區(qū)域計(jì)數(shù)點(diǎn)較分散,而正方形就不存在這一情況,這可能與正方形邊界上所有點(diǎn)都是直線,而圖形模擬都是按照直角坐標(biāo)系來對(duì)圖形中的點(diǎn)進(jìn)行標(biāo)注,即便是放大后也不會(huì)出現(xiàn)“鋸齒”有關(guān)。

    將圓形、三角形和正方形的相關(guān)特性進(jìn)行匯總,結(jié)果如表1所示。

    表1 簡單圖形距邊界距離的特性

    表1中斜率用周長與距邊界最長距離的比值來計(jì)算??梢宰⒁獾?這幾種簡單圖形具有一定的共性,如3種形狀的距邊界最長距離都可以表示為:

    (6)

    式中:Lm為距邊界最長距離;S為截面面積;Cb為邊界長度。

    由表1可見,形狀越偏離圓形,距邊界距離與等距點(diǎn)計(jì)數(shù)關(guān)系曲線的斜率也越大。為了提高催化劑的效率因子,希望距邊界最長距離更小,周長更大,相應(yīng)的斜率也更大。又由于距邊界最長距離可以用面積和周長來計(jì)算,也就是希望邊界周長更長,面積更小。

    上述3種簡單形狀,距邊界最長距離都是“一個(gè)點(diǎn)”,究其原因,圓形是等距向內(nèi),三角形和正方形具有角平分線相交于一點(diǎn)的特性,決定了只能匯聚到一個(gè)點(diǎn)上。而復(fù)雜圖形就不一定會(huì)匯聚到一個(gè)點(diǎn)上,另外,距邊界距離與等距點(diǎn)計(jì)數(shù)關(guān)系曲線也不再是簡單的直線。

    2.2 復(fù)雜形狀

    簡單形狀的幾何計(jì)算和計(jì)算機(jī)模擬結(jié)果對(duì)比說明,計(jì)算機(jī)模擬結(jié)果具有較高的可靠性,因此可以對(duì)一些相對(duì)復(fù)雜的圖形進(jìn)行模擬計(jì)算。

    2.2.1圓環(huán)形

    圓環(huán)形是在圓形的基礎(chǔ)上,中心部分加一個(gè)中空的孔道。從圓環(huán)形的幾何形狀可知,兩個(gè)極限點(diǎn)分別為:①距邊界距離為0的點(diǎn),計(jì)數(shù)值為內(nèi)環(huán)周長和外環(huán)周長之和;②距邊界距離最長的點(diǎn),位置在距外環(huán)和內(nèi)環(huán)相等的圓周上,計(jì)數(shù)值為內(nèi)環(huán)周長和外環(huán)周長和的一半。距邊界最近和最遠(yuǎn)點(diǎn)間等距點(diǎn)計(jì)數(shù)曲線為直線,如圖5所示,計(jì)算機(jī)輔助模擬結(jié)果如圖6所示。

    圖5 圓環(huán)形等距點(diǎn)幾何分布

    圖6 計(jì)算機(jī)輔助模擬得到的圓環(huán)形等距點(diǎn)幾何分布

    比較圖5和圖6可見,模擬得到的等距點(diǎn)分布與幾何分析特征一致;距邊界距離與等距點(diǎn)計(jì)數(shù)曲線也高度一致,模擬計(jì)算得到的分布也是與x軸平行的一條線段,在超過距邊界距離最長點(diǎn)后等距點(diǎn)計(jì)數(shù)降為0。

    2.2.2橢圓形

    橢圓形有不少的定義方式,取其中的一種形式:

    (7)

    式中:橢圓的短半軸長度為a;長半軸長度為b。橢圓形的面積S=πab,橢圓的周長等于該橢圓短半軸與長半軸長度之和與該橢圓系數(shù)的積,計(jì)算過程較復(fù)雜,不在此贅述。推測距邊界最長距離為短半軸長度。

    圖7為計(jì)算機(jī)輔助模擬得到的橢圓等距點(diǎn)幾何分布。由圖7可見,距邊界距離和等距點(diǎn)計(jì)數(shù)關(guān)系曲線并非一條直線,而是類似于橢圓右上部分的邊界曲線。這種分布用幾何方法計(jì)算的難度較大。此外,模擬結(jié)果證實(shí),距邊界距離最大值為短軸半徑。

    圖7 計(jì)算機(jī)輔助模擬得到的橢圓形等距點(diǎn)幾何分布

    2.2.3長方形

    長方形相對(duì)比較簡單。對(duì)一個(gè)長邊長度為L、短邊長度為H的長方形,長方形中任一點(diǎn)距離邊界的距離有兩個(gè)極值:①長方形周長上,距離邊界距離為0,計(jì)數(shù)數(shù)量為周長,即2(L+H);②中線,距邊界距離最大值為短邊長度的一半,即H/2,計(jì)數(shù)為2(L-H);長方形內(nèi)其他點(diǎn)的距邊界距離與計(jì)數(shù)點(diǎn)均在這兩點(diǎn)的連線上,如圖8所示。

    圖8 長方形等距點(diǎn)幾何分布

    長方形內(nèi)部距邊界距離與等距點(diǎn)計(jì)數(shù)的計(jì)算機(jī)輔助模擬結(jié)果如圖9所示。圖8和圖9結(jié)果顯示,計(jì)算機(jī)模擬結(jié)果和幾何法得到的結(jié)果具有很好的一致性。在距邊界距離較短區(qū),沒有出現(xiàn)計(jì)數(shù)點(diǎn)分散的情況。

    圖9 計(jì)算機(jī)輔助模擬得到的長方形等距點(diǎn)幾何分布

    2.2.4任意四邊形

    與長方形不同,任意四邊形無法用幾何方法確定距邊界最大距離,但仍然可以進(jìn)行計(jì)算機(jī)輔助模擬計(jì)算,結(jié)果見圖10。

    圖10 計(jì)算機(jī)輔助模擬得到的任意四邊形等距點(diǎn)幾何分布

    需要特別指出的是,任意四邊形距邊界最長距離可能不是一個(gè)點(diǎn),有可能是一條線或者其他形狀。在圖10所示的四邊形中,距邊界距離與等距點(diǎn)計(jì)數(shù)曲線基本呈直線。

    2.2.5心形

    催化劑截面很少出現(xiàn)心形,但是心形比較特殊,用極坐標(biāo)很容易繪制,在直角坐標(biāo)系也可以實(shí)現(xiàn),在此僅作為一個(gè)案例進(jìn)行探討。圖11為計(jì)算機(jī)輔助模擬得到的心形等距點(diǎn)幾何分布。由圖11可見,距邊界距離與等距點(diǎn)計(jì)數(shù)曲線在臨界點(diǎn)之前基本是一條趨勢向下的直線,過了臨界點(diǎn)則變成上凹曲線,最終與x軸相交。

    圖11 計(jì)算機(jī)輔助模擬得到的心形等距點(diǎn)幾何分布

    2.2.6齒輪形

    齒輪形偶見作為催化劑截面形狀,基本由圓環(huán)外套一個(gè)等間隔空缺的圓環(huán)組成,該形狀以函數(shù)方式表達(dá)難度較大,在此作為一個(gè)復(fù)雜形狀的案例進(jìn)行探討。圖12為計(jì)算機(jī)輔助模擬得到的八齒齒輪的等距點(diǎn)幾何分布。由圖12可見,齒輪形的距邊界距離與等距點(diǎn)計(jì)數(shù)曲線類似于心形和環(huán)形的組合,在臨界點(diǎn)之前計(jì)數(shù)值基本相等,過了臨界點(diǎn)后變?yōu)樯习记€,最終與x軸相交。

    圖12 計(jì)算機(jī)輔助模擬得到的齒輪形等距點(diǎn)幾何分布

    2.3 三維形狀

    三維圖形計(jì)算的難點(diǎn)是如何轉(zhuǎn)化為二維的距邊界距離與等距點(diǎn)計(jì)數(shù)曲線。

    所有催化劑都是以某種三維形狀出現(xiàn)的,純粹的三維形狀的模擬計(jì)算難度大,需要將三維形狀進(jìn)行降維處理,以方便計(jì)算和模擬。一種簡化方法,如前文采用的,將柱狀條的條長假定為無限長,不考慮長度的影響,僅考慮截面形狀;另外一種簡化方法,對(duì)一些“規(guī)范”的形狀,則可基于其幾何特性,采用一些簡化處理。

    (1)簡單形狀(球體、立方體、正四面體)。對(duì)于一些簡單形狀,比如球體,一個(gè)半徑為R的球體,圓球中任一點(diǎn)距球表面的距離有兩個(gè)極值:①球面上,距球面距離為0,計(jì)數(shù)數(shù)量為球的表面積,即4πR2;②球心,數(shù)量為1,距離圓球外表面的距離為R;若圓球中其他點(diǎn)距表面的距離為R1,則具有相同距離點(diǎn)的數(shù)量相當(dāng)于過該點(diǎn)的球面積,數(shù)值為 4π(R-R1)2計(jì)數(shù)。由圓球的表面積公式可知,距球表面距離與等距點(diǎn)計(jì)數(shù)的平方根呈線性關(guān)系。立方體和正四面體具有類似的情況。

    (2)柱狀體。在此假定柱狀體為圓柱體,如果橫截面為其他形狀,特性基本類似。柱狀體的模擬結(jié)果與柱狀體的長度有關(guān),要分為兩種情況,一種是柱體半徑小于柱長,一種是柱體半徑大于柱長,其中前一種情況比較多見。沿軸向過圓柱軸心的橫截面是一個(gè)長方形,從幾何特性分析,距邊界距離分布圖和長方形類似(見圖9)。

    一個(gè)長度為h、半徑為r的圓柱體,極限點(diǎn)分別為:①圓柱表面上,距離邊界距離為0,等距點(diǎn)計(jì)數(shù)為外表面積S=2πr2+2πrh;②中線,距離圓柱表面的最大距離點(diǎn)為圓柱形的半徑r。圓柱體內(nèi)的一點(diǎn),如果距圓柱體表面的距離為d,則等距點(diǎn)計(jì)數(shù)為S內(nèi)=2π(r-d)(h-d)+2π(r-d)2。由圓柱體的表面積計(jì)算公式可知,距圓柱表面的距離與等距點(diǎn)計(jì)數(shù)的平方根呈線性關(guān)系。

    2.4 催化劑形狀對(duì)效率因子的影響

    對(duì)于常見的催化劑外形,趙毅等[9]采用有限差分方法及松弛迭代法對(duì)三維體系的擴(kuò)散反應(yīng)方程進(jìn)行了求解,得到了幾種具有典型特征形狀催化劑的效率因子。

    對(duì)于復(fù)雜截面的條柱形,其三維結(jié)構(gòu)的模擬較復(fù)雜。假設(shè)催化劑條為無限長,則只需模擬計(jì)算截面形狀不同的催化劑的效率因子即可,結(jié)果如表2所示。

    表2 模擬計(jì)算值與文獻(xiàn)值比較

    表2中圓柱形、三葉形、環(huán)形、五葉形的形狀定義與文獻(xiàn)[9]一致,包括帶孔三葉形和心形的形狀,如圖13所示。

    圖13 催化劑截面形狀

    由表2可見,采用計(jì)算機(jī)輔助模擬計(jì)算得到的效率因子與文獻(xiàn)值相比,最大絕對(duì)誤差為0.036,最大相對(duì)誤差為8.3%,表明兩者具有較好的一致性。表2也提供了文獻(xiàn)未提及形狀及尺寸下的效率因子,同時(shí)提供了距邊界最大距離的數(shù)值。

    比較帶孔和不帶孔三葉形催化劑的效率因子,在同樣尺寸條件下,中心帶孔催化劑的效率因子比不帶孔催化劑提高11.8%(外徑1.8 mm),距界面最大距離縮短幅度更大。由前文分析可知:在距邊界最長距離的點(diǎn)處開孔,或者將該點(diǎn)轉(zhuǎn)化為向外的邊界,對(duì)效率因子的提升具有較大的幫助;反之,對(duì)原本效率因子較高的靠近邊界的點(diǎn)進(jìn)行改進(jìn),效率因子的提升幅度就比較有限。

    表2中提到的三葉形、帶孔三葉形和五葉形的距邊界距離與等距點(diǎn)分布如圖14所示。

    圖14 多葉形形狀及距表面距離分布

    在距邊界最長距離相等的條件下,不同形狀催化劑的效率因子如表3所示。由表3可見,在距邊界最長距離等于0.3 mm的情況下,除環(huán)形外,其他幾種形狀催化劑的效率因子很接近,最大值為0.755,最小值為0.719,這些外形與環(huán)形的最大差別是距邊界距離和等距點(diǎn)計(jì)數(shù)曲線是一條隨距離增加單調(diào)減少的曲線,而環(huán)形的是一條與x軸平行的直線。

    表3 距邊界最長距離相同時(shí)不同形狀催化劑的效率因子

    由表3還可見,在距邊界最長距離相等的條件下,催化劑最小外徑尺寸由大到小的順序?yàn)?環(huán)形>帶孔三葉形>三葉形>五葉形>心形>圓柱形,這一順序也代表催化劑形狀的效率。通過改變催化劑的形狀,可以用更大外徑實(shí)現(xiàn)小外徑圓柱形同樣的性能,對(duì)催化劑的工業(yè)生產(chǎn)以及使用均有益。

    距邊界最長距離與不同形狀催化劑效率因子的關(guān)系如圖15所示。由圖15可見,距邊界最長距離與效率因子之間的相關(guān)性較強(qiáng),包括環(huán)形在內(nèi),均存在距邊界最長距離越小效率因子越高這一基本規(guī)律。

    圖15 距邊界最長距離與不同形狀催化劑效率因子的關(guān)系

    如果將環(huán)形的效率因子進(jìn)行適當(dāng)校正,校正值取10%~20%,則可以基于距邊界最長距離與效率因子之間的關(guān)系,推算環(huán)形催化劑的效率因子。在實(shí)際應(yīng)用中,如果是相對(duì)較簡單的形狀,距邊界最長距離可以采用式(6)來估算,但形狀與規(guī)范形狀差別越大,相應(yīng)的偏差也會(huì)越大。

    由于采用了計(jì)算機(jī)輔助模擬圖像處理法,催化劑可以具有任意形狀,工業(yè)生產(chǎn)中即便是出現(xiàn)了不規(guī)則形狀,包括出現(xiàn)扭曲、變形、裂縫和起泡等情況,也都可以快速、便捷地計(jì)算出效率因子,為非規(guī)范形狀催化劑效率因子的計(jì)算開辟了一條可行的通道。

    3 結(jié) 論

    催化劑的效率因子與其外形有關(guān),特別是與距邊界等距點(diǎn)的分布有關(guān),研究了幾種簡單形狀催化劑的幾何特性,得到了相應(yīng)的距邊界距離與等距點(diǎn)計(jì)數(shù)的分布曲線。

    提出了一種計(jì)算機(jī)輔助模擬方法,結(jié)合圖像識(shí)別模擬計(jì)算法,可獲得一個(gè)圖中所有點(diǎn)距邊界距離與等距點(diǎn)計(jì)數(shù)的關(guān)系,經(jīng)過驗(yàn)證,模擬計(jì)算獲得的結(jié)果和幾何分析得到的結(jié)果一致,表明計(jì)算機(jī)輔助模擬方法具有較高的可靠性。

    采用計(jì)算機(jī)輔助方法,模擬并給出了環(huán)形、橢圓形、長方形、任意四邊形、心形和齒輪形的等距點(diǎn)幾何分布,基于計(jì)算機(jī)輔助模擬方法,可以得到更復(fù)雜形狀的等距點(diǎn)分布曲線。

    采用計(jì)算機(jī)輔助模擬計(jì)算得到的效率因子與文獻(xiàn)結(jié)果基本一致,相對(duì)誤差不大于10%,具有較好的一致性。結(jié)果證實(shí),除環(huán)形外,距界面最大距離與催化劑的效率因子之間具有較高的相關(guān)性。

    猜你喜歡
    等距模擬計(jì)算長距離
    R1234ze PVTx熱物性模擬計(jì)算
    能源工程(2022年1期)2022-03-29 01:06:26
    擬凸Hartogs域到復(fù)空間形式的全純等距嵌入映射的存在性
    水錘防護(hù)措施在某長距離供水系統(tǒng)中的應(yīng)用
    支持長距離4K HDR傳輸 AudioQuest Pearl、 Forest、 Cinnamon HDMI線
    保持算子束部分等距的映射
    擠出發(fā)泡片材褶皺分析及模擬計(jì)算
    中國塑料(2016年9期)2016-06-13 03:18:54
    等距延拓以及相關(guān)問題
    長距離多起伏輸水管道爆管分析與防護(hù)探討
    我國最長距離特高壓輸電工程開工
    實(shí)際發(fā)射工況下底排藥柱結(jié)構(gòu)完整性的模擬計(jì)算
    欧美日韩av久久| 亚洲av日韩精品久久久久久密| 999久久久国产精品视频| 超碰成人久久| 正在播放国产对白刺激| 黄色片一级片一级黄色片| a级毛片在线看网站| 成年人免费黄色播放视频| 国产真人三级小视频在线观看| 不卡一级毛片| 99国产精品一区二区三区| 国产亚洲欧美精品永久| 欧美色视频一区免费| 中文字幕另类日韩欧美亚洲嫩草| 涩涩av久久男人的天堂| 咕卡用的链子| 高清av免费在线| 欧美精品一区二区免费开放| 亚洲av日韩精品久久久久久密| 久久久久精品人妻al黑| 天天躁狠狠躁夜夜躁狠狠躁| 91麻豆精品激情在线观看国产 | 男人舔女人的私密视频| 丝袜在线中文字幕| 老司机靠b影院| 国产成人系列免费观看| 亚洲男人天堂网一区| 男女床上黄色一级片免费看| 狠狠狠狠99中文字幕| 亚洲熟女精品中文字幕| 国产成人一区二区三区免费视频网站| 三级毛片av免费| 久久人妻熟女aⅴ| 精品乱码久久久久久99久播| 亚洲成人免费av在线播放| 欧美日韩黄片免| 水蜜桃什么品种好| netflix在线观看网站| 9色porny在线观看| 777米奇影视久久| 成人黄色视频免费在线看| а√天堂www在线а√下载 | 午夜成年电影在线免费观看| www.精华液| 中文字幕另类日韩欧美亚洲嫩草| 久久久久精品人妻al黑| 日韩欧美一区二区三区在线观看 | 老汉色av国产亚洲站长工具| 日韩精品免费视频一区二区三区| 欧美+亚洲+日韩+国产| 欧美日韩精品网址| 精品国产国语对白av| 日日爽夜夜爽网站| 搡老乐熟女国产| 亚洲综合色网址| 少妇的丰满在线观看| 99国产综合亚洲精品| 美女午夜性视频免费| 国产欧美日韩一区二区三区在线| 久久香蕉精品热| 9热在线视频观看99| 女人爽到高潮嗷嗷叫在线视频| 久久青草综合色| 黑人猛操日本美女一级片| 大码成人一级视频| 国产成人av教育| 18禁裸乳无遮挡免费网站照片 | 91字幕亚洲| 亚洲国产欧美日韩在线播放| 黄网站色视频无遮挡免费观看| 亚洲国产精品合色在线| 久久精品国产综合久久久| 少妇猛男粗大的猛烈进出视频| 一级毛片精品| 看片在线看免费视频| 99国产精品99久久久久| 亚洲免费av在线视频| 国产三级黄色录像| 国产精品二区激情视频| 国产精品免费一区二区三区在线 | 久久中文字幕人妻熟女| 国产精品av久久久久免费| 美女高潮到喷水免费观看| 精品人妻熟女毛片av久久网站| 人人妻人人爽人人添夜夜欢视频| 精品乱码久久久久久99久播| 1024香蕉在线观看| 日韩精品免费视频一区二区三区| 亚洲第一av免费看| tube8黄色片| 国产亚洲精品第一综合不卡| 久久人人97超碰香蕉20202| 成人av一区二区三区在线看| 黑人巨大精品欧美一区二区蜜桃| www.自偷自拍.com| 999久久久精品免费观看国产| 欧美日韩av久久| 美女国产高潮福利片在线看| 亚洲午夜精品一区,二区,三区| 免费在线观看视频国产中文字幕亚洲| 免费一级毛片在线播放高清视频 | 亚洲少妇的诱惑av| 少妇粗大呻吟视频| 亚洲欧美日韩另类电影网站| 变态另类成人亚洲欧美熟女 | 国产免费现黄频在线看| 国产精品免费一区二区三区在线 | 亚洲欧美一区二区三区黑人| 久久中文字幕一级| 久久香蕉国产精品| 欧美日本中文国产一区发布| 亚洲熟妇中文字幕五十中出 | 在线视频色国产色| 亚洲国产毛片av蜜桃av| 日韩制服丝袜自拍偷拍| 超碰97精品在线观看| 身体一侧抽搐| 亚洲熟女毛片儿| 一边摸一边做爽爽视频免费| 中文亚洲av片在线观看爽 | 国内毛片毛片毛片毛片毛片| 亚洲av日韩精品久久久久久密| 亚洲熟妇熟女久久| 成年人免费黄色播放视频| 91精品三级在线观看| 一本一本久久a久久精品综合妖精| 亚洲精品一二三| 99热只有精品国产| 99精品欧美一区二区三区四区| 亚洲精品国产精品久久久不卡| 婷婷精品国产亚洲av在线 | 十八禁人妻一区二区| 久久久久国内视频| 超色免费av| 18禁裸乳无遮挡免费网站照片 | 亚洲精品一二三| 亚洲男人天堂网一区| 午夜福利,免费看| 日韩欧美三级三区| 精品人妻熟女毛片av久久网站| 色94色欧美一区二区| 欧美黄色片欧美黄色片| 精品电影一区二区在线| 久久人人97超碰香蕉20202| 乱人伦中国视频| 久久久国产欧美日韩av| 亚洲成人免费av在线播放| 啦啦啦免费观看视频1| 午夜福利欧美成人| 日韩 欧美 亚洲 中文字幕| 在线观看一区二区三区激情| 天天影视国产精品| avwww免费| 国产蜜桃级精品一区二区三区 | 成人国语在线视频| 精品欧美一区二区三区在线| 中文字幕精品免费在线观看视频| 久久久久久久国产电影| 欧美精品高潮呻吟av久久| 国产在线观看jvid| 午夜福利一区二区在线看| 国产av一区二区精品久久| 国产精品国产av在线观看| 99国产精品免费福利视频| 国产一区二区激情短视频| 国产av精品麻豆| 18禁美女被吸乳视频| 18禁美女被吸乳视频| cao死你这个sao货| 成年人黄色毛片网站| 精品熟女少妇八av免费久了| 国产成人av教育| 国产成人免费观看mmmm| 国产成人免费观看mmmm| 99久久综合精品五月天人人| 欧美日韩亚洲国产一区二区在线观看 | 国产精品香港三级国产av潘金莲| 男女午夜视频在线观看| 国产精品香港三级国产av潘金莲| 欧美黑人精品巨大| 深夜精品福利| 亚洲中文av在线| 亚洲欧美一区二区三区黑人| 1024视频免费在线观看| 国产精华一区二区三区| 69精品国产乱码久久久| 国产精品秋霞免费鲁丝片| 日本vs欧美在线观看视频| 在线播放国产精品三级| 亚洲色图综合在线观看| 91在线观看av| 亚洲自偷自拍图片 自拍| 国产乱人伦免费视频| 黄色 视频免费看| 黑人巨大精品欧美一区二区蜜桃| 无遮挡黄片免费观看| 国产又色又爽无遮挡免费看| 国产精品久久久人人做人人爽| 亚洲男人天堂网一区| 成人亚洲精品一区在线观看| 香蕉丝袜av| 久久人妻熟女aⅴ| 视频在线观看一区二区三区| 欧美成狂野欧美在线观看| 久久久精品国产亚洲av高清涩受| 亚洲成a人片在线一区二区| 久久久国产一区二区| 黄色a级毛片大全视频| 少妇猛男粗大的猛烈进出视频| 高清欧美精品videossex| 国产精品亚洲一级av第二区| 精品少妇久久久久久888优播| 亚洲国产精品一区二区三区在线| 十八禁人妻一区二区| 欧美成狂野欧美在线观看| 国产成人免费观看mmmm| 国产精品久久电影中文字幕 | 男男h啪啪无遮挡| 国产在视频线精品| 国产精品亚洲一级av第二区| 91字幕亚洲| 亚洲国产毛片av蜜桃av| 美女高潮到喷水免费观看| 五月开心婷婷网| 老鸭窝网址在线观看| 后天国语完整版免费观看| 一进一出抽搐动态| 91麻豆av在线| 黑人欧美特级aaaaaa片| 搡老熟女国产l中国老女人| 欧美精品亚洲一区二区| 午夜视频精品福利| 亚洲美女黄片视频| 日韩欧美一区二区三区在线观看 | 久久久久久久久久久久大奶| √禁漫天堂资源中文www| 日本黄色视频三级网站网址 | 国产亚洲欧美精品永久| 欧美精品啪啪一区二区三区| 亚洲专区中文字幕在线| 在线天堂中文资源库| 国产欧美日韩综合在线一区二区| 久久人人爽av亚洲精品天堂| 国产亚洲一区二区精品| 精品国产一区二区三区四区第35| 另类亚洲欧美激情| 久久人妻福利社区极品人妻图片| 久久精品亚洲av国产电影网| 亚洲人成伊人成综合网2020| 两个人看的免费小视频| 热re99久久国产66热| 成人av一区二区三区在线看| 一进一出好大好爽视频| 亚洲av第一区精品v没综合| aaaaa片日本免费| 人妻丰满熟妇av一区二区三区 | 色在线成人网| 亚洲av日韩精品久久久久久密| 久9热在线精品视频| 婷婷丁香在线五月| 黑人操中国人逼视频| 欧美激情久久久久久爽电影 | 国产国语露脸激情在线看| 中文字幕最新亚洲高清| 国产乱人伦免费视频| 欧美在线黄色| 51午夜福利影视在线观看| 久久这里只有精品19| 午夜91福利影院| 欧美日韩亚洲高清精品| 国产欧美日韩一区二区三| 国产成人精品在线电影| a在线观看视频网站| 国产男女内射视频| 精品久久久精品久久久| 老汉色∧v一级毛片| 日韩欧美一区视频在线观看| 黄色视频,在线免费观看| 久久香蕉激情| 亚洲欧洲精品一区二区精品久久久| 岛国在线观看网站| 国产真人三级小视频在线观看| 黑丝袜美女国产一区| 成年女人毛片免费观看观看9 | 亚洲专区国产一区二区| 国产精品影院久久| 日韩欧美一区视频在线观看| 国产av精品麻豆| 国产欧美日韩综合在线一区二区| 黄色怎么调成土黄色| 黄色成人免费大全| 人人妻,人人澡人人爽秒播| 天堂√8在线中文| 欧美精品一区二区免费开放| 欧美精品av麻豆av| 亚洲国产精品合色在线| videos熟女内射| 黑丝袜美女国产一区| 在线观看日韩欧美| 午夜免费鲁丝| 亚洲精品在线美女| 深夜精品福利| 99riav亚洲国产免费| 国产片内射在线| 黄色成人免费大全| 国产91精品成人一区二区三区| 成人黄色视频免费在线看| 人人澡人人妻人| 日韩中文字幕欧美一区二区| 日韩欧美一区视频在线观看| 满18在线观看网站| 99久久综合精品五月天人人| 亚洲va日本ⅴa欧美va伊人久久| 大型黄色视频在线免费观看| 黄片小视频在线播放| 香蕉丝袜av| av福利片在线| 啦啦啦视频在线资源免费观看| 国产视频一区二区在线看| 一区福利在线观看| 欧美精品人与动牲交sv欧美| 午夜成年电影在线免费观看| 亚洲精品在线美女| 人人妻人人爽人人添夜夜欢视频| 亚洲av欧美aⅴ国产| 免费观看人在逋| 久久这里只有精品19| 中文字幕av电影在线播放| 国产成人欧美| 两个人免费观看高清视频| 精品国产一区二区三区久久久樱花| bbb黄色大片| 午夜免费鲁丝| 国产一区二区三区视频了| 久久久久精品人妻al黑| 后天国语完整版免费观看| 日日爽夜夜爽网站| 丝瓜视频免费看黄片| 国产精品.久久久| 一级作爱视频免费观看| 免费一级毛片在线播放高清视频 | 高清毛片免费观看视频网站 | 久久久水蜜桃国产精品网| 久久性视频一级片| 正在播放国产对白刺激| 中文亚洲av片在线观看爽 | 亚洲综合色网址| 757午夜福利合集在线观看| 啦啦啦免费观看视频1| а√天堂www在线а√下载 | 亚洲美女黄片视频| 男女之事视频高清在线观看| 女人精品久久久久毛片| 久久精品亚洲熟妇少妇任你| 国产av一区二区精品久久| 日韩人妻精品一区2区三区| 欧美大码av| 一区二区三区激情视频| 宅男免费午夜| 成年人黄色毛片网站| 久久人妻av系列| 老汉色∧v一级毛片| 国产精品免费一区二区三区在线 | 免费久久久久久久精品成人欧美视频| 欧美乱色亚洲激情| av在线播放免费不卡| 母亲3免费完整高清在线观看| 我的亚洲天堂| 老司机靠b影院| 亚洲精品美女久久久久99蜜臀| 久久久国产成人免费| 欧美最黄视频在线播放免费 | 50天的宝宝边吃奶边哭怎么回事| 日本黄色视频三级网站网址 | 成人18禁在线播放| 亚洲专区字幕在线| 老熟女久久久| 国产男女超爽视频在线观看| 色婷婷av一区二区三区视频| 国产在视频线精品| 国产xxxxx性猛交| 交换朋友夫妻互换小说| 欧美大码av| 亚洲精品国产区一区二| 天堂√8在线中文| 很黄的视频免费| 欧美国产精品va在线观看不卡| 18禁国产床啪视频网站| 一区二区三区精品91| 欧美久久黑人一区二区| 亚洲全国av大片| 欧美日韩乱码在线| 中文字幕最新亚洲高清| 午夜精品国产一区二区电影| 日韩熟女老妇一区二区性免费视频| 亚洲成a人片在线一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 人人妻人人爽人人添夜夜欢视频| 国产精品 国内视频| www.熟女人妻精品国产| 日韩中文字幕欧美一区二区| 女性生殖器流出的白浆| 三级毛片av免费| 国产成+人综合+亚洲专区| 成年动漫av网址| 久久天躁狠狠躁夜夜2o2o| 久久影院123| 日韩制服丝袜自拍偷拍| 99精国产麻豆久久婷婷| 国产成人欧美在线观看 | 天堂动漫精品| x7x7x7水蜜桃| 国产精品99久久99久久久不卡| 捣出白浆h1v1| 在线av久久热| 一边摸一边抽搐一进一出视频| 精品国产美女av久久久久小说| 中文欧美无线码| 韩国av一区二区三区四区| 国产av一区二区精品久久| 国产av又大| 国产精品乱码一区二三区的特点 | 麻豆国产av国片精品| 亚洲午夜理论影院| 午夜福利视频在线观看免费| 久久精品91无色码中文字幕| 日日爽夜夜爽网站| 国产一区在线观看成人免费| 亚洲黑人精品在线| 正在播放国产对白刺激| ponron亚洲| 交换朋友夫妻互换小说| 黄片播放在线免费| 大型黄色视频在线免费观看| 久久中文字幕人妻熟女| 国产欧美日韩一区二区精品| 免费久久久久久久精品成人欧美视频| 久久九九热精品免费| 极品人妻少妇av视频| 一本大道久久a久久精品| 人妻 亚洲 视频| 国产极品粉嫩免费观看在线| 在线十欧美十亚洲十日本专区| 天堂√8在线中文| 欧美成人免费av一区二区三区 | 女人高潮潮喷娇喘18禁视频| 狠狠狠狠99中文字幕| 多毛熟女@视频| 男女午夜视频在线观看| 亚洲人成电影观看| 黑人操中国人逼视频| av在线播放免费不卡| 如日韩欧美国产精品一区二区三区| 天天躁日日躁夜夜躁夜夜| 99re6热这里在线精品视频| 国产成人av教育| 国产精品98久久久久久宅男小说| 天堂√8在线中文| 欧美一级毛片孕妇| 999久久久国产精品视频| 亚洲国产精品一区二区三区在线| 精品久久久久久久久久免费视频 | 精品久久久久久久久久免费视频 | 脱女人内裤的视频| 91成人精品电影| av超薄肉色丝袜交足视频| 男男h啪啪无遮挡| 宅男免费午夜| 亚洲人成伊人成综合网2020| 亚洲中文字幕日韩| 成在线人永久免费视频| 国产人伦9x9x在线观看| 天天躁夜夜躁狠狠躁躁| 久久久久久久精品吃奶| 久久人人爽av亚洲精品天堂| 十八禁高潮呻吟视频| 国产精品久久久久久人妻精品电影| 午夜免费鲁丝| 亚洲欧美色中文字幕在线| 久久九九热精品免费| 欧美国产精品va在线观看不卡| 亚洲精品美女久久av网站| 成人三级做爰电影| 久久久精品区二区三区| 亚洲av日韩在线播放| 日日爽夜夜爽网站| 精品视频人人做人人爽| 久9热在线精品视频| 亚洲va日本ⅴa欧美va伊人久久| 精品人妻熟女毛片av久久网站| 大片电影免费在线观看免费| 悠悠久久av| 交换朋友夫妻互换小说| 精品卡一卡二卡四卡免费| 99re6热这里在线精品视频| 亚洲成a人片在线一区二区| 每晚都被弄得嗷嗷叫到高潮| 午夜福利在线观看吧| 在线观看免费高清a一片| av视频免费观看在线观看| 性少妇av在线| 在线观看66精品国产| 亚洲中文av在线| 多毛熟女@视频| 亚洲欧美色中文字幕在线| 国产欧美日韩综合在线一区二区| 日韩中文字幕欧美一区二区| 电影成人av| 757午夜福利合集在线观看| 午夜福利免费观看在线| 免费av中文字幕在线| 国产精品影院久久| 国产精华一区二区三区| 老司机深夜福利视频在线观看| 天堂动漫精品| 成人18禁在线播放| 大型黄色视频在线免费观看| 人成视频在线观看免费观看| 亚洲av片天天在线观看| 人妻一区二区av| 国产97色在线日韩免费| 午夜精品国产一区二区电影| 中国美女看黄片| 欧美人与性动交α欧美软件| 岛国毛片在线播放| 制服人妻中文乱码| 久久狼人影院| 亚洲一区高清亚洲精品| av线在线观看网站| 嫁个100分男人电影在线观看| 亚洲午夜精品一区,二区,三区| 成在线人永久免费视频| 亚洲色图综合在线观看| 国产精品免费一区二区三区在线 | 丰满人妻熟妇乱又伦精品不卡| 国产在线一区二区三区精| 久99久视频精品免费| 丰满饥渴人妻一区二区三| av网站免费在线观看视频| 午夜91福利影院| 日韩视频一区二区在线观看| 免费观看精品视频网站| 亚洲第一青青草原| 欧美成人免费av一区二区三区 | 精品人妻熟女毛片av久久网站| 夫妻午夜视频| av欧美777| 欧美一级毛片孕妇| 免费黄频网站在线观看国产| 国产亚洲欧美98| 欧美 日韩 精品 国产| 大香蕉久久成人网| 亚洲成人免费电影在线观看| 久久天躁狠狠躁夜夜2o2o| 久久久久国产精品人妻aⅴ院 | 啦啦啦视频在线资源免费观看| 亚洲精品美女久久av网站| 亚洲情色 制服丝袜| 久久中文看片网| 搡老乐熟女国产| a级片在线免费高清观看视频| 久久久久精品国产欧美久久久| 亚洲精品久久成人aⅴ小说| 精品欧美一区二区三区在线| 91成人精品电影| 亚洲欧美激情综合另类| 午夜精品国产一区二区电影| 后天国语完整版免费观看| 亚洲中文字幕日韩| 国产精品成人在线| 麻豆乱淫一区二区| 嫁个100分男人电影在线观看| 在线观看免费视频日本深夜| 一级片免费观看大全| 国产精品国产高清国产av | 少妇猛男粗大的猛烈进出视频| 在线天堂中文资源库| 窝窝影院91人妻| 19禁男女啪啪无遮挡网站| 我的亚洲天堂| 女人久久www免费人成看片| 一级毛片女人18水好多| 99精品欧美一区二区三区四区| 黑人欧美特级aaaaaa片| 91麻豆精品激情在线观看国产 | 国产精品 欧美亚洲| 精品国产美女av久久久久小说| 午夜成年电影在线免费观看| 两个人免费观看高清视频| av一本久久久久| 亚洲精品一二三| 人人澡人人妻人| 久久久久精品人妻al黑| 窝窝影院91人妻| a级毛片黄视频| 日本a在线网址| 九色亚洲精品在线播放| 色精品久久人妻99蜜桃| 亚洲精品成人av观看孕妇| 欧美日韩成人在线一区二区| 国产精品98久久久久久宅男小说| 亚洲人成电影观看| 超碰97精品在线观看| 在线观看日韩欧美| 乱人伦中国视频| 亚洲成国产人片在线观看| 成年动漫av网址| 久久精品国产a三级三级三级| 亚洲欧美一区二区三区久久| 丰满的人妻完整版| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成人免费av在线播放| 美女扒开内裤让男人捅视频| 亚洲第一青青草原| 99久久99久久久精品蜜桃|