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ödinger 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.
IntroductionAspects of the theory of two-layer stratified flows in an infinite 2D channel have been the subject of intense recent studies. Layer models are widely used in a variety of geophysical applications (going back to early references such as that by Long [20] in the framework of meteorology), and are of conceptual value for illustrating many fundamental properties of stratified fluid dynamics. A typical configuration is depicted in Figure 1, with an interface between the two fluid representing the sharp pycnocline between superficial ("fresh") water, labeled by the index 1, and deep ("salty") water, labeled by the 2-index. (Other relevant notation used throughout the paper is defined by the figure). Long internal waves in such systems were studied in, e.g., [10,11], by deriving the two-layer models (including dispersive terms) by the layer-averaging method (see e.g., [30]). Their dispersionless counterparts were more recently reconsidered in papers by Milewski, Tabak and collaborators [25,12]. These papers were mainly interested in studying the Kelvin-Helmholtz (KH) instability (see also [4]) viewed as hyperbolic vs. elliptic transition for the resulting quasi-linear equations of motion, and its relation to the well-posedness of the initial value problem for these equations. In particular, for the so-called Boussinesq approximation, which in this context consists of disregarding density differences in the inertial terms while retaining them in the buoyancy terms, the following conditions for shear-flow stability are all equivalent:i) The "standard" stability criterion expressed by the Richardson number for the linearized two-layer equations ([23, 3]) around a velocity jump.ii) The hyperbolicity (nonlinear) criterion for the reduced quasi-linear system of PDEs in the variables ξ, w (the difference between ...