numerical simulations of the bondi  hoyle
lyttleton ( hereafter bhl ) accretion flow , in 2d and 3d , have not only enabled the determination of the mean rate of accretion of mass , linear and angular momentum , but have also significantly modified our understanding of its dynamics .
they have revealed that the accretion column foreseen by hoyle & lyttleton ( 1939 ) widens to form a detached bow shock if the adiabatic index of the accreted gas is @xmath0 ( hunt  1971 ) or @xmath3 ( hunt  1979 ) .
the shock structure observed in simulations was used by eadie ( 1975 ) to interpret the x  ray observations of vel x-1 .
+ real situations like wind accretion in x - ray binaries require taking into account the density and velocity gradients in the incoming flow ( see ishii  1993 and references therein ) , the effects of rotation ( see theuns  1996 and references therein ) as well as the intricacies of radiative transfer ( see an overview in blondin  1994 ) . for the sake of a better understanding , we restrict our present investigation to the axisymmetric meridional accretion of a gas with uniform adiabatic index onto a gravitating accretor with constant velocity .
+ although this configuration is highly simplified , a completely unexpected non  axisymmetric instability was discovered in numerical simulations for supersonic flows ( fryxell & taam  1988 ) .
this instability could explain the observed variability of the radio or x - ray luminosity in some astrophysical systems ( ruffert & melia  1994 ) .
+ however , its is not easy to establish the consistency of the results of the few published 3d numerical simulations which depend on the details of the numerical method : no instability appeared in the 3d sph simulations by boffin ( 1991 ) whereas it was observed in the 3d eulerian simulations by matsuda ( 1992 and references therein ) and ruffert ( 1995 and references therein ) . dome
shaped structures frequently appeared with roe s method ( approximate riemann solver ) of matsuda ( 1991 ) , but not in the nested grid piecewise parabolic method ( ppm ) of ruffert ( 1991 ) .
since numerical viscosity appears to play an important role , a sufficient resolution of the simulation is crucial . in this respect ,
numerical simulations encounter the following difficulty : the smallness of the accretor in simulations is limited by the computational power , because the shortest timestep of the simulation is usually determined by the region of high velocities in the vicinity of the accretor .
the ratio of scales between the accretion radius and the radius @xmath4 of a compact star , moving with a supersonic velocity @xmath5 is typically of order : @xmath6 where @xmath7 is the schwarzschild radius .
note that in the case of accretion onto a neutron star , the relevant accretor boundary is the boundary of its magnetosphere , typically a factor @xmath8 larger than the schwarzschild radius .
some of the most recent 3d simulations ( ruffert  1992 , 1994a ) could allow a ratio up to @xmath9 ( the full domain of the simulation extends up to @xmath10 accretion radii ) , by using a ppm code with multiple nested grids .
such an accreting sphere is still more than 100 times larger than @xmath7 , and typically @xmath11 times larger than the magnetosphere of a neutron star .
one would like to be able to predict , on an analytical basis , which changes are likely to occur if the size of the accretor is further decreased .
the role of the boundary conditions at the surface of the accretor is closely related to the position of the sonic surface , and we shall use analytical arguments to obtain some insights on this question .
an analytical description of axisymmetric stationary flows is also a first step towards the analysis of their stability , which we postpone to a future paper .
bisnovatyi - kogan ( 1979 ) found self - similar stationary solutions , and one would like to know to what extent any generic stationary solution would behave in the same way as their solution .
more generality is gained by studying , in the same spirit as theuns & david ( 1992 ) for the spherical case , the first order behaviour of all quantities near the accretor for an axisymmetric flow , and in particular their departure from sphericity .
some analytic formulae have been proposed to fit the mass accretion rate observed in simulations ( hunt  1979 ; ruffert  1994b ) , but none of them filled the gap between the well studied cases of axisymmetric accretion of dust ( hoyle & lyttleton  1939 ; bondi & hoyle  1944 ) on the one hand , and spherical accretion of a gas ( bondi  1952 ) on the other , for an arbitrary value of the adiabatic index .
continuing in this direction , we add to these reference models the spherical accretion of a shocked gas with uniform adiabatic index , in order to see the effect of both the kinetic and the thermal energies , and extend this approach to axisymmetric meridional stationary flows .
+ we first recall in sect .
[ sgeneral ] the general equations governing a stationary flow and formulate the axisymmetric problem with a single equation in cylindrical coordinates .
the shape of the sonic surface is analysed in sect .
[ ssonic ] .
we establish some properties of the flow by performing local expansions in the vicinity of the accretor in sect .
[ sexpansions ] .
we present in sect .
[ sspherical ] the analytical solution of the spherical accretion from a gas with constant kinetic and thermal energies at infinity .
an interpolation formula for the mass accretion rate of axisymmetric meridional flows is proposed in sect .
[ sinterp ] .
we summarize our conclusions in sect .
[ sconclusion ] .
the continuity equation for a gas of velocity @xmath12 and density @xmath13 is : @xmath14 we restrict our study to the case of a gas of uniform adiabatic index @xmath15 .
the entropy @xmath16 of a gas with pressure @xmath17 is defined as : @xmath18 in which @xmath19 is the gas constant and @xmath20 is the mean molecular weight .
the entropy is constant along each flow line ( @xmath21 ) , except when matter passes through a shock .
+ the temperature @xmath22 and the square of the sound velocity @xmath23 are simply related through a multiplicative constant : @xmath24 for the sake of simplicity , we shall set the ratio @xmath25 everywhere in what follows . since the mean molecular weight is constant throughout our study
, this convention is equivalent to changing the units of temperature and entropy , without loosing any generality .
+ using eq .
( [ defentro ] ) , the pressure force can be written as : @xmath26 the euler equation for a stationary flow with velocity @xmath12 around an accretor of mass @xmath27 is then : @xmath28 in which @xmath29 is the distance to the centre of the accretor and @xmath30 is the gravitational constant .
+ the bernoulli equation is obtained by integrating the euler equation along a flow line .
it states that the quantity @xmath31 @xmath32 is constant along each flow line .
the three terms of this sum are respectively the kinetic , the thermal , and the gravitational energies .
+ the vorticity of the flow is simply defined as @xmath33 by using the definition of the bernoulli constant ( [ bernoulli ] ) in the euler equation ( [ euler ] ) , we obtain the following equation for any stationary flow , @xmath34 this formula is particularly useful because it expresses the vorticity in terms of invariant quantities @xmath31 and @xmath16 . + the bernoulli constant @xmath31 remains constant along a flow line , even through a shock .
by contrast , the entropy @xmath16 increases through a shock . + in the classic bhl case of a flow with constant and uniform velocity at infinity , the bernoulli constant is uniformly constant ( @xmath35 ) , and so does not contribute to the vorticity : @xmath36 if , in addition , the flow is isentropic at infinity , the only possible source of vorticity gradients is the shock front .
note that eq .
( [ vortentro ] ) implies the conservation of entropy ( @xmath21 ) .
thus the only remaining differential equations are the continuity eq .
( [ continuity ] ) and the vorticity eq .
( [ vortentro ] ) .
we consider the axisymmetric stationary flow of a gas with uniform density and constant velocity @xmath37 with respect to the accretor .
the geometry of the flow lines is described by a single function , the angle @xmath38 between the radial direction and the flow line , defined by : @xmath39 we define by @xmath40 the direction perpendicular to the velocity , in the plane containing the symmetry axis . denoting by @xmath41 the unit vectors in cylindrical coordinates , we can write : @xmath42 when useful throughout the paper , we denote by @xmath43 the dimensionless distance normalized to the modified accretion radius @xmath44 , defined as : @xmath45 it coincides with the classic accretion radius @xmath46 for large mach numbers , when the kinetic energy largely exceeds the thermal energy : @xmath47 we shall denote distances normalized to @xmath44 by a bar .
+ using eqs .
( [ defentro ] ) and ( [ bernoulli ] ) , all relevant flow quantities can be rewritten as functions of @xmath48 : @xmath49 the azimuthal component of the euler eq .
( [ euler ] ) , multiplied by @xmath29 , is : @xmath50 we denote by @xmath51 the third and only non - zero component of the vorticity vector : @xmath52 by combining the vectorial vorticity eq .
( [ vortentro ] ) and the scalar continuity eq .
( [ continuity ] ) , we obtain a third order differential system in terms of @xmath53 only , with variables @xmath54 ( see appendix  [ apartial3 ] ) .
this system is hyperbolic in regions where @xmath55 , and elliptic where @xmath56 .
it requires _ three _ functions @xmath57 as boundary conditions on the surface of the accretor .
it should be noted that if the entropy is uniform , the flow is potential ( @xmath58 ) , and the bernoulli and the continuity equations provide a second order partial differential equation on @xmath59 . in this case , the boundary conditions at the surface of the accretor consist of _ two _ functions of @xmath60 , namely @xmath61 and its radial derivative .
we can simplify the system by using the modified stream function introduced by crocco ( 1936 ) , defined as follows : @xmath62 where there is no shock , @xmath63 is conserved along flow lines ( because @xmath64 ) .
it is also conserved across shocks , since the derivatives of @xmath63 have only a finite discontinuity there .
the entropy can be expressed as a function of @xmath63 , and thus becomes a function of a single variable @xmath65 , defined once and for all along the shock .
+ let us index the flow lines with their distance @xmath66 to the symmetry axis ( their impact parameter ) at infinity .
only the gas within a cylinder of radius @xmath67 is accreted , so that the total mass accretion rate is @xmath68 . at infinity , the mass flux inside the cylinder of radius @xmath66 is directly proportional to @xmath69 through the following product : @xmath70 for the sake of clarity of the equations , let us define : @xmath71 because the total mass accretion rate is finite , we deduce from eq .
( [ massflux ] ) that @xmath63 is bounded in the vicinity of the accretor , and @xmath72 . + using the definition of @xmath63 and the bernoulli equation , we can formally write the mach number @xmath73 as a function of the gradient of @xmath63 through the following algebraic equation : @xmath74 at a given position ( @xmath54 )
, the maximum of the right hand side term is reached at @xmath75 .
this defines a maximum value for the gradient of the stream function , for the local mass flux according to eq .
( [ massflux ] ) . for any value of @xmath76 below this maximum , two mach numbers are solutions of the eq .
( [ implicit ] ) . as for spherically symmetric flows classified by bondi ( 1952 ) ,
one solution is subsonic , while the other is supersonic .
+ equation  ( [ implicit ] ) can be derived with respect to @xmath29 or @xmath60 to express the partial derivatives of @xmath73 with respect to @xmath29 or @xmath60 as function of @xmath73 and the derivatives of @xmath63 .
the vorticity eq .
( [ vortentro ] ) can therefore be expressed with derivatives of @xmath63 only , as follows : @xmath77 this equation can be viewed as an explicit partial differential equation of second order on @xmath63 , with @xmath73 implicitly depending on @xmath76 through the algebraic eq .
( [ implicit ] ) . in an equivalent way ,
( [ diffpsi ] ) can also be considered as a second order polynom in @xmath78 , whose two roots explicitly determine @xmath78 as a function of @xmath63 and its derivatives . inserting one of the two roots into eq .
( [ implicit ] ) , we obtain an implicit differential equation of second order containing explicitely only @xmath63 and its derivatives of first and second order .
+ this compact formulation will allow us to perform in sect .
[ sexpansions ] a direct analysis of the behaviour of the solution in the vicinity of @xmath79 .
88.9 mm by 70.9 mm ( topology )    numerical simulations of axisymmetric accretion with different mach numbers , adiabatic indices and accretor sizes have exhibited three different topologies of the sonic surface .
these three topologies are illustrated in fig .
[ topologies ] .    - type f ( fast accretion ) : the sonic surface is detached from the accretor .
this occurs for subsonic flows for small enough accretors ( ruffert  1994b , 1995 , 1996 )
.    - type sf ( mixed slow / fast accretion ) : the sonic surface reaches the surface of the accretor .
it subtends an angle @xmath80 relative to the axis of symmetry at the position of the accretor .
this is the case for subsonic flows with large accretors , and supersonic flows with @xmath0 ( matsuda  1992 ; ruffert & arnett  1994 ; ruffert  1994a ) .
- type fsf ( mixed fast / slow / fast accretion ) : both the sonic surface and the shock surface touch the surface of the accretor .
they respectively subtend the angles @xmath80 and @xmath81 relative to the axis of symmetry .
this occurs for all supersonic flows when the accretor is large enough , for all the known supersonic 3d simulations with @xmath82 close to one ( ishii  1993 ; ruffert  1996 ) .
more complicated topologies were observed by hunt ( 1979 ) and ruffert ( 1995 ) in simulations of supersonic flows with @xmath3 , where the bow shock is linked to the accretor by a conical shock . these can also be classified as type fsf , with the difference that the entropy in the supersonic flow ahead of the attached shock is not uniform .
+ the topology of the flow in numerical simulations appears to depend strongly on the size of the accretor : all the supersonic flows are of type fsf if the accretor is big enough , but some become type sf or f when the accretor size is decreased . in order to spare computing time , which increases with decreasing accretor size
, one would like to use the largest accretor size leading to a physically relevant accretion process ( cf .
eq .  [ ratio5 ] ) .
type f flows are numerically easier to handle , because unphysical effects due to a possible reflection of waves at the boundary of the too large accretor can be avoided .
let us consider a flow tube of cross section @xmath83 , parametrized along the tube by the curvilinear variable @xmath84 .
we denote by @xmath85 the entropy that this flow tube gains passing through a shock , and by @xmath86 and @xmath87 the density and the sound speed at large distances from the accretor , before the shock . the density @xmath13 and sound speed @xmath88 after the shock satisfy the following equation : @xmath89 using the bernoulli equation ( [ bernoulli ] ) and the entropy equation ( [ entrop ] ) , we can write the momentum as a function of the mach number @xmath73 and the position @xmath29 : @xmath90    86.7 mm by 72.3 mm ( flowtube )    the continuity equation ( [ continuity ] ) states that the mass flux is conserved along the flow tube : @xmath91 by substituting eq .
( [ momentum ] ) into eq .
( [ dmoment ] ) , we obtain , along a given flow line : @xmath92 where the derivatives are calculated along a flow line , @xmath93
. the regularity of @xmath94 at the sonic point @xmath95 implies : @xmath96 we have denoted by @xmath97 the radius of the sonic sphere in the spherical case : @xmath98 it is remarkable that the spherical sonic radius is a function of the total energy @xmath99 only , and is independent of the entropy of the flow .
+ the regularity condition ( [ sonic ] ) is a constraint on the shape of the flux tube at the sonic radius .
note that this condition is valid for any stationary flow , with no hypothesis on its symmetry .
let us illustrate it for an inflowing axisymmetric flow tube : the spherical case corresponds to radial flow lines , and thus to a conical flow tube .
if the sonic radius is smaller than the spherical one , the cross section of the flow tube increases faster than the square of the distance to the accretor , thus the azimuthal component @xmath100 of the velocity at the sonic radius must be oriented towards the axis of symmetry .
conversely , if the sonic radius is larger than the spherical one , the cross section of the flow tube increases more slowly than the square of the distance to the accretor , and @xmath100 must be oriented away from the axis of symmetry ( see fig .  [ convexe ] ) .
86.7 mm by 72.3 mm ( property )    we prove in appendix  [ aproof ] the following property @xmath101 : + _ if the sonic surface is detached from the accretor , it must intersect at least once the sonic sphere of radius @xmath97 of the corresponding spherically symmetry flow with same energy_. + @xmath97 is defined in sect .
[ ssr0 ] by eq .
( [ rspherique ] ) .
this property can be understood intuitively from fig .
[ propriete ] , by using the geometrical interpretation of eq .
( [ sonic ] ) in sect .
[ ssr0 ] , in the simple case where the azimuthal velocity is of constant sign along the sonic surface in the interval @xmath102 .
if the azimuthal velocity is directed away from the axis on one side of the accretor ( for example , @xmath103 ) , it will be oriented towards the axis on the other side ( @xmath104 ) , thus implying by continuity an intersection of the sonic surface with the sonic sphere .
we show in appendix  [ aproof ] that this property is true for any axisymmetric flow , without any hypothesis on the distribution of azimuthal velocities .
+ former studies of the bhl accretion have dealt with two set of length scales , one associated to the hoyle - lyttleton accretion of dust ( the accretion radius @xmath46 ) , and the other associated to the spherical accretion ( the bondi radius @xmath105 and the corresponding sonic radius @xmath106 ) .
none of these scales took into account both the kinetic energy and the thermal energy of the gas .
in particular , the significance of the sonic radius @xmath106 for axisymmetric flows has always been dubious .
the property @xmath101 establishes the fundamental meaning of the spherical radius @xmath107 , in relation with the sonic surface , for axisymmetric flows of type f. + the geometrical property @xmath101 is useful for numerical simulations of the bhl flow , in order to fix the maximum size of the totally absorbing accretor . as discussed in sect .
[ sstypes ] , since the sonic surface must cross the sonic sphere , the accretor size @xmath4 should at least be smaller than the spherical sonic radius @xmath97 . however , although the sonic sphere is larger than the accretor size used in numerical simulations with @xmath3 and @xmath108 ( see ruffert ... ) , the sonic surface and the shock surface appear to be attached to the accretor .
although necessary , the condition @xmath109 is far from being sufficient .
+      since @xmath110 , the sonic surface of flows with @xmath0 must always reach the surface of the accretor .
following the nomenclature introduced in sect .
[ sstypes ] , _ the topology of the accretion with @xmath0 is either of type sf or fsf , for any accretor size and any mach number at infinity_.      the hoyle - lyttleton limit for the mass accretion rate at high mach numbers requires that the accreted matter comes from a cylinder with a radius scaling like the accretion radius ( hoyle & lyttleton  1939 ) . if @xmath82 is close to unity , the sonic spherical radius @xmath97 is @xmath111 larger than the accretion radius @xmath46 : @xmath112 if such a flow were of type f , the geometrical property @xmath101 would imply that the sonic surface goes out as far as @xmath111 times the accretion radius .
this leaves only two possibilities for a stationary flow :    \(i ) the accretor is able to keep material bound up to distances comparable to @xmath113 , and the mass accretion rate is infinitely larger than the hoyle - lyttleton limit .
this is the case of the spherical flows we shall study in sect .
[ sspherical ] .
\(ii ) the accretor can not retain material at such large distances , and the flow is of type sf or fsf , as observed in numerical simulations with @xmath82 close to one by ishii ( 1993 ) and ruffert ( 1996 ) .
+ let us note that the artificially constructed type f geometry , where the distance of the stagnation point exceeds @xmath97 , while the accreted matter collects within a cylinder of radius @xmath114 , would be strongly unstable to the kelvin - helmholtz instability : the growth rate would scale like @xmath115 while the time needed to reach the accretor would scale like @xmath116 .
+    86.7 mm by 72.3 mm ( 9/7 )    consequently , we rule out the possibility that the sonic surface is detached in the case of an isothermal stationary flow . according to eq .
( [ discrep ] ) , the spherical sonic radius is larger than the accretion radius at high mach numbers if ( see also fig .  [ r0ra ] )
spherically symmetric accretion flows with @xmath15 are separated into two distinct classes , depending on whether the flow is subsonic or supersonic at the accretor boundary ( bondi  1952 ) .
we learnt from eq .
( [ implicit ] ) in sect .
[ sscrocco ] that these two branches of solutions also exist for axisymmetric flows , for a given value of the local mass flux at the surface of the accretor .
the azimuthal parameter @xmath60 , however , introduces an additional degree of freedom , such that the solution might be subsonic in an angular sector , and supersonic in an other ( for a sonic surface attached to the accretor ) .
we must also distinguish between angular directions towards which the flow lines converge ( direction of  accreting " flow lines ) , and directions along which no gas is accreted ( direction of  passing " flow lines ) . by conducting a local analysis , we aim to clarify the various possible configurations .      for each flow line ( indexed by @xmath66 ) reaching the accretor ,
we denote by @xmath118 the angle between the velocity vector @xmath12 and the axis of symmetry , when @xmath119 .
our notations are illustrated in fig .
we define the direction @xmath120 along the axis of symmetry on the downstream side of the accretor , so that the function @xmath118 decreases monotonically in the interval @xmath121 $ ] , with @xmath122 and @xmath123 .
we define @xmath124 as the largest azimuthal angle , at @xmath125 , among all the converging flow lines , apart from the axis of symmetry : @xmath1260,\varpi_0[)\}\;.\ ] ] this angle can be smaller than @xmath127 .
in particular , @xmath128 for flows of type fsf . without loosing any significant physical generality
, we can assume that @xmath1290,\theta_{\rm c}[\subset \{\theta_0(]0,\varpi_0[)\}$ ] in the bhl flow .
+    86.7 mm by 35.2 mm ( nota )    near the surface of the accretor , let us denote by @xmath130 the local flux of mass between the azimuthal angles @xmath60 and @xmath131 : @xmath132 since the stream function is bounded in the flow ( eq .  [ massflux ] ) , integrating eq .
( [ tanbeta ] ) with respect to @xmath29 implies : @xmath133\;. \label{vrvt}\ ] ] thus , if the mass accretion rate @xmath134 is locally not zero , the flow lines become locally radial when @xmath119 ( no  spiraling " of the flow lines ) . on the other hand ,
if the flow lines are not radial , eq .  ( [ mpsi ] ) implies that the mass accretion rate is zero in this direction ( passing flow lines ) .
+ note that along a direction of accretion , we can consider a flux tube of cross section @xmath135 around the flow line such that @xmath136 , and write eq .
( [ locflux ] ) as follows : @xmath137 the quantity @xmath134 is consequently finite only if the flow tubes are asymptotically conical ( @xmath138 ) .
this leads us to call _ regular _ the accretion onto a point  like accretor such that @xmath139 remains finite , and _ singular _ when it is locally infinite .
+ we can also express the local mass accretion rate in terms of the mass flux between the cylinders of radii @xmath66 and @xmath140 at infinity : @xmath141 thus singular accretion corresponds to a localized azimuthal angle @xmath142 where @xmath143 .
an example of singular accretion is the self - similar solution found by bisnovatyi - kogan ( 1979 ) , where the whole accretion is along a single azimuthal angle @xmath1440,\varpi_0])$ ] .
+ for a regular flow ( @xmath145 ) , we take the limit of eq .
( [ consflux ] ) when @xmath119 and obtain , along a given flow line , @xmath146 this equation shows that for @xmath1 , the mach number along regular flow lines either tends to zero or to infinity when approaching the accretor . as in the case of spherical accretion ( bondi  1952 ) no intermediate finite mach number is possible .
moreover , this equation proves , by continuity , that a mixed subsonic / supersonic accretion requires a singular accretion at the transition if @xmath1 .
the sonic surface may still reach the accretor in a regular flow , by separating a region of supersonic accretion from a subsonic region of passing flow lines , as we shall see in sect .
[ ssinterface ] .
+ we will now examine the different regions of passing flow lines , supersonic and subsonic regular accretion successively , for @xmath1 , and treat the case @xmath0 separately .
we shall use the generic name @xmath147 for a function which tends to zero with @xmath29 , and the subscript @xmath148 for the limit of functions when @xmath125 .
since the flux of mass is zero , @xmath151 can not vanish in the interval @xmath129\theta_{\rm c},\pi[$ ] . according to our sign convention in sect .
4.3.3 and the definition of @xmath152 in eq .
( [ defbeta ] ) , @xmath153\theta_{\rm c},\pi[\;.\ ] ] since @xmath154 , the function @xmath155 must have a maximum in the interval @xmath129\theta_{\rm c},\pi[$ ] .
+    \(i ) if the shock is detached , the entropy is finite in the flow . according to eq .
( [ mm ] ) , the mach number must either diverge or tend to zero in order to satisfy @xmath156 .
+ let us suppose that @xmath157 .
we use eq .
( [ vm ] ) and obtain : @xmath158 using the vorticity equation ( [ vortentro ] ) , the entropy gradient can be expressed as : @xmath159 integrating this equation over a segment between @xmath160 and @xmath161 , for @xmath119 , contradicts the fact that the entropy is bounded for a flow with a detached shock .
consequently the mach number tends to 0 when @xmath125 , and @xmath162 .
_ we conclude that all the flows with a detached shock and @xmath163 are of type sf . _
+    \(ii ) if the shock is attached , the incoming supersonic velocity diverges near the accretor like @xmath164 when @xmath119 .
we decompose it into a parallel and a perpendicular component with respect to the shock .
the mach number @xmath165 associated to the component of the velocity , after and perpendicular to the shock , is subsonic according to the rankine  hugoniot jump conditions ( see landau & lifschitz  1987 ) : @xmath166 the mach number @xmath167 associated to the component parallel to the shock can be supersonic , but not arbitrarily large : @xmath168 _ the total mach number immediately after the attached shock is consequently finite when @xmath119_. this result will be important in order to study the local stability of such flows to the rayleigh  taylor instability .      if @xmath157 and @xmath134 is finite , eqs .
( [ vm])([pm ] ) can be expanded as follows : @xmath169 the first order azimuthal velocity can be estimated using eq .
( [ vrvt ] ) in the azimuthal euler eq .
( [ eulert ] ) .
if the pressure is not spherically symmetric to first order , @xmath170 @xmath171 being negative , the azimuthal velocity is oriented towards the region of increasing pressure .
the azimuthal pressure force is balanced by this dynamical term .
we deduce from eqs .
( [ defbeta ] ) and ( [ vtheta ] ) an upper bound for the departure from spherical symmetry of the function @xmath152 : @xmath172 by studying the leading terms of eq .
( [ diffpsi ] ) when @xmath119 , we derive in appendix  [ ainf53 ] a lower bound for the departure from spherical symmetry of the function @xmath152 , if the entropy of the flow is not uniform : @xmath173 on the other hand , the function @xmath152 may approach 0 much more rapidly than @xmath174 if the entropy is uniform . to phrase it more qualitatively
, the presence of an entropy gradient sets a lower bound on the  bending " of the flow lines near the accretor .
if a parallel can be drawn with the simpler case of two - dimensional supersonic flow past a concave profile ( landau & lifshitz 1987 ) , this  bending " could be related to the formation of long  lived radial shocks in the supersonic region of flows with a non  uniform entropy , as observed in numerical simulations with @xmath3 in 2d ( hunt 1979 ) and 3d ( ruffert  1995 ) .
+ because the radial velocity ( [ vrsup ] ) is spherically symmetric to first order , the azimuthal dependence of the accreted momentum @xmath175 comes mainly from the azimuthal dependence of the density @xmath13 , which we decompose into the contributions of pressure and entropy , @xmath176 thus we expect that for a flow with uniform entropy ( no shock ) , the mass is mainly accreted from the hemisphere where the pressure is the highest .
the geometrical interpretation of eq .
( [ sonic ] ) implies that near the axis of symmetry ( for @xmath177 and @xmath178 ) , the azimuthal speed is oriented away from the axis where @xmath179 , and towards the axis where @xmath180 . if one assumes that the sign of the azimuthal velocity does not change in the supersonic region between the sonic surface and the accretor , eq .
( [ prevt ] ) indicates that the highest pressure is on the side of the accretor where @xmath179 . + for an isentropic subsonic flow , the natural orientation of the azimuthal velocity is set by the sign of the angular momentum of the fluid element at infinity .
a subsonic medium would flow from the _ left _ of the accretor in fig .
[ propriete ] .
therefore the sonic surface is displaced opposite to the direction of the flow at infinity , and the maximum flux of mass is expected on the back hemisphere of the accretor . indeed , in the subsonic numerical simulations of a flow with @xmath3 by hunt ( 1979 ) in 2d , and by ruffert ( 1995 ) in 3d , the sign of the accreted linear momentum is negative for the smallest accretor considered .
+ if the flow is supersonic at infinity , the position of the shock along the axis is displaced in the direction of the flow , compared to its position for a spherical flow with same kinetic and thermal energies at infinity , because there is now a leakage of gas away from the accretor and the pressure ahead of the accretor is lower than in the corresponding spherical case .
the sonic point along the axis is displaced accordingly .
the sonic surface is consequently shifted in the direction of the flow , leading to a highest pressure on the front hemisphere of the accretor .
the simplest configuration is the one in fig .
[ propriete ] , where the supersonic flow now comes from the _ right _ of the accretor . + nevertheless , the entropy generated by the shock decreases with increasing distance from the axis of symmetry ( see appendix  [ aaxis ] ) , so that entropy and pressure contribute in opposite directions to the density in eq .
( [ densentr ] ) , precluding a direct conclusion about its azimuthal distribution near the accretor . however , the accreted linear momentum is negative in the numerical simulations by ruffert ( 1996 and references therein ) indicating that the entropy effect is stronger than the pressure effect .      if @xmath181 and @xmath134 is finite , eqs .
( [ vm])([pm ] ) can be written as follows : @xmath182 we use eqs .
( [ vrsub ] ) and  ( [ vrvt ] ) to estimate the dynamical terms balancing the azimuthal pressure force in eq .
( [ eulert ] ) and obtain : @xmath183 comparing this to the first order of @xmath184 given by eq .
( [ csub ] ) , we conclude for @xmath1 that _ the pressure in the regular subsonic region is spherically symmetric to first order when @xmath125_. + thus eq .
( [ psub ] ) implies that the subsonic flow for @xmath1 , with its first order expansions ( [ vrsub])([machsub ] ) , is possible _ only if the entropy is uniform_. + in particular , the spherical solution proposed by theuns & david ( 1992 ) , connecting a supersonic solution to a subsonic solution through a shock , can not resist the slightest asphericity if the accretor is point  like . the azimuthal pressure force due to
the convergence of flow lines of different entropies is too big to be balanced by dynamical forces in the regular subsonic region close to the accretor .
nevertheless , this constraint of the azimuthal euler equation is less stringent for large accretor sizes .
one can also argue that the radiative processes near the surface of the star may homogenize the entropy , and allow for a spherically symmetric subsonic accretion .
+ the entropy in flows of types sf and fsf is a _
strictly _ decreasing quantity along the shock :    \(i ) in type fsf flows , the velocity diverges to infinity when @xmath185 , thus the entropy created along the shock decreases strictly with distance from the accretor for type fsf flows , precluding a subsonic accretion with a uniform entropy .
\(ii ) for type sf flows , we show in appendix  [ aaxis ] that the entropy created along the shock strictly decreases with distance from the axis of symmetry .
+ we conclude that there can not be a region of regular subsonic accretion in shocked flows with @xmath1 , for a point  like accretor .
the subsonic region may reach the surface of the accretor for any finite size , but _
the mass accreted in the shocked regular subsonic region @xmath186 tends to zero when the accretor size decreases to zero_.      since no regular subsonic accretion is possible with a non  uniform entropy , all the regular flow lines are ultimately supersonic when they reach the point
like accretor @xmath187 .
note that in the singular solutions by bisnovatyi - kogan ( 1979 ) , the mach number tends to infinity along each flow line . + together with the results of sect .
[ noaccretion ] , we conclude that for a flow with a detached shock , and @xmath1 , @xmath188 according to eq .
( [ sonic ] ) , the cross section of the flow tube reaching the sound speed at the surface of the accretor varies like : @xmath189 comparing this equation to eq .
( [ regular ] ) , we conclude that the azimuthal angle @xmath80 can not be a direction of regular accretion . however , @xmath80 is not necessarily a direction of singular accretion .
a regular flow simply requires that @xmath1900,\varpi_0[)\}$ ] .
we denote by @xmath191 the limit of the mach number along a converging flow line .
for @xmath1490,\theta_{\rm c}[$ ] , eq .
( [ mm ] ) implies that the local mass accretion rate @xmath134 is always finite , its maximum value being reached for @xmath192 : @xmath193 the pressure forces for @xmath0 are strong enough to impede an arbitrarily high mass flux , even locally .
thus accretion flows with @xmath0 are always regular , and the maximum mass accretion rate is the spherical bondi value .
this result is particularly important to understand the instability observed in numerical simulations .
since the most unstable simulations were observed with @xmath0 , we can restrict our future stability analysis to the simpler case of regular flows .
+ using eq .
( [ vm ] ) and  ( [ vrvt ] ) in the azimuthal euler eq .
( [ eulert ] ) , we find that the azimuthal pressure force is bounded by eq .
( [ maxip ] ) . comparing it to eqs .
( [ rhom ] ) and  ( [ pm ] ) , we conclude that the pressure is spherically symmetric to first order for @xmath1490,\theta_{\rm c}[$ ] . according to eq .
( [ densentr ] ) , the density is spherically symmetric to first order if there is no shock ( uniform entropy ) , but increases towards the back hemisphere if there is a shock . + with eqs .
( [ vm ] ) to  ( [ pm ] ) , we conclude that :    \(i ) in the absence of a shock , the pressure , density , temperature , velocity and mach number are spherically symmetric to first order when @xmath194 .
the accreted momentum is therefore spherically symmetric as well .
this is confirmed by the numerical simulations by ruffert ( 1994b ) , with @xmath195 and @xmath0 .
the accreted linear momentum decreases to zero with the size of the accretor ( see fig .  2 , model
sl " ,  sm " ,  ss " ) .
moreover , the uniform mach number @xmath196 must be lower than or equal to one because of the property @xmath101 .
\(ii ) if a shock is present , the pressure is spherically symmetric to first order when @xmath125 , the temperature and entropy increase towards the upstream hemisphere , and the density , the velocity and the mach number increase towards the downstream hemisphere
. consequently , more mass is accreted from the downstream hemisphere than from the upstream hemisphere .
this is also confirmed by the numerical simulations with @xmath0 by ruffert ( 1994b ) , at @xmath197 .
this trend persists in the unstable cases , simulated by ruffert & arnett ( 1994 ) at @xmath198 , and ruffert ( 1994b ) , at @xmath199 .
+ thus we disagree with hunt ( 1971 ) , who chose the mach number at the surface of the accretor @xmath196 in order to maximize the mass accretion rate . as can be checked from eq .
( [ m53 ] ) , the highest mass accretion rate is obtained for @xmath192 .
the mach number at the surface of the accretor is , however , not a free parameter , but is constrained by the distribution of entropy .
let @xmath200 be the entropy of the flow line reaching the accretor with @xmath192 .
the pressure being uniform , eq .
( [ pm ] ) implies that the final mach number of the flow line indexed by @xmath66 depends on its entropy @xmath201 according to : @xmath202 the erroneous conclusion of hunt ( 1971 ) , that the axisymmetric flow with @xmath0 never becomes supersonic close to the accretor is in contradiction with the 3d numerical simulations ( see fig .
17 in ruffert  1994b ) .
it is even in contradiction with the mach number @xmath203 that one can deduce from the velocity , entropy and density displayed in figs .  5 , 6 and 8 in hunt ( 1971 ) for @xmath204 .
+ in a similar way as eq .
( [ sbmlies ] ) , obtained for @xmath1 , eq .
( [ machaccr ] ) shows that the functions @xmath205 also cease to be independent at the boundary of the accretor for @xmath0 if @xmath206 . for a given distribution of entropy @xmath207 , the only degree of freedom left by eqs .
( [ mpsi ] ) , ( [ m53 ] ) and ( [ machaccr ] ) is the value of @xmath200 , the amount of mass accreted subsonically .
the value of @xmath208 which maximizes the total mass accretion rate satisfies an integral equation depending only on @xmath207 , which can be solved numerically .
thus if the shape of the shock is known , the distribution @xmath207 is fixed by the rankine - hugoniot jump conditions along the shock ( see appendix  [ aaxis ] ) , and the total mass accretion rate @xmath209 is bounded by this maximum value @xmath210 .
+ since the subsonic accretion of material with a non  uniform entropy is possible for @xmath0 , we only have the inequality @xmath211 .
nevertheless , the difference of entropies in the region of subsonic accretion is limited by the eq .
( [ machaccr ] ) to : @xmath212 an angular sector of passing flow lines may exist , and is described in sect .
[ noaccretion ] . +
an estimate of the first order departures from spherical symmetry is obtained by writing the differential eq .
( [ diffpsi ] ) to first order in the vicinity of @xmath79 ( see appendix  [ adiff ] ) .
this provides us with a lower bound on the departure from sphericity for the function @xmath152 , when the entropy is not uniform : @xmath213 as in sect .
[ sssuper ] , this  bending " of the flow lines enforced by the entropy gradient could be related to the formation of long  lived radial shocks in the supersonic region of flows with a non  uniform entropy , as observed in numerical simulations ( ruffert  1994b ) .
only two situations are known to have analytical solutions :    \(i ) spherically symmetric accretion from a gas at rest at infinity ( bondi  1952 ) .
the mass accretion rate of the unique transsonic solution is : @xmath214    \(ii ) axisymmetric accretion from a gas with negligible temperature ( hoyle and lyttleton  1939 ) : @xmath215 other normalizations were also proposed to account for the pressure effects with more realism , a factor @xmath216 lower ( bondi & hoyle  1944 ) . + in order to bridge the analytical gap between the case of spherical accretion described by @xmath217 ( eq .  [ mb ] ) and
the case of supersonic accretion from a gas of negligible temperature described by @xmath218 ( eq .
[ mhl ] ) , bondi ( 1952 ) proposed a simple interpolation formula , meant to give the correct order of magnitude .
here we use the normalization adopted later by shima ( 1985 ) and ruffert & arnett  ( 1994 ) on the basis of numerical simulations , which is a factor 2 larger : @xmath219 it has the advantage of being very simple and giving a roughly correct limiting behaviour for @xmath220 or @xmath221 .
however , because it is independent of @xmath82 , it is bound to fail to reproduce quantitatively the analytical spherical accretion rate @xmath222 , which varies by a factor @xmath223 between @xmath0 and @xmath3 . moreover , numerical simulations have shown typical discrepancies of a factor 23 with this formula for various mach numbers ( see fig .
[ g5/3 ] and fig .
[ g4/3 ] , and also hunt  1979 ) .
an intermediate situation , also analytically tractable , is the case of a spherical accretion with a constant kinetic energy at infinity : the gas is  thrown " towards the accretor , from a reservoir at infinity where the flow speed is @xmath224 , the density @xmath86 and the sound speed @xmath87 are uniform and non  vanishing , the kinetic energy and the thermal energies are constant at infinity .
86.7 mm by 37.0 mm ( pipes )    to this end , we imagine converging radial  pipes " of constant cross section ( fig .
[ pipes ] ) .
they remain separate up to a minimum radius @xmath225 at which they merge . for a given distance
@xmath225 , if a stationary solution exists , the mass accretion rate imposed at infinity would be : @xmath226 thus @xmath227 is also the radius of a cylinder containing the accreted material in an axisymmetric flow where the transition between the parallel flow at infinity and the spherical flow near the accretor ( and also the decrease of angular momentum along each flow line ) would be artificially driven through such pipes .
this situation allows us to understand how the kinetic energy influences the mass accretion rate , for a gas with an adiabatic index @xmath82 .
+ for spherical flows , the mass accretion rate of the stationary solution is determined by the regularity condition at the sonic radius .
we obtain exactly the same dimensionless quartic equation as in the classic spherical bondi accretion , determining , in the notations of theuns & david ( 1992 ) , @xmath228 as a function of @xmath229 ( see for example their eq .
2.9 ) , except that the velocities @xmath230 and distance @xmath29 are now normalized differently : @xmath231 for spherically symmetric flows , the regularity condition at the sonic radius fixes the mass accretion rate @xmath232 for the unique transsonic solution : by integrating the mass flux ( [ momentum ] ) over the sonic sphere ( @xmath233 ) , we obtain : @xmath234 where @xmath85 is the entropy gained when passing through a spherical shock , if the flow at infinity is supersonic .
+ introducing the entropy @xmath235 of the gas at infinity defined by eq .
( [ defentro ] ) , the spherical mass accretion rate can also be written : @xmath236 _ the mass accretion rate is clearly an increasing function of energy _
@xmath99 , _ and a decreasing function of the total entropy _ @xmath237 .
+ we now consider the two possible cases of flows with and without a shock separately .
according to eq .
( [ mdotsph ] ) when @xmath238 , we expect the mass accretion rate of isentropic flows to increase with the mach number : @xmath239 this increase of the mass accretion rate with the mach number for spherical subsonic flows , is not obvious _ a priori_. an increase of the speed at infinity acts as an additional external pressure on the flow : according to eq .
( [ rspherique ] ) , the sonic radius @xmath97 decreases , and eq .
( [ momentum ] ) shows that the momentum increases .
however , the decrease of the area of the sonic surface , over which the momentum must be integrated to obtain the mass accretion rate , is more than compensated by the increase of momentum , for @xmath1 .
the mass accretion rate is independent of the mach number for @xmath0 .      increasing the velocity of the flow at infinity increases its total energy @xmath99 , but passing through a shock also increases its entropy . according to eq .
( [ mdotsph2 ] ) , these two effects act in opposite directions on the mass accretion rate . whether the mass accretion is eventually increased or decreased
depends on the position of the shock , which can be determined numerically as follows : for high mach numbers , we choose the distance of the pipes , so that the sound velocity prior to the shock @xmath240 equals the sound velocity at infinity @xmath241 .
this corresponds to the position @xmath242 such that the dilution due to gravitational acceleration is balanced by the concentration of flow lines .
we make this particular choice by analogy with the case of axisymmetric accretion with a detached bow shock : we show in appendix  [ ahyperbolae ] that the density enhancement , along the axis of symmetry , before the shock , is smaller than a factor 2 if @xmath243 . + for each radius @xmath29
, the bernoulli equation  ( [ bernoulli ] ) determines the supersonic speed @xmath244 and mach number @xmath245 prior to the shock . using the rankine
hugoniot jump conditions , we can compute the entropy jump @xmath246 and the corresponding subsonic velocity @xmath247 after the shock , if the shock occurs at that radius @xmath29 : @xmath248 for a given entropy @xmath85 and bernoulli constant @xmath99 , we denote the velocity profile of the unique transsonic solution , by @xmath249(r)$ ] .
the unique shock position @xmath250 corresponds to the solution of the equation @xmath251(r_{\rm sh})\;.\ ] ]    86.7 mm by 72.3 mm ( shock )    the position of the shock determined by this procedure is shown in fig .
[ solnum ] , as a function of @xmath82 and @xmath252 . except for values of @xmath82 close to @xmath253 ,
the spherical shock distance is larger than the accretion radius , and thus larger than was observed in numerical simulations of the axisymmetric bhl flow ( typically @xmath254 ) .
+ the spherical mass accretion rate is given by the eq .
( [ mdotsub ] ) diminished by a factor @xmath255 defined in eq .
( [ entrorh ] ) .
+    86.7 mm by 142.0 mm ( bosse )    increasing the kinetic energy of the supersonic flow always decreases ultimately the mass accretion rate below the value @xmath217 obtained for a gas at rest .
between the two opposite tendencies ( entropy and energy ) , the decrease due to entropy will dominate for large mach numbers , since the effective mach number @xmath256 at the shock radius is larger than the mach number at infinity . for @xmath257 , eqs .
( [ mdotsph2 ] ) and ( [ entrorh ] ) imply : @xmath258 moreover , fig .
[ solnum ] and fig .
[ bosse ] show that the shock radius is proportional to the accretion radius for high mach numbers , in other words the ratio @xmath259 is a non  vanishing constant for high mach numbers .
apart from a multiplicative factor @xmath260 depending on @xmath82 , the mass accretion rate in this spherical geometry scales like the mass accretion rate @xmath261 derived by hoyle & lyttleton ( 1939 ) for an axisymmetric flow of negligible pressure .
it is remarkable that the mass accretion rate is asymptotically independent of the value of the sound speed at infinity .
we see from eq .
( [ mdotsph2 ] ) that @xmath217 is the maximum spherical mass accretion rate for @xmath0 .
as we shall proceed to show , the situation is very different for @xmath82 close to 1 .
the value of @xmath82 defines the strength of a pressure change in response to a density change . since the pressure gradient is the only force limiting the mass accretion rate , the higher the value of @xmath82 , the higher the pressure forces and thus the slower the accretion ( see eq .
+ an isothermal gas is much more easily  swallowed " by the accretor : a significant increase of the mass accretion rate occurs for @xmath262 . in this range
, the thermal energy dominates the kinetic energy , and eq .
( [ rspherique ] ) implies that the sonic radius scales like the bondi radius @xmath263 .
more accurately , taking the rankine
hugoniot jump condition ( [ entrorh ] ) in the limit @xmath264 and combining with eq .  ( [ mdotsph2 ] ) gives : @xmath265 we can estimate the contribution of the exponential factor using the bernoulli equation ( [ bernouspher ] ) , for @xmath266 : @xmath267 consequently , the mass accretion rate scales like : @xmath268 using the geometrical constraint that @xmath269 , we can also write the mass accretion rate @xmath270 as the flux of mass through the shock surface , for @xmath266 . for high mach numbers , eq .
( [ bernouspher ] ) ensures that @xmath271 , and the density before the shock @xmath272 because the pipes are of constant cross section .
thus , @xmath273 comparing this equation to eq .
( [ entroenerg ] ) , we conclude from @xmath274 that the shock radius is a factor @xmath275 larger than the sonic radius . inserting this result back into eq .
( [ m1minf ] ) , we further conclude that the mass accretion rate is : @xmath276 which is a factor @xmath277 larger than the spherical mass accretion rate @xmath217 from a gas at rest at infinity , and a factor @xmath278 higher than the hoyle  lyttleton accretion rate @xmath261 .
+ we conclude that for values of @xmath82 approaching 1 , the largest spherical mass accretion rate is : @xmath279 and is reached when the kinetic and thermal energies at infinity are comparable , : @xmath280 this scaling obtained in the isothermal limit for high mach numbers is confirmed by the exact calculation of the mass accretion rate , for any value of @xmath82 , through the numerical derivation of the position of the shock , as shown by fig .
[ bosse ] .
+ the entropy  energy approach for a spherical flow has shown us that the mass accretion rate ultimately decreases like the hoyle
lyttleton rate when the kinetic energy dominates the thermal energy .
+ nevertheless , the mass accretion rate increases considerably above its value at @xmath281 , when the kinetic energy is comparable to the thermal energy , @xmath282 .
this effect is most pronounced for @xmath283 .
thus the value of the adiabatic index plays a very important role in the determination of the mass accretion rate in the case of spherical accretion with a shock .
similarly , we expect the mass accretion rate of an axisymmetric bhl flow with constant parallel velocities at infinity to depend on the adiabatic index .
the axisymmetric mass accretion rate @xmath285 can be formally calculated from the integral of the mass flux along any surface surrounding the accretor .
we can decompose the mass accretion into three components , @xmath286 where @xmath287 is the mass flux through the sonic surface , @xmath288 is the mass flux in the subsonic region ( types sf and fsf ) , and @xmath289 is the mass flux in the supersonic region ahead of the shock ( type fsf only ) . + the trajectories of the gas in the supersonic region ahead of the accretor are well approximated by hyperbolae , obtained by neglecting the temperature of the gas .
@xmath289 decreases to zero when the accretor size @xmath4 decreases to zero ( eddington  1926 ) : @xmath290 according to sect
.  [ ssubsonic ] , @xmath288 also decreases to zero with decreasing size of the accretor , for shocked flows with @xmath1 .
consequently , the mass flux through the sonic surface gives a good estimate of the total mass accretion rate for small accretors : @xmath291 using eq .
( [ momentum ] ) with @xmath75 , we can write @xmath287 as follows : @xmath292 as in the case of spherical accretion , the mass accretion rate is an increasing function of the energy and a decreasing function of the entropy .
however , the geometry of the sonic surface also influences the mass accretion rate . we know from eq .
( [ sonic ] ) that the shape of the sonic surface @xmath293 is directly related to the shape of the flow lines .
in contrast to the energy and the entropy , the geometrical factor depends on the shape of the flow lines at infinity , and is more difficult to estimate since it does not correspond to a conserved quantity .
one can try to obtain , for the axisymmetric case , an interpolation formula for the mass accretion rate @xmath294 , based on the spherical entropy  energy approach .
+ with some trigonometry we write : @xmath295 the direction @xmath120 being in the direction of the flow , the axisymmetric sonic surface is described by @xmath296 , with @xmath297 . for
type f flows , @xmath80 must be replaced by @xmath127 in what follows : @xmath298 first note by differentiating @xmath299 with respect to @xmath43 that : @xmath300 the integral in eq .
( [ mdotfr22 ] ) can be seen as an average of @xmath301 times the entropy weighted by the geometrical function @xmath302 , thus defining an effective entropy corresponding to an effective mach number @xmath303 : @xmath304 using eq .
( [ entrorh ] ) and the definition  ( [ mb ] ) of @xmath305 , eq .
( [ mdotfr22 ] ) becomes : @xmath306 one can check that for isentropic flows ( @xmath238 ) and small deviations from sphericity ( @xmath307 ) , the mass accretion rate is accurately described by the spherical model of sect .
[ sspherical ] : @xmath308 in order to approximate the mass accretion rate @xmath287 , one must estimate the value of the effective mach number @xmath309 , the effective sonic radius @xmath310 and the integral of @xmath311 .
if the shock is detached , we can make simple guesses if the spherical sonic radius @xmath97 is smaller than the accretion radius @xmath312 .
we then use the spherical approximation for @xmath83 and @xmath313 , and use an ad - hoc expression for @xmath309 only , to account for the asymmetry .
we see from fig .
[ r0ra ] that @xmath314 for @xmath315 .
the range @xmath316 , and especially the vicinity of @xmath317 , requires a more careful treatment because @xmath97 can be orders of magnitude larger than @xmath46 .
consequently , we restrict the present approach to the range @xmath315 , and make the simple approximation @xmath318 .
+ since we would like to recover the hoyle
lyttleton accretion rate @xmath261 for high values of the mach number , we see from eq .
( [ asymptot ] ) that this requires a constant value for the ratio @xmath319 when @xmath320 . as a first approximation , we propose the ad  hoc prescription that @xmath321 when @xmath320 , and retain the same ratio @xmath319 for intermediate values of @xmath252 . neglecting the variations of the sound velocity before the shock , we can write the effective mach number in terms of the angle @xmath322 of the  effective flow line " of entropy @xmath323 and the normal to the shock surface , and the distance @xmath324 between the accretor and the intersection of this flow line with the shock surface : @xmath325 a geometrical interpretation is the following .
let the accretion radius @xmath46 act as a scaling factor on the topology of the accreted flow , and keep the angle @xmath322 and the ratio @xmath326 independent of the mach number .
( [ varphie ] ) then implies that the ratio @xmath327 is constant .
this does not contradict the fact that the opening angle of the mach cone , well after the accretor , depends on the mach number , since our approximation deals with the region of the bow shock which is mainly ahead of the accretor .
+ we obtain , for @xmath328 : @xmath329 with @xmath309 defined by this relation if @xmath328 , and @xmath330 otherwise , our interpolation formula follows from eq .
( [ maxme ] ) : @xmath331    86.7 mm by 72.3 mm ( i5/3 )    86.7 mm by 72.3 mm ( i4/3 )    we extend our rough approximation to the case @xmath0 , neglecting the mass accreted in the subsonic region . we know from sect .
[ s53 ] that this corresponds to a small ( but not zero ) fraction of the total flux of mass . + the curves displayed in fig .  [ g5/3 ] and fig
.  [ g4/3 ] show the mass accretion rates obtained with our interpolation formula @xmath294 , the bondi
hoyle interpolation formula @xmath332 and the results of numerical simulations ( ruffert  1994b , 1995 ; ruffert & arnett  1994 ) . for @xmath0 our interpolation eq .
( [ interpol ] ) , in @xmath332 units , rises monotonously with mach number ( fig .
[ g5/3 ] ) .
the numerical results fit only marginally well .
note that the numerical simulations with @xmath0 give only an upper bound for the mass accretion rate , since the accretor is always larger than the spherical sonic radius . for @xmath3 ( fig .
[ g4/3 ] ) a prominent maximum is apparent both on the interpolation formula and on the numerical results . comparing fig .
[ g4/3 ] with the corresponding curve in the lower part of fig .
[ bosse ] gives an indication on the role played by non  sphericity in reducing the mass accretion rate .
+ despite the considerable simplification of ignoring the geometrical term @xmath333 and keeping the ratio @xmath319 independent of @xmath252 , these curves are in reasonable agreement with the numerical simulations for @xmath3 and @xmath0 . thus we expect this analytic formula to be relevant for the whole range @xmath334 .
+ note that we do not solve the question of the dependence of the mass accretion rate with respect to the adiabatic index for high mach numbers .
for this purpose , and following bondi & hoyle ( 1944 ) , we introduced the constant @xmath335 , assumed to be between @xmath336 and @xmath337 .
our interpolation formula is nothing more than a natural way to link the bondi mass accretion rate at @xmath281 with the hoyle - lyttleton mass accretion rate at @xmath320 .
+ one should also keep in mind that the mass accretion rate obtained in the numerical simulations is for most models a time average over the fluctuations or quasi periodic oscillations due to the instability of the flow , and therefore differs from the exact mass accretion rate of the solution of the stationary equations .
this instability is stronger for @xmath82 close to @xmath253 .
by using a geometrical interpretation of the shape of the flux tube at the sonic radius , we have shown the particular significance of the radius @xmath97 ( eq .  [ rspherique ] ) , which depends only on the adiabatic index and the energy of the flow .
it is not only the sonic radius of all spherical flows , but is also related through the property @xmath101 ( sect .  [ ssprop ] ) to the sonic surface of axisymmetric flows .
+ by reducing the stationary flow equations to a single partial differential equation , we have performed a local analysis near the accretor and have obtained some general results which aim at clarifying the diversity of possible configurations .
we have extended the classification introduced by bondi ( 1952 ) to axisymmetric stationary flows .
in addition to the subsonic and supersonic types of accretion , angular sectors without any accretion are also possible , as well as singular directions of infinitely high accretion flux .
this is illustrated by the self - similar solution found by bisnovatyi - kogan ( 1979 ) , which is a very particular singular case ( @xmath338 ) .
+ the angular sectors without any accretion is subsonic for flows with a detached shock , so that @xmath339 in the notation introduced in sect .  [ notations ] .
+ among some of the properties that we established for flows with @xmath1 are the following :    \(i ) the sonic surface is likely to be attached to the accretor for nearly isothermal flows with high mach numbers ,    \(ii ) isentropic accretion proceeds mostly from the downstream hemisphere of the accretor ,    \(iii ) subsonic accretion must be isentropic .
+ we have also found a series of properties for flows with @xmath0 :    \(i ) accretion is always regular ,    \(ii ) the sonic surface is always attached to the accretor ,    \(iii ) the pressure distribution is spherically symmetric to first order ,    \(iv ) the mach number and the entropy at a point like accretor are algebraically related to the local mass flux ,    \(v ) shocked matter is accreted mostly from the downstream hemisphere , where the density , velocity and mach number are maximum , and the temperature and entropy are minimum ,    \(vi ) isentropic accretion is spherically symmetric to first order for the velocity , temperature , density , mass flux , and mach number .
+ 3d numerical simulations have shown that the strongest instabilities occur with the adiabatic index @xmath0 .
the fact that subsonic regions invariably reach the accretor for @xmath0 might be considered as a warning about the possible role , in the instability mechanism , of the boundary conditions at the surface of the accretor . for @xmath1 and sufficiently small accretors , it would be interesting to check for the presence of the instability in cases where the accretor is surrounded by a supersonic region ( type  f in the nomenclature of sect .  [ sstypes ] ) .
+ we have also determined the departure from sphericity , and constrained the radial variation of the angle @xmath152 for @xmath119 . in particular , we have stressed the influence of the entropy gradient on the bending of the flow lines for supersonic flows with @xmath1 and @xmath0 .
+ we have stressed the role of the geometrical shape of the shock in the determination of the entropy of the accreted matter , and contrasted it with the simpler case of shocked spherical flows . in this latter case , we have shown the relevance of the entropy - energy approach , and noted that the mass accretion rate scales like the hoyle
lyttleton formula for high mach numbers .
moreover , we have outlined the particularity of nearly isothermal flows for which the mass accretion rate can be a factor @xmath340 higher than the bondi mass accretion rate , when the kinetic energy at infinity is comparable to the thermal energy .
this demonstrates the fundamental difference of the bhl flows that have @xmath82 close to 1 and those with @xmath82 close to 5/3 .
+ we have shown that the total mass accretion rate equals the mass accretion rate through the sonic surface , for shocked regular flows with @xmath1 .
we have derived an interpolation formula for the mass accretion rate valid in the range @xmath2 .
it links spherical accretion to the case of highly supersonic accretion , with reasonable agreement with numerical simulations .
+ analytical estimates of the following quantities , however , are still missing :    \(i ) the distance of the shock as a function of @xmath82 and @xmath252 ,    \(ii ) the value of @xmath80 and @xmath124 as a function of @xmath82 and @xmath252 ,
\(iii ) the mass accretion rate in the limit of infinite mach numbers , as a function of @xmath82 .
nevertheless , our classification will help perform , in a forthcoming paper , a local stability analysis of all the configurations considered .
tf was supported by the ec grant erb - chrx - ct93 - 0329 , as part of the research network ` accretion onto compact objects and protostars ' .
we acknowledge useful comments by dr .
h. spruit and dr .
the combined vorticity equation , continuity equation and entropy equation are written as functions of @xmath341 as follows : @xmath342
using the axisymmetry condition and some simple trigonometry , the variation of the cross  section @xmath83 along the flow line can be written : @xmath343 where @xmath38 is defined by eq .
( [ defbeta ] ) , and its derivative is taken in the direction @xmath344 perpendicular to the flow line ( oriented towards increasing @xmath60 ) .
for @xmath345 we recover the spherical case of conical flow tubes , with a cross section proportional to the square of the distance .
+ we can use eq .
( [ trigo ] ) along the detached sonic surface , together with eq .
( [ sonic ] ) , for @xmath346 : @xmath347 let us now prove that the sonic radius reaches the value of the spherical radius at least once , that @xmath348 has to vanish .
+ axisymmetry imposes @xmath349 , and by eq .
( [ deff ] ) : @xmath350 suppose , for example , that the sonic surface is inside the sphere of radius @xmath97 and detached from the accretor , @xmath351 for @xmath352 $ ] . equations ( [ f0 ] ) and ( [ fpi ] ) imply that @xmath38 is strictly positive near @xmath353 , and strictly negative near @xmath354
. a line must exist in the ( @xmath54 ) plane separating @xmath355 from @xmath356 , along which @xmath152 vanishes , when decreasing from a region of positive values to the region of negative values .
we denote by @xmath357 the intersection of this line with the sonic surface : @xmath358 which contradicts the hypothesis .
the same contradiction arises if we suppose that @xmath359 .
consequently @xmath360 must vanish for some value of @xmath60 .
by using eq .  ( [ tanbeta ] ) , we can relate the second radial derivative of @xmath63 to its first radial derivative , at @xmath79 , through the following equation : @xmath362 the local mass flux is finite for regular flows , and therefore the limit when @xmath119 of ( @xmath363 ) is zero in the angular sectors where the flow is regular .
+ let us define the new radial variable @xmath364 .
equation  ( [ implicit ] ) implies to first order : @xmath365    \(i ) if eq .
( [ betamax ] ) is a strict inequality , eq .
( [ eulert ] ) implies that the pressure is spherically symmetric to first order .
using eqs .
( [ mpsi ] ) and ( [ psup ] ) , the stream function @xmath366 at the surface of the accretor is directly related to the entropy @xmath367 by : @xmath368 the leading terms of the differential eq .
( [ diffpsi ] ) satisfy the following equation : @xmath369 where the six functions @xmath370 converge to zero when @xmath371 .
the functions @xmath372 depend on @xmath60 only , and are defined as follows : @xmath373 assuming @xmath374 and @xmath375 contradicts eq .
( [ linearinf53 ] ) .
thus we conclude that : @xmath376 by derivating eq .
( [ diffpsi ] ) with respect to @xmath29 and considering the leading terms , we find that @xmath377 satisfies the same eq .
( [ linearinf53 ] ) as @xmath378 does .
using eq .
( [ tanbeta ] ) , we obtain eq .
( [ betamin ] ) for regular sectors with a non  uniform entropy .
\(ii ) if eq .
( [ betamax ] ) is an equality , we use eq .
( [ secondprem ] ) in order to write the first order terms of eq .
( [ diffpsi ] ) as follows : @xmath379 in both cases ( i ) and ( ii ) , @xmath380 and the radial derivative @xmath381 cease to be a set of three free functions of @xmath60 on the surface of the accretor , when the accretor size decreases to zero . let us recall that for spherical flows , one parameter ( the mass flux ) is required as boundary condition on the surface @xmath382 of the accretor , with no particular change if @xmath383 . in this respect ,
the differential system ( [ diffpsi ] ) is more singular at @xmath79 for axisymmetric flows than for spherical flows .
because the stream function @xmath63 is finite everywhere , @xmath384 and @xmath385 tend to zero with decreasing @xmath43 .
( [ implicit ] ) implies that when @xmath0 , the mach number tends to a finite limit everywhere where the local mass flux ( @xmath386 according to eq .
( [ mpsi ] ) ) is not zero . taking the limit of eq .
( [ diffpsi ] ) when @xmath387 , we find that its last term must vanish at @xmath388 : @xmath389 where the subscript @xmath148 denotes the value of the function at @xmath79 .
this equation is equivalent to the condition that the pressure is spherically symmetric to first order .
+ in order to find the behaviour of the quantities @xmath390 near @xmath388 , we consider only the leading terms of eqs .
( [ implicit ] ) and ( [ diffpsi ] ) .
we first define @xmath256 and @xmath378 as follows : @xmath391 by linearizing eq .
( [ implicit ] ) near @xmath392 , we obtain to first order : @xmath393 introducing this into eq .
( [ diffpsi ] ) , we obtain the following leading terms : @xmath394 where the six functions @xmath395 converge to zero when @xmath387 . the functions @xmath396 depend on @xmath60 only , and are defined as follows : @xmath397 let us suppose that @xmath84 is a positive number such that @xmath398 if the entropy is not uniform ( @xmath399 ) , the first two terms of eq .
( [ linear53 ] ) are negligible compared to the rhs term . dividing the remaining equation by @xmath43 , we obtain : @xmath400 according to eq .
( [ hypolim ] ) , @xmath401 .
thus eq .
( [ psi1r ] ) contradicts our hypothesis ( [ hypolim ] ) , and we conclude that if @xmath399 , @xmath402 by derivating eq .
( [ diffpsi ] ) with respect to @xmath43 , and keeping the leading terms , @xmath403 satisfies exactly the same equation ( [ linear53 ] ) as @xmath378 does . from the same reasoning
, we conclude that @xmath404 this lower bound on the scaling of @xmath403 is also a lower bound on the scaling of @xmath152 : we obtain eq .
( [ dlogbeta ] ) by using eqs .
( [ tanbeta ] ) , ( [ conclulim2 ] ) and the regularity of the flow for @xmath0 .
we check the validity of neglecting the pressure forces in the supersonic flow ahead of the shock , and assuming that the sound velocity is uniformly constant to its value at infinity @xmath87 .
the velocity of the supersonic incoming gas is then defined by the unique hyperbolic trajectory passing trough this point , with velocity at infinity @xmath224 . + for each point ( @xmath54 ) of the shock surface , we can compute numerically the dimensionless impact parameter @xmath405 ( normalized to the accretion radius @xmath46 ) of this hyperbola : @xmath406 the half angle @xmath407 of the hyperbola is defined by @xmath408 let us calculate the variation of density along the axis of symmetry , due to both the compression by gravitational convergence of the hyperbolic trajectories and the dilation due to acceleration .
this can be done using the continuity equation  ( [ continuity ] ) in cylindrical coordinates , along the axis , for @xmath409 : @xmath410 the azimuthal derivative of the azimuthal component of the velocity is directly related to the impact parameter @xmath411 through the conservation of angular momentum , while the ( negative ) radial component of the velocity , along the axis , is given by the conservation of energy : @xmath412 in order to determine the azimuthal dependence of the impact parameter , let us consider the trajectories defined by eq .
( [ hyperbole ] ) for @xmath413 . since @xmath414
, we write : @xmath415 from first order expansions of eq .
( [ demiangle ] ) and eq .
( [ hyperbole ] ) , we find , for @xmath413 : @xmath416 choosing the negative solution of this equation ( @xmath405 is positive and decreases to zero when @xmath60 increases to @xmath127 ) , leads to : @xmath417 using eqs .
( [ dvtheta ] ) and  ( [ enervr ] ) in eq .
( [ contax ] ) we obtain : @xmath418 this can be integrated to obtain the radial dependence of the density along the axis of symmetry : @xmath419 the compression effect due to the convergence of the flow lines is consequently stronger than the dilatation due to the acceleration of the flow .
+ if the shock is not too close to the accretor , gravitational focusing is negligible and the density and temperature can be approximated as constant in the upstream supersonic flow .
one can check from eq .
( [ gravcomp ] ) that this density enhancement is smaller than a factor of two for @xmath420 . note that when the shock is detached , its distance is generally of order @xmath421 in numerical simulations .
we deduce from eq .
( [ gravcomp ] ) that the sound velocity scales like @xmath422 when @xmath423 , and is therefore always negligible compared to the gravitational potential scaling like @xmath424 .
let @xmath425 be local coordinates in a frame chosen in a way such that @xmath426 is perpendicular to the shock surface , and @xmath427 is parallel to it .
let us call @xmath256 the mach number relative to the component of the velocity , ahead of and perpendicular to the shock surface . in our local system of coordinates , @xmath428 we shall use the classical convention of using the index  1 " for quantities ahead of the shock , and the index  2 " for quantities just after the shock .
+ we write the entropy gradient after the shock as : @xmath429 we can use the classical jump conditions  ( [ entrorh ] ) to express the entropy gradient as a function of @xmath256 only : @xmath430    we obtain @xmath431 by the vector product of the velocity and the vector tangent to the shock surface .
let the supersonic flow before the shock be described in cylindrical coordinates by the flow velocity @xmath432 .
for the sake of clarity , we denote the velocity along the shock surface @xmath433 by @xmath434 and a dot stands for the derivative with respect to the angle @xmath60 : @xmath435 we now prove that a regular axisymmetric shock surface can not produce a local region of uniform entropy behind it , and more precisely , that the gradient of entropy created by the shock scales like the square of the azimuthal angle @xmath436 in the vicinity of the axis , for any shape of the shock which is regular and symmetric .
+ when not specified , the quantities are taken at @xmath160 in what follows .
we also define @xmath437 and @xmath438 .
the symmetry condition imposes : @xmath439 using the symmetry of the problem , we make the following expansions for small values of @xmath436 : @xmath440    using the expansions ( [ expvr ] ) , ( [ expvt ] ) , and ( [ exprs ] ) , we obtain : @xmath441 in order to evaluate this expression , we use the approximation of the supersonic flow @xmath432 by hyperbolic trajectories . derivating the bernoulli eq .
( [ bernoulli ] ) with respect to @xmath442 and @xmath60 , we obtain for @xmath160 : @xmath443 replacing these expressions , together with eqs .
( [ enervr ] ) and ( [ dvtheta ] ) into eq  ( [ vnbrut ] ) , the effective mach number is : @xmath444 where the positive contribution @xmath445 depends on the second derivative @xmath446 of the shock surface on the symmetry axis , and is defined as : @xmath447 we conclude that _ for any shape of the detached shock surface _ , the effective mach number near the symmetry axis decreases at least as the square of the azimuthal angle @xmath448 .
+ this result is particularly useful in order to prove that the mass accretion rate in the subsonic region of type sf flows decreases to zero with the size of the accretor .