Higher order shadow waves and delta shock blow up in the Chaplygin gas

The introductory part of this paper contains an overview of known results about elementary and delta shock solutions to Riemann problem for well known Chaplygin gas model (nowadays used in cosmological theories for dark energy) in terms of entropic shadow waves. Shadow waves are introduced in 17] and they are represented by shocks depending on a small parameter ε with unbounded amplitudes having a distributional limit involving the Dirac delta function. In a search for admissible solutions to all possible cases of mutual interactions of waves arising from double Riemann initial data we found same cases that cannot be resolved with already known types of elementary or shadow wave solutions. These cases are resolved by introducing a sequence of higher order shadow waves depending on integer powers of ε. It is shown that such waves have a distributional limit but only until some finite time T.