Issue 
A&A
Volume 521, October 2010



Article Number  A25  
Number of page(s)  7  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201015210  
Published online  15 October 2010 
Potential vorticity dynamics in the framework of disk shallowwater theory
I. The Rossby wave instability
O. M. Umurhan^{1,2}
1  School of Mathematical Sciences, Queen Mary
University of London, London E1 4NS, UK
2 
Astronomy Department, City College of San Francisco,
San Francisco, CA 94112, USA
Received 14 June 2010 / Accepted 13 July 2010
Abstract
Context. The Rossby wave instability in astrophysical disks
is a potentially important mechanism for driving angular momentum
transport in disks.
Aims. We attempt to more clearly understand this instability in
an approximate threedimensional disk model environment which we assume
to be a single homentropic annular layer we analyze using disk
shallowwater theory.
Methods. We consider the normal mode stability analysis of two
kinds of radial profiles of the mean potential vorticity: the first
type is a single step and the second kind is a symmetrical step of
finite width describing either a localized depression or peak of the
mean potential vorticity.
Results. For single potential vorticity steps we find there is
no instability. There is no instability when the symmetric step is a
localized peak. However, the Rossby wave instability occurs when the
symmetrical step profile is a depression, which, in turn, corresponds
to localized peaks in the mean enthalpy profile. This is in qualitative
agreement with previous twodimensional investigations of the
instability. For all potential vorticity depressions, instability
occurs for regions narrower than some maximum radial length scale. We
interpret the instability as resulting from the interaction of at least
two Rossby edgewaves.
Conclusions. We identify the Rossby wave instability in the
restricted threedimensional framework of disk shallow water theory.
Additional examinations of generalized barotropic flows are needed.
Viewing disk vortical instabilities from the conceptual perspective of
interacting edgewaves can be useful.
Key words: hydrodynamics  accretion, accretion disks  instabilities
1 Introduction
The Rossby wave instability (hereafter RWI) is a promising candidate mechanism to account for the observed anomalous transport of cold astrophysical disks^{}. The process involves the instability of waves in an environment in which there are radial variations in a potential vorticity quantity (PV hereafter). This effect is the disk analog of the Rayleigh instability of stably stratified barotropic shear flows that are familiar in meteorology and geophysical fluid dynamics. The original works wherein the RWI was proposed for disks (Lovelace et al. 1999; Li et al. 2000; Li et al. 2001) demonstrated the existence and evolution of the instability in vertically integrated disk models. The earliest known (to this author) identification of the relevance of PV dynamics to astrophysical settings is found in Lovelace & Hohlfeld (1978) containing a study of the stability of cold selfgravitating galactic disks subjected to radially localized perturbations. For the galactic disk study as well as the RWI, instability requires the existence of a minimum or maximum in the associated radial PVprofile of the disk^{}.Because of the potential importance of this mechanism for disks, it is worthwhile to consider and understand this instability in disk settings, which are three dimensional, at least in part. The goal of this study is to expose the machinery of this process as clearly as possible in the context of a threedimensional theory.
Disk shallow water theory (Umurhan 2008, and hereafter DSWtheory) is a model reduction of the disk equations that is threedimensional but asymptotic in the sense that azimuthal scales are much larger than the corresponding vertical and radial scales. It is essentially a model environment representing vortex dynamics occurring on thin annular sections of disks over timescales which are much longer than the local disk rotation time. The original motivation for this approximation was to develop a threedimensional framework to describe the dynamics of elongated vortices known to emerge in twodimensional studies such as that reported in Godon & Livio (1999). The benefit of using the DSWtheory is that it allows a transparent analysis of vortex dynamics free of other physical processes that are likely not to play a significant role  in particular, free of both acoustic and gravity mode oscillations. As such, by means of this framework one may develop a mechanical understanding of diskrelated vortex dynamics.
In Umurhan (2008), DSWtheory was applied to understand the emergence of one form of the stratorotational instability, which is an instability of a stably stratified Rayleigh stable shear profile in a channel (Yavneh 2001; Dubrulle et al. 2005). The nonormal flow boundary conditions at the walls of the domain bring into existence edgewaves (e.g., Goldreich et al. 1986) that propagate along the walls in opposite directions with respect to one another. However, because the edgewaves may also interact with one another across the domain, if the separation of the walls fall within a specific allowed range, then the phase speeds of the edgewaves can ``resonate'' with each other and become unstable (Umurhan 2006; Umurhan 2008).
This edgewave dynamical picture has been used to understand the development of many forms of potential vorticity driven geophysical fluid instabilities including, among others, the original problem considered by Rayleigh (Hoskins et al. 1985; Baines & Mitsudera 1994; Heifetz et al. 1999) as well as astrophysical processes such as the PapaloizouPringle instability (Papaloizou & Pringle 1984; Goldreich et al. 1985). In the geophysical fluid dynamics context, this is frequently referred to as the instability of counterpropagating Rossby waves (e.g., Heifetz et al. 1999). The use of edgewaves as an interpretive tool has precedence in the analysis of plasma instabilities. For example, the diocotron effect, which is an instability associated with finite thickness charge sheets subjected to drift, is rationalized to arise from the interaction of edgewaves counterpropagating along the surfaces of the charge sheet (Knauer 1966).
Figure 1: Sketches of profile plots considered in this work. a) Single defect/jump located at x=0. b) Double defect (symmetric step) with jumps located at . Double defect depicted is a depression whose mean (midplane) enthalpy profile corresponds to a peak. 

