*Panel methods are inviscid solutions. You will not capture viscous effects except via user "modeling" by changing the geometry.
*Solutions are invalid as soon as the flow changes locally from subsonic to supersonic (i.e. the critical Mach number has been exceeded) or vice versa.