We investigate the dynamics of linear perturbations in Keplerian flow under external stochastic force. To abstract from the details of flow structure and boundary conditions, we consider the problem in the shearing box approximation. An external force is assumed to have zero mean, even so, induced perturbations form a steady-state, which provides angular momentum transfer to the periphery of the flow. The most effective scenario is based on the transient amplification of induced vortices with the following emission of shearing sound wave, wherein the maximum of the flux linearly depends on Reynolds number. Thus such a mechanism is significant for astrophysical flows, for which enormous Reynolds numbers are typical. At the same time, addressing the problem analytically, we found that for incompressible fluid in the shearing box approximation stochastic forcing does not lead to average angular momentum transfer. Thus the compressibility of the fluid plays an important role here, and one cannot neglect it.

Enhanced angular momentum transfer through the boundary layer near the surface of weakly magnetised accreting star is required in order to explain the observed accretion timescales in low-mass X-ray binaries, cataclysmic variables or young stars with massive protoplanetary discs. Accretion disc boundary layer is locally represented by incompressible homogeneous and boundless flow of the cyclonic type, which is linearly stable. Its non-linear instability at the shear rates of order of the rotational frequency remains an issue. We argue that hydrodynamical subcritical turbulence in such a flow is sustained by the non-linear feedback from essentially three-dimensional vortices, which are generated by quasi-two-dimensional trailing shearing spirals grown to high amplitude via the swing amplification. We refer to those three-dimensional vortices as cross-rolls, since they are aligned in the shearwise direction in contrast to streamwise rolls generated by the anti-lift-up mechanism in rotating shear flow on the Rayleigh line. Transient growth of cross-rolls is studied analytically and further confronted with direct numerical simulations (DNS) of dynamics of non-linear erturbations. DNS performed in a tall box show that transition Reynolds number $R_T$ as function of shear rate accords with the line of constant maximum transient growth of cross-rolls. The transition in the tall box has been observed until the shear rate three times higher than the rotational frequency, when $R_T\sim 50000$. Assuming that the cross-rolls are also responsible for turbulence in the Keplerian flow, we estimate $R_T\lesssim 10^8$ in this case. The most favourable shear rate to test the existence of turbulence in the quasi-Keplerian regime may be sub-Keplerian and equal approximately to $1/2$.

In this paper analytical method for constructing baroclinic flow is presented. Presented method allows the construction of baroclinic tori with significant vertical shear, which are stable by the Solberg--H{\o}iland criteria. The presented method can be useful in the investigations of the influence of background structure on vertical-shear- and entropy-driven instabilities.

Turbulent state of spectrally stable shear flows may be developed and sustained according to the bypass scenario of transition. If it works in non-magnetised boundless and homogeneous quasi-Keplerian flow, transiently growing shearing vortices should supply turbulence with energy. Employing the large shearing box approximation, as well as a set of global disc models, we study the optimal growth of the shearing vortices in such a flow in the whole range of azimuthal length-scales, $\lambda_y$, as compared to the flow scale-height, $H$. It is shown that with the account of the viscosity the highest possible amplification of shearing vortices, $G_{max}$, attains maximum at $\lambda_y\lesssim H$ and declines towards both the large scales $\lambda_y\gg H$ and the small scales $\lambda_y\ll H$ in a good agreement with analytical estimations based on balanced solutions. We pay main attention to the large-scale vortices $\lambda_y\gg H$, which produce $G_{max}\propto (\Omega/\kappa)^4$, where $\Omega$ and $\kappa$ denote local rotational and epicyclic frequencies, respectively. It is demonstrated that the large-scale vortices acquire high density perturbation as they approach the instant of swing. At the same time, their growth is not affected by bulk viscosity. We check that $G_{max}$ obtained globally is comparable to its local counterpart and the shape and localisation of global optimal vortices can be explained in terms of the local approach. The obtained results allow us to suggest that the critical Reynolds number of subcritical transition to turbulence in quasi-Keplerian flow, as well as the corresponding turbulent effective azimuthal stress should substantially depend on shear rate.

