12.07.2015 Views

A Magnetic α ω Dynamo in Active Galactic Nuclei Disks - NMT Physics

A Magnetic α ω Dynamo in Active Galactic Nuclei Disks - NMT Physics

A Magnetic α ω Dynamo in Active Galactic Nuclei Disks - NMT Physics

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

– 21 –rotat<strong>in</strong>g frame with the orig<strong>in</strong> at the po<strong>in</strong>t r = r 0 such that the x-axis is directed radially outward,the y-axis is directed <strong>in</strong> the positive azimuthal direction, and the z-axis is perpendicular to the diskplane. Then, the effective gravitational and centrifugal potential <strong>in</strong> the Keplerian rotat<strong>in</strong>g frame<strong>in</strong> the neighborhood of the po<strong>in</strong>t r = r 0 is∆Φ = GM2r03 (z 2 − 3x 2 ). (21)The thermal energy of the hot column of gas is a fraction of the loss of k<strong>in</strong>etic energy of the stardue to the hydrodynamic collision with the disk. This latter energy loss dur<strong>in</strong>g one passage isF drag 2H = 2HπR 2 ⊙ρv 2 ∗ = πR 2 ⊙Σv 2 ∗.Without a hydrodynamic simulation <strong>in</strong> 3-dimensions an accurate description is miss<strong>in</strong>g. Neverthelessit is sufficient to approximate the solution as that fraction of the disk matter that hasan <strong>in</strong>ternal energy density, ǫ shk greater than that of the ambient disk by that factor such that itwill rise to a height, z, determ<strong>in</strong>ed by its potential energy, or ∆Φ = GMz22r 3 (equation (21)). S<strong>in</strong>ce∆S/S o = ǫ shk /ǫ o , where ǫ o ≃ ∆Φ(H) = GMH22r 3 , then <strong>in</strong> order for a plume to rise above the diskmid-plane to a height, l,( lH) 2≃(ǫshk)≃ǫ o( ) 2 ( ) 2 vK R⊙( r) ( ) 2 2 R⊙≃ . (22)c s R shk H R shkWe are concerned with plumes that rise well above the disk so that they can expand horizontallyby a factor several times the plume’s orig<strong>in</strong>al radius. In this case the moment of <strong>in</strong>ertia of theplume about its own axis will be <strong>in</strong>creased by several times before fall<strong>in</strong>g back to the disk. Thiscauses the plume to reduce its own rotation rate relative to the frame of the disk, that is to untwistrelative to that frame. For this expansion to take place, the plume must rise roughly ∼ 2H abovethe disk, or l ≃ 3H. At this height the pressure of the hydrostatic isothermal atmosphere with thedensity profile as ∝ exp(−z 2 /H 2 ) becomes negligible compared to that of the plume, and so thehot gas of the plume can expand both vertically and horizontally as a free expansion. With this lwe getR shk≃ 1 rR ⊙ 3 H . (23)Us<strong>in</strong>g expression (8) for H/r at r ≤ r ab we obta<strong>in</strong>R shk ≃ 0.24H( αss) ( ) (2/21−26/21lE(ǫ) 26/21−19/21 r M 8 1 −0.01 0.1 0.1 r ab√ ) −23rg, (24)Thus a plume, start<strong>in</strong>g from a size, R shk < H will expand to a size ≃ 2H both vertically andhorizontally, thus produc<strong>in</strong>g a near spherical bubble with radius r p = H above the disk. Post shockexpansion will <strong>in</strong>crease the estimate of R shk somewhat. For simplicity we will use R shk = H/2 forestimates of the toroidal flux entra<strong>in</strong>ed <strong>in</strong> the plumes <strong>in</strong> paper II. This is our standard plume.r

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!