Open with DEXTER 
DSWtheory is an appropriate platform to examine the RWI. Indeed, an examination of the vortex structures emerging from nonlinear calculations of the instability (see Fig. 1 of Li et al. 2000, hereafter LFLC2000) shows that the coherent vortex structure is severely elongated in the azimuthal direction while being radially tightly confined. This structural quality is similar to the vortices reported in Godon & Livio (1999). The advantage of DSWtheory is that it can capture some essence of threedimensionality, while its disadvantage is that, because of its construction, small azimuthal scales cannot be represented. By contrast, the vertically integrated model environment used in the original RWI studies are able to capture small azimuthal scales while coming at the expense of representing threedimensionality.
In this study, the RWI is examined in the framework of DSWtheory. The stability of two analytically tractable mean PVprofiles is studied: (i) a single defect, or ``singlestep''; and (ii) a double defect, or ``symmetricstep'' (see Fig. 1). As in the case of the stratorotational instability, the RWI is shown to result from the interaction of Rossby edgewaves that propagate along the locations of the steps of the mean PVprofile. However, the instability manifests itself only in the symmetric step case. In Sect. 2, the equations of motion and their perturbations are presented. In Sect. 3, the equations are analyzed for two kinds of mean potential vorticity profiles including the one just described. Analytic solutions are determined throughout by considering the quasilinear problem posed by considering the special case , where is the ratio of specific heats. Section 4 concludes with comments and a discussion about recent threedimensional RWI calculations (Meheut et al. 2010).
2 Equations and potential vorticity conservation
The asymptotic equations upon which DSWtheory is based, describing the evolution of annular disk sections, was developed in Umurhan (2008). In that study, the equations of motion for an annular section in which the vertical/radial extent is much smaller that the azimuthal scale R was developed as an asymptotic series expansion in powers of . The timescales of the dynamics () is comparatively long with respect to the local rotation time ( ) of the annular disk centered on R_{0} i.e., . To ensure a nontrivial asymptotic balance, the radial and vertical velocities viewed in the reference frame of the annulus rotating with must scale as where is the sound speed of the gas, while the corresponding deviations of the azimuthal velocities must scale as . This disparity in aspect and velocity ratios is an asymptotic description of azimuthally elongated slowly overturning vortices. Furthermore, the reduced equations of motion in the annular disk section possess quasigeostrophic characteristics familiar in meteorological studies.
Single constant entropy layers are considered when the equation of state is polytropic, i.e., it is assumed that everywhere at all times. This corresponds to a situation in which the specific entropy inside, defined by , is everywhere a constant. This homentropic layer is considered in the following. If there are radial variations of the mean height of the disk (h), then the vertically integrated entropy of the disk also exhibits radial variations even though the specific entropy is constant. This point is made to keep in mind how the equations analyzed here compares to the vertically integrated twodimensional equations studied in Li et al. (2000).
After nondimensionalizing all the quantities according to the above,
the leading order equations
describing the evolution of varicose disturbances
of a single constant specific entropy layer are
where x,y represent the nondimensional radial and azimuthal coordinates and t is the nondimensionalized time (Umurhan 2008). In this construction, all quantities appearing are functions of these two coordinates and time. The radial velocity is u. The total time derivative is
The background shear , is linear in the radial coordinate. For a Keplerian disk, q=3/2. Thus, the azimuthal velocity v appearing here is understood to be a deviation from this basic background state as viewed in the local rotating frame. The quantity may be understood to be the enthalpy of the disk midplane and is given by
(5) 
where h=h(x,y,t) corresponds to the height of the disk measured from the disk midplane. The height h is the position at which the enthalpy goes to zero in the constant (specific) entropy environment. Although is technically the midplane enthalpy, it is referred to hereafter as simply the enthalpy with full understanding that it is really its value at the midplane.
Equation (1) is the radial momentum balance equation but in this asymptotic limit it says that all perturbations are radially geostrophic. Equation (2) is the azimuthal momentum balance and, unlike the radial version, no geostrophic balance is implied. Equations (3) and (4) represent, respectively, the motion of the surface h and an equation of state for the gas (see below). In the derivation contained in Umurhan (2008), these two equations result from (i) exploiting the hydrostasy of perturbations at all times and (ii) an analysis of the energy equation, which makes use of the linear independence of powers of z.
The disk vertical velocity, though not explicitly appearing, is odd symmetric with respect to the disk midplane because of the assumed varicose symmetry of the disk response. In the development of the DSWtheory its functional form is shown to be . In the reduction leading to these equations only appears as well as the explicit evaluation of the velocity of the layer at the height z=h. Additional details are found in Umurhan (2008).
The combination of these equations of motion leads to the conservation of potential vorticity Qgiven by
(6) 
The next section considers the fate of normal modes for potential vorticity (or ``PV'' for short) profiles in which there is either a single or a double defect in the steady PV profile .
3 Analytical study
The analysis in this work considers the fate of disturbances where
.
This
is chosen because the resulting equations, both steady and normal mode disturances, may
be analyzed without intensive use of numerical methods. General more realistic values of
require numerical evaluation of both the steady and timedepenedent state and
this will be reserved for a future study. For
,
the PV expression takes on
the quasilinear form
(7) 
The following sections consider the two analytically tractable defect profiles and their normal mode responses.
3.1 Single defect
In this study, the term defect is used to indicate places where steps
occur in the mean PVprofile. The following simple single defect,
(8) 
is studied. The value of Q_{0} is constrained to be greater than zero. Solutions are sought that are bounded as . This means solving for the mean enthalpy
(9) 
on either side of x=0(and where q = 3/2 has been explicitly used). The particular boundary conditions for the profiles at infinity are (i) as and (ii) as . At the interface x=0, it is also required that and its radial gradient are matched, the latter being identical to requiring the matching of the azimuthal velocities from both sides. The solution for the mean enthalphy is therefore
(10) 
where
(11) 
This result may be read to mean that the mean height decreases when Q_{0} > 1 and increases when 0<Q_{0}<1. It follows from Eq. (1) that
(12) 
which means that the deviation velocity at x=0 is given by . As , the mean velocity deviation goes to zero. The mean states developed here are qualitatively sketched in Fig. 1.
Infinitesimal perturbations of this steady state are
introduced by writing for each dependent quantity
,
where F'(x,y,t) is small
compared to the basic state .
This means that
the potential vorticity equation is
(13) 
where Q' is given as
Because entails a deltafunction at x=0, special wellknown treatments must be utilized to solve these equations (for example, see Drazin & Reid 1982). For this type of problem, the procedure is detailed, for example, in Umurhan (2008), and may be summarized as follows: (i) solve Q' = 0 on either side of the defect; (ii) match coming in from both sides; and (iii) make sure that the normal velocities u' are matched across the boundary. The last of these may be straightforwardly enforced by inspecting the linearized version of Eq. (2),
(14) 
and requiring that the two solutions coming in from either side of x=0 satisfy the requirement that the normal velocities u' are the same there.
Normal mode solutions of the form
are assumed.
Solutions with c>0 correspond to waves that propagate in the forward y direction
and, correspondingly, in the negative y direction for c<0.
It follows that the solution for ,
away from the two boundaries,
satisfying the condition that their amplitudes are equal at x=0and,
as
,
is
(15) 
where A_{0} is an arbitrary constant. The requirement that the normal perturbation velocities, u', match at the defect location means ensuring that
(16) 
Satisfaction of this condition leads to an eigenvalue condition on the wavespeed c which, after some algebra, reduces to the simple form
For given above, it follows that that . The wave here is a Rossby edgewave phenomenon (e.g. Baines Mitsudera 1994) and the eigenvalue determined here describes the speed with which the edgewave propagates. If the defect were placed at some position other than x=0, then the observed speed of the edgewave would be Dopplershifted according to the Keplerian velocity at the position of the defect. Thus, if the defect were instead at x=x_{0}, the wavespeed would then be given by .
3.2 Double defect: the symmetric step
This subsection is an analysis of a symmetrical defect pair. The previous section demonstrates how a single defect produces a single Rossby edgewave along the defect. For two defects, therefore, two Rossby edgewaves are expected and it is under these circumstances that instability may emerge according to the principles of interacting waves (Hayashi & Young 1985; Baines & Mitsudera 1994; Heifetz et al. 1999).
The mean PVprofile appropriate for a double symmetric defect
is
(18) 
where . The edges of the defect are located at and separated by a distance , where is the region's halfwidth. Following the same procedures as outlined in the previous section and requiring that as , the mean enthalpy is found to be
(19) 
The value of in this case is
(20) 
It means that for 0<Q_{0}<1, the enthalpy deviation is a peak or bump, while for Q_{0}>1the enthalpy structure is a depression. The corresponding mean velocity deviations are
(21) 
Defining , it follows that at the defects . It is useful to define , i.e., the ratio of the height of the peak mean enthalpy to its value as ,
A qualitative sketch of the mean states are depicted in Fig. 1.
Because there are two defects, individual edgewaves
propagate along each of them, which interact with
each other as well. A normalmode analysis, as
executed in the previous section, is performed. The details
are presented in the Appendix. The resulting wavespeed csatisfies the relationship
where
and is defined for convenience. The results of this are plotted in Fig. 3. The first clear result is that instability ( ) may only occur for enthalpy bumps, which is to say that instability occurs for PVdepressions, 0<Q_{0}<1. The boundary between stable and unstable parameter regimes as a function of Q_{0} is given by . Thus, for PVdepressions, instability occurs for . For very small but nonzero values of Q_{0}, the results show that . Defining instead as the marginally stable value of the PVdepression then its asymptotic form has the functional dependence for . Figure 4 shows plots of the wavespeeds as a function of for representative values of a PV depression and peak.
4 Comments and discussion
It has been shown that the RWI is recovered in the framework of DSWtheory. The RWI emerges when the idealized symmetricstep mean (radial) potentialvorticity profile is a depression (or equivalently, when the mean radial enthalpy profile is a peak). However, single stepped mean PV profiles do not support the instability. The RWI is interpreted as being caused by the interaction of two Rossbyedgewaves propagating azimuthally along each defect of the symmetric step PVprofile. These types of dynamical processes are a generic feature of many vortical dynamics in geophysical flows. Below is a collection of comments about these results including comparisons between the results of this work and the original RWI studies:
1. A weakness of the analysis performed here is that an unrealistic
value of the ratio of specific heats was used: a value of
places the fluid in a thermodynamic regime somewhere between a normal gas
and an incompressible fluid (the latter being
,
Salby 1989).
The main reason for this choice
was to develop the analysis analytically as most other values of require numerical methods to study both the steady states and the perturbation evolution.
Among other things, developing analytical solutions are important for benchmarking future analyses for
more reasonable values of .
But most of all, the analytical solutions offer
the most transparency in interpreting the results.
Preliminary calculations for realistic values of
show that the main mechanical interpretations
and qualitative solutions, as well as conditions for onset, are not considerably altered.
This is the subject
of a study in preparation.
2. The limiting form of the DSWtheory used to study the RWI considers the response of a single homentropic disk layer. In this picture, the mean vertical height of the disk is a function of the radial coordinate, i.e., . Instability only occurs when the enthalpy of the disk layer (which is ) has a peaked structure. In LFLC2000, the RWI is demonstrated in a vertically integrated model of a disk. Similar in quality to the results here, instability occurs when there is a peak in the mean radial pressure profile (or, relatedly, the surface density profile). In contrast, however, instability occurs for constant specific entropy profiles. This might imply that these results conflict with one another.
Because the disk model considered in LFLC2000 is vertically integrated, all
of the fluid quantities considered there are also understood to be vertically integrated,
including the entropy. In the DSWtheory framework studied here, despite
the specific entropy being constant throughout the layer, its vertical integration
is not constant if the mean height of the disk is not uniform. That is to say
if one defines
as the vertically integrated specific entropy of the disk, then it is selfevident that is not uniform if varies with radius x. As such, the results of this study should be compared to the nonhomentropic incarnations of the RWI discussed in the original studies, for example, as in Li et al. (2001).
3. With respect to the last point, the version of the RWI uncovered in the analysis of the symmetric step is qualitatively similar to the homentropic Gaussian bump profile (``HGB'') considered in LFLC2000, where HGB describes the quality of the radial dependence of the mean pressure. However, as stated above, the symmetric step profile studied in this work is effectively nonhomentropic when viewed in terms of the vertically integrated disk model considered in Li et al. (2000). The symmetric step profile might, therefore, most closely correspond to a nonhomentropic Gaussian bump profile (``NGB'') in the framework of the vertically integrated model of LFLC2000. A NGB profile was examined in the nonlinear regime in Li et al. (2001), and, in a qualitative sense, its evolution was not found to differ considerably.
One may, therefore, reasonably compare the conditions in which the RWI arises in the HGB setting to the instability arising here for the symmetricstep profile. In the analysis presented in LFLC2000, it was shown that for the HGB profile a minimum amplitude in the pressure bump was required for instability to occur. It was indicated in Fig. 9. of that study that the minimum variation required in the associated surface mass density for instability to occur, on a radial scale similar to the thickness of the disk, must be between approximately 10% to 20% (LFLC2000). According to this study, for any given value of 0 <Q_{0} < 1, instability occurs for . This means that for a given peak value of the midplane enthalpy there is a maximum width size beyond which instability ceases. In this sense, the results here imply that there is no minimum required amplitude of the midplane enthalpy peak (and, consequently, the corresponding associated surface density peak) for instability to occur.
This does not necessarily mean that these results disagree with each other. Indeed in Fig. 3, it can be seen that for fixed values of , there is a value of Q_{0} smaller than some critical value for the system to be unstable. For example, for , instability occurs for Q_{0} < 0.5. According to Fig. 2, this means that for these parameter values the minimum midplane enthalpy peak variation over the mean must be at least 25% for instability to occur.
Figure 2: Plot of , the ratio of the minimum or maximum of the mean enthalpy to its value as (symmetric step profile only). Three values of Q_{0} are shown: Q_{0} = 0.25,0.5 (depression) and Q_{0} = 3.25 (peak). 

