arxiv: v1 [astro-ph] 7 Mar 2008

Transcription

1 Mon. Not.. Aston. Soc., 5 () Pinted 7th Decembe 8 (MN LATEX style file v.) Analytic solutions to the accetion of a otating finite cloud towads a cental object I. Newtonian appoach. axiv:83.v [asto-ph] 7 Ma 8 S. Mendoa & E. Tejeda Instituto de Astonomía, Univesidad Nacional Autónoma de México, AP 7-64, Distito Fedeal 45, México 7th Decembe 8 INTODUCTION The fist model of a spheical symmetic accetion flow towads a cental object was made by Bondi (95). With time, it has become the key to undestand diffeent accetion pocesses in the univese (see e.g. Fank et al. ()), despite the fact that it was ceated fo mathematical cuiosity, athe than fo astophysical applications (Bondi 5). Bondi thought of an infinite gas cloud with a non elativistic gavitational cental object located at the oigin of coodinates. He assumed the flow to have eached a steady state. Unde these assumptions Bondi was able to integate the hydodynamic equations fo this spheical symmetic flow. Yeas late, Ulich (976) modified Bondi s ideas by assuming all fluid paticles at the bode of the cloud to have a cetain amount of angula momentum, following a distibution given by a igid body otating about the axis. In his model, Ulich took no account of pessue effects on the infalling gas. In othe wods, his analysis was ballistic, which is appoximately tue if the initial specific angula momentum of an infalling paticle is small and if heating by adiation and viscosity effects ae negligible (cf. Mendoa et al. 4). The bounday conditions in Ulich s model imply that all fluid paticles move in paabolas whose focus coespond to the oigin of coodinates. As soon as paticles aive to the equatoial plane they collide with thei symmetic countepats and stay on it, foming an equatoial disc of adius d made of fluid paticles which otate about the cental object. Ulich s model is usually taken as the basic model ABSTACT We constuct a steady analytic accetion flow model fo a finite otating gas cloud that accetes mateial to a cental gavitational object. The pessue gadients of the flow ae consideed to be negligible and so, the flow is ballistic. We also assume a steady flow and conside the paticles at the bounday of the spheical cloud to be otating as a igid body, with a fixed amount of inwads adial velocity. This epesents a genealisation of the taditional infinite gas cloud model descibed by Ulich (976). We show that the fluid paticles ae no longe paabolic, instead they could be any conical section which depends on the specific initial bounday conditions of the flow. Natually, the model pedicts the fomation of an equatoial accetion disc about the cental object with a adius diffeent fom the one calculated by Ulich (976). Key wods: hydodynamics accetion, accetion discs which pedicts the fomation of an accetion disc of paticles obiting a cental object due to accetion with otation. Both models, Bondi (95) and Ulich (976) have been taken to the elativistic egime by assuming the cental object is a Schwaschild black hole (Michel 97; Hueta & Mendoa 7). They have wide applications in high enegy phenomena (see fo example Lee & amie-ui 6; Beloboodov & Illaionov ). In this aticle (and in a futue elativistic one) we assume that the otating gas cloud about a cental object has a finite adius. Hee we assume that the cental object is Newtonian. We show that unde this assumption it is possible to build a cylindical symmetic flow which also pedicts the existence of an accetion disc about the cental object. As it will be shown, the tajectoies ae no longe sections of paabolas. They ae geneal conical sections and once bounday conditions ae fixed they can be of diffeent families fo each fluid paticle. In Sections and 3 we calculate the velocity and density fields of the flow espectively. Sections 4 and 5 pesent two impotant limiting cases of the model, the last one giving ise to Ulich s (976) model. Finally, Section 6 descibes the behaviou of the flow fo it s diffeent paametes and discusses its limitations. Though all the aticle we use a system of units in which Newton s gavitational constant G =.

