TY - JOUR

T1 - Nonlinear stability of two-layer shallow water flows with a free surface

AU - De Melo Viríssimo, Francisco

AU - Milewski, Paul A.

N1 - Publisher Copyright:
© 2020 The Author(s) Published by the Royal Society. All rights reserved.

PY - 2020/4/1

Y1 - 2020/4/1

N2 - The problem of two layers of immiscible fluid, bordered above by an unbounded layer of passive fluid and below by a flat bed, is formulated and discussed. The resulting equations are given by a first-order, four-dimensional system of PDEs of mixed-type. The relevant physical parameters in the problem are presented and used to write the equations in a non-dimensional form. The conservation laws for the problem, which are known to be only six, are explicitly written and discussed in both non- Boussinesq and Boussinesq cases. Both dynamics and nonlinear stability of the Cauchy problem are discussed, with focus on the case where the upper unbounded passive layer has zero density, also called the free surface case. We prove that the stability of a solution depends only on two 'baroclinic' parameters (the shear and the difference of layer thickness, the former being the most important one) and give a precise criterion for the system to be well-posed. It is also numerically shown that the system is nonlinearly unstable, as hyperbolic initial data evolves into the elliptic region before the formation of shocks. We also discuss the use of simple waves as a tool to bound solutions and preventing a hyperbolic initial data to become elliptic and use this idea to give a mathematical proof for the nonlinear instability.

AB - The problem of two layers of immiscible fluid, bordered above by an unbounded layer of passive fluid and below by a flat bed, is formulated and discussed. The resulting equations are given by a first-order, four-dimensional system of PDEs of mixed-type. The relevant physical parameters in the problem are presented and used to write the equations in a non-dimensional form. The conservation laws for the problem, which are known to be only six, are explicitly written and discussed in both non- Boussinesq and Boussinesq cases. Both dynamics and nonlinear stability of the Cauchy problem are discussed, with focus on the case where the upper unbounded passive layer has zero density, also called the free surface case. We prove that the stability of a solution depends only on two 'baroclinic' parameters (the shear and the difference of layer thickness, the former being the most important one) and give a precise criterion for the system to be well-posed. It is also numerically shown that the system is nonlinearly unstable, as hyperbolic initial data evolves into the elliptic region before the formation of shocks. We also discuss the use of simple waves as a tool to bound solutions and preventing a hyperbolic initial data to become elliptic and use this idea to give a mathematical proof for the nonlinear instability.

UR - http://www.scopus.com/inward/record.url?scp=85084927265&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85084927265&partnerID=8YFLogxK

U2 - 10.1098/rspa.2019.0594

DO - 10.1098/rspa.2019.0594

M3 - Article

AN - SCOPUS:85084927265

SN - 1364-5021

VL - 476

JO - Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences

JF - Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences

IS - 2236

M1 - 20190594

ER -