The solar chromosphere is very dynamic, due to the presence of large amplitude hydrodynamic waves. Their propagation is affected by NLTE radiative transport in strong spectral lines, which can in turn be used to diagnose the dynamics of the chromosphere. We give a basic introduction into the equations of NLTE radiation hydrodynamics and describe how they are solved in current numerical simulations. The comparison with observation shows that one-dimensional codes can describe strong brightenings quite well, but the overall chromospheric dynamics appears to be governed by three-dimensional shock propagation.

We aim to identify the physical processes which occur in the magnetic network of the chromosphere and which contribute to its dynamics and heating. Specifically, we study the propagation of transverse (kink) MHD waves which are impulsively excited in flux tubes through footpoint motions. When these waves travel upwards, they get partially converted to longitudinal waves through nonlinear effects (mode coupling). By solving the nonlinear, time-dependent MHD equations we find that significant longitudinal wave generation occurs in the photosphere typically for Mach numbers as low as 0.2 and that the onset of shock formation occurs at heights of about 600 km above the photospheric base. We also investigate the compressional heating due to longitudinal waves and the efficiency of mode coupling for various values of the plasma $\beta$, that parameterises the magnetic field strength in the network. We find that this efficiency is maximum for field strengths corresponding to $\beta\approx 0.2$, when the kink and tube wave speeds are almost identical. This can have interesting observational implications. Furthermore, we find that even when the two speeds are different, once shock formation occurs, the longitudinal and transverse shocks exhibit strong mode coupling.

We study the evolution of spectra of acoustic waves that are generated in the convection zone and propagate upward into the photosphere, where we compare the simulated acoustic spectra with the spectrum observed in an Fe I line. Although there is no pronounced 3 min component in the spectra generated in the convection zone, there are dominant 3 min features in the theoretical spectra, in agreement with the observed spectrum. We interpret the occurrence of the 3 min features as the response of the solar atmosphere to the acoustic waves which shifts high frequency wave energy to low frequencies. We also find qualitative agreement for the acoustic power between the wave simulations and the observations.

We study the response of the solar atmosphere to excitations by large amplitude acoustic waves with radiation damping now included. Monochromatic adiabatic waves, due to unbalanced heating, generate continuously rising chromospheric temperature plateaus in which the low frequency resonances quickly die out. All non-adiabatic calculations lead to stable mean chromospheric temperature distributions determined by shock dissipation and radiative cooling. For non-adiabatic monochromatic wave excitation, a critical frequency f_cr ~ 1/25 Hz is confirmed, which separates domains of different resonance behaviour. For waves of frequency f < f_cr, the resonances decay, while for waves of f > f_cr persistent resonance oscillations occur, which are perpetuated by shock merging. Excitation with acoustic frequency spectra produces distinct dynamical mean chromosphere models where the detailed temperature distributions depend on the shape of the assumed spectra. The stochasticity of the spectra and the ongoing shock merging lead to a persistent resonance behaviour of the atmosphere. The acoustic spectra show a distinct shape evolution with height such that at great height a pure 3 min band becomes increasingly dominant. With our Eulerian code we did not find appreciable mass flows at the top boundary.