First-principles quantum simulations of many-mode open interacting Bose gases using stochastic gauge methods