We study the effects of heat flows and velocity shear on the parallel firehose instability in weakly collisional plasma flow. For this purpose we apply an anisotropic 16-moments MHD fluid closure model that takes into account the pressure and temperature anisotropy, as well as the effect of anisotropic heat flux. The linear stability analysis of the firehose modes is carried out in the incompressible limit, where the MHD flow is parallel to the background magnetic field, while the velocity is sheared in the direction transverse to the flow direction. It seems that an increase of the velocity shear parameter leads to higher growth rates of the firehose instability. The increase of the instability growth rate is most profound for perturbations with oblique wave-numbers $k_{\perp}/k_{\parallel} < 1$. The heat flux parameter introduces an asymmetry of the instability growth in the shear plane: perturbations with wave-vectors with a component in the direction of the velocity shear grow significantly stronger as compared to those with components in the opposite direction. We discuss the implications of the presented study on the observable features of the solar wind and possible measurements of local parameters of the solar wind based on the stability constraints set by the firehose instability.

We develop the theoretical model for the analytic description of hydrodynamic jets from protostellar disks employing the Beltrami-Bernoulli flow configuration of disk-jet structure. For this purpose we extend the standard turbulent viscosity prescription and derive several classes of analytic solutions using the flow parametrization in self-similar variables. Derived solutions describe the disk-jet structure, where for the first time jet properties are analytically linked with the properties of the accretion disk flow. The ratio of the jet ejection and disk accretion velocities is controlled by the turbulence parameter, while the ejection velocity increases with the decrease of local sound velocity and the jet launching radius. Derived solutions can be used to analyze the astrophysical jets from protostellar accretion disks and link the properties of outflows with the local observational properties of accretion disk flows.

We show that in decaying hydromagnetic turbulence with initial kinetic helicity, a weak magnetic field eventually becomes fully helical. The sign of magnetic helicity is opposite to that of the kinetic helicity - regardless of whether or not the initial magnetic field was helical. The magnetic field undergoes inverse cascading with the magnetic energy decaying approximately like t^{-1/2}. This is even slower than in the fully helical case, where it decays like t^{-2/3}. In this parameter range, the product of magnetic energy and correlation length raised to a certain power slightly larger than unity, is approximately constant. This scaling of magnetic energy persists over long time scales. At very late times and for domain sizes large enough to accommodate the growing spatial scales, we expect a cross-over to the t^{-2/3} decay law that is commonly observed for fully helical magnetic fields. Regardless of the presence or absence of initial kinetic helicity, the magnetic field experiences exponential growth during the first few turnover times, which is suggestive of small-scale dynamo action. Our results have applications to a wide range of experimental dynamos and astrophysical time-dependent plasmas, including primordial turbulence in the early universe.

We consider a scale-invariant helical magnetic field generated during inflation. We show that, if the mean magnetic helicity density of such a field is measured, it can be used to determine a lower bound on the duration of inflation. Upper bounds can be used to derive constraints on the minimal duration of inflation if one assumes that the magnetic field generated during inflation is helical. Using three-dimensional simulations, we show that an initially scale-invariant field develops, which is similar both with and without magnetic helicity. In the fully helical case, however, the magnetic field appears to have a more pronounced folded structure.

We present new simulations of decaying hydromagnetic turbulence for a relativistic equation of state relevant to the early universe. We compare helical and nonhelical cases either with kinetically or magnetically dominated initial fields. Both kinetic and magnetic initial helicities lead to maximally helical magnetic fields after some time, but with different temporal decay laws. Both are relevant to the early universe, although no mechanisms have yet been identified that produce magnetic helicity with strengths comparable to the big bang nucleosynthesis limit at scales comparable to the Hubble horizon at the electroweak phase transition. Nonhelical magnetically dominated fields could still produce picoGauss magnetic fields under most optimistic conditions. Only helical magnetic fields can potentially have nanoGauss strengths at scales up to 30 kpc today.