Open with DEXTER 
Figure 3: Stability boundaries for the double defect problem in the plane. Unstable region depicted by shaded zone. Instability occurs for all values of enthalpy bumps (i.e., PVdepressions: 0<Q_{0}<1). No instability expected for enthalpy depressions. The boundary separates stable and unstable regions. Critical values of Q_{0} corresponding to instability goes like for . 

Open with DEXTER 
This still leaves open the question about the origin of the initial equilibrium. In the general scenario proposed in LFLC2000, matter is accumulated in a piecemeal way, either by means of the coupling to largescale magnetic field processes or due to smallscale turbulence (possibility MHD driven). These possibilities are not disputed here. Instead attention should be paid to the way in which matter is accumulated. Reference to Fig. 3 best illustrates this point. The way that matter accumulates depends upn the path within the Q_{0} plane along which the accumulation process takes place. For example, if the process involves an increase in (the region size) for fixed values of 0<Q_{0}<1, then instability will occur before any sizeable enthalpy peak is reached because any such path (starting from ) is unstable to the RWI. On the other hand, if the initial equilibrium Q_{0} is constructed on a path of constant , then the system may be brought into the unstable regime with a finite large amplitude of the enthalpy peak variation. Thus, the way in which the equilibrium is constructed is crucial to assessing how and on what scale the instability will operate. This point deserves further attention in the future.
4. The single step profile in the mean PV considered in this analysis shows no instability. In contrast, the single jump profiles studied in LFLC2000 do show the emergence of instability. The reason for this difference is not entirely clear. One possibility is that the single step profile considered here may not be a faithful analog of the single jump profile studied in LFLC2000. Indeed, in the generalized picture of edgewave dynamics of Baines & Mitsudera (1994), one needs a minimum of two edgewaves for instability to happen. The singledefect model of this study can only support a single edgewave and, therefore, no instability is possible. A corresponding analog representation of the single jump profile considered in LFLC2000 probably needs to support at least two edgewaves in order to see the instability expressed. This is, however, conjecture at this stage and more analysis is needed.
5. Instability can be best rationalized in terms of the interaction of the two individual edgewaves propagating along each defect of the symmetric step profile. This picture is consistent with the HayashiYoung criterion for wave instability (Hayashi & Young 1987), which can be summarized as follows: instability of two waves separated by an evanescent layer occurs when the two waves (i) propagate opposite to each other; (ii) have almost the same Dopplershifted frequency; and (iii) can interact with one another. Interaction in this case means that the edgewave at one defect contributes to the perturbation velocity profile at the other defect.
All of these conditions for wave instability
are satisfied in the problem encountered here.
In particular,
one can view the Rossbyedgewaves propagating along each
defect as though they were disturbances
of isolated defects. Thus,
with respect
to the wavespeed calculation worked out for single isolated steps in PV
(Sect. 3.1), one can write down the wavespeeds associated
with each defect as
for the wave at and
at . The wavespeeds and the values of c appropriate for the doubledefect problem are plotted in Fig. 4^{}. It can be seen in Fig. 4b that instability occurs when the possibility exists for the two isolated waves to have nearly the same (and opposite) wavespeeds. In the rest frame x=0, this occurs when each of the wavespeeds are nearly equal to zero. In Fig. 4a, where instability does not occur, the wavespeeds are not nearly equal enough (they are not close to zero), hence instability cannot occur.
Figure 4: wavespeeds as a function of for symmetric step profile: a) Q_{0} = 2.5, b) Q_{0} = 0.5. Also shown are the speeds (discussed in the text, Sect. 4) corresponding to the edgewaves that would propagate along each defect as if each defect were treated in isolation. Instability occurs in b) only when these two wavespeeds are nearly equal, which is consistent with the HayashiYoung criterion for wave instability. 

