Numerical superposition of Gaussian beams over propagating domain for high frequency waves and high-order invariant-preserving methods for dispersive waves