A high-order hybrid flow solver for the simulation of nonlinear internal waves in long complex domains