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

    The Study on Precise and Quantitative Measurement of Flame OH Concentration by CRDS-CARS-PLIF Techniques

    2023-12-13 06:20:38BAIBingCHENGuozhuYANGWenbinCHEQingfengWANGLinsenSUNWeiminCHENShuang
    光譜學與光譜分析 2023年12期

    BAI Bing, CHEN Guo-zhu, YANG Wen-bin, CHE Qing-feng, WANG Lin-sen, SUN Wei-min, CHEN Shuang*

    1. College of Physics and Optoelectronic Engineering, Harbin Engineering University, Harbin 150001, China 2. State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang 621000, China 3. Facility Design and Instrumentation Institute, China Aerodynamics Research and Development Center, Mianyang 621000, China

    Abstract In this paper, a pulsed cavity-ring down spectroscopy (CRDS) is employed to measure the quantitative concentration of the OH radical in a plane flame burner with premixed methane/air. By analyzing the cavity ring-down absorption spectrum theory, we select the P1(2) absorption line spectrum of the electronic transition band OH A2Σ+-X2Π(0,0) and build a set of the pulsed CRDS experimental device with a laser wavelength of 308.6 nm. The device of the pulsed CRDS is composed of a pair of mirrors with a reflectance of 99.7%, the cavity length of the ring-down cavity is 270 cm, and the ring-down time of the empty cavity (without a flame in the optical cavity) is 2.33 μs. By analyzing the experimental parameters that affect the precise measurement of concentration, we use Planar Laser Induced Fluorescence (PLIF), Coherent Anti-Stokes Raman Scattering (CARS), and the pulsed CRDS to measure the effective absorption length of OH, high temperature of the flame, and cavity ring-down time. When the premixed methane (1.1 L·min-1) and air (15 L·min-1) are burned in a flat flame burner, and at the height of 6 mm from the burner surface, the precisely measured effective absorption length by PLIF is 7.1% higher than that of directly choosing the diameter of the burner surface as the absorption length, the measured precision of the temperature by CARS is increased by 45% than that measured by the thermocouple under room temperature, the measured precision of the optical cavity ring-down time with flame in the cavity and non-OH absorption wavelength is improved by 21.6% than that measured time of cavity ring-down without a flame in the cavity. By combining the above measurement techniques to measure all experimental parameters precisely, we obtain that the number density of OH molecules (3.59×1013 molecules·cm-3) can reach the maximum value when the height from the furnace burner is 6 mm, and the precision of OH concentration is 35.6% higher than that of the unmodified OH concentration. Under different equivalence ratios (Φ=0.7~1.1), with the increase of the height from the burner surface, the number of OH particles gradually decreases, and the curve fitting shows that the OH concentration decreases in an e-exponential decay. At the same combustion height, the concentration of OH increases with the increase of equivalent ratios. When the methane flow rate is kept constant, the OH concentration in the oxygen-rich combustion condition is higher than in the low-oxygen combustion condition. In the combustion field, the precise measurement method with the multi-spectral technology (CRDS-CARS-PLIF) can achieve the precise quantitative measurement of OH concentration and provide technical support for the quantitative measurement of the concentration of other combustion product molecules, which plays a crucial role in the study of combustion chemical reactions.

    Keywords CRDS (cavity ring-down spectroscopy); Effective absorption length; High-temperature testing;Ring-down time; OH concentration

    Introduction

    Quantitative measurement of substance concentration in the combustion field is necessary to understand the chemical reaction of combustion. As an important intermediate combustion product, OH plays a vital role in studying the chemical combustion reaction. Laser-induced fluorescence (LIF) has been used to measure the number density of OH particles[1]. However, it is difficult to achieve the quantitative measurement of OH concentration by the LIF technique, and the main reasons are constant calibration, quenching of excited particles, and low one-way absorption intensity. Chen Shuang’s group[2]obtained the real OH fluorescence signal in experience by subtracting the difference between the original LIF signal and pure Rayleigh scattering signal to calibrate the undesired background from scattered laser light. However, the uncertainty of the obtained OH concentration is still as high as ±25%. Among them, the largest error comes from the uncertainty of LIF measurements, the measured temperature, and OH concentration distribution at the edges. In addition, the most serious difficulty in achieving quantitative measurement for LIF is first defining the energy transfer of excited states such as the total impact quenching rate of electronically excited states, and the energy transfer of rotational and vibration[3].

    As a direct absorption spectroscopy technique, cavityring-down spectroscopy (CRDS) has gained widespread attention for the quantitative measurement of substance concentration[4]. Compared with the LIF technique, the CRDS has the advantages of high absorption intensity without being affected by laser energy fluctuations and the increased length of the absorption path to improve test sensitivity[5-7]. CRDS can be briefly summarized as a multi-pass absorption technique whereby the sample is placed in a high-finesse cavity, and the intensity decay of the laser is detected at the output mirror when the laser traverses in the cavity. In 1994, Gerard Meijer et al.[8]applied the cavity ring-down spectroscopy technique to combustion diagnosis and measured the absorption spectrum of OH. In recent years, CRDS was also employed in microwave plasma-assisted combustion[9]to measure the absolute concentration of the ground state OH radical at the ignition region. In addition, CRDS was also suitable for the quantitative measurement of other combustion intermediate products, such as CH4[10], CN[11], NO[12], HCO[13], and CO[14]. Although the CRDS technique has been widely used, the analysis of the measured concentration error is relatively rare. At present, Tu Xiaobo et al.[15]have carried out concentration error analyses caused by experimental equipment such as laser linewidth and line-type, laser frequency stability, and the response of photomultiplier tubes. However, by analyzing the CRDS theories, we find that the precise measurement of experimental parameters will also cause great errors in CRDS concentration measurement, such as replacing effective absorption length with burner surface diameter[15], the measured temperature errors, and the influence of flame for cavity ring-down time.

    In this paper, the quantitative measurement of OH concentration is carried out by a pulsed CRDS technique on a McKenna burner with premixed CH4/air gases. Due to heat radiation and heat transfer[16], it is difficult for the thermocouple to achieve accurate temperature measurement at room temperature, so the high-precision CARS technique is used to measure the real temperature for the high-temperature flame. The PLIF technique obtains the flame structure distribution, and then the effective absorption length at different heights from the burner surface is measured. The P1(2) absorption line spectrum of the electronic transition band of OH A2Σ+-X2Π(0,0) of OH is selected to build a pulsed CRDS experimental device with the 308.6 nm laser wavelength to measure OH concentration. By analyzing the absorption spectrum theories and using the accurate experiment measurement, we realize the precise and quantitative measurement of OH concentration on the flat flame burner.

    1 Theory of CRDS

    The ring-down time of the optical cavity is an important parameter for CRDS, and its expression is as follows[17].

    (1)

    whereLis the cavity length of the optical cavity,cis the speed of light,Ris the reflectivity of the cavity mirror,αis the absorption coefficient of the sample,lis the absorption length,ξis the additional loss of the optical cavity.

    When the optical cavity is empty (without flame), the ring-down timeτ0can be expressed as

    (2)

    By substituting Eq.(2) into Eq.(1), the absorption coefficient expression of the sample can be obtained.

    (3)

    It can be seen from Eq.(3) that the sample absorption coefficient is mainly determined by the ring-down time and absorption length and is not affected by the laser intensity[5,18].

    In addition, the absorption coefficient of combustion products can be expressed as[15]

    (4)

    WherefBrepresents the Boltzmann coefficient corresponding to the ground state of combustion products,his the Planck constant,νis laser frequency,ηis the overlap factor, andB12is the absorption coefficient of Einstein.

    Eq.(3) is substituted into Eq.(4), and the concentration of the combustion product OH can be expressed as

    (5)

    It can be seen from Eq.(5) that the main factors affecting OH concentration are absorption length, decay time, and Boltzmann coefficient. However, the decay time is also affected by the additional loss of the optical cavity, so keeping the consistency of the combustion state can effectively eliminate the error caused by the additional loss. The LIFBASE database finds the Boltzmann coefficientfBcorresponding to each temperature under the steady-state condition. As shown in Fig.1, the absorption intensity of OH gradually decreases with the increase in temperature, which can indicate that temperature impacts the accuracy of OH concentration measurement.

    Fig.1 The absorption intensity of OH at different temperatures

    2 Experimental Setup

    In the complex spectral cackground, the absorption wavelength of OH is selected to avoid the highly overlapping absorption lines in the spectral bandwidth as much as possible, which can improve the signal-to-noise ratio and reduce the interference of adjacent spectral lines. The laser linewidth in the experiment is 0.3 cm-1, so the P1(2) absorption line in the electron transition band of A2Σ+-X2Π(0,0) of OH is selected as the absorption wavelength of OH.

    The pulsed CRDS experimental device is shown in Fig.2. The laser system consists of a solid laser, a dye laser, and a frequency-doubling device. The laser system outputs a 308.6 nm UV laser with a pulse width of 10 ns, a repetition frequency of 10 Hz, and a pulse energy of 3 mJ. The plane mirror M1 fully reflects the 617.2 nm laser and transmits the 308.6 nm laser. The wavelength meter monitors the reflected 617.2 nm laser to ensure that the wavelength of the 308.6 nm laser can transmit through M1 without shifting. The transmitter 308.6 nm UV laser first passes through the pinhole to filter stray light and then enters the ring-down cavity. The two mirrors of the ring-down cavity (M4 and M5) use concave mirrors with a curvatuse radius of -6 000 mm and have 99.7% reflectance near 308 nm. The cavity length is 2.7 m to avoid interference between two adjacent pulses in the cavity. A filter plate (M6) is placed behind the optical resonator, which filters the background noise and only lets the light around 308 nm pass through. The filtered ring-down signal is received by a photomultiplier tube, and then the data is collected by a high-speed Oscilloscope.

    Fig.2 A pulsed CRDS experimental device

    In the experiment, a flowmeter is used to control the flow rates of CH4/air/N2, N2is used to protect the flame stability, and the premixed CH4/air is completed before reaching the flat flame burner. The parameters of experimental cases are shown in Table 1. Five cases with different equivalence ratios are selected in the experiment, which can realize high-temperature combustion[19]. Under each case,the interval of 3 mm height is performed from close to the burner surface, and a total of 9 groups of ring-down signal data are measured.

    Table 1 Parameters of testing cases

    3 Results and Discussion

    3.1 Temperature measured by CARS

    The temperature must be precisely measured to get the precise Boltzmann fractionfB. In the combustion environment of case 1, the CARS technique and thermocouple are used to measure the temperature of different heights at the same step, and the temperature difference and the relative error offBare shown in Fig.3. As shown in Fig.3, the larger the temperature difference between CARS and the thermocouple, the larger the relative error offB. With the increased distance from the burner surface, the temperature difference decreases gradually. When the distance from the burner surface is 6 mm, the temperature difference between the measuring point of CARS and the thermocouple is the highest, and the relative error offBcan reach 45%, which indicates that the deviation of the measured temperature by the thermocouple in the high-temperature environment is large. The thermocouple cannot meet the requirements of temperature measurement.

    Fig.3 The temperature difference between CARS and the thermocouple (red dot) and fB relative error (blue dot) at different heights for case 1

    To illustrate the accuracy of temperature measurement of the CARS system in the laboratory[19], five cases in the DLR[20](Deutsches Zentrum fur Luft-und Raumfahrt) are selected and measured for comparison. The results of the measured temperature by CARS are shown in Fig.4(a), and the results indicate that the measured temperature by CARS in the laboratory is completely consistent with those in the literature. Therefore, the CARS system is used to measure the temperature of the remaining casesat different heights, and the temperature results are shown in Fig.4(b).

    Fig.4 The measured temperatures by CARS(a): The relationship between DLR literature temperature and CARS temperature; (b): CARS temperatures of different cases

    3.2 Effective absorption length measured by PLIF

    To accurately measure the effective absorption length of OH, the PLIF technique is first used to obtain the flame structure under different cases, and then the absorption lengths of OH at different heights are measured. As shown in Fig.5(a), the PLIF technique is used to obtain the flame structure of case 1, and the red line in Fig.5(a) is the effective absorption length of OH when the height from the furnace surface is 15 mm. It can be seen from Fig.5(b) that the effective absorption length of OH decreases with the increase of distance from the burner surface. The relative error of absorption length increases with height increasing. If the effective absorption length is replaced by the furnace face diameter, it will cause a 7.1%~39.5% relative error, and the higher the height, the greater the error. Therefore, the effective absorption length of different cases must be accurately measured by the PLIF technique.

    Fig.5 (a) The measured active absorption length (red line) by PLIF at the height of 15 mm; (b) Active absorption lengths and relative errors at different heights for case 1

    Figure 6 shows the measured effective absorption lengths by the PLIF technique at different heights from the burner surface for the remaining cases. As can be seen from case 4 in Fig.6, when the gas flow rate is very high, only a part of the effective absorption length near the burner surface is equal to the diameter of the burner surface.

    Fig.6 Active absorption lengthsindifferent cases

    3.3 OH concentration measured by CRDS

    When the optical cavity is empty, the signal of the ring-down cavity is shown in the upper right corner of Fig.7(a). In Fig.7(a), the round-trip propagation time of the optical pulse in the cavity is 18 ns,which is larger than the laser pulse width of 10 ns, and the interference of two adjacent pulses in the cavity can be effectively avoided. As shown in Fig.7(b), the ring-down timeτ0of the empty cavity is 2.33 μs by selecting the maximum value of the ring-down signal and e-index fitting.

    Fig.7 Cavity ring-down signal without a flame in the cavity(a): Raw data; (b): Signal fitting

    To verify the influence withouta flame in the cavity for the ring-down timeτ0, it is also necessary to measureτ0when there is a flame in the optical cavity. When the flat flame burner is ignited in the cavity, and the laser wavelength is adjusted to the non-absorbing wavelength of OH, the ring-down timeτ0of case 1 is obtained at different heights from the burner surface, as shown in red star in Fig.8(a).It can be seen from Fig.8 (a) thatτ0of the empty cavity is greater than that of the flame inthe cavity, which indicates that the flame in the cavity will affect the value ofτ0by the additional loss of the optical cavity. As shown in Fig.8(b), the relative error of cavity ring-down timeτ0at different heights increases by 12.8%~21.6% compared with the ring-down time with flame in the combustion environment for case 1. It is mainly due to the Rayleigh scattering[15]caused by impurity particles in the flame combustion process, which causes a change in the additional lossξ. Moreover, it can be seen from Eq.(2) that the additional lossξwill affect the value of the ring-down timeτ0, but the ring-down timeτ0of the empty cavity ignores the loss caused by the flame impurity particles.

    Fig.8 (a) The influence for τ0 with and without a flame in the cavity, (b) τ0 relative error for case 1

    The overlap factorηis 2.188 cm andB12is 2.278×109m2·J-1·s-1, which can be obtained by checking the LIFBASE database. Then the OH concentration will be obtained by substituting the above parameters into Eq. (5). As shown in Fig.9, the OH concentration at different height positions is obtained by accurate measurement of experimental parameters in the combustion environment of case 1. The OH concentration increases by 35.6%~47.9% compared with the uncorrected one.

    Fig.9 The concentration of OH and relative error for case 1

    The OH concentration and fitting curves at different cases are shown in Fig.10. As can be seen from Fig.10, with the increase of heights from the burner surface, the concentration of OH gradually decreases, which is the same trend as the tested results in the literature[15,21]. It is found that the concentration of OH decreases exponentially with the increase of heights from the burner surface by fitting the curve. By comparing cases 1, 2, 3, and 4 at the same combustion height, we find that the concentration of OH increases with the increase of the equivalent ratio under the condition of oxygen-rich combustion. The comparison of cases 2 and 5 shows that the OH concentration in the oxygen-rich state is higher than that in the low-oxygen state when the CH4flow rate remains constant.

    Fig.10 OH concentrations and curve fittings at different cases

    4 Conclusion

    The experimental results show that CARS, PLIF, and the pulsed CRDS have successfully measured the OH concentration in a flat flame burner with the premixed CH4/air. A set of the pulsed CRDS experimental device with the 308.6 nm laser wavelength is constructed by selecting the P1(2) absorption line spectrum of the electronic transition band OH A2Σ+-X2Π(0,0). By analyzing the spectrum theory,we find that the experimental parameters that affect the precise measurement of concentration are the effective absorption length of OH, temperature, and decay time. Planar laser-induced fluorescence (PLIF) is used to measure the changed effective absorption length of the flame structure. Coherent anti-Stokes Raman scattering (CARS) is used to measure the temperature field of the high-temperature flame accurately. The pulsed CRDS is used to measure the ring-down time by igniting the flame inthe optical cavity, which is affected by the additional loss. When the premixed CH4(1.1 L·min-1) and air (15 L·min-1) are burned on the flat flame burner, and the measured height from the burner surface is 6 mm, the measured accuracy of the effective absorption length by the PLIF technique is 7.1% higher than that by directly choosing the diameter of the burner surface as the absorption length, the measured accuracy of thetemperature by CARS technique is 45% higher than that measured by a thermocouple at room temperature, the measured precision of the cavity ring-down time is improved by 21.6% than that measured ring-down time of the empty cavity when the flame is ignited in the cavity, and non-OH absorption wavelength is selected. By using a combination of the above technologies to achieve precise measuring experiment parameters, we find that the number density of OH molecules (3.59×1013molecules·cm-3) can reach the maximum value when the height from the burner surface is 6 mm, and the precision of OH concentration is 35.6% higher than that of the unmodified OH concentration. Under the different equivalence ratios (Φ=0.7~1.1), with the increase of the heights from the burner surface, the number of OH particles gradually decreases,and the curve fitting shows that the OH concentration decreases in an e-exponential decay. Under the condition of oxygen-rich combustion, the concentration of OH at the same height increases with the increasing equivalent ratios. When the CH4flow rate remains constant, the OH concentration in the oxygen-rich combustion condition will be higher than in the low-oxygen combustion condition.The precise measurement of OH concentration is improved by combining the CRDS-CARS-PLIF multi-technique, which can provide strong data support for chemical combustion analysis.

    久9热在线精品视频| 国产一区二区在线观看av| 日韩免费av在线播放| 久久久久视频综合| 99re6热这里在线精品视频| 亚洲专区字幕在线| 9色porny在线观看| 久久青草综合色| 丁香六月欧美| 免费一级毛片在线播放高清视频 | 最新在线观看一区二区三区| 啪啪无遮挡十八禁网站| 色精品久久人妻99蜜桃| cao死你这个sao货| 黄片大片在线免费观看| 桃红色精品国产亚洲av| 99香蕉大伊视频| 成人国语在线视频| 国产精品 国内视频| 亚洲一码二码三码区别大吗| 91国产中文字幕| 黄频高清免费视频| 电影成人av| 国产精品亚洲av一区麻豆| 久久 成人 亚洲| 建设人人有责人人尽责人人享有的| 搡老岳熟女国产| 国产精品av久久久久免费| 亚洲人成电影免费在线| 一区二区三区国产精品乱码| 99久久人妻综合| 国产在线免费精品| 女警被强在线播放| 精品国产国语对白av| 狠狠精品人妻久久久久久综合| 国产精品二区激情视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产一区二区 视频在线| 麻豆成人av在线观看| 国产又爽黄色视频| 男女下面插进去视频免费观看| 女人精品久久久久毛片| 天天躁日日躁夜夜躁夜夜| 日韩欧美国产一区二区入口| 亚洲av日韩精品久久久久久密| 一区在线观看完整版| 高潮久久久久久久久久久不卡| 亚洲中文字幕日韩| 高清视频免费观看一区二区| a在线观看视频网站| 夫妻午夜视频| 999久久久精品免费观看国产| kizo精华| 日韩大码丰满熟妇| 久久ye,这里只有精品| 水蜜桃什么品种好| 两个人免费观看高清视频| 国产成人精品无人区| 久久久久国产一级毛片高清牌| 国产精品免费大片| 久久久久网色| 天天躁狠狠躁夜夜躁狠狠躁| 国产欧美日韩精品亚洲av| 在线观看免费日韩欧美大片| 日日夜夜操网爽| 午夜视频精品福利| 一本一本久久a久久精品综合妖精| 亚洲,欧美精品.| 国产男靠女视频免费网站| 两性午夜刺激爽爽歪歪视频在线观看 | www.999成人在线观看| 欧美精品人与动牲交sv欧美| 最近最新中文字幕大全电影3 | 天天躁夜夜躁狠狠躁躁| 人人妻人人澡人人爽人人夜夜| 欧美精品av麻豆av| 啦啦啦在线免费观看视频4| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成77777在线视频| 老司机午夜十八禁免费视频| 老司机亚洲免费影院| 在线观看一区二区三区激情| 一级,二级,三级黄色视频| 久久久精品94久久精品| 国产一区二区三区视频了| 日本精品一区二区三区蜜桃| 日本vs欧美在线观看视频| 国产精品国产高清国产av | 欧美黄色片欧美黄色片| 日本av免费视频播放| 我的亚洲天堂| 成人av一区二区三区在线看| 亚洲精品一二三| 电影成人av| 国产aⅴ精品一区二区三区波| 成人特级黄色片久久久久久久 | 久久久久久久精品吃奶| 两个人免费观看高清视频| 国产人伦9x9x在线观看| 久久人人爽av亚洲精品天堂| 午夜91福利影院| 久久婷婷成人综合色麻豆| 国产又色又爽无遮挡免费看| 国产三级黄色录像| av网站免费在线观看视频| 国产欧美日韩综合在线一区二区| svipshipincom国产片| 热re99久久国产66热| 欧美日韩精品网址| 高清视频免费观看一区二区| 十分钟在线观看高清视频www| 在线十欧美十亚洲十日本专区| 日韩 欧美 亚洲 中文字幕| 天堂中文最新版在线下载| 国产成人免费观看mmmm| 女同久久另类99精品国产91| 国产精品九九99| 精品国产一区二区久久| a级毛片黄视频| 少妇被粗大的猛进出69影院| 制服诱惑二区| av不卡在线播放| 一区福利在线观看| 亚洲一区二区三区欧美精品| 精品少妇内射三级| 精品视频人人做人人爽| 十八禁人妻一区二区| 91老司机精品| 丁香六月天网| 99riav亚洲国产免费| 久久人人爽av亚洲精品天堂| 一个人免费看片子| 亚洲av欧美aⅴ国产| 国产精品一区二区在线不卡| 日本五十路高清| 亚洲av第一区精品v没综合| 多毛熟女@视频| 亚洲中文字幕日韩| 中文字幕人妻丝袜制服| 久久热在线av| 黑人巨大精品欧美一区二区蜜桃| 麻豆av在线久日| 国产aⅴ精品一区二区三区波| 成人av一区二区三区在线看| 欧美在线一区亚洲| 9色porny在线观看| 国产成人系列免费观看| 欧美成狂野欧美在线观看| 窝窝影院91人妻| 日本五十路高清| av网站免费在线观看视频| 国产精品一区二区免费欧美| 亚洲成人免费av在线播放| 日韩人妻精品一区2区三区| 国产日韩一区二区三区精品不卡| 女人高潮潮喷娇喘18禁视频| 亚洲一区二区三区欧美精品| 99九九在线精品视频| 日韩精品免费视频一区二区三区| 久久性视频一级片| 国产1区2区3区精品| 国产野战对白在线观看| 丝袜人妻中文字幕| 日日爽夜夜爽网站| 午夜福利一区二区在线看| 在线观看免费日韩欧美大片| 一本综合久久免费| 久久久欧美国产精品| 国产熟女午夜一区二区三区| 国产亚洲精品一区二区www | 啪啪无遮挡十八禁网站| 丁香六月欧美| 国产免费av片在线观看野外av| 麻豆国产av国片精品| 精品卡一卡二卡四卡免费| av欧美777| tocl精华| www日本在线高清视频| 国产欧美日韩一区二区三| 久久香蕉激情| 亚洲专区国产一区二区| 精品久久久久久久毛片微露脸| 精品欧美一区二区三区在线| 久久久精品区二区三区| 高清黄色对白视频在线免费看| 黄片播放在线免费| 免费少妇av软件| 欧美成狂野欧美在线观看| 国产亚洲精品一区二区www | 成年人黄色毛片网站| 91精品国产国语对白视频| 日韩视频在线欧美| 国产熟女午夜一区二区三区| 国产精品1区2区在线观看. | 中亚洲国语对白在线视频| 亚洲欧美日韩高清在线视频 | 久久午夜综合久久蜜桃| 亚洲伊人色综图| 啪啪无遮挡十八禁网站| 搡老熟女国产l中国老女人| 人妻 亚洲 视频| 99精品久久久久人妻精品| 热99国产精品久久久久久7| 18禁国产床啪视频网站| 亚洲成国产人片在线观看| 亚洲欧美日韩另类电影网站| 天堂中文最新版在线下载| 人人澡人人妻人| 在线观看人妻少妇| 最近最新免费中文字幕在线| 成人三级做爰电影| 久久久久久久久免费视频了| 在线观看舔阴道视频| 侵犯人妻中文字幕一二三四区| 日韩一区二区三区影片| 久久国产精品男人的天堂亚洲| 亚洲性夜色夜夜综合| 极品人妻少妇av视频| 91精品国产国语对白视频| 亚洲 国产 在线| 国产欧美亚洲国产| 日本av免费视频播放| 青青草视频在线视频观看| 中文亚洲av片在线观看爽 | 性高湖久久久久久久久免费观看| 一区在线观看完整版| 色视频在线一区二区三区| 亚洲成人免费av在线播放| 黄片大片在线免费观看| 久久精品亚洲熟妇少妇任你| av有码第一页| 午夜免费鲁丝| 1024香蕉在线观看| 老司机福利观看| 日韩免费av在线播放| 电影成人av| 亚洲午夜理论影院| 久久久久精品人妻al黑| 丰满人妻熟妇乱又伦精品不卡| 国产精品影院久久| 精品国产乱码久久久久久小说| 久久久久网色| 午夜福利免费观看在线| 91老司机精品| 蜜桃在线观看..| 亚洲成a人片在线一区二区| 在线观看www视频免费| 天堂中文最新版在线下载| 看免费av毛片| 国产免费av片在线观看野外av| 老司机午夜十八禁免费视频| 三上悠亚av全集在线观看| 麻豆av在线久日| 无限看片的www在线观看| 国产一区二区三区综合在线观看| 欧美精品av麻豆av| 久久精品aⅴ一区二区三区四区| 视频区图区小说| 亚洲专区国产一区二区| 肉色欧美久久久久久久蜜桃| 国产成人免费无遮挡视频| 久久天堂一区二区三区四区| 国产精品亚洲av一区麻豆| 少妇粗大呻吟视频| 久久午夜亚洲精品久久| 变态另类成人亚洲欧美熟女 | 一二三四在线观看免费中文在| 色94色欧美一区二区| 成人影院久久| 久久狼人影院| 国产精品亚洲av一区麻豆| 老熟妇仑乱视频hdxx| 国产精品 国内视频| 一本久久精品| 国产精品一区二区在线观看99| 国产高清国产精品国产三级| 精品视频人人做人人爽| 91麻豆av在线| 在线观看一区二区三区激情| 最新在线观看一区二区三区| 99国产极品粉嫩在线观看| 国内毛片毛片毛片毛片毛片| 久久精品国产综合久久久| 黄片小视频在线播放| av又黄又爽大尺度在线免费看| 嫩草影视91久久| 欧美乱妇无乱码| 女性被躁到高潮视频| 国产一区二区三区在线臀色熟女 | 在线播放国产精品三级| 啦啦啦免费观看视频1| 极品人妻少妇av视频| 老司机靠b影院| 久久久久国内视频| av不卡在线播放| 老司机影院毛片| 天堂8中文在线网| 日韩欧美一区二区三区在线观看 | 免费av中文字幕在线| 精品亚洲成a人片在线观看| 男人操女人黄网站| 中文亚洲av片在线观看爽 | 国产国语露脸激情在线看| 色播在线永久视频| 中文字幕制服av| 国产av精品麻豆| 99久久国产精品久久久| 久久影院123| 一级a爱视频在线免费观看| 亚洲视频免费观看视频| 欧美成狂野欧美在线观看| 9191精品国产免费久久| 久热这里只有精品99| tocl精华| 我要看黄色一级片免费的| av视频免费观看在线观看| 99精品久久久久人妻精品| 免费在线观看视频国产中文字幕亚洲| 日韩视频在线欧美| 欧美精品一区二区免费开放| 丰满少妇做爰视频| 国产真人三级小视频在线观看| 青草久久国产| 午夜日韩欧美国产| h视频一区二区三区| 男人操女人黄网站| 日本黄色视频三级网站网址 | 亚洲欧洲日产国产| 久久精品熟女亚洲av麻豆精品| 亚洲成a人片在线一区二区| 亚洲专区中文字幕在线| 人人妻人人澡人人看| 少妇粗大呻吟视频| 日韩有码中文字幕| 三上悠亚av全集在线观看| 男人操女人黄网站| 午夜视频精品福利| 日本一区二区免费在线视频| 伊人久久大香线蕉亚洲五| 亚洲色图 男人天堂 中文字幕| 亚洲五月色婷婷综合| 男女免费视频国产| 一边摸一边抽搐一进一出视频| 黄色怎么调成土黄色| 老司机靠b影院| 啦啦啦中文免费视频观看日本| 一级毛片女人18水好多| 午夜免费成人在线视频| 一级,二级,三级黄色视频| 中文字幕av电影在线播放| 好男人电影高清在线观看| 电影成人av| 99国产精品99久久久久| 日韩精品免费视频一区二区三区| 欧美在线一区亚洲| 狠狠精品人妻久久久久久综合| 中亚洲国语对白在线视频| 国产精品国产高清国产av | 久久久久久久久久久久大奶| 欧美精品亚洲一区二区| 这个男人来自地球电影免费观看| 女同久久另类99精品国产91| 桃红色精品国产亚洲av| 国产精品1区2区在线观看. | 麻豆av在线久日| 国产成人免费无遮挡视频| 日本精品一区二区三区蜜桃| 色尼玛亚洲综合影院| 国产高清视频在线播放一区| 免费人妻精品一区二区三区视频| 天堂俺去俺来也www色官网| 一区二区三区精品91| 建设人人有责人人尽责人人享有的| 久久免费观看电影| 亚洲欧美日韩另类电影网站| 久久性视频一级片| 亚洲精品国产区一区二| 大香蕉久久成人网| 999久久久精品免费观看国产| 两性夫妻黄色片| 欧美日韩黄片免| 久久精品亚洲av国产电影网| 欧美日韩亚洲高清精品| 国产一区二区在线观看av| 日韩中文字幕欧美一区二区| 亚洲九九香蕉| 老汉色av国产亚洲站长工具| 欧美在线一区亚洲| 丰满迷人的少妇在线观看| 欧美成人午夜精品| 亚洲精品美女久久av网站| 成人国语在线视频| 亚洲中文字幕日韩| 国内毛片毛片毛片毛片毛片| 人人澡人人妻人| 亚洲国产毛片av蜜桃av| 亚洲一卡2卡3卡4卡5卡精品中文| 成年人午夜在线观看视频| www.999成人在线观看| 大片电影免费在线观看免费| 每晚都被弄得嗷嗷叫到高潮| 国产免费福利视频在线观看| xxxhd国产人妻xxx| 欧美日韩黄片免| 51午夜福利影视在线观看| 久久热在线av| avwww免费| 欧美人与性动交α欧美精品济南到| 精品国产亚洲在线| 成年人黄色毛片网站| 国产高清视频在线播放一区| 99久久精品国产亚洲精品| 免费看a级黄色片| 日韩中文字幕视频在线看片| 久久久精品免费免费高清| 中文字幕最新亚洲高清| 日韩三级视频一区二区三区| 在线观看舔阴道视频| 中文字幕另类日韩欧美亚洲嫩草| 岛国在线观看网站| 最近最新免费中文字幕在线| 国产精品久久久久久精品古装| 波多野结衣一区麻豆| 好男人电影高清在线观看| e午夜精品久久久久久久| tube8黄色片| 亚洲国产中文字幕在线视频| 99精品欧美一区二区三区四区| 黄色怎么调成土黄色| 极品少妇高潮喷水抽搐| 大香蕉久久成人网| 国产精品一区二区免费欧美| 法律面前人人平等表现在哪些方面| 丰满饥渴人妻一区二区三| 久久国产精品大桥未久av| 精品欧美一区二区三区在线| 亚洲精品在线美女| 亚洲黑人精品在线| 50天的宝宝边吃奶边哭怎么回事| 亚洲美女黄片视频| 99香蕉大伊视频| 黑人巨大精品欧美一区二区mp4| 国产有黄有色有爽视频| 亚洲色图综合在线观看| kizo精华| 国产片内射在线| av在线播放免费不卡| 99热国产这里只有精品6| 亚洲三区欧美一区| 免费少妇av软件| 国产欧美亚洲国产| 正在播放国产对白刺激| tube8黄色片| a级毛片在线看网站| 91九色精品人成在线观看| 在线天堂中文资源库| 国产一区二区 视频在线| 午夜福利乱码中文字幕| 欧美一级毛片孕妇| www.自偷自拍.com| 亚洲熟女精品中文字幕| 热re99久久精品国产66热6| 国产色视频综合| 亚洲专区中文字幕在线| 777米奇影视久久| 人妻 亚洲 视频| 日本av手机在线免费观看| 黑丝袜美女国产一区| 岛国在线观看网站| 两性夫妻黄色片| 91字幕亚洲| 首页视频小说图片口味搜索| 在线永久观看黄色视频| 一级毛片女人18水好多| 国内毛片毛片毛片毛片毛片| 亚洲五月色婷婷综合| 亚洲人成电影免费在线| 久久国产精品男人的天堂亚洲| 亚洲欧洲精品一区二区精品久久久| 久久久水蜜桃国产精品网| 黄色毛片三级朝国网站| 少妇裸体淫交视频免费看高清 | 69精品国产乱码久久久| 精品亚洲成a人片在线观看| 午夜成年电影在线免费观看| 国产欧美亚洲国产| 久久久久网色| 男人操女人黄网站| 国产精品久久久av美女十八| 亚洲情色 制服丝袜| 午夜福利免费观看在线| 手机成人av网站| 青草久久国产| 女人高潮潮喷娇喘18禁视频| 亚洲精品国产色婷婷电影| 久久久久精品人妻al黑| 欧美性长视频在线观看| 香蕉国产在线看| 色在线成人网| 久久热在线av| 91av网站免费观看| 精品国产亚洲在线| 一本—道久久a久久精品蜜桃钙片| 精品一品国产午夜福利视频| 国产av国产精品国产| 少妇裸体淫交视频免费看高清 | 亚洲情色 制服丝袜| 亚洲性夜色夜夜综合| 国产精品98久久久久久宅男小说| av福利片在线| 欧美精品高潮呻吟av久久| 两个人看的免费小视频| 欧美激情高清一区二区三区| 女同久久另类99精品国产91| 菩萨蛮人人尽说江南好唐韦庄| 国产淫语在线视频| 啦啦啦中文免费视频观看日本| 精品免费久久久久久久清纯 | 中文字幕人妻丝袜一区二区| 欧美日韩国产mv在线观看视频| av网站免费在线观看视频| 中文字幕制服av| 2018国产大陆天天弄谢| 99re6热这里在线精品视频| 夜夜骑夜夜射夜夜干| 高清毛片免费观看视频网站 | 丰满少妇做爰视频| 又大又爽又粗| 国产深夜福利视频在线观看| 国产片内射在线| 精品久久久久久久毛片微露脸| 亚洲专区国产一区二区| 精品久久久久久电影网| 久久久精品免费免费高清| 成人精品一区二区免费| av网站免费在线观看视频| 亚洲专区中文字幕在线| 一级毛片电影观看| videosex国产| 777久久人妻少妇嫩草av网站| 99精品久久久久人妻精品| 日韩 欧美 亚洲 中文字幕| 久久99热这里只频精品6学生| 免费观看av网站的网址| 国产精品久久久久久精品古装| 日韩中文字幕视频在线看片| 日本黄色日本黄色录像| 五月天丁香电影| 一区在线观看完整版| 午夜福利欧美成人| 久久 成人 亚洲| 男女边摸边吃奶| 亚洲伊人色综图| 成人18禁高潮啪啪吃奶动态图| 精品国产乱码久久久久久男人| 999久久久精品免费观看国产| 欧美一级毛片孕妇| 大香蕉久久网| 国产精品一区二区精品视频观看| 精品熟女少妇八av免费久了| av在线播放免费不卡| 99国产精品99久久久久| 一级,二级,三级黄色视频| 99久久人妻综合| 成人影院久久| 男女之事视频高清在线观看| 女人被躁到高潮嗷嗷叫费观| 高清av免费在线| 国产一区二区激情短视频| 老司机在亚洲福利影院| 18禁黄网站禁片午夜丰满| 少妇被粗大的猛进出69影院| 香蕉久久夜色| 免费观看a级毛片全部| 久久亚洲精品不卡| 亚洲欧美一区二区三区久久| 欧美日韩成人在线一区二区| 男女午夜视频在线观看| 我的亚洲天堂| 亚洲av欧美aⅴ国产| 国产在线一区二区三区精| 久久影院123| 精品免费久久久久久久清纯 | 国产三级黄色录像| 女同久久另类99精品国产91| 国产成人欧美在线观看 | 无人区码免费观看不卡 | 王馨瑶露胸无遮挡在线观看| 久久精品亚洲av国产电影网| 男女无遮挡免费网站观看| 色播在线永久视频| 欧美乱妇无乱码| 久久狼人影院| 亚洲av成人一区二区三| 桃红色精品国产亚洲av| 国产又爽黄色视频| 亚洲中文日韩欧美视频| 免费人妻精品一区二区三区视频| 日韩视频在线欧美| 王馨瑶露胸无遮挡在线观看| 欧美成狂野欧美在线观看| h视频一区二区三区| 搡老熟女国产l中国老女人| 一区在线观看完整版| 在线永久观看黄色视频| 一本大道久久a久久精品| 国产熟女午夜一区二区三区| 日韩欧美一区二区三区在线观看 | 天堂俺去俺来也www色官网| 精品午夜福利视频在线观看一区 | 国产精品一区二区免费欧美| 欧美在线一区亚洲| 亚洲少妇的诱惑av| 91老司机精品| 色精品久久人妻99蜜桃|