Two-layer interfacial flows beyond the Boussinesq approximation: a Hamiltonian approach

The theory of integrable systems of Hamiltonian PDEs and their near-integrable deformations is used to study evolution equations resulting from vertical-averages of the Euler system for two-layer stratified flows in an infinite 2D channel. The Hamiltonian structure of the averaged equations is obtained directly from that of the Euler equations through the process of Hamiltonian reduction. Long-wave asymptotics together with the Boussinesq approximation of neglecting the fluids' inertia is then applied to reduce the leading order vertically averaged equations to the shallow-water Airy system, and thence, in a non-trivial way, to the dispersionless non-linear Schr\"odinger equation. The full non-Boussinesq system for the dispersionless limit can then be viewed as a deformation of this well known equation. In a perturbative study of this deformation, it is shown that at first order the deformed system possesses an infinite sequence of constants of the motion, thus casting this system within the framework of completely integrable equations. The Riemann invariants of the deformed model are then constructed, and some local solutions found by hodograph-like formulae for completely integrable systems are obtained.
Comment: 26 pages; 8 figures