The initial data of a general Riemann problem is constant along radial directions from an origin and it is piecewise constant as a function of angle. In this thesis, the special case with initial data of piecewise constant solutions joined by four forward rarefaction waves is considered. A transonic shock is discovered through refined numerical simulations and generalized characteristic analysis, and the mathematical mechanism of shock formation for four rarefactions cases is clarified. To our best knowledge, shock formation in a rarefaction wave-only 2D Riemann problem, had not been...