• <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
    国内精品美女久久久久久| 欧美绝顶高潮抽搐喷水| 午夜久久久久精精品| 欧美3d第一页| 国产av一区在线观看免费| 中文资源天堂在线| 国产成+人综合+亚洲专区| 夜夜躁狠狠躁天天躁| 亚洲国产精品久久男人天堂| 免费av观看视频| 麻豆一二三区av精品| 亚洲一区二区三区不卡视频| 99久久99久久久精品蜜桃| 国产熟女xx| 国产精品野战在线观看| 九色国产91popny在线| 国产av一区在线观看免费| 国产久久久一区二区三区| 欧美国产日韩亚洲一区| 直男gayav资源| 在线播放国产精品三级| 久久人妻av系列| 国产探花在线观看一区二区| 欧美潮喷喷水| 又爽又黄a免费视频| 国产在视频线在精品| 国产成人a区在线观看| 麻豆一二三区av精品| 婷婷亚洲欧美| 国产免费一级a男人的天堂| 人人妻人人澡欧美一区二区| 高清日韩中文字幕在线| 国产精品自产拍在线观看55亚洲| 午夜福利在线观看吧| 久久久久性生活片| 免费在线观看亚洲国产| 久久精品国产亚洲av天美| 日日摸夜夜添夜夜添小说| 国产伦在线观看视频一区| 日韩中文字幕欧美一区二区| 午夜福利高清视频| 久久精品国产自在天天线| 欧美精品国产亚洲| 国产一区二区三区在线臀色熟女| 性插视频无遮挡在线免费观看| 精品欧美国产一区二区三| 老女人水多毛片| 成年女人毛片免费观看观看9| 性插视频无遮挡在线免费观看| 午夜亚洲福利在线播放| 99久久无色码亚洲精品果冻| 中亚洲国语对白在线视频| 午夜免费激情av| 男人和女人高潮做爰伦理| 久久久国产成人免费| 一级黄片播放器| 日本黄色片子视频| 老司机福利观看| 亚洲天堂国产精品一区在线| 国产精品一区二区免费欧美| 国产探花在线观看一区二区| 天堂av国产一区二区熟女人妻| 久久久久免费精品人妻一区二区| 老熟妇乱子伦视频在线观看| 精品不卡国产一区二区三区| 精品一区二区三区人妻视频| 美女高潮喷水抽搐中文字幕| 十八禁人妻一区二区| 尤物成人国产欧美一区二区三区| 精品福利观看| 丰满人妻一区二区三区视频av| 日本 欧美在线| 中文字幕免费在线视频6| 99国产精品一区二区三区| 自拍偷自拍亚洲精品老妇| 亚洲av电影在线进入| 国产精品嫩草影院av在线观看 | 亚洲国产欧美人成| 精品国产三级普通话版| 国产高清激情床上av| 久久亚洲精品不卡| 在线观看一区二区三区| 丰满的人妻完整版| 日韩欧美三级三区| 欧美性感艳星| 国模一区二区三区四区视频| 成人性生交大片免费视频hd| 女人被狂操c到高潮| 黄色丝袜av网址大全| 99热这里只有是精品50| 成人高潮视频无遮挡免费网站| 白带黄色成豆腐渣| 全区人妻精品视频| 丰满人妻一区二区三区视频av| 国产亚洲精品久久久com| 真实男女啪啪啪动态图| 中文字幕av成人在线电影| 国产一区二区激情短视频| 欧美xxxx性猛交bbbb| 成人精品一区二区免费| 精华霜和精华液先用哪个| 岛国在线免费视频观看| 色尼玛亚洲综合影院| 一级黄色大片毛片| 丰满人妻一区二区三区视频av| 午夜精品在线福利| 99热这里只有是精品50| 欧美精品啪啪一区二区三区| 日本熟妇午夜| 久久精品人妻少妇| 日韩中字成人| 国产精品久久视频播放| 丁香六月欧美| 久久久久亚洲av毛片大全| 乱人视频在线观看| 特大巨黑吊av在线直播| 成人特级黄色片久久久久久久| 精品午夜福利在线看| 一级作爱视频免费观看| 91九色精品人成在线观看| 久久国产乱子伦精品免费另类| 欧美成人性av电影在线观看| 十八禁国产超污无遮挡网站| 欧美+亚洲+日韩+国产| 久久久久久久午夜电影| 欧美成人免费av一区二区三区| h日本视频在线播放| 色综合亚洲欧美另类图片| 99riav亚洲国产免费| 成人美女网站在线观看视频| 精品久久久久久,| 欧美黑人巨大hd| 真实男女啪啪啪动态图| 亚洲精品一区av在线观看| 午夜亚洲福利在线播放| 欧美一级a爱片免费观看看| 成人国产一区最新在线观看| 欧美xxxx黑人xx丫x性爽| av在线老鸭窝| 久久人人精品亚洲av| 欧美日韩乱码在线| 极品教师在线视频| 日本黄大片高清| 久久久久久大精品| 中文字幕人妻熟人妻熟丝袜美| 12—13女人毛片做爰片一| 淫秽高清视频在线观看| 亚洲美女搞黄在线观看 | 国产aⅴ精品一区二区三区波| 日韩欧美 国产精品| 九九热线精品视视频播放| 国产人妻一区二区三区在| 国产亚洲精品综合一区在线观看| 精品国产三级普通话版| 非洲黑人性xxxx精品又粗又长| 老司机深夜福利视频在线观看| 男女做爰动态图高潮gif福利片| 精品欧美国产一区二区三| 日韩欧美在线乱码| 久久99热这里只有精品18| 男人和女人高潮做爰伦理| 韩国av一区二区三区四区| 精品久久久久久久久亚洲 | 国产 一区 欧美 日韩| 我的老师免费观看完整版| or卡值多少钱| 国产精品一区二区免费欧美| 婷婷色综合大香蕉| 毛片女人毛片| 亚洲欧美日韩东京热| 热99在线观看视频| 在线免费观看不下载黄p国产 | 男人舔女人下体高潮全视频| 国产麻豆成人av免费视频| 老司机深夜福利视频在线观看| 国产91精品成人一区二区三区| 三级毛片av免费| 欧美日韩瑟瑟在线播放| 成年女人永久免费观看视频| 亚洲第一欧美日韩一区二区三区| 国产高清视频在线播放一区| 日本a在线网址| 国产综合懂色| 日韩免费av在线播放| 亚洲av.av天堂| 在线观看午夜福利视频| 欧美日本视频| 国产高清视频在线观看网站| 久9热在线精品视频| 亚洲最大成人手机在线| 神马国产精品三级电影在线观看| 最新中文字幕久久久久| 亚洲狠狠婷婷综合久久图片| 别揉我奶头~嗯~啊~动态视频| 亚洲第一欧美日韩一区二区三区| 午夜视频国产福利| 国产精品一区二区性色av| 久久久久久久久大av| 亚洲av熟女| 日本撒尿小便嘘嘘汇集6| 亚洲av成人精品一区久久| 国内精品美女久久久久久| 欧美黑人欧美精品刺激| 欧美成人性av电影在线观看| x7x7x7水蜜桃| 夜夜躁狠狠躁天天躁| 国产精品99久久久久久久久| 亚洲专区中文字幕在线| 国产高清视频在线观看网站| 偷拍熟女少妇极品色| 久久久色成人| 国产精品久久电影中文字幕| 欧美高清成人免费视频www| 国内少妇人妻偷人精品xxx网站| 黄色女人牲交| 男女之事视频高清在线观看| 51国产日韩欧美| 在线播放国产精品三级| 免费看光身美女| 亚洲av成人不卡在线观看播放网| 久久精品影院6| 成人无遮挡网站| 亚洲不卡免费看| 国产伦在线观看视频一区| av在线观看视频网站免费| 亚洲精品影视一区二区三区av| 波多野结衣巨乳人妻| 国产蜜桃级精品一区二区三区| 午夜激情福利司机影院| 亚洲专区国产一区二区| 亚洲av熟女| a级一级毛片免费在线观看| 国产主播在线观看一区二区| 国产精品美女特级片免费视频播放器| 夜夜看夜夜爽夜夜摸| a级毛片a级免费在线| 国产69精品久久久久777片| 国产免费男女视频| 精品久久久久久,| а√天堂www在线а√下载| 亚洲成a人片在线一区二区| 久久久久免费精品人妻一区二区| 有码 亚洲区| 色尼玛亚洲综合影院| 免费看a级黄色片| 亚洲最大成人手机在线| 在线观看66精品国产| 嫩草影院精品99| 观看免费一级毛片| 草草在线视频免费看| 桃色一区二区三区在线观看| 91久久精品电影网| 男插女下体视频免费在线播放| 天天一区二区日本电影三级| 丰满乱子伦码专区| 国内精品久久久久精免费| 婷婷色综合大香蕉| 首页视频小说图片口味搜索| 又粗又爽又猛毛片免费看| 十八禁人妻一区二区| 亚洲最大成人中文| 欧美日韩乱码在线| 在线观看舔阴道视频| 亚洲av成人av| 国产成人a区在线观看| 最新中文字幕久久久久| 老司机午夜福利在线观看视频| 国产爱豆传媒在线观看| 麻豆国产av国片精品| 亚洲av二区三区四区| 日韩欧美免费精品| 国产精品不卡视频一区二区 | 欧美另类亚洲清纯唯美| 无遮挡黄片免费观看| av女优亚洲男人天堂| 国产精品1区2区在线观看.| 久久久精品欧美日韩精品| 人妻丰满熟妇av一区二区三区| 草草在线视频免费看| 欧美黑人巨大hd| 看片在线看免费视频| 国产欧美日韩一区二区三| 免费看日本二区| 中文字幕av成人在线电影| 97热精品久久久久久| 色综合婷婷激情| 成人特级黄色片久久久久久久| 又紧又爽又黄一区二区| 亚洲av美国av| 99热这里只有是精品在线观看 | 波多野结衣高清作品| 国产精品美女特级片免费视频播放器| 午夜日韩欧美国产| 久久国产乱子免费精品| 亚洲黑人精品在线| 观看免费一级毛片| 精华霜和精华液先用哪个| 免费在线观看成人毛片| 好男人电影高清在线观看| 毛片女人毛片| 国产极品精品免费视频能看的| 亚洲精品影视一区二区三区av| 人妻丰满熟妇av一区二区三区| 亚洲精品亚洲一区二区| 偷拍熟女少妇极品色| 在线免费观看不下载黄p国产 | 动漫黄色视频在线观看| 成人无遮挡网站| 婷婷亚洲欧美| 久久亚洲真实| 欧美最新免费一区二区三区 | 人人妻人人看人人澡| 亚洲18禁久久av| 国内久久婷婷六月综合欲色啪| 少妇人妻一区二区三区视频| 97热精品久久久久久| 精品久久久久久久久久免费视频| 免费在线观看影片大全网站| 色哟哟·www| 免费在线观看成人毛片| 日韩欧美国产一区二区入口| 国产真实乱freesex| 动漫黄色视频在线观看| 欧美成狂野欧美在线观看| 国产精品亚洲一级av第二区| www.www免费av| 免费无遮挡裸体视频| 最近中文字幕高清免费大全6 | 亚洲国产日韩欧美精品在线观看| av欧美777| 久久6这里有精品| 亚洲一区高清亚洲精品| 成人高潮视频无遮挡免费网站| 久久国产乱子伦精品免费另类| 久久6这里有精品| av黄色大香蕉| 欧美+亚洲+日韩+国产| 欧美成人性av电影在线观看| 亚洲人成网站高清观看| 自拍偷自拍亚洲精品老妇| 老鸭窝网址在线观看| 国产精品野战在线观看| 欧美日本亚洲视频在线播放| 亚洲精华国产精华精| 成人欧美大片| 国产真实伦视频高清在线观看 | 可以在线观看的亚洲视频| 国内精品一区二区在线观看| 深夜精品福利| 国产成人aa在线观看| 韩国av一区二区三区四区| 亚洲精品成人久久久久久| 欧美xxxx黑人xx丫x性爽| 夜夜躁狠狠躁天天躁| 久久人人精品亚洲av| 亚洲激情在线av| 少妇人妻精品综合一区二区 | 永久网站在线| 深爱激情五月婷婷| 最近最新中文字幕大全电影3| 日韩欧美 国产精品| 色播亚洲综合网| 国产精品电影一区二区三区| 欧美又色又爽又黄视频| 91久久精品国产一区二区成人| 亚洲精品乱码久久久v下载方式| 亚洲国产高清在线一区二区三| 露出奶头的视频| av黄色大香蕉| 午夜免费成人在线视频| 日韩欧美精品免费久久 | 国产一区二区三区在线臀色熟女| 国产精品不卡视频一区二区 | 人人妻人人看人人澡| 国产高清激情床上av| 国产精品嫩草影院av在线观看 | 欧美高清性xxxxhd video| 搡老岳熟女国产| 欧美成人免费av一区二区三区| eeuss影院久久| 99久国产av精品| 美女免费视频网站| 直男gayav资源| 久久久久久久久大av| 国产综合懂色| 久久精品影院6| 国产av麻豆久久久久久久| 久久精品综合一区二区三区| 99热这里只有精品一区| 日韩中文字幕欧美一区二区| 国产一区二区激情短视频| 精品99又大又爽又粗少妇毛片 | 婷婷精品国产亚洲av在线| 午夜老司机福利剧场| 三级男女做爰猛烈吃奶摸视频| 精品久久久久久久末码| 国产极品精品免费视频能看的| 精品欧美国产一区二区三| 一区二区三区高清视频在线| 在线看三级毛片| 无遮挡黄片免费观看| 91麻豆av在线| 亚洲av成人不卡在线观看播放网| 一边摸一边抽搐一进一小说| 小蜜桃在线观看免费完整版高清| 一级a爱片免费观看的视频| 免费大片18禁| 直男gayav资源| 亚洲美女黄片视频| 一本一本综合久久| 婷婷精品国产亚洲av| 极品教师在线免费播放| 久久伊人香网站| 久久久成人免费电影| 国产 一区 欧美 日韩| 精品午夜福利在线看| av专区在线播放| 欧美精品国产亚洲| 又爽又黄a免费视频| 在线观看一区二区三区| 精品久久久久久久久av| www.www免费av| 精品福利观看| 丝袜美腿在线中文| 日本熟妇午夜| 日日干狠狠操夜夜爽| 久久久久久久久久成人| 舔av片在线| 18+在线观看网站| 亚洲中文字幕一区二区三区有码在线看| 深爱激情五月婷婷| 成人欧美大片| 精品乱码久久久久久99久播| 亚洲国产精品久久男人天堂| eeuss影院久久| 国产69精品久久久久777片| 国产精品影院久久| 好男人在线观看高清免费视频| 日韩中文字幕欧美一区二区| 国产精品自产拍在线观看55亚洲| 日日夜夜操网爽| 长腿黑丝高跟| 国产成人a区在线观看| 丰满人妻一区二区三区视频av| 色尼玛亚洲综合影院| www日本黄色视频网| 三级毛片av免费| 日韩中文字幕欧美一区二区| 在线观看午夜福利视频| 在线a可以看的网站| 久久午夜福利片| 99久久九九国产精品国产免费| 亚洲av一区综合| 少妇的逼水好多| 国产一区二区三区在线臀色熟女| 亚洲欧美激情综合另类| 精品人妻视频免费看| 欧洲精品卡2卡3卡4卡5卡区| 国产高清视频在线观看网站| 在线观看舔阴道视频| 亚洲aⅴ乱码一区二区在线播放| 麻豆成人av在线观看| 国产真实乱freesex| 亚洲乱码一区二区免费版| 亚洲av电影不卡..在线观看| 天堂动漫精品| 亚洲成av人片在线播放无| 69av精品久久久久久| 动漫黄色视频在线观看| 在线看三级毛片| 亚洲片人在线观看| 欧美最黄视频在线播放免费| 搡老妇女老女人老熟妇| 99视频精品全部免费 在线| 欧美日韩福利视频一区二区| 少妇人妻一区二区三区视频| 国产 一区 欧美 日韩| 久久久久久久午夜电影| 婷婷亚洲欧美| 亚洲国产色片| av国产免费在线观看| 真实男女啪啪啪动态图| 欧美黄色片欧美黄色片| 夜夜躁狠狠躁天天躁| 国产成年人精品一区二区| 老女人水多毛片| 亚洲av美国av| 精品久久久久久久人妻蜜臀av| 国产高清有码在线观看视频| 久久精品91蜜桃| 欧美色欧美亚洲另类二区| 久久亚洲精品不卡| 免费一级毛片在线播放高清视频| 精品免费久久久久久久清纯| 在线天堂最新版资源| x7x7x7水蜜桃| www.色视频.com| 真人做人爱边吃奶动态| 啪啪无遮挡十八禁网站| 成人av一区二区三区在线看| 色综合欧美亚洲国产小说| 国产成年人精品一区二区| 好看av亚洲va欧美ⅴa在| 国模一区二区三区四区视频| 97碰自拍视频| 国内久久婷婷六月综合欲色啪| 国产一区二区激情短视频| 久久久久国产精品人妻aⅴ院| 久久午夜亚洲精品久久| 亚洲欧美日韩无卡精品| 精品久久久久久久久久免费视频| 少妇裸体淫交视频免费看高清| 国产亚洲精品久久久久久毛片| x7x7x7水蜜桃| 丰满的人妻完整版| 91在线观看av| www.www免费av| 精品人妻1区二区| 丰满人妻熟妇乱又伦精品不卡| 床上黄色一级片| 成人高潮视频无遮挡免费网站| 成年人黄色毛片网站| 深夜a级毛片| 精品久久久久久久久久免费视频| 午夜亚洲福利在线播放| 色哟哟哟哟哟哟| 99久国产av精品| 色综合亚洲欧美另类图片| 国产aⅴ精品一区二区三区波| 亚洲av免费在线观看| 在线观看美女被高潮喷水网站 | 成人精品一区二区免费| 精品久久国产蜜桃| 看片在线看免费视频| 亚洲性夜色夜夜综合| 女同久久另类99精品国产91| 国产亚洲欧美98| 免费电影在线观看免费观看| 成人永久免费在线观看视频| 欧美成狂野欧美在线观看| 91av网一区二区| 怎么达到女性高潮| 91在线精品国自产拍蜜月| av在线蜜桃| 色av中文字幕| 国产三级黄色录像| 亚洲专区国产一区二区| 欧美色欧美亚洲另类二区| 久久中文看片网| 亚洲无线观看免费| 成年人黄色毛片网站| 18禁黄网站禁片免费观看直播| 啦啦啦观看免费观看视频高清| 嫩草影院新地址| 12—13女人毛片做爰片一| 蜜桃久久精品国产亚洲av| 亚洲av日韩精品久久久久久密| av天堂中文字幕网| 成年免费大片在线观看| 亚洲熟妇中文字幕五十中出| 国产毛片a区久久久久| 能在线免费观看的黄片| 久久99热这里只有精品18| 亚洲国产日韩欧美精品在线观看| 禁无遮挡网站| eeuss影院久久| 看片在线看免费视频| 成人无遮挡网站| 色在线成人网| 天堂网av新在线| 69av精品久久久久久| 日韩精品青青久久久久久| 又爽又黄a免费视频| 一级毛片久久久久久久久女| 国产免费一级a男人的天堂| 内地一区二区视频在线| 舔av片在线| 国产一级毛片七仙女欲春2| 成人亚洲精品av一区二区| 琪琪午夜伦伦电影理论片6080| 国产精品98久久久久久宅男小说| 国产视频内射| 亚洲自偷自拍三级| 亚洲国产精品合色在线| 少妇人妻精品综合一区二区 | 国产麻豆成人av免费视频| 日韩成人在线观看一区二区三区| 久久久久久久久久黄片| 久久草成人影院| 日韩欧美一区二区三区在线观看| 九色国产91popny在线| 色哟哟·www| 国产伦在线观看视频一区| 午夜福利在线观看吧| 91久久精品电影网| 免费看日本二区| 精品乱码久久久久久99久播| 午夜日韩欧美国产| 色播亚洲综合网| 两个人的视频大全免费| 亚洲欧美日韩高清专用| 91av网一区二区| 免费在线观看成人毛片| www.色视频.com| 精品国内亚洲2022精品成人| 亚洲av一区综合| 在线a可以看的网站| 亚洲人成电影免费在线| 亚洲真实伦在线观看| 欧美高清性xxxxhd video| 亚洲综合色惰| 99热精品在线国产| 国内毛片毛片毛片毛片毛片| 欧美日韩瑟瑟在线播放| 俺也久久电影网|