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

    Multi-band polarization switch based on magnetic fluid filled dual-core photonic crystal fiber

    2023-03-13 09:18:52LianzhenZhang張連震XuedianZhang張學典XiantongYu俞憲同XuejingLiu劉學靜JunZhou周軍MinChang常敏NaYang楊娜andJiaDu杜嘉
    Chinese Physics B 2023年2期
    關(guān)鍵詞:周軍

    Lianzhen Zhang(張連震), Xuedian Zhang(張學典), Xiantong Yu(俞憲同), Xuejing Liu(劉學靜),Jun Zhou(周軍), Min Chang(常敏), Na Yang(楊娜), and Jia Du(杜嘉)

    Key Laboratory of Optical Technology and Instrument for Medicine,Ministry of Education,University of Shanghai for Science and Technology,Shanghai 200093,China

    Keywords: polarization switch,magnetic fluid,dual-core photonic crystal fiber

    1.Introduction

    As an electromagnetic wave, light possesses several basic characteristics, which include phase, intensity, polarization,and wavelength.[1]Since the polarization characteristic is a vector with multi-dimensional characteristics, its regulation requires more advanced devices.[2]In recent years, with the improvement of the ability to control the polarization characteristics of light,the unique advantages of polarization optical methods are attracting more attention.Nowadays, polarized optics has been widely used in the fields of optical communication, optical sensing and biomedicine.[2-4]In order to further tap the potential of polarized light methods,the control of light polarization state is still a current research hotspot.At present,the polarization state of light in space is mainly controlled by birefringent elements,such as polarization splitting prism, thin-film polarization beam splitters (PBS), and wave plates.[5,6]For the development of integrated photonics, controlling the polarization state of light in silicon-based waveguides is of great significance.Many silicon-based waveguide PBS and rotators with excellent performance have been reported.[7,8]With the increasing application of optical fibers in modern optical systems, it is of great significance to study the polarization controllers based on optical fibers.At present,there have been many reports on the research of fiber PBS and fiber polarization filter,[9-15]but there are still some deficiencies in these reports.Firstly, the cost of introducing elliptical holes or inconsistent round holes to increase birefringence is that it will be a great challenge to fabricate fibers.The coexistence of different shapes holes is challenging in photonic crystal fiber (PCF) manufacturing.[16]Then, in pursuit of structural miniaturization,its compatibility with conventional fiber systems is ignored.Finally,the device cannot convert the polarization state of the output modes.In addition,there are few reports on all-fiber polarization switch.

    Recently, we have witnessed the rapid development of high-performance all-fiber devices in the PCF platforms, combining with different optoelectronic materials,such as metals,[17-19]two-dimensional materials,[20-22]and liquids,[23-25]to achieve different functions.The cladding of PCF is composed of regularly arranged air holes,which facilitate the combination of PCF with different optoelectronic materials.PCF,filled with functional materials, will exhibit different light guiding properties.Optical fiber sensing,[26]light modulators,[25,27]etc.can be designed according to the principles.The combination of magnetic fluid(MF)and PCF,by full or selective filling, has attracted widespread attention.Many related devices have been reported, such as magnetic field sensors,[28-30]magneto-optical modulators,[31]and magnetooptical isolators.[32]

    In this paper,we proposed a polarization switch based on dual-core photonic crystal fiber(DCPCF)with same diameter cladding air holes, and analyzed its characteristics with magnetic field or not,respectively.The result shows that,in a free magnetic condition or its strength is less than a critical value,it can be used as a PBS.However,when the value is more than the critical value,the power of the received polarization mode decreases or increases at out ports.Therefore, the extinction ratio(ER)of the output polarization mode can be ameliorated by adjusting the magnitude of the magnetic field,and the polarization state can be rotated under a certain magnetic field magnitude.In addition,we analyzed the influence of the fiber length interception on the performance of the device.The results show that the switch has a good fault tolerance for the interception of the fiber length.

    2.Model and theory

    Figure 1 shows the schematics illustration of the polarization controller based on DCPCF with uniform circular holes.The air holes in the cladding of the PCF are periodically arranged in a hexagonal shape, the diameter of all the air holes isd,and the distance between the two holes isΛ.Two cores,C1 and C2, are formed by removing two air holes from the periodically arranged air holes.The base material of the PCF is silica,and its dispersion model can be described by the Sellmeier equation.[33]The eight holes, painted in blue in Fig.1,are filled with MF.When there is no magnetic field around the MF or the magnetic field magnitude is lower than itsHc, the refractive index of the MF does not change.However, when the surrounding magnetic field magnitude exceedsHc,the refractive index can be changed by adjusting the magnitude of the magnetic field.The relationship between the refractive index of the MF and the strength of the magnetic field can be described by the Langevin function[34,35]

    wherenois the refractive index of the MF when the ambient magnetic field strength is less thanHc,nsis the saturation refractive index of the MF.HandTare the magnetic field strength(in unit Oe,1 Oe=79.5775 A·m-1)and temperature(in unit K)of the environment,respectively,andαrepresents the fitting parameters.In this paper,a water-based MF with a particle volume concentration of 0.53 emu/g was used,andHcandTis assumed 43 Oe and 294 K,respectively.[15]

    Combining mode coupling theory and super-mode theory,[36]when the DCPCF is viewed as a monolithic waveguide,each mode in the waveguide will form two supermodes,namely the symmetric mode (y-even,x-even) and the asymmetric mode(y-odd,x-odd).The coupling law of DCPCF can be analyzed by the superposition of supermodes when the light with the energy ofPinis input from the core C1.The most energy is transferred to C2 when the phase difference between the symmetric and antisymmetric modes isπ.The two cores in the DCPCF are identical,so the phase matching condition is satisfied.Theoretically, the complete coupling of energy can be achieved,and the coupling distance isLx,y,

    whereβx,ysandβx,yaare the propagation constants of the four supermodes,β=2πneff/λ,neffis the effective refractive index,andλis the wavelength in free space.Assuming the light is input in C1,then the normalized energy ofxorypolarization mode output by C1 is given by

    To achieve complete separation ofx- andy-polarized light, it is necessary to satisfy the propagation distancez=mLx=nLy,wheremandnare positive integers with different parities.It can be concluded that the device length is determined by the ratio ofLxtoLy,which is defined as the coupling length ratio(CLR):

    Fig.1.Schematic diagram and cross section of the polarization switch based on DCPCF.When the magnetic field is none or its strength is less than a critical value, the two polarization (pol) modes, y-pol (red) and x-pol(blue), are split into the two cores at the output ports.When the magnetic field strength is more than the critical value, the ER of the output polarization mode,y-pol(red)and x-pol(blue),can not only be adjusted with the change of the magnetic field,but also the polarization states can be rotated.

    For conventional PBS,the ER is a key parameter to measure its performance.ER is the ratio of the energy of the polarized light to be eliminated in a fiber core at the output port to the obtained.The calculation method is as follows:

    The separation of the two polarization modes in the fiber is achieved because of the difference in their coupling distance from one core to the other.The reason for the different coupling distances, however, is their different propagation constants in the fiber.Therefore, it can be known that the polarization state of the beam at the output end can be changed by changing the length of the device or the propagation constant of the polarization modes.

    3.Analysis and results

    3.1.Optimization of the fiber structure

    In this article, a finite element mode solver (COMSOL)is adopted to optimize the parameters and performance of the DCPCF.Most of the current reports on PBS have introduced irregular air holes,such as elliptical air holes,circular air holes with different diameters, in the cladding of DCPCF.As a result, more parameters need to be analyzed when optimizing the device.The diameters of the air holes in the cladding of the DCPCF designed in this paper are consistent.In this case,we only need to analyze the influence of the air hole arrangement periodΛand the duty cycled/Λon its performance.In addition,this approach can reduce the complexity of the fiber fabrication process compared to using non-uniform pores.

    Fig.2.At the wavelength of 1.55μm: (a)the effect of duty ratio of the air holes on CL and CLR with the parameters of Λ =4μm,(b)the influence of Λ on CL and CLR with the parameter of d/Λ =0.44.

    Figure 2 depicts the influence of duty ratio and the hole spacing on the coupling lengths and CLR.Figure 2(a) shows the effect of duty ratio of the air holes on CL and CLR with the parameter ofΛ=4μm.The results show thatLxandLyincrease with the increase of duty cycle,while CLR decreases.Because the coupling lengthLx&lt;Ly, the smaller the CLR,the shorter the length of the device will be.[10]When the duty cycle in the PCF is less than 0.45,[37]the single-mode transmission with infinite cutoff can be realized,sod/Λ=0.44 is set.Figure 2(b)shows that the increase of the air hole spacing also leads to the increase ofLx,ywhen the duty cycle is fixed atd/Λ=0.44.Although the CLR also increases,the change is not significant.The above results first show that this method can be used for the modulation of the polarization state.In addition,the results show that after we determine the duty cycle of the air holes, the change of the air hole spacing has little effect on the CLR, which will be beneficial for us to design polarization controllers of various sizes.In order to be compatible with the traditional single-mode fibre communication system,Λ=5.25μm.

    According to the optimized parameters,the device is subjected to mode analysis, and the distribution of transmission modes (y-even,y-odd,x-even, andx-odd) of the device at a wavelength of 1550 nm is obtained as shown in Fig.3.The electric field of the even mode is in the same direction (symmetric mode), and the direction of the odd mode is opposite(asymmetric mode).After analysis, both the even and odd modes of thexandypolarization modes have energy leakage into the MF-filled holes,but more than 85%of the energy is still confined in the core.Since the refractive index of the MF is higher than that of air, it is more likely to cause energy leakage into the fiber cladding.Confinement loss can be used to measure the PCF’s ability to confine light.A perfectly matched layer is introduced into the numerical model to analyze the confinement loss of DCPCF.The results show that the confinement loss of the modes supporting transmission in the DCPCF of the entire communication band is less than 10-8dB/m,which is negligible for millimeter-level device.

    Fig.3.Mode field distribution of supermodes in DCPCF, at λ =1550 nm,Λ =5.25μm,d/Λ =0.44.

    3.2.Characteristics of the polarization beam splitter

    The PBS performance of DCPCF in the absence of magnetic field was analyzed.When the incident light is incident from the fiber core C1, the energy of the two polarization modes in the two fiber cores changes as the propagation distance increases, and finally only one polarization state of light can be obtained in each fiber core.According to formula(3),the normalized power of the two polarization states in C1 can be obtained as a function of the transmission distance,atΛ=5.25μm,d/Λ=0.44,λ=1550 nm.As shown in Fig.4, when incident at C1, the normalized powers of the two polarization modes at the output of C1 vary with different periods.When the length of the fiber is 38.638 mm, the main energy at the output end of C1 isypolarization mode,and the energy in thexpolarization mode almost does not exist, thus realizing the function of polarization beam splitting.The ER is-88.6 dB,which is another important indicator for PBS.The results show that the device has excellent splitting function when no external magnetic field is applied.

    Fig.4.The normalized power of x and y polarizations versus the transmission distance,at λ =1550 nm,Λ =5.25μm,d/Λ =0.44.The red dotted line is magnetic field magnitude at the first order separate length of 38.638 mm.

    3.3.Analysis of magnetic tuning characteristics

    After the dimensions of the device are determined, the factors that affectLxandLyare mainly wavelength and magnetic field magnitude.The propagation constantsβof the four supermodes in the DCPCF will change with the wavelength, and the coupling distance will also change according to Eq.(2), as shown in Fig.5(a).In the wavelength range of 1260 nm-1675 nm,LxandLydecrease as the wavelength increases.In addition, the refractive index of the MF increases with the increase of the magnetic field magnitude, so the change of the magnetic field magnitude will also cause the change of the propagation constant of the four supermodes in the DCPCF.As shown in Fig.5(b), when the magnetic field strength is lower than the critical value,the coupling distance of the two polarization modes does not change with the change of the magnetic field, but when the magnetic field magnitude exceeds the critical value,it will decrease with the increase of the magnetic field magnitude.Furthermore,the rate of change of the coupling distance decreases with increasing magnetic field magnitude.This is because the rate of change of the refractive index of the magnetic fluid decreases as the magnetic field magnitude increases, eventually reaching its saturation value.

    The normalized power and ER for thexandypolarization modes at the outports vary with wavelength.Combining Eqs.(3)and(5), the function curve of the normalized energy and ER of thexandypolarization modes of the output port of the fiber core C1 as a function of wavelength can be drawn without magnetic field, as shown in Fig.6.As described in the figure,the change of wavelength will cause the change of the energy and ER of the two polarization modes at the same side output port under the condition of fixed device structure.This is because a change in wavelength causes a change in the coupling lengthsLxandLy, which have been analyzed in detail above.As shown in Fig.6,at the wavelengths of 1514 nm,1550 nm,1585 nm,and 1618 nm,the normalized energy and ER are relatively high,and their values are(98.97%,37.5 dB),(99.87%,-88.7 dB),(99.72%,37.9 dB),(98.95%,-39.7 dB),respectively,the positive ER means that the main energy isxpolarization mode, and the negative ER means that the main energy is in they-polarization mode.The magnetic tuning characteristics of the device are analyzed below.

    Fig.5.Coupling length versus the magnitude of(a)wavelength and(b)magnetic intensity,at λ =1550 nm,Λ =5.25μm,d/Λ =0.44.

    Fig.6.The extinction ratio and normalized power versus the wavelength without magnitude of magnetic field, at z=38.638 mm, Λ =5.25μm,d/Λ =0.44.

    Fig.7.Normalized power and extinction ratio versus the magnitude of magnetic field at different wavelengths (a) 1514 nm, (b) 1550 nm,(c)1585 nm,(d)1618 nm.

    Figure 7 depicts the normalized power and ER for the two polarization modes of the C1 output port as a function of magnetic field amplitude.According to formula(1),it can be known that the refractive index of MF does not change when the magnetic fieldH&lt;44 Oe,which does not affect the change of the intensity of the two polarization modes of the outports.As shown in Fig.7,the intensity and ER of the two polarization modes at the output port are stable in the range of 0 Oe-43 Oe.When the magnetic field intensityH&gt;44 Oe,the refractive index of MF changes with the magnetic field intensity, which will affect the intensity of the two polarization modes of the output port.Figure 7(a) illustrates the normalized energy and ER at the wavelength of 1514 nm.The results show that the normalized energy and ER, at magnetic fields of 103 Oe, 184 Oe, 596 Oe, are (99.5%, 38.03 dB),(99.52%, 38.84 dB), (99.85%, 67.92 dB), respectively.At the wavelength of 1550nm,the results are shown in Fig.7(b).The normalized energy and ER of the two polarization modes,at magnetic fields of 101 Oe, 187 Oe, 618 Oe, are (99.98%,46.6 dB),(99.98%,54.04 dB),(99.89%,67.29 dB),separately.Figure 7(c)demonstrates the normalized energy and ER at the wavelength of 1585 nm.The results show that the normalized energy and ER, at magnetic fields of 100 Oe, 181 Oe,are (99.62%, 47.57 dB), (96.4%,49.46 dB), respectively.At the wavelength of 1618 nm,the results are shown in Fig.7(d).The normalized energy and ER of the two polarization modes,at magnetic fields of 108 Oe,191 Oe,are(98.07%,76.56 dB),(97.51%,-59.44 dB),separately.

    The positive ER means that the main energy is in thexpolarization mode, and the negative ER means that the main energy is in they-polarization mode.It can be concluded that the rotation of the polarization mode of the output port can be realized by changing the amplitude of the magnetic field,and the ER of the output polarization mode can also be tuned by the magnetic field.The device can achieve polarization mode rotation at the first-order polarization separation length.When the magnetic field magnitude is 103 Oe,101 Oe,100 Oe,and 108 Oe, respectively, the polarization rotation at wavelengths of 1514 nm,1550 nm,1585 nm,and 1618 nm can be achieved.In this case, figure 8 depicts the wavelength as a function of ER when outputtingxorypolarization modes for different wavelengths.Since the device needs to output different polarization modes at the same wavelength,the ER when outputting thexandypolarization modes separately should be more than 10 dB.In this case,when the central wavelengths are 1514 nm,1550 nm,1585 nm,and 1618 nm,the wavelength ranges that satisfy the above conditions are 1509 nm-1520 nm,1544 nm-1556 nm,1578 nm-1591 nm,and 1611 nm-1624 nm,respectively.

    The above analysis results are obtained at the temperature of 21°C.However, when the ambient temperature changes,the characteristics of the MF will also be affected.Thus, the influence of temperature on the device performance should be considered.The variation of temperature not only affects the refractive index of the MF, but also changes its critical magnetic field strengthHc.[35,38]In the literature,the influence of temperature on the refractive index andHcof water-based MF with a concentration of 0.85 emu/g was analyzed.The criticalHcwas found to increase from 22 Oe to 50 Oe upon increasing the temperature from 8.0°C to 60.0°C, and it can be seen thatHcand temperature are basically linear.In addition, they found that the refractive index of the MF increases linearly and continuously with the particle concentration under zero field, and the magnitude of thermo-optic coefficient dnMF/dTis less than 10-4°C-1.The concentration of the water-based MF used in this paper is 0.53 emu/g.[39]Therefore, its characteristics will have a similar variation trend as the MF with a concentration of 0.85 emu/g.The above analysis indicates that the alteration of temperature will change the magnetic field strength of the polarization switch for polarization conversion to a certain extent.

    Fig.8.Extinction ratio of x or y polarization mode versus the wavelength,at the centre wavelength of(a)1514 nm;(b)1550 nm;(c)1585 nm;(d)1618 nm.

    Fig.9.Extinction ratio of x or y polarization mode versus the intercepted fibre length,at different wavelengths of(a)1514 nm,(b)1550 nm,(c)1585 nm,(d)1618 nm.

    The fiber length of the polarization controller designed in this paper is 38.638 mm.At the central wavelength,as a PBS,the ER are 24.79 dB,-88.65 dB,37.89 dB,and 39.69 dB,respectively, and the ER when the polarization rotation occurs are-38.03 dB, 46.67 dB,-47.57 dB, and 76.56 dB.This length can be precisely controlled when analyzing its characteristics by numerical simulation.However, there is an error when cutting the fiber in the experiment.Figure 9 depicts the curves of the ER as a function of the cutoff length of the fiber when the wavelengths are 1514 nm, 1550 nm, 1585 nm, and 1618 nm, respectively.As shown in the figure, the absolute value of the extinction ratio between the lengths of 38.5 mm and 38.7 mm is greater than 20 dB, which indicates that the device has a good fault tolerance for the interception of the fibre length.

    4.Conclusion

    In conclusion, we have conducted a detailed study of a novel polarization switch based on DCPCF, whose cladding air holes are uniform in diameter.The DCPCF is infiltrated with MF in the four holes longitudinal symmetrically of each core.Firstly,the influence of structural parameters on the performance of the polarization controller is analyzed by the full vector finite element method.The results show that the working length of the controller is 38.638 mm.The operating bands with ER greater than 10 dB when outputtingxorypolarization modes include 1509 nm to 1520 nm,1544 nm to 1556 nm,1578 nm to 1591 nm, and 1611 nm to 1624 nm.The magnetic magnitudes for polarization rotation in these bands are 103 Oe,101 Oe,100 Oe,and 108 Oe,respectively.The ER of the outputxorypolarization is greater than 20 dB in a range of 38.5 mm to 38.7 mm for the proposed polarization controller.It shows that the device has a good tolerance for the fiber cut length.In addition, the polarization controller has the characteristics of simple manufacture and fair compatibility with traditional optical fiber communication and sensing systems.

    Acknowledgement

    Project supported by the National Key Research and Development Program of China “National Quality Infrastructure”(Grant No.2021YFF0600902).

    Data availability

    Data underlying the results presented in this paper may be available from the corresponding author upon reasonable request.

    猜你喜歡
    周軍
    《拓荒者》
    民族藝林(2022年4期)2023-01-19 04:06:16
    周軍:手握“最準的尺”,從“0”到“∞”篤信前行
    華人時刊(2021年15期)2021-11-27 09:16:24
    27年后與恩人再聚首
    新天地(2018年8期)2018-10-31 10:31:30
    “軸對稱”復習點睛
    對一道行程問題解法的反思
    對一個還原坐標系問題的反思
    對一道直角三角形相似課本習題的反思
    《整式的乘法與因式分解》專題復習
    月光唇印
    延河(2014年7期)2014-07-25 17:18:59
    水中幽靈
    小小說月刊(2011年6期)2011-05-14 14:55:08
    中文天堂在线官网| 亚洲人成网站在线观看播放| 97在线视频观看| 制服诱惑二区| 纵有疾风起免费观看全集完整版| 亚洲视频免费观看视频| 亚洲少妇的诱惑av| 久久久久久人妻| 日韩一区二区三区影片| 国产成人免费无遮挡视频| 2018国产大陆天天弄谢| 少妇人妻 视频| 欧美亚洲 丝袜 人妻 在线| 99国产综合亚洲精品| 亚洲精品一二三| 亚洲精品一区蜜桃| 成人二区视频| 日韩一区二区三区影片| 国产精品嫩草影院av在线观看| 最近中文字幕高清免费大全6| 久久久国产一区二区| 99久久中文字幕三级久久日本| av免费观看日本| 麻豆av在线久日| 成人免费观看视频高清| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 在线观看国产h片| 色吧在线观看| 日韩制服骚丝袜av| 26uuu在线亚洲综合色| 少妇被粗大猛烈的视频| 如日韩欧美国产精品一区二区三区| 母亲3免费完整高清在线观看 | 国产97色在线日韩免费| 免费播放大片免费观看视频在线观看| 久久99蜜桃精品久久| 老鸭窝网址在线观看| 99精国产麻豆久久婷婷| 人人妻人人澡人人爽人人夜夜| 久久久精品免费免费高清| 狠狠精品人妻久久久久久综合| 亚洲成人一二三区av| 中文字幕亚洲精品专区| 美女视频免费永久观看网站| 一个人免费看片子| 国产精品不卡视频一区二区| 黄色一级大片看看| 亚洲,欧美精品.| 熟女av电影| 热99久久久久精品小说推荐| 国产xxxxx性猛交| 国产精品久久久av美女十八| 日本欧美视频一区| 精品国产一区二区三区四区第35| 天天躁狠狠躁夜夜躁狠狠躁| 新久久久久国产一级毛片| 黑人巨大精品欧美一区二区蜜桃| 91在线精品国自产拍蜜月| 日韩三级伦理在线观看| 欧美激情高清一区二区三区 | 精品少妇黑人巨大在线播放| 国产又色又爽无遮挡免| 国产成人午夜福利电影在线观看| 考比视频在线观看| 久久这里有精品视频免费| 一区福利在线观看| 成人午夜精彩视频在线观看| 日韩一区二区视频免费看| 少妇被粗大的猛进出69影院| 国产精品国产av在线观看| 久热这里只有精品99| 午夜福利乱码中文字幕| 精品人妻在线不人妻| 午夜福利一区二区在线看| 久久毛片免费看一区二区三区| 热99久久久久精品小说推荐| 亚洲,欧美精品.| 九草在线视频观看| 最近中文字幕2019免费版| 熟女电影av网| 欧美xxⅹ黑人| 亚洲国产日韩一区二区| 在线天堂中文资源库| 啦啦啦啦在线视频资源| 亚洲精品成人av观看孕妇| 一区二区三区激情视频| 免费黄色在线免费观看| 成年人午夜在线观看视频| 一级a爱视频在线免费观看| 日韩一区二区三区影片| 亚洲国产欧美日韩在线播放| 欧美成人精品欧美一级黄| 韩国精品一区二区三区| videos熟女内射| 美女中出高潮动态图| 宅男免费午夜| 99re6热这里在线精品视频| 黄色一级大片看看| 肉色欧美久久久久久久蜜桃| 欧美精品人与动牲交sv欧美| 久久久久久久久久久久大奶| 日产精品乱码卡一卡2卡三| 国产成人精品久久二区二区91 | 久久久a久久爽久久v久久| 精品国产一区二区三区四区第35| h视频一区二区三区| 免费看不卡的av| 久久鲁丝午夜福利片| 国产av码专区亚洲av| 满18在线观看网站| 久久精品亚洲av国产电影网| 最新中文字幕久久久久| 成人国产av品久久久| 黄网站色视频无遮挡免费观看| 性色av一级| 日本欧美视频一区| 综合色丁香网| 日本av手机在线免费观看| 国产在视频线精品| 在现免费观看毛片| 丝袜在线中文字幕| 精品视频人人做人人爽| 国产在线免费精品| 男人操女人黄网站| 国产片内射在线| 五月天丁香电影| 国产野战对白在线观看| 国产精品二区激情视频| 亚洲av日韩在线播放| 99热网站在线观看| 人人妻人人添人人爽欧美一区卜| 在现免费观看毛片| 男人舔女人的私密视频| 久久毛片免费看一区二区三区| 久久精品国产亚洲av高清一级| 精品久久久久久电影网| 国产av国产精品国产| 少妇 在线观看| 日本91视频免费播放| 日韩熟女老妇一区二区性免费视频| 亚洲av欧美aⅴ国产| 国产一区二区在线观看av| 美女xxoo啪啪120秒动态图| 在线 av 中文字幕| 日韩一区二区三区影片| 极品少妇高潮喷水抽搐| 91精品三级在线观看| 五月开心婷婷网| 如日韩欧美国产精品一区二区三区| 777米奇影视久久| 美女福利国产在线| 久久精品久久久久久久性| 亚洲精品中文字幕在线视频| 午夜老司机福利剧场| 午夜激情久久久久久久| 久久热在线av| 嫩草影院入口| 一二三四在线观看免费中文在| 国产精品久久久久久久久免| 日本91视频免费播放| a级片在线免费高清观看视频| 女人久久www免费人成看片| 免费女性裸体啪啪无遮挡网站| 91国产中文字幕| √禁漫天堂资源中文www| 啦啦啦中文免费视频观看日本| 9热在线视频观看99| 18禁动态无遮挡网站| 日本欧美视频一区| 99re6热这里在线精品视频| 中文字幕人妻丝袜制服| 母亲3免费完整高清在线观看 | 亚洲欧洲日产国产| av国产久精品久网站免费入址| 三级国产精品片| 久久久久久久大尺度免费视频| 巨乳人妻的诱惑在线观看| 日本黄色日本黄色录像| 久久人妻熟女aⅴ| 街头女战士在线观看网站| 女人精品久久久久毛片| 成人18禁高潮啪啪吃奶动态图| 免费日韩欧美在线观看| 在线观看www视频免费| av在线app专区| 亚洲在久久综合| 亚洲国产精品一区二区三区在线| 亚洲av日韩在线播放| 亚洲国产av新网站| 黑人欧美特级aaaaaa片| 亚洲四区av| 国产欧美日韩一区二区三区在线| 久久这里有精品视频免费| 国产精品.久久久| 日产精品乱码卡一卡2卡三| 国产一区亚洲一区在线观看| 日日爽夜夜爽网站| 五月天丁香电影| 久久精品国产亚洲av天美| 午夜福利一区二区在线看| 亚洲欧美精品自产自拍| 午夜av观看不卡| 国产精品欧美亚洲77777| 色播在线永久视频| 免费看av在线观看网站| h视频一区二区三区| 亚洲人成77777在线视频| 中文欧美无线码| 欧美bdsm另类| 国产精品 欧美亚洲| 美国免费a级毛片| 午夜免费男女啪啪视频观看| 日韩 亚洲 欧美在线| 亚洲美女搞黄在线观看| 久久亚洲国产成人精品v| 狠狠婷婷综合久久久久久88av| 男人操女人黄网站| 亚洲国产精品一区三区| 最近2019中文字幕mv第一页| 91成人精品电影| 美女脱内裤让男人舔精品视频| 日韩精品有码人妻一区| 91在线精品国自产拍蜜月| 韩国高清视频一区二区三区| 国产av一区二区精品久久| 亚洲av电影在线进入| 国产国语露脸激情在线看| 综合色丁香网| 国产 精品1| 精品99又大又爽又粗少妇毛片| 欧美变态另类bdsm刘玥| 精品国产一区二区三区久久久樱花| 只有这里有精品99| 成年女人毛片免费观看观看9 | 国产一区二区在线观看av| 亚洲情色 制服丝袜| 男人爽女人下面视频在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲成色77777| 岛国毛片在线播放| 亚洲欧美一区二区三区久久| 久久99精品国语久久久| 丰满少妇做爰视频| 久久免费观看电影| 热99国产精品久久久久久7| 亚洲精品日本国产第一区| 伊人久久大香线蕉亚洲五| 99re6热这里在线精品视频| 欧美精品一区二区大全| 成人国产av品久久久| 久久久欧美国产精品| 电影成人av| 中文精品一卡2卡3卡4更新| 卡戴珊不雅视频在线播放| 成人二区视频| 黄片小视频在线播放| 男女免费视频国产| 成人手机av| 在线观看免费高清a一片| 免费黄频网站在线观看国产| 久久人妻熟女aⅴ| 国产亚洲午夜精品一区二区久久| 久久久久精品性色| 欧美变态另类bdsm刘玥| 麻豆av在线久日| 欧美中文综合在线视频| 水蜜桃什么品种好| 男人添女人高潮全过程视频| 超色免费av| 啦啦啦中文免费视频观看日本| 中文字幕亚洲精品专区| 90打野战视频偷拍视频| 中文字幕色久视频| 国产精品国产三级国产专区5o| 国产又色又爽无遮挡免| 两性夫妻黄色片| xxxhd国产人妻xxx| 九九爱精品视频在线观看| 黄色配什么色好看| 99香蕉大伊视频| 91精品国产国语对白视频| 九九爱精品视频在线观看| 成人免费观看视频高清| 久久精品国产自在天天线| 久久精品久久精品一区二区三区| 中文字幕亚洲精品专区| 看十八女毛片水多多多| 看非洲黑人一级黄片| 免费观看无遮挡的男女| 亚洲图色成人| 欧美激情高清一区二区三区 | av不卡在线播放| 亚洲少妇的诱惑av| 国产精品av久久久久免费| 一本—道久久a久久精品蜜桃钙片| 亚洲国产欧美网| 国产在线一区二区三区精| 色94色欧美一区二区| 亚洲av.av天堂| 一本久久精品| 亚洲国产毛片av蜜桃av| 亚洲第一区二区三区不卡| 国产精品 欧美亚洲| 欧美成人午夜精品| av有码第一页| 99热国产这里只有精品6| 99九九在线精品视频| 国产亚洲精品第一综合不卡| 香蕉丝袜av| 日本猛色少妇xxxxx猛交久久| 久久人人爽人人片av| 国产精品二区激情视频| 久久精品久久久久久噜噜老黄| 春色校园在线视频观看| 巨乳人妻的诱惑在线观看| 只有这里有精品99| 日韩欧美精品免费久久| 少妇精品久久久久久久| 一级毛片黄色毛片免费观看视频| 日韩三级伦理在线观看| www日本在线高清视频| 国产午夜精品一二区理论片| 婷婷色麻豆天堂久久| 黄色怎么调成土黄色| 亚洲国产成人一精品久久久| 伦理电影大哥的女人| 日韩一卡2卡3卡4卡2021年| 大片免费播放器 马上看| 久久精品国产亚洲av天美| 午夜福利在线观看免费完整高清在| 大片电影免费在线观看免费| a级毛片黄视频| 18禁国产床啪视频网站| 欧美日韩精品成人综合77777| 色播在线永久视频| 国产一区有黄有色的免费视频| 黄片无遮挡物在线观看| 最近中文字幕高清免费大全6| 亚洲精品一区蜜桃| 久久人妻熟女aⅴ| 欧美国产精品一级二级三级| 久久精品久久久久久久性| 国产一区二区激情短视频 | 2018国产大陆天天弄谢| 人妻系列 视频| 美女高潮到喷水免费观看| 色吧在线观看| 99热全是精品| 波野结衣二区三区在线| 国产精品亚洲av一区麻豆 | 久久久久久人人人人人| 免费观看在线日韩| 日本爱情动作片www.在线观看| 久久精品国产a三级三级三级| 欧美精品一区二区大全| 日韩成人av中文字幕在线观看| 欧美97在线视频| 观看av在线不卡| 女性被躁到高潮视频| 久久久久精品人妻al黑| 十分钟在线观看高清视频www| 少妇熟女欧美另类| 国产成人精品福利久久| 中文精品一卡2卡3卡4更新| 国产精品久久久久久精品古装| 欧美激情极品国产一区二区三区| 亚洲av男天堂| 国产 精品1| 夫妻性生交免费视频一级片| 中文字幕色久视频| 麻豆精品久久久久久蜜桃| 波野结衣二区三区在线| 国产男女内射视频| 亚洲国产成人一精品久久久| 欧美xxⅹ黑人| 国产精品女同一区二区软件| 韩国精品一区二区三区| 纵有疾风起免费观看全集完整版| 视频区图区小说| 日本av手机在线免费观看| 国产精品99久久99久久久不卡 | 毛片一级片免费看久久久久| 啦啦啦啦在线视频资源| 在线观看免费视频网站a站| 国产熟女欧美一区二区| 久久精品国产自在天天线| xxx大片免费视频| 美女福利国产在线| 99热全是精品| 国产色婷婷99| 国产av码专区亚洲av| 2022亚洲国产成人精品| 亚洲成国产人片在线观看| 热re99久久国产66热| 亚洲精品中文字幕在线视频| 一区二区日韩欧美中文字幕| 成年美女黄网站色视频大全免费| 爱豆传媒免费全集在线观看| 天堂俺去俺来也www色官网| 欧美日韩综合久久久久久| 免费久久久久久久精品成人欧美视频| 黄色 视频免费看| 中文字幕另类日韩欧美亚洲嫩草| 欧美少妇被猛烈插入视频| 麻豆精品久久久久久蜜桃| 99精国产麻豆久久婷婷| 国产亚洲最大av| 国产亚洲av片在线观看秒播厂| 亚洲精品,欧美精品| 欧美日韩国产mv在线观看视频| 国产免费又黄又爽又色| 国产熟女午夜一区二区三区| 在现免费观看毛片| 美女高潮到喷水免费观看| 有码 亚洲区| 久久久久久人人人人人| 精品福利永久在线观看| 王馨瑶露胸无遮挡在线观看| kizo精华| 自拍欧美九色日韩亚洲蝌蚪91| 热re99久久精品国产66热6| 国产成人精品在线电影| 夫妻午夜视频| 91成人精品电影| 欧美日韩精品成人综合77777| 9热在线视频观看99| 最近中文字幕2019免费版| 黄色怎么调成土黄色| 国产精品 欧美亚洲| 韩国av在线不卡| 亚洲经典国产精华液单| 精品一区在线观看国产| 久久久a久久爽久久v久久| 中文字幕人妻熟女乱码| 国产精品99久久99久久久不卡 | 老司机影院毛片| 久久综合国产亚洲精品| 久热这里只有精品99| 一本色道久久久久久精品综合| 精品国产超薄肉色丝袜足j| 欧美xxⅹ黑人| 美女福利国产在线| 欧美日本中文国产一区发布| 色哟哟·www| 亚洲少妇的诱惑av| av免费观看日本| 九色亚洲精品在线播放| 人体艺术视频欧美日本| 中文字幕av电影在线播放| 欧美+日韩+精品| 91久久精品国产一区二区三区| 国产一区有黄有色的免费视频| 国产片特级美女逼逼视频| 国产av国产精品国产| 欧美精品高潮呻吟av久久| 男人操女人黄网站| 人妻系列 视频| 18在线观看网站| 亚洲欧美精品综合一区二区三区 | 婷婷色综合大香蕉| 午夜福利影视在线免费观看| 久久久久国产网址| 成人漫画全彩无遮挡| 男人添女人高潮全过程视频| 亚洲国产av影院在线观看| 国产一级毛片在线| 久久精品亚洲av国产电影网| 狂野欧美激情性bbbbbb| 免费女性裸体啪啪无遮挡网站| 精品少妇久久久久久888优播| 午夜福利一区二区在线看| 少妇的丰满在线观看| 乱人伦中国视频| 中文欧美无线码| 亚洲精品第二区| 欧美精品人与动牲交sv欧美| 国产免费一区二区三区四区乱码| 宅男免费午夜| 日韩一卡2卡3卡4卡2021年| 亚洲精品国产av蜜桃| 亚洲av中文av极速乱| 不卡av一区二区三区| 熟女电影av网| 在线观看免费高清a一片| 91精品三级在线观看| 欧美国产精品一级二级三级| av在线老鸭窝| 日本wwww免费看| 丝袜在线中文字幕| 久久99精品国语久久久| 热99国产精品久久久久久7| 免费高清在线观看日韩| 丰满迷人的少妇在线观看| 久久毛片免费看一区二区三区| 精品午夜福利在线看| 亚洲少妇的诱惑av| 精品国产乱码久久久久久小说| 亚洲成人一二三区av| 这个男人来自地球电影免费观看 | 国产精品熟女久久久久浪| 亚洲av电影在线进入| 在线观看免费高清a一片| 久久久久精品人妻al黑| 美女xxoo啪啪120秒动态图| 日韩欧美一区视频在线观看| 中文字幕色久视频| 久久精品国产亚洲av天美| 欧美人与性动交α欧美精品济南到 | 高清av免费在线| 国产免费又黄又爽又色| 精品一区二区免费观看| 国产欧美日韩一区二区三区在线| 亚洲天堂av无毛| 少妇熟女欧美另类| 麻豆av在线久日| 女性被躁到高潮视频| 国产高清不卡午夜福利| 国产免费视频播放在线视频| 成年人免费黄色播放视频| 日韩成人av中文字幕在线观看| 国产成人一区二区在线| 国产亚洲最大av| 日韩欧美一区视频在线观看| 免费久久久久久久精品成人欧美视频| 久久久久久久亚洲中文字幕| 精品久久久久久电影网| 免费黄网站久久成人精品| 欧美日韩视频精品一区| 久久久精品区二区三区| 超色免费av| 少妇 在线观看| 成人毛片a级毛片在线播放| 久久久久网色| 肉色欧美久久久久久久蜜桃| 亚洲欧美清纯卡通| 久久国产精品大桥未久av| 精品国产国语对白av| av有码第一页| 男男h啪啪无遮挡| 亚洲中文av在线| 午夜福利视频精品| 伦理电影大哥的女人| 免费观看无遮挡的男女| 日韩,欧美,国产一区二区三区| 青青草视频在线视频观看| 欧美激情 高清一区二区三区| 国产高清国产精品国产三级| 赤兔流量卡办理| 国产亚洲最大av| 蜜桃在线观看..| 制服人妻中文乱码| a级片在线免费高清观看视频| av天堂久久9| 欧美老熟妇乱子伦牲交| 两个人看的免费小视频| 熟女av电影| 人人澡人人妻人| 国产在线视频一区二区| 国产av国产精品国产| 国产精品.久久久| 亚洲国产欧美网| 99国产综合亚洲精品| 午夜老司机福利剧场| 欧美最新免费一区二区三区| 最近手机中文字幕大全| 天堂中文最新版在线下载| 亚洲精品日本国产第一区| 久久免费观看电影| 少妇 在线观看| 国产精品无大码| 中文字幕制服av| av片东京热男人的天堂| 午夜影院在线不卡| 亚洲欧洲日产国产| 国产精品国产av在线观看| 天堂8中文在线网| videosex国产| 一区二区三区乱码不卡18| 最近中文字幕2019免费版| 黄色配什么色好看| 9191精品国产免费久久| 久久久久精品人妻al黑| 少妇猛男粗大的猛烈进出视频| 久久久久久久久久人人人人人人| 丝袜美足系列| 中文字幕另类日韩欧美亚洲嫩草| 欧美国产精品va在线观看不卡| 亚洲一码二码三码区别大吗| 成年人午夜在线观看视频| 日韩人妻精品一区2区三区| 麻豆精品久久久久久蜜桃| 国产精品一国产av| 久久精品国产鲁丝片午夜精品| 中国三级夫妇交换| 亚洲三级黄色毛片| 有码 亚洲区| 亚洲精华国产精华液的使用体验| 校园人妻丝袜中文字幕| 久久精品人人爽人人爽视色| 午夜老司机福利剧场| 亚洲国产精品一区二区三区在线| 黄色视频在线播放观看不卡| 搡老乐熟女国产| tube8黄色片| 男女边吃奶边做爰视频| 久久久精品免费免费高清| 黄片播放在线免费| 国产白丝娇喘喷水9色精品| 美国免费a级毛片| 母亲3免费完整高清在线观看 | 欧美亚洲日本最大视频资源| 青草久久国产| 亚洲精品在线美女| 日韩不卡一区二区三区视频在线| 91成人精品电影|