English   español  
Please use this identifier to cite or link to this item: http://hdl.handle.net/10261/209589
logo share SHARE logo core CORE   Add this article to your Mendeley library MendeleyBASE

Visualizar otros formatos: MARC | Dublin Core | RDF | ORE | MODS | METS | DIDL | DATACITE
Exportar a otros formatos:


A 1D numerical model for the simulation of unsteady and highly erosive flows in rivers

AuthorsMartínez-Aranda, S.; Murillo, J.; García-Navarro, Pilar
KeywordsFinite volume method
Shallow water equations
Exner equation
Coupled method
Augmented Roe’s scheme
Variable geometry
Issue Date15-Mar-2019
CitationComputers and Fluids 181: 8-34 (2019)
AbstractThis work is focused on a numerical finite volume scheme for the coupled shallow water-Exner system in 1D applications with arbitrary geometry. The mathematical expressions modeling the hydrodynamic and morphodynamic components of the physical phenomenon are treated to deal with cross-section shape variations and empirical solid discharge estimations. The resulting coupled equations can be rewritten as a non-conservative hyperbolic system with three moving waves and one stationary wave to account for the source terms discretization. Moreover, the wave celerities for the coupled morpho-hydrodyamical system depend on the erosion-deposition mechanism selected to update the channel cross-section profile. This influence is incorporated into the system solution by means of a new parameter related to the channel bottom variation celerity. Special interest is put to show that, even for the simplest solid transport models as the Grass law, to find a linearized Jacobian matrix of the system can be a challenge in presence of arbitrary shape channels. In this paper a numerical finite volume scheme is proposed, based on an augmented Roe solver, first order accurate in time and space, dealing with solid transport flux variations caused by the channel geometry changes. Channel cross-section variations lead to the appearance of a new solid flux source term which should be discretized properly. The stability region is controlled by wave celerities together with a proper reconstruction of the approximate local Riemann problem solution, enforcing positive values for the intermediate states of the conserved variables. Comparison of the numerical results for several analytical and experimental cases demonstrates the effectiveness, exact well-balancedness and accuracy of the scheme.
Publisher version (URL)https://doi.org/10.1016/j.compfluid.2019.01.011
Appears in Collections:(LIFTEC) Artículos
Files in This Item:
File Description SizeFormat 
accesoRestringido.pdf59,24 kBAdobe PDFThumbnail
Show full item record
Review this work

Related articles:

WARNING: Items in Digital.CSIC are protected by copyright, with all rights reserved, unless otherwise indicated.