In this article, we have developed a new time-marching method to simulate sound and soliton propagation in homogeneous and heterogeneous media. The proposed time integration scheme can numerically preserve physical dispersion over a wide wavenumber range and conserve energy while solving wave propagation problems. Here, the Fourier stability analysis

has been used to assess the numerical properties of the developed method. The proposed numerical method's dispersion and dissipation properties have also been compared with the classical fourth-order Runge-Kutta (RK4) method. Stability property contours for the newly proposed method display that the maximum allowable time step is at least five times higher than the RK4 method. The Fourier stability analysis also explains the dispersion error associated with the used spatio-temporal discretization schemes. It is observed that the dispersion error is significantly small for the proposed time integration schemes compared with the RK4 method. The proposed

methods simulate sound propagation problems with fewer computational resources that otherwise demand high computational costs. The efficacy of the proposed time integration methods has been demonstrated by solving

benchmark sound wave propagation problems. Moreover, to test the developed method's efficiency and robustness, we have performed simulations of the sound wave propagation in a layered media, corner-edge model, and damped sine-Gordon equation.