The full potential of the hybrid large-particle method using the known and new Riemann problems in two-dimensional domains has been studied. The method includes a space-time second-order approximation for smooth solutions. Using the artificial viscosity nonlinear correction and the convective fluxes hybridization maintained monotonicity of solutions. Centrally symmetrical problems with a complex shock-wave structure and with the development of instability on the contact boundary were studied in details. The test calculations demonstrated high resolution, low dissipation, and stability of the method.