Open with DEXTER 
In this particular doubledefect setup it follows that one can roughly predict
the onset of instability to be when the isolated wavespeeds are
nearly zero in the reference frame of an observer at x=0. Using the
formula for
given above,
happens very nearly when
Clearly this can never happen when Q_{0} > 1 as must be positive. Furthermore, when this rough criterion means that the critical value of the separation scale is approximately
which is the same trend observed in Fig. 3.
The importance of this perspective is that one may approximately assess whether or not a more complicated looking PVprofile can become unstable to the RWI instability by appealing to this conceptually simple mechanistic view. For example, if a more sophisticated PVprofile were treated as if it were a series of PVsteps, then an approximate prediction for the onset of instability would only require calculating edgewave speeds along each defect taken as though each were in isolation. Onset conditions would be predicted for values of the parameters in which any two of these wavespeeds are nearly equal.
6. The asymptotic limit that the DSWtheory exploits is one in which the azimuthal length scales are much longer than the radial ones. In this limit there is no information as to how wavespeeds depend upon the horizontal wavenumber . In the theory developed here the wavespeed c is independent of . However, it is known that in the vertically integrated configurations investigated in the original RWI studies there is a clear dependence of c with .
On more concrete terms, in the numerical solution shown in LFLC2000 there is obviously a fastest growing azimuthal wavenumber of the instability, which leads to the final steady pattern state reached by the system. In contrast, in the approximate study done here for the symmetricstep profile, once stability is breached all wavenumbers become unstable and there is no highwavenumber cut off. The growth rates, , obviously increase without bound in proportion to  signaling the breakdown of the theory.
To correct this issue, the DSWtheory would
have to be appropriately extended by allowing for
vertical structure in the horizontal velocities (i.e., inclusion of the
dependence z^{2}) and some deviations from the linear dependenc
of the vertical velocity
with the z coordinate (i.e., an additional z^{3} dependency).
The resulting
next order correction to the
wavespeed in the extended DSWtheory
will show that the general dispersion relationship takes the form
where is a function of the parameters Q_{0} and . Preliminary calculations show that is greater than zero, which means there is a highwavenumber cutoff beyond which there would be no instability in the extended theory. The DSWtheory may be extended to higher orders in by asymptotically developing corrections for c near values of , which occur for values of near . Further investigation is needed in this regard.
7. Threedimensional numerical simulations of the RWI were announced during the preparation of this manuscript. Meheut et al. (2010) report on threedimensional simulations of a constant entropy disk with midplane symmetric perturbations with the main goal of following the evolution of a Gaussian bump profile in the initial pressure field. Notable features of the results show that there is a significant amount of threedimensional structure in the horizontal velocities as well as deviations from a linear z dependence of the vertical velocities. The resultant structures show interesting and complex threedimensional structure in the streamlines of the flow. In contrast, in the asymptotic threedimensional theory used here the horizontal velocities have no vertical variation and the vertical velocity is linear in z. Thus, the asymptotic strategy used here is far from matching the overall structure contained in the results of these particular simulations. The DSWtheory (or an analogous approach) may be extended to capture the quality of the dynamics contained in these numerical simulations much in line with the previous comments made in this regard. Moreover, it is important to understand how the threedimensional structure emerges in these simulations, i.e., whether it arises from other secondary barotropic processes or if they are dynamics that are tied to the primary quasitwo dimensional dynamics driven by the primary instability. Rationalization of the threedimensional dynamics in terms of the Rossby edgewave perspective utilized in this work could be a natural first step to this end.
AcknowledgementsThe author thanks Yiannis Tsapras, Oliver Gressel, Joe Barranco, and James Cho for discussions relating to this work. The author also thanks the anonymous referee for pointing out the analogy between the diocotron effect and the RWI. The author is indebted to the Kavli Institute of Theoretical Physics and the program Exoplanets Rising: Astronomy and Planetary Science at the Crossroads, where portions of this work were performed.
Appendix A: Doubledefect normalmode calculation
Following the procedure performed for the single defect, the equation Q' = 0 is solved separately in the three zones subject to the condition that as . For and , the equation to be solved is(A.1) 
while for , it is
(A.2) 
The solution for in each of the three zones is
(A.3) 
The solution as presented satisfies the continuity of at the boundaries . The remaining conditions are the continuity of the normal velocities at the location of the defects, which amounts to simultaneously satisfying
(A.4) 
and
(A.5) 
These two equations result in the matrix equation
where is a matrix and . Nontrivial satisfaction of this matrix equation requires enforcing . This results in the wavespeed relationship found in Eq. (23).
References
 Baines, P. G., & Mitsudera, H. 1994, J. Fluid Mech., 276, 327 [NASA ADS] [CrossRef] [Google Scholar]
 Dubrulle, B., Marie, L. , Normand, Ch., et al. 2004, A&A, 429, 1 [Google Scholar]
 Friedman, B. 1956, Principles and Techniques of Applied Mathematics (New York: Wiley) [Google Scholar]
 Godon, P., & Livio, M. 1999, ApJ, 523, 350 [NASA ADS] [CrossRef] [Google Scholar]
 Goldreich, P., & LyndenBell, D. 1965, MNRAS, 130, 125 [NASA ADS] [CrossRef] [Google Scholar]
 Goldreich, P., Goodman, J., & Narayan, R. 1986, MNRAS, 221, 339 [NASA ADS] [Google Scholar]
 Hayashi, Y.Y., & Young, W. R. 1987, J. Fluid Mech., 184, 477 [NASA ADS] [CrossRef] [Google Scholar]
 Heifetz, E., Bishop, C. H., & Alpert 1999, Q. J. R. Meteorol. Soc., 125, 2835 [NASA ADS] [CrossRef] [Google Scholar]
 Lovelace, R. V. E., & Hohlfeld, R. G. 1978, ApJ, 221, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Hoskins, B. J., McIntyre, M. E., & Robertson, A. W. 1985, Q. J. R. Meteorol. Soc., 111, 877 [Google Scholar]
 Knauer, W. 1966, J. App. Phys., 37, 602 [Google Scholar]
 Li, H., Finn, J. M., Lovelace, R. V. E., & Colgate, S. A. 2000, LFLC2000, ApJ, 533, 1023 [Google Scholar]
 Li, H., Colgate, S. A., Wendroff, B., & Liska R. 2001, ApJ, 551, 874 [NASA ADS] [CrossRef] [Google Scholar]
 Lovelace, R. V. E., Li, H., Colgate, S. A., & Nelson, A. F. 1999, ApJ, 513, 805 [NASA ADS] [CrossRef] [Google Scholar]
 Meheut, H., Casse, F., Varniere, P., & Tagger, M. 2010, A&A, 516, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Papaloizou, J. C. B., & Pringle, J. E. 1984, MNRAS, 208, 721 [NASA ADS] [CrossRef] [Google Scholar]
 Pedlosky, J. 1987, Geophysical Fluid Dynamics, 2nd edn. (New York: SpringerVerlag) [Google Scholar]
 Lord Rayleigh 1880, Proc. Royal Math. Soc., 9, 57 [Google Scholar]
 Salby, M. L. 1989, Tellus, 41A, 48 [NASA ADS] [CrossRef] [Google Scholar]
 Umurhan, O. M. 2006, MNRAS, 365, 85 [NASA ADS] [CrossRef] [Google Scholar]
 Umurhan, O. M. 2008, A&A, 489, 953 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Yavneh I., McWilliams J. C., & Molemaker M. J. 2001 J. Fluid Mech., 448, 1 [Google Scholar]
