Linear dynamics of the solar convection zone: excitation of waves in unstably stratified shear flows

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.Comment: 8 pages. To appear in the proceedings of the conference "Waves in Dusty, Solar and Space Plasmas", Leuven, Belgium 21-26 May 200

Fast magnetohydrodynamic oscillation of longitudinally inhomogeneous prominence threads: an analogue with quantum harmonic oscillator

Previous works indicate that the frequency ratio of second and first harmonics of kink oscillations has tendency towards 3 in the case of prominence threads. We aim to study the magnetohydrodynamic oscillations of longitudinally inhomogeneous prominence threads and to shed light on the problem of frequency ratio. Classical Sturm--Liouville problem is used for the threads with longitudinally inhomogeneous plasma density. We show that the spatial variation of total pressure perturbations along the thread is governed by the stationary Schr\"{o}dinger equation, where the longitudinal inhomogeneity of plasma density stands for the potential energy. Consequently, the equation has bounded solutions in terms of Hermite polynomials. Boundary conditions at the thread surface lead to transcendental dispersion equation with Bessel functions. Thin flux tube approximation of the dispersion equation shows that the frequency of kink waves is proportional to the expression \alpha(2n+1), where \alpha is the density inhomogeneity parameter and n is the longitudinal mode number. Consequently, the ratio of the frequencies of second and first harmonics tends to 3 in prominence threads. Numerical solution of the dispersion equation shows that the ratio only slightly decreases for thicker tubes in the case of smaller longitudinal inhomogeneity of external density, therefore the thin flux tube limit is a good approximation for prominence oscillations. However, stronger longitudinal inhomogeneity of external density may lead to the significant shift of frequency ratio for wider tubes and therefore the thin tube approximation may fail. The tendency of frequency ratio of second and first harmonics towards 3 in prominence threads is explained by the analogy of the oscillations with quantum harmonic oscillator, where the density inhomogeneity of the threads plays a role of potential energy.Comment: 8 pages, 7 figures (accepted in A&A

Linear coupling of modes in 2D radially stratified astrophysical discs

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.Comment: 10 pages, 10 figures, MNRAS (accepted

Hydrodynamic stability and mode coupling in Keplerian flows: local strato-rotational analysis

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.Comment: A&A (accepted

Numerical simulations of the decay of primordial magnetic turbulence

We perform direct numerical simulations of forced and freely decaying 3D magnetohydrodynamic turbulence in order to model magnetic field evolution during cosmological phase transitions in the early Universe. Our approach assumes the existence of a magnetic field generated either by a process during inflation or shortly thereafter, or by bubble collisions during a phase transition. We show that the final configuration of the magnetic field depends on the initial conditions, while the velocity field is nearly independent of initial conditions.Comment: 10 pages, 6 figures, references added, PRD accepte

Stability and nonlinear adjustment of vortices in Keplerian flows

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.Comment: 12 pages, 10 figure