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

    厄米-高斯光束在熱非局域介質(zhì)中傳輸?shù)臄?shù)值模擬研究*

    2011-10-23 12:13:52李少華楊振軍陸大全
    物理學(xué)報(bào) 2011年2期
    關(guān)鍵詞:孤子局域光束

    李少華 楊振軍 陸大全 胡 巍

    (華南師范大學(xué)信息光子技術(shù)廣東省高校重點(diǎn)實(shí)驗(yàn)室,廣州 510631)(2010年1月27日收到;2010年5月6日收到修改稿)

    厄米-高斯光束在熱非局域介質(zhì)中傳輸?shù)臄?shù)值模擬研究*

    李少華 楊振軍 陸大全 胡 巍?

    (華南師范大學(xué)信息光子技術(shù)廣東省高校重點(diǎn)實(shí)驗(yàn)室,廣州 510631)(2010年1月27日收到;2010年5月6日收到修改稿)

    基于非線性薛定諤方程和熱擴(kuò)散的泊松方程,采用分步傅里葉算法以及多重網(wǎng)格法對厄米-高斯光束在不同形狀的熱非局域介質(zhì)鉛玻璃中的傳輸進(jìn)行了數(shù)值模擬研究.結(jié)果表明,低階厄米-高斯光束可以較為穩(wěn)定地在鉛玻璃中傳輸.高階厄米-高斯光束在鉛玻璃中傳輸變得不穩(wěn)定,并且階數(shù)越高,穩(wěn)定性越差.樣品的形狀對于厄米-高斯光束的影響很大.在正方形樣品中,厄米-高斯光束的傳輸與Snyder-Mitchell模型符合得相對較好.在矩形樣品中厄米-高斯光束在傳輸過程中的強(qiáng)度分布將發(fā)生較大的變化.

    非局域非線性介質(zhì),厄米-高斯光束,光束傳輸

    PACS:42.65.Tg,42.65.Jx

    1.引 言

    光孤子自從 1973年被 Hasegawa和 Tappert[1]提出以后,引起了人們廣泛研究.空間光孤子是光束在傳播過程中由非線性效應(yīng)平衡衍射效應(yīng)的結(jié)果,空間光孤子一直是非線性光學(xué)研究前沿.光孤子在全光網(wǎng)絡(luò),光通信以及光邏輯器件方面有著非常重要的應(yīng)用[2—4]. 在 Snyder 和 Mitchell[5]開創(chuàng)性地提出強(qiáng)非局域下空間光孤子模型后,有關(guān)強(qiáng)非局域非線性介質(zhì)中的孤子研究在近幾年一直是熱點(diǎn).Guo等[6]在理論上提出和論證了在強(qiáng)非局域下孤子傳輸會(huì)出現(xiàn)大相移現(xiàn)象,為邏輯門和光開關(guān)的實(shí)現(xiàn)提供了一個(gè)強(qiáng)有力的理論指導(dǎo).但是Snyder模型為簡化的理想模型,在真實(shí)的物理系統(tǒng)中,光束的傳輸將變得更加復(fù)雜,而且理論和數(shù)值模擬表明,高階孤子不能夠穩(wěn)定傳輸.在不同非局域程度下,不同非線性介質(zhì)中的光孤子的傳輸性質(zhì)以及孤子間的相互作用的研究已取得了很大的成果[7—12].

    最近有人展開一系列關(guān)于有界非局域介質(zhì)的理論和實(shí)驗(yàn)研究[13—15],其中包括液晶和鉛玻璃.不像在弱非局域和高斯響應(yīng)的介質(zhì),在熱非局域介質(zhì)中,熱效應(yīng)的響應(yīng)函數(shù)和介質(zhì)邊界的熱擴(kuò)散有關(guān).因此,在熱非局域介質(zhì)中,樣品的幾何形狀對于孤子的結(jié)構(gòu)和形狀具有巨大的影響.鉛玻璃是一種目前被證實(shí)比較典型的熱非局域介質(zhì),研究在鉛玻璃中的孤子傳輸性質(zhì)也引起了廣泛的關(guān)注.Segev等[16]在鉛玻璃試驗(yàn)中發(fā)現(xiàn)了橢圓孤子以及單孤子偏轉(zhuǎn)振蕩的現(xiàn)象,并且還發(fā)現(xiàn)了非局域表面孤子.我們通過理論研究推算出了有界的鉛玻璃的響應(yīng)函數(shù)以及光束中心的振蕩軌跡,并且給出了光束振蕩周期和輸入功率的理論關(guān)系[17].

    本文在以上工作的基礎(chǔ)上做了針對厄米-高斯在鉛玻璃里傳輸性質(zhì)的數(shù)值研究.在前期的工作中,只是研究了光束中心軌跡的傳輸,未涉及光束光強(qiáng)分布的演化和穩(wěn)定性問題.本文進(jìn)一步討論了在熱非局域介質(zhì)中,不同階數(shù)厄米-高斯光束傳輸過程中所表現(xiàn)出來的光強(qiáng)分布的演化和高階厄米-高斯光束傳輸?shù)姆€(wěn)定性,以及在具有正方形和矩形邊界的樣品中,厄米-高斯光束傳輸過程中展現(xiàn)出的不同性質(zhì).同時(shí)對于不同入射方式下厄米-高斯光束的傳輸性質(zhì)也進(jìn)行了比較.

    2.理論模型

    在熱非局域介質(zhì)中,光束的光強(qiáng)是介質(zhì)中的熱源,熱量向邊界擴(kuò)散,引起介質(zhì)溫度的重新分布,從而導(dǎo)致折射率變化,所以在熱非局域介質(zhì)中的折射率分布不僅由光束的光強(qiáng)決定,而且介質(zhì)的邊界對折射率分布也有著重要的影響.

    傍軸光束在鉛玻璃中的傳輸可以由非線性薛定諤方程[16]來描述.非線性薛定諤方程由下式給出

    其中E代表傍軸光束,k表示僅考慮線性折射率時(shí)的波數(shù),n0表示線性折射率,Δn表示折射率的改變量.鉛玻璃中的熱擴(kuò)散過程可以用泊松方程[18]來描述

    其中 κ 是熱傳導(dǎo)系數(shù),T(x,y,z)表示溫度分布,I(x,y,z)表示光束的光強(qiáng),α是介質(zhì)的熱吸收系數(shù).折射率的變化可以表示為 Δn=β(T-T0),其中 β是介質(zhì)的熱光系數(shù),T0表示邊界處的溫度(本文假設(shè)邊界處為恒溫).

    為數(shù)值模擬的方便,采用下述歸一化對(1)和(2)式進(jìn)行簡化.得到歸一化后的(1)和(2)式為:

    3.數(shù)值分析

    本文主要模擬和討論厄米-高斯光束在方形和矩形鉛玻璃里中的傳輸,分析邊界對厄米-高斯光束的影響,重點(diǎn)討論厄米-高斯光束在傳輸中的強(qiáng)度分布的變化以及其在傳輸中的穩(wěn)定性.厄米-高斯光束可以表示為其中A是與功率有關(guān)的常數(shù).在數(shù)值模擬中,正方形樣品的邊長用D表示;矩形樣品的長邊和短邊分別用DL和DS分別表示,矩形樣品的長寬比為τ=DL/Ds=2,其中短邊為X方向,長邊為Y方向.光束都是采用中心入射的方式.

    3.1.基模高斯光束

    首先研究(0,0)階厄米-高斯光束,即基模高斯光束在鉛玻璃中的傳輸性質(zhì).在液晶介質(zhì)中,不存在嚴(yán)格的孤子功率[21],即保持光束寬度始終不變的功率.同樣鉛玻璃中由于邊界的作用使熱擴(kuò)散不均勻,從而邊界對于光束的作用也不均勻,所以也不存在光束寬度保持嚴(yán)格不變的孤子.事實(shí)上我們經(jīng)過大量的模擬證明了這一點(diǎn),無論如何調(diào)整輸入功率,光束的束寬總是保持著振蕩的,所以我們定義能夠保持光束寬度的變化小于2%的功率為臨界功率PC.在方形樣品和矩形樣品中,高斯光束都可以存在形成孤子的臨界功率.臨界功率下,從束寬W的變化圖(如圖1)發(fā)現(xiàn),由于樣品的形狀不同,孤子的束寬變化不同(但束寬變化都小于2%).在方形樣品中(圖1(a)),X方向和Y方向的束寬變化是一致的,所以光束總是呈圓形;然而在矩形樣品中(圖1(b)),束寬在X方向和Y方向上有微小的差別:短邊X方向的束寬基本上總是小于長邊Y方向的束寬,但光束的橢圓度很小,所以很難分辨出光束是橢圓型的.我們認(rèn)為其原因是,在矩形樣品中,熱量沿樣品的短邊即X方向擴(kuò)散時(shí)先受到邊界的作用,所以對光束在X方向上有一個(gè)比Y方向更大的壓縮作用.

    其次,我們模擬和研究了非臨界功率時(shí)高斯光束在鉛玻璃中的傳輸性質(zhì).經(jīng)過大量的數(shù)值模擬可知,非臨界功率下,高斯光束在方形樣品中在非臨界功率下傳輸時(shí)可以形成穩(wěn)定的呼吸子,而在矩形樣品中傳輸時(shí)是一個(gè)準(zhǔn)周期橢圓呼吸的過程.圖2給出了3倍臨界功率下高斯光束的束寬變化圖.其中圖2(a)是方形樣品中高斯光束在大功率傳輸時(shí)X方向和Y方向的束寬變化曲線,虛線表示X方向的束寬變化,實(shí)線表示Y方向的束寬變化.可以看到虛線和實(shí)線基本重合,很難分辨.這說明兩個(gè)方向上的束寬變化是完全一致的,光束在傳輸過程中的呼吸效應(yīng)沒有改變光束的形狀,光束一直是保持一個(gè)圓形.圖2(b)是高斯光束在大功率下矩形樣品中傳輸時(shí)X方向和Y方向的束寬變化曲線,從結(jié)果可以看到X方向和Y方向的束寬變化及其包絡(luò)變化都是一個(gè)準(zhǔn)周期性呼吸變化.當(dāng)Y方向的束寬振蕩到包絡(luò)的最頂點(diǎn)時(shí),X方向的束寬正好振蕩到束寬包絡(luò)的最低谷,此時(shí)光束的形狀是一個(gè)橢圓.如Z=7處的光束強(qiáng)圖分布圖(圖3(b)).當(dāng)X方向和Y方向的束寬變化曲線重合時(shí),此時(shí)光束的形狀又回到了圓(圖3(c)).這是由于光束在傳輸過程中,光束強(qiáng)度的變化是一個(gè)分?jǐn)?shù)傅里葉變換的過程[22].當(dāng)光束成一個(gè)橢圓形狀時(shí),分?jǐn)?shù)傅里葉變換在X軸和Y軸方向的積分范圍不同,使得光束光強(qiáng)變化在X軸和Y軸方向上分?jǐn)?shù)傅里葉變換的周期不相等,導(dǎo)致在光束的傳輸過程中,光束的兩個(gè)方向上的束寬變化步調(diào)不一致,也就是說橢圓的長短軸一直在改變,即模擬結(jié)果中光束的束寬一直在變化.同時(shí)從圖2(b)可以看出,當(dāng)光束束寬越大時(shí),包絡(luò)的振蕩越劇烈,光束束寬越小時(shí),包絡(luò)的振蕩越平緩.由此可知光束束寬越大,邊界對于光束的影響越大.比較圖2(a)和(b)的實(shí)線可知,在X方向上光束矩形樣品中的束寬總是小于方形樣品中的束寬,所以光束在矩形短邊X方向上受到了壓縮作用,這和臨界功率下高斯光束的性質(zhì)是一致的,只是在大功率下,光束橢圓度的較高,邊界作用較明顯.

    圖2 3倍臨界功率下高斯光束的光束束寬圖 (a)方形樣品中光束束寬圖,其中模擬參數(shù)是A3pc=12.0,D=20;(b)矩形樣品中光束束寬圖,其中模擬參數(shù)是A3pc=12.0,DL=20

    為了更直觀地理解樣品形狀對于高斯光束的影響,我們模擬了光束在不同傳輸距離時(shí)的強(qiáng)度分布圖,見圖3和4.從光束傳輸強(qiáng)度分布圖可以明顯看出:在大功率下,高斯光束在不同形狀的樣品中呈現(xiàn)不同的傳輸性質(zhì),在方形樣品中可以形成較穩(wěn)定的呼吸子,并且可以穩(wěn)定傳輸很遠(yuǎn)的距離.而在矩形樣品中卻出現(xiàn)準(zhǔn)周期性橢圓呼吸的形狀,這與上面的分析和結(jié)果符合得很好.為驗(yàn)證高斯光束在鉛玻璃中傳輸?shù)姆€(wěn)定性,我們模擬了加入微擾的情況下、高斯光束臨界功率時(shí)在方形樣品和矩形樣品中的傳輸.假設(shè)初始光場為(1+η·ρ(X,Y))ψ,其中 ρ(X,Y)為一個(gè)在(-1,1)之間隨機(jī)取值的函數(shù),η=0.1,模擬結(jié)果見圖5.從圖5可知,在加入微擾的情況下,無論在方形樣品還是矩形樣品中,高斯孤子都是穩(wěn)定的.但是經(jīng)過大量的模擬發(fā)現(xiàn),在矩形樣品中這種穩(wěn)定性會(huì)隨著樣品的長度變小而被破壞,隨著樣品邊界的減小,光束傳輸越來越不穩(wěn)定.

    3.2.低階厄米-高斯光束

    在體介質(zhì)中,由高斯光束取π的相位差形成的偶極孤子是可以線性穩(wěn)定傳輸?shù)模?3].對于非局域非線性介質(zhì),邊界對于偶極孤子的傳輸有著很大的影響.例如,在由液晶分子排列方向決定非局域程度的向列液晶介質(zhì)中,偶極孤子的振幅和相位的振蕩性已被充分地討論和研究[24],向列液晶分子排列方向?qū)ε紭O孤子的傳輸穩(wěn)定性有著很大的影響.在這節(jié)中我們討論厄米-高斯光束中的偶極孤子,即(0,1)厄米-高斯光束在鉛玻璃中的傳輸性質(zhì),研究鉛玻璃介質(zhì)邊界對(0,1)厄米-高斯光束傳輸性質(zhì)的影響以及光束傳輸?shù)姆€(wěn)定性.對于(0,1)階厄米-高斯光束入射矩形樣品時(shí),可采取兩種不同的入射方式:一種是平行入射方式,即光束的長邊平行于矩形樣品的長邊(如圖6(a));另一種是正交入射方式,即光束的長邊平行于矩形樣品的短邊(如圖6(b)).這兩種不同的入射方式使光束在樣品中產(chǎn)生不同的折射率分布,從而能夠使得邊界對光束在鉛玻璃中的傳輸產(chǎn)生不同的影響.

    和基模高斯光束一樣,先從束寬的角度來分析低階厄米-高斯光束的傳輸性質(zhì).經(jīng)過大量的模擬發(fā)現(xiàn),無論如何調(diào)整入射功率,都找不到光束在 X和Y方向同時(shí)保持不變的功率.當(dāng)光束束寬在Y方向上保持微小振蕩時(shí),在X方向上是明顯呼吸效應(yīng)(如圖7);而當(dāng)束寬保持X方向微小振蕩時(shí),在Y方向上則是明顯的呼吸效應(yīng),這和基模高斯光束在臨界功率下兩個(gè)方向都能保持孤子傳輸有著明顯不同.這是由于(0,1)階厄米-高斯光束是非對稱光束,X和Y方向的初始束寬不相等,而在鉛玻璃中,孤子的臨界功率只是和光束的束寬有關(guān)[17],所以非對稱光束在X和Y方向的臨界功率是不相等的.比較圖7(a)中的3條束寬變化曲線可知,在臨界功率下,平行入射矩形樣品時(shí)Y方向上的束寬基本上總是最大,方形樣品次之,正交入射時(shí)最小.這也就是說,矩形的短邊會(huì)給光束一個(gè)比長邊更大的壓縮力,使該方向上的光束被壓縮.

    為了深入理解邊界對低階厄米-高斯光束的影響,我們接著模擬了低階厄米-高斯光束以不同方式入射矩形樣品在不同傳輸距離時(shí)的強(qiáng)度分布(見圖8和9).對比圖8(b)和9(b)可以看出,矩形的短邊對于光束有壓縮的作用,這一點(diǎn)正好說明了上面束寬分析的正確性.并且從折射率等高線圖(圖10)可以看出,在矩形樣品中,短邊方向的溫度梯度大于長邊方向的溫度梯度,也就是短邊方向的折射率梯度大于長邊方向的折射率梯度,在等高線圖中就是短邊方向的等高線密度密一些,而長邊方向的等高線密度疏一些.這樣兩者相比較,光束在短邊方向上就如同被壓縮,相對在長邊方向上就如同被拉伸.

    圖9 矩形樣品中(1,0)階厄米-高斯光束(正交入射)的強(qiáng)度分布圖 (a)Z=0,(b)Z=1,(c)Z=2,(d)Z=3,(e)Z=4,(f)Z=5,(g)Z=64,(h)Z=100(模擬參數(shù)和光束圖相同)

    為了驗(yàn)證低階厄米-高斯光束在熱非局域介質(zhì)中傳輸?shù)姆€(wěn)定性,我們模擬了加入微擾后的低階厄米-高斯光束(1+ η·ρ(X,Y))ψ,在臨界功率下在方形樣品和矩形樣品中傳輸?shù)姆€(wěn)定性(圖11).比較圖11各分圖可知,在加入了微擾的情況下,光束平行入射矩形樣品的傳輸最穩(wěn)定,在方形樣品的傳輸穩(wěn)定性次之,正交入射矩形樣品的傳輸最不穩(wěn)定,這和偶極孤子[25]的穩(wěn)定性相匹配.

    圖11 在加入10%微擾的下(0,1)階厄米-高斯光束的強(qiáng)度圖(a)平行入射矩形樣品,Z=0;(b)平行入射矩形樣品,Z=390;(c)正交入射矩形邊界,Z=0;(d)正交入射矩形邊界,Z=64;(e)方形樣品,Z=0;(f)方形樣品,Z=64

    3.3.高階厄米-高斯光束

    在這一部分將討論高階厄米-高斯光束在鉛玻璃中的傳輸性質(zhì).在Snyder-Mitchell模型下高階厄米-高斯光束可以形成穩(wěn)定的孤子的傳輸.而在高斯響應(yīng)下,非局域程度對于高階厄米-高斯光束傳輸有著很大的影響[26].但Snyder-Mitchell模型與高斯響應(yīng)都是理想情況,實(shí)際中并不存在,所以研究在有界熱非局域介質(zhì)中的高階厄米-高斯光束的傳輸性質(zhì)更有實(shí)際意義.我們通過數(shù)值模擬發(fā)現(xiàn),高階厄米-高斯光束在鉛玻璃中的傳輸都是不穩(wěn)定的,但是光束在近距離傳輸過程中光強(qiáng)分布會(huì)出現(xiàn)周期性的演化.

    高階厄米-高斯光束有多個(gè)極,并且各極之間呈矩陣分布,所以我們把高階厄米-高斯光束稱為(m+1)×(n+1)矩陣光束,例如(1,1)階厄米-高斯光束稱為2×2矩陣光束,更高階的光束依次類推.首先討論(1,1)階厄米-高斯光束,圖12是2×2矩陣光束在鉛玻璃中傳輸?shù)讲煌嚯x的強(qiáng)度分布圖.從模擬結(jié)果可以看出,無論是在方形樣品還是矩形樣品中,2×2矩陣光束都可以保持4極狀態(tài)穩(wěn)定地傳輸一段距離.在穩(wěn)定傳輸一段距離后,我們觀察到了兩種光束花樣的相互轉(zhuǎn)化,光束從四極矩陣光束演化成一個(gè)類似圓環(huán)狀的光束,并且這種光束花樣演化是具有周期性的,盡管這種周期性是不穩(wěn)定的.這和高斯響應(yīng)下厄米-高斯和拉蓋爾高斯光束簇可以相互轉(zhuǎn)化的性質(zhì)是相一致的[27],這說明鉛玻璃作為一種典型的實(shí)際存在的強(qiáng)非局域介質(zhì)能體現(xiàn)著強(qiáng)非局域介質(zhì)的固有的性質(zhì).同時(shí)比較圖12(a)—(f)和(g)—(l)可以看出:2×2矩陣光束在方形樣品中可以穩(wěn)定傳輸?shù)木嚯x比在矩形樣品中穩(wěn)定傳輸?shù)木嚯x要長,這是由于方形樣品相比矩形樣品有著更加對稱的折射率分布,從而能夠?qū)哂芯鶆驅(qū)ΨQ性質(zhì)的2×2矩陣光束各極保持更加均衡的作用,因此在方形樣品中可以穩(wěn)定傳輸更遠(yuǎn)的距離.而且由于樣品形狀不一樣,矩陣光束的各極受到邊界的壓縮作用不一樣,相當(dāng)于對光束施加了一個(gè)力矩作用,因此光束發(fā)生了一個(gè)角度的偏轉(zhuǎn).

    為了更加全面地研究熱非局域介質(zhì)中高階厄米-高斯光束強(qiáng)度花樣演化的性質(zhì),我們繼續(xù)討論更高階厄米-高斯光束的傳輸性質(zhì).同樣,首先研究光束在方形樣品中的傳輸性質(zhì),與四極光束一樣,(1,2)階厄米-高斯光束在方形鉛玻璃中傳輸時(shí)也會(huì)發(fā)生類似的周期性的花樣演化.其中2×3矩陣光束轉(zhuǎn)化為一個(gè)半開的類圓環(huán)狀的光束,如圖13(a)—(d),這和體介質(zhì)中孤子態(tài)的相互轉(zhuǎn)化[28]是類似的;但是這種演化在經(jīng)過了兩個(gè)周期后,光束的花樣變化不會(huì)重現(xiàn),這是由于邊界的作用,這種周期性的花樣轉(zhuǎn)化是不穩(wěn)定的.接著模擬矩形樣品中高階厄米-高斯光束的傳輸性質(zhì),以研究樣品形狀對于高階厄米-高斯光束花樣演化的影響.從模擬結(jié)果發(fā)現(xiàn),2×3矩陣光束在傳輸一段距離后都會(huì)演化成一個(gè)橢圓形項(xiàng)鏈光束,如圖13(e)—(h).為了進(jìn)一步驗(yàn)證我們的結(jié)果,我們模擬了3×3和3×4矩陣光束在矩形樣品中的傳輸,同樣會(huì)出現(xiàn)矩陣光束向橢圓形項(xiàng)鏈光束的轉(zhuǎn)化現(xiàn)象.對比在方形樣品中矩陣光束向類圓環(huán)的光束演化,這說明樣品的形狀對于高階厄米-高斯光束的傳輸有很大的影響,樣品形狀不同,花樣演化的形狀也不同.并且可以看出由于矩形樣品的邊界對于對光束的影響不一致,橢圓形項(xiàng)鏈光束的長軸總是平行于矩形邊界的長邊,這和低階厄米-高斯光束在矩形樣品中被壓縮的性質(zhì)是一致的.

    圖12 (1,1)階厄米-高斯光束傳輸不同位置的強(qiáng)度圖 (a)—(f)是方形樣品,Z=0,44,84,92,100,108,花樣演化的周期 ZT=16,模擬參數(shù)是 Apc=3.00,D=20;(g)—(l)是矩形樣品,Z=0,38,76,84,92,100,花樣演化的周期 ZT=16,樣品的形狀沒有改變花樣演化的周期,模擬參數(shù)是 Apc=3.00,DL=20

    圖 13 (1,2)階厄米-高斯光束傳輸不同距離的強(qiáng)度圖 (a)—(d)是在方形樣品,Z=0,70,75,80,Apc=1.00,D=20;(e)—(h)是正交入射矩形樣品,Z=0,10,17,22,Apc=1.00,DL=20

    綜上可知,高階厄米-高斯光束在熱非局域介質(zhì)中不能穩(wěn)定傳輸,階數(shù)越高,光束傳輸越不穩(wěn)定,但是光束會(huì)出現(xiàn)周期性的花樣的演化.樣品的形狀對于光束花樣演化有很大影響.在方形樣品中,光束會(huì)從一個(gè)矩陣光束演化成一個(gè)環(huán)狀的光束,再回到矩陣光束;而在矩形樣品中,高階厄米-高斯光束會(huì)演化成一個(gè)橢圓形項(xiàng)鏈光束.

    4.結(jié) 論

    本文用數(shù)值模擬的方法研究了在熱非局域介質(zhì)中不同階數(shù)厄米-高斯光束的傳輸性質(zhì).通過模擬驗(yàn)證了熱非局域介質(zhì)中低階厄米-高斯孤子的存在,發(fā)現(xiàn)在熱非局域介質(zhì)中,隨著階數(shù)的增加,厄米-高斯光束越來越不穩(wěn)定.對于基模高斯光束而言,在臨界功率下,高斯光束能夠形成穩(wěn)定的孤子傳輸,邊界的影響幾乎可以忽略;在大功率下,邊界對高斯光束的傳輸性質(zhì)有很大的影響.通過比較非對稱厄米-高斯光束不同方式入射矩形樣品的模擬結(jié)果,發(fā)現(xiàn)矩形的短邊對光束具有壓縮作用.對高階厄米-高斯光束而言,光束并不會(huì)像在理想的Snyder模型中可以保持孤子傳輸,而是會(huì)出現(xiàn)有趣的周期性花樣演化,樣品的形狀對于光束花樣演化有著很大的影響,而且高階厄米-高斯光束在鉛玻璃中的傳輸都是不穩(wěn)定的.本文的結(jié)果揭示了鉛玻璃在光信息無失真?zhèn)鬏敽蛨D像處理當(dāng)中的應(yīng)用,同時(shí)對于熱非局域介質(zhì)實(shí)驗(yàn)具有一定的參考意義.在此基礎(chǔ)上我們將繼續(xù)研究在實(shí)際的熱非局域介質(zhì)中其他孤子解的傳輸特性.

    [1]Hasegawa A,Tappert F 1973Appl.Phys.Lett.23 142

    [2]Suter D,Blasberg T 1983Phys.Rev.A 48 4583

    [3]Desyatnikov A S,Sukhorukov A A,Kivshar Y S 2005Phys.Rev.Lett.95 203904

    [4]Kip D,Soljajic M,Segev M,Eugenieva E,Christodoulides D N 2000Science290 495

    [5]Snyder A W,Mitchell D J 1997Science276 1538

    [6]Guo Q,Luo B,Yi F,Chi S,Xie Y 2004Phys.Rev.E 69 016602

    [7]Zhu Y Q,Long X W,Hu W,Cao L G,Yang P B 2008Acta Phys.Sin.57 2260(in Chinese)[朱葉青、龍學(xué)文、胡 巍、曹龍貴、楊平保2008物理學(xué)報(bào)57 2260]

    [8]Yang P B,Cao L G,Hu W,Zhu Y Q,Guo Q,Yang X B 2008Acta Phys.Sin.57 285(in Chinese)[楊平保、曹龍貴、胡巍、朱葉青、郭 旗、楊湘波2008物理學(xué)報(bào)57 285]

    [9]Shao Y Q,Guo Q 2006Acta Phys.Sin.55 2751(in Chinese)[邵毅全、郭 旗2006物理學(xué)報(bào)55 2751]

    [10]Guo Q,Zhang X P,Hu W,Shou Q 2006Acta Phys.Sin.55 1832(in Chinese)[郭 旗、張霞萍、胡 巍、壽 倩 2006物理學(xué)報(bào) 55 1832]

    [11]Zheng Y J,Xuan W T,Lu D Q,OuYang S G,Hu W,Guo Q 2010Acta Phys.Sin.59 1075(in Chinese)[鄭亞建、宣文濤、陸大全、歐陽世根、胡 巍、郭 旗2010物理學(xué)報(bào)59 1075]

    [12]Chen L X,Lu D Q,Hu W,Yang Z J,Cao W W,Zheng R,Guo Q 2010Acta Phys.Sin.59 327(in Chinese)[陳利霞、陸大全、胡 巍、楊振軍、曹偉文、鄭 睿、郭 旗2010物理學(xué)報(bào) 59 327]

    [13]Rotschild C,Segev M,Xu Z,Kartashov Y V,Torner L,Cohen O 2006Opt.Lett.31 3312

    [14]Rotschild C,Alfassi B,Cohen O,Segev M 2006Nature Physics2 769

    [15]Dreischuh A,Neshev D N,Petersen D E,Bang O,Krolikowski W 2006Phys.Rev.Lett.96 043901

    [16]Rotschild C,Cohen O,Manela O,Segev M,Carmon T 2005Phys.Rev.Lett.95 213904

    [17]Liang Y B,Zheng Y J,Yang P B,Cao L G,Lu D Q,Hu W,Guo Q 2008Acta Phys.Sin.57 5690(in Chinese)[梁炎斌、鄭亞建、楊平保、曹龍貴、陸大全、胡 巍、郭 旗2008物理學(xué)報(bào) 57 5690]

    [18]Iturbe-Castillo M D,Stepanov S,Sanchz-Mondragon J J 1996Opt.Lett.21 1622

    [19]Agrawal G P Translated by Jia D F,Yu Z H,Tan B,Hu Z Y 2002Nonlinear Fiber Optics and Application of Nonlinear Fiber Optic(Beijing:Electronic Industry Press)pp34,35(in Chinese)[Agrawal G P著賈東方、余震虹、談 斌、胡智勇譯2002非線性光纖光學(xué)原理及應(yīng)用 (北京:電子工業(yè)出版社)第34,35頁]

    [20]Teukolsky S A,Vetterling W T,F(xiàn)lannery B P Translated by Fu Z Y,Zhao M N,Ding Y 1995Numerical Recips in C:The Art of Scientific Computing(Beijing:Electronic Industry Press)pp734—749(in Chinese)[Teukosky S A,Vetterling W T,F(xiàn)lannery B P著傅祖蕓、趙海娜、丁 巖譯1995 C語言數(shù)值算法程序大全 (北京:電子工業(yè)出版社)第 734—749頁]

    [21]Strinic A I,Petrovic M,Timotijevic D V,Aleksic N B,Belic M R 2009Opt.Express17 11698

    [22]Lu D Q,Hu W,Zheng Y J,Liang Y B,Cao L G,Lan S,Guo Q 2008Phys.Rev.A 78 043815

    [23]Garcia-Ripoll J J,Perez-Garcia V M,Ostrovskaya E A,Kivshar Y S 2000Phys.Rev.Lett.85 82

    [24]Garcia-Reimbert C,Minzoni A A,Marchant T R,Smyth N F,Worthy A L 2008PhysicaD 237 1088

    [25]Ye F W,Kartashov Y V,Torner L 2008Phys.Rev.A 77 043821

    [26]Deng D M,Zhao X,Guo Q,Lan S 2007J.Opt.Soc.Am.B 24 2537

    [27]Buccoliero D,Desyatnikov A S,Krolikowski W,Kivshar Y S 2007Phys.Rev.Lett.98 053901

    [28]Buccoliero D,Desyatnikov A S 2009Opt.Express17 9609

    PACS:42.65.Tg,42.65.Jx

    Numerical study of Hermite-Gaussian beams in nonlocal thermal media*

    Li Shao-Hua Yang Zhen-Jun Lu Da-Quan Hu Wei?
    (Laboratory of Photonic Information Technology,South China Normal University,Guangzhou 510631,China)
    (Received 27 January 2010;revised manuscript received 6 May 2010)

    Based on the nonlocal nonlinear Schr?dinger equation and Poisson equation of thermal diffusion,using the slip-step Fourier algorithm and multi-grid method,we numerically investigated the propagation properties of Hermite-Guassian beams in the nonlocal thermal media.The results show that low-order Hermite-Gaussian beams can propagate stably,in contrast with the unstable propagation of high-order Hermite-Gaussian beams.The worse the stability is,the higher the order is.The effect of the boundary of the sample with different cross sections on the propagation properties of Hermite-Guassian beam is also discussed in detail.We found that propagation properties in square geometry are in agreement with those in Snyder-Mitchell model.However,in rectangular sample,the evolution of intensity distribution of Hermite-Gaussian beams differs seriously from that in the square sample.

    nonlocal nonlinear media,Hermite-Gaussian beams,beam propagation

    *國家自然科學(xué)基金(批準(zhǔn)號(hào):10674050,10804033)、教育部高等學(xué)校博士點(diǎn)專項(xiàng)科研基金(批準(zhǔn)號(hào):200805740002)和廣東省高等學(xué)??萍紕?chuàng)新團(tuán)隊(duì)計(jì)劃(批準(zhǔn)號(hào):06CXTD005)資助的課題.

    ?通訊聯(lián)系人.E-mail:huwei@scnu.edu.cn

    *Project supported by the National Natural Science Foundation of China(Grant Nos.10674050,10804033),the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.200805740002)and the Program for Innovative Research Team of the High Education in Guangdong Province,China(Grant No.06CXTD005).

    ?Corresponding author.E-mail:huwei@scnu.edu.cn

    猜你喜歡
    孤子局域光束
    2維Airy光束陣列強(qiáng)度的調(diào)控技術(shù)研究
    詭異的UFO光束
    奧秘(2021年3期)2021-04-12 15:10:26
    一個(gè)新的可積廣義超孤子族及其自相容源、守恒律
    (3+1)維Potential-Yu-Toda-Sasa-Fukuyama方程新的多周期孤子解
    局域積分散列最近鄰查找算法
    電子測試(2018年18期)2018-11-14 02:30:34
    兩個(gè)孤子方程的高階Painlevé截?cái)嗾归_
    激光共焦顯微光束的偏轉(zhuǎn)掃描
    PET成像的高分辨率快速局域重建算法的建立
    基于局域波法和LSSVM的短期負(fù)荷預(yù)測
    電測與儀表(2015年7期)2015-04-09 11:39:50
    激光探索
    久久人人爽人人爽人人片va| 国产欧美亚洲国产| 午夜视频国产福利| 天堂中文最新版在线下载 | av在线天堂中文字幕| 日韩亚洲欧美综合| 国产 一区 欧美 日韩| 日韩国内少妇激情av| 神马国产精品三级电影在线观看| 日本猛色少妇xxxxx猛交久久| 国语对白做爰xxxⅹ性视频网站| 久久精品久久久久久久性| 国产探花在线观看一区二区| 国产午夜福利久久久久久| 国产精品久久久久久精品电影| 欧美zozozo另类| 麻豆乱淫一区二区| 亚洲,欧美,日韩| 亚洲精品aⅴ在线观看| 日日撸夜夜添| 成人无遮挡网站| 内地一区二区视频在线| 亚洲成人av在线免费| 日韩精品有码人妻一区| 成人无遮挡网站| 日本三级黄在线观看| 九草在线视频观看| www.色视频.com| 午夜福利视频精品| 精品人妻一区二区三区麻豆| 欧美日韩一区二区视频在线观看视频在线 | 91aial.com中文字幕在线观看| 一级毛片电影观看| 麻豆精品久久久久久蜜桃| 久热这里只有精品99| 免费观看在线日韩| 777米奇影视久久| 男人和女人高潮做爰伦理| 97精品久久久久久久久久精品| 欧美3d第一页| 狠狠精品人妻久久久久久综合| 91精品国产九色| 免费黄色在线免费观看| 美女高潮的动态| 欧美高清性xxxxhd video| 精品一区二区三卡| 99re6热这里在线精品视频| av天堂中文字幕网| 亚洲成人中文字幕在线播放| 精品久久久久久久末码| 亚洲人成网站高清观看| 身体一侧抽搐| 有码 亚洲区| 亚洲国产色片| 欧美日本视频| 国产真实伦视频高清在线观看| 美女xxoo啪啪120秒动态图| 免费黄色在线免费观看| 中国美白少妇内射xxxbb| 只有这里有精品99| 一二三四中文在线观看免费高清| av黄色大香蕉| 边亲边吃奶的免费视频| 国产精品女同一区二区软件| 成人鲁丝片一二三区免费| 久久久久久久久久成人| 久久99热这里只有精品18| 亚洲精品影视一区二区三区av| 欧美高清性xxxxhd video| 黄色日韩在线| 亚洲av福利一区| 久久久久网色| 欧美老熟妇乱子伦牲交| 久久精品久久久久久噜噜老黄| 欧美zozozo另类| 欧美成人午夜免费资源| 制服丝袜香蕉在线| 91久久精品国产一区二区三区| 秋霞在线观看毛片| 欧美激情国产日韩精品一区| 3wmmmm亚洲av在线观看| 日本一二三区视频观看| 国产在线男女| 亚洲熟女精品中文字幕| 一区二区三区免费毛片| 国产大屁股一区二区在线视频| 色网站视频免费| 97超视频在线观看视频| 久久精品夜色国产| 搞女人的毛片| 久久久久久伊人网av| 国产在视频线精品| 精品99又大又爽又粗少妇毛片| 精品久久久久久久人妻蜜臀av| 免费大片黄手机在线观看| 夜夜看夜夜爽夜夜摸| 亚洲av日韩在线播放| 日本午夜av视频| 色哟哟·www| 六月丁香七月| 国产亚洲一区二区精品| av黄色大香蕉| 久久久久久久久久成人| 国产精品一区二区在线观看99| 插逼视频在线观看| 精华霜和精华液先用哪个| 欧美精品国产亚洲| 亚洲成色77777| 男人舔奶头视频| 国产成年人精品一区二区| 中文乱码字字幕精品一区二区三区| 国产伦精品一区二区三区四那| 汤姆久久久久久久影院中文字幕| 精品少妇久久久久久888优播| 狠狠精品人妻久久久久久综合| 老司机影院毛片| xxx大片免费视频| 最近2019中文字幕mv第一页| 91狼人影院| 看非洲黑人一级黄片| 性插视频无遮挡在线免费观看| 极品少妇高潮喷水抽搐| 汤姆久久久久久久影院中文字幕| 王馨瑶露胸无遮挡在线观看| a级一级毛片免费在线观看| 婷婷色麻豆天堂久久| 久久ye,这里只有精品| 在线 av 中文字幕| 高清毛片免费看| 青春草视频在线免费观看| 涩涩av久久男人的天堂| 精品一区二区三卡| 久久久久久久久久久免费av| 日日啪夜夜撸| 五月玫瑰六月丁香| 欧美日韩视频高清一区二区三区二| 亚洲人与动物交配视频| 天天一区二区日本电影三级| av又黄又爽大尺度在线免费看| 国产成人aa在线观看| 色5月婷婷丁香| 少妇丰满av| 一本久久精品| 亚洲性久久影院| 国产伦精品一区二区三区视频9| 2021少妇久久久久久久久久久| 午夜免费观看性视频| 国产黄频视频在线观看| 国产亚洲精品久久久com| 日韩欧美一区视频在线观看 | 成年版毛片免费区| 国国产精品蜜臀av免费| 亚洲精品456在线播放app| 黄色日韩在线| 精品久久国产蜜桃| 免费看av在线观看网站| 国产午夜福利久久久久久| 噜噜噜噜噜久久久久久91| 中文字幕亚洲精品专区| 日本免费在线观看一区| 国产黄片视频在线免费观看| 真实男女啪啪啪动态图| 街头女战士在线观看网站| 69av精品久久久久久| 亚洲av成人精品一区久久| 国产亚洲5aaaaa淫片| 亚洲精品影视一区二区三区av| 大又大粗又爽又黄少妇毛片口| 久久99精品国语久久久| 国产亚洲av嫩草精品影院| 亚洲精华国产精华液的使用体验| 亚洲综合色惰| 成人鲁丝片一二三区免费| 久久久久久久久久人人人人人人| 少妇的逼好多水| 干丝袜人妻中文字幕| 99久久九九国产精品国产免费| 国产成年人精品一区二区| 国产乱来视频区| 成人免费观看视频高清| 搡老乐熟女国产| 丝瓜视频免费看黄片| 我的老师免费观看完整版| 免费大片18禁| 日日啪夜夜爽| 91aial.com中文字幕在线观看| 精品久久久精品久久久| 国产毛片在线视频| 国产国拍精品亚洲av在线观看| 99热6这里只有精品| 日日啪夜夜爽| www.av在线官网国产| 国产成人精品久久久久久| 少妇人妻久久综合中文| 国精品久久久久久国模美| 欧美丝袜亚洲另类| 大又大粗又爽又黄少妇毛片口| 亚洲欧美一区二区三区黑人 | xxx大片免费视频| 免费看av在线观看网站| 边亲边吃奶的免费视频| 国模一区二区三区四区视频| 搞女人的毛片| 国产黄a三级三级三级人| 久热久热在线精品观看| 久久久久久久久久人人人人人人| 午夜亚洲福利在线播放| 国产伦精品一区二区三区四那| 又大又黄又爽视频免费| 久久精品久久久久久久性| 高清午夜精品一区二区三区| 啦啦啦中文免费视频观看日本| 天堂网av新在线| 哪个播放器可以免费观看大片| 99久国产av精品国产电影| 久久精品国产亚洲av涩爱| 国产精品国产三级国产av玫瑰| 超碰97精品在线观看| 亚洲欧美一区二区三区黑人 | av在线天堂中文字幕| 在线观看人妻少妇| 亚洲美女视频黄频| 超碰97精品在线观看| 身体一侧抽搐| 亚洲国产精品成人综合色| 国产淫语在线视频| av免费在线看不卡| 日本一二三区视频观看| 国产精品麻豆人妻色哟哟久久| 人妻一区二区av| 欧美精品一区二区大全| 天美传媒精品一区二区| 亚洲人成网站在线观看播放| 国产男人的电影天堂91| 免费黄网站久久成人精品| 久久久亚洲精品成人影院| 天天躁日日操中文字幕| 十八禁网站网址无遮挡 | 亚洲av.av天堂| 天堂俺去俺来也www色官网| 人体艺术视频欧美日本| 久久久久久久精品精品| 欧美日本视频| 欧美另类一区| 男女无遮挡免费网站观看| 欧美bdsm另类| 欧美激情在线99| 一级毛片aaaaaa免费看小| 一区二区av电影网| 国产精品蜜桃在线观看| 日日啪夜夜撸| 狂野欧美激情性bbbbbb| 国产熟女欧美一区二区| 嫩草影院入口| 麻豆成人av视频| 国产成人精品婷婷| 国产高清三级在线| 99热网站在线观看| 狠狠精品人妻久久久久久综合| 久久精品人妻少妇| 亚洲国产色片| 寂寞人妻少妇视频99o| 卡戴珊不雅视频在线播放| 亚洲欧美成人精品一区二区| 久久久精品94久久精品| av专区在线播放| 最近最新中文字幕免费大全7| 美女视频免费永久观看网站| freevideosex欧美| 国产精品熟女久久久久浪| 黄色一级大片看看| av国产精品久久久久影院| 国产在线男女| 男女边吃奶边做爰视频| 免费人成在线观看视频色| 男女啪啪激烈高潮av片| 成人漫画全彩无遮挡| 超碰97精品在线观看| 欧美日本视频| 亚洲成人精品中文字幕电影| 中国国产av一级| 久久久久国产精品人妻一区二区| 天天一区二区日本电影三级| 亚洲自偷自拍三级| tube8黄色片| 久久精品久久久久久噜噜老黄| 99热这里只有是精品在线观看| 制服丝袜香蕉在线| 免费观看的影片在线观看| 欧美激情久久久久久爽电影| 久久精品久久精品一区二区三区| 十八禁网站网址无遮挡 | xxx大片免费视频| 亚洲图色成人| 日韩不卡一区二区三区视频在线| 18禁动态无遮挡网站| 汤姆久久久久久久影院中文字幕| 亚洲欧美一区二区三区国产| 22中文网久久字幕| 久久久久精品性色| 久久久久久九九精品二区国产| 日韩欧美精品免费久久| 亚洲成人av在线免费| 丝袜美腿在线中文| 女人久久www免费人成看片| 少妇裸体淫交视频免费看高清| 久热这里只有精品99| 大码成人一级视频| 免费黄频网站在线观看国产| 国产一区二区在线观看日韩| 热re99久久精品国产66热6| 国产爱豆传媒在线观看| av女优亚洲男人天堂| 夜夜看夜夜爽夜夜摸| 王馨瑶露胸无遮挡在线观看| av天堂中文字幕网| 亚洲国产最新在线播放| h日本视频在线播放| 深爱激情五月婷婷| 91精品一卡2卡3卡4卡| 男女边吃奶边做爰视频| av专区在线播放| 国产男人的电影天堂91| 日韩电影二区| 国产一区二区亚洲精品在线观看| 亚洲av中文字字幕乱码综合| 亚洲精品国产成人久久av| 一级毛片aaaaaa免费看小| 五月玫瑰六月丁香| 亚洲三级黄色毛片| 一边亲一边摸免费视频| av专区在线播放| 久久久国产一区二区| 亚洲av成人精品一区久久| 高清日韩中文字幕在线| 国产精品麻豆人妻色哟哟久久| 国产成人免费观看mmmm| 又黄又爽又刺激的免费视频.| 久久久久久伊人网av| 亚洲不卡免费看| 久久精品久久久久久噜噜老黄| 国产精品蜜桃在线观看| 视频区图区小说| 香蕉精品网在线| 亚洲自拍偷在线| 免费在线观看成人毛片| 亚洲真实伦在线观看| 少妇人妻精品综合一区二区| 日韩av在线免费看完整版不卡| 亚洲最大成人手机在线| av.在线天堂| 插阴视频在线观看视频| www.av在线官网国产| 99久久精品热视频| 亚洲成人久久爱视频| 久久精品熟女亚洲av麻豆精品| 大话2 男鬼变身卡| 在线观看国产h片| 伊人久久国产一区二区| 身体一侧抽搐| 搡老乐熟女国产| 涩涩av久久男人的天堂| 免费看a级黄色片| 97热精品久久久久久| 久久亚洲国产成人精品v| 色哟哟·www| 国产伦精品一区二区三区四那| 色哟哟·www| 久久这里有精品视频免费| 一级a做视频免费观看| 青春草视频在线免费观看| 2018国产大陆天天弄谢| 小蜜桃在线观看免费完整版高清| 久久影院123| 国产毛片在线视频| 久久97久久精品| 欧美激情久久久久久爽电影| 好男人视频免费观看在线| 黄色日韩在线| 99久久精品国产国产毛片| 肉色欧美久久久久久久蜜桃 | 神马国产精品三级电影在线观看| 欧美3d第一页| 久久久久精品久久久久真实原创| 精品亚洲乱码少妇综合久久| 丰满人妻一区二区三区视频av| 观看免费一级毛片| 少妇的逼好多水| 80岁老熟妇乱子伦牲交| 在线观看国产h片| 亚洲性久久影院| 精品少妇久久久久久888优播| 天天一区二区日本电影三级| 18禁在线无遮挡免费观看视频| 亚洲成人av在线免费| 亚洲欧美清纯卡通| 欧美精品人与动牲交sv欧美| 亚洲国产精品成人久久小说| 亚洲av国产av综合av卡| 超碰97精品在线观看| 免费观看a级毛片全部| 亚洲国产成人一精品久久久| 久久99热这里只有精品18| 中文字幕人妻熟人妻熟丝袜美| 免费不卡的大黄色大毛片视频在线观看| 最近最新中文字幕大全电影3| 国产精品国产三级国产av玫瑰| 极品少妇高潮喷水抽搐| 97人妻精品一区二区三区麻豆| 熟女人妻精品中文字幕| 大香蕉久久网| 国产成人福利小说| 超碰97精品在线观看| 纵有疾风起免费观看全集完整版| 亚洲,一卡二卡三卡| 人人妻人人看人人澡| 国产欧美另类精品又又久久亚洲欧美| 中文字幕亚洲精品专区| 亚洲激情五月婷婷啪啪| 成人鲁丝片一二三区免费| 亚洲精品成人久久久久久| 国产久久久一区二区三区| 我的女老师完整版在线观看| 男人添女人高潮全过程视频| 久久影院123| 国产一区二区三区综合在线观看 | 别揉我奶头 嗯啊视频| 最近最新中文字幕免费大全7| 成人综合一区亚洲| 婷婷色综合大香蕉| 成人亚洲精品av一区二区| 精品久久久久久电影网| 成人综合一区亚洲| 欧美性感艳星| 女人久久www免费人成看片| 一区二区三区乱码不卡18| 日本-黄色视频高清免费观看| 人妻少妇偷人精品九色| 少妇人妻一区二区三区视频| 国产在线一区二区三区精| 黑人高潮一二区| 久久这里有精品视频免费| 18禁在线播放成人免费| 黄色日韩在线| av女优亚洲男人天堂| 汤姆久久久久久久影院中文字幕| 成人亚洲精品av一区二区| 亚洲av男天堂| 国产在线男女| 天美传媒精品一区二区| 久久久久久久久久久免费av| 日韩亚洲欧美综合| 自拍偷自拍亚洲精品老妇| 狠狠精品人妻久久久久久综合| av免费在线看不卡| 少妇的逼水好多| 天美传媒精品一区二区| 美女脱内裤让男人舔精品视频| 亚洲国产精品成人久久小说| 国产午夜精品久久久久久一区二区三区| 国产91av在线免费观看| av又黄又爽大尺度在线免费看| 国产午夜精品一二区理论片| 丝袜美腿在线中文| 国产精品一区二区性色av| 人妻系列 视频| 男女下面进入的视频免费午夜| 极品少妇高潮喷水抽搐| 久久久久久久午夜电影| 岛国毛片在线播放| 日韩亚洲欧美综合| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美精品一区二区大全| 国产视频内射| 中文欧美无线码| 九九爱精品视频在线观看| 七月丁香在线播放| 亚洲人成网站在线观看播放| 日日摸夜夜添夜夜爱| 欧美xxⅹ黑人| 干丝袜人妻中文字幕| 最新中文字幕久久久久| 国产一区有黄有色的免费视频| 国国产精品蜜臀av免费| 精品酒店卫生间| 国模一区二区三区四区视频| 人妻夜夜爽99麻豆av| 日本黄色片子视频| 美女国产视频在线观看| 国产色爽女视频免费观看| 欧美少妇被猛烈插入视频| 亚洲最大成人中文| 国产亚洲av嫩草精品影院| 汤姆久久久久久久影院中文字幕| 99久久人妻综合| 精品久久国产蜜桃| 亚洲欧美日韩东京热| 美女cb高潮喷水在线观看| 午夜日本视频在线| 国产爽快片一区二区三区| 久久久国产一区二区| 色5月婷婷丁香| 午夜老司机福利剧场| 久久久亚洲精品成人影院| 国产一级毛片在线| 日韩中字成人| 大又大粗又爽又黄少妇毛片口| 青青草视频在线视频观看| av福利片在线观看| 一级黄片播放器| 亚洲怡红院男人天堂| a级一级毛片免费在线观看| 欧美日韩视频高清一区二区三区二| 国产av码专区亚洲av| 国产成人aa在线观看| 在线观看一区二区三区| 成人一区二区视频在线观看| 亚洲色图综合在线观看| 欧美精品一区二区大全| 精品99又大又爽又粗少妇毛片| 高清av免费在线| 国产欧美亚洲国产| 久久久久精品性色| 日韩精品有码人妻一区| 日韩欧美一区视频在线观看 | 亚洲最大成人av| 三级经典国产精品| 亚洲成人久久爱视频| 日韩精品有码人妻一区| 80岁老熟妇乱子伦牲交| 久久午夜福利片| 国产黄色免费在线视频| 97在线人人人人妻| 青春草亚洲视频在线观看| .国产精品久久| 一级黄片播放器| 免费av毛片视频| 午夜免费男女啪啪视频观看| 欧美成人a在线观看| 在线观看av片永久免费下载| 亚洲精品乱码久久久v下载方式| 国产精品三级大全| 网址你懂的国产日韩在线| 午夜爱爱视频在线播放| 视频区图区小说| 波多野结衣巨乳人妻| 久久午夜福利片| av在线观看视频网站免费| 观看免费一级毛片| 久久久久久久精品精品| 交换朋友夫妻互换小说| 亚洲不卡免费看| 噜噜噜噜噜久久久久久91| 熟女电影av网| 美女视频免费永久观看网站| 国产av不卡久久| 99久久九九国产精品国产免费| 国产成人aa在线观看| 久久久成人免费电影| 尾随美女入室| 久久久色成人| 99热国产这里只有精品6| 99久久精品一区二区三区| 国产成人a∨麻豆精品| 婷婷色麻豆天堂久久| 在线观看人妻少妇| 最近中文字幕高清免费大全6| 亚洲精品aⅴ在线观看| 亚洲美女搞黄在线观看| 久久影院123| 永久网站在线| 久久久久久久国产电影| 狂野欧美激情性xxxx在线观看| 最近中文字幕2019免费版| 欧美丝袜亚洲另类| 亚洲av福利一区| 成人免费观看视频高清| 久久精品国产亚洲网站| 爱豆传媒免费全集在线观看| 人妻制服诱惑在线中文字幕| 欧美xxxx性猛交bbbb| 日日撸夜夜添| 伊人久久精品亚洲午夜| 狂野欧美激情性bbbbbb| 日韩大片免费观看网站| 国产一区二区三区综合在线观看 | 18+在线观看网站| 欧美日韩国产mv在线观看视频 | 91在线精品国自产拍蜜月| 边亲边吃奶的免费视频| 国产乱人偷精品视频| 综合色丁香网| 午夜福利高清视频| 国产大屁股一区二区在线视频| 久久综合国产亚洲精品| 交换朋友夫妻互换小说| 高清在线视频一区二区三区| 黄色怎么调成土黄色| 久久热精品热| 免费黄频网站在线观看国产| 亚洲,一卡二卡三卡| 日日啪夜夜撸| 精品一区二区三卡| videossex国产| 亚洲天堂av无毛| 日韩一区二区视频免费看| 欧美成人一区二区免费高清观看| 深夜a级毛片| 日韩欧美一区视频在线观看 | 久久韩国三级中文字幕| 自拍偷自拍亚洲精品老妇| 一级毛片aaaaaa免费看小| 国语对白做爰xxxⅹ性视频网站| 精品久久久久久久久亚洲| 免费av观看视频|