The linear stability of viscous Keplerian flow around a gravitating center is studied using the rheological granular fluid model. The linear rheological instability triggered by the interplay of the shear rheology and Keplerian differential rotation of incompressible dense granular fluids is found. Instability sets in granular fluids, where the viscosity parameter grows faster than the square of the local shear rate (strain rate) at constant pressure. Found instability can play a crucial role in the dynamics of dense planetary rings and granular flows in protoplanetary disks.

We study the evolution of primordial magnetic fields in an expanding cosmic plasma. For this purpose we present a comprehensive theoretical model to consider the evolution of MHD turbulence that can be used over a wide range of physical conditions, including cosmological and astrophysical applications. We model different types of decaying cosmic MHD turbulence in the expanding universe and characterize the large-scale magnetic fields in such a medium. Direct numerical simulations of freely decaying MHD turbulence are performed for different magnetogenesis scenarios: magnetic fields generated during cosmic inflation as well as electroweak and QCD phase transitions in the early universe. Magnetic fields and fluid motions are strongly coupled due to the high Reynolds number in the early universe. Hence, we abandon the simple adiabatic dilution model to estimate magnetic field amplitudes in the expanding universe and include turbulent mixing effects on the large-scale magnetic field evolution. Numerical simulations have been carried out for non-helical and helical magnetic field configurations. The numerical results show the possibility of inverse transfer of energy in magnetically dominated non-helical MHD turbulence. On the other hand, decay properties of helical turbulence depend on whether the turbulent magnetic field is in a weakly or a fully helical state. Our results show that primordial magnetic fields can be considered as a seed for the observed large-scale magnetic fields in galaxies and clusters. Bounds on the magnetic field strength are obtained and are consistent with the upper and lower limits set by observations of extragalactic magnetic fields.

We present a linear stability analysis of the perturbation modes in anisotropic MHD flows with velocity shear and strong magnetic field. Collisionless or weakly collisional plasma is described within the 16-momentum MHD fluid closure model, that takes into account not only the effect of pressure anisotropy, but also the effect of anisotropic heat fluxes. In this model the low frequency acoustic wave is revealed into a standard acoustic mode and higher frequency fast thermo-acoustic and lower frequency slow thermo-acoustic waves. It is shown that thermo-acoustic waves become unstable and grow exponentially when the heat flux parameter exceeds some critical value. It seems that velocity shear makes thermo-acoustic waves overstable even at subcritical heat flux parameters. Thus, when the effect of heat fluxes is not profound acoustic waves will grow due to the velocity shear, while at supercritical heat fluxes the flow reveals compressible thermal instability. Anisotropic thermal instability should be also important in astrophysical environments, where it will limit the maximal value of magnetic field that a low density ionized anisotropic flow can sustain.

In the presence of magnetic helicity, inverse transfer from small to large scales is well known in magnetohydrodynamic (MHD) turbulence and has applications in astrophysics, cosmology, and fusion plasmas. Using high resolution direct numerical simulations of magnetically dominated self-similarly decaying MHD turbulence, we report a similar inverse transfer even in the absence of magnetic helicity. We compute for the first time spectral energy transfer rates to show that this inverse transfer is about half as strong as with helicity, but in both cases the magnetic gain at large scales results from velocity at similar scales interacting with smaller-scale magnetic fields. This suggests that both inverse transfers are a consequence of a universal mechanisms for magnetically dominated turbulence. Possible explanations include inverse cascading of the mean squared vector potential associated with local near two-dimensionality and the shallower $k^2$ subinertial range spectrum of kinetic energy forcing the magnetic field with a $k^4$ subinertial range to attain larger-scale coherence. The inertial range shows a clear $k^{-2}$ spectrum and is the first example of fully isotropic magnetically dominated MHD turbulence exhibiting weak turbulence scaling.

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.