2 S. Mendoa, E. Tejeda THE MODEL AND ITS VELOCITY FIELD Let us assume that a spheical cloud with adius accetes matte in a steady way towads a cental object located at the oigin of coodinates. All fluid paticles located at ae taken to follow igid body motion in the aimuthal diection, i.e. they otate about the axis of coodinates. At this position, paticles also have a adial velocity component v (hee and in what follows, due to the symmety of the poblem we use spheical coodinates, θ, and φ fo the adial coodinate, and the pola and aimuthal angles). The paticle numbe density n at is taken to be constant with a value n. To simplify the poblem, we assume the self gavity of the gas cloud to be small as compaed to the gavitational potential poduced by the cental object. In addition, if the gadients of pessue ae small compaed to the kinetic and potential enegies, then all fluid paticle tajectoies can be appoximated as ballistic. The total angula momentum of the gas cloud points in the diection of the axis, which contains the cental object. Note that when and v, the model conveges to the accetion poblem descibed by Ulich (976), which is well descibed by Mendoa et al. (4). In fact, Ulich s convegence is mathematically poved in Section 5. Since the gas is being acceted fom, it s accetion flow ate pe unit mass N at any fixed adial distance is constant and given by Ṅ = 4π n v = const. () At, each paticle has a constant specific angula momentum h given by h = h sin θ, whee θ is the initial pola angle of the paticle s tajectoy. Unde all the above assumptions, the specific enegy E and the specific angula momentum h ae constants of motion along each paticula tajectoy. The enegy E is given by E = v + h M = v + h sin θ M. () We now intoduce two dimensionless paametes µ and ν given by the paticle is located at and its pola and aimuthal angles ae given by θ and φ. Ove the plane, the paticle tajectoy is defined by an aimuthal angle ϕ, which at the initial position has the value ϕ. The tajectoy of a given paticle is defined by the solution to Keple s poblem (Landau & Lifshit 989): = µ sin θ ecos ϕ. (5) In the pevious equation and in what follows, unless stated othewise, distances ae measued in units of. The eccenticity e of the obit is given by e = q + 4εµ sin θ. (6) Substitution of = in equation (5) gives the following condition fo ϕ : cos ϕ = e ( µ sin θ ). (7) Pefoming standad spatial otations its easy to obtain the following fomulae between the angles ϕ, ϕ, θ, θ, φ and φ : whee cos(ϕ ϕ ) = cos θ cos θ, cos(φ φ ) = tanθ tan θ. (8) Using equation (8) to ewite equation (5) yields = µ sin θ ecos ξ, (9) «ξ = cos cos θ + ϕ. () cos θ We now use the fact that in spheical coodinates v φ = sin θdφ/dt = h sin θ / sin θ, v θ = dθ/dt = (dθ/dφ)(v φ / sin θ), v = d/dt = (d/dθ)(v θ /). With this and using equation (8) the expessions fo the velocity field ae given by: µ := h E, ν := v E, (3) whee E := M/ is the gavitational potential enegy of a fluid paticle evaluated at. The quantity µ epesents the atio of centifugal enegy to the gavitational potential enegy of a fluid paticle located at the equato of the cloud s bode (i.e. =, θ = π/). On the othe hand, ν epesents the atio of adial kinetic enegy to the gavitational potential enegy evaluated at. With these paametes, equation () can be witten into dimensionless fom as v θ = µsin θ sin θ v φ = µsin θ sin θ, () `cos θ cos θ /, () eµsin ξ sin θ v = ( e cos ξ). (3) In these equations, the velocity has been measued in units of the Kepleian velocity v k at position multiplied by, i.e. vk = E. ε = ν + µ sin θ, (4) whee the dimensionless enegy ε is given by ε := E/E. The tajectoy of each fluid paticle is contained on a plane P and is given by a conic section. The oigin epesents one of the foci of the obit. At a specific initial position, 3 PATICLE NUMBE DENSITY In ode to obtain the paticle numbe density field we stat fom the continuity equation fo a stationay situation, i.e.

3 Analytic solutions to the accetion of a otating finite cloud towads a cental object I. Newtonian appoach. 3 4 A SPECIAL CASE: ν. (nv) =. (4) We now integate the pevious equation ove a volume consisting on a collection of steamlines. In othe wods, this tube is built in such a way that it s lateal suface is bound by steamlines and its lids by spheical sections, the uppe one at = and the lowe at any abitay distance such that < <. With this selection, equation (4) can be integated to obtain: nv da = = nv da. (5) Taking into account the fact that the diffeential aea element is given by da = sin θ dθ dφ e + sin θ d dφ e θ + d dθ e φ, (6) then, equation(5) tansfoms to n ν sin θ dθ dφ = n v sin θ dθ dφ. (7) Using the symmety of the quantities in the aimuthal angle involved in the pevious equation, we can integate equation (7) with espect to φ fom to π, to obtain n = Ṅ» «sin θ θ v 4π. (8) sin θ θ Fom equation (9), it follows that esin ξ sin θ sin θ «θ θ «cos θ cos θ + = + µ 3 cos θ µν cos θ sin θ + cos θ cos θ cos θ cos θ «. (9) Substitution of equations (9) and (3) into (8) leads to the equied paticle numbe density field: Nµsin n = j» θ cos θ 3 cos + µ θ π cos θ cos θ«`cos θ cos θ / + + µν cot θ ` + cos θ cos θ ff. () As a esult of the collision at the equatoial plane between the steamlines coming fom the nothen hemisphee with those coming fom the southen one, a shock font is expected to fom, as well as the fomation of an equatoial disc shape stuctue (see e.g. Ulich 976; Mendoa et al. 4, and efeences theein). Taking θ π/ in equation (9) followed by the substitution θ π/ we obtain the following expession fo the accetion disc adius d = µ + µν () In this section we conside the case when ν. If we assume the accetion ate N in equation () to be constant even at this limit, then necessaily the paticle numbe density at the cloud s bode diveges, i.e. n. Also, using equation (3) in this limit, it follows fom equation () that the angle ϕ and d µ fom equation (). In the same limit, equation (9) tansfoms to: = µ sin θ e cos θ cos θ, () Fo the velocity field, it follows that equations () and () emain unchanged wheeas equation (3) yields eµsin θ `cos θ cos θ / v =. (3) cos θ e cos θ In this paticula case, the paticle numbe density field n is given by: Nµsin n = j» θ cos θ 3 cos + µ θ π cos θ cos θ «`cos θ cos θ / ). 5 CONVEGENCE TO ULICH S 976 MODEL (4) In the oiginal wok of Ulich (976) the gas cloud is taken to be of infinite extent. To ecove Ulich s model fom ous, we fist have to abandon as unit of distance, since in hee. If we conside h to emain finite as, then fom equation () it follows that E v /. Let us take = h /M as a new unit of distance. Unde these cicumstances, we must abandon the paamete ν since it is defined in tems of, and use a moe appopiate one η given by: η := v E (5) whee E = M/. Pefoming those changes equation (9) tansfoms to: = sin θ e cos ξ (6) whee ξ is the same that in equation () and ϕ now satisfies: while e is given by: e = cos ϕ = e Note that equation () becomes (7) q + η sin θ (8)

4 4 S. Mendoa, E. Tejeda d = + η (9).5.5 If we now take the Kepleian velocity v k* evaluated at, i.e. v k* := h M =, (3) * as unit of velocity, then the velocity field equations ()- (3) and the paticle numbe density () take the following fom: esin θ sin ξ v = ( ecos ξ), (3) v θ = sin θ `cos θ cos θ /, (3) sin θ n = v φ = sin θ sin θ, (33) j» Ṅ sin θ cos θ + 4π (3cos θ ) `cos θ cos θ / + + η cot θ ` + cos θ cos θ ff, (34) espectively. Anothe assumption in Ulich s model is that v = and so, η in such a way that Ṅ emains finite. Fom equation (9) it follows that d =, e and ϕ. All the equations obtained within this limit coespond to the ones calculated by Ulich (976). 6 DISCUSSION Figue shows pojected steamlines fo diffeent values of the paametes µ and ν. Each dot on a fixed cuve coespond to the intesection of the eal steamline and the plane φ = const as it is being swept fom φ = to φ = π/. In othe wods, each line is the intesection of a family of eal steamlines all having the same θ, with a fixed plane φ = const. The eason of doing this is because diect pojections of steamlines on a plane φ = const show false intesections of the steamlines on the obtained plots. Figue shows the pojected steamlines fo a fixed µ =.8 and values of ν =.75,.5,, 8 fom left to ight, top to bottom panels. Figue shows also pojected steamlines of the flow fo a fixed value of the adius of the disc d = /. The bottom panels in Figue have tiangles dawn on them. These ae ones out-of-ange in the model due to eal intesections of the accetion steamlines. Note that this behaviou appeas because the model is ballistic. We have exploed numeically in which cases these intesections occu and it tuns out that fo d = /, the steamlines show intesections fo ν.7. Fom values of ν 5 intesections Figue. The figue shows pojected steamlines of the flow fo diffeent values of the paamete ν (see text) at an abitay aimuthal angle φ = const in such a way that the paamete µ =.8. Fom left to ight, top to bottom, these values ae ν =.75,.5,, 8. These pojections do not coespond to the eal steamlines of the flow (see section 6). appea but the out-of-ange egion is kept fixed. Fo the same value of d, plots with ν.5 seem to be about the same. Fom ν.5 show a lot of vaiation as ν inceases. Howeve, as soon as ν eaches a value of 5, the plots ae vey simila as ν inceases without limit. Figues 3 and 4 show density isocontous and density pofiles fo a fixed height measued fom the equatoial plane espectively. The out-of-ange egion appeas only in the bottom panels of Figue 3. Figues 5-7 ae equivalent plots as those of Figues -4 fo the case, which can be thought of as the obvious modifications to Ulich s model, all of them depending on a single dimensionless paamete η defined in equation (5). Note that in all these cases, thee is no intesection of the flow steamlines. This appeas to indicate that as long as the cloud eaches a sufficiently lage adius, with espect to the accetion disc adius d, the intesections disappea. Indeed, we have tested diffeent accetion models vaying d fom. to. and the pevious statement tuns out to be coect. As an example, a value of.8 epesents the bounday of intesections of steamlines fo values of ν =.. We have constucted an analytic accetion flow that expands Ulich s (976) model in the sense that the adius of the igid body otating cloud is finite and an allowance fo a adial accetion velocity at the cloud s bode is assumed. When the adius of the cloud tends to infinity and the adial velocity of the input flow at infinity goes to eo, the model conveges to that of Ulich. We ae developing a elativistic model with all these popeties. These esults will be the subject of a subsequent aticle

