• <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
    中文乱码字字幕精品一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | tube8黄色片| 联通29元200g的流量卡| 亚洲经典国产精华液单| 成人国产麻豆网| 91久久精品国产一区二区三区| 老司机影院成人| 麻豆乱淫一区二区| 国产欧美日韩精品一区二区| 日韩一本色道免费dvd| 中文精品一卡2卡3卡4更新| 精品少妇久久久久久888优播| 久久久久网色| 老女人水多毛片| 一个人看视频在线观看www免费| 免费av不卡在线播放| 少妇被粗大猛烈的视频| 国产成人a∨麻豆精品| 亚洲精品成人av观看孕妇| 成人一区二区视频在线观看| 中文字幕av成人在线电影| 久久久国产一区二区| 免费观看的影片在线观看| 如何舔出高潮| 高清视频免费观看一区二区| 婷婷色综合大香蕉| 欧美zozozo另类| 国产精品伦人一区二区| 国产毛片a区久久久久| 久久国产乱子免费精品| 色视频在线一区二区三区| a级毛色黄片| 午夜视频国产福利| 日韩一本色道免费dvd| 精品久久国产蜜桃| 国产美女午夜福利| 97超碰精品成人国产| 在线观看av片永久免费下载| 精品一区在线观看国产| 又爽又黄无遮挡网站| 久久亚洲国产成人精品v| 婷婷色麻豆天堂久久| 一级毛片电影观看| 久久精品久久久久久久性| 亚洲精华国产精华液的使用体验| 深夜a级毛片| 久久人人爽av亚洲精品天堂 | av在线亚洲专区| 亚洲精品视频女| 国产亚洲午夜精品一区二区久久 | 2021天堂中文幕一二区在线观| 欧美另类一区| 国产69精品久久久久777片| 中国美白少妇内射xxxbb| 午夜福利在线在线| 校园人妻丝袜中文字幕| 99久久人妻综合| av福利片在线观看| av又黄又爽大尺度在线免费看| 国产在视频线精品| 国产成人精品一,二区| 永久网站在线| 综合色丁香网| 亚洲成人久久爱视频| 午夜免费观看性视频| 亚洲av欧美aⅴ国产| 九色成人免费人妻av| 亚洲国产精品成人久久小说| 国产毛片在线视频| 亚洲欧美成人综合另类久久久| 亚洲国产精品成人综合色| 男人狂女人下面高潮的视频| 成人亚洲精品一区在线观看 | 久久人人爽人人爽人人片va| 韩国av在线不卡| 日日摸夜夜添夜夜添av毛片| 特大巨黑吊av在线直播| 精品国产一区二区三区久久久樱花 | 欧美激情国产日韩精品一区| 中文欧美无线码| 视频中文字幕在线观看| 国产91av在线免费观看| 免费av不卡在线播放| 色吧在线观看| av女优亚洲男人天堂| 男女无遮挡免费网站观看| 纵有疾风起免费观看全集完整版| av一本久久久久| 久久精品国产亚洲网站| 久久久久久久久久成人| 精品人妻视频免费看| 91在线精品国自产拍蜜月| 欧美 日韩 精品 国产| 九草在线视频观看| 国产伦精品一区二区三区视频9| 亚洲婷婷狠狠爱综合网| 成人欧美大片| 久久久午夜欧美精品| 欧美性猛交╳xxx乱大交人| 97精品久久久久久久久久精品| 看黄色毛片网站| 只有这里有精品99| 午夜福利视频精品| 国产探花在线观看一区二区| 街头女战士在线观看网站| 中国国产av一级| 欧美3d第一页| 大陆偷拍与自拍| 成人国产av品久久久| 亚洲精品aⅴ在线观看| 亚洲内射少妇av| 久久久久国产精品人妻一区二区| 特级一级黄色大片| 久久久久久久久大av| 热99国产精品久久久久久7| 免费人成在线观看视频色| 男人舔奶头视频| 一级毛片久久久久久久久女| 伊人久久国产一区二区| 三级国产精品欧美在线观看| 深爱激情五月婷婷| 极品少妇高潮喷水抽搐| 日韩一本色道免费dvd| 在线观看美女被高潮喷水网站| 18禁动态无遮挡网站| 看十八女毛片水多多多| 少妇的逼好多水| 狂野欧美白嫩少妇大欣赏| 午夜视频国产福利| 国产精品久久久久久久电影| 最新中文字幕久久久久| 久久久精品欧美日韩精品| 亚洲av免费高清在线观看| 亚洲精品中文字幕在线视频 | 黄色一级大片看看| 亚洲精品视频女| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品无大码| 精品久久久久久久久av| 欧美成人一区二区免费高清观看| 久久99热6这里只有精品| 一级毛片aaaaaa免费看小| 国产老妇女一区| 婷婷色综合大香蕉| 国产精品精品国产色婷婷| 男人爽女人下面视频在线观看| 欧美高清性xxxxhd video| 特级一级黄色大片| 亚洲精品久久久久久婷婷小说| 国产一区亚洲一区在线观看| videossex国产| 中文字幕亚洲精品专区| 秋霞伦理黄片| 免费看光身美女| 国产在线一区二区三区精| www.av在线官网国产| 少妇被粗大猛烈的视频| 久久99热这里只频精品6学生| 青青草视频在线视频观看| 亚洲va在线va天堂va国产| 国产欧美日韩一区二区三区在线 | 精品国产露脸久久av麻豆| 一级毛片 在线播放| 插阴视频在线观看视频| 最近最新中文字幕大全电影3| 亚洲国产高清在线一区二区三| 国产一区二区在线观看日韩| 日本av手机在线免费观看| 国产精品一二三区在线看| 免费av观看视频| 天天躁夜夜躁狠狠久久av| av又黄又爽大尺度在线免费看| 日韩成人av中文字幕在线观看| 亚洲国产色片| 中文在线观看免费www的网站| 国产精品一二三区在线看| 春色校园在线视频观看| 97精品久久久久久久久久精品| 国产乱来视频区| 乱码一卡2卡4卡精品| 少妇高潮的动态图| 国产精品偷伦视频观看了| 97在线视频观看| 日韩亚洲欧美综合| 新久久久久国产一级毛片| 人人妻人人澡人人爽人人夜夜| 日本与韩国留学比较| 国产亚洲最大av| 亚洲成人久久爱视频| 成年女人在线观看亚洲视频 | 亚洲精品国产色婷婷电影| 99热6这里只有精品| 国内揄拍国产精品人妻在线| 有码 亚洲区| 欧美xxxx性猛交bbbb| 五月开心婷婷网| 97超视频在线观看视频| 少妇熟女欧美另类| 一个人看的www免费观看视频| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久精品古装| 日本免费在线观看一区| 国产综合精华液| 中文字幕人妻熟人妻熟丝袜美| 久久久亚洲精品成人影院| 女人十人毛片免费观看3o分钟| 成人美女网站在线观看视频| 视频区图区小说| 国产精品一区www在线观看| 日韩国内少妇激情av| 美女国产视频在线观看| 美女视频免费永久观看网站| 成人毛片a级毛片在线播放| 欧美一区二区亚洲| 一级毛片aaaaaa免费看小| 黄片wwwwww| 丝袜脚勾引网站| 看十八女毛片水多多多| 少妇裸体淫交视频免费看高清| 国产高清三级在线| 亚洲精品日韩av片在线观看| 国产黄色视频一区二区在线观看| 你懂的网址亚洲精品在线观看| 美女xxoo啪啪120秒动态图| 中文在线观看免费www的网站| 欧美亚洲 丝袜 人妻 在线| 午夜老司机福利剧场| 尾随美女入室| 日韩av免费高清视频| 女的被弄到高潮叫床怎么办| 成人毛片60女人毛片免费| 超碰av人人做人人爽久久| 成人美女网站在线观看视频| 亚洲国产精品专区欧美| 一级av片app| 亚洲av日韩在线播放| 亚洲欧美日韩卡通动漫| 国产伦在线观看视频一区| 各种免费的搞黄视频| 一级av片app| 成人毛片60女人毛片免费| 丝袜美腿在线中文| 亚洲无线观看免费| 国产亚洲午夜精品一区二区久久 | 国产精品三级大全| 国产一区有黄有色的免费视频| 欧美日韩在线观看h| 99久国产av精品国产电影| 亚洲欧美日韩东京热| 国产av不卡久久| 午夜精品国产一区二区电影 | 美女cb高潮喷水在线观看| 天天一区二区日本电影三级| 99九九线精品视频在线观看视频| 免费人成在线观看视频色| 亚洲怡红院男人天堂| 国产视频内射| 亚洲av中文av极速乱| 亚洲精品影视一区二区三区av| 97超碰精品成人国产| av又黄又爽大尺度在线免费看| 国产欧美日韩一区二区三区在线 | 99久久人妻综合| 黄片无遮挡物在线观看| 少妇熟女欧美另类| 男女边摸边吃奶| 高清欧美精品videossex| 日韩国内少妇激情av| 97在线视频观看| 丝袜喷水一区| 亚洲av.av天堂| 极品教师在线视频| 国产精品偷伦视频观看了| 中文乱码字字幕精品一区二区三区| 成人亚洲精品一区在线观看 | 日本一二三区视频观看| 亚洲经典国产精华液单| 欧美激情在线99| 97超碰精品成人国产| 国产男人的电影天堂91| 伊人久久国产一区二区| 亚洲怡红院男人天堂| 嘟嘟电影网在线观看| 亚洲精品亚洲一区二区| 水蜜桃什么品种好| 免费观看在线日韩| 男插女下体视频免费在线播放| 亚洲,欧美,日韩| 亚洲av.av天堂| 国产成人aa在线观看| av在线亚洲专区| 亚洲美女视频黄频| 久久国内精品自在自线图片| 久久精品熟女亚洲av麻豆精品| av在线观看视频网站免费| 男男h啪啪无遮挡| 男人添女人高潮全过程视频| 亚洲成人中文字幕在线播放| 欧美性感艳星| av黄色大香蕉| av网站免费在线观看视频| 精品一区在线观看国产| 熟妇人妻不卡中文字幕| 国产高清不卡午夜福利| 午夜激情久久久久久久| 精品人妻熟女av久视频| 免费av观看视频| 97人妻精品一区二区三区麻豆| 欧美日韩国产mv在线观看视频 | 国产一区有黄有色的免费视频| 三级国产精品欧美在线观看| 麻豆国产97在线/欧美| 青春草国产在线视频| 看免费成人av毛片| 亚洲国产精品专区欧美| 禁无遮挡网站| 国产色婷婷99| 成人毛片a级毛片在线播放| 中国国产av一级| 熟妇人妻不卡中文字幕| 九九在线视频观看精品| 久久久久久久午夜电影| 亚洲精品国产av蜜桃| 国产男女超爽视频在线观看| 国产成人福利小说| 天堂俺去俺来也www色官网| 色5月婷婷丁香| 99热国产这里只有精品6| videossex国产| 联通29元200g的流量卡| 亚洲精品,欧美精品| 欧美日韩国产mv在线观看视频 | 国产午夜精品一二区理论片| 成人毛片a级毛片在线播放| 亚洲aⅴ乱码一区二区在线播放| 国产伦精品一区二区三区四那| 久久久久精品久久久久真实原创| 国产人妻一区二区三区在| 白带黄色成豆腐渣| a级毛片免费高清观看在线播放| 久久久久久久国产电影| 真实男女啪啪啪动态图| 日本色播在线视频| 久久精品久久久久久噜噜老黄| 大香蕉97超碰在线| 国产色爽女视频免费观看| eeuss影院久久| 人妻少妇偷人精品九色| 波野结衣二区三区在线| 中文精品一卡2卡3卡4更新| 亚洲精品成人av观看孕妇| 99视频精品全部免费 在线| 久久久精品欧美日韩精品| 香蕉精品网在线| 色网站视频免费| 国产亚洲最大av| 亚洲美女搞黄在线观看| 汤姆久久久久久久影院中文字幕| 丰满少妇做爰视频| 中国美白少妇内射xxxbb| 精品久久久久久久末码| 在线播放无遮挡| 国产一区亚洲一区在线观看| 国产v大片淫在线免费观看| 1000部很黄的大片| 九色成人免费人妻av| 啦啦啦啦在线视频资源| 精品酒店卫生间| 丰满乱子伦码专区| 新久久久久国产一级毛片| 午夜福利在线在线| 2021少妇久久久久久久久久久| 搞女人的毛片| 一级a做视频免费观看| 九九爱精品视频在线观看| 国产亚洲5aaaaa淫片| 色播亚洲综合网| 久热这里只有精品99| av免费在线看不卡| 又爽又黄无遮挡网站| 高清午夜精品一区二区三区| 国产精品精品国产色婷婷| 亚洲在线观看片| 久久久久精品久久久久真实原创| 少妇 在线观看| 人体艺术视频欧美日本| 青青草视频在线视频观看| 秋霞伦理黄片| 亚洲国产av新网站| 禁无遮挡网站| 国产精品不卡视频一区二区| 精品国产三级普通话版| 男人和女人高潮做爰伦理| 午夜免费观看性视频| 少妇丰满av| 男人爽女人下面视频在线观看| 亚洲色图av天堂| 色视频在线一区二区三区| 美女被艹到高潮喷水动态| 日韩精品有码人妻一区| 精品久久久精品久久久| 99热这里只有是精品50| 亚洲精品视频女| 日日摸夜夜添夜夜爱| 国产久久久一区二区三区| 国产午夜福利久久久久久| 国产精品精品国产色婷婷| 夜夜爽夜夜爽视频| 少妇人妻一区二区三区视频| 免费观看性生交大片5| 大又大粗又爽又黄少妇毛片口| 黄片wwwwww| 成年人午夜在线观看视频| 免费av不卡在线播放| 啦啦啦中文免费视频观看日本| 免费av观看视频| 大片电影免费在线观看免费| 蜜桃久久精品国产亚洲av| 国产黄频视频在线观看| 国产片特级美女逼逼视频| 三级经典国产精品| 我要看日韩黄色一级片| 中文字幕久久专区| 成人国产麻豆网| 超碰97精品在线观看| 中文天堂在线官网| 亚洲欧美清纯卡通| 国产成人精品福利久久| 国产精品女同一区二区软件| 日韩视频在线欧美| 高清日韩中文字幕在线| 国产69精品久久久久777片| 大话2 男鬼变身卡| 丰满少妇做爰视频| 久久久久国产网址| 又粗又硬又长又爽又黄的视频| 欧美潮喷喷水| 少妇被粗大猛烈的视频| 国产美女午夜福利| 一个人看的www免费观看视频| 嫩草影院新地址| 高清视频免费观看一区二区| 亚洲最大成人av| 国产日韩欧美在线精品| 亚洲综合色惰| 干丝袜人妻中文字幕| 国产永久视频网站| 国产在视频线精品| 国产又色又爽无遮挡免| 欧美精品国产亚洲| 97在线人人人人妻| 国产一区二区三区av在线| 久久精品国产亚洲网站| 午夜福利网站1000一区二区三区| 亚洲精品国产av成人精品| 国产成人精品福利久久| 中国国产av一级| 国产色爽女视频免费观看| 免费观看a级毛片全部| 少妇熟女欧美另类| 欧美xxⅹ黑人| 男男h啪啪无遮挡| 在线精品无人区一区二区三 | 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区乱码不卡18| 秋霞在线观看毛片| 中文天堂在线官网| 国产成人a∨麻豆精品| 99re6热这里在线精品视频| 欧美另类一区| 国产精品麻豆人妻色哟哟久久| 亚洲欧美日韩无卡精品| 免费黄色在线免费观看| 亚洲性久久影院| 国产精品三级大全| 涩涩av久久男人的天堂| 国模一区二区三区四区视频| 99re6热这里在线精品视频| 日韩三级伦理在线观看| 亚洲成人一二三区av| 午夜福利在线观看免费完整高清在| 亚洲,一卡二卡三卡| 亚洲人成网站在线观看播放| av免费在线看不卡| 18禁动态无遮挡网站| 在线观看av片永久免费下载| 一本久久精品| 亚洲av二区三区四区| 18禁在线播放成人免费| 午夜免费男女啪啪视频观看| 18禁裸乳无遮挡动漫免费视频 | 国产毛片a区久久久久| 一级毛片电影观看| 精品一区二区三卡| 亚洲美女视频黄频| 久久久久久久大尺度免费视频| 日日啪夜夜撸| 在线 av 中文字幕| 免费观看av网站的网址| 国产成人a∨麻豆精品| 久久久久久久久久成人| 精品久久国产蜜桃| 国产成人a区在线观看| 久久久久精品久久久久真实原创| 性色av一级| 2022亚洲国产成人精品| 99热6这里只有精品| 午夜福利高清视频| 男女国产视频网站| 午夜爱爱视频在线播放| 人妻夜夜爽99麻豆av| 国产亚洲av嫩草精品影院| 久久女婷五月综合色啪小说 | 在线免费十八禁| 白带黄色成豆腐渣| av国产精品久久久久影院| 久久久欧美国产精品| 亚洲欧洲日产国产| 一级爰片在线观看| 国产精品av视频在线免费观看| 免费在线观看成人毛片| 中文字幕av成人在线电影| 国产精品熟女久久久久浪| 日本-黄色视频高清免费观看| 91精品伊人久久大香线蕉| 免费看光身美女| 亚洲av电影在线观看一区二区三区 | 欧美高清成人免费视频www| 久久ye,这里只有精品| 国产成人午夜福利电影在线观看| 啦啦啦在线观看免费高清www| av黄色大香蕉| 18禁裸乳无遮挡动漫免费视频 | 看非洲黑人一级黄片| 婷婷色av中文字幕| 欧美日韩视频精品一区| 久久精品夜色国产| 欧美精品人与动牲交sv欧美| 美女cb高潮喷水在线观看| 国产精品一区二区在线观看99| 亚洲aⅴ乱码一区二区在线播放| av在线亚洲专区| 亚洲精品国产成人久久av| kizo精华| a级毛片免费高清观看在线播放| 精品久久久噜噜| 国产精品国产三级国产av玫瑰| 又爽又黄a免费视频| 久久久久久伊人网av| 久久久久久久亚洲中文字幕| 亚洲av男天堂| 亚洲精品乱码久久久v下载方式| 草草在线视频免费看| 婷婷色av中文字幕| 国产黄频视频在线观看| 午夜激情福利司机影院| 成人毛片a级毛片在线播放| 人妻少妇偷人精品九色| 国内少妇人妻偷人精品xxx网站| 久久精品夜色国产| 久久精品人妻少妇| 最近的中文字幕免费完整| 国产永久视频网站| 精品人妻熟女av久视频| 精品人妻偷拍中文字幕| 亚洲精品成人av观看孕妇| 亚洲av国产av综合av卡| 亚洲国产精品成人久久小说| 在线看a的网站| 男人和女人高潮做爰伦理| 欧美xxⅹ黑人| 精品一区二区三区视频在线| 国产免费一级a男人的天堂| 在现免费观看毛片| 久久热精品热| 在线免费十八禁| 亚洲精品亚洲一区二区| 免费黄色在线免费观看| 久久久午夜欧美精品| 少妇 在线观看| 黄片无遮挡物在线观看| 嘟嘟电影网在线观看| 久久精品久久久久久噜噜老黄| 亚洲美女搞黄在线观看| 草草在线视频免费看| 国产亚洲av片在线观看秒播厂| 亚洲欧洲国产日韩| 欧美高清性xxxxhd video| 免费电影在线观看免费观看| 在线观看人妻少妇| 日韩制服骚丝袜av| 久久精品国产鲁丝片午夜精品| 少妇人妻 视频| 97在线人人人人妻| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲av中文字字幕乱码综合| 国产视频首页在线观看| 国产精品人妻久久久久久| 日韩成人av中文字幕在线观看| 黄色配什么色好看| 精品久久久久久久末码| 精品国产一区二区三区久久久樱花 | 国产精品一区二区三区四区免费观看| 久久韩国三级中文字幕| 在线a可以看的网站| 欧美成人a在线观看| 免费观看在线日韩| 国产在视频线精品| 两个人的视频大全免费| 亚洲精品第二区| 免费观看a级毛片全部| av线在线观看网站| 黑人高潮一二区| 国产午夜精品一二区理论片|