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

    Accretion onto Charged Black Holes in Einstein and Massive Theories of Gravity?

    2019-07-16 12:29:34AbdulJawadandShahidChaudhary
    Communications in Theoretical Physics 2019年6期

    Abdul Jawad and Shahid Chaudhary

    1Department of Mathematics,COMSATS University,Islamabad,Lahore Campus-54000,Pakistan

    2Department of Mathematics,Sharif College of Engineering and Technology,Lahore-54000,Pakistan

    AbstractThe accretion process is being investigated onto some important black holes such as Born-Infeld-AdS black hole,non-linear charged black hole solution in AdS space-time and Einstein-Yang-Mills massive gravity in the presence of Born-Infeld nonlinear electrodynamics.We find out the relations of radial velocity,energy density and change of mass for mention black holes and analyze their behavior graphically for different values of equation of state parameters ω.We also examine the relations for critical speed for these black holes.It is observed that for different state parameters different fluids exhibit different evolutions in black holes backgrounds.The energy density of some fluids is negative or positive near the black hole while other fluids become cause to increase or decrease in black hole mass.

    Key words:accretion,black holes,Einstein gravity,massive gravity

    1 Introduction

    The process in which black holes(BHs)or stars like massive astrophysical objects[1]can take particles from its vicinity and increase their mass is known as accretion.[2]It is one of the most worldwide processes in the universe.In fact accretion in some inhomogeneous regions of gas and dust form stars and planets.[3?5]It can be possible that presence of super-massive BHs at centers of some galaxies may be due to accretion process.Due to the accretion of dust or matter from nearby area for long period of time results the formation of huge massive giant like BH.[6]The small compact objects form the astrophysical jets and this shows the presence of huge hot dust in shape of accretion disk around jet formation.[7?10]It is not necessary that accretion always increases the mass of compact objects in some cases in falling matter move away in the form of jet or cosmic rays.[11]Most probably accretion is a dynamic process,velocity and the energy density depend upon time and position.It is really very important to learn about the accretion of different types of dark energies onto BHs because the universe is dominated by dark energy(DE).[12]

    Astronomical results such as cosmic microwave background radiations,Supernova Type Ia[13]and large scale structure data[14?16]provide evidence that the universe expanded exponentially which is known as in fl ation.Infl ation has remarkable importance in cosmology and it is found that dark energy is responsible for this accelerated expansion.Repulsive gravitational force produces due to dark energy and it breaks the null energy conditions and weak energy conditions(WEC).[17]Approximately 71.4 percent energy of the universe is DE.In order to understand the behavior of this DE is still an open problem and it has become one of the most challenging problems in theoretical physics.[18]Many theorists have been trying to handel the challenging problems.Many ideas like cosmological constant,phantom energy,quintessence,kessence,dynamic scalar fields[19]and others have been provided for this problem.DE is modeled by using equation of state,p= ωρ where p and ρ are the pressure and energy density of perfect fluid.Here ω is the equation of state parameter with the following cases of models ω= ?1,?1< ω< ?1/3 and ω< ?1 represents cosmological constant,quintessence and phantom models.[20]

    In the solution of Einstein field equations,presence of essential singularities in different BHs is a prime issue in general relativity.[21]In order to solve this issue regular BHs(RBHs)have been constructed.The metrics of these BHs have no essential singularity hence are regular everywhere.[22]Some of the RBHs satisfy weak energy conditions[23]while strong energy conditions fails to hold.The necessary conditions for RBHs is to satisfy the weak energy conditions having de sitter center.[24]RBHs have an important role to understand gravitational collapse.Bardeen[25]introduced the first RBH known as Bardeen BH,and WEC hold for it.By utilizing Newtonian gravity Bondi[26]studied accretion for some important compact objects.After the discovery of general relativity,the first person who discuss the accretion process for the Schwarzschild BH was Michel.[27]Generalizing the work of Babichev et al.,Debnath[28]studied static accretion onto general static spherically symmetric BH.Bahamonde1 and Jamil[29]studied more general ansatz for a static spherically symmetric spacetime.Jawad and Umair[30]studied accretion onto well known BHs.Kim et al.[31]worked on the accretion on static BH and Kerr-Newman BH of DE.Sharif and Abbas[32]contributed their work for accretion on stringy charged BHs due to phantom energy.

    Recently,lot of work is done on BHs in massive gravity.Hendi and Momennia[33]discussed thermodynamic description of AdS BHs in Born-Infeld(BI)massive gravity with a non-abelian hair.Hendi et al.[34]studied Einstein-BI-massive gravity AdS-BH solutions and their thermodynamical properties.Cai et al.[35]worked on thermodynamics of BH in massive gravity.Myung et al.[35]contributed their work on thermodynamics and phase transitions in the BI-AdS BH.Cembranos et al.[36]studied thermodynamic analysis of non-linear Reissner-Nordstrm BH.Meng et al.[37]found out the BH solution of Einstein-BIYang-Mills theory.Sarioglu[38]studied stationary lifshitz BH of new massive gravity.Sumeet et al.[39]discussed BTZ-BH in massive gravity.

    We discuss the accretion process on BI-AdS BH in massive gravity,non-linear charged BH solution in AdS spacetime and BH solution in Einstein-Yang-Mills(EYM)massive gravity in the presence of Binonlinear electrodynamics.This paper is organized as follows:In Sec.2,we provide the generalized relations for spherically static accretion process.In Sec.3,we discuss accretion onto BI-AdS BH in massive gravity.In Sec.4,we find out different parameters of accretion for non-linear charged BH solution in AdS spacetime.In Sec.5,we discuss accretion for EYM-massive gravity in the presence of Binonlinear electrodynamics and finally in the end we conclude our results.

    2 Generalized Relations for Accretion

    The general relation for static spherical symmetry line element can be written as

    here E(r)>0,F(r)>0,and G(r)>0 are function of r only.The equation of energy momentum tensor is given by

    here ρ,p,and uμare energy density,pressure and four velocity respectively.However,the four velocity can be defined as follows

    Due to restrictions of spherical symmetry,other components of velocity like uθand u?become zero.Here τ is proper time and one can see pressure,energy density,and four velocity are functions of r only.[40]The condition of normalization satisfies uμuμ= ?1,and one can get

    where u=dr/dτ=urand one can see that due to the presence of square root the term utcan be negative or positive which is backward or forward in time conditions.[41]With the help of Michels theory,[42]the proper dark energy model for spherical symmetric BH,[43]can be obtained.Here the equation of continuity leads to

    C0is after integration.

    By projecting the conservation law onto the 4-velocity uμTμν;ν=0,thereby contracting all indices,we can find the equation of continuity

    Using equation of state p=p(ρ),the above equation will take the following form

    here prime is derivative w.r.t r.Taking into account the last equation,we can get

    where C1appears as integration constant.Using Eqs.(5)and(8)one can get

    Here C3is another constant term obtained from C0and C1.However,the mass flux can be defined as

    where C2is also constant of integration.Using Eqs.(5)and(10),we can get

    where C4depends on C1and C2.The differentials of Eqs.(10)and(11)lead to

    The points where the flow speed is equal to speed of sound during the accretion of fluid on the BH is called critical points.Along the particle trajectories the fluid moves toward the BH with increasing velocity.For any critical point r=rc,we have the following possibilities:

    (i)u2=Vs2at r=rc,u2rcand u2>Vs2for r

    (ii)u2>Vs2and u2

    (iii)u2=Vs2for all the values of rrc.

    It is not possible that flow speed inside and outside the critical points is same.It is also a non-physical case.Thus the solution(1)is only valid physical case.In order to obtain turn-around point,one or the other of the bracketed factors in Eq.(12)is terminated,for this case,we get two values in either r or u.The solutions are passing through a critical point that assembles the material falling out(or lf owing into)and along with the particle trajectory the object has monotonically increasing velocity.In order to obtain r=rcwhere the critical point of accretion is located the two bracketed terms(the coefficients of dr and du)in Eq.(12)need to be set zero.

    Now we define the new variable[27]

    As we know critical point at which speed of sound equal to flow speed located at r=rc,thus we have

    Eq.(12)takes the following form

    Here ucrepresents the critical speed of flow at r=rc.[43]Initially the speed of the fluid that moves towards the BH is less than the speed of sound but its speed may transit to a supersonic level as it comes closer to the BH horizons.The region where flow speed is equal to the speed of sound is called a sound horizon.Inside(outside)the sound horizons the flow speed is supersonic(subsonic).The flow speed is supersonic and less than the speed of light inside the sound horizon r

    For E(r)=F(r),Eq.(16)leads to

    The speed of sound at r=rccan be given as

    as both u2cand V2ccan not be negative so following condition must be hold

    The rate of change of mass of BH is given by[44]

    Here derivative w.r.t time is represented by dot and one can see that for the condition ρ+p>0 the mass of BH will increase and accretion happens at exterior of BH.[32]On the other hand,for the fluid with condition ρ+p<0,the mass will decrease.As in case of accretion,mass of BH increases and it decreases in case of Hawking radiations[45]so mass cannot remain constant.

    To discuss the accretion for different barotropic fluids,we take C1=1+ω and C2=1.Using this we must have C3=C2C4/(1+ω).

    3 Born-Infeld-AdS Black Hole in Massive Gravity

    We start by considering a 4-dimensional BI-AdS BH,which has the following action

    where F=(1/4)Fμνμν,R is scalar curvature,G and Λ are gravitational and cosmological constants respectively.Cosmological constant is related to AdS radius as Λ=?3/l2and b represents BI parameter.[46]The coefficients of line element of the BI-AdS BH are given by[47]

    and the value of speed of light c and gravitational constant G are set to be 1,2F1is the hypergeometric function.Also,M and Q represent the mass and charge of BH,respectively.The best thing about the above relation is that one can see in the limit b→ ∞,Q=0,the solution reduces to the Reissner-Nordstrom-AdS BH[48]and the limit Q ?→ 0,it becomes Schwarzchild AdS BH.The electrostatic potential difference between the horizon and infinity is given as

    The relations for radial velocity and the energy density for BI-AdS BH for the case of barotropic fluid[49]having equation p(r)= ωρ(r)are given by

    Thus,the relation for change of the mass of BI-AdS BH take the form

    Figure 1 shows the velocity trajectories w.r.t x=r/M.As ω =1,0,?1 represents sti ff,dust and cosmological constant respectively.We discuss six different cases for ω = ?2,?1.5,?0.5,0,0.5,1.The trajectories of velocity profile are negative for ω = ?2,?1.5 and positive for remaining values of ω.The fluid is at rest at x=4.2,2.2 for ω = ?1.5,?0.5,and ω = ?2,0 respectively.

    Figure 2 shows the energy density trajectories of fluid in surrounding areas of BH w.r.t x=r/M.We plot for six different values of ω = ?2,?1.5,?0.5,0,0.5,1.The trajectories of energy density are negative for ω = ?2,?1.5 and positive for other values of ω.The WEC and dominated energy conditions(DEC)satisfied by dust,sti ffand quintessence fluids.When phantom fluid(ω = ?1.5,?2)moves towards BH then energy density decreases and reverse will happen for dust,sti ffand quintessence fluids(ω =?0.5,0,0.5,1).

    Fig.1 (Color online)Plot of velocity versus x=r/M for BIAdS BH in massive gravity.Other constants are β=1,q=1.055M,M=1,and C4=0.45.

    Fig.2 (Color online)Plot of energy density versus x=r/M for BI-AdS BH in massive gravity.The values of other constants are β=1,q=1.055M,M=1,and C4=0.45.

    Figure 3 shows the change of mass trajectories w.r.t x=r/M.For small values of x,near the BH,0 showing that mass of BH increases for ω =0 and ω = ?0.5 being matter and quintessence respectively.In other words,mass of the BH will increase for matter and quintessence accretions.<0 for phantom-like equations of state such as ω=?1.5,?2 thus mass of BH will decrease by the accretion of phantom like fluids.Figure 4 shows the critical speed of flow against r+.u2cis zero for r+=1 and positive for other values of r+.The critical speed of flow is increasing function of r+as expected.The critical speed of flow is also satisfied the condition mentioned in Eq.(19).

    Fig.3(Color online)Rate of change of mass of BI-AdS BH in massive gravity versus x=r/M.We use β=1,M=1,q=1.055M,and C4=0.45.

    Fig.4(Color online)Plot of critical speed of flow for BI-AdS BH in massive gravity versus r+.We use M=1;b=?1;q=1.055;l=1.

    4 Non-linear Charged BH in AdS Spacetime

    The Einstein gravity coupled to a non-linear electromagnetic field[50]in the AdS space time is described by the action

    where R is the scalar curvature,l is the curvature radius of the AdS spacetime,andL (F)is the non-linear electrodynamic term,which is a function of the invariant Fμν=?μAν??νAμto be the field strength of the non linear electromagnetic field.The non-linear electrodynamic is explicitly defined as[51]

    where M and Q are mass and charge of the system.For this BH,the metric coefficients are defined as

    Similarly for this BH the relations for radial velocity and energy density[50]become

    In this case,the change of mass takes the following form

    Fig.5(Color online)Plot of velocity versus x=r/M for nonlinear BH in AdS spacetime.Other constants are β=1,q=1.055M,M=1,and C4=0.45.

    Fig.6 (Color online)Plot of energy density versus x=r/M for non-linear BH in AdS spacetime.We use β=1,q=1.055M,M=1,and C4=0.45.

    Figure 5 shows the velocity trajectories w.r.t x=r/M.We discuss six different cases for ω = ?2,?1.5,?0.5,0,0.5,1.Trajectories of velocity profile are negative for ω = ?2,?1.5 and positive for remaining valves of ω.Fluid is at rest at x=6.8,3.4 for ω = ?1.5,?0.5 and ω =2,0 respectively.Figure 6 shows the energy density trajectories of fluid in surrounding areas of BH w.r.t x=r/M.Trajectories of energy density are negative for ω = ?2,?1.5 and positive for other values of ω.The WEC and DEC satisfied by dust,sti ffand quintessence fluids.When phantom fluid(ω = ?1.5,?2)approaches to BH,energy density decreases and reverse will happen for dust,sti ffand quintessence fluids(ω=?0.5,0,0.5,1).

    Fig.7 (Color online)Plot of rate of change of mass versus x=r/M for non-linear BH in AdS spacetime.Other values are β=1,q=1.055M,M=1,and C4=0.45.

    Fig.8 (Color online)Plot of critical speed of flow versus r+for non-linear BH in AdS spacetime.Other values are M=1;b=?1;q=1.055;l=1.

    Figure 7 shows the change of mass trajectories w.r.t x=r/M.In short we find that mass of the BH will increase where energy conditions hold while it decreases for fluids violating the same energy conditions.Figure 8 shows the critical speed of flow against r+.u2cis zero at r+=0.While it is positive and increasing function of r+as expected.The critical values for this case are given as

    5 Black Holes Solution in Einstein-Yang-Mills Massive Gravity in the Presence of Born-Infeld Nonlinear Electrodynamics

    We consider the following(3+1)-dimensional action of EYM-massive gravity[50]with BI-NED for the model

    whereLBI(zM)and zYM=Tr(Fμ(aν)F(a)μν)are,respectively,the Lagrangian of BI-NED and the YM invariant.[52]In addition,m is related to the graviton mass while f refers to an auxiliary reference me√tric.Moreover,ciare some free constants andμiare symmetric polynomials of the eigenvalues of 4×4 matrix Kνμ=gμσfσν,we can obtain the following metric functions[53]

    Here M is the total mass of BH,β is linear electrodynamics,magnetic parameter ν is a non-vanishing integer,c is a constant,2F1is hypergeometric function,Q is an integration constant,which is related to the total electric charge of the BH.The fourth term in above equation is related to the magnetic charge(hair),the fi fth term is related to the massive gravitons,and finally,the last term comes from the nonlinearity of electric charge.The important relations of radial velocity and energy density for the model of EYM-massive gravity with BI-NED are given as

    The change of mass by using above general relation yields

    Fig.9 (Color online)Plot of velocity versus x=r/M for EYM-massive gravity in the presence of Binonlinear electrodynamics.Other constant are β=1,q=1.055M,M=1,and C4=0.45.

    Fig.10 (Color online)Plot of energy density versus x=r/M for EYM-massive gravity in the presence of Binonlinear electrodynamics for β=1,M=1,q=1.055M,and C4=0.45.

    Fig.11(Color online)Rate of change of mass of EYM-massive gravity in the presence of Binonlinear electrodynamics versus x=r/M for β=1,M=1,q=1.055M,and C4=0.45.

    Fig.12 (Color online)Plot of critical speed of flow versus r+for EYM-massive gravity in the presence of Binonlinear electrodynamics.Other values are M=1;b=?1;q=1.055;l=1.

    Figure 9 shows the velocity trajectories w.r.t x=r/M.We discuss six different cases for ω = ?2,?1.5,?0.5,0,0.5,1.Trajectory of velocity profile is negative for ω = ?2 and positive for ω = ?1.5,?0.5,0,0.5,1.Figure 10 shows the energy density trajectories of fluid in surrounding areas of BH w.r.t x=r/M.Trajectories of energy density are negative for ω = ?2,?1.5 and positive for other values of ω.The WEC and DEC satisfied by dust,sti ffand quintessence fluids.When phantom fluid(ω = ?1.5,?2)approaches to BH,energy density decreases and reverse will happen for dust,sti ff and quintessence fluids(ω = ?0.5,0,0.5,1).Figure 11 shows the change of mass trajectories w.r.t x=r/M.Mass of the BH increases for matter and quintessence accretions.When˙M<0,the mass of BH will decrease by the accretion of phantom like fluids for ω = ?1.5,?2.Figure 12 shows the critical speed of flow against r+.u2cis positive and increasing function of r+as expected.

    Finally the critical relations for the above mention massive gravity become

    6 Concluding Remarks

    In this paper,we have discussed the general formalism for the study of spherical accretion onto well known BHs in massive gravity.We studied various features of BI-AdS BH in massive gravity,non-linear charged BH in AdS spacetime,EYM-massive gravity in the presence of Binonlinear electrodynamics.We have followed the work of Jawad and Umair[30]and find out the relations of radial velocity,energy density,rate of change of mass for the above mention BHs.We have plotted these parameters for six different values of the equation of state parameter.In order to find out these quantities we took barotropic EoS and obtained correspondence in law of conservation and the barotropic EoS.We established result from different trajectories of u that radial velocity is negative for phantom-like fluids and it is positive for other three fluids mention above.

    The velocity profile for BI-AdS BH in massive gravity are negative for ω = ?2,?1.5 and positive for remaining values of ω.The fluid is at rest at x=4.2,2.2 for ω =?1.5,?0.5 and ω =?2,0 respectively.For non-linear charged BH in AdS spacetime the velocity profile are negative for ω = ?2,?1.5 and positive for remaining values of ω.The fluid is at rest at x=6.8,3.4 for ω = ?1.5,?0.5 and ω=2,0 respectively.While for EYM-massive gravity velocity profile is negative for ω = ?2 and positive for ω = ?1.5,?0.5,0,0.5,1.Other same constants are β =1,q=1.055M,M=1,and C4=0.45.

    For the energy density of BI-AdS BH the WEC and DEC satisfied by dust,sti ffand quintessence fluids.When phantom fluid(ω = ?1.5,?2)moves towards BH then energy density decreases and reverse will happen for dust,sti ffand quintessence fluids(ω = ?0.5,0,0.5,1). For the case of non-linear charged BH in AdS spacetime the WEC and DEC satisfied by dust,sti ffand quintessence fluids.When phantom fluid(ω = ?1.5,?2)approaches to BH,energy density decreases and reverse will happen for dust,sti ffand quintessence fluids(ω = ?0.5,0,0.5,1).

    While for EYM-massive gravity the WEC and DEC satisfied by dust,sti ffand quintessence fluids.When phantom fluid(ω = ?1.5,?2)approaches to BH,energy density decreases and reverse will happen for dust,sti ffand quintessence fluids(ω = ?0.5,0,0.5,1).

    The rate of change of mass for small values of x,near the BI-AdS BH,˙M>0 showing that mass of BH increases for ω =0 and ω = ?0.5 being matter and quintessence respectively.˙M<0 for phantom-like equations of state such as ω = ?1.5,?2 thus mass of BH will decrease by the accretion of phantom like fluids.For the case of non-linear charged BH we find that mass of the BH will increase where energy conditions hold while it decreases for fluids violating the same energy conditions.The mass of EYM-massive gravity increases for matter and quintessence accretions.When˙M<0,the mass of BH will decrease by the accretion of phantom like fluids for ω = ?1.5,?2.Other constants we use are β=1,M=1,q=1.055M,and C4=0.45.

    For the inward and outward flows the conditions that are not allowed are u<0 and u>0 respectively.Similarly we evaluated that the value of energy density is negative for phantom-like fluid near RBHs while it is positive for other models of fluids.Increase and decrease in mass of BH depend upon fluids nature which accretes onto it.We discussed dust,sti ffmatter,quintessence and phantom dark energy accretion onto BH.It is found that different fluids with distinct state parameters have different evolutions in the BH backgrounds.Certain fluids acquire positive or negative energy density near the BH while some fluids become the cause of BH masses to increase or decrease.We have found that the rate of change of mass of our considered BHs increases for dust and sti ffmatter,quintessence-like fluid since these fluids do not have enough repulsive force.However,the mass decreases in the presence of phantom-like fluid(and the corresponding energy density and radial velocity becomes negative)because it has strong negative pressure.

    一级毛片精品| 女人高潮潮喷娇喘18禁视频| a在线观看视频网站| 久久天堂一区二区三区四区| 亚洲久久久国产精品| 欧美日本中文国产一区发布| 国产精品九九99| 岛国毛片在线播放| 亚洲国产毛片av蜜桃av| 欧美精品亚洲一区二区| 少妇猛男粗大的猛烈进出视频| 精品国产乱子伦一区二区三区| 在线播放国产精品三级| 国产成人欧美在线观看 | 色综合婷婷激情| 久久国产精品人妻蜜桃| 国产高清国产精品国产三级| 亚洲精品乱久久久久久| 久久人妻福利社区极品人妻图片| 热99久久久久精品小说推荐| 人人妻,人人澡人人爽秒播| 不卡av一区二区三区| 久久青草综合色| 久久久久久久久免费视频了| 一级片'在线观看视频| 成人三级做爰电影| 成年人黄色毛片网站| 中亚洲国语对白在线视频| 高清毛片免费观看视频网站 | 在线观看66精品国产| 脱女人内裤的视频| 高清在线国产一区| 国产成人啪精品午夜网站| 制服人妻中文乱码| 成人影院久久| 黑人欧美特级aaaaaa片| 亚洲成a人片在线一区二区| 在线观看免费日韩欧美大片| 天天添夜夜摸| 亚洲av日韩精品久久久久久密| 精品电影一区二区在线| 亚洲色图av天堂| 9色porny在线观看| 精品欧美一区二区三区在线| 777米奇影视久久| aaaaa片日本免费| 国产熟女午夜一区二区三区| 亚洲一区二区三区欧美精品| 大型av网站在线播放| 国产黄色免费在线视频| 日韩欧美免费精品| 99国产综合亚洲精品| 91精品三级在线观看| 亚洲av成人av| 五月开心婷婷网| 老汉色av国产亚洲站长工具| 日韩视频一区二区在线观看| 在线观看一区二区三区激情| 国产高清视频在线播放一区| 国产乱人伦免费视频| 咕卡用的链子| 久久久久精品国产欧美久久久| 人妻久久中文字幕网| 欧美日韩亚洲国产一区二区在线观看 | 怎么达到女性高潮| 两人在一起打扑克的视频| 久久久久国内视频| 久久人妻熟女aⅴ| 中文字幕人妻丝袜一区二区| 日韩精品免费视频一区二区三区| 韩国精品一区二区三区| 麻豆乱淫一区二区| 亚洲va日本ⅴa欧美va伊人久久| 色综合婷婷激情| 怎么达到女性高潮| 久久热在线av| 成人18禁在线播放| 国产激情久久老熟女| 国产97色在线日韩免费| 麻豆国产av国片精品| a级毛片在线看网站| 久久人妻福利社区极品人妻图片| videosex国产| 亚洲成av片中文字幕在线观看| 国产高清激情床上av| 涩涩av久久男人的天堂| 亚洲专区中文字幕在线| 麻豆成人av在线观看| 一个人免费在线观看的高清视频| 国产高清激情床上av| 精品亚洲成国产av| 韩国精品一区二区三区| 国产99白浆流出| 久热爱精品视频在线9| 老司机深夜福利视频在线观看| 真人做人爱边吃奶动态| 久久国产精品人妻蜜桃| 免费av中文字幕在线| 女人被狂操c到高潮| 中文字幕人妻熟女乱码| 男人操女人黄网站| 1024视频免费在线观看| 人成视频在线观看免费观看| 黄色视频,在线免费观看| 亚洲成av片中文字幕在线观看| 久久精品国产清高在天天线| 精品乱码久久久久久99久播| 天堂动漫精品| 99在线人妻在线中文字幕 | 国产片内射在线| 久久久久精品国产欧美久久久| 制服人妻中文乱码| 亚洲欧美日韩高清在线视频| 日日夜夜操网爽| 色播在线永久视频| 欧美国产精品va在线观看不卡| 久热爱精品视频在线9| 免费看a级黄色片| 国产精品二区激情视频| 久久中文字幕人妻熟女| 黑丝袜美女国产一区| 电影成人av| 欧美亚洲 丝袜 人妻 在线| 国产成人精品无人区| 欧美精品一区二区免费开放| 国产av又大| 国产欧美日韩一区二区精品| 亚洲人成伊人成综合网2020| 最新在线观看一区二区三区| 中出人妻视频一区二区| 国产精品免费大片| 国产精品一区二区在线不卡| 人人澡人人妻人| 久久精品国产综合久久久| 9热在线视频观看99| 妹子高潮喷水视频| 一进一出抽搐gif免费好疼 | 欧美性长视频在线观看| 水蜜桃什么品种好| 少妇 在线观看| 12—13女人毛片做爰片一| 乱人伦中国视频| 久久精品aⅴ一区二区三区四区| 两个人免费观看高清视频| 国产精品久久久av美女十八| 欧美av亚洲av综合av国产av| av在线播放免费不卡| 欧美老熟妇乱子伦牲交| 亚洲视频免费观看视频| 亚洲精品粉嫩美女一区| 日本a在线网址| 久久精品成人免费网站| 久久人妻熟女aⅴ| 国内毛片毛片毛片毛片毛片| 亚洲性夜色夜夜综合| 国产精品自产拍在线观看55亚洲 | 亚洲五月婷婷丁香| a级毛片在线看网站| 18在线观看网站| 国产成人一区二区三区免费视频网站| 久久久久国产一级毛片高清牌| 久久久国产一区二区| 国产男女内射视频| 国产片内射在线| 免费观看精品视频网站| 黄色毛片三级朝国网站| 香蕉久久夜色| 激情在线观看视频在线高清 | 国产精品久久视频播放| 女警被强在线播放| 十分钟在线观看高清视频www| 欧美一级毛片孕妇| 久久国产精品人妻蜜桃| 黄色毛片三级朝国网站| 亚洲一区二区三区不卡视频| 成人永久免费在线观看视频| 又黄又粗又硬又大视频| 欧美激情久久久久久爽电影 | av视频免费观看在线观看| 国产主播在线观看一区二区| 精品国产一区二区三区四区第35| 久久国产乱子伦精品免费另类| 在线视频色国产色| 巨乳人妻的诱惑在线观看| 国产又色又爽无遮挡免费看| 日韩成人在线观看一区二区三区| 少妇 在线观看| 久久久精品区二区三区| 久久久久久久国产电影| 亚洲美女黄片视频| 一级a爱片免费观看的视频| 国产成人系列免费观看| 欧美在线黄色| 一区二区日韩欧美中文字幕| 国产精品免费一区二区三区在线 | 天堂中文最新版在线下载| 亚洲国产精品一区二区三区在线| 51午夜福利影视在线观看| 亚洲av熟女| 超色免费av| 国产av又大| 18禁裸乳无遮挡动漫免费视频| 高清av免费在线| 最近最新中文字幕大全电影3 | 这个男人来自地球电影免费观看| 欧美精品一区二区免费开放| 别揉我奶头~嗯~啊~动态视频| 99精品久久久久人妻精品| 亚洲成a人片在线一区二区| 色综合欧美亚洲国产小说| 国产精品欧美亚洲77777| 国产日韩一区二区三区精品不卡| 脱女人内裤的视频| 欧美乱码精品一区二区三区| 美女视频免费永久观看网站| 午夜福利欧美成人| 欧美激情 高清一区二区三区| 黄色视频不卡| 大香蕉久久网| 日韩视频一区二区在线观看| 性色av乱码一区二区三区2| 国产成人精品无人区| 18禁裸乳无遮挡免费网站照片 | 他把我摸到了高潮在线观看| 一二三四社区在线视频社区8| 国产精品.久久久| 交换朋友夫妻互换小说| 国产黄色免费在线视频| 免费观看a级毛片全部| 日本撒尿小便嘘嘘汇集6| 免费一级毛片在线播放高清视频 | 在线av久久热| 欧美日韩瑟瑟在线播放| 久久精品aⅴ一区二区三区四区| www.999成人在线观看| 日韩熟女老妇一区二区性免费视频| 精品国产一区二区三区久久久樱花| 夜夜爽天天搞| 一级作爱视频免费观看| 99久久综合精品五月天人人| av天堂久久9| 在线观看免费高清a一片| 热99国产精品久久久久久7| 国产成人系列免费观看| 捣出白浆h1v1| 91av网站免费观看| 午夜福利在线观看吧| 一区二区三区精品91| 99久久国产精品久久久| 亚洲第一欧美日韩一区二区三区| 欧美久久黑人一区二区| 在线国产一区二区在线| 亚洲一区中文字幕在线| 女人被躁到高潮嗷嗷叫费观| 精品一品国产午夜福利视频| 国产午夜精品久久久久久| 在线看a的网站| 后天国语完整版免费观看| 日本撒尿小便嘘嘘汇集6| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品久久午夜乱码| 国产精品乱码一区二三区的特点 | 精品福利永久在线观看| 老司机午夜十八禁免费视频| 一级毛片高清免费大全| 不卡一级毛片| 久久久精品免费免费高清| 99国产综合亚洲精品| 69精品国产乱码久久久| 久久久久国产一级毛片高清牌| 国产人伦9x9x在线观看| 亚洲va日本ⅴa欧美va伊人久久| 成人国产一区最新在线观看| 久热这里只有精品99| 久久亚洲真实| 很黄的视频免费| 99精品在免费线老司机午夜| 精品国产亚洲在线| 国产不卡一卡二| 69av精品久久久久久| 中出人妻视频一区二区| 岛国在线观看网站| 精品久久久久久,| 午夜福利在线免费观看网站| x7x7x7水蜜桃| 久久国产亚洲av麻豆专区| 亚洲三区欧美一区| 日韩 欧美 亚洲 中文字幕| 99精国产麻豆久久婷婷| 欧美激情高清一区二区三区| 日本一区二区免费在线视频| 久久久久精品国产欧美久久久| 色老头精品视频在线观看| 91字幕亚洲| 亚洲熟妇熟女久久| 1024视频免费在线观看| 欧美一级毛片孕妇| 不卡av一区二区三区| 国产精品影院久久| 在线天堂中文资源库| a在线观看视频网站| 男女高潮啪啪啪动态图| 老司机深夜福利视频在线观看| 嫩草影视91久久| 国产精品二区激情视频| 这个男人来自地球电影免费观看| 正在播放国产对白刺激| 手机成人av网站| 在线视频色国产色| 黄色视频,在线免费观看| 国产精品美女特级片免费视频播放器 | 国产男女超爽视频在线观看| 欧美亚洲 丝袜 人妻 在线| av网站在线播放免费| 热re99久久国产66热| 日韩中文字幕欧美一区二区| 国产精品二区激情视频| 99riav亚洲国产免费| 欧美精品av麻豆av| 在线观看www视频免费| 男女床上黄色一级片免费看| 亚洲人成电影免费在线| 日日爽夜夜爽网站| 日本五十路高清| 日韩熟女老妇一区二区性免费视频| 18禁裸乳无遮挡免费网站照片 | 亚洲国产看品久久| 老鸭窝网址在线观看| www.熟女人妻精品国产| 欧美日韩av久久| 亚洲视频免费观看视频| 成人黄色视频免费在线看| 人人澡人人妻人| 每晚都被弄得嗷嗷叫到高潮| 国产精品九九99| 热99re8久久精品国产| 黄片播放在线免费| 亚洲精品国产区一区二| 岛国毛片在线播放| 久99久视频精品免费| 色播在线永久视频| 搡老岳熟女国产| 狠狠狠狠99中文字幕| 亚洲第一欧美日韩一区二区三区| netflix在线观看网站| 亚洲精品在线观看二区| 精品国产乱码久久久久久男人| 大码成人一级视频| 水蜜桃什么品种好| 欧美乱妇无乱码| 黄色a级毛片大全视频| 飞空精品影院首页| 国产一区二区三区在线臀色熟女 | 久久性视频一级片| 9色porny在线观看| 黑人巨大精品欧美一区二区mp4| 亚洲美女黄片视频| 中文字幕另类日韩欧美亚洲嫩草| √禁漫天堂资源中文www| 老司机在亚洲福利影院| 国产91精品成人一区二区三区| 亚洲人成电影免费在线| 天天躁日日躁夜夜躁夜夜| 亚洲片人在线观看| 亚洲精品自拍成人| 日本vs欧美在线观看视频| 一级毛片女人18水好多| 国产一区二区三区视频了| 精品亚洲成a人片在线观看| 久久精品人人爽人人爽视色| 操出白浆在线播放| 黑人欧美特级aaaaaa片| 看片在线看免费视频| 在线十欧美十亚洲十日本专区| 亚洲国产毛片av蜜桃av| 久久亚洲真实| 国产aⅴ精品一区二区三区波| 老熟妇仑乱视频hdxx| 久久99一区二区三区| 欧美一级毛片孕妇| 精品福利观看| 久久久精品免费免费高清| 最近最新免费中文字幕在线| 精品国产乱子伦一区二区三区| 欧美在线黄色| 精品无人区乱码1区二区| 国产亚洲一区二区精品| 中文字幕另类日韩欧美亚洲嫩草| 法律面前人人平等表现在哪些方面| 每晚都被弄得嗷嗷叫到高潮| 一个人免费在线观看的高清视频| 精品一区二区三区四区五区乱码| 黄片播放在线免费| 91字幕亚洲| 69av精品久久久久久| 久久ye,这里只有精品| 999精品在线视频| 90打野战视频偷拍视频| 精品久久蜜臀av无| 97人妻天天添夜夜摸| 人人妻人人爽人人添夜夜欢视频| 欧美色视频一区免费| 黑人操中国人逼视频| 午夜福利一区二区在线看| 中文欧美无线码| 色综合婷婷激情| 欧美精品高潮呻吟av久久| 最近最新中文字幕大全免费视频| 美女扒开内裤让男人捅视频| 少妇 在线观看| 99国产精品免费福利视频| 亚洲人成电影免费在线| 丰满人妻熟妇乱又伦精品不卡| 啦啦啦视频在线资源免费观看| 日韩欧美国产一区二区入口| 日韩欧美免费精品| 美女扒开内裤让男人捅视频| 少妇 在线观看| 十八禁网站免费在线| 在线观看日韩欧美| 精品卡一卡二卡四卡免费| 欧美av亚洲av综合av国产av| 亚洲黑人精品在线| 久久人妻av系列| 亚洲少妇的诱惑av| 久久精品成人免费网站| 久久青草综合色| 波多野结衣av一区二区av| 18禁裸乳无遮挡免费网站照片 | 18在线观看网站| netflix在线观看网站| 欧美午夜高清在线| 久久这里只有精品19| 亚洲欧美日韩另类电影网站| 亚洲人成77777在线视频| 大码成人一级视频| 欧美日韩黄片免| 中文字幕另类日韩欧美亚洲嫩草| av一本久久久久| 亚洲av成人不卡在线观看播放网| 大型av网站在线播放| 中文字幕最新亚洲高清| 日韩成人在线观看一区二区三区| 91国产中文字幕| 中文字幕精品免费在线观看视频| 人妻一区二区av| 成年女人毛片免费观看观看9 | av视频免费观看在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久人人97超碰香蕉20202| 99精国产麻豆久久婷婷| 国产xxxxx性猛交| 高清视频免费观看一区二区| 国产麻豆69| a级毛片黄视频| 欧美日本中文国产一区发布| 日本vs欧美在线观看视频| 777久久人妻少妇嫩草av网站| 亚洲第一青青草原| av片东京热男人的天堂| 在线观看66精品国产| 在线观看美女被高潮喷水网站 | 在线观看午夜福利视频| 波多野结衣巨乳人妻| 99久久精品热视频| 欧美成狂野欧美在线观看| 香蕉久久夜色| 五月玫瑰六月丁香| 蜜桃亚洲精品一区二区三区| 亚洲精品在线美女| 女人十人毛片免费观看3o分钟| 人妻丰满熟妇av一区二区三区| 亚洲国产精品999在线| 最近在线观看免费完整版| 亚洲va日本ⅴa欧美va伊人久久| 久久99热这里只有精品18| 一进一出抽搐动态| 久久久久国产精品人妻aⅴ院| 麻豆国产97在线/欧美| 黄色丝袜av网址大全| 国产野战对白在线观看| 精品一区二区三区视频在线 | 亚洲内射少妇av| 丰满人妻熟妇乱又伦精品不卡| 神马国产精品三级电影在线观看| 亚洲精品久久国产高清桃花| 欧美绝顶高潮抽搐喷水| 天堂√8在线中文| svipshipincom国产片| 午夜免费成人在线视频| 国产亚洲av嫩草精品影院| 色老头精品视频在线观看| 久久久久久国产a免费观看| 国产aⅴ精品一区二区三区波| 午夜福利在线在线| 男人的好看免费观看在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 18禁在线播放成人免费| 国产熟女xx| or卡值多少钱| 中文字幕久久专区| 国产欧美日韩精品亚洲av| 久久久久久久亚洲中文字幕 | 少妇丰满av| 免费在线观看成人毛片| 日本黄大片高清| 中亚洲国语对白在线视频| 欧美+日韩+精品| 露出奶头的视频| 中文字幕久久专区| 久久久久久久久大av| 国产成人啪精品午夜网站| 无遮挡黄片免费观看| 久久性视频一级片| 精品福利观看| 午夜福利免费观看在线| 99久久成人亚洲精品观看| 国产野战对白在线观看| 一级作爱视频免费观看| 国产一区二区三区视频了| 嫁个100分男人电影在线观看| 激情在线观看视频在线高清| 国产高清有码在线观看视频| 韩国av一区二区三区四区| 国产一区二区亚洲精品在线观看| 高潮久久久久久久久久久不卡| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av免费高清在线观看| 国产乱人伦免费视频| 69人妻影院| 两性午夜刺激爽爽歪歪视频在线观看| 国产高清有码在线观看视频| 免费av观看视频| 国产国拍精品亚洲av在线观看 | 欧美日本视频| 麻豆国产97在线/欧美| 国产在线精品亚洲第一网站| 草草在线视频免费看| 久久精品国产99精品国产亚洲性色| 日韩欧美 国产精品| 男人的好看免费观看在线视频| 变态另类成人亚洲欧美熟女| 欧美bdsm另类| 又爽又黄无遮挡网站| 男女午夜视频在线观看| 国产精品美女特级片免费视频播放器| 精品一区二区三区人妻视频| 男插女下体视频免费在线播放| 在线十欧美十亚洲十日本专区| 性色av乱码一区二区三区2| 香蕉久久夜色| 综合色av麻豆| 午夜日韩欧美国产| 国产av不卡久久| 亚洲在线观看片| 他把我摸到了高潮在线观看| 熟女少妇亚洲综合色aaa.| 亚洲自拍偷在线| 美女高潮的动态| 亚洲欧美日韩无卡精品| 国产单亲对白刺激| 男女下面进入的视频免费午夜| 老司机福利观看| 在线视频色国产色| 亚洲人成伊人成综合网2020| 男人舔奶头视频| 日韩亚洲欧美综合| 欧美日韩黄片免| 亚洲成av人片免费观看| 国产午夜福利久久久久久| 免费观看的影片在线观看| 国产精品亚洲一级av第二区| 国产精品 国内视频| 白带黄色成豆腐渣| 免费人成视频x8x8入口观看| 成人国产综合亚洲| 欧美zozozo另类| 亚洲av二区三区四区| 在线观看午夜福利视频| 亚洲成人免费电影在线观看| 啪啪无遮挡十八禁网站| 国产乱人伦免费视频| 老鸭窝网址在线观看| 国产乱人伦免费视频| 久久天躁狠狠躁夜夜2o2o| 波野结衣二区三区在线 | 两人在一起打扑克的视频| 国产野战对白在线观看| 国产亚洲欧美98| 在线视频色国产色| 99久国产av精品| 3wmmmm亚洲av在线观看| 亚洲男人的天堂狠狠| 草草在线视频免费看| 丰满人妻熟妇乱又伦精品不卡| 久久久久久久久大av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 夜夜躁狠狠躁天天躁| 亚洲熟妇中文字幕五十中出| 亚洲av中文字字幕乱码综合| 我要搜黄色片| 老司机在亚洲福利影院| 午夜老司机福利剧场| 两个人的视频大全免费| 欧美日韩福利视频一区二区| 91字幕亚洲| 老司机午夜十八禁免费视频| 午夜福利18| 日韩欧美在线二视频| 精品乱码久久久久久99久播| 欧美最黄视频在线播放免费| 乱人视频在线观看| 欧美绝顶高潮抽搐喷水| 91久久精品国产一区二区成人 | 欧美日本亚洲视频在线播放|