A characteristics approach to shock formation in 2D Euler with azimuthal symmetry and entropy

We provide a detailed analysis of the shock formation process for the non-isentropic 2d Euler equations in azimuthal symmetry. We prove that from an open set of smooth and generic initial data, solutions of Euler form a first singularity or gradient blow-up or shock. This first singularity is termed a H\"{o}lder $C^{\frac{1}{3}}$ pre-shock, and our analysis provides the first detailed description of this cusp solution. The novelty of this work relative to [Buckmaster-Drivas-Shkoller-Vicol, 2022] is that we herein consider a much larger class of initial data, allow for a non-constant initial entropy, allow for a non-trivial sub-dominant Riemann variable, and introduce a host of new identities to avoid apparent derivative loss due to entropy gradients. The method of proof is also new and robust, exploring the transversality of the three different characteristic families to transform space derivatives into time derivatives. Our main result provides a fractional series expansion of the Euler solution about the pre-shock, whose coefficients are computed from the initial data.