The mechanism of transition from laminar state to turbulent state in Keplerian disks is still unknown. The most popular version today is generation of turbulence due to magnetorotational instability (MRI). However magnetohydrodynamic simulations give the value of Shakura-Sunyaev parameter more then an order of magnitude smaller rather than that found from observations. One way to solve this problem is the existence of an alternative or additional mechanism for generating turbulence. It can be the bypass mechanism, which is responsible for transition to turbulence in Couette and Poiseuille flows. This mechanism is based on the transient growth of linear perturbations in the flow with the subsequent transition to the nonlinear stage. In order to clarify the role of this mechanism in astrophysical disks first of all it is necessary to calculate the maximal possible growth factor of linear perturbations in the flow. In this paper the results of such calculations are presented for perturbations on different scales compared with the disk thickness. Qualitative description of mechanisms responsible for the growth will also be presented. It was found that the most rapidly growing shear harmonics have azimuthal wavelength of the order of the disk thickness. In addition, their initial form is always similar to the vortex perturbations with the same potential vorticity. It was shown that the vortices with azimuthal wavelength more than an order of magnitude in excess of the disc thickness, are still able to grow by dozens of times.

This paper reviews some aspects of one of the major unsolved problems in understanding astrophysical (in particular, accretion) disks: whether the disk interiors may be effectively viscous in spite of the absence of marnetorotational instability? In this case a rotational homogeneous inviscid flow with a Keplerian angular velocity profile is spectrally stable, making the transient growth of perturbations a candidate mechanism for energy transfer from the regular motion to perturbations. Transient perturbations differ qualitatively from perturbation modes and can grow substantially in shear flows due to the nonnormality of their dynamical evolution operator. Since the eigenvectors of this operator, alias perturbation modes, are mutually nonorthogonal, they can mutually interfere, resulting in the transient growth of their linear combinations. Physically, a growing transient perturbation is a leading spiral whose branches are shrunk as a result of the differential rotation of the flow. This paper discusses in detail the transient growth of vortex shear harmonics in the spatially local limit as well as methods for identifying the optimal (fastest growth) perturbations. Special attention is given to obtaining such solutions variationally, by integrating the direct and adjoint equations forward and backward in time, respectively. The material is presented in a newcomer-friendly style.

We study linear transient dynamics in a thin Keplerian disc employing a method based on variational formulation of optimisation problem. It is shown that in a shearing sheet approximation due to a prominent excitation of density waves by vortices the most rapidly growing shearing harmonic has azimuthal wavelength, $\lambda_y$, of order of the disc thickness, $H$, and its initial shape is always nearly identical to a vortex having the same potential vorticity. Also, in the limit $\lambda_y\gg H$ the optimal growth $G\propto (\Omega/\kappa)^4$, where $\Omega$ and $\kappa$ stand for local rotational and epicyclic frequencies, respectively, what suggests that transient growth of large scale vortices can be much stronger in areas with non-Keplerian rotation, e.g. in the inner parts of relativistic discs around the black holes. We estimate that if disc is already in a turbulent state with effective viscosity given by the Shakura parameter $\alpha<1$, the considered large scale vortices with wavelengths $H/\alpha>\lambda_y>H$ have the most favourable conditions to be transiently amplified before they are damped. At the same time, turbulence is a natural source of the potential vorticity for this transient activity. We extend our study to a global spatial scale showing that global perturbations with azimuthal wavelengths more than an order of magnitude greater than the disc thickness still are able to attain the growth of dozens of times in a few Keplerian periods at the inner boundary of disc.

A thin gaseous disc with an almost keplerian angular velocity profile, bounded by a free surface and rotating around point-mass gravitating object is nearly spectrally stable. Despite that the substantial transient growth of linear perturbations measured by the evolution of their acoustic energy is possible. This fact is demonstrated for the simple model of a non-viscous polytropic thin disc of a finite radial size where the small adiabatic perturbations are considered as a linear combination of neutral modes with a corotational radius located beyond the outer boundary of the flow.