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

    一種基于TDLAS的高分辨率二維溫度場重建算法及數(shù)值仿真

    2014-07-10 13:13:26李金義朱可柯杜振輝姚宏寶
    實驗流體力學 2014年3期
    關(guān)鍵詞:平均偏差光路溫度場

    李金義,朱可柯,杜振輝,周 濤,姚宏寶

    (1.天津大學精密測試技術(shù)及儀器國家重點實驗室,天津 300072;2.天津津航技術(shù)物理研究所,天津 300308)

    一種基于TDLAS的高分辨率二維溫度場重建算法及數(shù)值仿真

    李金義1,朱可柯1,杜振輝1,周 濤2,姚宏寶2

    (1.天津大學精密測試技術(shù)及儀器國家重點實驗室,天津 300072;2.天津津航技術(shù)物理研究所,天津 300308)

    發(fā)展了一種基于TDLAS正交光路的二維燃燒場溫度重建算法,為了評價該重建方法,數(shù)值仿真了2N(N≥3)條視線(Light-of-sight,LOS)測量路徑呈N×N均勻排布時對兩種不同的溫度場的二維重建結(jié)果,并定義了最大偏差和平均偏差兩個物理量來描述N值對重建結(jié)果的影響。結(jié)果表明,對于對稱的單峰非均勻溫度場,重建的溫度場最大偏差在50K(2.5%)以內(nèi),最大偏差隨N的奇偶而波動,其總體趨勢隨N增大而減小;平均偏差隨著N的增大而減小,N≥9時最大偏差和平均偏差的改善均已不明顯;對于雙峰溫度場,用該方法二維重建的結(jié)果最大偏差超過350K(15.2%),平均偏差大于0.03,并出現(xiàn)嚴重失真,此時可通過額外布置斜穿待測場的測量路徑來改善測量結(jié)果。研究結(jié)果對實際的二維測量系統(tǒng)的搭建和應用有指導意義。

    二維溫度場重建;調(diào)諧激光吸收光譜;正交光路;數(shù)值仿真

    0 引 言

    可調(diào)諧二極管激光吸收光譜(TDLAS)技術(shù)用于燃燒流場的溫度診斷起源于上世紀70年代末[1-3],相比于傳統(tǒng)的探針式接觸測溫法,TDLAS不干擾待測流場,具有高靈敏度、高的譜分辨率、快速的響應時間、大的動態(tài)測量范圍、多參量同時測量等優(yōu)點。

    傳統(tǒng)的TDLAS技術(shù)是對光路上的吸收氣體進行一維路徑積分的視線(Line of Sight,LOS)測量[4-7]。測量的量是該條路徑的積分吸光度,但這種測量方法對測量路徑上的溫度等參量的分布無能為力。將TDLAS技術(shù)與計算機斷層成像技術(shù)(CT)結(jié)合起來(稱為Tunable Diode Laser Absorption Tomography,TDLAT技術(shù)),可以達到獲取燃燒流場的空間二維(2D)分布的目的。TDLAT屬于透射斷層成像,從光源發(fā)出激光在通過待測燃燒流場時被目標氣體吸收,透射的部分被探測器接收,得到氣體的吸收光譜。由于流場的非均勻特性,其各部分對激光的吸收不同,探測器獲得的氣體吸收光譜實際上反映了流場各部分對激光的吸收情況,然后利用各種反演算法從測得的激光吸收光譜數(shù)據(jù)重建出待測流場內(nèi)部某種物理量(如濃度、溫度、流速等)的分布圖。

    TDLAT最早出現(xiàn)在20世紀90年代中期[8-9],近年來愈來愈得到國內(nèi)外研究者的重視,發(fā)展很快。從TDLAT成像光路系統(tǒng)構(gòu)成角度而言,TDLAT可分為掃描式和固定光路式。掃描式一般是將光源和探測器固定在旋轉(zhuǎn)臺上,通過機械轉(zhuǎn)動使得測量光路掃描整個待測場,然后利用CT成像技術(shù)進行二維重建。弗吉尼亞大學的Martin等人[10-11]使用二維計算機層析掃描技術(shù)對超燃沖壓發(fā)動機氣體管進行水蒸汽濃度測量,激光發(fā)射與接收在同一端,另一端安裝反射鏡,二者在同一圓形軌道上運動,掃描得到的氣體濃度二維幾何分布數(shù)據(jù)用MATLAB進行重建分析。另外,該研究小組與NASA合作,為了實現(xiàn)對超聲速燃燒風洞流場特性的測量,對TDLAT進行了方法和實驗研究,采用光源和探測器旋轉(zhuǎn)+扇形掃描的方式來獲取投影數(shù)據(jù),每15°進行旋轉(zhuǎn)步進,每個位置點收集45條扇形投影,共用時29min[12]。Kasyutich等人[13]采用波長1391.67nm的可調(diào)諧激光器測通過將兩運動平臺安裝在直徑97cm的圓形軌道的一條直徑兩端,其中一個運動平臺上裝有激光器和反光鏡,另一平臺上裝有光接收器和反光鏡,每個運動平臺自身也可轉(zhuǎn)動11°的轉(zhuǎn)角,這些部件的運動均由程序控制,由4個伺服電機完成,通過掃描獲取5個角度共55條光線對水的濃度和溫度二維分布重建。王飛等人[14]利用4個不同位置的激光發(fā)射端旋轉(zhuǎn)的方式實現(xiàn)每個測量點100條投影的扇形掃描。李飛等人[15]設(shè)計了一套基于6平行光束—旋轉(zhuǎn)測量的TDLAT系統(tǒng),使用代數(shù)重建算法(ART)對超燃直連臺燃燒室出口的溫度和濃度分布進行了重建。宋俊玲等人[16-17]將燃燒爐固定在步進電機上,通過移動和旋轉(zhuǎn)燃燒爐的位置獲取多個投影,并通過在相鄰光線間設(shè)置虛擬光線,從而提高重建圖像的精確度。掃描法的優(yōu)點是通過很密集的掃描可以得到大量的沿各個角度的投影,從而可實現(xiàn)較高的空間分辨率和重建精確度,但由于使用了機械旋轉(zhuǎn)結(jié)構(gòu),系統(tǒng)復雜并且時間分辨率低,不能滿足瞬態(tài)燃燒過程測量的需求。

    固定光路式通過布置多條LOS正交測量光路將待測區(qū)域分成許多網(wǎng)格,在測量過程中光路和被測燃燒器都不發(fā)生運動,然后利用多路測量的結(jié)果對待測場進行二維重建。美國空軍實驗室的Givens等人[18]采用的激光光源有1391、1392和1397nm 3種波長,以適應水分子在不同狀態(tài)時吸收線的不同,共有16條固定的吸收光路,交叉分布,應用計算流體動力學(CFD)模擬與TDLAT結(jié)合的方法對測量數(shù)據(jù)進行計算分析,重建的二維溫度結(jié)果與熱電偶所測結(jié)果相比較,誤差在5%以內(nèi)。姜治深等人[19]利用24路光束從4個不同的角度穿過待測燃燒區(qū)域,利用ART算法重建了燃燒場中H2O的溫度和濃度二維分布。闞瑞峰等人[20]利用4×4的光路布置對發(fā)動機出口的溫度進行了二維重建。由于固定光路法得到的投影數(shù)較少,因此空間分辨率和重建精度不如掃描法,但其系統(tǒng)簡單,可實現(xiàn)燃燒流場的高時間分辨率的二維測量,故在實際的燃燒流場測量中更具實用性,其中又以正交光路的使用最為普遍[21-23]。對于正交光路二維重建方法,據(jù)我們所知,目前并沒有對LOS測量路徑數(shù)目對重建結(jié)果影響分析的研究,在實際測量時選擇光路數(shù)目仍比較盲目,因此,對正交光路重建方法進行更深入的研究評價是非常必要的。

    本文發(fā)展了一種正交光路重建算法,數(shù)值仿真研究了不同視線測量路徑數(shù)目分別對單峰分布和雙峰分布溫度場的2D重建結(jié)果,并給出定量評價,為實際的2D測量系統(tǒng)的設(shè)計和應用奠定基礎(chǔ)。

    1 二維溫度重建方法

    1.1 雙線測溫法

    基于LOS吸收測量的2D溫度場重建時,通過對激光器的輻射波長進行調(diào)諧,從而覆蓋所選的一對氣體吸收譜線,然后根據(jù)測得的積分吸光度反演氣體的溫度、濃度等信息。許多文獻已對雙線測溫法的原理進行了詳細介紹[5-7],簡言之,由Beer-Lambert定律,每條吸收譜線的積分吸光度A1和A2(cm-1)可通過一個可調(diào)諧、窄線寬的激光器測得,如(1)式所示:

    其中,I0和I分別表示入射和透射的激光光強,積分對整個頻率v(cm-1)域進行。而每個積分吸光度還可以表示成:

    其中,T(K)和x分別表示氣體的溫度和摩爾濃度,總壓強p(atm)和光程長L(cm)通常事先測得。分子吸收線強Si(cm-2atm-1)由下式給出:

    其中,h(J·s)為普朗克常量,c(cm/s)為光速,k(J/K)為波爾茲曼常數(shù),T0(K)是參考溫度,Q為吸收氣體的配分函數(shù),v0,i(cm-1)和分別為第i個躍遷譜線的中心頻率和低能級能量。

    因此,吸收譜線的面積比A1/A2可以化簡為只包含T的函數(shù):

    利用測得的譜線A1和A2,通過求解方程(4)可以得到溫度T,然后利用(2)求得摩爾濃度x。

    1.2 基于正交光路2D溫度重建

    對于一般的2D分布測量,可以將非均勻特性的待測區(qū)域劃分成M×N的網(wǎng)格,并假設(shè)每個網(wǎng)格內(nèi)的分布是均勻的,每個網(wǎng)格可稱為一個像素,最簡單的情況是等距均勻劃分網(wǎng)格,如圖1中虛線所示,將待測區(qū)域均勻劃分成M×N的網(wǎng)格,并由M×N條正交排布的LOS測量路徑對感興趣的場進行測量,如圖1中實線所示。

    當吸收氣體的特性(溫度或/和濃度)沿測量路徑為非均勻時,方程(2)可修改為:

    這里,假設(shè)壓強p為均勻的,l為沿LOS的坐標積分元,ai(l)是第i個吸收譜線的吸收系數(shù),a和b待測區(qū)域的積分限。在2D成像測量中,x和T是空間坐標X和Y的函數(shù),如圖1所示。

    圖1 正交網(wǎng)格法2D測量示意圖Fig.1 Schematic of the orthogonal grid for 2D measurement

    在利用圖1所示的LOS路徑排布進行2D圖像重建時,我們基于如下兩個假設(shè):(1)LOS光束的寬度比像素的寬度小得多,因此每條光束不能同時穿過多個像素;(2)每個像素中溫度和濃度的分布都是均勻的。因此,第m條橫向測量路徑與第n條縱向測量路徑的交點的溫度值Tmn和濃度值xmn即代表整個[m,n]像素的溫度和濃度。這樣,第i個躍遷譜線在[m,n]像素的吸收吸收可用表征。

    因此,(5)式的積分可以寫成離散求和的形式:

    值得注意的是式(7)為統(tǒng)一的表達式,不受LOS測量路徑排布和網(wǎng)格劃分的限制。其中,F(xiàn)為向量ai的權(quán)重系數(shù)矩陣。其元素fmn(m=1,2…M;n=1,2…N)為

    利用(9)式為目標函數(shù)對方程(7)進行求解時,通常使用基于Kacmarz迭代的代數(shù)重建算法(Algebra Reconstruction Technique,ART),但在本方案中,由于采用正交交叉測量路徑,路徑數(shù)較少,方程(7)為病態(tài)性很嚴重的不適定反問題,并且在實際投影獲取過程中,由于燃燒流場的惡劣環(huán)境等因素,測量路徑積分吸光度會存在較大程度的噪聲,因此,采用Tikhonov正則化算法[24]對其進行求解,可以克服上述兩個問題,得到每個交叉點(即每個像素)處的吸收系數(shù)。

    圖2 L曲線法求最優(yōu)正則化參數(shù)示意圖Fig.2 Schematic of L-curve method for optimal regularization parameter

    使用Tikhonov正則化方法時,(7)和(9)式的求解變?yōu)榍蠼庖韵聠栴}:

    其中aη為Tikhonov正則化解,a*為a的初始估計值,在這里可設(shè)置為0。η為正則化參數(shù),其作用是平衡正則化解的殘差范數(shù)與正則化解的范數(shù)在優(yōu)化時的權(quán)重。正則化參數(shù)η的選取異常重要,目前最廣泛采用的是L曲線法[24]。首先以正則化參數(shù)η為變量,做正則化解的范數(shù)與殘差范數(shù)關(guān)系曲線,并且橫縱坐標都使用對數(shù)形式,即得到L曲線,如圖2所示,然后尋找該曲線上曲率最大的點,對應的參數(shù)即最優(yōu)正則化參數(shù)。

    將正則化解aη即像素內(nèi)的吸收系數(shù)代入式(10)即可得到每個像素的溫度值Tmn。然后,通過(2)式可以求得該像素內(nèi)的摩爾濃度xmn。因此,就獲得了具有M×N像素分辨率的2D溫度和濃度場。

    更進一步,假設(shè)溫度和濃度平滑分布,即在整個的重建2D平面內(nèi)其1階和2階導數(shù)連續(xù)?;贛×N個網(wǎng)格的重建值,利用自然雙三次樣條插值的方法可以得到網(wǎng)格交叉點之間的溫度和濃度值,從而獲得更高分辨率的2D圖像。

    2 數(shù)值仿真

    在TLDAS技術(shù)用于燃燒診斷時,最常用的檢測氣體是水蒸汽。因為水蒸汽是燃燒的產(chǎn)物,濃度較高,并且在近紅外波段有著很強的吸收,易于測量。我們使用H2O氣體位于1391.67和1468.89nm的兩條吸收譜線,其線強隨溫度的變化關(guān)系可由Hitran譜庫獲得,如圖3(a)所示;由此可推導出線強比及其對溫度檢測的靈敏度隨溫度的變化關(guān)系,如圖3(b)所示;線強比和檢測靈敏度都隨溫度的升高而降低,但在所感興趣的溫度范圍(1000~2600K)內(nèi)都有較高的靈敏度(>1)。在實際的TDLAS燃燒診斷系統(tǒng)中,這兩條譜線可由兩個可調(diào)諧激光器的時分復用(TDM)或波分復用(WDM)獲得。

    圖3 所選譜線線強隨溫度的變化及線強比和靈敏度隨溫度的變化Fig.3 Line strengthen of selected line pair versus temperature and integration absorption ratio of selected line pairs and its sensitivity versus temperature

    為了計算方便并不失一般性,假設(shè)待測場為15cm ×15cm的正方形2D平面,在此區(qū)域內(nèi)H2O氣體摩爾濃度均勻分布(x=15%)。并假設(shè)在實際燃燒器中可能遇到的兩種溫度場的分布模型[22]如圖4所示,模型都用1501×1501個像素點描述,圖4(a)為對稱的單峰溫度分布情況,為拋物面中心疊加高斯分布;圖4(b)為雙峰溫度情況,為拋物面的不同位置疊加不同的高斯分布。對假設(shè)的方形溫度場進行2D重建時,使用2N條LOS測量路徑,按圖1所示的方式進行等間距的正交排布,并穿過N×N網(wǎng)格組成的測量區(qū)域。

    2.1 對稱的單峰非均勻溫度場

    根據(jù)上述2D重建算法,首先利用Tikhonov正則化方法得到兩條譜線分別在N×N個像素內(nèi)的吸收系數(shù),然后利用(10)式得到每個像素內(nèi)的溫度,最后利用自然雙三次樣條插值對N×N個交叉點處的溫度進行平滑,從而得到與模型同一分辨率(1501× 1501)的2D分布。對圖4(a)所示的對稱單峰非均勻溫度場模型進行重建。不同的LOS測量路徑數(shù)N(N≥3)情況下的重建結(jié)果及偏差如圖5所示。需要說明的是,由于從最邊緣的測量路徑到待測區(qū)域邊界處的溫度值獲得是依賴外延法獲取而非內(nèi)插值法,因此,誤差相對較大[21],為了更好地評價重建結(jié)果,只對最邊緣的測量路徑所包圍的測量區(qū)域的重建結(jié)果及誤差進行分析。

    由圖5可見,該算法對單峰溫度場的重建效果較好,最大絕對偏差基本都小于50K(N=4時例外,為57K),最大相對偏差小于2.5%。重建的絕對偏差存在相似的規(guī)律,即峰值和4個角位置處重建結(jié)果偏小,而靠中間的部分重建結(jié)果偏大,并且負偏差遠大于正偏差,這是由于算法本身在計算時的峰值位置處的吸光度比偏大造成的。

    圖4 待重建的溫度場模型Fig.4 The model of reconstruction temperature field

    圖5 不同N值時單峰溫度場2D重建的絕對偏差分布圖Fig.5 2D reconstructionresults and absolute deviations for single peak temperature field with differentN

    2.2 雙峰溫度場

    類似地,對圖4(b)所示的雙峰溫度場模型進行重建。不同的LOS測量路徑數(shù)N(N≥3)情況下的重建結(jié)果及偏差如圖6所示。

    圖6 不同N值時雙峰溫度場2D重建溫度及絕對偏差分布圖Fig.6 2D reconstruction results and absolute deviations for the bimodal peak temperature field temperature field with differentN

    由圖6可知,對于雙峰溫度場重建,該算法效果較差,最大絕對偏差都超過了350K(對應相對偏差為15.2%),N=3和4時,不能分辨雙峰,N>4時,重建結(jié)果雖然出現(xiàn)了兩個峰,但峰的位置在(3,6)和(6,3)附近,與原始溫度場模型中的峰值坐標(4,9.5)和(10,7)不相符。

    3 討 論

    為了更好地評價上述2D重建結(jié)果,定義最大偏差eM和平均偏差eT兩個物理量,分別描述重建溫度場與模型對應像素點偏差的最大值和綜合平均值。

    兩種待測場模型重建的偏差隨N的變化曲線如圖7所示。兩種情況下,重建的最大偏差和平均偏差與N的關(guān)系曲線呈“L”型,隨著N的增大,eT和eM都是先迅速減小然后再緩慢變小,拐點都在N=5處。因此,綜合考慮2D網(wǎng)格成像系統(tǒng)復雜性和成本及誤差因素,選擇5×5的路徑排布是最佳方案。另外,由圖7(a)可以看出,eM隨著N的奇偶而波動,這是由于N為奇數(shù)時,中間測量路徑經(jīng)過待測溫度場峰值點,因此重建結(jié)果的最大偏差稍小一些,而N為偶數(shù)時,測量路徑不能經(jīng)過峰值點,偏差稍大。N≥9時最大偏差和平均偏差的改善不明顯。由圖7(a)可以看出,N>6時,最大偏差eM不再像圖7(a)中隨N的增大而減小,而是基本保持不變,而平均偏差eT雖隨著N增大有減小趨勢,但已非常平緩,甚至N=100時eT和eM都未能有較大改善。其原因在于布置的LOS測量路徑?jīng)]有通過溫度峰值處,導致溫度場的峰值信息未能被捕捉到。

    圖8 對圖4(b)中雙峰溫度場模型進行重建時的LOS光路布置示意圖Fig.8 Diagrammatic arrangement of LOSoptical path for reconstruction of bimodal temperature phantom shown as Figure 4(b)

    表1 新增LOS測量光路前后的2D溫度重建結(jié)果比較Table 1 Comparison of reconstruction results without and with added LOS measurement optical path

    因此,為了改善雙峰溫度場的2D重建結(jié)果,對N=5時的正交光路排布進行改進,新增兩條LOS測量路徑,按對角線交叉排布,如圖8中的點畫線所示,依然利用上述算法,重建的結(jié)果如圖9(a)所示??梢姡黾訙y量光路后的重建結(jié)果有了較大改善,具體如表1所示,eM和eT的值比較可見,重建結(jié)果比N=5,6,7時都有了改善。特別地,若再增加一條穿過兩個峰值點的路徑,重建結(jié)果改善更為明顯,如圖9(b)所示,eM和eT相比與之前N=5時分別改善了34.5%和25.9%。

    圖9 對N=5情形下新增斜向測量光路對雙峰場的重建結(jié)果和偏差Fig.9 Reconstruction results and deviations for bimodal peak phantom with additional oblique measurement path atN=5

    4 結(jié) 論

    發(fā)展了一種基于TDLAS正交光路的二維燃燒場溫度重建算法,并利用數(shù)值仿真的方法評價了LOS測量路徑數(shù)N值對重建結(jié)果的影響。結(jié)果表明,對于對稱的單峰非均勻溫度場,重建的溫度場最大偏差在50K以內(nèi),相對偏差在2.5%以內(nèi),最大偏差隨N的奇偶而波動,其總體趨勢隨N增大而減??;平均偏差隨著N的增大而減小,N≥9時最大偏差和平均偏差的改善不明顯;對于雙峰溫度場,用該方法二維重建的結(jié)果最大偏差超過350K,最大相對偏差大于10%,平均偏差大于0.03,并出現(xiàn)嚴重失真。因此,該方法適合待測對象為單峰分布的燃燒場,并且綜合考慮2D網(wǎng)格成像系統(tǒng)復雜性、成本及誤差因素,選擇5×5的路徑排布是最佳方案。對雙峰分布的燃燒場,該重建方法誤差較大,但可通過增加斜穿過待測場的路徑提高重建精度。研究結(jié)果對實際的二維測量系統(tǒng)的搭建和應用有指導意義。

    [1]Hanson R K,Kuntz P A,Kruger C H.High-resolution spectroscopy of combustion gases using a tunable IR diode laser[J].Applied Optics,1977,16:2045-2048.

    [2]Lück K C,Thielen W.Measurements of temperatures and OH-concentrations in a lean methane-air flame using high-resolution laser-absorption spectroscopy[J].Journal of Quantitative Spectroscopy and Radiative Transfer,1978,20(1):71-79.

    [3]Hanson R K,F(xiàn)alcone P K.Temperature measurement technique for high-temperature gases using a tunable diode laser[J].Applied Optics,1978,17(16):2477-2480.

    [4]Arroyo M P,Hanson R K.Absorption measurements of watervapor concentration,temperature,and line-shape parameters using a tunable In Ga AsP diode laser[J].Applied Optics,1993,32(30):6104-6116.

    [5]Nagali V,Chou S I,Baer D S,et al.Diode-laser measurements of temperature-dependent half-widths of H2O transitions in the 1.4μm region[J].Journal of Quantitative Spectroscopy and Radiative Transfer,1997,57(6):795-809.

    [6]李飛,余西龍,陳立紅,等.TDLAS測量甲烷/空氣預混平面火焰溫度和H2O濃度[J].實驗流體力學,2009,23(2):40-44.Li Fei,Yu Xilong,Chen Lihong,et al.Temperature and water vapour concentration measurements of CH4/air premixed flat flame based on TDLAS[J].Journal of Experiments in Fluid Mechanics,2009,23(2):40-44.

    [7]許振宇,劉文清,劉建國,等.基于可調(diào)諧半導體激光器吸收光譜的溫度測量方法研究[J].物理學報,2012,61(23):234204.Xu Zhenyu,Liu Wenqing,Liu Jianguo,et al.Temperature measuremtns based on tunable diode laser absorption spectroscopy[J].Acta Phys Sin,2012,61(23):234204.

    [8]Kauranen P,Hertz H M,Svanberg S.Tomographic imaging of fluid flows by the use of two-tone frequency-modulation spectroscopy[J].Optics Letters,1994,19(18):1489-1491.

    [9]Kessler W J,Allen M G,Lo E Y M,et al.Tomographic reconstruction of air temperature and density profiles using tunable diode laser absorption measurements on O2[C].26th AIAA Plasmadynamics and Lasers Conference.San Diego,USA,1995:1-11.

    [10]Martin E F,Goyne C P,Diskin G S.Analysis of a tomography technique for a scramjet wind tunnel[J].International Journal of Hypersionics,2010,1(3):173-180.

    [11]Bryner E,Sharma M G,Goyne C P,et al.Tunable diode laser absorption technique development for determination of spatially resolved water concentration and temperature[C].48th AIAA Aerospace Sciences Meeting,2010.

    [12]Martin E F,Goyne C P,Diskin G S.Analysis of a tomography technique for a scramjet wind tunnel[J].International Journal of Hypersonics,2010,1(3):173-180

    [13]Kasyutich V L,Martin P A.Towards a two-dimensional concentration and temperature laser absorption tomography sensor system[J].Applied Physics B,2011,102(1):149-162

    [14]Wang F,Cen K F,Li N,et al.Two-dimensional tomography for gas concentration and temperature distributions based on tunable diode laser absorption spectroscopy[J].Measurement Science and Technology,2010,21(4):045301.

    [15]李飛,余西龍,林鑫,等.基于TDLAS的層析成像技術(shù)[C].高超聲速專題研討會暨第五屆全國高超聲速科學技術(shù)會議,2012.Li Fei,Yu Xilong,Lin Xin,et al.Computer tomography technique based on TDLAS[C].Hypersonic Symposium and 5th National Hypersonic Science and Technology Conference,2012.

    [16]宋俊玲,洪延姬,王廣宇,等.基于激光吸收光譜技術(shù)的燃燒場氣體溫度和濃度二維分布重建研究[J].物理學報,2012,61(24):240702.Song Junling,Hong Yanji,Wang Guangyu,et al.Two-dimensional resconstructions of gas temperature and concentration in combustion based on tunable diode laser absorption spectroscopy[J].Acta Phys Sin,2012,61(24):240702.

    [17]Song J,Hong Y,Wang G,et al.Algebraic tomographic reconstruction of two-dimensional gas temperature based on tunable diode laser absorption spectroscopy[J].Applied Physics B,2013:1-9.

    [18]Givens R N.Demonstration of a strategy to perform two-dimensional diode laser tomography[D].Air Force Institute of Technology Air University,2008

    [19]姜治深,王飛,邢大偉,等.可調(diào)諧半導體激光吸收光譜技術(shù)應用于平面火焰中氣體濃度二維分布重建的研究[J].光譜學與光譜分析,2012,32(11):2891-2896.Jiang Zhishen,Wang Fei,Xing Dawei,et al.The reconstruction of two-dimensional distributions of gas concentration in the flat flame based on tunable laser absorption spectroscopy[J].Spectroscopy and Spectral Analysis,2012,32(11):2891-2896.

    [20]闞瑞峰.TDLAS在大氣監(jiān)測及燃燒診斷中的應用方法研究[C].光譜技術(shù)在燃燒診斷及大氣檢測中的應用學術(shù)研討會,2013,4.Kan Ruifeng.Application and method research of TDLAS on atmospheric monitoring and combustion diagnosis[C].Spectroscopy Symposium on Combustion Diagnosis and Atmospheric Detection,2013.

    [21]Li J,Du Z,Zhou T,et al.Numerical investigation of two-dimensional imaging for temperature and species concentration using tunable diode laser absorption spectroscopy[C].2012 International Conference on Measurement,Information and Control(MIC),IEEE,2012,1:234-238.

    [22]Ma L,Cai W.Numerical investigation of hyperspectral tomography for simultaneous temperature and concentration imaging[J].Applied Optics,2008,47(21):3751-3759.

    [23]Deguchi Y,Yasui D,Adachi A.Development of 2D temperature and concentration measurement method using tunable diode laser absorption spectroscopy[J].Journal of Mechanics Engineering and Automation,2012,(2):543-549

    [24]Hansen P C,O′Leary D P.The use of the L-curve in the regularization of discrete ill-posed problems[J].SIAM Journal on Scientific Computing,1993,14(6):1487-1503.

    A high-resolution algorithm for 2D temperature reconstruction using TDLAS and numerical simulation

    Li Jinyi1,Zhu Keke1,Du Zhenhui1,Zhou Tao2,Yao Hongbao2
    (1.State Key Lab of Precision Measuring Technology and Instruments,Tianjin University,Tianjin 300072,China;2.Tianjin Jinhang Institute of Technical Physics,Tianjin 300308,China)

    We develop a tunable diode laser absorption tomography algorithm for two-dimensional(2D)temperature reconstruction based on orthogonal optical paths.In order to evaluate the reconstruction method,we analyzed the 2D reconstruction results of two different temperature phantoms by means of numerical simulation ofN×Nuniform grid arrangement consisting of 2N(N≥3)light of sight(LOS)measurement path.The maximum deviation and the mean deviation were defined to describe the effect of the value ofNon the reconstruction results.The results show that for the single symmetrical peak non-uniform temperature phantom,the maximum deviation of the reconstruction temperature is less than 50K(2.5%),the maximum deviation fluctuates withNparity and the overall trend is decrease withNincreasing.The mean deviation is reduced asNincreases,which means the reconstruction is better withNincreasing.But whenN≥9,the improvements are not distinct any more.For the bimodal temperature phantom,the maximum deviations of the results using this 2D reconstruction method are more than 350K(15.2%)and the mean deviations are greater than 0.03.What’s worse,there is a serious distortion.However,the reconstruction results can be improved by added measurement path diagonally across the tested temperature field.The study is helpful for the real 2D measurement system setup and actual application.

    two-dimensional temperature reconstruction;tunable diode laser absorption spectroscopy;orthogonal optical path;numerical simulation

    O433.5;TH811.2

    :A

    1672-9897(2014)03-0063-09doi:10.11729/syltlx20130059

    (編輯:張巧蕓)

    2013-07-05;

    :2013-09-28

    精密測試技術(shù)及儀器國家重點實驗室探索課題(PILT1107)和天津市應用基礎(chǔ)及前沿技術(shù)項目(11JCYBJC04900)

    杜振輝,E-mail:duzhenhui@tju.edu.cn

    LiJY,ZhuKK,DuZH,etal.Ahigh-resolutionalgorithmfor2DtemperaturereconstructionusingTDLASandnumericalsimulation.

    JournalofExperimentsinFluidMechanics,2014,28(3):63-71.李金義,朱可柯,杜振輝,等.一種基于TDLAS的高分辨率二維溫度場重建算法及數(shù)值仿真.實驗流體力學,2014,28(3):63-71.

    李金義(1986-),男,河北吳橋人,博士后。研究方向:調(diào)諧激光吸收光譜及燃燒診斷。通信地址:天津大學精密儀器與光電子學院17樓403室(300072)。E-mail:jinyilee@tju.edu.cn

    猜你喜歡
    平均偏差光路溫度場
    河北地方性震級量規(guī)函數(shù)與方位角校正值研究1
    鋁合金加筋板焊接溫度場和殘余應力數(shù)值模擬
    FY-3C/VIRR西北太平洋區(qū)域海表溫度精度評估?
    基于紋影法的溫度場分布測量方法
    MJS工法與凍結(jié)法結(jié)合加固區(qū)溫度場研究
    建筑科技(2018年6期)2018-08-30 03:41:08
    自制立體光路顯示儀
    通天之光路
    脛前動脈穿刺可行性及心肺流轉(zhuǎn)下脛前動脈與橈動脈壓力監(jiān)測的一致性研究
    X80鋼層流冷卻溫度場的有限元模擬
    數(shù)控機床自動編程技術(shù)中尺寸公差的研究
    久久久久久久久中文| 丁香欧美五月| 黑人欧美特级aaaaaa片| 久久性视频一级片| 啦啦啦免费观看视频1| 黑人操中国人逼视频| 日本撒尿小便嘘嘘汇集6| 久久久水蜜桃国产精品网| 国产成人欧美| 美女扒开内裤让男人捅视频| 人成视频在线观看免费观看| 欧美av亚洲av综合av国产av| 国产成+人综合+亚洲专区| 69精品国产乱码久久久| 欧美日本中文国产一区发布| 身体一侧抽搐| 怎么达到女性高潮| 国产精品秋霞免费鲁丝片| 久久久水蜜桃国产精品网| 久久欧美精品欧美久久欧美| 欧美性长视频在线观看| 校园春色视频在线观看| 久久亚洲真实| 免费高清视频大片| 99国产极品粉嫩在线观看| 欧美乱码精品一区二区三区| 久久伊人香网站| 午夜免费激情av| or卡值多少钱| av视频免费观看在线观看| 每晚都被弄得嗷嗷叫到高潮| 黄色丝袜av网址大全| 精品国产一区二区久久| 狠狠狠狠99中文字幕| 国产又色又爽无遮挡免费看| 国产精品免费视频内射| 精品高清国产在线一区| 国产激情久久老熟女| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲欧美在线一区二区| 成人三级黄色视频| 色哟哟哟哟哟哟| 成人精品一区二区免费| 桃色一区二区三区在线观看| 免费女性裸体啪啪无遮挡网站| avwww免费| 侵犯人妻中文字幕一二三四区| 精品欧美一区二区三区在线| 一二三四社区在线视频社区8| 人人妻人人澡欧美一区二区 | 美国免费a级毛片| 法律面前人人平等表现在哪些方面| 满18在线观看网站| 91麻豆精品激情在线观看国产| 12—13女人毛片做爰片一| 一级毛片高清免费大全| 免费不卡黄色视频| а√天堂www在线а√下载| 精品国产一区二区三区四区第35| 黑人巨大精品欧美一区二区mp4| 看免费av毛片| 色尼玛亚洲综合影院| 最新美女视频免费是黄的| 国产麻豆成人av免费视频| 亚洲七黄色美女视频| 免费在线观看完整版高清| 午夜福利欧美成人| 身体一侧抽搐| 欧美激情久久久久久爽电影 | 亚洲欧美日韩无卡精品| 国产亚洲av高清不卡| 在线观看舔阴道视频| 高清在线国产一区| 亚洲成a人片在线一区二区| 午夜免费激情av| 男女之事视频高清在线观看| 天堂√8在线中文| 亚洲欧美日韩无卡精品| 欧美一区二区精品小视频在线| 久久精品亚洲熟妇少妇任你| 黄色毛片三级朝国网站| 国产精品久久电影中文字幕| 麻豆av在线久日| 大陆偷拍与自拍| 欧美激情 高清一区二区三区| 日本欧美视频一区| 午夜视频精品福利| 一级黄色大片毛片| 一本大道久久a久久精品| 亚洲欧美日韩高清在线视频| 国产日韩一区二区三区精品不卡| 午夜福利成人在线免费观看| 两人在一起打扑克的视频| 给我免费播放毛片高清在线观看| 一夜夜www| 成人免费观看视频高清| 亚洲熟妇中文字幕五十中出| 国产又爽黄色视频| 男女午夜视频在线观看| 午夜福利高清视频| 老司机在亚洲福利影院| 欧美乱色亚洲激情| 脱女人内裤的视频| 黄色丝袜av网址大全| 欧美成人免费av一区二区三区| 精品国产一区二区久久| 18美女黄网站色大片免费观看| 大型av网站在线播放| 亚洲欧美一区二区三区黑人| 免费搜索国产男女视频| 嫩草影视91久久| 国产精品香港三级国产av潘金莲| 黄色视频,在线免费观看| 日韩高清综合在线| 午夜福利视频1000在线观看 | 波多野结衣av一区二区av| 午夜视频精品福利| 一区二区三区激情视频| 久热爱精品视频在线9| 一区二区三区激情视频| 一区二区三区激情视频| 中文亚洲av片在线观看爽| 777久久人妻少妇嫩草av网站| 国产精品国产高清国产av| 亚洲色图av天堂| 国产精品99久久99久久久不卡| 免费无遮挡裸体视频| 手机成人av网站| 女性被躁到高潮视频| 欧美激情久久久久久爽电影 | 国产私拍福利视频在线观看| 91字幕亚洲| 国产一级毛片七仙女欲春2 | 欧美大码av| 最新美女视频免费是黄的| 亚洲欧美日韩高清在线视频| 免费在线观看日本一区| 色播在线永久视频| 成人永久免费在线观看视频| netflix在线观看网站| 国产伦一二天堂av在线观看| 国产1区2区3区精品| 国产成人精品久久二区二区免费| 人妻丰满熟妇av一区二区三区| 精品第一国产精品| 日日夜夜操网爽| 国产激情欧美一区二区| 精品久久久久久久人妻蜜臀av | 人妻丰满熟妇av一区二区三区| 91大片在线观看| svipshipincom国产片| 欧美另类亚洲清纯唯美| 热re99久久国产66热| 久久精品国产综合久久久| 国产人伦9x9x在线观看| 精品欧美一区二区三区在线| 国产亚洲欧美98| 老司机午夜十八禁免费视频| 午夜精品久久久久久毛片777| 成人手机av| 欧美成人性av电影在线观看| 国产av精品麻豆| 亚洲精品粉嫩美女一区| 人人妻人人爽人人添夜夜欢视频| 久久精品国产清高在天天线| 日韩欧美一区视频在线观看| 在线视频色国产色| 免费看a级黄色片| 亚洲成人精品中文字幕电影| 男人舔女人的私密视频| 伊人久久大香线蕉亚洲五| 动漫黄色视频在线观看| 欧美精品亚洲一区二区| 曰老女人黄片| 97人妻天天添夜夜摸| 精品熟女少妇八av免费久了| 午夜福利一区二区在线看| 啦啦啦免费观看视频1| 真人一进一出gif抽搐免费| 免费女性裸体啪啪无遮挡网站| 亚洲av第一区精品v没综合| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美一区二区三区黑人| 亚洲七黄色美女视频| 国产精品久久久久久人妻精品电影| 日本免费a在线| 日日干狠狠操夜夜爽| 免费观看人在逋| 亚洲av熟女| 日本 欧美在线| 麻豆av在线久日| 亚洲人成电影观看| 男人操女人黄网站| 99精品久久久久人妻精品| 十八禁人妻一区二区| 一本久久中文字幕| 成年人黄色毛片网站| 欧美日韩一级在线毛片| 黑人欧美特级aaaaaa片| 国产国语露脸激情在线看| 18禁黄网站禁片午夜丰满| 很黄的视频免费| 欧美另类亚洲清纯唯美| 日韩大尺度精品在线看网址 | 午夜视频精品福利| 男人舔女人的私密视频| 亚洲第一av免费看| 婷婷六月久久综合丁香| 99精品欧美一区二区三区四区| 久久久水蜜桃国产精品网| 色综合站精品国产| 叶爱在线成人免费视频播放| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精华国产精华精| 欧美人与性动交α欧美精品济南到| 国产精品美女特级片免费视频播放器 | 久久久久国产精品人妻aⅴ院| 国产乱人伦免费视频| 悠悠久久av| 涩涩av久久男人的天堂| 久久亚洲精品不卡| 中文字幕人妻丝袜一区二区| 国产精品精品国产色婷婷| 黄片大片在线免费观看| 精品无人区乱码1区二区| 一卡2卡三卡四卡精品乱码亚洲| av天堂久久9| 一区二区三区国产精品乱码| 一级毛片女人18水好多| 中文字幕最新亚洲高清| 午夜免费鲁丝| 日日摸夜夜添夜夜添小说| 国产亚洲欧美精品永久| 欧美成人免费av一区二区三区| 好男人在线观看高清免费视频 | 久久久久久久久久久久大奶| 男人的好看免费观看在线视频 | 亚洲最大成人中文| 久久国产精品影院| 免费在线观看影片大全网站| 国产男靠女视频免费网站| 啦啦啦 在线观看视频| 亚洲性夜色夜夜综合| 在线观看免费午夜福利视频| 国产97色在线日韩免费| 久久国产精品人妻蜜桃| 久久精品国产综合久久久| 亚洲自偷自拍图片 自拍| 亚洲欧美精品综合一区二区三区| 丝袜在线中文字幕| 国产在线观看jvid| 国产精品,欧美在线| 日韩视频一区二区在线观看| 十八禁人妻一区二区| 日本vs欧美在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| 50天的宝宝边吃奶边哭怎么回事| 欧美日本视频| 亚洲熟女毛片儿| 亚洲 国产 在线| 日韩欧美一区二区三区在线观看| 99在线人妻在线中文字幕| 欧美成人性av电影在线观看| av福利片在线| 国内久久婷婷六月综合欲色啪| 男男h啪啪无遮挡| 亚洲 欧美 日韩 在线 免费| 国产成人系列免费观看| 国产97色在线日韩免费| 高清在线国产一区| 精品国产国语对白av| 亚洲精品国产区一区二| 国产精品香港三级国产av潘金莲| 久久中文字幕一级| 日本一区二区免费在线视频| 久久久久精品国产欧美久久久| 欧美乱码精品一区二区三区| 国产高清视频在线播放一区| 日韩精品免费视频一区二区三区| 国产av精品麻豆| 国产欧美日韩综合在线一区二区| 制服人妻中文乱码| 99国产精品99久久久久| 动漫黄色视频在线观看| 男女下面进入的视频免费午夜 | 亚洲国产高清在线一区二区三 | 欧美绝顶高潮抽搐喷水| 亚洲国产高清在线一区二区三 | 欧美性长视频在线观看| 波多野结衣高清无吗| 黄片小视频在线播放| 亚洲五月天丁香| 国产成人欧美在线观看| 在线国产一区二区在线| 12—13女人毛片做爰片一| 国产精品香港三级国产av潘金莲| 亚洲av电影在线进入| 操出白浆在线播放| 人人澡人人妻人| 国语自产精品视频在线第100页| 村上凉子中文字幕在线| 琪琪午夜伦伦电影理论片6080| 欧美成狂野欧美在线观看| 国产高清videossex| 日韩精品中文字幕看吧| av中文乱码字幕在线| av有码第一页| 午夜精品在线福利| 伦理电影免费视频| 一边摸一边抽搐一进一出视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲性夜色夜夜综合| 国产麻豆成人av免费视频| 黄色毛片三级朝国网站| 麻豆久久精品国产亚洲av| 国产91精品成人一区二区三区| 国产精品久久视频播放| 欧美日韩亚洲国产一区二区在线观看| 一本综合久久免费| 成人国产一区最新在线观看| 在线天堂中文资源库| 如日韩欧美国产精品一区二区三区| 多毛熟女@视频| 午夜激情av网站| 久久精品国产亚洲av香蕉五月| 日本免费a在线| 热99re8久久精品国产| 一级毛片精品| 久久香蕉精品热| 欧美老熟妇乱子伦牲交| a级毛片在线看网站| 大码成人一级视频| 精品久久久精品久久久| 999久久久国产精品视频| 他把我摸到了高潮在线观看| 香蕉丝袜av| 国产欧美日韩一区二区精品| 91大片在线观看| 亚洲一区中文字幕在线| 亚洲国产精品sss在线观看| 免费观看人在逋| 久久久国产成人精品二区| 欧美日本视频| 午夜影院日韩av| 一个人免费在线观看的高清视频| 久久香蕉激情| 亚洲欧美日韩高清在线视频| 亚洲三区欧美一区| 婷婷六月久久综合丁香| 最近最新免费中文字幕在线| 欧美中文综合在线视频| 国产精品一区二区精品视频观看| 国产在线观看jvid| 最近最新中文字幕大全免费视频| 国产单亲对白刺激| 久久精品亚洲精品国产色婷小说| 午夜影院日韩av| 丝袜美腿诱惑在线| www国产在线视频色| 国产亚洲精品久久久久久毛片| 成人亚洲精品一区在线观看| 夜夜夜夜夜久久久久| 亚洲国产欧美网| 亚洲精品一区av在线观看| 少妇裸体淫交视频免费看高清 | 亚洲精品国产色婷婷电影| 一级a爱视频在线免费观看| 99re在线观看精品视频| 亚洲精品中文字幕一二三四区| 久99久视频精品免费| 日日夜夜操网爽| 12—13女人毛片做爰片一| 国产精品一区二区三区四区久久 | 999久久久国产精品视频| 国产一区二区三区综合在线观看| 真人一进一出gif抽搐免费| 亚洲avbb在线观看| 97人妻天天添夜夜摸| 久久精品aⅴ一区二区三区四区| 久久亚洲真实| 日本在线视频免费播放| 国产精品免费视频内射| 波多野结衣av一区二区av| 国产伦人伦偷精品视频| 热99re8久久精品国产| 热99re8久久精品国产| 熟女少妇亚洲综合色aaa.| 岛国在线观看网站| 亚洲少妇的诱惑av| 亚洲熟妇熟女久久| 国产国语露脸激情在线看| 精品少妇一区二区三区视频日本电影| 国产精品av久久久久免费| 中文字幕最新亚洲高清| 欧美在线黄色| 国产亚洲av嫩草精品影院| 色综合欧美亚洲国产小说| 咕卡用的链子| 午夜影院日韩av| 日韩视频一区二区在线观看| 久久久久久亚洲精品国产蜜桃av| 国产亚洲精品av在线| 黑人巨大精品欧美一区二区蜜桃| 久久中文字幕人妻熟女| 亚洲成国产人片在线观看| 色婷婷久久久亚洲欧美| 在线观看免费视频网站a站| 丝袜人妻中文字幕| 日本在线视频免费播放| 精品国产超薄肉色丝袜足j| 在线观看一区二区三区| 成年版毛片免费区| 精品人妻在线不人妻| 欧美成人一区二区免费高清观看 | 青草久久国产| 我的亚洲天堂| 亚洲专区中文字幕在线| 少妇裸体淫交视频免费看高清 | 一级黄色大片毛片| 老汉色av国产亚洲站长工具| 高清毛片免费观看视频网站| 日本在线视频免费播放| 国产亚洲精品久久久久久毛片| 啦啦啦 在线观看视频| 99精品久久久久人妻精品| 韩国av一区二区三区四区| 午夜福利一区二区在线看| 国产成人av教育| 在线十欧美十亚洲十日本专区| 欧美老熟妇乱子伦牲交| 国产极品粉嫩免费观看在线| 老司机深夜福利视频在线观看| 午夜福利一区二区在线看| 一级作爱视频免费观看| 亚洲自偷自拍图片 自拍| 丰满的人妻完整版| 亚洲中文日韩欧美视频| 亚洲 国产 在线| 免费观看人在逋| 亚洲天堂国产精品一区在线| 亚洲av成人不卡在线观看播放网| 日本a在线网址| 国产成人av激情在线播放| 欧美日本亚洲视频在线播放| 最近最新免费中文字幕在线| 日日夜夜操网爽| 97人妻精品一区二区三区麻豆 | 国产又色又爽无遮挡免费看| 国产亚洲av嫩草精品影院| 亚洲国产日韩欧美精品在线观看 | 丰满人妻熟妇乱又伦精品不卡| 91字幕亚洲| 欧美在线一区亚洲| 波多野结衣巨乳人妻| 亚洲精品久久成人aⅴ小说| 午夜精品久久久久久毛片777| 不卡av一区二区三区| 黄片小视频在线播放| 成人国语在线视频| 国产一级毛片七仙女欲春2 | 女生性感内裤真人,穿戴方法视频| 精品欧美一区二区三区在线| 亚洲精品国产区一区二| 最近最新中文字幕大全免费视频| 人妻丰满熟妇av一区二区三区| 一级a爱片免费观看的视频| 波多野结衣一区麻豆| 搞女人的毛片| 天堂√8在线中文| 欧美另类亚洲清纯唯美| 97碰自拍视频| a在线观看视频网站| 国产伦一二天堂av在线观看| www.999成人在线观看| 国产亚洲欧美精品永久| 国产欧美日韩一区二区三区在线| 欧美在线一区亚洲| 亚洲中文日韩欧美视频| 午夜免费观看网址| www国产在线视频色| 99久久久亚洲精品蜜臀av| 国产成人精品久久二区二区91| 日韩大尺度精品在线看网址 | av片东京热男人的天堂| 日本黄色视频三级网站网址| 免费女性裸体啪啪无遮挡网站| 日韩大尺度精品在线看网址 | 久久久国产成人精品二区| 在线观看www视频免费| 久久人妻熟女aⅴ| 精品久久久久久久久久免费视频| 亚洲自拍偷在线| 久久亚洲真实| 国产人伦9x9x在线观看| 午夜福利18| 欧美日韩福利视频一区二区| 天天一区二区日本电影三级 | 丰满人妻熟妇乱又伦精品不卡| 日韩欧美国产一区二区入口| 一边摸一边抽搐一进一小说| 国产亚洲精品av在线| 男女床上黄色一级片免费看| 日韩大尺度精品在线看网址 | 一级毛片女人18水好多| 久久天堂一区二区三区四区| 成人18禁高潮啪啪吃奶动态图| 免费少妇av软件| 一级作爱视频免费观看| 久久国产精品影院| 国产精品自产拍在线观看55亚洲| 国产av一区在线观看免费| 精品久久久久久,| 国产成人av教育| 国产亚洲精品久久久久久毛片| 国产伦人伦偷精品视频| 99精品在免费线老司机午夜| 一进一出好大好爽视频| 亚洲国产精品久久男人天堂| 亚洲精品久久成人aⅴ小说| 精品国产一区二区久久| 动漫黄色视频在线观看| 亚洲精华国产精华精| 咕卡用的链子| 99在线人妻在线中文字幕| av天堂久久9| 搡老岳熟女国产| 亚洲欧美激情综合另类| 亚洲成人久久性| 久久久久九九精品影院| 黑人操中国人逼视频| 免费高清视频大片| √禁漫天堂资源中文www| 男女午夜视频在线观看| 日本一区二区免费在线视频| 在线观看www视频免费| 国产精品久久久久久人妻精品电影| 久久久久久久久中文| 亚洲av成人av| 久久人妻av系列| 久久精品国产清高在天天线| 欧美另类亚洲清纯唯美| 男女午夜视频在线观看| 人人妻人人澡人人看| 黄网站色视频无遮挡免费观看| 亚洲一码二码三码区别大吗| 亚洲国产日韩欧美精品在线观看 | 色综合欧美亚洲国产小说| av天堂在线播放| 一进一出抽搐动态| 欧美成人午夜精品| 啪啪无遮挡十八禁网站| 亚洲一区中文字幕在线| 久久人妻熟女aⅴ| 一本综合久久免费| 人人妻,人人澡人人爽秒播| 一级,二级,三级黄色视频| 老司机在亚洲福利影院| 欧美在线一区亚洲| 搞女人的毛片| 极品教师在线免费播放| 99久久99久久久精品蜜桃| 久久久国产精品麻豆| 免费一级毛片在线播放高清视频 | 中文字幕最新亚洲高清| 国产精品久久久av美女十八| 十分钟在线观看高清视频www| 此物有八面人人有两片| 热99re8久久精品国产| 12—13女人毛片做爰片一| 男人操女人黄网站| 在线观看www视频免费| 国产色视频综合| 国产午夜福利久久久久久| 青草久久国产| 侵犯人妻中文字幕一二三四区| 一级,二级,三级黄色视频| 亚洲欧美精品综合一区二区三区| 亚洲黑人精品在线| 色婷婷久久久亚洲欧美| 国产成人精品久久二区二区91| 757午夜福利合集在线观看| 18禁美女被吸乳视频| 中文字幕高清在线视频| 美女高潮到喷水免费观看| 一级毛片精品| 视频区欧美日本亚洲| 日韩大尺度精品在线看网址 | 亚洲国产中文字幕在线视频| 国产又色又爽无遮挡免费看| 免费久久久久久久精品成人欧美视频| 精品一品国产午夜福利视频| 亚洲avbb在线观看| 久久婷婷人人爽人人干人人爱 | 欧美在线黄色| 校园春色视频在线观看| 日本 av在线| 国产精品 欧美亚洲| 国产欧美日韩一区二区三| 国产精品一区二区免费欧美| 免费高清视频大片| 精品久久久久久久人妻蜜臀av | 国产欧美日韩综合在线一区二区| www.自偷自拍.com| av片东京热男人的天堂| 国产在线观看jvid| www.精华液| 亚洲成人国产一区在线观看| 亚洲 欧美一区二区三区| 久久久国产精品麻豆| 欧美日本视频| 久久人妻熟女aⅴ| 亚洲一区中文字幕在线| 中亚洲国语对白在线视频| 精品国产亚洲在线| 亚洲一卡2卡3卡4卡5卡精品中文|