With the use of our JOANNA code, which solves radiative equations for ion + electron and neutral fluids, we perform realistic 2.5D numerical simulations of plasma outflows associated with the solar granulation. These outflows exhibit physical quantities consistent to the order of magnitude with the observational findings for mass and energy losses in the upper chromosphere, transition region and inner corona, and they may originate the fast solar wind.

Using high-resolution numerical simulations we investigate the plasma heating driven by periodic two-fluid acoustic waves that originate at the bottom of the photosphere and propagate into the gravitationally stratified and partially ionized solar atmosphere. We consider ions+electrons and neutrals as separate fluids that interact between themselves via collision forces. The latter play an important role in the chromosphere, leading to significant damping of short-period waves. Long-period waves do not essentially alter the photospheric temperatures, but they exhibit the capability of depositing a part of their energy in the chromosphere. This results in up about a five times increase of ion temperature that takes place there on a time-scale of a few minutes. The most effective heating corresponds to waveperiods within the range of about 30-200 s with a peak value located at 80 s. However, we conclude that for the amplitude of the driver chosen to be equal to 0.1 km s$^{-1}$, this heating is too low to balance the radiative losses in the chromosphere.

We aim to study the formation and evolution of solar spicules by means of numerical simulations of the solar atmosphere. With the use of newly developed JOANNA code, we numerically solve two-fluid (for ions + electrons and neutrals) equations in 2D Cartesian geometry. We follow the evolution of a spicule triggered by the time-dependent signal in ion and neutral components of gas pressure launched in the upper chromosphere. We use the potential magnetic field, which evolves self-consistently, but mainly plays a passive role in the dynamics. Our numerical results reveal that the signal is steepened into a shock that propagates upward into the corona. The chromospheric cold and dense plasma lags behind this shock and rises into the corona with a mean speed of 20-25 km s$^{-1}$. The formed spicule exhibits the upflow/downfall of plasma during its total lifetime of around 3-4 minutes, and it follows the typical characteristics of a classical spicule, which is modeled by magnetohydrodynamics. The simulated spicule consists of a dense and cold core that is dominated by neutrals. The general dynamics of ion and neutral spicules are very similar to each other. Minor differences in those dynamics result in different widths of both spicules with increasing rarefaction of the ion spicule in time.

We investigate numerically Alfv\'en waves propagating along an axisymmetric and non-isothermal solar flux tube embedded in the solar atmosphere. The tube magnetic field is current-free and diverges with height, and the waves are excited by a periodic driver along the tube magnetic field lines. The main results are that the two wave variables, the velocity and magnetic field perturbations in the azimuthal direction, behave differently as a result of gradients of physical parameters along the tube. To explain these differences in the wave behavior, the time evolution of the wave variables and the resulting cutoff period for each wave variable are calculated, and used to determine regions in the solar chromosphere where strong wave reflection may occur.