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

    航空發(fā)動機燃燒室光學可視模型試驗件及其流場測量研究進展

    2021-03-20 03:37:10王于藍范雄杰劉存喜楊金虎劉富強
    實驗流體力學 2021年1期
    關鍵詞:視窗旋流器旋流

    王于藍, 范雄杰, 高 偉, 劉存喜,3,*, 楊金虎,3,劉富強,3, 穆 勇,3, 徐 綱,3,*

    1. 中國科學院工程熱物理研究所 中國科學院輕型動力重點實驗室, 北京 100190;2. 中國科學院大學, 北京 100049; 3. 中國科學院 輕型動力創(chuàng)新研究院, 北京 100190

    0 引 言

    航空發(fā)動機燃燒室將燃料的化學能轉化為熱能從而對渦輪做功[1]。航空發(fā)動機是飛機的“心臟”,而燃燒室可以說是“心臟”的“心臟”。隨著發(fā)動機燃燒室由傳統(tǒng)富油頭部設計向多級旋流燃油分級的貧油頭部設計發(fā)展[2],燃燒室頭部進氣量大幅增加,流場組織、燃油霧化、油-氣混合、燃燒組織等都在燃燒室頭部完成,燃燒室頭部的流場組織在燃油霧化、蒸發(fā)、混合和燃燒組織中的作用越來越重要,決定著燃燒室的點/熄火穩(wěn)定工作邊界、出口溫度分布質量、污染物排放水平和火焰筒壽命等燃燒性能。

    為滿足發(fā)動機推重比、壽命、污染物排放和安全性等性能要求,軍用發(fā)動機的高溫升燃燒室和民用發(fā)動機的低污染燃燒室普遍采用先進的分區(qū)分級耦合燃燒組織方法,如貧油預混預蒸發(fā)燃燒(LPP)、富燃-淬熄-貧燃燃燒(RQL)以及貧油直接噴射燃燒(LDI)等[3-5]。這些先進的燃燒組織思想已經部分成功運用于在役軍用和商用發(fā)動機燃燒室中,如美國通用電氣(GE)的雙環(huán)預混旋流器燃燒室(TAPS),英國羅爾斯·羅伊斯(RR)的Lean-burn燃燒室,普惠(PW)公司的TALON燃燒室,日本宇航局(JAXA)的Lean-stage燃燒室等。這些先進燃燒室采用多級旋流器的燃燒室頭部結構,流場組織過程和流場結構特征復雜,湍流強度高,且流場特性受旋流器氣動結構參數(shù)和燃燒室結構影響大。國內外均對燃燒室內的流場特性進行了大量的數(shù)值仿真和試驗研究[6-7],數(shù)值仿真方法雖然可以得到較完整的流場細節(jié),但由于燃燒室內流場的高湍流度,并且旋流流動與火焰筒受限空間及其形狀耦合,同時旋流流動過程受旋流器及其出口位置的微小結構影響較大,數(shù)值模擬方法尚不能得到完全真實的流場信息,尤其是在燃燒狀態(tài)下,對帶化學反應的湍流模擬更加困難。目前,發(fā)動機燃燒室內的流場組織機理及流場特性研究仍然以流場測量試驗為主。

    先進發(fā)動機或工業(yè)燃氣輪機燃燒室普遍采用環(huán)形燃燒室,具有機匣和火焰筒兩層結構,火焰筒上設計主燃孔、補燃孔、摻混孔等大孔結構和多種類型的氣膜或發(fā)散冷卻小孔等,幾何結構復雜,形狀各異,徑向尺寸達到米量級,幾何尺度范圍寬,無論是接觸式測點布置還是光學診斷的光路布置都相對困難[4]。為更好地捕捉燃燒室內流場特征,流場測試手段不斷發(fā)展,激光診斷技術由于不干擾流場,已被廣泛運用于燃燒室流場測量中,并隨著計算機技術、激光器技術、圖像采集及處理技術的快速發(fā)展,其流場測量精度逐漸提高。如單點測量的熱線風速儀(HWA)、激光多普勒測速儀(LDV)和平面測量的粒子圖像測速儀(PIV)等,國內外大多借助這些成熟的流場診斷技術對燃燒室流場開展研究,并且向高頻和三維體測量方向發(fā)展,同時也在積極探索新的測量方法[6]。美國密歇根大學、辛辛那提大學和普渡大學等與GE的聯(lián)合團隊、歐洲各個高校與RR公司、日本宇航局(JAXA)、德國宇航局(DLR)和法國宇航局(ONERA)等均建立了較為完備的燃燒室流場測量試驗平臺,借助光學模型燃燒室進行了大量的流場測量試驗。國內中國航空發(fā)動機集團公司下屬的動力所聯(lián)合北京航空航天大學、南京航空航天大學、西北工業(yè)大學等航空航天院校、中國科學院、中國空氣動力研究與發(fā)展中心等基于發(fā)動機燃燒室的研制需求,在光學模型燃燒室試驗件內開展了大量的流場測量試驗研究。這些光學模型燃燒室內的流場研究在發(fā)動機燃燒室研制過程中發(fā)揮了重要作用,已經引起了國內外研究院所和發(fā)動機公司的重視。

    光學模型燃燒室是研究發(fā)動機燃燒室流場組織機理和流場結構特征的關鍵載體,但由于目前在光學模型燃燒室試驗件設計方面缺少共識,設計什么樣的光學模型燃燒室試驗件困擾著發(fā)動機燃燒室燃燒組織基礎問題或應用基礎研究的科研工作者和燃燒室工程研制的工程師。即使上文提到的成熟流場測量技術應用于發(fā)動機燃燒室時,也受到燃燒室內的高溫高壓燃燒惡劣環(huán)境和復雜結構引起的光路布局等限制,影響精度和測試空間。本文重點總結光學模型燃燒室的發(fā)展趨勢,詳細分析在光學模型燃燒室試驗件設計過程中針對流場測量需要考慮的要求和存在的問題,進一步,分析兩類典型燃燒室的流場特征,從而把燃燒室本身的流場特征與模型試驗件設計結合起來,以期為發(fā)動機燃燒室研制過程中燃燒室光學模型試驗件的設計提供支撐。

    1 光學可視模型燃燒室

    光學可視模型燃燒室的作用是模擬燃燒室內的真實氣動、霧化和燃燒特征,并且滿足光學測量裝置的光路布置要求。為深入認識燃燒室內的流場組織機理和流場特性,國內外學者基于不同的研究目標,開展了大量的相關研究[1,4]。為了模擬發(fā)動機燃燒室內的真實流動組織過程和流場特征,其模型燃燒室試驗件的發(fā)展趨勢為從開放空間到受限空間、從常溫常壓到中溫中壓甚至高溫高壓、從冷態(tài)流場到燃燒條件下的熱態(tài)流場和從單頭部到多頭部甚至全環(huán)燃燒室等(見圖1),為發(fā)動機燃燒室研制過程中相關燃燒問題的解決和新型燃燒組織技術的發(fā)展提供了支撐,有效促進了發(fā)動機燃燒室的研制。

    根據(jù)圖1,光學可視模型燃燒室設計的最終目標是模擬真實燃燒室的真實工作條件并捕獲燃燒室內的氣動、霧化和燃燒等物理參數(shù)。受經濟因素、試驗條件、測量技術和測點及光路空間布局限制,在發(fā)動機燃燒室研制過程中,光學可視模型燃燒室的設計采用對測量參數(shù)數(shù)量及其空間、燃燒室結構和進氣參數(shù)等進行折中的方法。光學可視模型燃燒室試驗件向3個方向發(fā)展:1) 高溫高壓模型燃燒室試驗件。其更加接近燃燒室的真實工作參數(shù),忽略燃燒室頭部之間相互作用,并且對燃燒室結構進行簡化,只能采用少量能承受燃燒室內高溫高壓惡劣環(huán)境的測量技術,測量參數(shù)受限,試驗費用高,適用于技術成熟度達到TRL3的燃燒室技術研究。2) 含真實火焰筒并且工作在真實環(huán)境下的燃燒室試驗件。該試驗件充分考慮了進氣參數(shù)、燃燒室的結構和頭部之間相互作用,試驗件結構復雜,受測量方法在高溫高壓下的環(huán)境適應性和空間布局等限制,只能測量極少量參數(shù),試驗費用較高,適用于技術成熟度達到TRL5-TRL6的燃燒室技術研制。3) 全環(huán)燃燒室模型試驗件。該試驗件充分考慮了燃燒室頭部之間耦合相互作用,采用縮放或真實燃燒室頭部結構,簡化火焰筒結構,大多在常溫常壓條件下開展實驗,測量參數(shù)受環(huán)境適應性和空間布局等因素的限制較小,試驗成本相對較低,適用于技術成熟度達到1級的燃燒室技術研制和燃燒組織機理基礎研究。

    圖1 光學可視模型燃燒室試驗件發(fā)展趨勢Fig. 1 The development of optically accessible model combustor

    1.1 旋流杯光學模型燃燒室

    旋流杯光學模型燃燒室(以下簡稱“旋流杯模型燃燒室”)是美國GE公司發(fā)展的用于燃油霧化和油氣混合的燃燒組織技術,由于其良好的燃油霧化質量和燃燒室穩(wěn)定性,在發(fā)動機上得到廣泛應用,如CFM56[8]。國內外許多學者針對旋流杯模型燃燒室開展了大量的基礎及應用研究,辛辛那提大學在旋流杯模型燃燒室的流場組織機理研究和模型燃燒室試驗件設計方面都作出了較大貢獻[9-12],在旋流杯模型燃燒室流場研究過程中依次發(fā)展了受限空間單頭部模型燃燒室、真實火焰筒結構的單頭部模型燃燒室和多頭部模型燃燒室。首先,借助受限空間單頭部模型燃燒室在常溫常壓條件下開展了限制域有/無、限制域尺寸和旋流器參數(shù)等對流場結構的影響。其次,對于采用RQL燃燒技術的旋流杯模型燃燒室,頭部的旋流與主燃孔、補燃孔和摻混孔射流相互作用對燃燒穩(wěn)定邊界、污染物排放和出口溫度分布質量等有重要影響,Hassa等[13]和Mohammad等[14]借助帶主燃孔、補燃孔和摻混孔、更接近真實燃燒室結構的光學模型燃燒室開展了旋流與射流混合過程及其流場特征的研究。最近,燃燒室頭部之間及其與側壁相互作用引起研究者的重視,Kao等[10]和Kwong等[15]借助直線型多頭部模型燃燒室開展了頭部個數(shù)、頭部間距、軸線安裝位置和側壁距離等對流場結構影響研究。

    圖2 辛辛那提大學旋流杯模型燃燒室試驗件[9-15]Fig. 2 Swirl cupmodel combustor rig in University of Cincinnati [9-15]

    為進一步模擬燃燒室的真實進氣條件,旋流杯模型燃燒室試驗件向高溫高壓方向發(fā)展。德國宇航局(DLR)[16]和法國宇航局(ONERA)[17]在高溫高壓模型燃燒室設計及測量方面一直處于領先地位,國內的研究機構近年來也在積極探索高溫高壓旋流杯模型燃燒室及其對應的測量方法。Willert等[16]在2001年已完成了2 MPa和850 K進氣條件的高溫高壓光學可視模型燃燒室試驗件設計,并采用PIV開展了流場測量實驗,通過燃燒室出口的臨界節(jié)流孔板控制空氣流量和壓力,主流空氣按2∶1的比例分別進入燃燒室頭部和內側可視窗的冷卻狹縫。西北工業(yè)大學的Xiao等[18]采用真實火焰筒結構,利用從燃燒室試驗件尾部進光側面拍攝的方法,在0.55 MPa和300 K的進氣條件下分析了主燃孔和摻混孔與旋流流場的耦合過程及其對流場結構的影響。

    1.2 分區(qū)分級耦合燃燒光學模型燃燒室

    分區(qū)分級耦合燃燒技術是低排放燃燒室和高溫升燃燒室普遍采用的燃燒組織方法,也是目前國內外航空發(fā)動機燃燒室研制致力突破和創(chuàng)新的關鍵技術之一[1-2,5]。下面針對目前國內外在新型燃燒室研制過程中普遍采用的典型光學模型燃燒室試驗件進行分析和總結。

    1) TAPS模型燃燒室。美國密歇根大學為研究TAPS燃燒室在不同狀態(tài)下的流場特征(剪切層、回流區(qū)等),設計了高壓燃燒室模型[19-20],如圖5所示。模型燃燒室的圓筒形高壓機匣內安裝了TAPS燃燒室的燃油噴嘴和火焰筒,燃燒室機匣耐壓能力為1 MPa。為便于光學測量,燃燒室壁面設計了3個可視窗,其中上側可視窗用于圖像采集,水平方向成180°的2個窗口分別用于讓激光進入燃燒室和激光經過測試區(qū)后透過燃燒室,從而減小激光在燃燒室可視窗之間反射而產生的信號噪聲。為保證進入燃燒室的空氣分布均勻,在主燃燒室的上游設計了一個整流段,依次采用多孔孔板、玻璃珠和蜂窩結構對來流空氣進行整流。同時,燃燒區(qū)下游的火焰筒上分布有一排冷卻孔,這種冷卻布置使火焰筒的冷卻過程與真實的燃氣輪機燃燒室相同,進一步保證了模型燃燒室內流場與真實燃燒室內流場的相似性。試驗過程中,進入TAPS燃燒室頭部的空氣量和冷卻空氣量比例沒有進行測試,而是通過旋流器和冷卻孔的有效面積計算得到。

    圖5 密歇根大學TAPS燃燒室模型[19]Fig. 5 TAPS model combustor of University of Michigan[19]

    2) Lean-burn燃燒室模型。德國宇航局(DLR)在研究RR公司的Lean-burn低污染燃燒室過程中設計了如圖6所示的光學研究試驗臺和燃燒室[21-22]。為獲取燃燒室在實際工況下的燃燒流場數(shù)據(jù),模型燃燒室耐壓能力高達2 MPa,進口空氣溫度為850 K,空氣質量流量為10 kg/s。為滿足光學可視化測量光路布置的需要,整個試驗件測量段高壓機匣設計成矩形而不是圓形,方便相干反斯托克斯拉曼散射(CARS)、LDV、PIV和OH-PLIF等光學診斷設備在此試驗件上的使用。可視窗玻璃可以沿流動方向調整位置,以適應不同軸向位置速度場測量過程中發(fā)射光源光路和接受信號光路的布置。為滿足測試區(qū)尺寸要求,德國宇航局開發(fā)了具有較大可視窗口的BOSS光學可視模型燃燒室,如圖7所示[23]。

    圖6 德國宇航局(DLR)燃燒室模型[22]Fig. 6 Model combustor of DLR[22]

    圖7 BOSS光學可視模型燃燒室(DLR)[23]Fig. 7 Big optical single sector model combustor (DLR)[23]

    3) Lean-stage燃燒室。日本宇航局(JAXA)建立了單頭部燃燒室高壓噴霧試驗裝置,如圖8所示,最大工作狀態(tài)參數(shù)為常溫條件下壓力1 MPa[24-29]。其試驗段內徑為310.5 mm,高壓可視窗直徑分別約為150和200 mm,可視窗的石英玻璃厚度為40 mm。在可視窗的內側設計了輔助空氣氣簾系統(tǒng)用于清潔可視窗,速度場和粒徑測量采用PDPA方法,霧錐幾何特性和燃油分布測量采用PMie方法。此模型燃燒室可以開展高壓試驗,但沒有考慮火焰筒限制域對流場特性的影響。隨后,在進行常溫常壓點火和熄火特性試驗過程中,火焰筒采用直徑為125 mm的石英玻璃筒,用于試驗過程中的火焰探測。

    4) 南京航空航天大學TAPS/MIDL燃燒室模型。南京航空航天大學在研制一種帶多點燃油直接噴射雙環(huán)預混旋流頭部的貧油預混預蒸發(fā)(LPP)燃燒室時,忽略燃燒室機匣、火焰筒形狀及其冷卻,將模型燃燒室設計成一個單層的圓筒形燃燒室[30-32]。為滿足氣動霧化場光學測量過程中光路布置的需求,如測量燃燒室不同橫向和縱向截面的速度分布,燃燒室筒體采用全透明的石英玻璃筒體,并在測量段末端設置石英玻璃觀察窗。在中心截面的速度場測量過程中,激光通過測量段末端的可視窗進入燃燒室,利用CCD相機通過圓筒形可視窗采集粒子圖像;而在橫截面的速度場測量過程中,激光通過圓筒形可視窗進入燃燒室,在測量段末端的可視窗進行粒子圖像采集。在進行燃燒室熱態(tài)流場試驗時,為防止燃燒火焰發(fā)光對測量結果的影響,利用波長為532 nm、帶寬為5 nm的濾光片對火焰光干擾信號進行屏蔽。

    5) 北京航空航天大學TeLESS模型燃燒室。北京航空航天大學在TeLESS中心分級低排放燃燒室的研制過程中開發(fā)了光學可視模型燃燒室試驗件[33],如圖10所示。TeLESS模型燃燒室采用四面開窗的方形結構,上下側的火焰筒保留了原始火焰筒的初始結構特征和氣膜冷卻結構,并且預留了相對較大的側壁可視窗,用以采集盡可能大的區(qū)域的氣動霧化場和燃燒場信息。此模型燃燒室考慮了限制域、火焰筒幾何形狀及氣膜冷卻等因素,在保持燃燒室內氣動霧化和燃燒場真實性方面具有一定的優(yōu)勢。

    圖8 日本JAXA高壓試驗裝置[24]Fig. 8 High pressure combustor facility of Japan Aerospace Exploration Agency[24]

    圖9 南航LPP燃燒室試驗方案[30-32]Fig. 9 Lean premixed prevaporized combustor facility of Nanjing University of Aeronautics and Astronautics[30-32]

    圖10 TeLESS光學可視模型燃燒室[33]Fig. 10 Optically accessible model combustor of TeLESS[33]

    6) 中國科學院工程熱物理研究所分層部分預混模型燃燒室。針對組合式燃油霧化噴嘴對高溫高壓氣動霧化特性的測量需求,中國科學院工程熱物理研究所研制了燃油噴嘴高壓高溫霧化特性平面激光測量系統(tǒng),開發(fā)了高溫高壓光學可視模型燃燒室試驗件[34],如圖11所示。采用內層承溫外層承壓的設計理念,火焰筒內側設計氣簾用于吹掃火焰筒可視窗,減小示蹤粒子及液滴對可視窗的污染;模型試驗件中采用方形火焰筒,火焰筒橫截面積保持與燃燒室單頭部橫截面積相同,同時保持了火焰筒收縮段長度及收縮比例。應用于采用中心分級燃燒技術的分層部分預混燃燒室時,由于分層部分預混燃燒室取消了主燃孔和補燃孔,約50%~70%的空氣由燃燒室頭部的旋流器進入,燃燒室內的流場結構主要由旋流器氣動結構參數(shù)決定,火焰筒壁面的氣膜冷卻孔和摻混孔對燃燒室內的流場結構影響較小,并考慮到火焰筒的加工難度和光學可視需求,光學火焰筒省略了氣膜冷卻孔和摻混孔。目前已經利用此模型燃燒室試驗件對分層部分預混燃燒室進行了大量的氣動霧化場試驗,也為從事發(fā)動機研究相關的高等院校、科研院所和工業(yè)部門的新型高溫升和低排放燃燒室研制提供了氣動霧化場測量試驗服務和數(shù)據(jù)支撐。

    圖11 高溫高壓光學可視模型燃燒室試驗件[34]Fig. 11 High temperature and high pressure model combustor with optical access[34]

    7) 中國空氣動力研究與發(fā)展中心三級軸向旋流模型燃燒室。中國空氣動力研究與發(fā)展中心和西北工業(yè)大學合作,針對采用三級軸向旋流器的航空發(fā)動機燃燒室設計了高溫高壓光學可視模型燃燒室試驗件[35]。供氣系統(tǒng)提供常溫或加熱的壓縮空氣來流,經燃燒室后水平向側后方45°轉向流入排氣管路,管路轉向是為了在正對燃燒室出口的方向設置石英玻璃窗,以進行光學測量,如圖12所示。該試驗件可以承受真實燃燒室來流條件,目前已經開展了全狀態(tài)進氣壓力2.78 MPa和進氣溫度813 K條件下的PIV流場測量。

    圖12 三級軸向旋流模型燃燒室[35]Fig. 12 Three-staged axial swirler model combustor[35]

    以上典型的光學模型燃燒室為發(fā)動機燃燒室研制過程中燃燒室頭部方案設計、分區(qū)分級氣動霧化場和燃燒組織機理認識、霧化性能和部分燃燒性能評估提供了重要支撐。目前已服役的GE航空公司的TAPS系列低排放燃燒室和RR航空公司的Trent系列發(fā)動機上的低排放燃燒室等都離不開研制前期在光學模型燃燒室上的大量基礎研究試驗。我國在研制民用發(fā)動機低排放燃燒室和高推重比發(fā)動機高溫升燃燒室的過程中,已在光學模型燃燒室上開展了部分基礎試驗,從支撐在研發(fā)動機燃燒室和探索未來新型燃燒室的角度考慮,還需要進行更加全面的氣動霧化和燃燒組織機理試驗、數(shù)據(jù)積累和針對具體工程問題的應用試驗等。

    1.3 流場診斷技術在光學模型燃燒室上應用的難點分析

    根據(jù)前面對光學模型燃燒室發(fā)展趨勢的分析,從光學模型燃燒室設計的角度考慮,流場診斷技術在光學模型燃燒室上應用的主要挑戰(zhàn)是可視窗空間布局、可視窗的清潔和可視窗的燒蝕等問題??梢暣翱臻g布局局限性的分析已在典型的光學模型燃燒室中介紹。燃燒室內流場測量試驗過程中,需要向燃燒室入口空氣中添加示蹤粒子,包括液滴和固體示蹤粒子,示蹤粒子材質和顆粒粒徑等選擇可參考文獻[36],而在高溫冷態(tài)和燃燒條件下,由于液滴示蹤粒子會快速蒸發(fā),故只能采用耐高溫的固體示蹤粒子。無論采用液滴或固體示蹤粒子,進入燃燒室后都可能部分附著于可視窗表面,對光學模型燃燒室的可視窗產生污染,而發(fā)動機燃燒室內的旋流環(huán)境進一步加劇了示蹤粒子在可視窗表面的附著幾率。目前普遍通過在可視窗內側設置氣簾來降低示蹤粒子對可視窗的污染。進一步,燃燒條件下的流場測量過程中,在相對較小的發(fā)動機工況時,可視窗壁面易產生積炭;同時,對于采用液體燃料的燃燒室,存在部分液滴撞擊到可視窗表面的問題,一方面影響可視窗光學通透性,另一方面也會加劇可視窗表面的積炭問題。此時氣簾設計需考慮對可視窗的加熱作用,加速可視窗表面附著液滴的蒸發(fā),減小液滴及隨后的積炭對可視窗的污染。而在發(fā)動機大工況條件下,火焰局部溫度達到2000 ℃以上,可視窗表面易被高溫燃氣燒蝕,可視窗表面的積炭和燒蝕都影響流場測量過程中的光學信號輸入/采集,給流場測量帶來誤差,可視窗壁面污染嚴重時甚至難以開展流場測量試驗。此時吹掃氣簾設計需考慮對可視窗的冷卻作用。綜上所述,故在可視窗內外兩層分別設置吹掃氣簾。

    2 限制域對流場特性的影響

    發(fā)動機燃燒室內的流場組織是受限空間內的旋流流動過程,因此,限制域是影響旋流流場結構的重要因素,是在進行模型燃燒室試驗件設計過程中需要考慮的重要因素之一[37]。為深入認識旋流杯模型燃燒室內流場組織過程和驗證燃燒模型,辛辛那提大學的Fu和Jeng等[9]系統(tǒng)開展了限制域對旋流杯模型燃燒室內的流場特性影響研究。在旋流器壓降為4%的常溫常壓條件下,采用LDV測量燃燒室內的速度場,方形燃燒室寬度分別為76.2 mm(3.0 inch)、101.6 mm (4.0 inch)、114.4 mm (4.5 inch)、127.0 mm (5.0 inch)、152.4 mm (6.0 inch)和開放空間。如圖13所示,研究結果表明,在受限空間工況下,流場中都存在角渦回流區(qū)(Corner Recirculation Zone,CORZ),角渦回流區(qū)的尺寸和旋流在火焰筒上的附著點受火焰筒寬度與旋流器直徑的比值影響。在限制域尺寸較大時,旋流在火焰筒上的附著點靠近下游,這也導致燃燒室內存在2個中心回流區(qū)(Center Toroidal Recirculation Zone,CTRZ)。而在限制域尺寸較小時(3.0和4.0 inch),旋流產生的中心低壓區(qū)形成一個中心回流區(qū);限制域由4.0 inch縮小到3.0 inch時,CTRZ尺寸減小。流場模式由2個CTRZ轉換成1個CTRZ時,經過詳細的試驗發(fā)現(xiàn)流場模式轉換的臨界點為限制域尺寸4.3 inch,并且在流場模式轉換臨界點工況下,流場隨機呈現(xiàn)2種不同的速度分布,如圖14所示。目前,尚未見限制域對中心分級燃燒室流場特性研究的論文發(fā)表。

    圖13 限制域對流場的影響[9]Fig. 13 Confinement effects on flow field[9]

    圖14 流場模式轉換臨界點時的兩種流場結構[9]Fig. 14 Flow fields with two metastable modes[9]

    3 燃燒室頭部之間流場耦合作用

    燃燒室由單管燃燒室發(fā)展到環(huán)管燃燒室,進一步發(fā)展到目前在燃氣輪機上普遍采用的環(huán)形燃燒室,燃燒室頭部之間的耦合作用不斷增強[4]。最近,高推重比軍用發(fā)動機的高溫升燃燒室和大涵道比民用發(fā)動機的低排放燃燒室頭部進氣量比例由富油頭部設計理念的傳統(tǒng)燃燒室的15%~25%增加到60%~80%,旋流器直徑大幅增加,進一步增強了燃燒室頭部之間的耦合作用,因此,認識燃燒室頭部之間耦合作用對燃燒室模型試驗件設計及其試驗研究均至關重要。

    辛辛那提大學的Jeng等對旋流杯模型燃燒室直線陣列布置的頭部之間旋流相互作用對流場的影響開展了研究[10,38-39]。在旋流器壓降為4%的常溫常壓條件下,采用LDV測量燃燒室內的速度場。研究結果表明,在三頭部和五頭部試驗件中各個頭部的CTRZ尺寸和強度并不相同(見圖15)。相對于兩側的頭部,中心頭部的CTRZ強度相對較小,而CTRZ尺寸相對較大;中心頭部的流場組織過程類似于開放空間中的旋流器。產生這種現(xiàn)象的原因是頭部之間旋流空氣的相互作用,中心頭部的旋流強度減弱,其旋流方向主要受中心第一級旋流器決定,而兩側頭部中的旋流更類似于單頭部受限空間中的流場組織過程。

    圖15 旋流杯模型燃燒室頭部之間耦合相互作用對流場結構的影響(軸向速度)[38]Fig. 15 The effect of sector-to-sector interactions on flow field of swirl cup combustor (axial velocity)[38]

    為分析燃燒室頭部之間的氣動作用,Jeng研究了頭部間距(S)和側邊旋流器中心到火焰筒壁面的距離(Dw)對直線型五頭部旋流杯模型燃燒室試驗件(見圖16)冷態(tài)流場的影響[39],相對于S=2D(D為旋流器直徑)的流場(從中心頭部向外呈現(xiàn)“大-小-大”周期性的非均勻CTRZ分布),當頭部間距減小到1.75D時,回流區(qū)空間布局不變,而當頭部間距增加到2.50D,回流區(qū)空間布局發(fā)生明顯變化,中心頭部向外呈現(xiàn)出 “小-大-小”周期性的非均勻CTRZ分布,繼續(xù)增加頭部間距到2.75D時,回流區(qū)空間布局與2.50D時相同,如圖17和18所示。為進一步深入分析五頭部旋流杯模型燃燒室試驗件中周期性非均勻回流區(qū)形成的原因,Jeng結合限制域和旋流器深入火焰筒深度2個參數(shù)詮釋了2種旋流流場組織模式——貼壁旋流模式和自由射流旋流模式的形成機理,如圖19所示。

    圖16 旋流杯五頭部模型燃燒室試驗件[39]Fig. 16 Linear-arranged 5-swirler array of swirl cup combustor[39]

    圖17 旋流杯模型燃燒室頭部間距對流場結構的影響[39]Fig. 17 Effect of inter-swirler spacing on flow field of swirl cup combustor[39]

    圖18 中心回流區(qū)空間布局示意圖[39]Fig. 18 Periodic flow patterns of center recirculation zone[39]

    圖19 旋流杯模型燃燒室內流動模式分析[39]Fig. 19 Flow pattern sketches in swirl cup combustor[39]

    北京航空航天大學的林宇震等[40]通過對比單頭部受限空間和三頭部旋流杯模型燃燒室內的流場結構,分析了邊界條件對燃燒室內渦分布和耗散速率的影響,結果表明,在三頭部燃燒室中,渦結構在旋流器下游較長距離內都更能保持連貫性。中國科學院工程熱物理研究所的徐綱和北京航空航天大學的林宇震合作研究了分區(qū)分級新型燃燒室頭部之間的流場耦合作用(見圖20)[41]。借助三頭部光學模型試驗件開展了氣動霧化特性試驗,中心頭部的回流區(qū)尺寸相對于側邊頭部較小,不同于旋流杯模型燃燒室三頭部試驗件的流場分布結果,這說明燃燒室頭部之間的相互作用不僅受到燃燒室頭部之間間距的影響,也取決于燃燒室內的旋流器及其流場特性。

    圖20 某分級分區(qū)低排放燃燒室三頭部試驗件內回流區(qū)(Plane 1:中心頭部,Plane 5:側邊頭部)[41]Fig. 20 Recirculation zone in the three-sector rig of a internal staged low emission combustor (Plane 1: center dome, Plane 5: side dome)[41]

    4 發(fā)動機燃燒室流場研究進展

    除發(fā)動機燃燒室的全環(huán)、火焰筒及其多孔等結構特征外,燃燒室內的流場組織過程和流場特征也是模型燃燒室試驗件設計過程中需要考慮的重要因素。目前在役和在研的發(fā)動機燃燒室主要包括2種類型:基于頭部富油設計的傳統(tǒng)旋流杯燃燒室和采用分區(qū)分級擴散/預混耦合燃燒技術的中心分級燃燒室[1],因此,本文主要針對這2種類型的發(fā)動機燃燒室流場特性研究進展進行分析。

    4.1 傳統(tǒng)旋流杯燃燒室

    旋流杯燃燒室的頭部油/氣混合裝置由壓力霧化噴嘴和周圍兩層旋向相反的旋流空氣組成,由于其良好的霧化和燃燒性能,在發(fā)動機上得到廣泛應用,如CF6-80、F404、CFM56、GE90、F101、F110、F414及T700等[42]。自20世紀80年代以來,國內外學者針對旋流杯燃燒室的流場組織開展了大量研究,國外主要代表有GE公司的Mongia等[8]、辛辛那提大學Jeng等[9,42]、加州大學歐文分校的McDonell等[43]和德國宇航中心(DLR)的Hassa等[13],國內的主要代表有北京航空航天大學、南京航空航天大學、西北工業(yè)大學3個航空航天院校、中國航空發(fā)動機集團公司下屬的動力所和中國科學院等。傳統(tǒng)旋流杯燃燒室的典型流場特征如圖21所示。在燃燒室中心形成中心回流區(qū)(CTRZ),在燃燒室火焰筒與頭部安裝邊交接處形成角渦回流區(qū)(CORZ)。影響旋流杯燃燒室流場的因素包括旋流器的氣量分配、旋流強度、旋流器軸向/徑向位置和出口套筒結構等氣動結構參數(shù),在此方面已有大量成果發(fā)表,本文不再贅述。

    圖21 旋流杯燃燒室典型流場結構Fig. 21 The typical flow field structure of swirl cup combustor

    模型燃燒室設計過程中還需要考慮的另一個因素是燃燒化學反應釋熱對流場結構的影響(見圖22)[11,44]。目前得到的一致結論是:燃燒化學反應釋熱引起主燃區(qū)空氣膨脹,導致回流區(qū)尺寸減小。因此,根據(jù)需要試驗目標,確定是在冷態(tài)條件下還是在熱態(tài)燃燒化學反應條件下開展流場測試試驗。

    圖22 旋流杯燃燒室冷/熱態(tài)流場對比[11]Fig. 22 Comparison between non-reacting and reacting flow fields in swirl cup combustor[11]

    進一步的,在模擬燃燒室真實結構方面,導流板、主燃孔和摻混孔等也是模型燃燒室設計過程中需要考慮的因素,Mohammad等和Hassa等[12-14]在帶導流板、主燃孔和摻混孔的真實燃燒室火焰筒內開展了旋流杯燃燒室流場特性研究,其流場結構(見圖23)明顯不同于簡化火焰筒的模型燃燒室流場,角渦回流區(qū)消失,中心回流區(qū)被主燃孔射流截斷,中心回流區(qū)尺寸明顯減小,在主燃孔和摻混孔下游形成尾流低速區(qū)。

    圖23 旋流杯燃燒室真實火焰筒結構的流場特征[12]Fig. 23 The flow field in swirl cup combustor with real liner[12]

    4.2 中心分級燃燒室

    中心分級燃燒室是目前先進燃燒室的代表和未來的發(fā)展方向,目前的高溫升燃燒室和低排放燃燒室普遍采用中心分級燃燒室技術,借助燃油分級和多級旋流器實現(xiàn)分區(qū)分級擴散/預混耦合燃燒,分為主燃級和預燃級,以適應發(fā)動機全工況范圍內燃燒性能的需求。國內外學者均對中心分級燃燒室流場開展了大量基礎研究,促進了對中心分級燃燒室內流場組織過程和流場特性的認識。圖24為中心分級主燃燒室的典型流場結構示意圖[45],燃燒室內存在主回流區(qū)(PRZ)、角回流區(qū)(CRZ)、臺階回流區(qū)(LRZ)和剪切層(Shear layer)。

    圖24 中心分級燃燒室典型流場結構[45]Fig. 24 The typical flow field of internal staged combustor[45]

    GE公司的TAPS燃燒室是最早在發(fā)動機上成功應用的中心分級燃燒室方案,目前已經發(fā)展到第三代(TAPS I,TAPS II和TAPS III)[46]。燃燒室頭部的燃燒組織采用預燃級與主燃級同心共軸設計,預燃級為兩級旋流器,主燃級為一級徑向旋流器,預燃級采用富油擴散燃燒,主燃級采用貧油預混預蒸發(fā)燃燒。

    圖25為Sulabh等[19-20]利用PIV得到的TAPS燃燒室在某工況下的冷態(tài)、熱態(tài)流場??梢钥闯?在冷態(tài)工況下,燃燒室內形成了尺寸相對較大的主回流區(qū)(PRZ),主回流區(qū)幾乎充滿了整個燃燒室。而在預燃級單獨工作和主/預燃級同時工作的2種熱態(tài)工況下,燃燒室中心線上的軸向速度由負變?yōu)檎?,主回流區(qū)尺寸明顯減小,并且主回流區(qū)被分為2個獨立的環(huán)形區(qū)域。冷/熱態(tài)流場對比結果表明,熱態(tài)條件下的熱釋放改變了回流區(qū)形狀和速度場結構。圖26為TAPS燃燒室冷、熱態(tài)工況下的湍流強度和雷諾應力分布。冷態(tài)工況下,相對于主回流區(qū),主/預旋流空氣出口位置湍流強度明顯較大,在回流區(qū)邊界處存在較大的剪切應力;熱態(tài)工況下,預燃級火焰的熱釋放使燃燒室中心的湍流強度大幅增加,預燃級旋流空氣與主回流區(qū)之間、主燃級旋流空氣與預燃級旋流空氣之間的區(qū)域都形成了剪切層。

    圖25 TAPS燃燒室冷態(tài)、熱態(tài)流場[19]Fig. 25 Non-reacting and reacting flow fields in TAPS combustor[19]

    圖26 TAPS燃燒室冷/熱態(tài)工況下的湍流強度和雷諾應力[19]Fig. 26 Turbulence intensity and shear strain rate in TAPS combustor under non-reacting and reacting conditions[19]

    RR公司在Engine 3E項目中研制了多種類型的中心分級燃油噴嘴用于Lean-burn燃燒室[47],預燃級方案采用預膜空氣霧化噴嘴或離心噴嘴同心旋流空氣霧化技術;主燃級燃油霧化方案都采用預膜式空氣霧化技術。Meier等[48]針對Lean-burn燃燒室研究了預燃級旋流強度對流場結構的影響,其中弱旋流的噴嘴下游未形成回流區(qū),而強旋流噴嘴下游形成了中心回流區(qū),如圖27所示。地面起動點/熄火試驗表明,強旋流噴嘴的燃燒室有較好的點火能力,但顯示出較差的貧油燃燒穩(wěn)定性,而弱旋流噴嘴的點/熄火燃燒穩(wěn)定性趨勢與此相反。

    圖27 RR公司Lean-burn燃燒室內流場[48]Fig. 27 Flow field in Lean-burn combustor of Rolls-Royce company[48]

    JAXA和川崎重工(KHI)聯(lián)合研制了Lean-Staged低排放燃燒室[28-29]。其采用四級旋流器,主燃級采用噴射成膜霧化方式,預燃級采用通道內預膜或旋流杯預膜空氣霧化方式,典型流場結構如圖28(a)所示。Kobayashi等在優(yōu)化Lean-Staged低排放燃燒室點火性能的過程中研究了預燃級內旋流器角度對燃燒穩(wěn)定性的影響,結果表明預燃級內級旋流器角度由0°增加到45°時,主回流區(qū)尺寸略有增加,預燃級內級旋流器通過改變燃油空間分布特性而影響點火性能。進一步的,為阻止主燃級旋流空氣對預燃級霧錐的耦合作用,在預燃級和主燃級旋流器周向布置了狹縫,狹縫對流場結構的影響如圖28(b)所示,通過增加狹縫改善了小推力工況下的燃燒穩(wěn)定性。

    圖28 Lean-Staged燃燒室流場[29]Fig. 28 Flow field of Lean-Staged combustor[29]

    中國科學院工程熱物理研究所針對分層部分預混低排放燃燒室和高溫升燃燒室的流場特性及組織機理開展了研究[49-51]。圖29對比了單頭部常壓下預燃級內級旋流器是否有旋對燃燒室流場結構的影響[49]??梢钥闯?,預燃級有旋時可以形成幾乎占滿頭部的回流區(qū),旋流器出口會形成小的臺階回流區(qū),同時燃油幾乎全部分布在中心回流區(qū)里;而預燃級無旋的回流區(qū)則被高速軸向射流分成2個小回流區(qū),燃油分布于中心高速射流區(qū)。進一步研究了主燃級旋流角度(40°、45°和50°)對流場結構的影響[50],獲得了如圖30所示的時均流場分布??梢钥闯?,隨著葉片角度增加,主回流區(qū)位置向上移動,尺寸增大,主燃級回流氣與預燃級進氣的相互作用增強,但主燃級回流區(qū)與預燃級回流區(qū)在中心線附近有重合。

    圖29 預燃級內級旋流器有旋/無旋對流場及燃油分布的影響[49]Fig. 29 Effect of inner swirler angle on flow field and spray pattern[49]

    圖30 主燃級旋流角度對流場的影響[50]Fig. 30 Effect of main swirler angle on flow field[50]

    5 結 論

    光學模型燃燒室試驗件是認識發(fā)動機燃燒室內流場組織機理和流場結構特征的基礎,在發(fā)動機燃燒室的研制過程中發(fā)揮了重要作用。目前,針對復雜旋流組織機理的光學模型燃燒室試驗件已經形成相對一致的認識,而針對發(fā)動機燃燒室工程研制的光學模型燃燒室試驗件設計尚未形成統(tǒng)一的意見,但完全保持燃燒室結構和進出口試驗條件,對于測量方法的環(huán)境適用性和光學診斷光路布置而言,在現(xiàn)實燃燒室流場測量中是不可行的。當前,對于發(fā)動機燃燒室內流場的測試,在測量方法、進氣條件、燃燒室結構和光路布置等方面都需要通過權衡、盡量得到真實的燃燒室流場信息,從而進一步認識燃燒室內復雜旋流組織過程。

    光學模型燃燒室的應用可以分為3個主要方面:1) 發(fā)動機燃燒室工程應用。深入認識發(fā)動機燃燒室的氣動、霧化和燃燒等過程,支撐新型燃燒室研發(fā)及燃燒相關問題的解決。2) 氣動熱力學模型和數(shù)值仿真工具開發(fā)。用于驗證湍流模型、霧化模型和燃燒模型等,利用試驗數(shù)據(jù)開發(fā)氣動霧化和燃燒性能預測模型,形成設計工具,支撐發(fā)動機燃燒室的研發(fā)。3) 解決燃燒科學問題。詮釋復雜旋流燃燒過程中的流場組織、霧化、蒸發(fā)、混合和燃燒組織機理。

    光學模型燃燒室以上3個方面的應用需求對模型燃燒室的進氣條件和燃燒室結構真實性的要求是逐次降低的,即對能否真實反映發(fā)動機燃燒室內流場組織過程的要求是逐次降低的,但是低工況下或簡化模型燃燒室上得到的氣動熱力學模型和數(shù)值仿真工具應用于發(fā)動機燃燒室設計時,需要考慮其適用范圍或進行外推驗證。受試驗條件、光路布置、測量方法和試驗成本等限制時,需要從以下幾個方面考慮:1) 所研究的燃燒室本身的特點。例如:對于采用分區(qū)分級耦合燃燒技術的先進燃燒室,其燃燒極大程度取決于頭部進氣的空氣動力學特性,故在模型燃燒室設計過程中,燃燒室頭部結構及其與頭部流場組織密切相關的因素應盡量與真實發(fā)動機燃燒室結構相同,此時火焰筒壁面的冷卻小孔(如氣膜孔)和位于主燃區(qū)下游的摻混孔等在開展流場組織機理研究時可忽略;對于傳統(tǒng)旋流杯燃燒室,通過主燃孔和補燃孔的高速射流都與主燃區(qū)的旋流流動強耦合,從而影響流場特性,因此,在開展旋流杯燃燒室頭部流場組織機理研究時,火焰筒壁面的主燃孔和補燃孔需要保留。2) 模型燃燒室的研發(fā)目標。如果模型燃燒室研發(fā)的目的是用于燃燒基礎科學問題或模型/數(shù)值仿真工具方法研究,在分析所研究的燃燒室特點基礎上,首要考慮的因素是模型燃燒室邊界條件清晰,其次是待研究的流場、霧化特性和燃燒特性等物理量對應的測量方法和光路布置要求,最后考慮進氣條件和模型燃燒室結構;如果模型燃燒室研發(fā)的目的是燃燒室工程研制,在分析所研究的燃燒室特點基礎上,首要考慮的因素是所處的研制階段或技術成熟度,例如在役發(fā)動機燃燒室使用過程的燃燒問題、在研發(fā)動機燃燒室技術成熟度達到TRL3及TRL3以上性能的驗證過程中的燃燒問題和預研新型發(fā)動機燃燒室在TRL3以下時借助模型燃燒室認識氣動、霧化和燃燒過程并評估初步燃燒性能等;其次需考慮的是具體的燃燒問題,例如點火、熄火、積炭、火焰筒燒蝕、出口溫度溫度分布質量和污染物排放指數(shù)等。

    國內外在傳統(tǒng)旋流杯燃燒室和分區(qū)分級耦合燃燒燃燒室的流場組織機理和流場特性方面都已經開展了大量的基礎和應用基礎研究,得到了旋流杯燃燒室和分區(qū)分級耦合燃燒室的典型流場結構,但由于對旋流、剪切層形成等過程中的渦核進動、渦形成和渦破碎等微觀機理認識不足,目前還難以通過燃燒室頭部結構和旋流器氣動參數(shù)等對流場結構進行精細化設計。

    猜你喜歡
    視窗旋流器旋流
    視窗
    工會博覽(2022年33期)2023-01-12 08:53:38
    重介質旋流器選煤技術在我國的創(chuàng)新發(fā)展與應用
    選煤技術(2022年3期)2022-08-20 08:38:50
    視窗
    工會博覽(2022年17期)2022-07-15 07:26:14
    雙級徑向旋流器對燃燒性能的影響
    燒嘴旋流器優(yōu)化設計計算
    科學家(2021年24期)2021-04-25 13:25:34
    視窗
    中國水利(2017年9期)2017-06-07 08:23:27
    上旋流厭氧反應器在造紙廢水處理中的應用
    一級旋流偏置對雙旋流杯下游流場的影響
    微氣泡旋流氣浮選設備的研究與應用
    旋流脫氣性能影響因素的CFD模擬
    化工進展(2015年6期)2015-11-13 00:30:04
    18禁在线播放成人免费| 黄片无遮挡物在线观看| 能在线免费看毛片的网站| 大话2 男鬼变身卡| 91精品国产国语对白视频| 日韩,欧美,国产一区二区三区| 嫩草影院入口| 国产高清不卡午夜福利| 亚洲,一卡二卡三卡| 久久精品国产自在天天线| 视频区图区小说| 亚洲第一av免费看| 51国产日韩欧美| 亚洲国产色片| 一本色道久久久久久精品综合| 久久免费观看电影| 97超碰精品成人国产| 乱人伦中国视频| 在线免费观看不下载黄p国产| 嫩草影院新地址| 国产精品国产三级国产av玫瑰| 男女啪啪激烈高潮av片| 久久久久久久久久久免费av| 久久国产精品男人的天堂亚洲 | 老熟女久久久| 国产 精品1| 日日摸夜夜添夜夜爱| 国产亚洲欧美精品永久| 天天躁夜夜躁狠狠久久av| 久久鲁丝午夜福利片| 亚洲av成人精品一区久久| 国产精品一二三区在线看| 国产免费一区二区三区四区乱码| 夫妻午夜视频| 欧美bdsm另类| 3wmmmm亚洲av在线观看| 五月伊人婷婷丁香| 制服丝袜香蕉在线| 十八禁网站网址无遮挡 | 亚洲天堂av无毛| 欧美三级亚洲精品| 人妻制服诱惑在线中文字幕| 久热这里只有精品99| 久久精品国产鲁丝片午夜精品| 国产中年淑女户外野战色| 亚洲av免费高清在线观看| 亚洲欧美日韩卡通动漫| 男女国产视频网站| 亚洲国产av新网站| 自拍偷自拍亚洲精品老妇| 免费大片黄手机在线观看| 亚洲不卡免费看| 国产老妇伦熟女老妇高清| 日本猛色少妇xxxxx猛交久久| 亚洲精品日韩av片在线观看| 秋霞伦理黄片| 亚洲av综合色区一区| 亚洲精品国产色婷婷电影| 街头女战士在线观看网站| 人人妻人人澡人人爽人人夜夜| 亚洲精品乱码久久久v下载方式| 一级毛片久久久久久久久女| 国产精品秋霞免费鲁丝片| 人人妻人人澡人人爽人人夜夜| 国产精品蜜桃在线观看| 久久久久人妻精品一区果冻| 日韩av不卡免费在线播放| 日本黄色日本黄色录像| 另类精品久久| 久久97久久精品| 亚洲av成人精品一二三区| h日本视频在线播放| 国产一区有黄有色的免费视频| 亚洲欧美清纯卡通| 中国三级夫妇交换| 男女边摸边吃奶| 看十八女毛片水多多多| 新久久久久国产一级毛片| 亚洲四区av| 国产精品女同一区二区软件| 国产黄频视频在线观看| av福利片在线观看| h视频一区二区三区| 亚洲成人手机| 美女国产视频在线观看| 亚洲精品亚洲一区二区| 午夜日本视频在线| 免费看光身美女| av专区在线播放| 精品亚洲成a人片在线观看| 水蜜桃什么品种好| 99九九线精品视频在线观看视频| 国产伦精品一区二区三区视频9| 亚洲国产精品一区二区三区在线| 在线观看免费视频网站a站| 国产精品久久久久久久久免| 六月丁香七月| 一级毛片电影观看| 国产午夜精品久久久久久一区二区三区| 国产高清三级在线| 韩国av在线不卡| 日韩成人av中文字幕在线观看| 久久久久精品久久久久真实原创| 国产精品久久久久久久电影| 丁香六月天网| 视频中文字幕在线观看| 少妇人妻久久综合中文| 亚洲精品aⅴ在线观看| 91成人精品电影| 99热6这里只有精品| 久久精品久久精品一区二区三区| 一区二区av电影网| 国产极品天堂在线| 国产有黄有色有爽视频| 欧美精品高潮呻吟av久久| 国产日韩欧美视频二区| 一二三四中文在线观看免费高清| 免费观看a级毛片全部| 午夜福利在线观看免费完整高清在| 国产在线免费精品| 免费人成在线观看视频色| 韩国高清视频一区二区三区| 日本欧美国产在线视频| 亚洲av电影在线观看一区二区三区| 亚洲av.av天堂| 国产精品无大码| 精品人妻偷拍中文字幕| 一级毛片黄色毛片免费观看视频| 高清午夜精品一区二区三区| 欧美最新免费一区二区三区| 亚洲,欧美,日韩| 免费人成在线观看视频色| 国产精品久久久久久久久免| 国产免费一区二区三区四区乱码| 亚洲精品国产色婷婷电影| 菩萨蛮人人尽说江南好唐韦庄| 蜜臀久久99精品久久宅男| 久久久国产精品麻豆| 日韩精品免费视频一区二区三区 | 91精品一卡2卡3卡4卡| 在现免费观看毛片| 精品少妇内射三级| 在线观看国产h片| 国产成人免费观看mmmm| 亚洲av欧美aⅴ国产| 日韩熟女老妇一区二区性免费视频| 国产一区有黄有色的免费视频| 久热久热在线精品观看| 久久久久精品性色| 亚洲,一卡二卡三卡| 夫妻午夜视频| 日韩中文字幕视频在线看片| 26uuu在线亚洲综合色| 亚洲av日韩在线播放| 国产精品99久久久久久久久| 亚洲国产毛片av蜜桃av| 国产精品蜜桃在线观看| 曰老女人黄片| 91精品伊人久久大香线蕉| 国产视频首页在线观看| 亚洲自偷自拍三级| 亚洲精品乱码久久久久久按摩| 久久久国产精品麻豆| av专区在线播放| 亚洲国产色片| 久久久国产一区二区| 日韩精品有码人妻一区| 国产黄频视频在线观看| 天天操日日干夜夜撸| av.在线天堂| 全区人妻精品视频| 久久精品国产亚洲av涩爱| 欧美变态另类bdsm刘玥| 国产免费又黄又爽又色| 亚洲av免费高清在线观看| 成人影院久久| 亚洲欧美日韩另类电影网站| 亚洲欧美清纯卡通| 最后的刺客免费高清国语| 欧美性感艳星| 精品久久国产蜜桃| 亚州av有码| 少妇人妻久久综合中文| 日本av免费视频播放| 婷婷色麻豆天堂久久| 久久免费观看电影| 亚洲欧美精品专区久久| 一区二区三区精品91| 男人舔奶头视频| 纵有疾风起免费观看全集完整版| 亚洲av日韩在线播放| 22中文网久久字幕| 日韩av不卡免费在线播放| 国产女主播在线喷水免费视频网站| 亚洲电影在线观看av| 久久久久久久久久人人人人人人| 三级国产精品欧美在线观看| 男女免费视频国产| 国产熟女欧美一区二区| 久久久久久久久久人人人人人人| 欧美三级亚洲精品| 交换朋友夫妻互换小说| 国产女主播在线喷水免费视频网站| 永久免费av网站大全| 欧美亚洲 丝袜 人妻 在线| 天天操日日干夜夜撸| 男女免费视频国产| 精品久久久久久久久亚洲| 最近最新中文字幕免费大全7| 一级毛片久久久久久久久女| 一本—道久久a久久精品蜜桃钙片| 五月玫瑰六月丁香| 亚洲自偷自拍三级| 亚洲精品国产色婷婷电影| 丰满乱子伦码专区| 国产男女内射视频| 国产亚洲最大av| 午夜久久久在线观看| 各种免费的搞黄视频| 国产一区亚洲一区在线观看| 丰满人妻一区二区三区视频av| 一区二区av电影网| 一区二区三区精品91| av播播在线观看一区| 只有这里有精品99| 国产精品成人在线| 汤姆久久久久久久影院中文字幕| 久久久亚洲精品成人影院| 亚洲真实伦在线观看| 91久久精品电影网| 亚洲一区二区三区欧美精品| 汤姆久久久久久久影院中文字幕| 免费高清在线观看视频在线观看| 寂寞人妻少妇视频99o| 国产成人精品福利久久| 免费av不卡在线播放| 国产精品一区二区三区四区免费观看| 国产成人一区二区在线| 亚洲精华国产精华液的使用体验| 国产av一区二区精品久久| 赤兔流量卡办理| 国产精品国产av在线观看| 亚洲国产精品999| 免费少妇av软件| 大又大粗又爽又黄少妇毛片口| 精品一品国产午夜福利视频| 精品一区二区三卡| 少妇被粗大的猛进出69影院 | 热re99久久精品国产66热6| 国产片特级美女逼逼视频| 国产精品久久久久久精品古装| 久久久国产欧美日韩av| 又爽又黄a免费视频| 观看免费一级毛片| 久热这里只有精品99| 亚洲欧美成人精品一区二区| 黄色一级大片看看| 国产免费福利视频在线观看| 精品人妻偷拍中文字幕| 街头女战士在线观看网站| 黄色日韩在线| 亚洲精品久久久久久婷婷小说| 最新的欧美精品一区二区| 夫妻性生交免费视频一级片| 精品一品国产午夜福利视频| 18禁动态无遮挡网站| 97在线人人人人妻| 亚洲精品国产av蜜桃| 精品99又大又爽又粗少妇毛片| 亚洲图色成人| 国产欧美日韩一区二区三区在线 | 日日啪夜夜撸| 少妇人妻精品综合一区二区| 午夜91福利影院| 日日爽夜夜爽网站| 国模一区二区三区四区视频| 五月天丁香电影| 亚洲精品成人av观看孕妇| 另类精品久久| 老女人水多毛片| 国产又色又爽无遮挡免| 精品酒店卫生间| 精品国产一区二区三区久久久樱花| 国内少妇人妻偷人精品xxx网站| 最新的欧美精品一区二区| 观看av在线不卡| 亚洲一级一片aⅴ在线观看| 色哟哟·www| 狂野欧美激情性xxxx在线观看| 一级毛片 在线播放| 午夜激情久久久久久久| 一级毛片电影观看| 久久国产亚洲av麻豆专区| 我要看黄色一级片免费的| 在现免费观看毛片| 人妻系列 视频| 狂野欧美白嫩少妇大欣赏| 99热6这里只有精品| 国产精品久久久久久精品古装| 成年av动漫网址| 日本欧美视频一区| 久久久精品免费免费高清| 一个人看视频在线观看www免费| 黄色配什么色好看| h日本视频在线播放| 如日韩欧美国产精品一区二区三区 | 国产精品欧美亚洲77777| 少妇裸体淫交视频免费看高清| 久久久久久久精品精品| 天堂8中文在线网| 久久99热这里只频精品6学生| 热99国产精品久久久久久7| 色视频www国产| 成人美女网站在线观看视频| 精品少妇黑人巨大在线播放| 亚洲成色77777| 亚洲av欧美aⅴ国产| av在线观看视频网站免费| 国内揄拍国产精品人妻在线| 久久人人爽人人爽人人片va| 国产精品女同一区二区软件| 国产精品嫩草影院av在线观看| 男人狂女人下面高潮的视频| 赤兔流量卡办理| 精品酒店卫生间| 十八禁网站网址无遮挡 | 亚洲精品日韩av片在线观看| 久久热精品热| 国产 一区精品| 人人妻人人澡人人看| 日韩强制内射视频| 99九九在线精品视频 | 亚洲欧美成人精品一区二区| 亚洲中文av在线| 亚洲精品自拍成人| 色吧在线观看| 日本与韩国留学比较| 国产国拍精品亚洲av在线观看| 91精品国产国语对白视频| 一级黄片播放器| 久久精品国产a三级三级三级| 99热这里只有精品一区| 一个人免费看片子| kizo精华| 熟女电影av网| 只有这里有精品99| 国产精品人妻久久久影院| av卡一久久| 久久人人爽人人片av| 99热这里只有是精品50| 成年人免费黄色播放视频 | 两个人的视频大全免费| 亚洲av中文av极速乱| 日本wwww免费看| 免费av中文字幕在线| 男男h啪啪无遮挡| 国产视频内射| a级一级毛片免费在线观看| 亚洲成色77777| 亚洲精品视频女| 国产成人91sexporn| 国产综合精华液| 亚洲国产色片| 国产免费又黄又爽又色| 99精国产麻豆久久婷婷| 亚洲国产精品专区欧美| 精品卡一卡二卡四卡免费| 国产精品国产三级国产专区5o| 人妻 亚洲 视频| 亚洲精品久久久久久婷婷小说| 亚洲精品国产色婷婷电影| 久久久久国产精品人妻一区二区| 国产色婷婷99| 老司机影院毛片| 纵有疾风起免费观看全集完整版| 99久久精品热视频| 啦啦啦视频在线资源免费观看| 97精品久久久久久久久久精品| 精品酒店卫生间| 欧美激情极品国产一区二区三区 | 婷婷色综合大香蕉| av在线播放精品| 香蕉精品网在线| 国产成人午夜福利电影在线观看| 亚洲激情五月婷婷啪啪| a级毛片免费高清观看在线播放| 男人爽女人下面视频在线观看| 成人美女网站在线观看视频| 久久女婷五月综合色啪小说| 亚洲国产av新网站| 国模一区二区三区四区视频| av在线老鸭窝| 久久国产精品大桥未久av | 最近最新中文字幕免费大全7| 久久青草综合色| 日本欧美国产在线视频| 亚洲不卡免费看| 一级,二级,三级黄色视频| 人人妻人人添人人爽欧美一区卜| 极品教师在线视频| 人妻制服诱惑在线中文字幕| 在线观看av片永久免费下载| 丰满乱子伦码专区| 亚洲精品成人av观看孕妇| 欧美一级a爱片免费观看看| 国产免费一级a男人的天堂| 国产一区二区在线观看av| 少妇裸体淫交视频免费看高清| 免费看av在线观看网站| 欧美少妇被猛烈插入视频| 国产视频首页在线观看| 男人爽女人下面视频在线观看| 观看免费一级毛片| 免费大片黄手机在线观看| 国产午夜精品久久久久久一区二区三区| 丰满乱子伦码专区| av女优亚洲男人天堂| 成人午夜精彩视频在线观看| 国产熟女欧美一区二区| 一区二区av电影网| 久久久久久久精品精品| 中文字幕免费在线视频6| 99国产精品免费福利视频| 亚洲精品视频女| 国产一区亚洲一区在线观看| 亚洲国产精品成人久久小说| 久久久国产欧美日韩av| 久久久久久久久久人人人人人人| 欧美bdsm另类| 91久久精品国产一区二区成人| 久久久久久久久久久丰满| 欧美人与善性xxx| 久久久午夜欧美精品| 男女边吃奶边做爰视频| 麻豆成人午夜福利视频| 久久精品国产a三级三级三级| 人妻制服诱惑在线中文字幕| 国产真实伦视频高清在线观看| 久久 成人 亚洲| 乱系列少妇在线播放| 免费人成在线观看视频色| 青春草亚洲视频在线观看| 成人无遮挡网站| 熟妇人妻不卡中文字幕| 国产探花极品一区二区| 亚洲欧美精品专区久久| 成人漫画全彩无遮挡| 日韩 亚洲 欧美在线| 成年女人在线观看亚洲视频| 最近中文字幕2019免费版| 日韩不卡一区二区三区视频在线| 国产精品三级大全| 99久国产av精品国产电影| 欧美精品高潮呻吟av久久| 全区人妻精品视频| 日韩大片免费观看网站| 一级毛片电影观看| 老司机亚洲免费影院| 免费少妇av软件| 精品熟女少妇av免费看| 亚洲精品日韩在线中文字幕| 中国国产av一级| 少妇熟女欧美另类| 国产精品人妻久久久久久| 欧美精品人与动牲交sv欧美| 少妇的逼好多水| 久久久久网色| 国产色婷婷99| 亚洲人与动物交配视频| a级一级毛片免费在线观看| 欧美bdsm另类| 国产精品熟女久久久久浪| 欧美97在线视频| 亚洲久久久国产精品| 午夜av观看不卡| 国产男女内射视频| 最近手机中文字幕大全| 免费少妇av软件| 22中文网久久字幕| 男女国产视频网站| 久久精品国产亚洲网站| 日韩成人伦理影院| 日本av免费视频播放| 亚洲色图综合在线观看| 国产91av在线免费观看| 伊人久久精品亚洲午夜| 精华霜和精华液先用哪个| 免费观看无遮挡的男女| 在线亚洲精品国产二区图片欧美 | av卡一久久| 极品少妇高潮喷水抽搐| 九草在线视频观看| 3wmmmm亚洲av在线观看| 18禁在线无遮挡免费观看视频| 精品亚洲成国产av| 亚洲第一区二区三区不卡| 亚洲欧洲精品一区二区精品久久久 | 熟女人妻精品中文字幕| av视频免费观看在线观看| 一区二区三区精品91| 日韩成人av中文字幕在线观看| 蜜桃久久精品国产亚洲av| 欧美区成人在线视频| 国产精品免费大片| 亚洲国产欧美日韩在线播放 | 免费观看性生交大片5| 久久精品国产鲁丝片午夜精品| 亚洲av中文av极速乱| 性色avwww在线观看| 亚洲在久久综合| 成人特级av手机在线观看| videos熟女内射| 观看av在线不卡| 亚洲精品成人av观看孕妇| 国产成人精品久久久久久| 国产成人免费无遮挡视频| 久久亚洲国产成人精品v| xxx大片免费视频| 久久久久久久国产电影| 国产黄色免费在线视频| 国产在视频线精品| 免费人成在线观看视频色| 日本与韩国留学比较| 欧美国产精品一级二级三级 | 午夜免费观看性视频| 日韩伦理黄色片| 麻豆成人午夜福利视频| 久久鲁丝午夜福利片| 我的女老师完整版在线观看| 水蜜桃什么品种好| 99国产精品免费福利视频| 人人妻人人澡人人看| 男人爽女人下面视频在线观看| 国产伦理片在线播放av一区| 久久精品国产自在天天线| 最近中文字幕高清免费大全6| 亚洲国产毛片av蜜桃av| 国产黄片美女视频| 成人国产av品久久久| 日韩一区二区三区影片| 色5月婷婷丁香| 国产精品福利在线免费观看| 美女内射精品一级片tv| 日韩欧美一区视频在线观看 | 久久99蜜桃精品久久| 午夜福利视频精品| 青春草视频在线免费观看| 久久99蜜桃精品久久| 日韩免费高清中文字幕av| 午夜视频国产福利| 五月伊人婷婷丁香| 欧美精品亚洲一区二区| 亚洲欧美日韩东京热| 少妇人妻一区二区三区视频| 精品午夜福利在线看| 妹子高潮喷水视频| 黄色配什么色好看| 亚洲自偷自拍三级| 妹子高潮喷水视频| 人体艺术视频欧美日本| 久久ye,这里只有精品| 十八禁高潮呻吟视频 | 久久久午夜欧美精品| 中文字幕av电影在线播放| 欧美日韩亚洲高清精品| 色网站视频免费| av又黄又爽大尺度在线免费看| 青春草国产在线视频| 欧美+日韩+精品| 午夜激情久久久久久久| 亚洲av综合色区一区| 永久免费av网站大全| 午夜福利视频精品| 国产日韩欧美在线精品| 国产极品粉嫩免费观看在线 | 精品亚洲成a人片在线观看| 2022亚洲国产成人精品| 一本久久精品| 亚洲人成网站在线观看播放| 久久99精品国语久久久| 麻豆成人午夜福利视频| 成人毛片a级毛片在线播放| 色婷婷久久久亚洲欧美| 午夜福利在线观看免费完整高清在| 插阴视频在线观看视频| 国产有黄有色有爽视频| 大香蕉97超碰在线| 寂寞人妻少妇视频99o| 老司机亚洲免费影院| 一级爰片在线观看| 在线 av 中文字幕| 免费大片黄手机在线观看| 国国产精品蜜臀av免费| 亚洲国产欧美日韩在线播放 | 亚洲欧美成人综合另类久久久| 天美传媒精品一区二区| 久久婷婷青草| 免费播放大片免费观看视频在线观看| 欧美日韩av久久| 欧美日本中文国产一区发布| 亚洲自偷自拍三级| 一级av片app| 国产免费一级a男人的天堂| 国产精品女同一区二区软件| 国产免费又黄又爽又色| 午夜激情久久久久久久| av专区在线播放| 久久久国产一区二区| 免费看不卡的av| av天堂中文字幕网| 亚洲综合精品二区| 色网站视频免费| 国产片特级美女逼逼视频| 亚洲精品视频女| 性色av一级|