Footnotes
 ... disks^{}
 I.e., nonmagnetized disks.
 ... disk^{}
 Note that in Lovelace & Hohlfeld (1978) the associated quantity f(r) corresponds roughly to the inverse of the PV.
 ...
^{}  The single wavespeed formula quoted above for c_{}is written in terms of the formula given in Eq. (17) augmented by the Dopplershift consistent with the the local Keplerian speed at . The isolated edge wavespeed c_{+} is correspondingly written taking into consideration the symmetry of the mean PVprofile as viewed at , which is the opposite of that at .
All Figures
Figure 1: Sketches of profile plots considered in this work. a) Single defect/jump located at x=0. b) Double defect (symmetric step) with jumps located at . Double defect depicted is a depression whose mean (midplane) enthalpy profile corresponds to a peak. 

Open with DEXTER  
In the text 
Figure 2: Plot of , the ratio of the minimum or maximum of the mean enthalpy to its value as (symmetric step profile only). Three values of Q_{0} are shown: Q_{0} = 0.25,0.5 (depression) and Q_{0} = 3.25 (peak). 

Open with DEXTER  
In the text 
Figure 3: Stability boundaries for the double defect problem in the plane. Unstable region depicted by shaded zone. Instability occurs for all values of enthalpy bumps (i.e., PVdepressions: 0<Q_{0}<1). No instability expected for enthalpy depressions. The boundary separates stable and unstable regions. Critical values of Q_{0} corresponding to instability goes like for . 

Open with DEXTER  
In the text 
Figure 4: wavespeeds as a function of for symmetric step profile: a) Q_{0} = 2.5, b) Q_{0} = 0.5. Also shown are the speeds (discussed in the text, Sect. 4) corresponding to the edgewaves that would propagate along each defect as if each defect were treated in isolation. Instability occurs in b) only when these two wavespeeds are nearly equal, which is consistent with the HayashiYoung criterion for wave instability. 

Open with DEXTER  
In the text 
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.