TY - JOUR

T1 - Multi-dimensional shear shallow water flows

T2 - Problems and solutions

AU - Gavrilyuk, S.

AU - Ivanova, K.

AU - Favrie, N.

N1 - Publisher Copyright: © 2018 Elsevier Inc.

PY - 2018/8/1

Y1 - 2018/8/1

N2 - The mathematical model of shear shallow water flows of constant density is studied. This is a 2D hyperbolic non-conservative system of equations that is mathematically equivalent to the Reynolds-averaged model of barotropic turbulent flows. The model has three families of characteristics corresponding to the propagation of surface waves, shear waves and average flow (contact characteristics). The system is non-conservative: for six unknowns (the fluid depth, two components of the depth averaged horizontal velocity, and three independent components of the symmetric Reynolds stress tensor) one has only five conservation laws (conservation of mass, momentum, energy and mathematical ‘entropy’). A splitting procedure for solving such a system is proposed allowing us to define a weak solution. Each split subsystem contains only one family of waves (either surface or shear waves) and contact characteristics. The accuracy of such an approach is tested on 2D analytical solutions describing the flow with linear with respect to the space variables velocity, and on the solutions describing 1D roll waves. The capacity of the model to describe the full transition scenario as commonly seen in the formation of roll waves: from uniform flow to 1D roll waves, and, finally, to 2D transverse ‘fingering’ of the wave profiles, is shown.

AB - The mathematical model of shear shallow water flows of constant density is studied. This is a 2D hyperbolic non-conservative system of equations that is mathematically equivalent to the Reynolds-averaged model of barotropic turbulent flows. The model has three families of characteristics corresponding to the propagation of surface waves, shear waves and average flow (contact characteristics). The system is non-conservative: for six unknowns (the fluid depth, two components of the depth averaged horizontal velocity, and three independent components of the symmetric Reynolds stress tensor) one has only five conservation laws (conservation of mass, momentum, energy and mathematical ‘entropy’). A splitting procedure for solving such a system is proposed allowing us to define a weak solution. Each split subsystem contains only one family of waves (either surface or shear waves) and contact characteristics. The accuracy of such an approach is tested on 2D analytical solutions describing the flow with linear with respect to the space variables velocity, and on the solutions describing 1D roll waves. The capacity of the model to describe the full transition scenario as commonly seen in the formation of roll waves: from uniform flow to 1D roll waves, and, finally, to 2D transverse ‘fingering’ of the wave profiles, is shown.

KW - Godunov-type methods

KW - Non-conservative hyperbolic equations

KW - Roll waves

KW - ROLL-WAVES

KW - TO-DETONATION TRANSITION

KW - MODEL

KW - JUMP

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

U2 - 10.1016/j.jcp.2018.04.011

DO - 10.1016/j.jcp.2018.04.011

M3 - Article

AN - SCOPUS:85045245826

VL - 366

SP - 252

EP - 280

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -