We studied dynamical balances in magnetorotational instability (MRI) turbulence with a net vertical field in the shearing box model of disks. Analyzing the turbulence dynamics in Fourier (${\bf k}$-)space, we identified three types of active modes that define turbulence characteristics. These modes have lengths similar to the box size, i.e., lie in the small wavenumber region in Fourier space labeled the vital area and are: (i) the channel mode - uniform in the disk plane with the smallest vertical wavenumber,(ii) the zonal flow mode - azimuthally and vertically uniform with the smallest radial wavenumber and (iii) the rest modes. The rest modes comprise those harmonics in the vital area whose energies reach more than $50 \%$ of the maximum spectral energy. The rest modes individually are not so significant compared to the channel and zonal flow modes, however, the combined action of their multitude is dominant over these two modes. These three mode types are governed by interplay of the linear and nonlinear processes, leading to their interdependent dynamics. The linear processes consist in disk flow nonmodality-modified classical MRI with a net vertical field. The main nonlinear process is transfer of modes over wavevector angles in Fourier space - the transverse cascade. The channel mode exhibits episodic bursts supplied by linear MRI growth, while the nonlinear processes mostly oppose this, draining the channel energy and redistributing it to the rest modes. As for the zonal flow, it does not have a linear source and is fed by nonlinear interactions of the rest modes.

We investigate magnetohydrodynamic turbulence driven by the magnetorotational instability (MRI) in Keplerian disks with a nonzero net azimuthal magnetic field using shearing box simulations. As distinct from most previous studies, we analyze turbulence dynamics in Fourier (${\bf k}$-) space to understand its sustenance. The linear growth of MRI with azimuthal field has a transient character and is anisotropic in Fourier space, leading to anisotropy of nonlinear processes in Fourier space. As a result, the main nonlinear process appears to be a new type of angular redistribution of modes in Fourier space -- the \emph{nonlinear transverse cascade} -- rather than usual direct/inverse cascade. We demonstrate that the turbulence is sustained by interplay of the linear transient growth of MRI (which is the only energy supply for the turbulence) and the transverse cascade. These two processes operate at large length scales, comparable to box size and the corresponding small wavenumber area, called \emph{vital area} in Fourier space is crucial for the sustenance, while outside the vital area direct cascade dominates. The interplay of the linear and nonlinear processes in Fourier space is generally too intertwined for a vivid schematization. Nevertheless, we reveal the \emph{basic subcycle} of the sustenance that clearly shows synergy of these processes in the self-organization of the magnetized flow system. This synergy is quite robust and persists for the considered different aspect ratios of the simulation boxes. The spectral characteristics of the dynamical processes in these boxes are qualitatively similar, indicating the universality of the sustenance mechanism of the MRI-turbulence.

To understand the self-sustenance of subcritical turbulence in spectrally stable shear flows, we performed direct numerical simulations of homogeneous shear turbulence for different aspect ratios of the flow domain and analyzed the dynamical processes in Fourier space. There are no exponentially growing modes in such flows and the turbulence is energetically supported only by the linear growth of perturbation harmonics due to the shear flow non-normality. This non-normality-induced, or nonmodal growth is anisotropic in spectral space, which, in turn, leads to anisotropy of nonlinear processes in this space. As a result, a transverse (angular) redistribution of harmonics in Fourier space appears to be the main nonlinear process in these flows, rather than direct or inverse cascades. We refer to this type of nonlinear redistribution as the nonlinear transverse cascade. It is demonstrated that the turbulence is sustained by a subtle interplay between the linear nonmodal growth and the nonlinear transverse cascade that exemplifies a well-known bypass scenario of subcritical turbulence. These two basic processes mainly operate at large length scales, comparable to the domain size. Therefore, this central, small wave number area of Fourier space is crucial in the self-sustenance; we defined its size and labeled it as the vital area of turbulence. Outside the vital area, the nonmodal growth and the transverse cascade are of secondary importance. Although the cascades and the self-sustaining process of turbulence are qualitatively the same at different aspect ratios, the number of harmonics actively participating in this process varies, but always remains quite large. This implies that the self-sustenance of subcritical turbulence cannot be described by low-order models.

We find and investigate via numerical simulations self-sustained two-dimensional turbulence in a magnetohydrodynamic flow with a maximally simple configuration: plane, noninflectional (with a constant shear of velocity) and threaded by a parallel uniform background magnetic field. This flow is spectrally stable, so the turbulence is subcritical by nature and hence it can be energetically supported just by transient growth mechanism due to shear flow nonnormality. This mechanism appears to be essentially anisotropic in spectral (wavenumber) plane and operates mainly for spatial Fourier harmonics with streamwise wavenumbers less than a ratio of flow shear to the Alfv\'{e}n speed, $k_y < S/u_A$ (i.e., the Alfv\'{e}n frequency is lower than the shear rate). We focused on the analysis of the character of nonlinear processes and underlying self-sustaining scheme of the turbulence, i.e., on the interplay between linear transient growth and nonlinear processes, in spectral plane. Our study, being concerned with a new type of the energy-injecting process for turbulence -- the transient growth, represents an alternative to the main trends of MHD turbulence research. We find similarity of the nonlinear dynamics to the related dynamics in hydrodynamic flows -- to the \emph{bypass} concept of subcritical turbulence. The essence of the analyzed nonlinear MHD processes appears to be a transverse redistribution of kinetic and magnetic spectral energies in wavenumber plane [as occurs in the related hydrodynamic flow, see Horton et al., Phys. Rev. E {\bf 81}, 066304 (2010)] and differs fundamentally from the existing concepts of (anisotropic direct and inverse) cascade processes in MHD shear flows.

We investigate linear dynamics of non-axisymmetric perturbations in incompressible, vertically stratified Keplerian discs with a weak vertical magnetic field in the shearing box approximation. Perturbations are decomposed into shearing waves whose evolution is followed via numerical integration of the linearized ideal MHD equations. There are two basic modes in the system -- inertia-gravity waves and magnetic mode, which displays the magnetorotational instability (MRI). As distinct from previous studies, we introduce `eigenvariables' characterizing each (counter-propagating) component of the inertia-gravity and magnetic modes, which are governed by a set of four first order coupled ordinary differential equations. This allowed us to identify a new process of linear coupling of the two above non-axisymmetric modes due to the disc's differential rotation. We did a comparative analysis of the dynamics of non-axisymmetric and axisymmetric magnetic mode perturbations. It is shown that the growth of optimal and close-to-optimal non-axisymmetric harmonics of this mode, having transient nature, can prevail over the exponential growth of axisymmetric ones (i.e., over the axisymmetric MRI) during dynamical time. A possible implication of this result for axisymmetric channel solutions is discussed. Specifically, the formation of the channel may be affected/impeded by non-axisymmetric modes already at the early linear stage leading to its untimely disruption -- the outcome strongly depends on the amplitude and spectrum of initial perturbation. So, this competition may result in an uncertainty in the magnetic mode's non-linear dynamics. It is also shown that a maximum non-axisymmetric growth is at vertical wavelengths close to the scale-height for which compressibility effects are important. This indirectly suggests that compressibility plays a role in the dynamics of the non-axisymmetric MRI.

We investigate mode coupling in a two dimensional compressible disc with radial stratification and differential rotation. We employ the global radial scaling of linear perturbations and study the linear modes in the local shearing sheet approximation. We employ a three-mode formalism and study the vorticity (W), entropy (S) and compressional (P) modes and their coupling properties. The system exhibits asymmetric three-mode coupling: these include mutual coupling of S and P-modes, S and W-modes, and asymmetric coupling between the W and P-modes. P-mode perturbations are able to generate potential vorticity through indirect three-mode coupling. This process indicates that compressional perturbations can lead to the development of vortical structures and influence the dynamics of radially stratified hydrodynamic accretion and protoplanetary discs.

Aims. Qualitative analysis of key (but yet unappreciated) linear phenomena in stratified hydrodynamic Keplerian flows: (i) the occurrence of a vortex mode, as a consequence of strato-rotational balance, with its transient dynamics; (ii) the generation of spiral-density waves (also called inertia-gravity or $g\Omega$ waves) by the vortex mode through linear mode coupling in shear flows. Methods. Non-modal analysis of linearized Boussinesq equations written in the shearing sheet approximation of accretion disk flows. Results. It is shown that the combined action of rotation and stratification introduces a new degree of freedom -- vortex mode perturbation -- which is linearly coupled with the spiral-density waves. These two modes are jointly able to extract energy from the background flow and they govern the disk dynamics in the small-scale range. The transient behavior of these modes is determined by the non-normality of the Keplerian shear flow. Tightly leading vortex mode perturbations undergo substantial transient growth, then, becoming trailing, inevitably generate trailing spiral-density waves by linear mode coupling. This course of events -- transient growth plus coupling -- is particularly pronounced for perturbation harmonics with comparable azimuthal and vertical scales and it renders the energy dynamics similar to the 3D unbounded plane Couette flow case. Conclusions. Our investigation strongly suggests that the so-called bypass concept of turbulence, which has been recently developed by the hydrodynamic community for spectrally stable shear flows, can also be applied to Keplerian disks. This conjecture may be confirmed by appropriate numerical simulations that take in account the vertical stratification and consequent mode coupling in the high Reynolds number regime.

Flow nonnormality induced linear transient phenomena in thin self-gravitating astrophysical discs are studied in the shearing sheet approximation. The considered system includes two modes of perturbations: vortex and (spiral density) wave. It is shown that self-gravity considerably alters the vortex mode dynamics -- its transient (swing) growth may be several orders of magnitude stronger than in the non-self-gravitating case and 2-3 times larger than the transient growth of the wave mode. Based on this finding, we comment on the role of vortex mode perturbations in a gravitoturbulent state. Also described is the linear coupling of the perturbation modes, caused by the differential character of disc rotation. The coupling is asymmetric -- vortex mode perturbations are able to excite wave mode ones, but not vice versa. This asymmetric coupling lends additional significance to the vortex mode as a participant in spiral density waves and shocks manifestations in astrophysical discs.

We investigate the stability, nonlinear development and equilibrium structure of vortices in a background shearing Keplerian flow. We make use of high-resolution global two-dimensional compressible hydrodynamic simulations. We introduce the concept of nonlinear adjustment to describe the transition of unbalanced vortical fields to a long-lived configuration. We discuss the conditions under which vortical perturbations evolve into long-lived persistent structures and we describe the properties of these equilibrium vortices. The properties of equilibrium vortices appear to be independent from the initial conditions and depend only on the local disk parameters. In particular we find that the ratio of the vortex size to the local disk scale height increases with the decrease of the sound speed, reaching values well above the unity. The process of spiral density wave generation by the vortex, discussed in our previous work, appear to maintain its efficiency also at nonlinear amplitudes and we observe the formation of spiral shocks attached to the vortex. The shocks may have important consequences on the long term vortex evolution and possibly on the global disk dynamics. Our study strengthens the arguments in favor of anticyclonic vortices as the candidates for the promotion of planetary formation. Hydrodynamic shocks that are an intrinsic property of persistent vortices in compressible Keplerian flows are an important contributor to the overall balance. These shocks support vortices against viscous dissipation by generating local potential vorticity and should be responsible for the eventual fate of the persistent anticyclonic vortices. Numerical codes have be able to resolve shock waves to describe the vortex dynamics correctly.

We analyze the scepticism on the hydrodynamic turbulence in Keplerian astrophysical disks expressed in Balbus and Hawley 2006 and show the failure of arguments of the paper.

Special features of magnetohydrodynamic waves linear dynamics in smooth shear flows are studied. Quantitative asymptotic and numerical analysis are performed for wide range of system parameters when basic flow has constant shear of velocity and uniform magnetic field is parallel to the basic flow. The special features consist of magnetohydrodynamic wave mutual transformation and over-reflection phenomena. The transformation takes place for arbitrary shear rates and involves all magnetohydrodynamic wave modes. While the over-reflection occurs only for slow magnetosonic and Alfv\'en waves at high shear rates. Studied phenomena should be decisive in the elaboration of the self-sustaining model of magnetohydrodynamic turbulence in the shear flows.

Special features of surface gravity waves in deep fluid flow with constant vertical shear of velocity is studied. It is found that the mean flow velocity shear leads to non-trivial modification of surface gravity wave modes dispersive characteristics. Moreover, the shear induces generation of surface gravity waves by internal vortex mode perturbations. The performed analytical and numerical study provides, that surface gravity waves are effectively generated by the internal perturbations at high shear rates. The generation is different for the waves propagating in the different directions. Generation of surface gravity waves propagating along the main flow considerably exceeds the generation of surface gravity waves in the opposite direction for relatively small shear rates, whereas the later wave is generated more effectively for the high shear rates. From the mathematical point of view the wave generation is caused by non self-adjointness of the linear operators that describe the shear flow.

We perform a detailed analytical and numerical study of the dynamics of perturbations (vortex/aperiodic mode, Rossby and spiral-density waves) in 2D compressible disks with a Keplerian law of rotation. We draw attention to the process of spiral-density wave generation from vortices, discussing, in particular, the initial, most peculiar stages of wave emission. We show that the linear phenomenon of wave generation by vortices in smooth (without inflection points) shear flows found by using the so-called non-modal approach, is directly applicable to the present case. After an analytical non-modal description of the physics and characteristics of the spiral-density wave generation/propagation in the local shearing-sheet model, we follow the process of wave generation by small amplitude coherent circular vortex structures, by direct global numerical simulation, describing the main features of the generated waves.

This paper deals with the problem of hydrodynamic shear turbulence in non-magnetized Keplerian disks. We wish to draw attention to a route to hydrodynamic turbulence which seems to be little known by the astrophysical community, but which has been intensively discussed among fluid dynamicists during the past decade. In this so-called `bypass' concept for the onset of turbulence, perturbations undergo a transient growth, and they may reach an amplitude that is sufficiently large to allow positive feedback through nonlinear interactions. This transient growth is linear in nature, and thus it differs in principle from the well-known nonlinear instability. We describe the type of perturbations that according to this process are the most likely to lead to turbulence, namely non-axisymmetric vortex mode perturbations in the two dimensional limit. We show that the apparently inhibiting action of the Coriolis force on the dynamics of such vortical perturbations is substantially diminished due to the pressure perturbations, contrary to current opinion. We stress the similarity of the turbulent processes in Keplerian disks and in Cartesian flows and conclude that the prevalent skepticism of the astrophysical community on the occurrence of hydrodynamic shear turbulence in such disks is not founded.

In this paper we report on the nonresonant conversion of convectively unstable linear gravity modes into acoustic oscillation modes in shear flows. The convectively unstable linear gravity modes can excite acoustic modes with similar wave-numbers. The frequencies of the excited oscillations may be qualitatively higher than the temporal variation scales of the source flow, while the frequency spectra of the generated oscillations should be intrinsically correlated to the velocity field of the source flow. We anticipate that this nonresonant phenomenon can significantly contribute to the production of sound waves in the solar convection zone.