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 perform numerical simulations of impulsively generated magnetic swirls in an isolated flux tube which is rooted in the solar photosphere. These swirls are triggered by an initial pulse in a horizontal component of the velocity. The initial pulse is launched either: (a) centrally, within the localized magnetic flux tube; or (b) off-central, in the ambient medium. The evolution and dynamics of the flux tube is described by three-dimensional, ideal magnetohydrodynamic equations. These equations are numerically solved to reveal that in case (a) dipole-like swirls associated with the fast magnetoacoustic kink and $m=1$ Alfv\'en waves are generated. In case (b), the fast magnetoacoustic kink and $m=0$ Alfv\'en modes are excited. In both these cases, the excited fast magnetoacoustic kink and Alfv\'en waves consist of similar flow pattern and magnetic shells are also generated with clockwise and counter-clockwise rotating plasma within them, which can be the proxy of dipole-shaped chromospheric swirls. The complex dynamics of vortices and wave perturbations reveals the channelling of sufficient amount of energy to fulfill energy losses in the chromosphere ($\sim$ 10$^{4}$ W m$^{-1}$) and in the corona ($\sim$ 10$^{2}$ W m$^{-1}$). Some of these numerical findings are reminiscent of signatures in recent observational data.