In order to study chromospheric magnetosonic wave propagation including, for the first time, the effects of ion-neutral interactions in the partially ionized solar chromosphere, we have developed a new multi-fluid computational model, accounting for ionization and recombination reactions in gravitationally stratified magnetized collisional media. The two-fluid model used in our 2D numerical simulations treats neutrals as a separate fluid and considers charged species (electrons and ions) within the resistive MHD approach with Coulomb collisions and anisotropic heat flux determined by Braginskii's transport coefficients. The electromagnetic fields are evolved according to the full Maxwell equations and the solenoidality of the magnetic field is enforced with a hyperbolic divergence cleaning scheme. The initial density and temperature profiles are similar to VAL III chromospheric model in which dynamical, thermal and chemical equilibrium are considered to ensure comparison to existing MHD models and avoid artificial numerical heating. In this initial setup we include simple homogeneous flux tube magnetic field configuration and an external photospheric velocity driver to simulate the propagation of MHD waves in the partially ionized reactive chromosphere. In particular, we investigate the loss of chemical equilibrium and the plasma heating related to the steepening of fast magnetosonic wave fronts in the gravitationally stratified medium.

We have derived the full set of MHD equations for incompressible shear flow of a magnetized fluid and considered their solution in the wave-vector space. The linearized equations give the famous amplification of slow magnetosonic waves and describe the magnetorotational instability. The nonlinear terms in our analysis are responsible for the creation of turbulence and self-sustained spectral density of the MHD (Alfven and pseudo-Alfven) waves. Perspectives for numerical simulations of weak turbulence and calculation of the effective viscosity of accretion disks are shortly discussed in k-space.

Propagation of three dimensional magnetosonic waves is considered for a homogeneous shear flow of an incompressible fluid. The analytical solutions for all magnetohydrodynamic variables are presented by confluent Heun functions. The problem is reduced to finding a solution of an effective Schroedinger equation. The amplification of slow magnetosonic waves is analyzed in great details. A simple formula for the amplification coefficient is derived. The velocity shear primarily affects the incompressible limit of slow magnetosonic waves. The amplification is very strong for slow magnetosonic waves in the long-wavelength limit. It is demonstrated that the amplification of those waves leads to amplification of turbulence. The phenomenology of Shakura-Sunyaev for the friction in accretion disks is derived in the framework of the Kolmogorov turbulence. The presented findings may be the key to explaining the anomalous plasma heating responsible for the luminosity of quasars. It is suggested that wave amplification is the keystone of the self-sustained turbulence in accretion disks.

The set of equations for magnetohydrodynamic (MHD) waves in a shear flow is consecutively derived. The proposed scenario involves the presence of a self-sustained turbulence and magnetic field. In the framework of Langevin--Burgers approach the influence of the turbulence is described by an additional external random force in the MHD system. Kinetic equation for the spectral density of the slow magnetosonic (Alfvenic) mode is derived in the short wavelength (WKB) approximation. The results show a pressing need for conduction of numerical Monte Carlo (MC) simulations with a random driver to take into account the influence of the long wavelength modes and to give a more precise analytical assessment of the short ones. Realistic MC calculations for the heating rate and shear stress tensor should give an answer to the perplexing problem for the missing viscosity in accretion disks and reveal why the quasars are the most powerful sources of light in the universe. The planned MC calculations can be incorporated in global models for accretion disks and also in all other physical conditions where there is a shear flow in a magnetized turbulent plasma. It is supposed that the heating mechanism by Alfven waves absorption is common for many kinds of space plasmas from solar corona to active galactic nuclei and the solution of these longstanding puzzles deserves active interdisciplinary research. The work is illustrated by numerical calculations and by exact solutions for the time dependence of the magnetic field given by the Heun function.

Absorbtion of Alfven waves is considered to be the main mechanism of heating in the solar corona. It is concluded that the sharp increase of the plasma temperature by two orders of magnitude is related to a self-induced opacity with respect to Alfven waves. The maximal frequency for propagation of Alfven waves is determined by the strongly temperature dependent kinematic viscosity. In such a way the temperature jump is due to absorption of high frequency Alfven waves in a narrow layer above the solar surface. It is calculated that the power per unit area dissipated in this layer due to damping of Alfven waves blows up the plasma and gives birth to the solar wind. A model short wave-length (WKB) evaluation takes into account the 1/f^2 frequency dependance of the transversal magnetic field and velocity spectral densities. Such spectral densities agree with old magnetometric data taken by Voyager 1 and recent theoretical calculations in the framework of Langevin-Burgers MHD. The present theory predicts existence of intensive high frequency MHD Alfven waves in the cold layer beneath the corona. It is briefly discussed how this statement can be checked experimentally. It is demonstrated that the magnitude of the Alfven waves generating random noise and the solar wind velocity can be expressed only in terms of satellite experimental data. It is advocated that investigation of properties of the solar surface as a random driver by optical methods is an important task for future solar physics. Jets of accretion disks are speculated as a special case of the wind from magnetized turbulent plasma.

The heating of magnetized plasma by propagation of Alfven waves is calculated as a function of the magnetic field spectral density. The results can be applied to evaluate the heating power of the solar corona at known data from satellites' magnetometers. This heating rate can be incorporated in global models for heating of the solar corona and creation of the solar wind. The final formula for the heating power is illustrated with a model spectral density of the magnetic field obtained by analysis of the Voyager 1 mission results. The influence of high frequency dissipative modes is also taken into account and it is concluded that for evaluation of the total coronal heating it is necessary to know the spectral density of the fluctuating component of the magnetic field up to the frequency of electron-proton collisions.

Heating of magnetized turbulent plasma is calculated in the framework of Burgers turbulence [A.M. Polyakov, Phys. Rev. E. 52, 6183 (1995)]. Explicit formula for the energy flux of Alfven waves along the magnetic field is presented. The Alfven waves are considered as intermediary between the turbulent energy and the heat. The derived results are related to a wave channel of heating of the solar corona. If we incorporate amplification of Alfven waves by shear flow the suggested model of heating can be applied to analysis of the missing viscosity of accretion discs and to reveal why the quasars are the most powerful sources of light in the universe. We suppose that the Langevin-Burgers approach to turbulence we have applied in the current work can be also helpful for other systems where we have intensive interaction between a stochastic turbulent system and waves and can be used in many multidisciplinary researches in hydrodynamics and MHD.