Hodograph Method for Solving the Problem of Shallow Water Under a Solid Lid

Abstract

One of the mathematical models describing the behavior of two horizontally infinite adjoining layers of an ideal incompressible liquid under a solid cover moving at different speeds is investigated. At a large difference in the layer velocities, the Kelvin-Helmholtz instability occurs, which leads to a distortion of the interface. At the initial point in time, the interface is not necessarily flat. From a mathematical point of view, the behavior of the liquid layers is described by a system of four quasilinear equations, either hyperbolic or elliptic, in partial derivatives of the first order. Some type shallow water equations are used to construct the model. In the simple version of the model considered in this paper, in the spatially one-dimensional case, the unknowns are the boundary between the liquid layers h(x,t) and the difference in their velocities γ(x,t). The main attention is paid to the case of elliptic equations when |h|<1 and γ>1. An evolutionary Cauchy problem with arbitrary sufficiently smooth initial data is set for the system of equations. The explicit dependence of the Riemann invariants on the initial variables of the problem is indicated. To solve the Cauchy problem formulated in terms of Riemann invariants, a variant of the hodograph method based on a certain conservation law is used. This method allows us to convert a system of two quasilinear partial differential equations of the first order to a single linear partial differential equation of the second order with variable coefficients. For a linear equation, the Riemann-Green function is specified, which is used to construct a two-parameter implicit solution to the original problem. The explicit solution of the problem is constructed on the level lines (isochrons) of the implicit solution by solving a certain Cauchy problem for a system of ordinary differential equations. As a result, the original Cauchy problem in partial derivatives of the first order is transformed to the Cauchy problem for a system of ordinary differential equations, which is solved by numerical methods. Due to the bulkiness of the expression for the Riemann-Green function, some asymptotic approximation of the problem is considered, and the results of calculations, and their analysis are presented.