High-precision numerical methods are utilized to study the shock waves interacting with an elliptical heavy bubble.The influence of different bubble gases(SF6 and R22)and shock intensities(Ma=1.21 and Ma=2.1)is analyzed qualitatively and quantitatively.The results show that the focusing position is farther from the right bubble interface in the SF6 bubble(case 1)than in the R22 bubble(case 2)when Ma=1.21;thus,case 2 exhibits an outward jet structure,while case 1 does not.When Ma=2.1(case 3),the shock wave propagates faster,and the shock focusing is nearer to the right bubble interface.Finally,outward jet structures form on the bubble interfaces.The maximum values of density and pressure of shock focusing are different in the three cases,which imply that different gas densities and shock intensities significantly affect the shock-bubble interaction.The effective bubble volume and gases mixing degree are both smaller in case 2 than in case 1,but the trends over time are essentially the same.The increased shock intensity in case 3 leads to a smaller effective bubble volume than in case 1,but a much greater gases mixing degree.In all three cases,the compression term plays a more important role in the vorticity development than the other terms.