We study the scaling laws and structure functions of stratified shear flows by performing high-resolution numerical simulations of inviscid compressible turbulence induced by Kelvin-Helmholtz instability. An implicit large eddy simulation approach is adapted to solve our conservation laws for both two-dimensional (with a spatial resolution of 16,384<sup>2</sup>) and three-dimensional (with a spatial resolution of 512<sup>3</sup>) configurations utilizing different compressibility characteristics such as shocks. For three-dimensional turbulence, we find that both kinetic energy and density-weighted energy spectra follow the classical Kolmogorov <i>k</i><sup>−5/3</sup> inertial scaling. This phenomenon is observed due to the fact that the power density spectrum of three-dimensional turbulence yields the same <i>k</i><sup>−5/3</sup> scaling. However, we demonstrate that there is a significant difference between these two spectra in two-dimensional turbulence since the power density spectrum flattens to <i>k</i><sup>−1/4</sup>. This flattening may be assumed to be a reason for the <i>k</i><sup>−7/3</sup> scaling observed in the two-dimensional density-weight kinetic every spectra for high compressibility as compared to the <i>k</i><sup>−3</sup> scaling traditionally assumed with incompressible flows. Further inquiries are made to validate the statistical behavior of the various configurations studied through the use of second and third order velocity structure functions where it is noticed that scaling behavior differs between the two- and three-dimensional cases wherein only the latter is seen to follow trends from K41 theory.