Large-amplitude internal solitary waves in a two-layer fluid system with a free surface are investigated in this paper. Laboratory experiments on strongly nonlinear internal solitary waves with a free surface for a deep configuration are conducted. After comparing the experimental data with the results of the Miyata–Choi–Camassa model that includes the free-surface effects (MCC-FS), we find that the MCC-FS model does not calculate accurately the internal solitary waves with a free surface. Thus, we develop a strongly nonlinear model for a deep configuration, namely the two-layer high-level Green–Naghdi (HLGN-FS) model that includes the free-surface effects. Numerical results of the HLGN-FS model, including the wave profile, velocity field and wave speed, are presented for three cases. The first case is a shallow configuration with . It is shown that the MCC-FS model can provide accurate results for the shallow configurations. Meanwhile, the HLGN-FS model is shown to be accurate for describing the internal solitary waves for both shallow and deep configurations.

The strongly nonlinear Miyata–Choi–Camassa model under the rigid lid approximation (MCC-RL model) can describe accurately the dynamics of large-amplitude internal waves in a two-layer fluid system for shallow configurations. In this paper, we apply the MCC-RL model to study the internal waves generated by a moving body on the bottom. For the case of the moving body speed $U=1.1c_{0}$ , where ${c_0}$ is the linear long-wave speed, the accuracy of the MCC-RL results is assessed by comparing with Euler's solutions, and very good agreement is observed. It is found that when the moving body speed increases from $U=0.8c_{0}$ to $U=1.241c_{0}$ , the amplitudes of the generated internal solitary waves in front of the moving body become larger. However, a critical moving body speed is found between $U=1.241c_{0}$ and $U=1.242c_{0}$ . After exceeding this critical speed, only one internal wave right above the body is generated. When the moving body speed increases from $U=1.242c_{0}$ to $U=1.5c_{0}$ , the amplitudes of the internal waves become smaller.