Numerical study of GaN-on-Si HEMT breakdown instability accounting for substrate and packaging interactions