5 Analytic solutions to the accetion of a otating finite cloud towads a cental object I. Newtonian appoach Figue. The figue shows pojected steamlines of the flow fo diffeent values of the paametes µ and ν (see text) at an abitay aimuthal angle φ = const in such a way that the adius of the disc has a value d = /, measued in units of the cloud s adius. Fom left to ight, top to bottom, these values ae µ = /,.64,.8, 5.5, ν =,.5,,. These pojections do not coespond to the eal steamlines of the flow (see section 6). The filled tiangles dawn in the bottom plots ae ones out-ofange whee steamlines of the flow intesect. This occus because the flow has been assumed to be ballistic. Figue 3. The figue shows logaithmic paticle numbe density n isocontous. Fom left to ight, top to bottom the values of the paametes of the flow ae the same as in Figue and the values of the isocontous fo each panel ae given by {.3,.4,.6,.75,.5,.75}, {.3,.34,.4,.5,.65,.9,.3,.}, {.,.3,.33,.45,.6,.,.5,.} and {.,.85,.65,.45,.3,.,.,.5} espectively. The dak egions on the bottom panels ae the out-of-ange ones. 7 ACKNOWLEDGEMENTS We thank the pofound comments made by Pof. E. Ley Koo duing the pepaation of this aticle. The authos gatefully acknowledge suppot fom DGAPA (IN93-3) at Univesidad Nacional Autónoma de México (UNAM). SM acknowledges financial suppot ganted by CONACyT (6344). Ulich. K., 976, ApJ,, 377 efeences Beloboodov A. M., Illaionov A. F.,, MNAS, 33, 67 Bondi H., 95, MNAS,, 95+ Bondi H., 5, Accetion. The Scientific Legacy of Fed Hoyle, pp Fank J., King A., aine D. J.,, Accetion Powe in Astophysics: Thid Edition. Accetion Powe in Astophysics, by Juhan Fank and Andew King and Deek aine, pp ISBN Cambidge, UK: Cambidge Univesity Pess, Febuay. Hueta E. A., Mendoa S., 7, evista Mexicana de Astonomia y Astofisica, 43, 9 Landau L., Lifshit E., 989, Mechanics, 3d ed. edn. Vol. of Couse of Theoetical Physics, Pegamon Lee W. H., amie-ui E., 6, ApJ, 64, 96 Mendoa S., Cantó J., aga A. C., 4, evista Mexicana de Astonomia y Astofisica, 4, 47 Michel F. C., 97, Astophysics and Space Science, 5, 53

