A&A 486, 341-345 (2008)
DOI: 10.1051/0004-6361:20077484
A. Sternberg^{1} - O. M. Umurhan^{1,2,3} - Y. Gil^{1} - O. Regev^{1,4}
1 - Department of Physics, Technion-Israel Institute of
Technology,
32000 Haifa, Israel
2 - Department of Geophysics and Space Sciences, Tel-Aviv
University, Tel-Aviv, Israel
3 - Department of Astronomy, City College of San Francisco,
San Francisco, CA 94112, USA
4 - Department of Astronomy, Columbia University,
New York, NY 10025, USA
Received 15 March 2007 / Accepted 1 April 2008
Abstract
We examine the hydrodynamic response of the inviscid small
shearing box model of a midplane section of a rotationally
supported astrophysical disk. We formulate
an energy functional for the general nonlinear problem. We find that the fate of disturbances is
related to the conservation of this quantity which, in turn, depends on the
boundary conditions utilized:
is conserved for
channel boundary conditions, while it is not conserved,
in general, for shearing box conditions.
Linearized disturbances
subject to channel boundary conditions have normal-modes described
by Bessel Functions and are qualitatively governed by a quantity ,
which is a measure of the ratio between the azimuthal and vertical
wavelengths.
Inertial oscillations ensue if
- otherwise
disturbances must be treated generally
as an initial value
problem. We reflect upon these results
and offer a speculation.
Key words: accretion, accretion disks - hydrodynamics
Far from answering whether or not such transitions do exist in these flows, we instead identify two clues that could be used to cast some light onto this question in the future. To this end, we reconsider the dynamics of inviscid anticylonic Rayleigh stable flows in the SSB model, i.e. rpCf (rotating plane Couette flow); subject to shearing boundary conditions (SBC, Goldreich & Lynden-Bell 1965; Rogallo 1981; Knobloch 1984; Korycansky 1992) or channel boundary conditions (CBC, Yecko 2004; Mukhopadhyay et al. 2004; Lesur & Longaretti 2005; Rincon et al. 2007).
The first result is a general nonlinear feature pertaining to the energetics of such flows. If one defines as the total energy of the flow in the SSB, whose dynamical constituent includes the total kinetic energy (i.e., the energy formed by considering the total velocity, which is to say the velocity fluctuations plus background shear flow), the fate of depends on the boundary conditions employed. The is conserved for (but not limited to) (i) disturbances subject to CBC; and (ii) azimuthally and vertically periodic, radially localized disturbances in radially unbounded domains. On the other hand, disturbances subject to SBC do not identically conserve and its temporal quality depends critically upon the dynamics within the domain. The nonconservation of for flows subject to SBC relates to the feature that the kinetic energy of a fluid parcel jumps when it crosses the (quasi) periodic radial boundaries.
For our second result, we find that there are qualitative differences in the
behavior of linearized disturbances subject to CBC,
depending on whether
or not the quantity
is less than or greater than 1.
This parameter is defined as
(1) |
In flows subject to CBC, with given vertical and azimuthal wavenumbers, if the system supports a countably infinite number of inertial normal modes described by Bessel functions oscillating without growth or decay. Otherwise, if , the system allows for only one eigenmode, meaning that in general the dynamical response cannot be characterized by discrete normal modes in which case such linearized disturbances must be treated as an initial value problem, a situation well-known in studies of plane-Couette flow (pCf, Schmid & Henningson 2000).
We would like to note that the two main results we report on here may be clues in explaining the results of Lesur & Longaretti (2005) who show that there is a subcritical transition beyond the Rayleigh line in rpCf when SBC conditions are applied, whereas no such transition was observed for such flows under CBC boundary conditions. On a related note, as well, Lithwick (2007) demonstrated the existence of steady vortex structures in 3D rpCf flow for Keplerian shear profiles when SBC are applied. We return to this in the discussion.
The flow on these scales is equivalently known as incompressible rpCf (e.g., Nagata 1986; Yecko 2004). In this form, these equations are the inviscid limit of those considered by Longaretti (2002); Mukhopadyay et al. (2005); Yecko (2004). Here, x corresponds to the radial (shearwise) coordinate, y corresponds to the azimuthal (streamwise) coordinate and z corresponds to the vertical (spanwise) coordinate - with the corresponding velocity disturbances . These velocities represent deviations over the steady Keplerian flow (as manifesting itself in this rotating frame) given by . The quantity is sometimes also referred to as the Coriolis parameter, and in these nondimensionalized units is 1. The local shear gradient is defined as the exponent of the general rotation law .
Boundary conditions. We first present and discuss the means by which one administers SBC in the SSB. The flow variables , and p are (a) periodic in the azimuthal and vertical directions on scales L_{y}=1 and L_{z}=1, respectively, and (b) are simply periodic (on scale L_{x}=1) in the shearwise (radial) direction with respect to a coordinate frame that moves with the local shear (e.g., Lynden-Bell & Goldreich 1965; Rogallo 1981). We refer to this coordinate frame as the shearing coordinate frame (SC). The coordinate frame of an observer in the rotating frame, will be referred to as the non-shearing coordinate frame (NSC). The fundamental Eqs. (2)-(5) have been expressed in the NSC frame.
If (X,Y,Z,T) represent the independent variables in the SC frame then we say they are related to the NSC variables by If f is any dynamical quantity then in the SC frame the periodic boundary conditions in the shearwise, azimuthal and vertical directions respectively are and this means, in general, that f(X,Y,Z) = f(X+L_{x},Y+L_{y},Z+L_{z}). As expressed in the NSC frame, the periodicity in the azimuthal and vertical directions appears like f(x,y,z) = f(x,y+L_{y},z); and f(x,y,z) = f(x,y,z+L_{z});while the periodicity in the shearwise direction is Because of the periodicity in the azimuthal direction this statement takes the form where . Note that this means that only after precisely integer values of , (i.e., when ) are the SC and NSC frames coincident, meaning every time units. Evidentally, the SBC are time dependent when viewed by an observer in the NSC frame.
When enforcing channel boundary conditions (CBC) on (2)-(5) we require that all quantities are periodic in the vertical and azimuthal directions (in the NSC frame) while requiring no normal-flow at the channel boundaries. This latter statement amounts to u = 0 at x = 0,L_{x}.
A pair of integral statements.
The dynamical Eqs. (2)-(5)
describe the evolution
of velocities that are disturbances within a steady
rotationally supported flow which, on the SSB scales, is
the linear shear.
Thus, the steady solution,
,
is
interpreted as the undisturbed state.
We now consider the total velocity ,
defined as,
As such the governing equations of motion
((2)-(5))
are more concisely written in vector form,
By contrast, we follow the same procedure as above but consider
an energy functional comprised only of the disturbance velocities.
By taking the inner product of (3)-(5)
with the disturbance velocity ,
followed by the use
of (2) and either the
SBC or
CBC boundary conditions, one can integrate the resulting
inner product over the domain
to find
Inspection of the Reynolds-Orr
relationship says that
as long as
there is some correlation between the radial and azimuthal
velocities, the disturbance energy must fluctuate in time.
However, from the definitions of
and E we have
To appreciate the significance of these relationships
we first apply CBC to (9) and find that
We note that if instead of the CBC we require that the flow be periodic in y and z (again) but that all flow quantities decay sufficiently fast as , then the boundary term in (9) also vanishes and (12) again follows. The same is true if the vertical domain is unbounded and disturbances are vertically localized as well, i.e. (12) also follows if disturbances similarly decay as .
Matters are less clear when SBC are employed because
is no longer
conserved, evolving according to
(with boundaries of the periodic box set at
),
(17) |
(18) |
We observe that if , as appearing, is a real quantity then there is no way for the quantization condition to be met by real values of (other than the single value ) according to the known behavior of the Bessel functions (Abramowitz & Stegun 1972). This means that if , then there is only one normal mode solution for this problem. For there exists a countably infinite set of normal-mode solutions (overtones) as per the theory of Bessel functions (see again Abramowitz & Stegun 1972).
Nevertheless, inferring
the dynamical content of the disturbances, including the shape of
the eigenfunctions and the general formulae for the dispersion
relationship, is not easy. One can more transparently study the
behavior of these features, including their character change, if an asymptotic solution for
is developed instead for
the small limiting values of
and .
We proceed
by assuming
and,
followed by the expansion procedure
and
to the solution of (15).
We find that the lowest order solution depends critically on the
parameter
defined by
which is really a reexpression of .
The leading order solution
to u_{0} is
In a more general sense, when , aside from the mode , the response of the flow must be considered as an initial value problem since the system supports no other discrete normal modes. When , then there are two branches of normal modes (i.e., overtones) clustering about and . In general, the frequencies around which the overtones cluster correspond to wavespeeds (i.e., ) equal to flow speeds at the two channel boundaries.
Modes with n>0 are the symmetric brethren of those with n<0. The reason is that (15) remains unchanged after applying the spatial reflection operation about the point x=1/2together with the frequency reflection/shift . A given eigensolution with eigenfrequency generates another eigensolution with eigenfrequency such that is the mirror reflection of about x=1/2. In general, we see that the approximate theory, including the predicted eigenfrequencies, are good even for order 1 values of and , but the theory begins to breakdown when they exceed order 1 values (Fig. 1).
Figure 1: Results from normal mode theory in which and q=3/2 and . a) The exact eigenfunctions for depicted for the first three overtones each at two values of . The eigensolutions developed via the asymptotic theory lie exactly on the eigenfunctions for where . b) The frequencies as a function of for the first three overtones are shown (open circles). The predicted frequencies using the asymptotic theory are shown with solid lines. The values corresponding to the eigenfunctions depicted in the left panel are designated with diamonds. Note that only the positive overtones (i.e., n>0) are shown here - the negative overtones are reflection symmetric about x=1/2 as discussed in the text. | |
Open with DEXTER |
With respect to the question of dynamics of wall-bounded rpCf, these results are quite suggestive. Rincon et al. (2007) showed that Rayleigh-stable anticyclonic rpCf does not appear to show a transition into a turbulent state in the manner in which transition is observed in pCf flow (e.g., Waleffe 2003). The strategy in those problems is to identify steady nonlinear solutions of the flow, usually by some continuation method from known solutions in otherwise unstable parameter regimes, (e.g., Nagata 1990; or forcing methods, Waleffe 2003) and study their stability properties. If there are a collection of steady nonlinear structures, some of which unstable, then the transition into a turbulent state may be approached. These states get easily triggered (bypass mechanism) because small scale disturbances in pCf get amplified due to the lift-up effect, a mechanism by which axisymmetric disturbances result in algebraic instability (Ellingsen & Palm 1975). Rincon et al. (2007) however show that there are two problems with this approach when employed to identify structures of rpCf near the Rayleigh-stable line (e.g., near q=2 for ). First, steady structures could not be identified, and that the antilift-up effect, which is the lift-up effect analog on the Rayleigh line, induces radial/vertical velocity fields that do not undergo the sorts of secondary instabilities that are characteristic of the induced azimuthal velocity fields generated by the lift-up effect in pCf flow.
Since the linear analysis of the inviscid Rayleigh-stable rpCf shows that for a given and there are conditions where there exists a countable number of oscillating modes, we openly speculate here that one might consider identifying oscillating nonlinear structures instead of steady ones. Perhaps one strategy would be to (i) identify oscillating nonlinear states in the infinite Reynolds number case ( ); and then (ii) ratchet down Re and follow/study the structures and their activity until; (iii) either one identifies a Re of transition, if a unique one exists free of hysteretic effects or determine if, as Rincon et al. (2007) suggest, rpCf exhibits the character of a chaotic saddle. However we note that identifying oscillating (travelling wave) nonlinear solutions is difficult in plane Poiseuille flow (pPf) as Toh & Itano (2003) have recently shown. However, it might be easier to identify such structures in rpCf because of the presence there of linear oscillations whereas such oscillations are absent in linearized pCf.
We note that for a general collection of perturbations utilizing SBC, (13) says that . The fluctuations in the total energy relate directly to the lateral exchange of kinetic energy from radially neighboring copies of the shearing box, where azimuthal velocities enter or exit boosted by the difference in the background shear velocities between the two boundaries. The quality of this exchange is a result of the type of correlations imposed upon the disturbances. For SBC, it means that these correlations are set a priori by the periodic length scale of disturbances which is set to be on the scale of the box L (e.g. see Regev & Umurhan 2008). By contrast, the imposition of no-normal flow boundary conditions removes this source of energy into and out of the domain for problems invoking CBC; in that case .
We feel that these observations about the differences between CBC and SBC, together with recent results, raises more questions that need to be resolved. For instance, linear disturbances of the CBC type are understood not to result in decay because of reflections off of the channel walls. SBC conditions, on the other hand, support modes that can result in fluctuations of the total energy depending on the correlations that exist across the radial boundaries. It is interesting, therefore, that Lesur & Longaretti (2005) demonstrate the persistence of activity for simulations of the SSB (slightly below the Rayleigh line) using SBC, while they do not report similar persistent activity for CBC. On the other hand, there is in Shen et al. (2006) a recent demonstration of the decaying quality of disturbances in the SBC for flows run with q=3/2^{}. Equally intriguing is that Lithwick (2007) shows that a single vortex may persist for a long time in a three-dimensional simulation of the SSB for q=3/2 using SBC which involves the complex interplay of three-wave mode interactions. While CBC conserves its total energy no activity has yet been observed under Keplerian-like conditions while for SBC, which does not conserve its , sustained activity seems to occur in some cases while decay characterizes others.
Given the preceeding, important questions to be asked:
Acknowledgements
We thank the anonymous referee who gave us critical comments which helped us focus the presentation of this work and pointed us to earlier references. This work was partially supported by BSF grant number 0603414082.