the current successes of cold dark matter cosmologies , and in particular the @xmath1cdm models , are well known .
these hierarchical models begin with the microwave background fluctuations and proceed forward to reproduce the structure in the ly@xmath2 forest at @xmath3 ( e.g. , mandelbaum et al . 2003 ) and continue to evolve the dark matter structure distribution to arrive at results consistent with the power spectrum of the distribution of galaxies at the present epoch ( e.g. , tegmark et al .
2004 ) , the fraction of baryons in rich clusters ( white et al .
1993 ) , and other notable results .
these models follow the evolution of dark matter , and the accompanying galaxies are assumed to form as gas flows into the ever growing dark matter condensations ( e.g. , white and rees 1978 ) .
the most detailed results describing this structure evolution have come from highly sophisticated numerical simulations that follow the growth of mass fluctuations in the early universe , such as the millennium simulation ( springel et al .
2005 ) . however , there are some well known difficulties with these models in that they predict too many faint and too many bright galaxies ( e.g. , benson et al .
2003 ) , and in addition they predict that the largest galaxies should be blue and star forming .
these are serious difficulties , and they have been addressed through the use of various assumptions about physical processes that work on scales not resolved by the numerical simulations .
these `` semi - analytic '' models have addressed the problems at both the dwarf galaxy and the brightest galaxy scales ; the issues addressed in this paper involve only the problems associated with the evolution of the brightest galaxies .
the problems there are a direct result of the hierarchical nature of the models ; minute fluctuations in the very early universe grow due to continuing dark matter and baryonic infall , and in the context of the models this infall should continue to the present epoch , possibly leading to significant star formation now in the most massive galaxies .
thus in this picture the most massive galaxies would be blue due to the large population of young stars , whereas in fact the largest galaxies are generally seen as `` red and dead '' at the present epoch , implying that the most active epoch of star formation occurred earlier ( e.g. , madau et al .
1996 , cowie et al .
1996 ) .    to address this problem ,
one of the most commonly used assumptions in semi - analytic models is that of `` feedback '' from the central regions of the galaxy , where energy in some form is produced in sufficient quantity that it can heat the inflowing gas , perhaps stopping or reversing the inflow , and in any case providing enough energy to suppress star formation .
once this has been done , the subject galaxy will then passively evolve , and its aging stellar population will become redder with time until the galaxy colors match the observations . motivated by the observation of active galactic nuclei ( agn ) in many massive red galaxies , several agn related driving mechanisms for this feedback have been suggested .
radiatively driven winds , agn driven shocks and mass outflows , qsos , and radio sources associated with agns have all been proposed ( e.g. , sijacki et al .
2007 ) . in some instances qsos
have been referred to as `` high luminosity '' agn feedback and radio - loud agn ( or `` radio agn '' ) have been called `` low luminosity '' agn feedback .
this may not be a meaningful distinction , since many qsos may well be examples of powerful radio sources that are directed toward the observer , with their collimated emission doppler boosted to very high luminosities . the use of agns , and in particular radio agns , as sources of the energy required to suppress star formation in feedback models has appeal because estimates for the total energies of these objects derived from energy equipartition arguments have existed for some time .
these minimum energies ( from @xmath4 ergs ) are , in an integral sense , adequate to provide the needed heating for agn feedback to work .
moreover , recent estimates of the total energies present in some radio sources tend to lie factors of 10 higher than the estimates based on the equipartition energy present in relativistic electrons and fields ( brzan et al .
2008 ) .
thus the total energies available in radio agn make these objects viable candidates for powering feedback and for suppressing star formation in massive galaxies , and radio agn have been incorporated into semi - analytic models that are based on cosmological simulations such as the millennium run ( e.g. , croton et al .
2006 , bower et al .
2006 , okamoto et al .
for example , croton et al .
( 2006 ) show that a form of radio agn energy input can not only reproduce the observed truncation of the bright end of the luminosity function , but that it can also possibly suppress accretion and give rise to an old , red population of elliptical galaxies .
radio agn feedback has also been used to explain the presence of old `` red and dead '' galaxies in voids ( croton & farrar 2008 ) , where environmental effects are much less pronounced . finally , radio agn feedback has also been suggested by several authors as a means of suppressing cooling or providing reheating in the cores of rich clusters of galaxies , thereby solving the `` cooling flow '' problem and at the same time accounting for the colors of brightest cluster galaxies ( e.g , fabian et al .
2000 , brggen & kaiser 2002 , ruszkowski et al .
2004 , mcnamara & nulsen 2007 ) .    while the integrated total energies of radio sources are adequate , and sometimes more than adequate ( e.g. , best et al .
2006 ) to provide the needed feedback in semi - analytic models , the key question of `` does it really work ? '' remains .
having enough total energy is only the first of several conditions that need to be met . in order for radio agn feedback to work ,
this energy contained in the outflowing radio sources must be transferred to the surrounding medium . in order to suppress star formation
, this energy transfer has to be fairly efficient , and equally important , for radio agn feedback it has to be turned into heat , and it has to be distributed over most of the volume of gas in and around the halo of the parent galaxy . finally , all this has to happen in a time less than the local cooling time or dynamical infall time .
until all of these conditions can be met , the effectiveness of radio agn feedback in solving significant difficulties with current cosmological models will remain in question .
this paper provides a calculation of one way this energy transfer process can work .
section ii presents an overview of the physical processes involved in radio agn feedback , while section iii describes a calculation of the energy transfer during the interaction of a radio source with its environment that covers up to 5 orders of magnitude in scale .
the results of this calculation are given in section iv , and the conclusions are found in section v.
a fundamental property of radio agn that has a major impact on the issue of feedback is the bipolar nature of the energy outflow . except for low energy `` compact radio sources '' , all radio agn show this morphology .
the degree of collimation can be moderate , with opening angles of the flow of order twenty degrees or so , or it can be extreme , as in the case of the very energetic and highly collimated flows associated with high luminosity radio sources and quasars .
this two - sided collimated outflow is the principal obstacle to radio agn feedback being an efficient and effective mechanism for energizing and heating ( or blowing away ) the gas in the central regions of the parent galaxy , since the heating and suppression of star formation must occur throughout the central region and not just along the two narrow cones in the vicinity of the radio jets .
thus for feedback to be effective , some processes must be found that transfers energy from the jet to the surrounding medium , converts that energy to heat , and distributes it over most of the volume in the central regions of the parent galaxy . some possible mechanisms for this energy transfer
have been previously suggested , such as shock waves from highly supersonic outflows ( brggen et al .
2007 , graham et al . 2008 ) or sound wave dissipation from less energetic events ( sanders & fabian 2007 , 2008 ) . these will be considered in more detail later , but in general shock waves are one - time events that rapidly slow and dissipate in the direction transverse to jet motion , while sound waves may or may not carry enough energy to heat large volumes of gas through their eventual dissipation .
more importantly , the process of energy deposition by the sound waves in a magnetoionic medium is not clearly understood in this case , especially with regard to the dissipation lengths .
clarification of possible feedback processes can be found from the luminosities , morphologies and demographics of the radio agn outflows .
most extended radio sources can be placed in one of two categories first proposed by fanaroff & riley ( fanaroff & riley 1974 ) .
these two classes , fr - i and fr - ii , are characterized by different luminosities and by very different radio morphologies .
high luminosity fr - ii objects are generally characterized by two very highly collimated jets of emission extending over large distances , typically hundreds of kpc , with the jets terminating in regions of high surface brightness ( `` hot spots '' ) surrounded by lower surface brightness lobes which often have a `` swept back '' or `` bow shock '' appearance ( e.g. ferrari 1998 , de young 2002 ) .
the morphologies of the lower luminosity fr - i objects are dramatically different , with the highest surface brightness near the nucleus of the parent galaxy and a gradual dimming of brightness along two moderately collimated jets with much wider opening angles than seen in the fr - ii objects .
these outflows are also often bent or distorted and display a meandering morphology as the flow moves outward .
this pattern of radio emission is very suggestive of the slow , subsonic flow associated with plumes and turbulent flows , and it was suggested many years ago ( e.g. , de young 1981 , bicknell 1984 ) that these radio sources were characterized by significant interaction with the surrounding ism or igm that resulted in mass entrainment and a deceleration of the agn outflow as it moves away from the nucleus
. it may also be that the original outflow speeds in the fr - i objects are less than those existing in fr - ii sources , but in any case the current modeling of these objects usually involves deceleration to subsonic flow , presumably mediated by mass entrainment from the ambient medium ( e.g. , laing & bridle 2002 , 2004 ) .
the second key point in considering radio agn feedback comes from the demographics of fr - i and fr - ii radio sources .
it has been known for some time that the local space density of extended radio sources is dominated by fr - i objects ( e.g. , ledlow & owen 1996 , owen & ledlow 1997 , jackson & wall 1999 , willot et al .
2001 ) , where the local ( @xmath5 ) density of fr - i objects is roughly @xmath6 times that of fr - ii radio galaxies . the local density of fr - i radio sources , and a calibration point for constructing radio luminosity functions ( rlf ) is about 400 gpc@xmath7 ( rigby , best & snellen 2008 ) . due to the rapid decline in the detection of the fainter fr - i objects with distance , the dependence of the fr - i / fr - ii ratio with redshift is less clear , but recent surveys and modeling of the rlf out to moderate redshifts of less than 1 suggest that the space density of fr - i objects clearly dominates in the local universe ( willot et al .
2001 , cara & lister 2008 , donoso et al .
thus to a very good approximation the `` generic '' radio agn outflow in the nearby universe can be taken to be an fr - i radio source .
this may also be the case at moderate redshifts as well , since the previous picture of a constant fr - i space density with cosmic epoch coupled with a rapidly evolving fr - ii space density may be overly simple ( e.g. , rigby , best & snellen 2008 ) . in any case
it is clear that `` radio agn feedback '' for redshifts less than one can be taken as `` fr - i agn feedback '' , and henceforth the focus here will be on fr - i objects as sources of radio agn feedback .    a third observational result that constrains the feedback process comes from radio and x - ray observations of radio sources in rich clusters of galaxies . for the first time since their discovery over thirty years ago , it is now possible to carry out calorimetry of some radio source outflows by using estimates of the @xmath8 work done by the outflows in inflating cavities observed in the icm .
( e.g. , brzan et al .
2008 , mcnamara & nulsen 2007 and references therein . )
the role of such cavities as being a signature of radio agn feedback in clusters has been discussed elsewhere , and the relation of this phenomenon to more general fr - i outflows will be discussed at the end of section 4 . the focus here is on the general problem of radio agn feedback , both in and out of clusters and groups of galaxies , and the relevant result from the radio source calorimetry in clusters is that it reveals outflow energies that are generally in excess of the equipartition values that were commonly assumed in the past ( brzan et al .
2008 ) . these @xmath8 values are often an order of magnitude greater than the total energies contained in relativistic electrons and magnetic fields in minimum energy equipartition , and they are consistent with most of the energy residing in the kinetic energy flux of the outflow ( e.g. , de young 2006 ) .    with these three observational constraints - the meandering plume - like morphology of fr - i radio sources , the dominance of fr - i objects in the nearby radio agn population , and the existence of much more energy in the outflows than is required to account for the synchrotron radio emission ,
the most likely paths for radio agn feedback become more clearly defined .
first , if these outflows are truly subsonic over most of their length , and the morphological data are very suggestive of this , then heating by shock waves is not likely to be a significant energy transfer mechanism in these objects , nor is momentum transfer likely to form a nearly spherical outflowing wind that removes the ambient gas .
second , because there is strong evidence that agn outflows are relativistic on very small ( @xmath9 parsec ) scales , the fr - i morphology strongly suggests significant deceleration of the flow .
the uniform appearance of the outflows from the nucleus outward without any clear discontinuities together with the spreading of the outflow from the nucleus outward implies that this deceleration process begins near the nucleus and continues more or less uniformly throughout the outflow region within the parent galaxy .
deceleration of the flow means a loss of momentum from the jet via a transfer of momentum from the jet to another medium , with the most obvious candidate being the ambient gas in the ism of the host galaxy .
the most natural and most obvious candidate for this deceleration process is thus entrainment of , and transfer of momentum to , the ambient medium at the interface between the outflowing agn material and the ism and galactic halo gas .
in addition to continuous deceleration from entrainment , deceleration of the outflow from its original speed near the galactic nucleus could entail the production of shock waves and subsequent heating , since vlbi observations of agn suggest very rapid outflow speeds on the parsec and subparsec scales in the innermost regions .
however , the propagation of a decelerating supersonic jet into the circumnuclear gas will produce a driven bow - shock configuration that , for uniform flow , is a one - time event .
in addition , except for the small angular region near the head of the jet , such shocks are generally weak , oblique , and decelerating .
nonetheless , they probably occur around the initial outflow regions of the jets , and as such they could produce localized heating in the innermost regions of the parent galaxy , especially if the agn outflow involves multiple outbursts .
this process may be at work in outflows in rich clusters such as 3c 84 , where deceleration is most dramatic due to the very high ambient densities .
because the total energy involved in these shocks is much less than the kinetic energy of the outflow , and because the focus here is not on agn in rich clusters , the emphasis here will be on the transport of energy by the outflow itself , with most of the flow deceleration being due to continuous entrainment of the ambient medium .
this entrainment process is mediated through the non - linear development of surface instabilities , most notably the kelvin - helmholtz ( kh ) instability , and it is driven by the large ( @xmath10 erg s@xmath7 ) kinetic energy flux of these outflows derived from the calorimetry described above .
analytic approximations , laboratory experiments , and numerical simulations have all shown that this instability in shearing flows is essentially inevitable .
incompressible flows , compressible flows , supersonic flows , mhd flows and relativistic shearing flows all show the onset and non - linear growth of the kh instability ( e.g. , chandrasekhar 1961 , brown & roshko 1974 , clemens & mungal 1995 , aloy et al , 1999 , ryu , jones & frank 2000 , perucho et al .
2004 , 2007 ) . in the fully non - linear stage of development
the instability evolves into a mixing layer dominated by large scale vortex structures that entrain material from both sides of the layer into the layer itself ( e.g. , dimotakis & brown 1976 ) .
the thickness of this layer grows as one moves down the jet from the onset of the non - linear phase of the instability , and the growth of the layer is nearly linear with distance along the jet until the jet is fully infiltrated by the layer . at that point
a different evolutionary phase begins , as is described in section 3 . for a plane mixing layer , which is a good approximation to the cylindrical mixing layer in a round jet in its initial stages ( e.g. , dimotakis et al .
1983 , freund et al . 2000 ) , the initial growth is characterized by a nearly constant opening angle that depends upon the relative flow speeds and the fluid densities as @xmath11 where @xmath12 and @xmath13 are the densities in the outflow and the ambient medium respectively , and @xmath14 is the mach number of the outflow relative to the sound speed in the ambient medium ( e.g. , brown & roshko 1974 , dimotakis & brown 1976 , de young 1993 ) .
empirical data suggest @xmath15 and @xmath16 . within this
growing mixing layer the entrainment proceeds through the initial production of large scale eddies whose diameter is roughly equal to the thickness of the layer , followed by the development of finer scale structures within the eddies as well as merging of the larger scale structures . in other words ,
the mixing layer becomes fully turbulent .
moreover , the thickness of this turbulent mixing layer grows as one moves downstream along the jet , and eventually the turbulence will penetrate throughout the entire volume of the outflow . from eq .
1 , the distance along the jet from the onset of the non - linear mixing layer to the point of complete infiltration of turbulence into the interior of the jet can be approximated by @xmath17 where @xmath18 is the radius of the jet at the onset of the mixing layer . beyond this point
the flow within the jet is fully turbulent ; many experimental examples of such turbulent , decelerated jets can be found ( e.g. , dimotakis et al . 1983 , mungal et al .
1992 ) , and within the limits of resolution the outflow from fr - i objects such as 3c 31 is indistinguishable from these experimental results .
thus the fr - i demographics , when combined with the evidence for uniform deceleration in their agn outflows , the universality of shear driven surface instabilities that entrain mass and provide momentum transfer through the growth of turbulent mixing layers , and the morphological similarities between fr - i outflows and subsonic turbulent plumes all provide evidence that the most common energy outflow associated with radio agn feedback is that of diverging , subsonic , fully turbulent plumes .    this argument is given further support from a comparison of high resolution radio observations of fr - i outflows with experimental results concerning fully turbulent round jets propagating into a uniform medium .
there is considerable evidence that such jets are self - similar flows ( e.g. , pope 2000 , hussein et al .
1994 ) , and that these flows all have the same opening angle once self - similarity is established .
this angle is about 23 - 24 degrees , and an examination of radio observations of the fr - i objects 3c 31 and 3c 296 ( hardcastle et al . 2005 , laing et al . 2008 ) shows that both the north and south jets in these objects are very well enclosed by an opening angle of about 23 degrees ( @xmath19 ) .
once the outflows show sudden divergence or bending this is no longer the case , but in the inner regions of the flow where the ambient conditions are likely to be fairly constant or slowly varying , the agreement with self - similar turbulent jets is remarkable .
elements of this picture have been suggested before , and the first numerical simulation of mass entrainment in astrophysical jets was performed many years ago ( de young 1986 ) ; there has since been considerable development of phenomenological models of these flows ( e.g. , laing et al .
1999 , laing & bridle 2002 , 2004 , wang et al .
2009 ) .
because magnetic fields are present in both the agn outflows and in the ambient medium , possible inhibition of the mixing layer growth by magnetic effects might occur , especially for `` favorable '' field geometries where @xmath20 is parallel to @xmath21 .
however , it is clear that in the context of mhd outflow models , the average kinetic energy density must be much greater than the average magnetic energy density in order for the agn outflow to occur .
( the opposite case of nearly mass - free `` poynting '' jets is discussed in the last section . ) in the mhd case the concern is then that very local b field amplification due to turbulence could lead to local field energies sufficient to damp the development of the mixing layer . in this case , detailed 3d numerical mhd simulations of the development of the k - h instability and its fully non - linear growth ( e.g. , ryu et al . 2000 ) show that the flow remains `` essentially hydrodynamic '' .
small scale field amplifications do occur , and at late times the amplitudes are comparable to the flow energy densities , but these areas are very localized and do not inhibit the growth of turbulent flow in the mixing layer on the larger scales of energy injection into the turbulence .
it is important to note the difference between these dynamic , driven outflows and the much slower , passive evolution of buoyant " bubbles in an intracluster medium , where large scale icm fields may have an influence on the bubble evolution .
this characterization of radio agn feedback as turbulent fr - i outflows provides the initial conditions for finding processes that can transform this outflow into widespread heating of the ambient medium .
a major advantage of turbulent flows is that the ultimate fate of fully developed turbulent flows is conversion of the turbulent energy into heat .
thus if mass entrainment and the generation of turbulence is the most efficient coupling of the agn outflow to the ambient medium , then the the natural evolution of that turbulence will provide the most effective heating of the gas in the galaxy by radio agn outflows .
a central issue is how long this conversion takes .
the next section provides a quantitative estimate of this time through a calculation of the evolution of a driven turbulent flow from its onset until its dissipation .
the turbulent flow that describes the radio agn outflow can be recast in a form common to many turbulent flows . on some large scale comparable to a defining scale of the problem ( e.g. ,
the radius of the outflow ) energy is injected in the form of turbulent eddies .
this process , as verified by experiment , involves a large scale ingestion of material on both sides of flow boundary in a `` gulping '' mode by the largest eddies ( e.g. , dimotakis & brown 1976 ) , followed by a rapid generation of finer and finer structure and mixing within this large scale as the flow evolves to smaller scales . in terms of the evolution of the energy spectrum of the turbulence , this can be described as energy injection over a small range of wavenumber @xmath22 followed by a loss free cascade of energy to larger wavenumbers . for homogeneous and isotropic hydrodynamic turbulence
this process results in the standard power - law kolmogorov spectrum with @xmath23 .
( for mhd turbulence the power - law exponent is less well defined but is generally somewhat smaller . )
the dissipation of the turbulent kinetic energy into heat occurs at and above some wavenumber @xmath24 which generally lies many orders of magnitude above the energy injection scale @xmath25 .
the kolmogorov formalism describes an equilibrium state where the energy input at @xmath25 is exactly balanced by the dissipation at @xmath24 .
the problem being addressed here is not as simple , since what is needed is a calculation of the time required for the flow to reach the dissipation range at @xmath24 after the injection has begun at @xmath25 . because the evolution of the turbulent flow involves the non - linear interaction and energy transfer among turbulent scales of many different sizes , both upward and downward in @xmath22 space , the evolution must be calculated explicitly .
in addition to the time dependent nature of the flow , a major hurdle is the many orders of magnitude in scale between the size of the large eddy injection range ( tens of parsecs or more ) to the energy dissipation scale due to effective viscosities in ionized plasma ( possibly meters or less ) .
some simplifying assumptions can be made .
if the overall outflow of the fr - i agns is subsonic or transonic relative to the ambient medium , then the flow can be regarded as generally incompressible .
in addition , in a frame comoving with the mean local outflow speed the turbulence can be treated as locally homogeneous and weakly isotropic ( in that helicity is allowed ) to a good approximation . under these conditions
a technique may be used that calculates the three dimensional time dependent evolution of turbulence , both hydrodynamic and magnetohydrodynamic , over several orders of magnitude in spatial scale .
the original method was devised to treat purely hydrodynamical turbulence ( e.g. , orszag 1970 , 1972 ) but was subsequently generalized to include magnetic fields in the mhd approximation ( pouquet et al .
1976 , pouquet et al . 1978 ,
de young 1980 ) .
this approach takes moments of the fourier transforms of the equations of mass , momentum and energy conservation in order to obtain equations that describe the time evolution of second moments of velocity and magnetic field , i.e. , the kinetic energy and magnetic energy , as a function of time and spatial scale .
the system of moments is closed by using a quasi - normal approximation that assumes the fourth order moments are related to the second order moments in the same way as is the case for a gaussian distribution of fourier modes .
this `` quasi - normal '' approximation permits closure of the system of equations at second order and results in differential equations that describe the temporal and spatial evolution in fourier space of the square of the fluid velocity ( i.e. , kinetic energy per unit mass ) and the magnetic energy per unit mass .
the calculation uses a complex form of eddy viscosity that employs a markov cutoff in `` memory '' of the eddy interactions and includes the non - linear transfer of energy among turbulent eddies of different sizes , including triads of eddies with differing wavenumbers .
some additional details are found in the appendix .
the principal advantages of this eddy damped , quasi - normal markov ( edqnm ) method are its ability to calculate the time evolution of the turbulence over very large spatial ranges and its relatively modest computational requirements , especially when compared to the more straightforward methods of direct numerical simulation .
the technique is is wide use , and its accuracy has been verified by comparison with both direct numerical simulations and with laboratory experiments ( e.g. , gomez et al .
2007 , vedula et al . 2005 ,
turner & pratt 2002 , lesieur & ossia 2000 , staquet & godeferd 1998 ) .
details of this technique can be found in orszag ( 1977 ) , lesieur ( 2008 ) , and pouquet et al .
( 1976 ) and de young ( 1980 ) for the mhd generalization . as discussed in the appendix , evolution of the turbulence
is followed by solving a series of integro - differential equations of the general form @xmath26dpdq.\ ] ] here @xmath27 is the fourier transformed magnetic or kinetic energy at wavenumber @xmath22 , @xmath28 is a dissipation constant ( either viscosity or magnetic diffusivity ) , @xmath29 is the forcing function or driving term , which in this case arises from the surface instabilities at the jet - ambient medium interface that create the large scale turbulent eddies driving the turbulent flow .
the term @xmath30 represents the non - linear transfer terms in the flow which mediate the transfer of magnetic and kinetic energies over all wavenumbers ; i.e. , @xmath30 describes the processes that create the cascade of energy from large scale to small scales , the amplification of magnetic fields by fluid turbulence , and the inverse cascades , if any , of energy from small scales to larger scales .
( which in general is significant only in the presence of flows with net helicity . )
in general all of @xmath27 , @xmath31 , and @xmath30 are time dependent . the explicit form for @xmath30 can be complex because it describes the interaction of hydrodynamic turbulent eddies with each other , the interaction of hydrodynamic turbulence with any magnetic fields and the resulting possible magnetic field amplification , the effects of eddy viscosity , and the self interaction of magnetic turbulent structures with each other .
the complete mhd forms of @xmath30 , including the effects of helicity ( * @xmath32 * and its magnetic counterpart * @xmath33 * ) are given in de young ( 1980 ) .
in the present case magnetic effects are not of major concern .
however , a very small , dynamically unimportant magnetic field is included for completeness ; its initial energy density is @xmath34 that of the turbulent kinetic energy density .
it could in principle play a role in the dissipation process at small scales , but to include these effects would require a knowledge of the detailed field reconnection and dissipation mechanisms , which in turn would require a specification of field geometries and strengths as a function of scale length .
since these are unknown , the magnetic diffusivity is chosen to be consistent with the hydrodynamic dissipation as described below . magnetic field amplification to equipartition with the hydrodynamic energies
will occur on very small scales , but in the absence of any helicity this will have no effect on the evolution of the turbulence on the scales of interest here .
the problem is specified by the injection of kinetic turbulent energy at some rate @xmath29 at a scale defined as the energy range , and the evolution of the turbulent energy spectrum is then followed as the energy cascades to smaller and smaller scales .
a natural unit for length is that corresponding to the energy injection scale , @xmath35 , and a corresponding time scale comes from the largest eddy turnover time ; i.e. , @xmath36 , where @xmath37 is the mean velocity of the turbulent flow . because the turbulence is driven by shear instabilities in the mixing layer between the agn outflow and the surrounding medium , @xmath37 is comparable to the local jet outflow speed .
the calculations here are basically scale free .
this is because the actual dissipation scale is very many orders of magnitude smaller than the energy injection scale , and inclusion of the energy scale , which is necessary , means that even with this technique the calculation can not span the scale lengths from energy injection to dissipation .
the dissipation scale is given by @xmath38 , where @xmath39 is the reynolds number appropriate to the large scale flow and @xmath25 is the wavenumber corresponding to the energy injection scale .
the hydrodynamic reynolds number is @xmath40 , where @xmath41 is the average velocity at scales @xmath25 ; i.e. , @xmath42 .
the magnetic reynolds number is just this expression with the magnetic diffusivity @xmath43 substituted for the kinetic viscosity @xmath28 . for a hot ( @xmath44k ) rarefied ( @xmath45 ) plasma @xmath46 ( e.g. , spitzer 1962 ) , and energy injection scales of tens of parsecs and sonic or subsonic flow speeds of @xmath47 km / sec imply ten or more orders of magnitude in scale between the energy injection scale and the dissipation scale . even though the equations of the form ( 2 ) above are solved numerically on a logarithmic scale in @xmath22 space , a solution that spans the complete scale from @xmath25 to @xmath24
is not feasible , especially in view of a causality condition discussed below .
thus the actual values of the dissipation constants @xmath28 and @xmath43 never enter into the calculation , and no inherent scales are present .
this means that application of the results can be applied to many different spatial and time scales , and this characteristic will be used in estimating the time for the turbulent energy flow to actually reach the regime where energy is dissipated into heating the ambient medium .    because the wavenumber range included in the calculations
is less than the complete range to @xmath24 , some mechanism has to be used at the upper end @xmath48 of the calculated wavenumber range to simulate the continuation of the loss free cascade of energy to wavenumbers beyond @xmath48 .
if this is not done , unphysical reflection of turbulent energy from the upper end of the calculated wavenumber range can occur . a standard method for treating
this is to use artificially large values of the dissipation constants @xmath28 and @xmath43 .
these are adjusted so that the resulting effective values of @xmath24 are somewhat less than @xmath48 so that no spurious energy reflections are seen .
the differential equations of the form given by eq . 2 are solved numerically on a discrete grid in @xmath22 space , and in order to avoid non - physical energy fluctuations , a time step control similar to a courant condition must be used .
the value of the time step in the integration is chosen to be less than the minimum of @xmath49 .
this ensures the correct behavior at high wavenumbers ; it also illustrates why extending the calculated wavenumber range beyond about 5 orders of magnitude is not practical .
because the calculation is scale free , the obvious scaled values to choose are @xmath50 and @xmath51 ; once the physical dimension and speed corresponding to these are chosen , the time step is determined and the problem is well defined .
the scale - free nature of the calculation allows a simple estimate of the time required to actually reach the true dissipation range once a calculation has reached equilibrium over the more limited wavenumber range described above .
suppose a calculation is made of the evolution of the turbulence spectrum between some @xmath52 and @xmath53 using a @xmath54-function energy input at @xmath52 .
let @xmath53 be some large multiple @xmath55 of @xmath52 , about 100 , and suppose the spectrum reaches a steady state at some time @xmath56 , where @xmath57 is measured in units of @xmath58 with @xmath59 ; @xmath60 is roughly the large scale eddy turnover time at @xmath52 , and @xmath61 is the mean flow speed of the turbulence on this scale .
hence @xmath56 can be measured in units of the turnover time and thus @xmath62 .
because it is scale free , the same calculation also applies to driving turbulence between @xmath53 and @xmath63 with a @xmath54-function at @xmath53 . in this case a steady state is reached in some time @xmath64 where now @xmath57 is measured in units of @xmath65 with @xmath66 .
similarly , the calculation can be made between @xmath67 and @xmath68 , and so on until some @xmath69 is reached where @xmath70 .
then the total time required for the turbulence to cascade from @xmath52 to the dissipation range is just @xmath71 . but
@xmath72 and similarly for the other @xmath73 .
but in this scale - free calculation , @xmath74 , and thus @xmath75 , and so @xmath76 .
to a first approximation all of the @xmath77 are equal , though it could be argued that as one goes to smaller and smaller scales the relevant values of @xmath77 decrease due to an increasing effective viscosity , especially in the presence of magnetic fields . by setting all @xmath78 a lower limit on the value of @xmath79
is obtained , and this is the most constraining value for the issue of concern here .
then @xmath80 , and since @xmath81 the time for the turbulent flow to reach the dissipation range and begin to heat the ambient medium is @xmath82 .
% $ $ \ ] ] thus a good approximation of the time for the turbulent cascade to reach the dissipation stage is just the number @xmath83 of large scale eddy turnover times that are required for the establishment of a steady state turbulent cascade over a factor of about @xmath84 in scale and beginning with the scale of the large scale energy injection .
i.e. , @xmath85 .      in this model
the energy injection occurs on the scale of the largest eddies when the outflowing jet first becomes fully turbulent from outer boundary to centerline . in the limiting case
the energy injection spectrum could be taken as a delta function centered on @xmath86 , and the energy injection function in eq .
2 would have the form @xmath87 ) . since the goal of the calculation is to determine the time required for the turbulence originating on a scale @xmath25 to cascade down to the dissipation range
, the injection spectrum should be some sharply peaked function centered on @xmath25
. a narrow gaussian form may be more cosmetically appealing but will be seen to have little effect on the outcome of the calculation .
the value of the constant @xmath88 controls the amplitude of the energy injection spectrum , and because the calculation is looking at the minimum time required for the turbulent cascade to reach the dissipation range , values of @xmath88 near the maximum allowed are of interest .
the units of @xmath88 are ( kinetic ) energy per unit mass per unit time , and in the dimensionless units used here , @xmath88 should be some fraction of @xmath89 , where @xmath90 is the value of the jet radius at which the flow becomes fully turbulent throughout its volume .
the maximum value of @xmath88 is then equal to one .
a `` reasonable '' choice might be to set an average large scale eddy size to half @xmath90 and the similarly for some average rotation speed of those eddies as @xmath91 , which gives a value of @xmath92 in dimensionless units where @xmath93 and @xmath94 .
in fact , it will be seen that values of @xmath88 between 1.0 and 0.1 yield roughly the same results . introducing a time variation in @xmath88
is easily done , but such variations could only be less than the maximal value and thus defeat the purpose of finding minimum equilibration times . hence @xmath88 will be kept constant .
a major virtue of the formalism used here is that the calculation is carried out with logarithmic intervals in @xmath22 space , and though this introduces some complications to ensure that energy transfer among structures in neighboring regions of wavenumber space is not incomplete , it allows the calculation to cover many orders of magnitude in scale and thus permits the development of a turbulent cascade through the loss - free inertial range . in general
the calculations here cover almost four orders of magnitude in scale in three dimensions , which would require a major computational effort in a direct numerical simulation .
most calculations used a wavenumber range of @xmath95 $ ] , while more extensive runs extended the upper limit to @xmath96 .
the kinetic energy in the edqnm formalism must be positive definite , so the initial energy spectra were set equal to a constant but negligible value of @xmath97 in scaled units at all wavenumbers .
this is three to four orders of magnitude below the injection energy and has no effect on the evolution of the energy spectrum resulting from the turbulent mixing layer .
the energy levels at the upper and lower ends of the wavenumber range were not fixed after initialization . in practice
there is no significant inverse cascade of energy in the absence of helicity , so the value of @xmath98 at the low end of the wavenumber range ( largest scale ) does not evolve from its initial value .
the energy at the largest wavenumber @xmath48 begins to fall immediately due to the very rapid energy transfer on the smallest scales , but this value always remains non - zero .
the major free parameters in these calculations are the shape and amplitude of the energy injection spectrum @xmath99 .
other parameters that can be varied are the wavenumber range , the length of time covered by the calculation , the value of the viscosity to test for any influence on the spectral evolution , and the density of points in wavenumber space . for most of the calculations
the shape of @xmath99 was taken to be a delta function centered at the energy injection scale @xmath25 , as discussed above . in
some runs this shape was changed to a gaussian distribution in the log of injection energy sharply peaked at @xmath100 , and with @xmath101 in log@xmath22 .
the values of @xmath88 ranged from the maximum of 1.0 to 0.01 , with the most common value being 0.25 as described previously .
the normalization constant for the gaussian runs was chosen so that the peak of the distribution had the same amplitude at @xmath102 as did the delta function energy injection , thus providing a comparable energy injection rate .
the scaled wavenumber range was in most cases from 0.015 to 128 , with a few runs being made over the range @xmath103 $ ] .
values of the time evolution of the turbulent energy spectra were calculated at equally spaced logarithmic intervals in the @xmath22 domain ; some runs were repeated with the density of points increased by a factor of two to see what effect , if any , would result from changing the size of @xmath104 . as described above , values of the viscosity @xmath28 were chosen so that an effective @xmath24 lay near the value of @xmath48 .
various small changes in @xmath28 were made in some runs to find the smallest value of @xmath28 that would not produce any `` feedback '' from the upper end of the @xmath22 range into the larger scales .
such unphysical inverse cascades were immediately apparent by the growth or `` bunching up '' of elevated @xmath98 values at the top end of the wavenumber range .
adequately high dissipation to ensure only a direct and unimpeded cascade to smaller scales ( large @xmath22 ) was always shown by a rolloff in values of @xmath98 from a power law at the highest wavenumbers .
the timestep was determined for each run via the effective courant condition described earlier ; increases in @xmath105 much above this value immediately result in zero or negative values of @xmath98 at the highest wavenumbers .
the objective of these calculations is a determination of the time elapsed from the onset of driven turbulence on some large scale to the transformation of the turbulent energy into heat at the dissipation scale .
this first occurs when the turbulent cascade reaches an equilibrium or steady state condition with energy flowing in at @xmath25 and being dissipated at @xmath24 after transit through a loss free cascade to ever smaller scales .
establishment of this equilibrium is an asymptotic process , hence the time required for equilibrium will be set at that time when any changes in the energy spectrum per unit time are less than a few percent .
this will give a minimum time for establishment of equilibrium and will result in the minimum time for feedback into heat to occur .
the approach to equilibrium is slightly complicated because there is no well defined `` front '' of energy marching down to smaller scales into a constant background of very low turbulent energy .
instead , as soon as the calculation is begun , energy begins to flow to higher wavenumbers , as it should .
this means that the initial constant and very small value of @xmath98 set at the beginning shows a decay at large @xmath22 before the energy flow from the injection region arrives at the smallest scales .
hence the shape of the spectrum evolves over the    full wavenumber range for all times greater than zero .
figure 1 , which shows the early evolution of the turbulent kinetic energy spectrum as a function of time , illustrates this behavior . in this case
the energy injection spectrum has a delta function form for @xmath29 at @xmath106 , and time is measured in units of the @xmath106 eddy turnover time , @xmath107 . at early times ( @xmath108 )
the effects of the energy drain at large wavenumber is very evident , since the spectrum in this range is reduced far below the initial constant value of @xmath97 described in section 3.2 .
nonetheless , at @xmath109 the flow of energy from the injection region at @xmath106 to smaller scales is clearly seen , with a `` front '' of the energy cascade having arrived as far as @xmath110 . by @xmath111
the effects of energy injection have reached to scales of order @xmath34 smaller than the injection region , but , as figure 1 shows , the spectrum has not yet reached an equilibrium state by @xmath112 .
figure 2 shows the late time evolution of the kinetic energy spectrum , and there it is clear that for scales smaller than the injection scale , there is very little evolution in the spectrum at times later than @xmath113 , with little or no evolution past @xmath114 . at larger scales
there continues to be a slow growth of turbulent energy on scales of up to 10 times the size of the injection region , consistent with previous calculations .
however , this region is of no concern for the present problem , and in addition the applicability of homogeneity on these scales is less certain for the jet structures considered here .
the region @xmath115 has no effect on small scale evolution , and in the absence of helicity almost all the injected energy cascades to the small scale regions ( e.g. , de young 1980 ) .
figure 3 shows the resultant spectrum at @xmath113 when the wavenumber space is extended to @xmath116 , thus spanning over 5 orders of magnitude in scale . as discussed above , this is a practical maximum for a single calculation due to the time step constraints at the smallest scales .
figure 3 also shows a comparison between a @xmath54 function energy injection at @xmath106 and a gaussian energy injection on the same scale with the normalization discussed in the previous section .
it is clear that there is no significant difference in the late time spectra that result from the two different energy injection functions .
also evident in figure 3 is the establishment of a clear power - law energy spectrum for times greater than    @xmath113 .
finally , figure 4 shows the time evolution of the total ( kinetic plus magnetic ) energy for the @xmath54 function injection case that covers the range @xmath117 . from this figure
it is clear that the turbulence has not reached an equilibrium state for times less than @xmath118 .
however , it does seem that equilibrium is established by @xmath114 .
the early rapid growth of @xmath119 for times less than ten is due to the constant inflow of injected energy with no loss at the high wavenumber end having yet occurred . the slower growth for @xmath120 reflects the evolution of the spectrum after energy loss begins but is not yet in equilibrium .
in addition , after @xmath111 the conversion of kinetic energy into magnetic energy begins to become noticeable at the few percent level . throughout all these calculations ,
the total turbulent energy is always dominated by the kinetic energy ; the maximum magnetic energy integrated over all scales never exceeds about 70 percent of the kinetic energy .
additional calculations were made with variations of other parameters mentioned in the previous section , such as variation of @xmath121 , @xmath104 , @xmath105 , @xmath122 , and @xmath28 .
however , these variations did not significantly alter the overall results presented in figures 1 though 4 .
in particular , the variation of @xmath121 , which might be thought to have the largest effect , was not particularly notable .
calculations were made with values of @xmath121 varying from 0.1 to 1.0 , and while the amplitude of the turbulent energy spectrum increased as @xmath121 was increased , the shapes of the equilibrium spectra were virtually invariant over an order of magnitude change in @xmath121 .
this is not surprising given the tendency for the turbulent mhd cascade to evolve toward a common power - law form .
moreover , the scaled time required to establish equilibrium was about @xmath123 in all cases , which is consistent with a scale free calculation .
variations in @xmath105 and @xmath28 have been discussed above ; these are controlled by the courant condition to preserve positive definite energies and to avoid unphysical inverse cascades of energy at large wavenumber .
the effects of changing the upper limit on the wavenumber range are shown in the figures ; a more clear definition of the loss - free power - law energy cascade is seen with larger @xmath122 , but the time required to reach equilibrium is not affected .
this is also an expected result , since equilibrium is reached most rapidly at the smallest scales ( largest @xmath22 ) .
decreasing the value of @xmath104 and increasing the density of computational points has the expected effect of making the energy transfer among triads of wavenumbers slightly more efficient . a doubling of the density of points in @xmath22 results in maximum change in @xmath124 of 7 percent at the upper end of the wavenumber range .
again , the time to reach equilibrium , which is governed by processes on the larger scales , is unaffected .
thus a robust conclusion that emerges from these calculations and which is illustrated in figures 2 and 4 is that the turbulent cascade reaches equilibrium over the wavenumber range @xmath125 in about 100 large scale eddy turnover times ; @xmath126 .
this result can then be combined with equation 3 for the time required after the onset of energy injection to reach turbulent dissipation of energy into heat ; @xmath127 \\   \approx ( n\ell_{1}/v).\ ] ] thus these results imply that a value of @xmath128 is appropriate for these turbulent jet flows . with this result in hand
it is now possible to estimate the time required for radio agn outflows to convert their outflow energy into heat . to do this
requires a conversion from the scale free calculation here into the relevant astrophysical quantities .
the first scale to establish is the spatial scale that corresponds to @xmath25 where the energy is injected .
this is basically the size of the large scale fully developed turbulent eddies which are driving the turbulence , powered by the vorticity generated in the shearing motions of the mixing layer .
as described earlier , the mixing layer grows in thickness until it permeates the entire jet , at which point the flow is fully turbulent , and it is at stage that the interior of the jet can be characterized as a homogeneous and isotropic ( in a comoving frame ) turbulent flow .
hence the size of the large scale eddies at the point of complete turbulent infiltration or saturation is an appropriate measure of the energy injection scale .
given the internal structure of mixing layers , a good approximation for this eddy size is just the radius of the outflowing jet at saturation , since that is where the linearly expanding mixing layers merge together and the potential flow in the core of the jet ceases . experimental data on gaseous turbulent jets also show that the large scale structures at this point are comparable to the jet radius ( dimotakis et al . 1983 ) .
in addition , it is at this point that the opening angle of the outflow increases slightly and takes on the common value of @xmath129 seen in fully turbulent round jets .
( some experimental data for turbulent jets also show the creation of larger scale structures that develop in the downstream flow , but these are often attributed to merging or pairing of large eddies formed earlier rather than an indication of an energy injection scale ( e.g. , brown & roshko 1974 , dimotakis et al .
in addition , these structures are most apparent in fluid turbulent jets and are less dominant in very high reynolds number gaseous jets ( dimotakis et al .
1983 , dahm & dimotakis 1990 ) . ) determination of the jet radius when it becomes fully turbulent can be estimated from analytic approximations to @xmath130 or from observations of the early stages of fr - i outflows .
the latter is probably the more robust method , given the uncertainty in the initial jet densities and outflow speeds .
radio observations of the inner regions of 3c 31 and 3c 296 ( hardcastle et al . 2002 , laing et al . 2008 , hardcastle et al . 2005 ) permit an estimation of the jet width in its earliest turbulent phase .
in particular , the southern jets in these objects show an increase of the jet opening angle to @xmath131 degrees that marks the transition to self - similar fully turbulent outflow .
estimates of the jet radii at these points give values in the the range @xmath132 kpc , and thus a reasonable estimate from these data for the scale corresponding to @xmath133 is about @xmath134 kpc .
a consistency check can be made by looking at the corresponding values of @xmath130 and the density ratios and outflow speeds that would be required .
for example , in 3c 31 the estimated point where the jets appear fully turbulent is at a distance of about 1.7 kpc from the nucleus , while in the case of 3c 296 this distance is about @xmath135 kpc . using this as a rough measure of @xmath130 , and assuming that the functional form for plane mixing layers given above can be generalized to mixing layers for round jets , given the overall similarity of the mixing layers in the two cases ( e.g. , freund et al .
2000 , dimotakis et al .
1983 ) , then it is possible to see if `` reasonable '' values of outflow densities and speeds near the nucleus result .
the jet becomes fully turbulent at a radius @xmath136 , and using @xmath137 and @xmath138 kpc as an average of the above numbers , together with the observed temperatures and inferred number densities for the ambient medium in 3c 31 and 3c 296 at about 1.5 kpc from the their nuclei ( hardcastle et al .
2002 , 2005 ) gives jet flow speeds in the innermost regions of @xmath139 cm s@xmath7 and jet number densities of @xmath140 @xmath141 , this is consistent with models of `` light '' jets which are mildly relativistic in the innermost regions ( e.g. , laing et al .
thus different and independent approaches indicate that a value of @xmath142 kpc is a consistent value to take for the scale of the energy injection range , and this value will be taken to be the size corresponding to @xmath25 .
once a value of @xmath90 is established , the remaining scales can be fixed by setting a velocity scale , and this is usually chosen to be the turbulent turnover speed of the large scale eddies in the energy range .
the flow in the mixing layer is subsonic relative to the ambient medium as argued previously , perhaps transsonic at most in the initial stages .
hence the value of the eddy turnover speed @xmath143 will be some fraction of the local sound speed in the ambient medium , and this in turn has been established for some fr - i objects via observations of x - ray emission in the central regions of the parent galaxies .
these observations are consistent with thermal bremsstrahlung radiation , and thus a temperature and sound speed for this gas can be established . in particular ,
chandra observations of 3c 31 and 3c 296 ( hardcastle et al . 2002 , 2005 ) are consistent with gas temperatures of order @xmath144k , with corresponding sound speeds of @xmath145 cm s@xmath7 .
thus a an average value of @xmath146 in the energy injection range is consistent with observations and will be used here .
this value is significantly less than the above estimate of @xmath147 in the innermost region of the jet , and so it is also consistent with significant deceleration having occurred , especially in the outer layers of the jet , due to mass entrainment via the growing turbulent mixing layer .
once these scales are established , a time unit of @xmath148 , the large eddy turnover time , follows naturally , as does the unit of turbulent kinetic energy density per unit mass , @xmath149 .
this last quantity is the unit ( per unit wavenumber ) used in the energy spectra shown in figs . 1 - 3 .
the previous two sections provide the information needed to establish an estimate of the time required from the establishment of a fully turbulent outflow in an fr - i jet to the onset of conversion of turbulent energy into heat . as argued above
, both observations and analytic approximations suggest a large scale eddy size at the energy injection region of @xmath142 kpc , with an accompanying turbulent eddy speed of @xmath150 cm s@xmath7 .
this gives a time unit for the turbulence calculations of @xmath151s or @xmath152 yr .
the nonlinear turbulence calculation shows the establishment of an equilibrium spectrum in about 100 large scale eddy turnover times , which gives an estimate of the time for the onset of conversion of outflow energy into heat of about 100 million years .
these times are comparable to or less than the bremsstrahlung cooling times deduced from x - ray observations of hot gas in and around the host galaxies of fr - i radio sources ( e.g. , hardcastle et al . 2002 , 2005 ; laing et al .
2008 ) .
where is this heat deposited ?
the above time estimate is a measure of the beginning of the conversion into heat ; it presumably continues for a time comparable to the lifetime of the agn outflow .
this injection of heat can not take place within the confines of the outflows that are observed to be well defined jets , since the energy density in these regions is so high that conversion of all the directed kinetic energy into heat at that point would result in pressures high enough to decollimate the jet .
hence there is a need to know how far the jets will propagate in the @xmath0 years after they become fully turbulent .
the mean outflow speeds in these objects are unknown , but an estimate can be gained from the arguments that the jets are fully turbulent flows . as mentioned before ,
such jets are self - similar flows , and experimental data show that the on - axis outflow speed in such jets is approximately 10 times the mean outflow speed on the surface of the jet ( e.g. , pope 2000 , hussein et al .
1994 ) . from previous considerations
the outflow speed on the surface of the jet is of the order of or less than the local sound speed , and this has been estimated to be @xmath44 cm s@xmath7 in the region of turbulent energy injection .
thus in this phase of the jet propagation a mean outflow speed of @xmath0 cm s@xmath7 may be used , which would place the onset of conversion into heat at a distance of 100 kpc or so from the onset of the fully turbulent phase .
this is well beyond the collimated outflow region of most fr - i objects , and hence deposition of heat in this region would not be in conflict with the observed properties of the jets .
moreover , this region may lie in that portion of the outflow where the agn flow has become much more broadly dispersed and less well collimated .
this decollimation process is necessary if radio agn feedback is to work , since the heating must take place over a large volume in order to provide significant suppression of star formation .
the 100 kpc estimate is somewhat uncertain , since beyond ten kpc or so from the nucleus the radio morphology of fr - i outflows shows a clear departure from simple self - similar flow .
however , it could be that there is very little deceleration beyond this point , and that most of the mass entrainment is completed by this time .
the jet morphologies are not inconsistent with this , in which case the outflow speeds near the centerline of the jet remain about the same as they were in the self - similar regime , and the above estimate of @xmath153 kpc retains some validity .
this is suggestive , but hardly definitive .
there is in addition the fact that centerline outflow speeds in self - similar jets decrease linearly with distance , though the relevant constant in this relation is unknown ; hence there will be some , possibly non - negligible , deceleration of the jet from the onset of self - similarity until the @xmath154 kpc scale is reached .
the late stages of the evolution of these flows raises the issue of their relation to the agn outflows seen in the cores of rich galaxy clusters .
much work has been done on the role of `` radio bubbles '' in the icm , relic and otherwise , and their possible influence on the evolution of the icm and brightest cluster galaxies ( e.g. , fabian et al . 2002 , brggen & kaiser 2001 , fabian et al .
2002 , de young 2003 , jones & de young 2005 , oneill et al .
2009 ) , and a natural question arises about the relation of these objects to the flows discussed here . in general , most fr - i agn do not show radio `` bubbles '' , nor are they usually found in the cores of rich clusters .
for those that are , it it likely that relatively low powered fr - i jets would be rapidly decelerated and confined by the high ambient pressures and could slowly inflate buoyant cavities in the hot icm , as has been suggested by many models .
such cavities are usually small for low radio powers , a few tens of kpc , and the jet structures feeding them are even smaller ( e.g. , mcnamara & nulsen 2007 ) .
by contrast , the fully turbulent fr - i jets discussed here are themselves tens of kpc in length before the onset of complex large scale flow patterns such as bends or meanders .
hence there is likely to be a continuum in radio source morphology for a given radio power as one moves from very high pressure regions in the centers of rich clusters to lower pressure environments found in small clusters , groups , and individual galaxies .
the onset of mixing layers , entrainment and deceleration could occur in all these environments , as has been argued here , with the process occurring on smaller scales in more dense and high pressure environments ( cf .
however , one important added factor in the dense icm in cluster cores is the presence of non - negligible external magnetic fields .
although the kh instability and resultant mixing occurs unabated in mhd flows with weak fields ( ryu , jones & frank 2000 ) , if the average magnetic energy density becomes comparable to the kinetic energy density , then the mixing layer picture may have to be modified , especially if the external fields have long coherence lengths .
in particular , expansion of an inflated `` bubble '' into a magnetized medium may not result in the late stage mixing that simple hydrodynamic flows would suggest ( de young 2003 , de young et al .
2008 ) , so in the case of cluster cores the final step in agn feedback may be more difficult to accomplish .
this issue is at present still unresolved .    though the creation of bubbles and cavities may be the end state of the flow in high pressure environments , in lower pressure regions the final state could be much more diffuse and distorted due to large scale external magnetic fields or to any flows in the circumgalactic gas caused by galaxy orbits or merging events .
the detection ( or not ) of such faint extended regions is very much a function of x - ray and radio telescope sensitivity and selection effects . in this context
is it very interesting to note that rudnick & brown ( 2009 ) have recently observed very large , diffuse polarized radio structures surrounding the fr - i object 3c 31 . though further observations are needed , such very large structures could be the signature of the decollimated and widely dispersed portions of radio agn outflows where the heat is deposited from fully turbulent jets .
a similar widely dispersed radio emitting region has been known for some time to surround the collimated outflows from the fr - i source associated with m87 ( owen et al .
thus the issue of where the turbulent energy is actually deposited remains an important unresolved problem .
it was argued above that a distance from the nucleus of @xmath153 kpc was not unreasonable for the onset of heat deposition , but also discussed there were the uncertainties associated with agn outflow speeds for distances greater than @xmath154 kpc . for isolated galaxies ,
heat deposition at @xmath153 kpc is clearly too distant for thermal suppression of star formation in the near galactic halo , yet the presence of coherent flows on scales @xmath155 10 kpc argues for conversion to heat beyond 10 kpc .
the picture is further complicated by low surface brightness radio data for some fr - i objects ( e.g. , 3c 296 , m87 ) ; these show the presence of faint radio emission very close to the stellar extent of the galaxy , and the implication here is that agn outflow material from earlier epochs may have been advected or diffused back into the near circumgalactic regions .
these data could imply that deposition of heat is possible into inner galactic halo at distances a few tens of kpc from the nucleus , probably on timescales in excess of @xmath156 years .
however , it seems clear that if energy deposition is needed for isolated galaxies in the inner regions of the ism , further investigation into the very late time evolution of outflows is needed to see if this process is viable .
the object of this paper is an exploration of the feasibility of radio agn feedback .
the constraints on this process are that it must heat the interstellar and circumgalactic medium around the host galaxy in a time less than or the order of the ambient gas cooling time and with enough volume coverage so that significant star formation can be suppressed and the host galaxy will appear to be undergoing passive evolution .
the issue of radio agn feedback in massive , rich clusters of galaxies is not specifically addressed , as this process has been treated in many other papers . instead
, the more general case of radio agn feedback for galaxies in small groups and clusters and even in isolated systems is considered here .
the emphasis is to proceed from the basic properties of agn outflows to see if a continuous chain of physical processes can be established that will lead to the required heating in the required time .
the analysis has a small number of fairly simple steps .
first , observations show that the most common class by far of radio agn in universe out to @xmath157 and perhaps beyond is the class of fr - i radio sources .
second , the morphology of these objects strongly suggests that their bipolar outflows are subsonic or transsonic beyond a few kpc from the nucleus .
third , given the evidence for relativistic outflows very close to the nucleus together with the nearly universal onset of global surface instabilities in such shearing flows , this slow outward motion strongly suggests that it is due to deceleration of the flow via the formation of non - linear mixing layers and their accompanying mass entrainment on the surface of the jet .
fourth , experimental evidence argues convincingly that such turbulent mixing layers will soon penetrate throughout the jet volume , resulting in a fully turbulent flow in the interior of the bipolar outflows .
this conclusion is reinforced by the opening angle of the radio jets observed in some nearby fr - i objects ; this angle is completely consistent with the experimental data on fully turbulent self - similar jet outflows .
fifth , it is well known that such turbulent flows eventually cascade to smaller and smaller structures until a dissipation range is reached where most of the turbulent kinetic energy is converted into heat . a series of 3d time dependent calculations of the evolution of such turbulent flows is then performed to see if this eventual conversion of collimated outflow energy into heat can take place in a time short enough to be of interest for agn feedback .
the answer appears to be yes , with the first onset of heating occurring roughly @xmath156 years after the development of fully turbulent flow ; this time is comparable to or less than the cooling times derived from x - ray observations of gas surrounding nearby fr - i radio galaxies . finally , the question of where this heat is deposited and whether or not the deposition is sufficiently isotropic is briefly addressed . to retain the observed radio source structure ,
heating probably occurs beyond the extent of the coherent jet outflow .
estimates of outflow speeds for fully turbulent jets indicate energy deposition on scales as large as @xmath9 100 kpc , which is well beyond the coherent jet structures for most fr - i objects .
isolated galaxies may require heating on smaller scales , and there is some evidence this might be occurring due to late time evolution of the flows .
recent observations of very large scale and diffuse radio structures around 3c 31 and earlier observations of m87 suggest that the suggested large scale heat deposition may be taking place .    these simple arguments , while suggestive , are not conclusive .
it seems likely that the idea of slow and turbulent outflows from fr - i radio agn is either almost certainly right or very wrong , the latter because we have missed the fundamental nature of jet outflows from agn .
for example , if these flows are `` poynting flows '' that are completely dominated at all times by magnetic fields and have very few particles ( e.g. , li et al .
2006 ) , then the nature of the interaction with the surrounding medium might be very different , with little or no role for mhd or hydrodynamic turbulence .
however , it is not clear how such magnetically dominated outflows can reproduce the morphology observed in the fr - i radio sources
. even if turbulent jets are the correct description , there are details of their evolution suggested by observations that may mean departures from the isotropic turbulence in their interiors that is used here .
deflections , bending , or propagation into decreasing density gradients or into inhomogeneities may all cause departures from this condition and perhaps changes in the time required for turbulence to decay into heat .
the magnitude of such departures is currently unknown , and they may be small enough that the above treatment remains a good approximation , especially in the interior of the jet . finally , the large scale distribution of energy from radio agn outflows remains a largely unsolved problem .
the observational hints are suggestive , but more observations of faint and diffuse radio emission around fr - i outflows are needed , as are additional calculations of the very late stages of these outflows .
portions of this work were completed during the aspen center for physics astrophysics programs in 2008 and 2009 .
thanks are due to bob rosner and to tom jones for illuminating discussions and to jessica moy for bibliographic assistance .
presented here are some additional concepts and relations that lie behind the derivation of equation ( 2 ) .
most of the basic framework of the edqnm method is given in the main body of the paper , and complete derivations for all aspects of the method , together with more detailed descriptions of the method and its variants , are provided in the references cited .
a particularly complete treatment is found in lesieur ( 2008 ) .
where @xmath161 is the magnetic field , @xmath162 is the density , @xmath163 is the fluid velocity , @xmath164 is the pressure and @xmath165 is the magnetic susceptibility . here
@xmath28 is the kinematic viscosity , @xmath43 is the magnetic diffusivity , and @xmath166 is a forcing term that represents the injection of kinetic energy into the system . in principal
the formalism also accommodates the injection of magnetic energy as well as kinetic and magnetic helicities ( @xmath167 and @xmath168 ) .
the treatment described here is restricted in nearly incompressible flows that are almost isotropic in that they admit helicity in the flow fields .
extensions of the method to anisotropic and quasi compressible flows are described in the references .    as mentioned in the text ,
the edqnm method derives equations that govern the second order moments of velocity and magnetic field which thus provide the evolution of the kinetic and magnetic energy densities .
these governing equations actually determine the behavior of the fourier transformed moments ; because fourier components are global averages over configuration space , the equations provide the evolution of energy on a particular scale over the entire flow and thus provide immediately the time evolution of the energy spectra .
in addition , operation in fourier space allows in this formalism the integration of the governing equations using logarithmic intervals in wavenumber space and so permits coverage over many orders of magnitude in spatial scale in three dimensions , a task that would be extremely difficult through direct numerical simulation .
the general form of equation ( 2 ) in the text can be obtained by first taking the fourier transform of the above navier - stokes equation , which gives , for no magnetic fields and no helicity , with @xmath162 explicit , the general form @xmath169\ ] ] where @xmath170 is the fourier transform of the velocity @xmath171 , @xmath31 is the fourier transform of the forcing ( energy injection ) term and f.t . is the fourier transform operator .
the transformed incompressibility condition becomes @xmath172 , and this condition allows elimination of the pressure from the equation of motion .
this is because the incompressibility equation shows that @xmath173 is in a plane perpendicular to @xmath174 , and thus so are @xmath175 and @xmath176 , while the pressure gradient @xmath177 is parallel to @xmath174 . as a result , the f.t .
of @xmath178 is the projection onto the plane containing @xmath173 of the f.t . of the non - linear term @xmath179 .
thus the form of the fluid equation in fourier space becomes , in general , @xmath180 where the vector notation has been dropped .
the exact form of this equation can be found in lesieur ( 2008 ) .
it can now be seen that a similar equation can be found for the fourier transform of the second moment of eq .
a1 , and it is also clear that the result of this will include a third order term in @xmath170 on the right hand side . similarly an equation for the third order moments will include an integral over a term of fourth order moments ; this difficulty is of course a result of the non - linear nature of the flow arising from the @xmath181 term in the momentum equation .
closure of this system of moments is accomplished by assuming that the fourth order moments are nearly gaussian ; hence the `` quasi - normal '' approximation .
this approximation allows the fourth order moments to be written as sums of products of the second order moments , and because the third order moments can already be written as a function of second order moments , the system is closed at this level .
a discussion of the accuracy of this approximation and tests of its rigor are found in lesieur ( 2008 ) and references therein .
thus the origin of the overall form of eq ( 2 ) in the text can be seen . in general
the integral in this equation will contain many terms that are quadratic in kinetic energy , helicity , and magnetic energy , and numerical solution of this integro - differential equation is necessary .
it can be shown ( lesieur 2008 ) that energy is strictly conserved when integration over wavenumber space includes all triads of turbulent eddies whose wavenumbers ( k , p , q ) satisfy the requirement that @xmath182 ; thus interactions with @xmath183 include non - local interactions in wavenumber space .
a complete tabulation of all the terms included on the right side of eq ( 1 ) is found in de young ( 1980 ) .
in addition to a closure approximation , two other constraints are added to the right side of eq .
( 2 ) . in its simplest form ,
the quasi - normal approximation can , at late times , lead to negative energies in some portions of the spectra .
this has been shown to be due to an eventual excess value building up in the third order moments .
these moments have been shown experimentally to saturate , and a modification to the theory incorporates this effect by adding a linear damping term in the equation governing the growth of the third order moments which represents the deformation rate of eddies of size @xmath184 by larger scale eddies .
this `` eddy damping '' stabilizes the growth rate of the third order moments and results in the eddy damped quasi - normal ( edqn ) version of the theory .
a final modification guarantees the positive definite character of the spectra in all situations , not just at late times .
this is the introduction of a markov process to remove excessive `` memory '' in the system on timescales of order of the large eddy turnover time , including removal of any time variations in the eddy damping term . with this step
the formalism takes on the final form of the edqnm theory .
specific forms for these modifications can be found in de young ( 1980 ) and lesieur ( 2008 ) .
aloy , m.a . ,
ibaez , j.m . ,
mart , j - m , gmez , j - l . ,
& mller , e. 1999 , apj , 523 , l125 benson , a.j .
, bower , r.g . ,
frenk , c.s . ,
lacey , c.g . ,
baugh , c.m .
, & cole , s. 2003 , apj , 599 , 38 best , p.n . , kaiser , c.r . ,
heckman , t.m .
, & kauffmann , g. 2006 , mnras , 368 , l67 bicknell , g.v .
1984 , apj , 286 , 68 brzan , l. , mcnamara , b.r . , nulsen , p. , carilli , c. , & wise , m.w .
2008 , apj , 686 , 859 bower , r.g .
2006 , mnras , 370 , 645 brown , g.l . , & roshko , a. 1974 , j.fl.mech , 64 , 775 brggen , m. , & kaiser , c.r 2001 , mnras , 325 , 676 brggen , m. , & kaiser , c.r 2002 , nature , 418 , 301 brggen , m. , heinz , s. , roediger , e. , ruszkowski , m. , & simionescu , a. 2007 , mnras , 380 , l67 cara , m. , & lister , m.l .
2008 , apj , 674 , 111 chandrasekhar , s. 1961 , hydrodynamic and hydromagnetic stability ( oxford : clarendon ) , chap .
11 clemens , n.t . , & mungal , m.g .
1995 , j.fl.mech , 284 , 171 cowie , l.l . ,
songaila , a. , hu , e.m . , & cohen , j.g .
1996 , aj , 112 , 839 croton , d.j .
2006 , mnras , 365 , 11 croton , d.j . , &
farrar , g. 2008 , mnras , 386 , 2285 dahm , w. , & dimotakis , p.e . 1990 , j .. fl.mech , 217 , 299 de young , d.s . 1980 ,
apj , 241 , 81 de young , d.s .
1981 , nature , 293 , 43 de young , d.s .
1986 , apj , 307 , 62 de young , d.s .
1993 , apj , 405 , l13 de young , d.s .
2002 , the physics of extragalactic radio sources ( chicago : univ .
chicago press ) de young , d.s .
2003 , mnras , 343 , 719 de young , d.s .
2006 , apj , 648 , 200 de young , d.s . ,
oneill , s.m . , & jones , t.w .
2008 , aspc , 388 , 343 dimotakis , p.e . , &
brown , g.l .
1976 , j.fl.mech , 78 , 535 dimotakis , p.e .
, miake - lye , r.c . , &
papantoniou , d.a .
1983 , phys.fl . , 26 , 3185 donoso , e. , best , p.n . , &
kauffmann , g. 2009 , mnras , 392 , 617 fabian , a.c .
2000 , mnras , 318 , l65 fabian , a.c .
, celotti , a. , blundell , k.m . ,
kassim , n.e .
, & perley , r.a .
2002 , mnras , 331 , 369 fanaroff , b.l . ,
& riley , j.m .
1974 , mnras , 167 , 31p ferrari , a. 1998 , araa , 36 , 539 freund , j.b . ,
lele , s.k . , & moin , p. 2000 , j.fl.mech , 421 , 229 gomez , t. , sagaut , p. , schilling , o. , & zhou , y. 2007 , phys .
19 , 48101 graham , j. , fabian , a.c . , & sanders , j.s .
2008 , mnras , 386 , 278 hardcastle , m.j . ,
alexander , p. , pooley , g.g .
, & riley , j.m .
1997 , mnras , 288 , l1 hardcastle , m.j . ,
worrall , d.m .
, birkinshaw , m. , laing , r.a . , & bridle , a.h .
2002 , mnras , 334 , 182 hardcastle , m.j . ,
worrall , d.m .
, birkinshaw , m. , laing , r.a . , & bridle , a.h .
2005 , mnras , 358 , 843 hussein , h.j . , capp , s. , & george w.k .
1994 , j.fl.mech 258 , 31 jackson , c.a .
, & wall , j.v .
1999 , mnras , 304 , 160 jones , t.w . , & de young , d.s .
2005 , apj , 624 , 586 laing , r.a . ,
parma , p. , de riter , h.r .
, & fanti , r. 1999 , mnras , 306 , 513 laing , r.a . , & bridle , a.h .
2002 , mnras , 336 , 328 laing , r.a . , & bridle , a.h .
2004 , mnras , 348 , 1459 laing , r.a . ,
bridle , a.h . ,
parma , p. , feretti , l. , giovannini , g. , murgia , m. , & perley , r.a .
2008 , mnras , 386 , 657 ledlow , m.j . , &
owen , f.n .
1996 , aj , 112 , 9 lesieur , m. , & ossia , s. 2000 , j. turb . , 1 ,
7 lesieur , m. 2008 , turbulence in fluids , ( springer ) , 4th ed .
li , h. , lapenta , g. , finn , j.m . , li , s. , & colgate , s.a .
2006 , apj .
643 , 92 madau , p. , ferguson , h.c . ,
dickinson , m.e . ,
giavalisco , m. , steidel , c.c . , & fruchter , a 1996 , mnras , 283 , 1388 mandelbaum , r. , mcdonald , p. , seljak , u. , & cen , r. 2003 , mnras , 344 , 776 mcnamara , b.r . , & nulsen , p.e.j .
2007 , araa , 45 , 117 mungal , m.g . , lozano , a. , & van cruyningen 1992 , exp.fl , 12 , 141 okamoto , t. , nemmen , r. , rodrigo , s. , & bower , r.g .
2008 , mnras , 385 , 161 oneill , s.m .
, de young , d.s . , &
jones , t.w .
2009 , apj , 694 , 1317 orszag , s.a . 1970 ,
j.fl.mech , 41 , 363 orszag , s.a . , & patterson , g.s .
1972 , phys.rev.lett , 28 , 76 owen , f.n .
, & ledlow , m.j .
1997 , apjs , 108 , 41 owen , f.n . ,
eilek , j.a . , & kassim , n.e .
2000 , apj , 543 , 611 perucho , m. , hanasz , m. , mart , j - m .
, & miralles , j - a .
2007 , phys.rev.e , 75 , 56312 pope , s.b .
2000 , turbulent flows ( cambridge : cambridge univ press ) pouquet , a. , frisch , u. , & lorat , j. 1976 , j.fl.mech 77 , 321 pouquet , a. , & patterson , g.s .
1978 , j.fl.mech , 85 , 305 rigby , e.e . , best , p.n . , & snellen , i. 2008 , mnras , 385 , 310 ruszkowski , m. , brggen , m. , & begelman , m.c .
2004 , apj , 615 , 675 rudnick , l. , & brown , s. 2009 , apj , 137 , 145 ryu , d. , jones , t.w .
, & frank , a. 2000 , apj , 545 , 475 sanders , j.s . , & fabian , a.c .
2007 , mnras , 381 , 1381 sanders , j.s . , &
fabian , a.c .
2008 , mnras , 390 , l93 sijacki , d. , springel , v. , di matteo , t. , & hernquist , l. 2007 , mnras , 380 , 877 spitzer , l. 1962 , the physics of fully ionized gases ( new york : wiley - interscience ) springel , v. et al .
2005 , nature , 435 , 629 staquet , c. , & godeferd , f.s .
1998 , j.fl.mech 360 , 295 tegmark , m. et al .
, 2004 , phys.rev.d , 69 , 103501 turner , l. , & pratt , j. 2002 , j. phys .
a , 35 , 781 vedula , p. , moser , r.d . ,
& zandonade , p.s .
2005 , phys .
17 , 55106 wang , w. , kaiser , c.r . ,
laing , r. , alexander , p. , pavlovski , g. , & knigge , c. 2009 , mnras , 397 , 1113 white , s.d.m .
, & rees , m.j .
1978 , mnras , 183 , 341 white , s.d.m . ,
navarro , j.f . ,
evrard , a.e . , & frenk , c. s. 1993 , nature , 366 , 429 willott , c.j . , rawlings , s. , blundell , k.m . ,
lacy , m. , & eales , s. 2001 , mnras , 322 , 536