A numerical calculation for internal waves over topography

The internal waves propagating from the deep to shallow, and the shallow to deep, areas in the two-layer fluid systems, have been numerically simulated by solving the set of nonlinear equations, based on the variational principle in consideration of both the strong nonlinearity and strong dispersion of internal waves. The incident wave in the deep area, is the BO-type downward convex internal wave, which is the numerical solution obtained for the present fundamental equations. In the cases where the interface elevation is below, or equal to, the critical level in the shallow area, the disintegration of the internal waves occurs remarkably, leading to a long wave train. The lowest elevation of the interface, increases after its gradual decrease in the shallow area, where the interface is above the critical level, while the lowest elevation of the interface, increases through the internal-wave propagation in the shallow area, where the interface elevation is below, or equal to, the critical level, after its steep decrease around the boundary between the area over the upslope, and the shallow region.