Numerical simulations of the flowfield structures and properties of underexpanded supersonic jet and planar shear layer are performed by solving the time-dependent, compressible Euler equations. The numerical code uses the high-order weighted essentially