143 research outputs found
A robust well-balanced scheme for multi-layer shallow water equations
International audienceThe numerical resolution of the multi-layer shallow water system encounters two additional difficulties with respect to the one-layer system. The first is that the system involves nonconservative terms, and the second is that it is not always hyperbolic. A splitting scheme has been proposed by Bouchut and Morales, that enables to ensure a discrete entropy inequality and the well-balanced property, without any theoretical difficulty related to the loss of hyperbolicity. However, this scheme has been shown to often give wrong solutions. We introduce here a variant of the splitting scheme, that has the overall property of being conservative in the total momentum. It is based on a source-centered hydrostatic scheme for the one-layer shallow water system, a variant of the hydrostatic scheme. The final method enables to treat an arbitrary number of layers, with arbitrary densities and arbitrary topography. It has no restriction concerning complex eigenvalues, it is well-balanced and it is able to treat vacuum, it satisfies a semi-discrete entropy inequality. The scheme is fast to execute, as is the one-layer hydrostatic method
On the rigid-lid approximation for two shallow layers of immiscible fluids with small density contrast
The rigid-lid approximation is a commonly used simplification in the study of
density-stratified fluids in oceanography. Roughly speaking, one assumes that
the displacements of the surface are negligible compared with interface
displacements. In this paper, we offer a rigorous justification of this
approximation in the case of two shallow layers of immiscible fluids with
constant and quasi-equal mass density. More precisely, we control the
difference between the solutions of the Cauchy problem predicted by the
shallow-water (Saint-Venant) system in the rigid-lid and free-surface
configuration. We show that in the limit of small density contrast, the flow
may be accurately described as the superposition of a baroclinic (or slow)
mode, which is well predicted by the rigid-lid approximation; and a barotropic
(or fast) mode, whose initial smallness persists for large time. We also
describe explicitly the first-order behavior of the deformation of the surface,
and discuss the case of non-small initial barotropic mode.Comment: Compared to version 2, typos have been corrected and additional
remarks/discussion added. To appear in Journal of Nonlinear Scienc
A consistent reduction of the two-layer shallow-water equations to an accurate one-layer spreading model
The gravity-driven spreading of one fluid in contact with another fluid is of
key importance to a range of topics. To describe these phenomena, the two-layer
shallow-water equations is commonly employed. When one layer is significantly
deeper than the other, it is common to approximate the system with the much
simpler one-layer shallow water equations. So far, it has been assumed that
this approximation is invalid near shocks, and one has applied additional front
conditions for the shock speed. In this paper, we prove mathematically that an
effective one-layer model can be derived from the two-layer equations that
correctly captures the behaviour of shocks and contact discontinuities without
any additional closure relations. The proof yields a novel formulation of an
effective one-layer shallow water model. The result shows that simplification
to an effective one-layer model is well justified mathematically and can be
made without additional knowledge of the shock behaviour. The shock speed in
the proposed model is consistent with empirical models and identical to the
front conditions that have been found theoretically by e.g. von K\'arm\'an and
by Benjamin. This suggests that the breakdown of the shallow-water equations in
the vicinity of shocks is less severe than previously thought. We further
investigate the applicability of the shallow water framework to shocks by
studying shocks in one-dimensional lock-exchange/lock-release. We derive
expressions for the Froude number that are in good agreement with the widely
employed expression by Benjamin. We then solve the equations numerically to
illustrate how quickly the proposed model converges to solutions of the full
two-layer shallow-water equations. We also compare numerical results using our
model with results from dam break experiments. Predictions from the one-layer
model are found to be in good agreement with experiments.Comment: 23 pages, 17 figure
- âŠ