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

    脈沖星磁輻射制動力矩對具有磁輻射的兩成分模型自旋的長期減速*1

    2016-08-01 11:44:22李林森
    天文研究與技術 2016年3期
    關鍵詞:脈沖星

    李林森

    (東北師范大學物理學院,吉林 長春 130024)

    ?

    脈沖星磁輻射制動力矩對具有磁輻射的兩成分模型自旋的長期減速*1

    李林森

    (東北師范大學物理學院,吉林 長春130024)

    摘要:利用分析法研究了具有磁輻射的兩成分模型(殼層和中子超流體)的脈沖星在可變的磁輻射制動力矩的作用下,兩成分自旋角速度隨時間的長期變化。給出了具有磁輻射兩成分模型的耦合方程組的分析解。理論結果給出兩成分模型在外力可變的磁輻射制動力矩的作用下,自旋角速度隨時間長期減慢。利用所得的分析解對具有磁輻射的兩成分模型蟹狀星云脈沖星(PSR0531+21)(Crab)在磁輻射力矩可變的情況下做了數(shù)值計算。并討論了所得的理論和數(shù)值結果。結果表明,蟹狀星云脈沖星(PSR0531+21)在磁輻射制動力矩的作用下,殼層自旋角速度隨時間的長期減速每年為-0.245 s。

    關鍵詞:脈沖星;具有磁輻射的兩成分模型;磁輻射力矩;自旋減速

    為了解釋脈沖星突然加速的現(xiàn)象,許多學者對脈沖星提出各式各樣的模型,其中較成功的模型如磁偶極輻射模型、四極彈性能模型以及兩成分模型等,兩成分模型最先由文[1-3]提出。兩成分模型主要是由導電較高的固體外殼和殼內的中子超流體的混合物組成。因此中子星的構造模型顯然分為兩部分:一部分是帶電的固體外殼;另一部分是殼內的中子超流體的混合體。此外,兩成分以不同角速度自轉。自轉使超流體成分同帶電成分在粘滯作用下相耦合,盡管這種耦合是弱的,但也反應出超流體的豐富度和耦合的程度。由于脈沖星磁偶極輻射在脈沖星表層產生磁輻射制動力矩,這種力矩使兩成分自旋角速度改變。因此,只要知道磁轉矩變化以及可觀測的宏觀馳豫時間τ和中子超流體的豐富度Q,兩成分的自旋耦合的微分方程可解。文[4]給出了對兩成分耦合方程組在磁輻射制動力矩不變的情況下耦合方程組的解并對解做了討論。文[5]研究了中子星的發(fā)射噪聲的兩成分模型, 但沒有研究磁輻射力矩對兩成分模型角速度改變的影響。文[6]研究中子星兩成分模型在廣義相對論構架內的自轉動力方程,并假定殼成分的轉動Ωc為常數(shù),而中子超流體層的角速度依賴于坐標,但沒有研究磁輻射力矩對兩成分模型角速度改變的影響。本文作者研究了在磁輻射制動力矩隨時間可變的情況下,具有磁輻射的兩成分角速度變化的規(guī)律,并給出方程組的解。最后將解應用于具有磁輻射兩成分模型的蟹狀星云脈沖星(Crab pulsar)的兩成分角速度隨時間變化的數(shù)值解,此解不同于文[4]給出的解。

    1脈沖星兩成分的耦合方程式及其在磁輻射力矩不變條件下的解

    根據兩成分模型理論,兩成分的自旋角速度并不一致。假定帶電的固體外殼的轉動慣量為IC并以角速度Ω繞自轉軸旋轉,而殼內中子超流體的轉動慣量In并以角速度Ωn繞自轉軸旋轉,但Ω≠Ωn,帶電外殼成分與殼內中子超流體成分之間的耦合由下式給出[4]:

    (1)

    (2)

    (1)式、(2)式是在星震后不存在躍變時的方程組。其中,N(t)為外部輻射制動力矩;τc是理論上計算的脈沖星弛豫時間或稱微觀弛豫時間, 它依賴于Ω。文[7]給出計算τc的公式,文[8]將τc用下式表示:

    (3)

    其中,Δ為超流體間隙參數(shù)(能量間隔);EF為電子的費米能量。

    微觀弛豫時間τc可用宏觀弛豫時間τ表示,宏觀弛豫時間τ是可以觀測的時間,兩者之間關系由下式表示[8]:

    (4)

    中子超流體的豐富度Q:

    (5)

    其中,ΔΩn為Ωn的初始躍變,而總轉動慣量I:

    (6)

    文[8]指出只要Ic?In,Q可以表示:

    (7)

    由于大多數(shù)脈沖星均滿足上述條件,故:

    (8)

    因此,(4)式可以寫成:

    (9)

    文[4]給出了在磁輻射制動力矩N不變的情況下,耦合方程組(1)~(2)的解:

    (10)

    以下本文給出在輻射制動力矩N可變的情況下,方程組(1)~(2)的解并應用于脈沖星。

    2可變磁輻射制動轉矩的形式

    本文研究具有磁輻射的兩成分模型的脈沖星。有磁輻射必然有磁輻射力矩作用在具有磁輻射的兩成分模型上,所以方程組(1)~(2)中N是磁輻射力矩。本文給出磁輻射力矩隨時間變化的公式。

    按磁偶極輻射模型,脈沖星的輻射功率Wd是由自轉能的變率E·轉化來的,即

    (11)

    輻射功率和自轉能的變率可寫成[ 9 ]:

    (12)

    由(11)和(12)式可得

    (13)

    其中,磁矩μ=R3Bssinα;R為脈沖星半徑;Bs為脈沖星的表面磁場。假定磁轉矩垂直于自轉軸,則α=900,故磁矩μ=R3Bs。如果在短時間內不考慮磁衰減或增長,則(13)式成為

    (14)

    積分式的下限t0是在星震后不存在躍變時具有磁輻射的兩成分模型統(tǒng)一體的初始時間。這個初始時間t0不是脈沖星誕生開始t=0時的初始時間,它是具有磁輻射的兩成分模型統(tǒng)一體從現(xiàn)在年齡開始的時間。J(t0)是t0時的角動量:

    所以:J(t)m=J(t0)m[1+K1(t-t0)]-1/2,

    (15a)

    (15b)

    利用(15a)和(15b)可以給出磁輻射制動力矩N(t)的表達式:

    將(15a)式代入上式后得

    (16)

    3磁輻射制動力矩可變情況下兩成分耦合方程組(1)~(2)的分析解

    下形式:

    (17)

    (18)

    這樣,方程組(1)~(2)可用觀測到的宏觀弛豫時間τ表示。將方程組(1)和(2)的兩端各項相加,可得

    轉動慣量Ic和In不變或為常數(shù)的情況下,將(16)式的N(t)代入上式積分后可得

    所以: J(t)m=IcΩ+InΩn=J(t0)m[1+K1(t-t0)]-1/2.

    在鐘揚排得密密麻麻的時間表里,學生的事總是優(yōu)先的。他還是很多中小學生喜愛的明星專家、“科學隊長”,心甘情愿將大量寶貴的時間分給科普。他說:“小時候家中那套殘缺不全的《十萬個為什么》讓我相信,科學能深入兒童心靈?!?/p>

    (19)

    此式正好是(15a)式,這說明角動量變化的(19)式左端的J(t)m是兩成分角動量之和。

    利用(7)~(8)式將(19)改成下列兩式:

    (20)

    (21)

    將(21)式的Ωn代入(17)式,方程組(17)~(18)變成下列一組一階非線性齊次方程組:

    (22)

    (23)

    解(22)式得到Ω后再代入(23)式得到Ωn的解。

    同樣,將(20)的Ω代入(18)式,方程組(17)~(18)式變成下列另一種形式:

    (24)

    (25)

    解的形式為

    令:Ω=e-t/τ(I1+I2)+C1e-t/τ.

    (26)

    積分式:

    I1和I2代入(26)式,得

    (27)

    (28)

    利用初始條件t=t0,Ω(t0)=Ω(t0)m,這表示磁輻射和兩成分統(tǒng)一體在星震和躍變后從現(xiàn)在輻射年齡t0開始, 則由上式可得

    所以殼層角速度Ω在磁輻射力矩作用下隨時間變化是

    (29)

    (30)

    在(29)式中,t-t0=(t0+Δt)-t0=Δt,Δt是以年為單位的時間間隔, 它與脈沖星年齡無關, 也與起始時間無關,所以 (29)式可以寫成用時間間隔表示的式子:

    (31)

    此式只與取的時間間隔有關,而與時間起點無關。故本文在計算時取時間間隔為一年。

    將(29)式代入(23)式后可得到Ωn的一階線性非齊次方程, 令其解為

    (32)

    利用前面的積分方法可得I1、I2、I3、I4各積分結果和所得常數(shù)C,將其代入上式,由此得到中子超流體的角速Ωn隨時間變化的解(推導式子較長,因篇幅關系略去):

    (33)

    (34)

    (29)式和(33)式就是本文給出的脈沖星磁輻射力矩對兩成分模型自轉角速隨時間長期變化的式子。

    4脈沖星PSR0531+21(Crab)的理論數(shù)值結果

    本文研究的脈沖星必須具有磁輻射的兩成分模型,根據文[3],蟹狀星云脈沖星(Crab, PSR0531+21)是兩成分模型,根據文[4]它又是磁輻射模型或是具有這兩種模型的統(tǒng)一體。本文選取Crab脈沖星(PSR0531+21)作為計算實例。首先給出這個脈沖星的物理參數(shù)如表1,其中Ω0、Q、τ引自文[3, 10],表面磁場Bs引自文[11],μ=R3Bs。一般假定中子星的半徑R=1.2 × 106cm,轉動慣量I=1.4 × 1045(cm3·g)[6]

    表1 蟹狀星云脈沖星(PSR0531+21)的物理數(shù)據

    將表1中的數(shù)據代入(28)式,得到PSR0531+21的K1=5.43 × 10-11,K2=4.33 × 10-4,K3=2.35 × 10-8,再取時間間隔Δt=1年。 將K1、K2、K3的數(shù)值代入(30)和(31)式,得到PSR0531+21(Crab)的殼層在磁輻射力矩作用下角速每年變化的數(shù)值如表2。

    表2蟹狀星云脈沖星(PSR0531+21)的數(shù)值的結果(每年變化的數(shù)值: Δt=1年)

    Table 2The numerical results for PSR0531+21 per year (Δt=1 year)

    脈沖星PulsarΩ(t)(rad/s)Ω/Ω0δΩ(rad/s)PSR0531+21(Crab)189.75500.999871-0.2450

    Ω0的數(shù)值在表1中給出。

    對于中子超流體的角速度變化情形由于初始角速度Ω(t0)n難以觀測到, 所以中子超流體的角速度Ω(t)n隨時間變化也同樣難以觀測到。故本文只能根據所推得的理論式子對殼層角速度變化做一計算,但給出中子超流體的角速度變化的理論式子仍有理論價值和意義。

    5討論和結論

    (1)對推得的理論結果式子的驗證。將初始時間t=t0代入理論式子(29)和(33),可得到右端的初始角速度Ω(t0)和Ω(t0)n。即(29)~(30)和(33)~(34)式變成Ω=Ω(t0),δΩ=0;Ωn=Ω(t0)nδΩn=0。這說明文中給理論結果是正確的。此外,(29)和(33)式的級數(shù)隨時間增大是收斂的。

    (2)根據表2給出的數(shù)值結果,脈沖星PSR0531+21在磁輻射外力矩的作用下,外殼角速度隨時間逐漸減慢,然而這種減速是長期的,而和由于殼震角速突然加速的躍變兩者迥然不同。前者是長期性的,后者是突然臨時性的。此外,前者長期減速不會影響后者突然躍變的加速,而后者的突然躍變加速對前者的長期變化也不產生影響。因為方程組(1)~(2)是在星震躍變不存在的情況下成立,因此它的解(長期變化)不受星震和躍變的影響。

    (3)脈沖星磁輻射的演化起點是從脈沖星誕生t=0開始,而兩成分模型的演化起點是從星震和躍變后為起點, 兩者并不一致。 本文致力于如何使兩個模型的演化起點合二為一,這是本文解決此問題的特點。本文研究的脈沖星具有磁輻射和兩成分的統(tǒng)一體(Crab脈沖星就是這兩種模型的統(tǒng)一體),因此,演化的起點t0必須在星震和躍變后統(tǒng)一體的初始時間。如果按磁輻射的初始時間是脈沖星誕生t=0為起始時間,可是脈沖星誕生時的物理參量是一個不確定的較大的物理量。另外也不符合兩成分模型的演化起點。 因此應該采用脈沖星的現(xiàn)在年齡t0為兩個模型統(tǒng)一體的初始時間。 在積分式子(29)~(30)中,取從t0到t為積分上下限, 可是t=t0+Δt,故在(29)~(30)式子中取t-t0=(t0+Δt)-t0=Δt。按本文計算取時間間隔Δt=1年,所以演化時間間隔與脈沖星演化的起始時間無關。這樣,具有磁輻射的兩成分模型,如Crab脈沖星的演化時間可選取在星震和躍變后磁輻射和兩成分的統(tǒng)一體的現(xiàn)在年齡為兩者的初始時間。這個方法解決了上述兩個模型演化起點不一致的問題。

    (4)本文假定(13)式中的磁偶極矩μ為常量,即不隨時間衰減或增長。實際上,如果從長期考慮磁矩μ是隨時間變化的。即μ=μ0e-ξt,ξ=2/τD。文[9]給出的τD=1.6 × 106年,ξ=1.3 × 10-6年。因此,磁衰減時間較長,磁衰減系數(shù)每年很小。本文只計算每年角速的變化值,因此,磁衰減的影響可以不必考慮。

    另一方面,有的脈沖星的磁場是增長的。例如文[12]給出的PSR0531+21(Crab脈沖星)屬于磁場增長的脈沖星,即μ=μ0e+ξt=R3B0e+ξt,但磁場增量很小,增量為+0.000 6 × 1 012 G/年。這樣小的增量對Crab脈沖星每年殼層角速度變化的影響也可不必考慮。

    (5)在第3節(jié)中I1、I2、I3、I4各積分式子利用了下面的積分公式:

    級數(shù)展開式應用于本文n=1/2,a=1/K1τ和x=[1+K(t-t0)]。表1中K=5.43 × 10-11,τ=7.7 d=665 280 s,K1τ=3.612 4 × 10-5,x>1隨時間t延長而增加。將這些代入上式, 此級數(shù)是收斂的。所以,I1、I2、I3、I4各積分式是收斂的,故此所得(29)式也是收斂的。另外,如果將第4節(jié)表1下面的K1、K2、K3的數(shù)值代入(29)式也可得到(29)式是收斂的,當[1+K1(t-t0)]隨時間t延長而增大。但由于在級數(shù)中略去了[1+K1(t-t0)]-7/2以上的高階項,所以,解(29)式是近似的。

    參考文獻:

    [1]Baym G, Pethick C, Pines D, et al. Spin up in neutron stars: the future of the Vela pulsar[J]. Nature, 1969, 224: 872-874.

    [2]Ruderman M. Pulsar: structure and dynamics[J]. Annual Review of Astronomy and Astrophysics, 1972, 10: 427.

    [3]Pines D. Observing neutron stars: information on stellar structure from pulsars and compact X-ray sources[C]// Proceedings of the Sixteenth Solvay Conference on Physics, Brussels, Belgium, September 24-28, 1973. 1974: 147-173.

    [4]Shapiro S L, Teukolsky S A. Black holes, white dwarfs, and neutron stars: the physics of compact objects[M]. New York: John Wiley & Sons, 1983: 247-295.

    [5]Baykal A, Alpar A, Kizilaglu U. A shot noise model for a two-component neutron star[J]. Astronomy & Astrophysics, 1991, 252(2): 664-668.

    [6]Sedrakian D M. Rotation of a two-component model neutron star in GTR[J]. Astrophysics, 1997, 40(3): 260-266.

    [7]Feibelman P J. Relaxation of electron velocity in a rotating neutron superfluid: application to the relaxation of a pulsar’s slowdown rate[J]. Physical Review D, 1971, 4(6): 1589-1597.

    [8]Pines D, Shaham J, Rudeman M A. Physics of dense matter[J]. Proceedings of the International Astronomical Union, 1972, 1053: 21-25.

    [9]曲欽岳, 汪珍茹, 陸埮, 等. 脈沖星的統(tǒng)計分析與JP 1953[J]. 科學通報, 1976(4): 176-177.

    [10]唐小英. 中子星星震與脈沖星加速[J]. 北京天文臺臺刊, 1975(4): 68-81.

    [11]Manchester R N, Hobbs G B, Teoh A, et al. The Australia Telescope National Facility pulsar catalogue[J]. The Astronomical Journal, 2005, 129(4): 1993-2006.

    [12]Li Linsen. A method for judging decay or growth of the magnetic field of pulsar[J]. Journal of Astrophysics & Astronomy, 2009, 30(3): 145-151.

    *收稿日期:2015-12-20;

    修訂日期:2016-02-06

    作者簡介:李林森,男,教授. 研究方向:天體自轉理論. Email: dbsd-lls@163.com

    中圖分類號:P1

    文獻標識碼:A

    文章編號:1672-7673(2016)03-0277-07

    The Impact of Magnetic Radiation Braking Torque on the Secular Retardation of Spin of Two-Components of Pulsar

    Li Linsen

    (School of Physics, Northeast Normal University, Changchun 130024, China, Email: dbsd-lls@163.com)

    Abstract:The secular influences of the variable magnetic radiation braking torque on the spin down of two-components of pulsar are studied by using analytical method. The analytical solutions of the coupled equation system are given with due consideration of the time-dependent magnetic braking torque by solving the system of non-homogeneous equations. The analytical solutions obtained are applied to the research for PSR0531+21 (Crab Pulsar). The numerical results are obtained in Table 2. The numerical results show that the spin of Crab pulsar speeds down -0.245rad/yr. Discussions are held on the obtained results.

    Key words:Pulsar; Two-component model; Magnetic radiation torque; Speed down

    CN 53-1189/PISSN 1672-7673

    猜你喜歡
    脈沖星
    “中國天眼”已發(fā)現(xiàn)740余顆新脈沖星
    軍事文摘(2023年12期)2023-06-12 07:51:00
    航空學報(2023年3期)2023-03-12 08:38:40
    發(fā)現(xiàn)脈沖星的女天文學家——貝爾
    科學(2022年4期)2022-10-25 02:43:42
    這些小眾的脈沖星可能有大作用
    脈沖星方位誤差估計的兩步卡爾曼濾波算法
    脈沖星角位置對脈沖模板的影響及其削弱策略
    宇宙時鐘——脈沖星
    基于虛擬觀測值的X射線單脈沖星星光組合導航
    長征十一號成功發(fā)射脈沖星試驗衛(wèi)星
    太空探索(2016年12期)2016-07-18 11:13:43
    中國即將發(fā)射全世界首顆脈沖星導航試驗衛(wèi)星
    電子世界(2016年19期)2016-03-17 14:16:59
    a级毛片黄视频| 777米奇影视久久| 亚洲精品,欧美精品| 国产有黄有色有爽视频| 欧美日本中文国产一区发布| 欧美97在线视频| 高清毛片免费看| 热re99久久精品国产66热6| a 毛片基地| 国产精品嫩草影院av在线观看| 99久久综合免费| 久久精品久久精品一区二区三区| 色哟哟·www| 国产国拍精品亚洲av在线观看| av免费在线看不卡| 最新中文字幕久久久久| 亚洲激情五月婷婷啪啪| 国产成人精品无人区| 熟女电影av网| 看十八女毛片水多多多| 久久精品国产亚洲av涩爱| 精品一区二区三卡| 老女人水多毛片| 亚洲图色成人| 中文字幕免费在线视频6| 蜜桃国产av成人99| 黄网站色视频无遮挡免费观看| 九色亚洲精品在线播放| 成年人午夜在线观看视频| 大香蕉97超碰在线| 一区二区三区乱码不卡18| 丰满迷人的少妇在线观看| 在线观看三级黄色| 最近最新中文字幕免费大全7| 成人影院久久| 免费av中文字幕在线| 亚洲av免费高清在线观看| 久久 成人 亚洲| 丰满迷人的少妇在线观看| 久久ye,这里只有精品| √禁漫天堂资源中文www| 国产男人的电影天堂91| av免费在线看不卡| 国产精品国产三级国产av玫瑰| 久久久久久伊人网av| 搡女人真爽免费视频火全软件| 国产一区二区在线观看日韩| 制服诱惑二区| 免费播放大片免费观看视频在线观看| 国产精品国产三级专区第一集| 久久免费观看电影| 亚洲av欧美aⅴ国产| 捣出白浆h1v1| 99国产精品免费福利视频| 欧美激情 高清一区二区三区| 亚洲精品一二三| 18+在线观看网站| 久久久精品免费免费高清| 日韩三级伦理在线观看| 国产极品粉嫩免费观看在线| 街头女战士在线观看网站| 亚洲精品视频女| freevideosex欧美| 高清黄色对白视频在线免费看| 丝袜人妻中文字幕| av网站免费在线观看视频| 永久免费av网站大全| 国产免费一级a男人的天堂| 亚洲欧美精品自产自拍| 亚洲欧美一区二区三区黑人 | 午夜精品国产一区二区电影| 亚洲精品456在线播放app| 亚洲精品456在线播放app| 999精品在线视频| 亚洲欧美日韩另类电影网站| 亚洲美女视频黄频| 亚洲国产av新网站| 亚洲国产精品专区欧美| 美女国产视频在线观看| 色吧在线观看| 国产欧美另类精品又又久久亚洲欧美| 色吧在线观看| 一级片'在线观看视频| 午夜久久久在线观看| 男女国产视频网站| 久久女婷五月综合色啪小说| 久久精品人人爽人人爽视色| 久久久久久久久久久免费av| 亚洲人成77777在线视频| 色网站视频免费| 国产成人精品久久久久久| 国产免费一级a男人的天堂| 久久久亚洲精品成人影院| 亚洲精品视频女| 亚洲精品久久久久久婷婷小说| 国产黄色免费在线视频| 亚洲精品自拍成人| 这个男人来自地球电影免费观看 | 国产 精品1| 麻豆精品久久久久久蜜桃| 人妻系列 视频| 国产淫语在线视频| 赤兔流量卡办理| 国产毛片在线视频| 九草在线视频观看| 欧美最新免费一区二区三区| 狠狠婷婷综合久久久久久88av| 日本免费在线观看一区| 成人毛片a级毛片在线播放| 国产黄色免费在线视频| 五月伊人婷婷丁香| 色婷婷av一区二区三区视频| 少妇的丰满在线观看| 国国产精品蜜臀av免费| 青春草亚洲视频在线观看| 国产av国产精品国产| 中文欧美无线码| 欧美人与性动交α欧美精品济南到 | 97在线视频观看| 久久人人爽av亚洲精品天堂| 亚洲,一卡二卡三卡| 人妻人人澡人人爽人人| 日韩免费高清中文字幕av| 制服人妻中文乱码| 久久99热这里只频精品6学生| 国产成人av激情在线播放| av免费在线看不卡| 亚洲国产av新网站| 亚洲精品久久久久久婷婷小说| 国产 一区精品| 亚洲av成人精品一二三区| 在线看a的网站| 只有这里有精品99| 国产欧美另类精品又又久久亚洲欧美| 亚洲综合色网址| 97在线视频观看| 免费观看性生交大片5| 黄网站色视频无遮挡免费观看| 国产成人欧美| 欧美成人午夜免费资源| 国产一区有黄有色的免费视频| 久久毛片免费看一区二区三区| 久久鲁丝午夜福利片| 免费观看av网站的网址| 亚洲av中文av极速乱| 中文精品一卡2卡3卡4更新| 国产精品一区二区在线不卡| 午夜免费鲁丝| 午夜福利影视在线免费观看| av女优亚洲男人天堂| 热re99久久精品国产66热6| 老熟女久久久| 日韩一区二区视频免费看| 亚洲久久久国产精品| 国产男女超爽视频在线观看| 母亲3免费完整高清在线观看 | 尾随美女入室| 国产白丝娇喘喷水9色精品| 男女下面插进去视频免费观看 | 国产精品一区www在线观看| 国产精品久久久久久精品古装| 欧美成人午夜免费资源| 亚洲丝袜综合中文字幕| 亚洲,欧美,日韩| 黑丝袜美女国产一区| 精品人妻偷拍中文字幕| 26uuu在线亚洲综合色| 有码 亚洲区| av在线播放精品| 在现免费观看毛片| 欧美3d第一页| 乱码一卡2卡4卡精品| 九草在线视频观看| 国产精品久久久久久久电影| 综合色丁香网| 午夜福利视频在线观看免费| 亚洲精品久久成人aⅴ小说| 在线精品无人区一区二区三| 精品国产一区二区三区久久久樱花| 久久久亚洲精品成人影院| 欧美日韩精品成人综合77777| 黑人高潮一二区| 亚洲情色 制服丝袜| 中文字幕精品免费在线观看视频 | 亚洲欧美精品自产自拍| 国产 一区精品| 久久久国产一区二区| 国产精品三级大全| 国精品久久久久久国模美| 免费日韩欧美在线观看| 午夜激情av网站| 婷婷色av中文字幕| 少妇精品久久久久久久| 久久亚洲国产成人精品v| 精品国产露脸久久av麻豆| 91国产中文字幕| 一边亲一边摸免费视频| 久久99热这里只频精品6学生| 欧美人与性动交α欧美软件 | 精品人妻偷拍中文字幕| 免费播放大片免费观看视频在线观看| 国产激情久久老熟女| 国产成人免费无遮挡视频| 国产成人aa在线观看| 伊人久久国产一区二区| 尾随美女入室| 成人毛片60女人毛片免费| 少妇被粗大猛烈的视频| 久久人人爽人人爽人人片va| 亚洲,一卡二卡三卡| 亚洲精品久久午夜乱码| 亚洲国产最新在线播放| 亚洲精品中文字幕在线视频| 精品国产一区二区三区四区第35| 久久久久国产网址| 亚洲精品aⅴ在线观看| 免费播放大片免费观看视频在线观看| 欧美少妇被猛烈插入视频| 18禁在线无遮挡免费观看视频| 免费黄色在线免费观看| 免费黄频网站在线观看国产| 99re6热这里在线精品视频| 天天操日日干夜夜撸| 国产成人91sexporn| 国产精品女同一区二区软件| 亚洲精品日本国产第一区| 久久人人97超碰香蕉20202| 人妻 亚洲 视频| 免费观看在线日韩| 26uuu在线亚洲综合色| 中文字幕最新亚洲高清| 国产老妇伦熟女老妇高清| 日韩制服丝袜自拍偷拍| 色5月婷婷丁香| 天天操日日干夜夜撸| 成人国语在线视频| 中文字幕av电影在线播放| 免费黄色在线免费观看| 国产高清三级在线| 久久久国产精品麻豆| 街头女战士在线观看网站| 国产国拍精品亚洲av在线观看| 咕卡用的链子| 欧美成人精品欧美一级黄| a级片在线免费高清观看视频| 99re6热这里在线精品视频| 久久久久久久久久久久大奶| 国产成人免费观看mmmm| 国精品久久久久久国模美| 一级黄片播放器| 午夜影院在线不卡| 久久久久久久久久人人人人人人| 99香蕉大伊视频| 天堂8中文在线网| 最近2019中文字幕mv第一页| 巨乳人妻的诱惑在线观看| 欧美日韩综合久久久久久| 看非洲黑人一级黄片| 天堂8中文在线网| 国产日韩欧美视频二区| 日韩电影二区| 免费观看在线日韩| 精品一区二区免费观看| 精品少妇久久久久久888优播| 秋霞在线观看毛片| 日韩电影二区| 国产免费福利视频在线观看| 桃花免费在线播放| 美女内射精品一级片tv| 熟妇人妻不卡中文字幕| 国产一区二区三区综合在线观看 | 美女国产高潮福利片在线看| 91成人精品电影| 欧美最新免费一区二区三区| 国产极品粉嫩免费观看在线| 久久99热这里只频精品6学生| 久久久精品区二区三区| 国产成人午夜福利电影在线观看| 天天操日日干夜夜撸| 中文精品一卡2卡3卡4更新| 又黄又粗又硬又大视频| 伊人久久国产一区二区| 看免费成人av毛片| 国产精品 国内视频| 一级a做视频免费观看| 国产老妇伦熟女老妇高清| 免费黄频网站在线观看国产| 国产熟女欧美一区二区| 99热网站在线观看| 欧美日韩国产mv在线观看视频| 精品99又大又爽又粗少妇毛片| 不卡视频在线观看欧美| 午夜av观看不卡| 国产精品女同一区二区软件| 国产日韩欧美视频二区| 蜜桃在线观看..| 男女边摸边吃奶| 有码 亚洲区| 九色亚洲精品在线播放| 人妻 亚洲 视频| 91午夜精品亚洲一区二区三区| 国产男女内射视频| 老司机影院成人| 国产xxxxx性猛交| 大话2 男鬼变身卡| 亚洲内射少妇av| 国产精品久久久久久精品电影小说| 精品国产一区二区三区四区第35| 满18在线观看网站| xxx大片免费视频| 在线 av 中文字幕| 成人综合一区亚洲| 成人毛片a级毛片在线播放| 大片电影免费在线观看免费| 色婷婷久久久亚洲欧美| 巨乳人妻的诱惑在线观看| 亚洲,一卡二卡三卡| 18禁动态无遮挡网站| 精品酒店卫生间| 亚洲国产av新网站| 丰满少妇做爰视频| 国产亚洲av片在线观看秒播厂| 夫妻性生交免费视频一级片| 美女国产视频在线观看| 亚洲欧美清纯卡通| 国产色婷婷99| 母亲3免费完整高清在线观看 | 午夜日本视频在线| 91久久精品国产一区二区三区| 狂野欧美激情性bbbbbb| 国产成人午夜福利电影在线观看| 久久99蜜桃精品久久| 综合色丁香网| 中文字幕免费在线视频6| 免费av不卡在线播放| 亚洲av.av天堂| 在现免费观看毛片| 80岁老熟妇乱子伦牲交| 亚洲精品av麻豆狂野| 国产乱来视频区| 精品人妻熟女毛片av久久网站| 99热国产这里只有精品6| 最后的刺客免费高清国语| 国产精品久久久av美女十八| 久久99一区二区三区| 黑人猛操日本美女一级片| 日韩一区二区三区影片| 久久av网站| 成人国产麻豆网| 一本大道久久a久久精品| 久久狼人影院| 丝袜人妻中文字幕| 精品少妇久久久久久888优播| 18禁观看日本| 熟女电影av网| 男女边吃奶边做爰视频| 啦啦啦啦在线视频资源| videosex国产| 成人国产av品久久久| 精品久久久精品久久久| 夜夜爽夜夜爽视频| 婷婷色麻豆天堂久久| 久久精品久久精品一区二区三区| 乱人伦中国视频| 你懂的网址亚洲精品在线观看| 一级a做视频免费观看| 韩国高清视频一区二区三区| 免费大片黄手机在线观看| 久久这里只有精品19| 欧美 亚洲 国产 日韩一| 国产伦理片在线播放av一区| 日日撸夜夜添| 最近2019中文字幕mv第一页| 五月开心婷婷网| xxx大片免费视频| 丰满乱子伦码专区| 一边摸一边做爽爽视频免费| 日本黄色日本黄色录像| 搡女人真爽免费视频火全软件| 日日啪夜夜爽| 日本黄大片高清| 制服人妻中文乱码| 中文字幕亚洲精品专区| 九九在线视频观看精品| 高清av免费在线| 国产无遮挡羞羞视频在线观看| 国产白丝娇喘喷水9色精品| 精品午夜福利在线看| 少妇被粗大的猛进出69影院 | 久久久久人妻精品一区果冻| 狠狠精品人妻久久久久久综合| 国产av码专区亚洲av| 国产亚洲最大av| 亚洲精品国产av蜜桃| 亚洲欧美一区二区三区黑人 | 伦理电影免费视频| 天天躁夜夜躁狠狠躁躁| 亚洲欧美精品自产自拍| 成人18禁高潮啪啪吃奶动态图| 女的被弄到高潮叫床怎么办| 在线天堂最新版资源| 男女边吃奶边做爰视频| 中文字幕免费在线视频6| 高清黄色对白视频在线免费看| 久久人人爽av亚洲精品天堂| 男女午夜视频在线观看 | 成年人午夜在线观看视频| 日韩成人伦理影院| 美女国产视频在线观看| 免费黄色在线免费观看| 久久久久久久久久久免费av| 九色成人免费人妻av| 欧美成人精品欧美一级黄| 日本-黄色视频高清免费观看| 久久ye,这里只有精品| 亚洲av在线观看美女高潮| 91午夜精品亚洲一区二区三区| 啦啦啦视频在线资源免费观看| 国产xxxxx性猛交| 男人操女人黄网站| 亚洲av免费高清在线观看| 免费在线观看完整版高清| 99久久中文字幕三级久久日本| 超碰97精品在线观看| 精品国产一区二区三区久久久樱花| 激情五月婷婷亚洲| 乱人伦中国视频| 亚洲一区二区三区欧美精品| 男男h啪啪无遮挡| av视频免费观看在线观看| 最新中文字幕久久久久| 欧美精品亚洲一区二区| 亚洲中文av在线| 日本猛色少妇xxxxx猛交久久| 国产精品蜜桃在线观看| 国产成人精品无人区| 免费观看性生交大片5| 欧美国产精品一级二级三级| 黑人巨大精品欧美一区二区蜜桃 | 又黄又粗又硬又大视频| 伦精品一区二区三区| 国产成人免费观看mmmm| 韩国高清视频一区二区三区| 一级毛片黄色毛片免费观看视频| 天堂8中文在线网| 日本欧美视频一区| 亚洲熟女精品中文字幕| 久久久久网色| 看十八女毛片水多多多| 国产亚洲最大av| 国产亚洲精品第一综合不卡 | 人妻系列 视频| 免费女性裸体啪啪无遮挡网站| 久久99热这里只频精品6学生| 成年动漫av网址| 菩萨蛮人人尽说江南好唐韦庄| 亚洲五月色婷婷综合| 精品午夜福利在线看| 中文精品一卡2卡3卡4更新| 国产精品久久久久久精品古装| 国产成人精品无人区| 卡戴珊不雅视频在线播放| 国产av一区二区精品久久| 久久久久久伊人网av| 成人二区视频| 男的添女的下面高潮视频| 精品卡一卡二卡四卡免费| 午夜福利,免费看| 又粗又硬又长又爽又黄的视频| 18禁动态无遮挡网站| 91午夜精品亚洲一区二区三区| 久久国内精品自在自线图片| 精品少妇黑人巨大在线播放| 人妻少妇偷人精品九色| 大陆偷拍与自拍| 亚洲,一卡二卡三卡| 国产成人午夜福利电影在线观看| 国产成人免费观看mmmm| 黄色 视频免费看| 男女啪啪激烈高潮av片| 最后的刺客免费高清国语| 国产一区二区三区综合在线观看 | 日韩视频在线欧美| 精品一区二区免费观看| 999精品在线视频| 亚洲精品日韩在线中文字幕| 国产亚洲精品久久久com| 日韩熟女老妇一区二区性免费视频| 青春草亚洲视频在线观看| 亚洲丝袜综合中文字幕| av片东京热男人的天堂| 欧美日韩亚洲高清精品| 亚洲av国产av综合av卡| 国产片特级美女逼逼视频| 亚洲成人av在线免费| xxxhd国产人妻xxx| 亚洲精华国产精华液的使用体验| 五月玫瑰六月丁香| 伊人久久国产一区二区| 美女主播在线视频| 青春草国产在线视频| a级毛片在线看网站| 最近手机中文字幕大全| videosex国产| 久久精品久久精品一区二区三区| 亚洲欧美成人精品一区二区| 日韩三级伦理在线观看| 久久精品国产a三级三级三级| 欧美日韩视频精品一区| 亚洲精品第二区| 国产精品一区二区在线观看99| 久久99热这里只频精品6学生| 亚洲精华国产精华液的使用体验| 午夜激情久久久久久久| 伊人久久国产一区二区| 又大又黄又爽视频免费| 国产精品久久久久久精品古装| 视频在线观看一区二区三区| 精品一区二区三区视频在线| 国产av码专区亚洲av| 毛片一级片免费看久久久久| 久久鲁丝午夜福利片| 国产精品女同一区二区软件| 久久久欧美国产精品| 亚洲综合色网址| 97在线人人人人妻| 99热国产这里只有精品6| 黄色 视频免费看| 黄片无遮挡物在线观看| 亚洲五月色婷婷综合| 天天躁夜夜躁狠狠久久av| 色婷婷久久久亚洲欧美| 免费观看av网站的网址| 人妻人人澡人人爽人人| 大陆偷拍与自拍| av国产久精品久网站免费入址| videosex国产| 美女福利国产在线| 黄色配什么色好看| 一区二区三区四区激情视频| 亚洲性久久影院| 久久午夜综合久久蜜桃| 午夜福利乱码中文字幕| 国产欧美另类精品又又久久亚洲欧美| 自线自在国产av| 少妇被粗大猛烈的视频| 亚洲精品中文字幕在线视频| 丝袜人妻中文字幕| 久久精品夜色国产| 交换朋友夫妻互换小说| 亚洲国产色片| 天天操日日干夜夜撸| 亚洲欧美中文字幕日韩二区| 精品福利永久在线观看| 丝袜脚勾引网站| 男女下面插进去视频免费观看 | 最近中文字幕2019免费版| 精品少妇内射三级| 90打野战视频偷拍视频| av国产久精品久网站免费入址| 波野结衣二区三区在线| 少妇猛男粗大的猛烈进出视频| 女性生殖器流出的白浆| 好男人视频免费观看在线| 久久久久久久久久人人人人人人| 成人无遮挡网站| 中文字幕另类日韩欧美亚洲嫩草| 蜜桃国产av成人99| 色哟哟·www| 日本与韩国留学比较| 韩国高清视频一区二区三区| 欧美精品一区二区免费开放| 免费看av在线观看网站| 91aial.com中文字幕在线观看| 建设人人有责人人尽责人人享有的| 日产精品乱码卡一卡2卡三| 久久人人爽人人片av| 久久99蜜桃精品久久| 亚洲精品中文字幕在线视频| 欧美人与性动交α欧美精品济南到 | 日韩一本色道免费dvd| 欧美精品亚洲一区二区| 国产精品久久久久久久久免| 亚洲天堂av无毛| 熟妇人妻不卡中文字幕| 日韩制服骚丝袜av| 极品少妇高潮喷水抽搐| 国产一区亚洲一区在线观看| 欧美97在线视频| 97在线视频观看| 久久av网站| av在线播放精品| 国产男女超爽视频在线观看| 国语对白做爰xxxⅹ性视频网站| 免费观看在线日韩| 国产片特级美女逼逼视频| 91午夜精品亚洲一区二区三区| 人妻系列 视频| 午夜av观看不卡| 免费黄网站久久成人精品| 久久久久久伊人网av| 国产日韩欧美在线精品| 国产亚洲午夜精品一区二区久久| 国产精品 国内视频| 黄网站色视频无遮挡免费观看| 91aial.com中文字幕在线观看| 18禁动态无遮挡网站| 菩萨蛮人人尽说江南好唐韦庄| videos熟女内射| 国产男人的电影天堂91| 亚洲伊人色综图| 日本色播在线视频| 美女视频免费永久观看网站| 国产国语露脸激情在线看| 亚洲国产精品国产精品|