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

    基于微米級顆粒臨界沉積/剝離標準的研究進展

    2024-01-16 11:29:52邵宏勛謝俊桂玉雙李潤東
    化工進展 2023年12期
    關鍵詞:恢復系數(shù)法向入射角

    邵宏勛,謝俊,桂玉雙,李潤東

    (沈陽航空航天大學能源與環(huán)境學院,遼寧 沈陽110136)

    顆粒沉積現(xiàn)象普遍存在于自然界與工業(yè)生產(chǎn)中,例如:當微米級顆粒流經(jīng)渦輪葉片的表面,葉片上氣膜孔的微小結構極易發(fā)生沉積現(xiàn)象,造成冷卻孔的堵塞,進而影響飛機空氣動力學特性,甚至引發(fā)動態(tài)失速的可能[1];此外,在鍋爐的運行中,爐膛內氣流攜帶的飛灰顆粒會沉積在換熱管面上誘發(fā)結渣,進一步增大鍋爐換熱熱阻,致使換熱效率降低,嚴重時還會引起鍋爐爆管及停機停產(chǎn)[2-3]。因此,為合理控制顆粒沉積帶來的危害,提出切實可行的沉積/剝離判據(jù)具有良好的現(xiàn)實意義。

    近年來,大批學者致力于開發(fā)應用更為深入的顆粒沉積及剝離的理論與依據(jù),Chen等[4]基于三維JKR的離散元法(DEM)來研究顆粒不同入射角度下的碰撞,得到一個廣義表達式作為顆粒的黏附/反彈判據(jù)。Zhang 等[5]借助高速相機記錄了微米級石墨顆粒碰撞的全過程,獲得了不同尺寸顆粒下的臨界捕集速度和恢復系數(shù)。Dong等[6]結合黏性、彈塑性以及滑移特性提出了一種新型微米級顆粒DEM 模型,能夠很好地對恢復系數(shù)和臨界捕集速度進行預測。Mu等[7]耦合非定常氣固湍流建立了一個數(shù)值沉積模型,以預測具有灰顆粒的黏附/反彈特性。Kleinhans 等[3]、Fang 等[8]也簡要匯總了一些顆粒沉積準則。盡管許多學者提出了不同情況下的顆粒沉積及剝離判據(jù),并得到了不同程度的應用,然而面對復雜多樣的沉積判據(jù),缺乏合理的總結歸納與分析,對相關判據(jù)的聯(lián)系與應用評判還有待進一步補充。

    本文以鍋爐為背景,闡述了相應的沉積機理與特性,然后基于歷年顆粒沉積判據(jù)的研究成果,從多種沉積因素(撞擊參數(shù)、顆粒粒徑、溫度)出發(fā),對微顆粒黏附判據(jù)之間的不足、聯(lián)系、應用與改進進行分析綜述,最后根據(jù)顆粒起動方式對剝離模型進行整合與展望,為進一步深入判據(jù)的數(shù)值推廣與工程應用提供方案參考。

    1 細顆粒沉積機理與特性

    積灰污染的本質為顆粒與壁面或沉積層碰撞后黏附在表面的一個過程,具體可以分為運輸過程、撞擊過程以及黏附過程。值得一提的是,并非所有顆粒都會發(fā)生黏附,這與背后的黏附機理息息相關,涉及燃燒、流體及材料化學等多種應用理論。此外,細顆粒的沉積機制也不盡相同,常見的有擴散、攔截、慣性碰撞及沉降等[9-11]。而對于鍋爐系統(tǒng)中的飛灰顆粒,按照顆粒輸運方式與積灰比重的不同,公認的沉積機制包含以下方面:慣性碰撞、湍流擴散、熱遷移擴散、冷凝機理以及多變的化學反應[12]。具體情況如圖1 所示[13],盡管沉積類型繁多,但從積灰份額來看,慣性碰撞無疑是絕對的主導,而如冷凝、擴散等機制無關緊要[14]。慣性顆粒碰撞是指:攜帶顆粒的氣流在運動中遇到障礙物發(fā)生繞流,其中粒徑大于10μm的顆粒因動量充足而未跟隨流線繞流,繼而保持初始方向運動并撞擊接觸表面的過程。碰撞概率取決于顆粒運動軌跡,即受流場條件控制,可通過斯托克斯數(shù)St(表征顆粒慣性作用與擴散作用的比值)進行表述[15]。而黏附概率取決于顆粒的入射速度、角度、粒徑、碰撞損失以及溫度等多個因素。

    作為煤粉燃燒后的破碎產(chǎn)物,飛灰顆粒的形成主要來源于煤粉中無機組分(可溶性鹽、與有機物結合的成灰化合物、原生礦物質、外來礦物質)在燃燒過程中發(fā)生的一系列物理化學反應[13]。參考表1中飛灰的組分分析[16],可以看到其成分復雜多樣。此外,實際條件下的飛灰顆粒大多為不規(guī)則球狀,為便于深入探究顆粒沉積機理及臨界沉積理論,大多數(shù)學者會采取選用一些簡化后的標準顆粒(如SiO2、玻璃、銅等)開展研究。因此在本文的論述分析過程中,除飛灰顆粒外,也引入了部分標準顆粒作為對象。

    表1 飛灰顆粒全元素分析[16]

    與標準顆粒相比,飛灰顆粒不僅沉積類型較多,且材料物性參數(shù)不一,但從顆粒沉積過程中的形態(tài)變化來看,除上述提到的化學反應、氣相擴散等沉積類型外,其他沉積方式幾乎都可以通過顆粒碰動力學理論進行解釋[17]。因此,可以在不同材料的標準顆?;A上引入必要的附加力和能量損失,從而描述飛灰顆粒的沉積過程。依托研究主體,對涉及微顆粒運輸過程的碰撞概率不再作討論。圖2顯示了顆粒(非熔融態(tài)、熔融態(tài)、半熔融態(tài))慣性撞擊壁面后可能存在的5 種現(xiàn)象[3],在此僅對顆粒的黏附與反彈行為做主要分析。顆粒碰撞作為一個能量轉化與耗散的過程,關鍵在于入射動能與能量損失的相對大小,當入射顆粒碰撞后的動能不足以對抗碰撞過程的能量損失,顆粒就會發(fā)生黏附,否則反彈。這樣的碰撞結果常常采用臨界捕集速度和恢復系數(shù)進行描述[4],不僅能夠直觀地反映顆粒是否發(fā)生黏附,其中恢復系數(shù)還能夠度量碰撞能量損失。而能量損失受界面黏附力、黏彈性、彈性/塑性變形、溫度等多重因素影響[4,18-20]。

    圖2 固體顆粒撞擊壁面:可能的現(xiàn)象[3]

    恢復系數(shù)是描述物體碰撞過程中能量損耗程度的重要參數(shù)[21],其應用廣泛,有以下三種定義方式[22-24]:第一種,以碰撞后速度與碰撞前速度之比的牛頓定義;第二種,以卸載與加載階段沖量之比的Poisson 定義;第三種,以恢復階段彈性應變釋放的能量與壓縮階段儲存的能量之比的Stronge 定義。本文恢復系數(shù)基于1687年被提出的牛頓定義,即認為領域內的恢復系數(shù)e為顆粒碰撞后的反彈速度Vr與初始入射速度Vi的比值。當e=0 時,顆粒黏附;當0<e<1時,顆粒反彈;而e=1時,顆粒反彈,但此時屬于完全彈性碰撞,碰撞過程無能量耗散。而臨界捕集速度Vcr是指當顆粒與壁面發(fā)生碰撞時顆粒剛好發(fā)生反彈而非黏附在壁面的速度。為了便于理解和區(qū)分不同的判據(jù)標準,引入圖3對顆粒碰撞過程中涉及的參數(shù)進行說明,并給出了恢復系數(shù)、法向恢復系數(shù)和切向恢復系數(shù)的表達式,特別指出法向臨界捕集速度對應上述描述臨界捕集速度的法向分量。

    圖3 顆粒動力學示意圖

    而針對與已經(jīng)被接觸壁面所捕獲的顆粒,決定顆粒最終是否沉積還要取決于顆粒所受流體的剪切力與壁面之間黏附力的相對大小。當流體剪切流作用于顆粒的剪切力小于等于顆粒與壁面的黏附力的時候,才可確保顆粒繼續(xù)沉積[25],而黏附力受到范德華力、靜電力和毛細力的影響[26],力的大小與顆粒與壁面兩者的材料特性、帶電性以及濕度等有關[27]。

    2 基于沉積因素對細顆粒沉積判據(jù)的研究

    如何針對性地提出有效合理的沉積判據(jù)對鍋爐設計及后處理具有重大意義,而針對顆粒臨界沉積標準的研究,體現(xiàn)在提供合適的臨界黏附/反彈判據(jù),主要方法集中于實驗測量和理論數(shù)值分析兩方面,其中實驗研究通常是搭建微米級顆粒撞擊裝置,利用相關設備記錄顆粒碰撞前后運動行為參數(shù);理論數(shù)值計算主要是基于牛頓第二定律,借助顆粒動力學來研究顆粒的運動、受力及二者的相互作用,以此來定量描述并驗證實驗結論。相反,實驗也為顆粒碰撞的數(shù)值分析與模擬提供相應的數(shù)據(jù)支撐。實際碰撞過程中影響顆粒慣性碰撞行為的因素有很多,要想找到確定單一公式來對臨界捕集速度Vcr或恢復系數(shù)e進行有效預測也并非易事,本章描述了各類沉積因素對細顆粒沉積標準的研究及影響。

    2.1 入射速度對細顆粒沉積模型的影響

    基于離散相的飛灰顆粒較連續(xù)相氣體的密度高很多,且異相顆粒數(shù)量較少[28]。因此在實驗研究中,研究者常常忽略顆粒-顆粒之間的相互碰撞,基于單顆粒碰撞的入射速度出發(fā),嘗試提出臨界捕集速度或恢復系數(shù)等黏附判據(jù)。作為最直觀的沉積因素,衍生出了大量臨界沉積準則。

    顆粒黏附效率能夠借助顆粒的速度函數(shù)來間接測評,Dahneke[29-30]設計了一款顆粒光束發(fā)射裝置,用來測量顆粒轟擊平板表面的黏附效率,實驗發(fā)現(xiàn)顆粒存在某一特定速度(臨界捕集速度),但并未獲得入射速度在臨界捕集速度附近時的詳細反彈數(shù)據(jù)。針對這一問題,Wall等[31]對顆粒法向撞擊壁面進行了全面的研究,其中利用激光多普勒測速裝置(LDV)獲得了不同材料壁面下的顆粒法向恢復系數(shù)隨入射速度的變化規(guī)律,依托能量守恒原理得到了不同條件下的臨界捕集速度,為后續(xù)研究提供了廣泛的數(shù)據(jù)支撐。

    同Wall 等[31]一樣,Dunn 等[32-33]也先后采用了LDV 及類設備-粒子軌跡成像系統(tǒng)(PTIS),并與視頻分析系統(tǒng)結合使用,來確定粒子法向碰撞與斜向碰撞時的入射與反彈速度及其速度分量;不同的是,將單一目標顆粒的入射和反彈階段進行了一體化考慮,從而測量碰撞顆粒的恢復系數(shù),而并非是平均值[32]。由于Dahneke[29]、Wall等[31]中提出的簡單顆粒動力學及捕獲模型在顆粒反彈和捕獲階段的數(shù)據(jù)不完全合理化,為賦予顆粒碰撞結果更直觀的判定標準,在此基礎上,Dahneke[34]詳細地推導了各自模型對應黏附能增量的適應范圍,假設球體完全彈性,測量了不同入射速度下單分散聚苯乙烯乳膠小球對石英壁面的影響,建立了黏附模型并推導了臨界捕集速度和恢復系數(shù),但礙于模型中部分參數(shù)測量的局限性,限制了該判據(jù)的推廣。在黏附撞擊動力學模型[35](基于赫茲理論)的基礎上,Cheng等[36]進行了開發(fā),通過耦合計算流體力學(CFD)軟件,拓寬了應用范圍。依托EA 模型,Kim 和Dunn[37]考慮了彈性和黏附耗散的影響,由于實驗入射速度遠低于材料屈服速度,因此無需考慮塑性變形損失,根據(jù)法向入射速度和材料特性進一步推進了法向恢復系數(shù)的研究,相關解析解如式(1)所示(CA和CH為阻尼系數(shù);ψA和ψH為取決于組合剛度、粒子特性等因素的物理量);而后借助高速攝影機對低速域EA 模型(實線)驗證[38],數(shù)據(jù)(圓點)吻合良好,置信度95%,詳見圖4。

    圖4 法向恢復系數(shù)與法向入射速度的數(shù)據(jù)[38]

    深入研究發(fā)現(xiàn),上述實驗數(shù)據(jù)[37]均處于理論模型下方,即對應實驗臨界捕集速度更大。一些學者常常將這類臨界捕集速度實驗值偏高歸因于靜電效應或微凸處塑性變形等[8]。但已有學者排除接觸面粗糙度和電荷的影響,認為實驗中的臨界捕集速度偏高的結果來源于對材料中微觀的能量耗散認知不足[39]。因此,推測這種未知的能量耗散可能是黏彈性阻尼和表面黏附之間的耦合,能量耗散增加,從而提高了臨界捕集速度。

    不同于純彈性碰撞,Thornton和Ning[40]考慮了顆粒加載階段界面黏附效應,依托JKR理論提出了與顆粒速度相關的恢復系數(shù)解析解,如式(2)~式(4)。

    恢復系數(shù)隨歸一化入射速度的預測數(shù)據(jù)見圖5,其中橫坐標為入射速度與屈服速度之比,分析發(fā)現(xiàn)顆粒突破臨界捕集速度Vs(實線),也突破了材料屈服速度,產(chǎn)生的塑性變形帶來了較大的能量耗散,從而導致恢復系數(shù)極大值隨著臨界捕集速度與屈服速度之比的增大而減小,且遠低于單位值1。此類恢復系數(shù)將入射速度與屈服速度歸一化處理,通過簡單沖擊實驗參數(shù)表述,避免了一些材料特性(如屈服應力、界面能等參數(shù))難以測量的問題,非常便于開發(fā)為用戶自定義函數(shù)(UDF)并納入數(shù)值模擬當中,依托入射速度提供一個可變的恢復系數(shù),可以更好地體現(xiàn)基于能量耗散下的顆粒沉積分布。

    圖5 恢復系數(shù)與入射速度的預測數(shù)據(jù)[41]

    從上述Kim等[37]實驗臨界捕集速度偏大的推測中可以看出,迫切需要一種更加微觀的、能夠解決顆粒局部變形的方法。鑒于此類需求,有限單元法(FEM)得到了廣泛發(fā)展,基于模型定義在網(wǎng)格尺度上的優(yōu)勢,可以將黏彈性耗散與局部應力、應變相關聯(lián),還能夠對卸載問題進行有效處理。從而更真實地反映顆粒碰撞的黏彈性行為。在研究入射速度對顆粒反彈的影響時,Wu 等[42]利用FEM 發(fā)現(xiàn)針對微小塑性變形,顆粒的恢復系數(shù)主要受顆粒入射速度Vi與屈服速度Vy的比值主導,通過對比,這一點與Thornton 等[40]的沖擊力學理論預測基本類似。相關有限塑性變形下的恢復系數(shù)滿足式(5)關系(E*和Y分別為楊氏模量與屈服應力)。

    在對顆粒模型簡化的過程,與Thornton 等[40]恰恰相反,清華大學Feng等[43]忽略掉顆粒加載階段的黏附能,采用有限元法(FEM)研究了彈性及彈塑性顆粒在法向方向的碰撞行為,在彈性模型結合赫茲接觸模型與能量守恒定律求解臨界捕集速度、恢復系數(shù),如式(6)、式(7)所示(K為彈性系數(shù);w2為反彈階段黏附能)。

    便于對比,將式(6)、式(7)中的理論恢復系數(shù)同有限元結果一同繪制在圖6中,觀察發(fā)現(xiàn),除了在接近捕獲速度的低入射沖擊速度下,解析解與數(shù)值結果具有良好的一致性(黑色曲線為理論模型曲線,離散點為有限元結果)。不難看出,彈性顆?;謴拖禂?shù)隨著入射速度增大而增大并最終趨于穩(wěn)定,這樣的結果一方面證實了界面黏附能對于入射速度較小的顆粒影響更大,另一方面也符合彈性顆粒在模型簡化后的結果(忽略加載階段黏附能)。值得注意的是,在接近臨界捕集速度的低速段存在一定差異,考慮為理論模型中對于顆粒入射階段黏附能的忽略,導致相同條件下的顆粒,理論恢復系數(shù)較大,對應的臨界捕集速度較小。

    圖6 3種半徑的彈性顆粒的恢復系數(shù)[43]

    結合圖7中彈塑性顆粒碰撞,恢復系數(shù)隨入射速度增加呈先上升后下降的趨勢。速度作為影響顆粒沉積的主導因素,很顯然,此類材料下的顆粒在較高入射速度下(突破屈服速度)產(chǎn)生的塑性變形損失對顆粒行為影響巨大。此外,對于彈塑性顆粒而言,存在雙向耦合作用:黏彈性會增大塑性變形,相反塑性變形也會引起黏彈性損失增加。相關數(shù)值數(shù)據(jù)誤差驗證的過程可參考圖8,同Wall 等[31]實驗數(shù)據(jù)吻合良好。

    圖7 3種半徑的彈塑性顆粒的恢復系數(shù)[43]

    圖8 基于恢復系數(shù)的數(shù)據(jù)誤差分析[43]

    本節(jié)判據(jù)準則中,不僅涉及了顆粒彈性沖擊和黏附力,而且涵蓋了塑性變形損失,并基于黏彈性材料本構關系復雜性給定了相應的解決途徑。鑒于飛灰的復雜程度,在沉積準則的探索中引入能量損失是十分必要的,可以顯著提高沉積模型的準確性與易用性。結合表2來看,一些顆粒模型并未引入黏彈性損失,這會導致判據(jù)中的臨界捕集速度低于飛灰實際值,具有局限性,因此需要尋求結合彈性、黏附力、黏彈性以及塑性變形損失等更為全面的顯式解析解更具切實意義。

    表2 臨界捕集速度和恢復系數(shù)的經(jīng)典解析式

    2.2 入射角度對細顆粒沉積模型的影響

    顆粒入射角度會影響顆粒的運動狀態(tài),從而進一步影響顆粒的沉積狀況,大多數(shù)研究主要集中于顆粒的法向碰撞,而實際情況下顆粒與顆粒之間或者顆粒與壁面的碰撞多為斜向碰撞,即入射角度與接觸表面會存在一個非零碰撞角,除正常黏附或反彈外,斜向碰撞可能會引起顆粒發(fā)生滑動或者滾動。然而,法向臨界捕集速度不考慮滑動、滾動行為,以此簡化標準,用來判定顆粒的黏附或反彈行為,但這樣的標準太過于片面,尤其針對高入射角的碰撞行為無法精準衡量[4]。與常規(guī)法向碰撞相比,斜向碰撞并不那么便于解析,其顆粒黏附/反彈標準也存在差異,因此很多學者基于入射角度選取恰當?shù)酿じ浇佑|力學模型對顆粒運動方程進行數(shù)值分析,旨在尋求一個切實的顆粒黏附/反彈判據(jù)標準。

    在研究氣固兩相流中顆粒對渦輪葉片的碰撞行為時,有人提出第一個顆粒碰撞回彈模型,認為入射角度是影響顆?;謴拖禂?shù)的重要因素[48]。在此基礎上,為進一步深入研究反彈特性,選取飛灰顆粒作為研究對象,基于3種不同壁面的反彈行為進行了實驗[49],利用LDV獲得各種入射角度下的恢復系數(shù),見圖9。

    圖9 顆粒速度恢復比的比較[49]

    以2024 鋁合金為例,對顆粒恢復系數(shù)平均值的最小二乘多項式曲線擬合出如式(8)關聯(lián)式。

    對應關聯(lián)式不同角度下的顆粒法向/切向標準差分析可參考表3。圖9 中顆粒入射角在0°~75°范圍內的恢復系數(shù)隨著入射角增大逐漸降低,但本文中是否存在臨界黏附角度還有待進一步深究。不過可以確定的是,除渦輪葉片上沉積的顆粒外,入射角對顆?;謴拖禂?shù)的影響巨大,這一結論對于鍋爐中飛灰同樣適用。

    表3 不同入射角下顆粒對應的恢復系數(shù)標準差[49]

    PIV是一種多點、瞬態(tài)流動測試技術,與上述LDV 相比,具備相當優(yōu)異的空間分辨率。為進一步表征顆粒碰撞規(guī)律,Cai 等[50]將恢復系數(shù)分為法向恢復系數(shù)與切向恢復系數(shù),基于PIV技術獲取了氧化鐵顆?;謴拖禂?shù)隨入射角度的變化規(guī)律,運用最小二乘法擬合了恢復系數(shù)分量的修正平均值和相應入射角(加權平均)的關聯(lián)式如式(9)、式(10)。

    上述研究工況及考慮因素的局限性也直接導致了該擬合判據(jù)的局限性,此外受流場加速的影響,法向恢復系數(shù)會降低3%~6%,測量可能會產(chǎn)生較大誤差。通過部分工況下的實驗數(shù)據(jù)總結實驗規(guī)律,從而擬合恢復系數(shù)與入射角度的關系,但是此類方案適用性有限,難以保證擬合模型的通用性。故Chen 等[4]利用DEM 結合黏性JKR 模型探究了微米顆粒在不同角度下的碰撞,分析了入射角對黏附判據(jù)的影響,具體結果見圖10。在此之前,Dominik等[51]證明了滾動與滑動行為在顆粒碰撞作用中可以忽略,認為法向臨界捕集速度VCN會存在固定值,很顯然,該結論與圖10中VCN曲線十分吻合,可以看到,在入射角低于65°時,法向臨界捕集速度VCN幾乎恒定;不過值得注意的是,當入射角度大于突變角度,VCN便迅速下降,這進一步證實了Dominik等[51]的結論存在適用范圍,即小入射角范圍內奏效,利用法向碰撞的臨界捕集速度VC0與回彈階段δN=0 處的臨界速度VC2之間的差值[式(11)]來表征黏彈性效應,見圖11,將入射速度Vin替換上式VC0,獲取Vin的函數(shù),采用圖解法求解法向碰撞下的臨界捕集速度VC0(雙線交點)。

    圖10 臨界捕集速度VC及法向臨界捕集速度VCN與入射角θ的函數(shù)關系[4]

    圖11 VC0 - VC2和HV i4n 5作為入射速度Vin的函數(shù)關系[4](插圖表示VC0作為耗散系數(shù)α的函數(shù)的變化)

    由于斜向碰撞求解方程非線性的性質,直接推導較大入射角下的臨界捕集速度判據(jù)不太可能,還會出現(xiàn)顆粒入射速度小于臨界捕集速度,但并未發(fā)生黏附的不合理狀況。不過可以結合圖10 的DEM結論進行數(shù)值轉換,獲得法向臨界捕集速度的冪律,結合圖10 插圖所示,采用雙曲線擬合相對偏差ΔV*CN和1 -b?,得到?jīng)_擊參數(shù)b?下法向臨界速度VCN的顯式表達式如式(12)、式(13)(其中b?90為特征角對應的沖擊參數(shù))。通過引入尺寸比,結合撞擊參數(shù)提出了全方位角度的顆粒黏附/反彈判據(jù)表達式,如式(14)。

    在探究SiO2顆粒慣性撞擊鋼板的實驗中,Xie等[52]基于動力學模型探究了入射角度對反彈特性的影響,根據(jù)歸一化能量損失(滑動損失)TKL和圖12來看,歸一化能量損失與入射角密切相關,且隨著入射角呈現(xiàn)先上升后下降的趨勢,但下降速度十分平緩,由此可以認為較高的入射角度更有利于顆粒發(fā)生黏附。

    圖12 不同速度下歸一化能量損失隨入射角的變化[52]

    切向反彈速度利用單/雙系數(shù)模型對兩個不同的動力學模型進行預測,依托單系數(shù)模型為例,將實驗數(shù)據(jù)計算出的阻尼系數(shù)與相應實驗恢復系數(shù)數(shù)據(jù)進行多項式擬合,如式(15)。

    當入射速度小于臨界捕集速度,法向恢復系數(shù)en=0,引入阻尼系數(shù)至離散運動方程[53],即可通過數(shù)值計算的方法獲得臨界捕集速度。通過對比單、雙系數(shù)模型結果,發(fā)現(xiàn)除極值外,兩個模型的切向力與位移計算結果相似,不過兩個矢量基于顆粒高入射角下的變化,均表現(xiàn)出了對入射角度的強烈依賴性。因此,為進一步預測不同入射角度下的顆粒碰撞行為,利用高速攝像機記錄了煤灰顆粒的撞擊過程,得到了法向恢復系數(shù)隨入射法向速度的變化規(guī)律[54],也進一步驗證了顆粒的臨界捕集速度隨著入射角的增大而減小。

    2.3 粒徑尺寸對細顆粒沉積模型的影響

    入射顆粒的尺寸會引發(fā)顆粒慣性和黏附力的相對重要性發(fā)生變化。隨著顆粒尺寸逐漸增大,影響顆粒行為的因素會從黏附主導轉移到慣性主導狀態(tài),從而進一步導致不同的碰撞結果[55]。此外,已有多篇研究[4,6,45-46,56]指出尺寸差異也會破壞原有機制引起的能量耗散的分數(shù)。

    許多學者已經(jīng)進行了不同尺寸下顆粒沉積行為的研究,浙江理工大學季浪宇[57]對大尺寸(3~7mm)顆粒進行碰撞反彈的研究中發(fā)現(xiàn),顆粒的尺寸不會對恢復系數(shù)產(chǎn)生影響。與之相反,Sommerfeld 等[58]研究發(fā)現(xiàn)當顆粒直徑較小時(100~500μm),顆粒尺寸對恢復系數(shù)影響顯著,且粒徑越大,恢復系數(shù)越小,但影響程度隨表面粗糙度增加而降低。

    上述結論中顆粒尺寸處于毫米或次毫米級,而對于微米級顆粒而言,其碰撞結果和判定標準截然不同。Rogers 等[59]針對玻璃、銅及不銹鋼3 種材料在不同組合下的碰撞過程進行了測量,并獲得相應材料不同粒徑下的臨界捕集速度,隨著顆粒尺寸增加,各類顆粒的臨界捕集速度均呈現(xiàn)不同情況的下跌,見圖13 的▲和★。另外,圖13 中的▼為韓健等[16,46]以電廠中飛灰顆粒為研究對象進行了慣性撞擊鋼板的實驗,基于JKR理論[40]推導的臨界捕集速度如式(16)所示(其中γ為界面能)。

    圖13 臨界速度與顆粒大小的可用數(shù)據(jù)的比較[46]

    與玻璃、銅顆粒相比,其臨界捕集速度隨尺寸的變化規(guī)律表現(xiàn)為良好一致性,這也進一步驗證了上述的猜想:微米級顆粒隨尺寸的判定標準與大尺寸顆粒(毫米級)不同。不過需要注意的是,對于尺寸越小的顆粒,首次碰撞能量損失對顆粒行為的主導作用越大。式(16)中臨界捕集速度的表達式與顆 粒 尺 寸 滿 足 冪 律Vcr∝Rα的 關 系,這 與Kim 和Dunn[38]的研究中簡化阻尼系數(shù)后提出的臨界捕集速度形式相似,均利用顆粒尺寸進行了表征,如式(17)所示(K為組合剛度、CR與黏附力相關)。

    正如預期那樣,Jassim[45]、Dong 等[46]也先后在研究中發(fā)現(xiàn),伴隨著顆粒尺寸增大,相應的臨界捕集速度呈指數(shù)形式下降,滿足上述冪律關系。

    僅探究微米級顆粒對光滑壁面的碰撞是遠遠不夠的,基于顆粒覆層的碰撞研究更具工程指導意義,但與此同時,撞擊覆層的機制也更為復雜。因此,Zhu等[60]通過DEM研究了顆粒和覆層碰撞后的黏附/反彈行為,模擬不同入射顆粒尺寸與覆層顆粒尺寸對平均恢復系數(shù)的影響,同Chen等[4]做法類似,將入射顆粒與覆層中顆粒的尺寸歸一化,即引入顆粒尺寸比(入射顆粒直徑與覆層中顆粒直徑的比值為dpi/dpp),結合圖14、圖15 發(fā)現(xiàn),顆?;謴拖禂?shù)隨尺寸比增大而增大,由于覆層結構能量耗散的復雜機制,結合3個可控參數(shù)的調整對恢復系數(shù)進行曲線擬合如式(18)所示。同DEM 模擬數(shù)據(jù)比較,擬合相關性良好,可參考文獻[60]。

    圖14 不同入射顆粒直徑與平均恢復系數(shù)作為粒子入射速度的函數(shù)[60]

    圖15 不同填料顆粒粒徑的平均恢復系數(shù)與粒子入射速度的函數(shù)[60]

    式中,a為穩(wěn)定恢復系數(shù)相關;b為當恢復系數(shù)達到其穩(wěn)定值的一半時的入射速度;c為與快速增長區(qū)域的斜率有關。

    基于顆粒碰撞接觸的時間與連續(xù)性劃分為離散模型和連續(xù)相模型[61],離散模型主要基于剛體材料,而連續(xù)相模型則能夠詳細地描述實際碰撞行為(接觸力與變形),如彈簧-阻尼模型[62]、Stronge 模型[24,63]以 及J-G 模 型[64]。Yu 等[65]基 于Stronge 模 型[24]建立了一種修正型的細顆粒碰撞模型,將黏附力以及粗糙度因子相結合考慮得到黏附功,利用動摩擦系數(shù)μ與法向接觸壓力建立了顆粒的切向恢復系數(shù)模型,如式(19)。

    在后續(xù)模型驗證過程中,基于4種不同案例對改進后的模型進行誤差分析,較初始模型[24,64]而言,改進后的模型對恢復系數(shù)的預測更為精確,平均誤差不高于9%,遠遠低于初始模型的平均誤差(詳細誤差分析可參考文獻[65])。結合圖16 分析(2024-鋁合金,粗糙度因子CR=1.0),不難發(fā)現(xiàn),尺寸小的顆粒更加具備發(fā)生黏附的條件,不過一旦發(fā)生反彈,考慮為小顆粒具備較大的屈服強度和楊氏模量,對應塑性變形引發(fā)的能量耗散也較小,因此在相同工況下,小顆粒往往比大顆粒具備更大的恢復系數(shù)。

    圖16 不同直徑顆粒碰撞的法向恢復系數(shù)曲線[65](粗糙因子為1.0)

    通過對基于顆粒尺寸對顆粒沉積模型影響的文獻總結分析,能夠發(fā)現(xiàn)微米級顆粒和大尺寸顆粒對應的沉積規(guī)律及判定標準不能夠相提并論,也進一步強調了顆粒尺寸對微米級顆粒黏附/反彈標準的影響不容忽略。結合表2,觀察發(fā)現(xiàn)微米級顆粒臨界捕集速度對顆粒尺寸具有冪律依賴性[38,40,43,45-46],此外,由于顆粒尺寸較小,不同工況下的顆粒受動力學及工況的影響,會使得顆粒碰撞行為差異顯著。覆層的存在使得顆粒在能量耗散方面占據(jù)重要因素,由于覆層顆粒尺寸與位置不一,還需要充分考慮到顆粒碰撞過程隨機性。此外,分析顆粒尺寸變化對碰撞結果造成的顯著影響可利用黏性微米級顆粒的三維離散元法[4],有助于探究涵蓋所有角度的顆粒黏附/反彈準則表達式。

    2.4 溫度對細顆粒沉積模型的影響

    在上述的論述研究中,針對微顆粒沉積準則的探討已經(jīng)做了大量工作,但此類沉積模型大多致力于恒溫條件下的應用,而飛灰顆粒的運行常常伴隨溫度變化,鍋爐中煙氣、飛灰及換熱表面都會彼此換熱[66],更重要的是,飛灰顆粒的物性參數(shù)也會隨溫度發(fā)生改變[67]。通過對上述顆粒沉積判據(jù)縱向分析,不難發(fā)現(xiàn),改變單一沉積因素(如改變顆粒尺寸、入射速度、入射角度),其黏附標準受單一沉積因素的影響效果幾乎相同。以入射角度為例,高溫下的顆粒往往不受撞擊角度的控制,此時溫度對于顆粒碰撞結果占據(jù)主導作用。因此,顆粒沉積判據(jù)的研究同溫度的相關性密不可分。

    在沉積過程中,環(huán)境溫度通過影響顆粒及接觸面材料的楊氏模量、表面能等參數(shù)進而影響碰撞結果。從這個角度來看,獲取不同溫度下材料的物性參數(shù)成為了探究沉積準則的首要任務。Ai[68]提出了一種計算飛灰彈性模量的方法,為了填補Brach 模型[44]的缺失,Ai 等[69]在El-Batsh 模型[70]的基礎上,進一步調整了對楊氏模量的計算方案。高溫下的顆粒,黏附行為通常會由有效黏度決定,常規(guī)的臨界速度模型或許不再適用,因此一些學者開拓了經(jīng)驗模型進行評估。Sreedharan 等[71]指定臨界黏附溫度為顆粒熔融溫度,當溫度大于臨界黏附溫度時,顆粒黏附概率為1,而溫度遠低于臨界溫度時黏附概率為0。定義臨界溫度下黏度與當前溫度下的黏度比值為顆粒沉積的概率,如式(20);其黏度μTp與溫度關系如式(21)[72]。

    其中顆?;瘜W成分決定A、B大小。從顆粒黏度的角度探究溫度對顆粒黏附的影響還有很多,但是否發(fā)生黏附的影響因素是材料的黏彈性變形和黏附力,不能直接判定這兩種因素和材料黏度存在某種關系。

    臨界黏度模型僅僅考慮到溫度的影響,并沒有把顆粒速度這一影響因素考慮進去;而臨界速度模型考慮到溫度對材料彈性模量的影響,并沒有把溫度對顆粒黏度這一影響因素考慮進去,兩者都存在一定的缺陷。圖17 分別顯示了葉片臨界黏度模型和臨界速度模型的沉積圖,通過云圖對比發(fā)現(xiàn),除前葉板外,兩種模型的沉積分布與沉積總量均差異較大[73]。

    圖17 沉積濃度云圖[73]

    針對上述模型缺點,在此基礎上,裴鈺[74]耦合出了新的沉積概率模型,綜合考慮溫度場與速度場對粒子沉積行為的影響。耦合思想為:當前溫度Tp大于等于a倍臨界黏附溫度Tcrit時,顆粒為熔融態(tài),黏附概率為1;當前溫度Tp小于等于b倍臨界黏附溫度Tcrit時,顆粒為固體,黏附概率為0;而當前溫度介于兩者之間時,顆粒沉積概率如式(22)所示(文獻[74]選取a=1、b=0.5,具體取值因實驗修正而異)。

    為深入了解溫度對微顆粒沉積的影響,一些學者往往會選取一些物性參數(shù)受溫度影響比較大的材料著手。孫奇[47]將聚苯乙烯顆粒作為研究對象,通過表征彈性位移與黏附力作用特征距離的比值,確定了JKR模型的適用性,并發(fā)現(xiàn)單分散聚苯乙烯顆粒首次碰撞能量損失隨著溫度升高而增大,臨界捕集速度也增大。值得一提的是,與飛灰顆粒相比,由于聚苯乙烯顆粒在沉積過程中的楊氏模量受溫度影響很大,便于探究溫度對顆粒黏附規(guī)律的影響,但從本質上來講,飛灰形貌不一,成分復雜,且不同煤種的組分含量差別巨大,對于物性參數(shù)的測量也是非常不易的,除此之外,飛灰顆粒沉積過程中楊氏模量隨溫度變化較小,且熔點較高。因此可以盡量選用低熔點飛灰顆粒,區(qū)別不同溫度區(qū)間下的黏附規(guī)律??偠灾侠聿牧系倪x取對于推進高溫沉積判據(jù)的發(fā)展是非常有利的。溫度作用于沉積判據(jù)的影響不僅基于顆粒,Abd-Elhady等[75]發(fā)現(xiàn)溫度梯度引起沉積層發(fā)生燒結,而燒結后沉積層的硬度更強(孔隙率降低),而對于較硬平面的顆粒撞擊,粒子的動能損失很小,傾向反彈,相應的顆粒沉積概率降低,直至沉積層不再生長,對應臨界捕集速度也隨之降低,然后趨于穩(wěn)定。

    截止目前,無論常溫還是高溫領域,JKR模型都是作為首選對象,通過引入阻尼系數(shù)來推測相應的黏附接觸模型,不過,在熔融溫度下的飛灰顆粒會發(fā)生相態(tài)變化,必要時需要引入液橋力,對黏附接觸模型進行更深度的改進。另外,顆粒覆層作為一個龐大的沉積體,會發(fā)生擠壓變形、燒結、脫離等行為,受沉積分布的影響,不同區(qū)域的沉積程度也不盡相同,壁面的受熱不均使得溫度梯度出現(xiàn)巨大差異,因此不易掌握沉積體的整體溫度分布,這一問題對于探究溫度下顆粒沉積準則意義重大。需要進一步優(yōu)化和改進模型,調整實驗方案,實現(xiàn)量化研究,以便提出恰當?shù)陌虢?jīng)驗沉積準則。

    3 基于起動方式對顆粒臨界剝離判據(jù)的研究

    前文已經(jīng)分析了多種因素(入射速度、入射角度、粒徑尺寸及溫度)對顆粒沉積模型的影響,然而對于一個完整的積灰模型,除沉積模型外,還應包括剝離模型,兩者相應共生,不可分割[76-77]。當目標顆粒被壁面捕獲,盡管在短期內處于黏附狀態(tài),但自身仍處于一個復雜環(huán)境下,隨時可能因氣流剪切力或外部顆粒注入等因素發(fā)生再次脫離[78],這一剝離行為對沉積質量的影響不容忽略。因此,積灰剝離是積灰過程中一個關鍵環(huán)節(jié),根據(jù)顆粒起動原理,可以分為非接觸起動和接觸起動。非接觸起動是指當顆粒受到流體的剪切力作用,使得顆粒發(fā)生移動并脫離沉積表面的過程,起動方式主要包括滑動、滾動及上升。而接觸起動是指顆粒依靠懸浮顆粒碰撞而被去除的過程。在非接觸起動的研究中發(fā)現(xiàn),滾動是顆粒去除的主導機制[79]。

    3.1 非接觸起動

    由于非接觸起動主要方式為滾動,滾動去除力來源于流體動力矩,當流體動力矩大于黏附力矩,顆粒就發(fā)生滾動去除。因此可以采用臨界力矩理論對顆粒的剝離模型進行分析。為了研究流體速度對廢氣再循環(huán)(EGR)冷卻器沉積顆粒的影響,Abd-Elhady 等[80]基于Zhang 等[81]定義的滾動力矩比對處于剪切流中的水平顆粒受力進行分析,如式(23)。

    式中,F(xiàn)a為范德瓦爾斯黏附力;Fl為升力。當RM>1,顆粒發(fā)生滾動;當RM=1,流體流速為臨界滾動流速。通過改變氣流方向來實現(xiàn)沉積顆粒的去除,朝下流動是通過外部注入顆粒去除沉積的最佳流向[82]。Tang等[83]基于能量守恒定律與臨界力矩理論,開發(fā)了顆粒沉積過程的去除算法,提出了流速和顆粒尺寸對沉積特性的影響,但并未進一步考慮結構幾何參數(shù)對換熱面沉積過程的影響[84]。圖18為處于流場中黏附顆粒的受力分析圖,能看出黏附顆粒受阻力FD、黏附力Fa、升力Fl和流體動力矩MD組成。

    圖18 黏附顆粒受力分析[83]

    結合臨界力矩理論,能夠很好地理解顆粒去除的滾動機理,以“0”點為滾動支點的顆粒所受流體動力矩大于顆粒黏附力力矩時(圖18),便發(fā)生滾動效應,用不等式(24)表示。

    其中升力對顆粒的影響可以忽略不計[76],利用JKR理論計算黏附力和接觸半徑,代入臨界力矩模型,獲得臨界剪切速度,如式(25)。

    同樣也是粒子能夠繼續(xù)保持黏附在壁面上的臨界捕集速度,其中C為修正系數(shù)。

    在Tang 等[83]研究的顆粒臨界力矩模型的基礎上,Tong等[85]添加了顆粒斜向碰撞下產(chǎn)生的沖擊力Fc,耦合有限體積法(FVM)和格子玻爾茲曼法(LBM)探索管道上顆粒沉積機制與剝離機制。這也就意味著該臨界剝離模型中的顆粒處于非靜止的動態(tài)過程,而非捕獲顆粒的剝離標準。結合圖19進行分析,當流體與沖擊力的合力矩大于黏附力時,顆粒發(fā)生滾動脫離,如式(26)。

    圖19 作用在附著粒子上的力和力矩[85]

    當阻力與沖擊切向力合力大于摩擦力時,顆粒發(fā)生滑動脫離,如式(27)。

    剪切速度定義為式(28)。

    上述表達式經(jīng)計算驗證,同Toscano 等[86]研究數(shù)據(jù)吻合良好。且臨界滑動分離的剪切速度大于臨界滾動分離的剪切速度,并隨著黏附顆粒和沖擊顆粒直徑的增大臨界剪切速度減小。表4歸納了部分顆粒起動脫離研究。

    表4 顆粒去除行為研究歸納

    3.2 接觸起動

    相較無顆粒注射,存在外部顆粒注射的臨界剪切速度要遠遠小于沒有顆粒的時候,進一步驗證了顆粒注射是主要去除機制[85]。不過,通過結合前面Abd-Elhady等[82]的研究內容分析發(fā)現(xiàn),雖然通過外部顆粒注射為沉積去除主要機制,但是若流動方向不佳,是無法單憑外部注入顆粒來完成沉積顆粒的去除,僅在非接觸力的流體加持下,外部注射顆粒的接觸力起動才占據(jù)主導機制。

    基于上述顆粒去除機制的研究分析,Tang等[83-84]和Tong等[85]應用了各種模型來預測積灰顆粒的去除,但并未考慮沉積層的表面特性對沉積的影響,如表面顆粒尺度[88]、表面粗糙度[89]、溫度以及沉積形態(tài)[87]等?;诖祟惾毕荩琙heng 等[88]提出了滑動和滾動顆粒分離的新準則,為全面深入地展示沉積細節(jié),選取了更適用常溫下的臨界速度模型,因此忽略了溫度影響,由于實際的顆粒撞擊會出現(xiàn)多粒子碰撞,因此引入了修正系數(shù)下的沖擊力公式,如式(29)。

    修正系數(shù)取決于沉積顆粒的位置與尺寸,根據(jù)方程發(fā)現(xiàn)沖擊力與顆粒速度成正比,即氣流速度越大,越容易導致沉積顆粒的去除。通過比較有無去除模型下的灰沉積情況,發(fā)現(xiàn)去除模式下的沉積總量與Tang等相似,且呈現(xiàn)漸進趨勢[82-83],但兩種情況下的沉積量隨時間推移差異逐漸顯著。造成這一現(xiàn)象的原因可能是由于隨著沉積的進行,沉積層表面的填充顆粒粒徑也隨之變化,相較于細小顆粒,大顆粒更傾向于剝離而無法繼續(xù)黏附,此時受入射顆粒碰撞影響,被捕獲的的大顆??赡芤蚺鲎捕l(fā)生去除。

    從理論上來看,通過控制氣流速度確實可以有效降低沉積速率,不過需要注意的是,氣流速度雖然影響了顆粒剝離機制,但同時也改變了顆粒注入的質量流量以及沉積水平。因此顆粒最終的沉積速率并非是氣流速度的單一線性關系,這一點在探究模型結構對沉積速率的數(shù)值模擬中應用較多[78,90],需要進行耦合分析,根據(jù)不同工況尋求一個合理的速度區(qū)間。此外,剝離模型耦合數(shù)值模擬對吹灰器效率的評估也得到了較為廣泛的發(fā)展[91]。

    4 結語

    顆粒沉積的黏附機理與沉積判據(jù)對于鍋爐參數(shù)的設計起著關鍵作用。本文以細顆粒慣性撞擊為主導機制,基于臨界捕集速度、恢復系數(shù)、能量守恒定律、臨界沖擊角以及臨界力矩理論等模型探討預測了顆粒的運動行為。然后通過對比不同沉積判據(jù)的數(shù)學描述,指出各類模型的聯(lián)系與改進方向,為有效抑制飛灰沉積提供了一定的理論參考?;诋斍把芯繜狳c及趨勢,認為顆粒沉積研究在以下方面仍需進一步深入。

    (1)基于慣性碰撞的顆粒沉積的分析是一個復雜的動力學過程,除材料本身屬性外,還涉及入射參數(shù)、溫度、流場等外來因素,因此針對不同工況下顆粒沉積的判定依據(jù)也不盡相同。在實際環(huán)境下,多為顆粒與沉積體(覆層)的碰撞,由此導致?lián)Q熱表面的熱阻會大大增加,使得換熱面的溫度梯度與溫度分布發(fā)生劇烈變化,進而影響材料的物性參數(shù)而改變沉積結果。因此,如何獲取一個合理的沉積體溫度場至關重要。與此同時,也對作用于溫度效應下的顆粒沉積行為判定提出了挑戰(zhàn)。

    (2)針對高溫顆粒沉積領域,顆粒黏附與黏度存在經(jīng)驗關系,這對于部分工程實踐有一定的指導性,但無法作為黏附接觸模型的本質。另外,溫度直接導致的材料楊氏模量、表面能等參數(shù)的變化也間接影響了顆粒碰撞黏附力和黏附功,需要對具體的顆粒黏附過程進行具體分析。若涉及顆粒液膜生成,有需要的情況下可引入液橋力等來進行修正。目前缺乏高溫顆粒沉積的量化研究,可基于能量守恒定律進一步完善與深入,結合溫度提出不同溫度區(qū)間內(固態(tài)、熔融態(tài))的沉積標準。

    (3)在研究顆粒沉積過程中,大多集中于利用高速攝相機結合理論數(shù)值計算來吻合實驗數(shù)據(jù),并且也通過搭建試驗臺進行了大量的顆粒碰撞實驗,但依然缺乏數(shù)值模擬方面的研究,目前CFD-DEM模塊以FLUENT-DEM 為主線,也可通過耦合開源軟 件OpenFOAM 與LAMMPS 編 譯 的SediFoam 求 解器,其可擴展性與適用性也十分優(yōu)異。

    (4)針對顆粒動力學的實驗研究中,通常將固態(tài)顆粒視為稀疏相,因此大多采用氣流輸送單顆粒進行入射碰撞分析,這也就忽略掉粒子運輸過程中的碰撞行為,使得推導的顆粒黏附判據(jù)在工程的應用中存在相對誤差。

    符號說明

    CR——接觸粗糙度因子

    D——特征長度,m

    dpi(dp) ——入射顆粒直徑,m

    dpp——覆層顆粒直徑,m

    dδ——接觸直徑,m

    E——彈性模量,Pa

    E*——等效彈性模量,Pa

    e——恢復系數(shù)

    eavg——平均恢復系數(shù)

    en——法向恢復系數(shù)

    et——切向恢復系數(shù)

    Fa、Fpo——黏附力,N

    Fb——浮力,N

    Fc——沖擊力,N

    Fg——重力,N

    Fl——升力,N

    Ft、FD——曳力,N

    K——彈性系數(shù)

    m——質量,kg

    m*——有效質量,kg

    Ps——沉積概率

    R(Rp) ——半徑,m

    R*——有效半徑,m

    Re——雷諾數(shù)

    Rs——無黏附恢復系數(shù)

    St——斯托克斯數(shù)

    T——溫度,K

    TKL——歸一化能量損失

    u* ——臨界剪切速度,m/s

    uˉ——平均速度,m/s

    Vcr——臨界捕集速度,m/s

    Vin(Vi) ——顆粒入射速度,m/s

    Vn——法向入射速度,m/s

    Vr——顆粒反彈速度,m/s

    Vt——切向入射速度,m/s

    Vy——屈服速度,m/s

    ν——運動黏度,m2/s

    β——尺寸比

    β1——入射角(面夾角),(°)

    βˉ——平均入射角,(°)

    γ——界面能,J/m2

    η——阻尼系數(shù),kg/s

    ρ——密度,kg/m3

    ω——黏附功,J/m2

    μ——摩擦系數(shù)

    μcrit——臨界黏度,Pa·s

    μTp——顆粒黏度,Pa·s

    δ——接觸位移,m

    τw——壁面剪切應力,Pa

    猜你喜歡
    恢復系數(shù)法向入射角
    剛體彈性碰撞中恢復系數(shù)的探討
    一般三棱鏡偏向角與入射角的關系
    落石法向恢復系數(shù)的多因素聯(lián)合影響研究
    利用恢復系數(shù)巧解碰撞問題
    預制圓柱形鎢破片斜穿甲鋼靶的破孔能力分析*
    用經(jīng)典定理證明各向異性巖石界面異常入射角的存在
    低溫狀態(tài)下的材料法向發(fā)射率測量
    落石碰撞法向恢復系數(shù)的模型試驗研究
    不透明材料波段法向發(fā)射率在線測量方法
    一維對心碰撞恢復系數(shù)的討論
    物理與工程(2014年5期)2014-02-27 11:23:18
    国产真实乱freesex| 国产精品永久免费网站| 久久久a久久爽久久v久久| 一级毛片aaaaaa免费看小| 中文亚洲av片在线观看爽| 日韩一本色道免费dvd| 免费av毛片视频| 天美传媒精品一区二区| 久久久国产成人免费| 国产在视频线在精品| 国产在线男女| 国产午夜福利久久久久久| 男女下面进入的视频免费午夜| 国产成人精品婷婷| 日日摸夜夜添夜夜爱| 最近的中文字幕免费完整| 久久精品国产鲁丝片午夜精品| 国产单亲对白刺激| 赤兔流量卡办理| 久久久久国产网址| 3wmmmm亚洲av在线观看| 18禁在线播放成人免费| 有码 亚洲区| 国产亚洲精品久久久com| 成人鲁丝片一二三区免费| 赤兔流量卡办理| av天堂中文字幕网| 亚洲精品日韩av片在线观看| 99九九线精品视频在线观看视频| 国产老妇伦熟女老妇高清| 免费看av在线观看网站| 黄色欧美视频在线观看| 久久久久久久久久久免费av| 亚洲欧美日韩无卡精品| 爱豆传媒免费全集在线观看| 国产真实乱freesex| 午夜精品国产一区二区电影 | 久久久精品大字幕| 国产伦精品一区二区三区四那| 少妇的逼水好多| 久久久久免费精品人妻一区二区| 亚洲一区二区三区色噜噜| 精品国产三级普通话版| 国模一区二区三区四区视频| 中文字幕免费在线视频6| 欧美+日韩+精品| 亚洲欧美精品自产自拍| 国产蜜桃级精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 国产麻豆成人av免费视频| 亚洲美女搞黄在线观看| 伦理电影大哥的女人| 男人舔奶头视频| 国产一区二区亚洲精品在线观看| 春色校园在线视频观看| 如何舔出高潮| 赤兔流量卡办理| 边亲边吃奶的免费视频| 国产成年人精品一区二区| 校园人妻丝袜中文字幕| 亚洲av不卡在线观看| 免费一级毛片在线播放高清视频| 熟女电影av网| 九色成人免费人妻av| 亚洲五月天丁香| 国产黄色视频一区二区在线观看 | 十八禁国产超污无遮挡网站| 人妻制服诱惑在线中文字幕| 久久欧美精品欧美久久欧美| 久久久精品94久久精品| 人人妻人人看人人澡| 国产老妇女一区| www.色视频.com| 在线观看美女被高潮喷水网站| 精品欧美国产一区二区三| 真实男女啪啪啪动态图| 国产熟女欧美一区二区| 婷婷精品国产亚洲av| 在线天堂最新版资源| 国产毛片a区久久久久| 大又大粗又爽又黄少妇毛片口| 亚洲丝袜综合中文字幕| 色5月婷婷丁香| 岛国毛片在线播放| 色5月婷婷丁香| 久久久久免费精品人妻一区二区| 伦理电影大哥的女人| h日本视频在线播放| 国产精品.久久久| 成年女人看的毛片在线观看| 亚洲自拍偷在线| 91麻豆精品激情在线观看国产| 看非洲黑人一级黄片| 欧美激情在线99| 欧美色视频一区免费| 亚洲三级黄色毛片| 天堂中文最新版在线下载 | 亚洲无线在线观看| 精品欧美国产一区二区三| 国产成人a区在线观看| 国产精品.久久久| 亚洲美女视频黄频| 亚洲国产精品成人综合色| 麻豆精品久久久久久蜜桃| 国产伦理片在线播放av一区 | 国产精品爽爽va在线观看网站| 好男人在线观看高清免费视频| 色尼玛亚洲综合影院| 三级男女做爰猛烈吃奶摸视频| 国内揄拍国产精品人妻在线| 欧美zozozo另类| 久久精品国产亚洲av天美| 在线观看66精品国产| 国产高清视频在线观看网站| 欧美另类亚洲清纯唯美| 亚洲高清免费不卡视频| 国产高清有码在线观看视频| 亚洲国产精品成人综合色| 久久久精品欧美日韩精品| 美女cb高潮喷水在线观看| 99热这里只有是精品在线观看| 午夜福利在线观看吧| 亚洲av成人av| 亚洲av第一区精品v没综合| 性插视频无遮挡在线免费观看| 91久久精品国产一区二区成人| 亚洲欧美清纯卡通| 亚洲欧美日韩高清专用| 日本成人三级电影网站| 亚洲在线观看片| 久久久久网色| 久久综合国产亚洲精品| 一边摸一边抽搐一进一小说| 久久久国产成人免费| 午夜免费激情av| 欧美成人a在线观看| 嫩草影院入口| 国产一区二区亚洲精品在线观看| 国产男人的电影天堂91| 国产精品一及| 亚洲av.av天堂| 村上凉子中文字幕在线| 亚洲精品久久久久久婷婷小说 | 国产午夜精品久久久久久一区二区三区| 亚洲激情五月婷婷啪啪| 哪个播放器可以免费观看大片| 日韩欧美国产在线观看| 国产麻豆成人av免费视频| 级片在线观看| 国产伦精品一区二区三区四那| 九九在线视频观看精品| 国产精品一区二区性色av| 国产精品人妻久久久影院| 美女内射精品一级片tv| 女人被狂操c到高潮| 久久久久久久久中文| 成人毛片60女人毛片免费| 大型黄色视频在线免费观看| 久久精品91蜜桃| 日韩国内少妇激情av| 国产亚洲av片在线观看秒播厂 | 能在线免费看毛片的网站| 一级二级三级毛片免费看| 国产成人福利小说| 狂野欧美激情性xxxx在线观看| 欧美性猛交黑人性爽| 国产视频首页在线观看| 偷拍熟女少妇极品色| 国产一区二区三区在线臀色熟女| 久久久久久大精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲真实伦在线观看| 亚洲国产高清在线一区二区三| 亚洲中文字幕一区二区三区有码在线看| 日韩一本色道免费dvd| h日本视频在线播放| 日韩国内少妇激情av| 内射极品少妇av片p| 国产精品麻豆人妻色哟哟久久 | 日韩一区二区视频免费看| 男女那种视频在线观看| 美女被艹到高潮喷水动态| 国产欧美日韩精品一区二区| 国产日韩欧美在线精品| 精品久久久久久久久久免费视频| 久久亚洲国产成人精品v| 亚洲在久久综合| 成年版毛片免费区| 亚洲激情五月婷婷啪啪| 欧美+亚洲+日韩+国产| 91aial.com中文字幕在线观看| 国产精品爽爽va在线观看网站| 中文字幕精品亚洲无线码一区| 久久精品国产亚洲网站| 麻豆成人av视频| 村上凉子中文字幕在线| 国产高清有码在线观看视频| 中文在线观看免费www的网站| av福利片在线观看| 国产精品爽爽va在线观看网站| 在线免费观看不下载黄p国产| 天堂av国产一区二区熟女人妻| 2021天堂中文幕一二区在线观| 中国国产av一级| 国产精品一二三区在线看| 少妇高潮的动态图| 国产精品国产三级国产av玫瑰| 亚洲欧美清纯卡通| 波多野结衣巨乳人妻| 悠悠久久av| 国产成人福利小说| 欧美3d第一页| h日本视频在线播放| 亚洲欧洲国产日韩| 亚洲国产精品成人久久小说 | avwww免费| 岛国在线免费视频观看| 国产av一区在线观看免费| 成年版毛片免费区| 久久99蜜桃精品久久| 99riav亚洲国产免费| 看黄色毛片网站| 国产伦精品一区二区三区视频9| 好男人视频免费观看在线| 联通29元200g的流量卡| 只有这里有精品99| 亚洲精品国产成人久久av| 亚洲av电影不卡..在线观看| www.av在线官网国产| 亚洲欧洲国产日韩| 久久99热6这里只有精品| 波多野结衣巨乳人妻| 久久久久久九九精品二区国产| 美女xxoo啪啪120秒动态图| 桃色一区二区三区在线观看| 国产精品1区2区在线观看.| 国产成人福利小说| 欧美成人免费av一区二区三区| 白带黄色成豆腐渣| 丰满的人妻完整版| 秋霞在线观看毛片| 人人妻人人澡人人爽人人夜夜 | 亚洲精品乱码久久久v下载方式| 婷婷精品国产亚洲av| 只有这里有精品99| 一本久久中文字幕| 国产精品不卡视频一区二区| 高清日韩中文字幕在线| 少妇人妻精品综合一区二区 | 国产精品国产三级国产av玫瑰| 久久久色成人| 国内精品一区二区在线观看| 亚洲自偷自拍三级| 桃色一区二区三区在线观看| 寂寞人妻少妇视频99o| АⅤ资源中文在线天堂| 午夜福利成人在线免费观看| 天堂中文最新版在线下载 | 99热6这里只有精品| 中文字幕av在线有码专区| 中文字幕av成人在线电影| 少妇猛男粗大的猛烈进出视频 | 国产真实伦视频高清在线观看| 少妇裸体淫交视频免费看高清| 晚上一个人看的免费电影| 男的添女的下面高潮视频| 国产精品女同一区二区软件| 日韩一区二区视频免费看| 欧美日韩在线观看h| 少妇猛男粗大的猛烈进出视频 | 神马国产精品三级电影在线观看| 日韩av在线大香蕉| 国产av麻豆久久久久久久| 亚洲乱码一区二区免费版| 国产黄色小视频在线观看| 人妻久久中文字幕网| 国产在线精品亚洲第一网站| 熟妇人妻久久中文字幕3abv| 久久久久国产网址| 99久久九九国产精品国产免费| 欧美又色又爽又黄视频| 日日撸夜夜添| 日日撸夜夜添| 嫩草影院新地址| 一本久久中文字幕| 成人亚洲精品av一区二区| 18+在线观看网站| 夜夜看夜夜爽夜夜摸| 3wmmmm亚洲av在线观看| 美女内射精品一级片tv| 精品日产1卡2卡| 在线播放无遮挡| 久久精品91蜜桃| 亚洲精品日韩av片在线观看| 床上黄色一级片| 一边摸一边抽搐一进一小说| 亚洲一区高清亚洲精品| 成人美女网站在线观看视频| av在线蜜桃| 亚洲欧美成人综合另类久久久 | 日日撸夜夜添| 精品免费久久久久久久清纯| 久久人人爽人人片av| 国产探花在线观看一区二区| 一区二区三区高清视频在线| 男人的好看免费观看在线视频| 国内精品一区二区在线观看| 国产黄片视频在线免费观看| 午夜免费男女啪啪视频观看| 91av网一区二区| 18禁裸乳无遮挡免费网站照片| 小蜜桃在线观看免费完整版高清| 亚州av有码| 久久人人爽人人片av| 国产午夜福利久久久久久| 国产精品精品国产色婷婷| 一本精品99久久精品77| 黄片wwwwww| 久久精品夜色国产| 国产亚洲精品久久久久久毛片| 一个人免费在线观看电影| 观看免费一级毛片| 日本与韩国留学比较| av卡一久久| 99精品在免费线老司机午夜| 91麻豆精品激情在线观看国产| 中文字幕精品亚洲无线码一区| 看十八女毛片水多多多| 欧美日韩乱码在线| 成人综合一区亚洲| 26uuu在线亚洲综合色| 欧美3d第一页| 成人毛片a级毛片在线播放| 亚洲va在线va天堂va国产| 免费av观看视频| 国产黄a三级三级三级人| 国产精品蜜桃在线观看 | 国产成人91sexporn| 日韩 亚洲 欧美在线| 国产精品.久久久| 一级毛片久久久久久久久女| 国产综合懂色| 国产极品天堂在线| 国产片特级美女逼逼视频| 亚洲欧美日韩东京热| 青春草亚洲视频在线观看| 亚洲精品乱码久久久v下载方式| 欧美+亚洲+日韩+国产| 麻豆成人午夜福利视频| 欧美区成人在线视频| 亚洲aⅴ乱码一区二区在线播放| 美女被艹到高潮喷水动态| 久久久久久九九精品二区国产| 毛片女人毛片| 欧美最新免费一区二区三区| 午夜福利在线观看吧| 婷婷精品国产亚洲av| 嫩草影院精品99| 一级毛片aaaaaa免费看小| 成人国产麻豆网| 精品无人区乱码1区二区| ponron亚洲| 国产精品久久久久久av不卡| 天天一区二区日本电影三级| 中文字幕熟女人妻在线| 国产 一区 欧美 日韩| 久久亚洲精品不卡| 99热这里只有是精品在线观看| 久久精品人妻少妇| 久久国内精品自在自线图片| 国产精品人妻久久久久久| 久久久成人免费电影| 在线免费观看不下载黄p国产| ponron亚洲| 中国美女看黄片| 亚洲欧美成人精品一区二区| 精品久久久久久久久av| 国产亚洲av片在线观看秒播厂 | 99热全是精品| 黄片无遮挡物在线观看| 国产片特级美女逼逼视频| 国产成人精品久久久久久| 日日摸夜夜添夜夜爱| 国产色爽女视频免费观看| 村上凉子中文字幕在线| 三级毛片av免费| 日韩一本色道免费dvd| 日韩欧美精品免费久久| 精品少妇黑人巨大在线播放 | 狂野欧美激情性xxxx在线观看| 秋霞在线观看毛片| 此物有八面人人有两片| 最近手机中文字幕大全| 国产成人精品久久久久久| 搡老妇女老女人老熟妇| 日韩成人伦理影院| 久久久久久国产a免费观看| 成人欧美大片| 亚洲av免费在线观看| 亚洲人成网站高清观看| 搡女人真爽免费视频火全软件| 国产精品久久电影中文字幕| 91午夜精品亚洲一区二区三区| 国产精华一区二区三区| 日韩,欧美,国产一区二区三区 | 免费人成视频x8x8入口观看| 天堂影院成人在线观看| 日日撸夜夜添| 99热这里只有精品一区| 在线观看一区二区三区| 美女高潮的动态| 亚洲精品成人久久久久久| 中文字幕精品亚洲无线码一区| 国产高清不卡午夜福利| 免费人成视频x8x8入口观看| 国产成人freesex在线| 一区二区三区四区激情视频 | 亚洲欧美日韩卡通动漫| 亚洲国产日韩欧美精品在线观看| 欧美日韩乱码在线| 国产成人一区二区在线| 国产男人的电影天堂91| 99热精品在线国产| 久久久久久久久久黄片| 99热精品在线国产| 美女黄网站色视频| 精品国内亚洲2022精品成人| 爱豆传媒免费全集在线观看| 午夜视频国产福利| 亚洲av一区综合| 免费av观看视频| av又黄又爽大尺度在线免费看 | 欧美色视频一区免费| 草草在线视频免费看| 精品久久久久久成人av| 精华霜和精华液先用哪个| 成熟少妇高潮喷水视频| 欧美日韩在线观看h| 亚洲欧美精品综合久久99| 久久欧美精品欧美久久欧美| 最近2019中文字幕mv第一页| 久久这里有精品视频免费| 欧美精品一区二区大全| 久久精品影院6| 校园春色视频在线观看| 亚洲人成网站在线播放欧美日韩| 国产一区二区激情短视频| 尤物成人国产欧美一区二区三区| 夜夜看夜夜爽夜夜摸| a级毛片免费高清观看在线播放| 免费看日本二区| 精品一区二区三区人妻视频| 国产一区二区激情短视频| 国产中年淑女户外野战色| 亚洲在久久综合| 国产精品日韩av在线免费观看| 久久久国产成人精品二区| 国产老妇女一区| 最近中文字幕高清免费大全6| 国产又黄又爽又无遮挡在线| 日本与韩国留学比较| 亚洲精品久久国产高清桃花| 夫妻性生交免费视频一级片| 性插视频无遮挡在线免费观看| 在线观看免费视频日本深夜| 欧美在线一区亚洲| 91精品一卡2卡3卡4卡| 97在线视频观看| 国产精品电影一区二区三区| 欧美一区二区亚洲| 青春草视频在线免费观看| 在线观看一区二区三区| 蜜桃久久精品国产亚洲av| 国产 一区精品| 欧美高清性xxxxhd video| 欧美bdsm另类| 天堂影院成人在线观看| 久久久久久伊人网av| 国产单亲对白刺激| 久久人人精品亚洲av| 午夜爱爱视频在线播放| 日产精品乱码卡一卡2卡三| 少妇熟女欧美另类| 亚洲精品成人久久久久久| 国产v大片淫在线免费观看| 热99在线观看视频| 九九爱精品视频在线观看| 亚州av有码| 午夜久久久久精精品| АⅤ资源中文在线天堂| 亚洲丝袜综合中文字幕| 我要搜黄色片| 欧美激情久久久久久爽电影| а√天堂www在线а√下载| 国产精品一区二区三区四区免费观看| 亚洲欧美日韩卡通动漫| 你懂的网址亚洲精品在线观看 | 久久99热6这里只有精品| 男人舔女人下体高潮全视频| 欧美成人免费av一区二区三区| av免费在线看不卡| 中文在线观看免费www的网站| 国产伦精品一区二区三区视频9| 精品久久久噜噜| 深夜精品福利| 日韩成人伦理影院| 一级毛片久久久久久久久女| 国产精品1区2区在线观看.| 好男人视频免费观看在线| 麻豆久久精品国产亚洲av| 欧美色视频一区免费| 欧美+亚洲+日韩+国产| 久久久久性生活片| 亚洲真实伦在线观看| 亚洲国产精品国产精品| 亚洲无线在线观看| 亚洲最大成人中文| 日韩高清综合在线| 99久久九九国产精品国产免费| 给我免费播放毛片高清在线观看| 寂寞人妻少妇视频99o| 免费看av在线观看网站| 精品久久久久久久久亚洲| 国产亚洲av片在线观看秒播厂 | 国产成人午夜福利电影在线观看| 麻豆国产97在线/欧美| 欧美人与善性xxx| а√天堂www在线а√下载| 老熟妇乱子伦视频在线观看| 欧美zozozo另类| 国内精品美女久久久久久| 一级黄片播放器| 91av网一区二区| 亚洲av电影不卡..在线观看| 精华霜和精华液先用哪个| 最近最新中文字幕大全电影3| 国产精品一及| 国产精品人妻久久久久久| 超碰av人人做人人爽久久| 美女脱内裤让男人舔精品视频 | 午夜精品国产一区二区电影 | 99在线视频只有这里精品首页| 欧美三级亚洲精品| a级毛色黄片| 亚洲精品成人久久久久久| 亚洲在线自拍视频| 国产精品日韩av在线免费观看| 黄色欧美视频在线观看| 欧美性猛交╳xxx乱大交人| 91狼人影院| 性色avwww在线观看| 精品熟女少妇av免费看| 日本黄色片子视频| 久久精品久久久久久噜噜老黄 | 中文字幕av在线有码专区| 亚洲欧洲日产国产| 国语自产精品视频在线第100页| 能在线免费观看的黄片| 精品人妻视频免费看| 男女那种视频在线观看| 国产一区二区激情短视频| 国产精品不卡视频一区二区| 观看美女的网站| 永久网站在线| 国产人妻一区二区三区在| 最新中文字幕久久久久| 亚洲国产精品国产精品| 舔av片在线| 欧美日本视频| 国模一区二区三区四区视频| 熟女电影av网| 国产黄色视频一区二区在线观看 | 亚洲欧美成人精品一区二区| 少妇裸体淫交视频免费看高清| 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产自在天天线| 日本黄色片子视频| 午夜福利高清视频| 久久国内精品自在自线图片| 久久99蜜桃精品久久| 啦啦啦啦在线视频资源| 亚洲av中文av极速乱| 日本黄色片子视频| 午夜a级毛片| 成人美女网站在线观看视频| 亚洲无线在线观看| 天天一区二区日本电影三级| 亚洲婷婷狠狠爱综合网| 国产免费一级a男人的天堂| 一级毛片我不卡| 在线观看免费视频日本深夜| 人妻系列 视频| 欧美bdsm另类| 亚洲av中文字字幕乱码综合| 免费看美女性在线毛片视频| 欧美激情久久久久久爽电影| 天天躁日日操中文字幕| 国产乱人视频| 99久国产av精品| 亚洲精品粉嫩美女一区| 国产伦在线观看视频一区| 一边摸一边抽搐一进一小说| 三级国产精品欧美在线观看| 久久精品夜色国产| www.色视频.com| 联通29元200g的流量卡| 国产成人a∨麻豆精品| .国产精品久久| 久久99热6这里只有精品| 国产一区二区在线观看日韩| 美女被艹到高潮喷水动态| 亚洲av一区综合| 婷婷精品国产亚洲av| 99在线视频只有这里精品首页| 少妇人妻一区二区三区视频|