A fully well-balanced hydrodynamic reconstruction

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Berthon, Christophe, Michel-Dansac, Victor
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866910776000774144
author Berthon, Christophe
Michel-Dansac, Victor
author_facet Berthon, Christophe
Michel-Dansac, Victor
contents The present work focuses on the numerical approximation of the weak solutions of the shallow water model over a non-flat topography. In particular, we pay close attention to steady solutions with nonzero velocity. The goal of this work is to derive a scheme that exactly preserves these stationary solutions, as well as the commonly preserved lake at rest steady solution. These moving steady states are solution to a nonlinear equation. We emphasize that the method proposed here never requires solving this nonlinear equation; instead, a suitable linearization is derived. To address this issue, we propose an extension of the well-known hydrostatic reconstruction. By appropriately defining the reconstructed states at the interfaces, any numerical flux function, combined with a relevant source term discretization, produces a well-balanced scheme that preserves both moving and non-moving steady solutions. This eliminates the need to construct specific numerical fluxes. Additionally, we prove that the resulting scheme is consistent with the homogeneous system on flat topographies, and that it reduces to the hydrostatic reconstruction when the velocity vanishes. To increase the accuracy of the simulations, we propose a well-balanced high-order procedure, which still does not require solving any nonlinear equation. Several numerical experiments demonstrate the effectiveness of the numerical scheme.
format Preprint
id arxiv_https___arxiv_org_abs_2501_03975
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle A fully well-balanced hydrodynamic reconstruction
Berthon, Christophe
Michel-Dansac, Victor
Numerical Analysis
65M08
The present work focuses on the numerical approximation of the weak solutions of the shallow water model over a non-flat topography. In particular, we pay close attention to steady solutions with nonzero velocity. The goal of this work is to derive a scheme that exactly preserves these stationary solutions, as well as the commonly preserved lake at rest steady solution. These moving steady states are solution to a nonlinear equation. We emphasize that the method proposed here never requires solving this nonlinear equation; instead, a suitable linearization is derived. To address this issue, we propose an extension of the well-known hydrostatic reconstruction. By appropriately defining the reconstructed states at the interfaces, any numerical flux function, combined with a relevant source term discretization, produces a well-balanced scheme that preserves both moving and non-moving steady solutions. This eliminates the need to construct specific numerical fluxes. Additionally, we prove that the resulting scheme is consistent with the homogeneous system on flat topographies, and that it reduces to the hydrostatic reconstruction when the velocity vanishes. To increase the accuracy of the simulations, we propose a well-balanced high-order procedure, which still does not require solving any nonlinear equation. Several numerical experiments demonstrate the effectiveness of the numerical scheme.
title A fully well-balanced hydrodynamic reconstruction
topic Numerical Analysis
65M08
url https://arxiv.org/abs/2501.03975