• <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
    淫妇啪啪啪对白视频| 最近最新中文字幕大全免费视频| 在线十欧美十亚洲十日本专区| 国产精华一区二区三区| 国产三级在线视频| 久久精品亚洲精品国产色婷小说| 桃红色精品国产亚洲av| 亚洲av免费在线观看| 一级毛片精品| 欧美乱妇无乱码| 大型黄色视频在线免费观看| 国产单亲对白刺激| 九九久久精品国产亚洲av麻豆 | 日本在线视频免费播放| 国产精品女同一区二区软件 | h日本视频在线播放| 在线观看日韩欧美| 亚洲乱码一区二区免费版| 脱女人内裤的视频| 国产亚洲精品久久久com| 天堂动漫精品| 成人av一区二区三区在线看| 欧美另类亚洲清纯唯美| 亚洲精品美女久久久久99蜜臀| 成人精品一区二区免费| 欧美乱妇无乱码| 国产午夜精品久久久久久| 1000部很黄的大片| 99久久99久久久精品蜜桃| 97超级碰碰碰精品色视频在线观看| 又大又爽又粗| 美女黄网站色视频| 国产亚洲av高清不卡| 久久人人精品亚洲av| 麻豆成人av在线观看| 欧美丝袜亚洲另类 | 99久久久亚洲精品蜜臀av| 国产欧美日韩精品一区二区| 免费在线观看成人毛片| 搡老妇女老女人老熟妇| 国产不卡一卡二| 成年女人看的毛片在线观看| 国产高清videossex| a级毛片在线看网站| 亚洲精品美女久久av网站| 国内少妇人妻偷人精品xxx网站 | 国产精品一及| 亚洲专区中文字幕在线| 国产精品98久久久久久宅男小说| 免费在线观看日本一区| 亚洲成a人片在线一区二区| 99久久无色码亚洲精品果冻| 99精品欧美一区二区三区四区| 精品国产超薄肉色丝袜足j| 国产人伦9x9x在线观看| 成人性生交大片免费视频hd| 亚洲欧美日韩卡通动漫| 亚洲国产欧洲综合997久久,| 村上凉子中文字幕在线| 日本成人三级电影网站| 黑人欧美特级aaaaaa片| 午夜福利成人在线免费观看| 操出白浆在线播放| 一个人观看的视频www高清免费观看 | 亚洲av片天天在线观看| 老司机午夜福利在线观看视频| 高清毛片免费观看视频网站| 久久精品夜夜夜夜夜久久蜜豆| 亚洲国产高清在线一区二区三| 精品国产乱子伦一区二区三区| 在线免费观看不下载黄p国产 | 久久久国产精品麻豆| 91字幕亚洲| 在线观看免费视频日本深夜| 国产欧美日韩精品一区二区| 久久人妻av系列| 亚洲天堂国产精品一区在线| 中出人妻视频一区二区| 三级毛片av免费| 精品久久久久久久末码| 久久久久久久午夜电影| 久久99热这里只有精品18| 午夜福利在线在线| 亚洲熟女毛片儿| 久久精品亚洲精品国产色婷小说| 国产精品一区二区免费欧美| 欧美另类亚洲清纯唯美| 亚洲七黄色美女视频| 亚洲精品美女久久久久99蜜臀| 精品人妻1区二区| 国产一级毛片七仙女欲春2| 午夜两性在线视频| 网址你懂的国产日韩在线| 免费看光身美女| 天堂av国产一区二区熟女人妻| 国产精品一区二区免费欧美| 在线播放国产精品三级| 精品无人区乱码1区二区| 可以在线观看的亚洲视频| 男人和女人高潮做爰伦理| a级毛片a级免费在线| 一个人免费在线观看的高清视频| 九色成人免费人妻av| 国产欧美日韩一区二区三| 国产高清激情床上av| 男女下面进入的视频免费午夜| 国产野战对白在线观看| 男女午夜视频在线观看| 一区二区三区激情视频| 又黄又爽又免费观看的视频| 天天一区二区日本电影三级| 一本精品99久久精品77| 欧美大码av| 伦理电影免费视频| 级片在线观看| 国产极品精品免费视频能看的| 亚洲精品456在线播放app | 国产精品自产拍在线观看55亚洲| 国产欧美日韩一区二区三| 2021天堂中文幕一二区在线观| 1000部很黄的大片| 亚洲成人免费电影在线观看| 中文字幕av在线有码专区| 啦啦啦观看免费观看视频高清| 一本精品99久久精品77| 蜜桃久久精品国产亚洲av| e午夜精品久久久久久久| 色综合欧美亚洲国产小说| 久久久久国产一级毛片高清牌| 可以在线观看的亚洲视频| 热99re8久久精品国产| 美女 人体艺术 gogo| 久99久视频精品免费| bbb黄色大片| 欧美日韩一级在线毛片| 亚洲人与动物交配视频| 狂野欧美白嫩少妇大欣赏| 99久久国产精品久久久| 99热6这里只有精品| 久9热在线精品视频| 久久亚洲真实| 久久亚洲真实| 老鸭窝网址在线观看| 91在线精品国自产拍蜜月 | 亚洲激情在线av| 国产精品香港三级国产av潘金莲| 日韩成人在线观看一区二区三区| 国产精品久久久久久亚洲av鲁大| 亚洲美女黄片视频| 日韩欧美一区二区三区在线观看| 亚洲美女视频黄频| 久久中文看片网| 久久久久久久久久黄片| 国产精品免费一区二区三区在线| 久久久成人免费电影| 美女高潮喷水抽搐中文字幕| 宅男免费午夜| 黄色女人牲交| 91麻豆精品激情在线观看国产| 香蕉av资源在线| 亚洲熟妇熟女久久| av在线天堂中文字幕| 听说在线观看完整版免费高清| 深夜精品福利| www.999成人在线观看| 久久精品综合一区二区三区| 偷拍熟女少妇极品色| 精品久久久久久久末码| 丝袜人妻中文字幕| 国产欧美日韩一区二区精品| 日韩有码中文字幕| 很黄的视频免费| 亚洲avbb在线观看| 99久久综合精品五月天人人| 欧美黄色淫秽网站| 国产久久久一区二区三区| 国产综合懂色| 黑人欧美特级aaaaaa片| 美女被艹到高潮喷水动态| 国产激情欧美一区二区| 国产高清三级在线| 最近最新中文字幕大全电影3| 伦理电影免费视频| 精品久久久久久久久久免费视频| 啦啦啦观看免费观看视频高清| 真人做人爱边吃奶动态| 亚洲无线观看免费| 少妇熟女aⅴ在线视频| 麻豆成人av在线观看| 色哟哟哟哟哟哟| 日韩精品青青久久久久久| 美女午夜性视频免费| 欧美午夜高清在线| 久久久久久久午夜电影| 国产精品自产拍在线观看55亚洲| 国产精品久久视频播放| 成年人黄色毛片网站| 欧美成人免费av一区二区三区| 国产又黄又爽又无遮挡在线| 99精品在免费线老司机午夜| 老司机在亚洲福利影院| 国产极品精品免费视频能看的| 欧美绝顶高潮抽搐喷水| 欧美日韩国产亚洲二区| 高潮久久久久久久久久久不卡| 精品一区二区三区av网在线观看| 日日夜夜操网爽| 天天一区二区日本电影三级| 三级国产精品欧美在线观看 | 一个人观看的视频www高清免费观看 | 午夜福利免费观看在线| 热99re8久久精品国产| 国产精品久久久久久人妻精品电影| 亚洲五月天丁香| 操出白浆在线播放| 亚洲一区高清亚洲精品| 精品久久久久久久久久免费视频| 美女被艹到高潮喷水动态| 天天一区二区日本电影三级| 亚洲 欧美 日韩 在线 免费| 99国产极品粉嫩在线观看| 久久亚洲精品不卡| 18禁国产床啪视频网站| 亚洲av电影在线进入| 免费搜索国产男女视频| 日本撒尿小便嘘嘘汇集6| 国产人伦9x9x在线观看| 欧美最黄视频在线播放免费| 国产亚洲精品一区二区www| 免费av不卡在线播放| 午夜免费成人在线视频| 精品99又大又爽又粗少妇毛片 | 免费av毛片视频| 90打野战视频偷拍视频| 国产精品av久久久久免费| 免费av毛片视频| 日韩免费av在线播放| 丁香欧美五月| 99热只有精品国产| 99在线人妻在线中文字幕| 日本免费a在线| 久久久久国内视频| xxxwww97欧美| 一进一出抽搐动态| 精品99又大又爽又粗少妇毛片 | 国产又黄又爽又无遮挡在线| 成人一区二区视频在线观看| 国语自产精品视频在线第100页| 色哟哟哟哟哟哟| av国产免费在线观看| 久久婷婷人人爽人人干人人爱| 欧美一区二区国产精品久久精品| 国产高清有码在线观看视频| 麻豆国产av国片精品| 狂野欧美激情性xxxx| 国产高清三级在线| 国产成+人综合+亚洲专区| 又黄又爽又免费观看的视频| 亚洲精品在线美女| 午夜免费成人在线视频| 午夜免费观看网址| 噜噜噜噜噜久久久久久91| 国产精品久久久久久人妻精品电影| 国产精华一区二区三区| 日韩欧美在线二视频| 韩国av一区二区三区四区| 中文亚洲av片在线观看爽| 国产亚洲av高清不卡| 婷婷精品国产亚洲av| 国产一级毛片七仙女欲春2| 亚洲国产精品999在线| 哪里可以看免费的av片| 亚洲精品久久国产高清桃花| 最新在线观看一区二区三区| 中文字幕久久专区| 麻豆久久精品国产亚洲av| 美女免费视频网站| 国产一区二区三区在线臀色熟女| 黄片大片在线免费观看| 又爽又黄无遮挡网站| 国产激情偷乱视频一区二区| 可以在线观看毛片的网站| 亚洲欧美精品综合一区二区三区| 国产精品影院久久| 黄片小视频在线播放| av在线天堂中文字幕| 国产真实乱freesex| 国产高清videossex| 欧美乱色亚洲激情| 人妻夜夜爽99麻豆av| 在线免费观看的www视频| 欧美成人性av电影在线观看| 操出白浆在线播放| 亚洲人成网站高清观看| 国产精品国产高清国产av| av中文乱码字幕在线| 亚洲无线在线观看| 操出白浆在线播放| 亚洲在线自拍视频| 黄色日韩在线| 少妇丰满av| 性色av乱码一区二区三区2| 亚洲欧美激情综合另类| 丁香六月欧美| 一进一出抽搐gif免费好疼| 国产高清激情床上av| 亚洲精品中文字幕一二三四区| 亚洲国产色片| 日韩欧美国产在线观看| 亚洲欧美精品综合一区二区三区| 日韩 欧美 亚洲 中文字幕| 久久精品综合一区二区三区| 国产精品美女特级片免费视频播放器 | 国产久久久一区二区三区| 亚洲专区字幕在线| 国产爱豆传媒在线观看| 色综合站精品国产| 99热6这里只有精品| 啦啦啦观看免费观看视频高清| 国内精品美女久久久久久| 97超视频在线观看视频| 最近最新中文字幕大全免费视频| 女人被狂操c到高潮| 成人欧美大片| 亚洲中文字幕日韩| avwww免费| 黄频高清免费视频| e午夜精品久久久久久久| 黑人操中国人逼视频| 欧美+亚洲+日韩+国产| 一二三四社区在线视频社区8| 免费人成视频x8x8入口观看| 欧美成人免费av一区二区三区| 好看av亚洲va欧美ⅴa在| 91字幕亚洲| 国产精品九九99| 国产成人aa在线观看| 亚洲第一电影网av| 网址你懂的国产日韩在线| 久久亚洲精品不卡| 精品国产乱码久久久久久男人| 国产精品美女特级片免费视频播放器 | 99久久99久久久精品蜜桃| 男人舔奶头视频| 色综合站精品国产| 非洲黑人性xxxx精品又粗又长| 亚洲av中文字字幕乱码综合| 夜夜看夜夜爽夜夜摸| 久久精品国产清高在天天线| 岛国在线免费视频观看| 99久国产av精品| 两性午夜刺激爽爽歪歪视频在线观看| 91av网一区二区| 99国产精品一区二区蜜桃av| 国产激情久久老熟女| 亚洲色图av天堂| 1024香蕉在线观看| 久久久久九九精品影院| 99热这里只有是精品50| 我的老师免费观看完整版| 老汉色av国产亚洲站长工具| 麻豆成人午夜福利视频| 免费在线观看日本一区| 成年女人看的毛片在线观看| 99久久精品热视频| 久久久久九九精品影院| 黄色丝袜av网址大全| 在线观看免费视频日本深夜| 亚洲国产欧洲综合997久久,| 别揉我奶头~嗯~啊~动态视频| 午夜精品一区二区三区免费看| 久久久久久九九精品二区国产| 国产一区二区激情短视频| 99热6这里只有精品| 曰老女人黄片| 亚洲激情在线av| 欧美日韩福利视频一区二区| 韩国av一区二区三区四区| 91av网一区二区| 一个人观看的视频www高清免费观看 | 女警被强在线播放| 香蕉久久夜色| 午夜福利欧美成人| 免费看十八禁软件| 国产aⅴ精品一区二区三区波| 精品国产超薄肉色丝袜足j| 黄色女人牲交| 精品一区二区三区视频在线 | 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品色激情综合| 日韩精品青青久久久久久| 97超视频在线观看视频| 麻豆国产av国片精品| 中文资源天堂在线| avwww免费| 午夜免费激情av| 真人做人爱边吃奶动态| 日韩欧美精品v在线| 嫩草影院入口| 中文字幕久久专区| av天堂在线播放| 色播亚洲综合网| 两个人看的免费小视频| 久久精品人妻少妇| 我要搜黄色片| 国产精品永久免费网站| 身体一侧抽搐| 淫秽高清视频在线观看| 色av中文字幕| 成年女人看的毛片在线观看| a级毛片a级免费在线| 熟女少妇亚洲综合色aaa.| 亚洲中文字幕一区二区三区有码在线看 | 国产精品九九99| 成人鲁丝片一二三区免费| 免费av不卡在线播放| 麻豆国产97在线/欧美| 首页视频小说图片口味搜索| 免费无遮挡裸体视频| 男女床上黄色一级片免费看| 99热6这里只有精品| 激情在线观看视频在线高清| 亚洲av免费在线观看| 国产成人av教育| 亚洲va日本ⅴa欧美va伊人久久| 免费av毛片视频| 日本一二三区视频观看| 精品99又大又爽又粗少妇毛片 | 国产av一区在线观看免费| 我的老师免费观看完整版| 搡老妇女老女人老熟妇| 久久久精品大字幕| 亚洲国产精品sss在线观看| 嫁个100分男人电影在线观看| netflix在线观看网站| 国内揄拍国产精品人妻在线| 欧美色视频一区免费| 国产午夜精品论理片| 国产精品久久久人人做人人爽| 久99久视频精品免费| 两性午夜刺激爽爽歪歪视频在线观看| 香蕉av资源在线| 国产精品精品国产色婷婷| 老司机午夜十八禁免费视频| 欧美日韩精品网址| 97超视频在线观看视频| 99久久国产精品久久久| 少妇熟女aⅴ在线视频| 国产精品免费一区二区三区在线| 成人av在线播放网站| 超碰成人久久| 一进一出好大好爽视频| 亚洲欧美激情综合另类| 少妇熟女aⅴ在线视频| 在线免费观看的www视频| 免费在线观看成人毛片| 国产视频一区二区在线看| 国产精品av久久久久免费| 久久人妻av系列| 精品久久久久久,| 免费观看的影片在线观看| 91九色精品人成在线观看| 国产精品免费一区二区三区在线| 国产日本99.免费观看| 日本成人三级电影网站| 国产精品一区二区三区四区久久| 两人在一起打扑克的视频| 久久香蕉精品热| 免费大片18禁| 精品国产乱子伦一区二区三区| 亚洲成人中文字幕在线播放| 国产欧美日韩一区二区精品| 欧美日韩瑟瑟在线播放| 久久精品影院6| 精品久久久久久成人av| 国产精品综合久久久久久久免费| 男女之事视频高清在线观看| 黑人欧美特级aaaaaa片| 女人高潮潮喷娇喘18禁视频| 99国产精品一区二区蜜桃av| 免费在线观看日本一区| av欧美777| 久久这里只有精品19| 成人特级黄色片久久久久久久| 少妇裸体淫交视频免费看高清| 九色国产91popny在线| 亚洲国产精品sss在线观看| 午夜福利视频1000在线观看| 亚洲一区二区三区不卡视频| 99精品欧美一区二区三区四区| 婷婷丁香在线五月| 欧美+亚洲+日韩+国产| 亚洲精品久久国产高清桃花| 国产精品99久久久久久久久| 天堂√8在线中文| 欧美极品一区二区三区四区| 日本免费一区二区三区高清不卡| 日韩av在线大香蕉| 可以在线观看毛片的网站| 日本 av在线| 色视频www国产| 国内少妇人妻偷人精品xxx网站 | 免费观看精品视频网站| 男人舔女人下体高潮全视频| 美女 人体艺术 gogo| 亚洲熟女毛片儿| 99久久无色码亚洲精品果冻| 好看av亚洲va欧美ⅴa在| 国产黄a三级三级三级人| 久久久久精品国产欧美久久久| 亚洲欧美日韩无卡精品| 国产成人aa在线观看| 夜夜看夜夜爽夜夜摸| 亚洲国产精品成人综合色| netflix在线观看网站| 老汉色av国产亚洲站长工具| 国产av一区在线观看免费| 无遮挡黄片免费观看| 国产美女午夜福利| 国产精品久久久久久久电影 | av在线天堂中文字幕| 国产欧美日韩一区二区三| aaaaa片日本免费| 12—13女人毛片做爰片一| 国产成人精品久久二区二区91| 国产成人精品久久二区二区免费| 免费电影在线观看免费观看| 亚洲国产中文字幕在线视频| 日韩欧美精品v在线| 国产美女午夜福利| 国产99白浆流出| 18禁裸乳无遮挡免费网站照片| 搞女人的毛片| 久久久久久久午夜电影| 亚洲在线自拍视频| 丰满的人妻完整版| 校园春色视频在线观看| 午夜影院日韩av| 久久这里只有精品19| 久久久久久九九精品二区国产| 十八禁人妻一区二区| 久久国产精品影院| 国产成人精品久久二区二区91| 国产亚洲精品综合一区在线观看| 亚洲精品美女久久久久99蜜臀| 天堂影院成人在线观看| 最新中文字幕久久久久 | 啦啦啦韩国在线观看视频| 国产精品,欧美在线| 欧美极品一区二区三区四区| 天天躁狠狠躁夜夜躁狠狠躁| 狂野欧美白嫩少妇大欣赏| 成年女人看的毛片在线观看| a级毛片在线看网站| 成年版毛片免费区| 母亲3免费完整高清在线观看| 熟妇人妻久久中文字幕3abv| 亚洲狠狠婷婷综合久久图片| 国产1区2区3区精品| 久久精品影院6| 国产在线精品亚洲第一网站| 亚洲激情在线av| 日本在线视频免费播放| 美女午夜性视频免费| 久久久久国产一级毛片高清牌| 免费看a级黄色片| 久久这里只有精品19| 此物有八面人人有两片| 在线观看日韩欧美| 国产一区二区三区视频了| 欧美xxxx黑人xx丫x性爽| 亚洲欧美精品综合一区二区三区| 国产成人精品久久二区二区91| 在线十欧美十亚洲十日本专区| 国产黄色小视频在线观看| 亚洲专区字幕在线| 在线观看美女被高潮喷水网站 | 男人舔奶头视频| 亚洲人与动物交配视频| 国内少妇人妻偷人精品xxx网站 | 99视频精品全部免费 在线 | 免费无遮挡裸体视频| 亚洲国产精品成人综合色| h日本视频在线播放| 日本 欧美在线| 在线观看66精品国产| a在线观看视频网站| 动漫黄色视频在线观看| 真人一进一出gif抽搐免费| av在线蜜桃| 在线观看日韩欧美| 91字幕亚洲| 亚洲人成伊人成综合网2020| 最好的美女福利视频网| xxxwww97欧美| 亚洲美女黄片视频| 亚洲精品456在线播放app | 一级a爱片免费观看的视频| 国产成人欧美在线观看| 一级a爱片免费观看的视频| 一卡2卡三卡四卡精品乱码亚洲| 最新中文字幕久久久久 | 欧美绝顶高潮抽搐喷水| 欧美日韩乱码在线| 免费看光身美女| 免费搜索国产男女视频| 婷婷精品国产亚洲av| 亚洲精品久久国产高清桃花| 亚洲人成网站在线播放欧美日韩| 丝袜人妻中文字幕| 一区二区三区激情视频| 后天国语完整版免费观看| 欧美一区二区精品小视频在线| bbb黄色大片|