Internal waves arise in a wide array of oceanographic problems of both theoretical and engineering interest. In this contribution we present a new model, valid in the weakly nonlinear regime, for the propagation of disturbances along the interface between two ideal fluid layers of infinite extent and different densities. Additionally, we present a novel high-order/spectral algorithm for its accurate and stable simulation. Numerical validation results and simulations of wave-packet evolution are provided.
