A fully well-balanced hydrodynamic reconstruction
-
Christophe Berthon
Abstract
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.
Funding statement: C. Berthon acknowledges the support of ANR MUFFIN ANR-19-CE46-0004, funded by Agence Nationale de la Recherche. The SHARK-FV conference has greatly contributed to this work.
References
[1] M. Abramowitz and I. A. Stegun (eds.), Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover Publications, Inc., New York, 1992.Suche in Google Scholar
[2] A. I. Aleksyuk, M. A. Malakhov, and V. V. Belikov, The exact Riemann solver for the shallow water equations with a discontinuous bottom, J. Comput. Phys. 450 (2022), 110801.10.1016/j.jcp.2021.110801Suche in Google Scholar
[3] E. Audusse, F. Bouchut, M.-O. Bristeau, R. Klein, and B. Perthame, A fast and stable well-balanced scheme with hydrostatic reconstruction for shallow water flows, SIAM J. Sci. Comput. 25 (2004), No. 6, 2050–2065.10.1137/S1064827503431090Suche in Google Scholar
[4] A. Bermúdez and M. E. Vázquez, Upwind methods for hyperbolic conservation laws with source terms, Comput. & Fluids 23 (1994), No. 8, 1049–1071.10.1016/0045-7930(94)90004-3Suche in Google Scholar
[5] C. Berthon, Stability of the MUSCL schemes for the Euler equations, Commun. Math. Sci. 3 (2005), No. 2, 133–157.10.4310/CMS.2005.v3.n2.a3Suche in Google Scholar
[6] C. Berthon, S. Bulteau, F. Foucher, M. M’Baye, and V. Michel-Dansac, A very easy high-order well-balanced reconstruction for hyperbolic systems with source terms, SIAM J. Sci. Comput. 44 (2022), No. 4, A2506–A2535.10.1137/21M1429230Suche in Google Scholar
[7] C. Berthon and C. Chalons, A fully well-balanced, positive and entropy-satisfying Godunov-type method for the shallow-water equations, Math. Comp. 85 (2016), No. 299, 1281–1307.10.1090/mcom3045Suche in Google Scholar
[8] C. Berthon, A. Duran, F. Foucher, K. Saleh, and J. D. D. Zabsonré, Improvement of the hydrostatic reconstruction scheme to get fully discrete entropy inequalities, J. Sci. Comput. 80 (2019), No. 2, 924–956.10.1007/s10915-019-00961-ySuche in Google Scholar
[9] C. Berthon and F. Foucher, Efficient well-balanced hydrostatic upwind schemes for shallow-water equations, J. Comput. Phys. 231 (2012), No. 15, 4993–5015.10.1016/j.jcp.2012.02.031Suche in Google Scholar
[10] C. Berthon, M. M’Baye, M. H. Le, and D. Seck, A well-defined moving steady states capturing Godunov-type scheme for shallow-water model, Int. J. Finite Vol. 15 (2021).Suche in Google Scholar
[11] F. Bouchut and T. M. de Luna, A subsonic-well-balanced reconstruction scheme for shallow water flows, SIAM J. Numer. Anal. 48 (2010), No. 5, 1733–1758.10.1137/090758416Suche in Google Scholar
[12] J. Britton and Y. Xing, High order still-water and moving-water equilibria preserving discontinuous Galerkin methods for the Ripa model, J. Sci. Comput. 82 (2020), No. 2, 30.10.1007/s10915-020-01134-ySuche in Google Scholar
[13] M. Castro, J. M. Gallardo, J. A. López-García, and C. Parés, Well-balanced high order extensions of Godunov’s method for semilinear balance laws, SIAM J. Numer. Anal. 46 (2008), No. 2, 1012–1039.10.1137/060674879Suche in Google Scholar
[14] M. J. Castro, A. P. Milanés, and C. Parés, Well-balanced numerical schemes based on a generalized hydrostatic reconstruction technique, Math. Models Methods Appl. Sci. 17 (2007), No. 12, 2055–2113.10.1142/S021820250700256XSuche in Google Scholar
[15] M. J. Castro and C. Parés, Well-balanced high-order finite volume methods for systems of balance laws, J. Sci. Comput. 82 (2020), No. 2, 48.10.1007/s10915-020-01149-5Suche in Google Scholar
[16] O. Castro-Orgaz and H. Chanson, Minimum specific energy and transcritical flow in unsteady open-channel flow, J. Irrig. Drain. Engrg. – ASCE 142 (2016), No. 1, 04015030.10.1061/(ASCE)IR.1943-4774.0000926Suche in Google Scholar
[17] G. Chen and S. Noelle, A new hydrostatic reconstruction scheme based on subcell reconstructions, SIAM J. Numer. Anal. 55 (2017), No. 2, 758–784.10.1137/15M1053074Suche in Google Scholar
[18] T. M. de Luna, M. J. Castro Díaz, and C. Parés, Reliability of first order numerical schemes for solving shallow water system over abrupt topography, Appl Math Comput 219 (2013), No. 17, 9012–9032.10.1016/j.amc.2013.03.033Suche in Google Scholar
[19] O. Delestre and P.-Y. Lagrée, A ‘well-balanced’ finite volume scheme for blood flow simulation, Int. J. Numer. Methods Fluids 72 (2012), No. 2, 177–205.10.1002/fld.3736Suche in Google Scholar
[20] S. Diot, S. Clain, and R. Loubère, Improved detection criteria for the multi-dimensional optimal order detection (MOOD) on unstructured meshes with very high-order polynomials, Comput. & Fluids 64 (2012), 43–63.10.1016/j.compfluid.2012.05.004Suche in Google Scholar
[21] S. Diot, R. Loubère, and S. Clain, The multidimensional optimal order detection method in the three-dimensional case: very high-order finite volume method for hyperbolic systems, Int. J. Numer. Methods Fluids 73 (2013), No. 4, 362–392.10.1002/fld.3804Suche in Google Scholar
[22] A. Duran, J.-P. Vila, and R. Baraille, Energy-stable staggered schemes for the shallow water equations, J. Comput. Phys. 401 (2020), 109051.10.1016/j.jcp.2019.109051Suche in Google Scholar
[23] I. Gómez-Bueno, M. J. Castro, and C. Parés, High-order well-balanced methods for systems of balance laws: a control-based approach, Appl. Math. Comput. 394 (2021), 125820.10.1016/j.amc.2020.125820Suche in Google Scholar
[24] I. Gómez-Bueno, M. J. Castro Díaz, C. Parés, and G. Russo, Collocation methods for high-order well-balanced methods for systems of balance laws, Mathematics 9 (2021), No. 15, 1799.10.3390/math9151799Suche in Google Scholar
[25] L. Gosse, A well-balanced flux-vector splitting scheme designed for hyperbolic systems of conservation laws with source terms, Comput. Math. Appl. 39 (2000), No. 9-10, 135–159.10.1016/S0898-1221(00)00093-6Suche in Google Scholar
[26] S. Gottlieb and C.-W. Shu, Total variation diminishing Runge–Kutta schemes, Math. Comp. 67 (1998), No. 221, 73–85.10.1090/S0025-5718-98-00913-2Suche in Google Scholar
[27] S. Gottlieb, C.-W. Shu, and E. Tadmor, Strong stability-preserving high-order time discretization methods, SIAM Rev. 43 (2001), No. 1, 89–112.10.1137/S003614450036757XSuche in Google Scholar
[28] N. Goutal and F. Maurel, Proc. of the 2nd Workshop on Dam-Break Wave Simulation, Groupe Hydraulique Fluviale, Département Laboratoire National d’Hydraulique, Electricité de France, Report, 1997.Suche in Google Scholar
[29] J. M. Greenberg and A.-Y. LeRoux, A well-balanced scheme for the numerical processing of source terms in hyperbolic equations, SIAM J. Numer. Anal. 33 (1996), No. 1, 1–16.10.1137/0733001Suche in Google Scholar
[30] A. Harten, P. D. Lax, and B. van Leer, On upstream differencing and Godunov-type schemes for hyperbolic conservation laws, SIAM Rev. 25 (1983), No. 1, 35–61.10.1137/1025002Suche in Google Scholar
[31] S. Jin, A steady-state capturing method for hyperbolic systems with geometrical source terms, M2AN Math. Model. Numer. Anal. 35 (2001), No. 4, 631–645.10.1051/m2an:2001130Suche in Google Scholar
[32] P. G. LeFloch and M. D. Thanh, The Riemann problem for the shallow water equations with discontinuous topography, Commun. Math. Sci. 5 (2007), No. 4, 865–885.10.4310/CMS.2007.v5.n4.a7Suche in Google Scholar
[33] P. G. LeFloch and M. D. Thanh, A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime, J. Comput. Phys. 230 (2011), No. 20, 7631–7660.10.1016/j.jcp.2011.06.017Suche in Google Scholar
[34] R. J. LeVeque, Finite volume methods for hyperbolic problems, Cambridge Texts in Applied Mathematics, Cambridge University Press, Cambridge, 2002.10.1017/CBO9780511791253Suche in Google Scholar
[35] G. Li and Y. Xing, Well-balanced discontinuous Galerkin methods with hydrostatic reconstruction for the Euler equations with gravitation, J. Comput. Phys. 352 (2018), 445–462.10.1016/j.jcp.2017.09.063Suche in Google Scholar
[36] V. Michel-Dansac, C. Berthon, S. Clain, and F. Foucher, A well-balanced scheme for the shallow-water equations with topography, Comput. Math. Appl. 72 (2016), No. 3, 568–593.10.1016/j.camwa.2016.05.015Suche in Google Scholar
[37] V. Michel-Dansac, C. Berthon, S. Clain, and F. Foucher, A well-balanced scheme for the shallow-water equations with topography or Manning friction, J. Comput. Phys. 335 (2017), 115–154.10.1016/j.jcp.2017.01.009Suche in Google Scholar
[38] V. Michel-Dansac, C. Berthon, S. Clain, and F. Foucher, A two-dimensional high-order well-balanced scheme for the shallow water equations with topography and Manning friction, Comput. & Fluids 230 (2021), 105152.10.1016/j.compfluid.2021.105152Suche in Google Scholar
[39] R. Natalini, M. Ribot, and M. Twarogowska, A well-balanced numerical scheme for a one dimensional quasilinear hyperbolic model of chemotaxis, Commun. Math. Sci. 12 (2014), No. 1, 13–39.10.4310/CMS.2014.v12.n1.a2Suche in Google Scholar
[40] S. Noelle, Y. Xing, and C.-W. Shu, High-order well-balanced finite volume WENO schemes for shallow water equation with moving water, J. Comput. Phys. 226 (2007), No. 1, 29–58.10.1016/j.jcp.2007.03.031Suche in Google Scholar
[41] B. Schmidtmann, B. Seibold, and M. Torrilhon, Relations between WENO3 and third-order limiting in finite volume methods, J. Sci. Comput. 68 (2015), No. 2, 624–652.10.1007/s10915-015-0151-zSuche in Google Scholar
[42] E. F. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics. A Practical Introduction, 3rd ed., Springer-Verlag, Berlin, 2009.10.1007/b79761Suche in Google Scholar
[43] B. van Leer, Towards the ultimate conservative difference scheme, V. A second order sequel to Godunov’s method, J. Comput. Phys. 32 (1979), 101–136.10.1016/0021-9991(79)90145-1Suche in Google Scholar
[44] J.-P. Vila, Simplified Godunov schemes for 2 × 2 systems of conservation laws, SIAM J. Numer. Anal. 23 (1986), No. 6, 1173–1192.10.1137/0723079Suche in Google Scholar
[45] Y. Xing, Exactly well-balanced discontinuous Galerkin methods for the shallow water equations with moving water equilibrium, J. Comput. Phys. 257 (2014), No. part A, 536–553.10.1016/j.jcp.2013.10.010Suche in Google Scholar
[46] Y. Xing and C.-W. Shu, High-order finite volume WENO schemes for the shallow water equations with dry states, Adv. Water Resour. 34 (2011), No. 8, 1026–1038.10.1016/j.advwatres.2011.05.008Suche in Google Scholar
[47] Y. Xing, C.-W. Shu, and S. Noelle, On the advantage of well-balanced schemes for moving-water equilibria of the shallow water equations, J. Sci. Comput. 48 (2011), No. 1-3, 339–349.10.1007/s10915-010-9377-ySuche in Google Scholar
[48] Y. Xing, X. Zhang, and C.-W. Shu, Positivity-preserving high order well-balanced discontinuous Galerkin methods for the shallow water equations, Adv. Water Resour. 33 (2010), No. 12, 1476–1493.10.1016/j.advwatres.2010.08.005Suche in Google Scholar
A Taylor expansions of 𝓗
The goal here is to provide a Taylor expansion of the function 𝓗 given by (4.5), in the case where ΔZ > 0, Δh > 0, and 1 – Fr2 > 0, when ΔZ approaches zero. The computations are performed below, where we have temporarily set 𝔽 = 1 – Fr2 in order to save some space:
B Taylor expansions of 𝓒
In this section, we give the Taylor expansions of the function 𝓒(hL, hR, qL, ΔZ), given by (4.7), when either hL or hR go to 0. The goal is to prove H3, i.e., prove that 𝓒 is continuous when either hL tends to 0 and hR ≠ 0, and when hL tends to 0 and hR = hL. Recall that
Since qL = hLuL, with u the velocity, we also note that
First, we consider the case where hL goes to 0 and hR ≠ 0. According to assumptions (1.3), in this case, uL also goes to zero. To model this phenomenon, we assume that uL = u(hL), where the function u is such that u(0) = 0. In this case, we get, again using symbolic computation software,
In the above expression, for the sake of clarity, the ± symbols correspond to sgn(ΔZ). In any case, since u(0) = 0, this Taylor expansion shows that
which is what we had set out to prove.
Second, we have to prove that
However, recall from (4.6) that 𝓗 = (hR – hL) 𝓑, with 𝓑 a bounded function. The above result is established by arguing the boundedness (1.3) of the Froude number.
Therefore, property H3 is satisfied by 𝓗, when 𝓗 is given by (4.5).
© 2024 Walter de Gruyter GmbH, Berlin/Boston
Artikel in diesem Heft
- Frontmatter
- Efficient numerical solution of the Fokker–Planck equation using physics-conforming finite element methods
- Stability and convergence of relaxed scalar auxiliary variable schemes for Cahn–Hilliard systems with bounded mass source
- Local parameter selection in the C0 interior penalty method for the biharmonic equation
- A fully well-balanced hydrodynamic reconstruction
Artikel in diesem Heft
- Frontmatter
- Efficient numerical solution of the Fokker–Planck equation using physics-conforming finite element methods
- Stability and convergence of relaxed scalar auxiliary variable schemes for Cahn–Hilliard systems with bounded mass source
- Local parameter selection in the C0 interior penalty method for the biharmonic equation
- A fully well-balanced hydrodynamic reconstruction