We study primordial magnetic field effects on the matter perturbations in the Universe. We assume magnetic field generation prior to the big bang nucleosynthesis (BBN), i.e. during the radiation dominated epoch of the Universe expansion, but do not limit analysis by considering a particular magnetogenesis scenario. Contrary to previous studies, we limit the total magnetic field energy density and not the smoothed amplitude of the magnetic field at large (order of 1 Mpc) scales. We review several cosmological signatures, such as halos abundance, thermal Sunyaev Zel'dovich (tSZ) effect, and Lyman-$\alpha$ data. For a cross check we compare our limits with that obtained through the CMB faraday rotation effect and BBN. The limits are ranging between 1.5 nG and 4.5 nG for $n_B \in (-3; -1.5)$.

We consider the evolution of primordial magnetic fields generated during cosmological, electroweak or QCD, phase transitions. We assume that the magnetic field generation can be described as an injection of magnetic energy to cosmological plasma at a given scale determined by the moment of magnetic field generation. A high Reynolds number ensures strong coupling between magnetic field and fluid motions. The subsequent evolution of the magnetic field is governed by decaying hydromagnetic turbulence. Both our numerical simulations and a phenomenological description allow us to recover "universal" laws for the decay of magnetic energy and the growth of magnetic correlation length in the turbulent (low viscosity) regime. In particular, we show that during the radiation dominated epoch, energy and correlation length of non-helical magnetic fields scale as conformal time to the powers -1/2 and +1/2, respectively. For helical magnetic fields, energy and correlation length scale as conformal time to the powers -1/3 and +2/3, respectively. The universal decay law of the magnetic field implies that the strength of magnetic field generated during the QCD phase transition could reach $\sim 10^{-9}$\,G with the present day correlation length $\sim 50$ kpc. The fields generated at the electroweak phase transition could be as strong as $\sim 10^{-10}$ G with correlation lengths reaching $\sim 0.3$\,kpc. These values of the magnetic fields are consistent with the lower bounds of the extragalactic magnetic fields.

We study the evolution of an inflation-generated magnetic field, due to its coupling to fluid motions, during cosmological phase transitions. We find that the magnetic field stays almost unchanged on large scales, while on small scales the spectrum is modified in such a way that power at small scales becomes progressively suppressed. We also show that the magnetic field generates turbulent motions in the initially turbulence-free plasma. On large scales, the slope of the resulting kinetic energy spectrum is consistent with that of white noise.

We study the evolution of QCD phase transition-generated magnetic fields in freely decaying MHD turbulence of the expanding Universe. We consider a magnetic field generation model that starts from basic non-perturbative QCD theory and predicts stochastic magnetic fields with an amplitude of the order of 0.02 $\mu$G and small magnetic helicity. We employ direct numerical simulations to model the MHD turbulence decay and identify two different regimes: "weakly helical" turbulence regime, when magnetic helicity increases during decay, and "fully helical" turbulence, when maximal magnetic helicity is reached and an inverse cascade develops. The results of our analysis show that in the most optimistic scenario the magnetic correlation length in the comoving frame can reach 10 kpc with the amplitude of the effective magnetic field being 0.007 nG. We demonstrate that the considered model of magneto-genesis can provide the seed magnetic field for galaxies and clusters.

We constrain a primordial magnetic field (PMF) generated during a phase transition (PT) using the big bang nucleosynthesis bound on the relativistic energy density. The amplitude of the PMF at large scales is determined by the shape of the PMF spectrum outside its maximal correlation length scale. Even if the amplitude of the PMF at 1 Mpc is small, PT-generated PMFs can leave observable signatures in the potentially detectable relic gravitational wave background if a large enough fraction ($1-10%$) of the thermal energy is converted into the PMF.

We study limits on a primordial magnetic field arising from cosmological data, including that from big bang nucleosynthesis, cosmic microwave background polarization plane Faraday rotation limits, and large-scale structure formation. We show that the physically-relevant quantity is the value of the effective magnetic field, and limits on it are independent of how the magnetic field was generated.

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.

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.

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.

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.