6 6 S. Mendoa, E. Tejeda Figue 4. Both panels show density pofiles fo a fixed aimuthal angle, i.e. φ = const and a fixed height =,.,.,.3,.5 measued fom the equatoial plane. The figues wee dawn using values of ν =.5,. espectively fo a fixed value of the adius of the disc d = /. Note that fo a height = the paticle numbe density diveges at =, ± Figue 5. The plots shows pojections of steamlines (in the same sense as Figue which coespond to an infinite value of the cloud s adius. Fom top to bottom the paamete η has values of,, espectively (see Section 5). The top panel coesponds to Ulich s accetion model. Distances ae measued in units of the adius of the disc d.

7 Analytic solutions to the accetion of a otating finite cloud towads a cental object I. Newtonian appoach Figue 6. The panels show the logaithmic paticle numbe density isocontous measued in units of Ṅ/4π. Distances ae measued in units of the adius of the disc d and fom top to bottom the sets of contous take the following values: {.55,.47,.35,.,.,.5,.}, {.5,.4,.6,.8,.,.5, }, and {.7,.8,.95,.,.5,.5, 3, 3.5}. Fom Figues 5 and 6 it is clea that the paticle numbe density does not accumulate that much in the oigin as η inceases Figue 7. Fom top to bottom, the figues show the logaithm of the paticle numbe density of the flow in units of N/4π, coesponding to the same values of η that Figue has. All plots have density as a function of the equatoial adial coodinate = sinθ fo diffeent fixed heights =,.5,.3,.5. In all figues, the plot at = has density divegences at =, ±.