Skip to main content Accessibility help
Hostname: page-component-55b6f6c457-4lvx9 Total loading time: 0.224 Render date: 2021-09-23T08:46:14.830Z Has data issue: true Feature Flags: { "shouldUseShareProductTool": true, "shouldUseHypothesis": true, "isUnsiloEnabled": true, "metricsAbstractViews": false, "figures": true, "newCiteModal": false, "newCitedByModal": true, "newEcommerce": true, "newUsageEvents": true }

IMEX Large Time Step Finite Volume Methods for Low Froude Number Shallow Water Flows

Published online by Cambridge University Press:  03 June 2015

Georgij Bispen*
Institute of Mathematics, University of Mainz, Germany
K. R. Arun*
School of Mathematics, Indian Institute of Science Education and Research Thiruvananthapuram, India
Mária Lukáčová-Medvid’ová*
Institute of Mathematics, University of Mainz, Germany
Sebastian Noelle*
IGPM, RWTH Aachen, Germany
Get access


We present new large time step methods for the shallow water flows in the low Froude number limit. In order to take into account multiscale phenomena that typically appear in geophysical flows nonlinear fluxes are split into a linear part governing the gravitational waves and the nonlinear advection. We propose to approximate fast linear waves implicitly in time and in space by means of a genuinely multidimensional evolution operator. On the other hand, we approximate nonlinear advection part explicitly in time and in space by means of the method of characteristics or some standard numerical flux function. Time integration is realized by the implicit-explicit (IMEX) method. We apply the IMEX Euler scheme, two step Runge Kutta Cranck Nicolson scheme, as well as the semi-implicit BDF scheme and prove their asymptotic preserving property in the low Froude number limit. Numerical experiments demonstrate stability, accuracy and robustness of these new large time step finite volume schemes with respect to small Froude number.

Research Article
Copyright © Global Science Press Limited 2014

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)


[1]Arun, K.R., Kraft, M., Lukáčová-Medvid’ová, M., and Prasad, Ph.Finite volume evolution Galerkin method for hyperbolic conservation laws with spatially varying flux functions. J. Comput. Phys., 228: 565590, 2009.CrossRefGoogle Scholar
[2]Arun, K.R. and Noelle, S.An asymptotic preserving scheme for low Froude number shallow flows. Accepted to the Proceedings of the 14th International Conference on Hyperbolic Problems. Theory, Numerics and Applications, American Institute of Mathematical Sciences, 2013.Google Scholar
[3]Arun, K.R., Noelle, S., Lukáčová-Medvid’ová, M., and Munz, C.-D.Asymptotic preserving all Mach number scheme for the Euler equations of gas dynamics. submitted to SIAM J. Sci. Comput. 2012.Google Scholar
[4]Audusse, E., Bouchut, F., Bristeau, M., Klein, R., and Perthame, B.A fast and stable well- balanced scheme with hydrostatic reconstruction for shallow water flows. SIAM J. Sci. Com- put., 25(6): 20502065, 2004.CrossRefGoogle Scholar
[5]Bermudez, A. and Vazquez, M.E.Upwind methods for hyperbolic conservation laws with source terms. Comput. Fluids, 23(8): 10491071, 1994.CrossRefGoogle Scholar
[6]Bijl, H. and Wesseling, P.A unified method for computing incompressible and compressible flows in boundary-fitted coordinates. J. Comput. Phys., 141(2): 153173, 1998.CrossRefGoogle Scholar
[7]Botta, N., Klein, R., Langenberg, S., and Lutzenkirchen, S.Well balanced finite volume methods for nearly hydrostatic flows. J. Comput. Phys., 196(2): 539565, 2004.CrossRefGoogle Scholar
[8]Bouchut, F., Le Sommer, J., and Zeitlin, V.Frontal geostrophic adjustment and nonlinear wave phenomena in one-dimensional rotating shallow water. II: High-resolution numerical simulations. J. Fluid Mech., 514: 3563, 2004.CrossRefGoogle Scholar
[9]Bresch, D., Klein, R., and Lucas, C.Multiscale analyses for the shallow water equations. Computational Science and High Performance Computing IV, Notes on Numerical Fluid Mechanics and Multidisciplinary Design 115, 149164, 2011.Google Scholar
[10]Canestrelli, A., Dumbser, M., Siviglia, A., and Toro, E.F.Well-balanced high-order centered schemes on unstructured meshes for shallow water equations with fixed and mobile bed. Advanced in Water Resources, 33: 291303, 2010.CrossRefGoogle Scholar
[11]Castro Diaz, Manuel J., Gonzalez-Vida, J. M., Macias, Jorge, and Pares, Carlos. Realistic application of a tidal 2D two-layer shallow water model to the Strait of Gibraltar. AIP Conf. Proceedings, 168, 14291432, 2009.CrossRefGoogle Scholar
[12]Degond, P. and Tang, M.All speed scheme for the low Mach number limit of the isentropic Euler equations. Commun. Comput. Phys., 10: 131, 2011.CrossRefGoogle Scholar
[13]Feistauer, M., Felcman, J., and Straskraba, I.Mathematical and Computational Methods for Compressible Flow. Oxford University Press, 2003.Google Scholar
[14]Giraldo, F.X. and Restelli, M.High-order semi-implicit time-integrators for a triangular discontinuous Galerkin oceanic shallow water model. Int. J. Numer. Methods Fluids, 63(9): 10771102, 2010.Google Scholar
[15]Giraldo, F.X., Restelli, M., and Lauter, M.Semi-implicit formulations of the Navier-Stokes equations: application to nonhydrostatic atmospheric modeling. J. Sci. Comput., 32(6): 33943425, 2010.Google Scholar
[16]Greenberg, J.M. and Le Roux, A.-Y.A well-balanced scheme for the numerical processing of source terms in hyperbolic equations. SIAM J. Numer. Anal., 33(1): 116, 1996.CrossRefGoogle Scholar
[17]Haack, J., Jin, S., and Liu, J.An all-speed asymptotic-preserving method for the isentropic Euler and Navier-Stokes equations. Commun. Comput. Phys., 12: 955980, 2012.CrossRefGoogle Scholar
[18]Hoffmann, L. B.Ein zeitlich selbstadaptives numerisches Verfahren zur berechnung von Strömungen aller Mach-Zahlen basierend auf Mehrskalenasymptotik und diskreter Datenanalys. PhD thesis, Universitat Hamburg, 2000.Google Scholar
[19]Jin, S.Efficient asymptotic-preserving (AP) schemes for some multiscale kinetic equations. SIAM J. Sci. Comput., 21(2): 441454, 1999.CrossRefGoogle Scholar
[20]Klainerman, S. and Majda, A.Singular limits of quasilinear hyperbolic systems with large parameters and the incompressible limit of compressible fluids. Commun. Pure Appl. Math., 34: 481524, 1981.CrossRefGoogle Scholar
[21]Klein, R., Botta, N., Schneider, T., Munz, C.-D., Roller, S., Meister, A., Hoffmann, L., and Sonar, T.Asymptotic adaptive methods for multi-scale problems in fluid mechanics. J. Eng. Math., 39(1-4): 261343, 2001.CrossRefGoogle Scholar
[22]Klein, R.An applied mathematical view of meteorological modelling. James M., Hill et al. (eds.), Applied mathematics entering the 21st century, Proceedings in Applied Mathematics, 116: 227269, 2004.Google Scholar
[23]Klemp, J.B., Skamarock, W.C., and Dudhia, J.Conservative split-explicit time integration methods for the compressible nonhydrostatic equations. Monthly Weather Rev., 135:28972913,2007.CrossRefGoogle Scholar
[24]Kurganov, A. and Levy, D.Central-upwind schemes for the Saint-Venant system. Math. Model. Numer. Anal., 36(3): 397425, 2002.CrossRefGoogle Scholar
[25]LeVeque, R.J.Balancing source terms and flux gradients in high-resolution Godunov methods: The quasi-steady wave-propagation algorithm. J. Comput. Phys., 146(1): 346365, 1998.CrossRefGoogle Scholar
[26]Liang, Q., Borthwick, A.G.L.Adaptive quadtree simulation of shallow flows with wet-dry fronts over complex topography. Comput. Fluids, 38(2): 221234, 2009.CrossRefGoogle Scholar
[27]Liang, Q., Marche, F.Numerical resolution of well-balanced shallow water equations with complex source terms. Adv. Water Resour., 32(6): 873884, 2009.CrossRefGoogle Scholar
[28]Lukáčová-Medvid’ová, M., Morton, K.W., and Warnecke, G.Finite volume evolution Galerkin methods for hyperbolic systems. J. Sci. Comput., 26(1): 130, 2004.Google Scholar
[29]Lukáčová-Medvid’ová, M., Morton, K.W., and Warnecke, Gerald. Evolution Galerkin meth-ods for hyperbolic systems in two space dimensions. Math. Comput., 69(232):13551384, 2000.CrossRefGoogle Scholar
[30]Lukáčová Medvid’ová, M., Müller, A., Wirth, V., and Yelash, L.Large time step discontinuous evolution Galerkin methods for multiscale atmospheric flow. submitted to J. Comput. Phys., 2013.Google Scholar
[31]Lukáčová-Medvid’ová, M. and Morton, K.W.Finite volume evolution Galerkin methods - A survey. Indian J. Pure Appl. Math., 41(2): 329361, 2010.CrossRefGoogle Scholar
[32]Lukáčová-Medvid’ová, M., Noelle, S., and Kraft, M.Well-balanced finite volume evolution Galerkin methods for the shallow water equations. J. Comput. Phys., 221(1): 122147, 2007.CrossRefGoogle Scholar
[33]Majda, A.Introduction to PDEs and Waves for the Atmosphere and Ocean, Lecture Notes, Courant Institute, 2003.CrossRefGoogle Scholar
[34]Müller, A., Behrens, J., Giraldo, F.X., and Wirth, V.Comparision between adaptive and uni-form discontinuous Galerkin simulations in dry 2d bubble experiments. J. Comput. Phys., 235: 371393, 2013.CrossRefGoogle Scholar
[35]Munz, C.-D., Roller, S., Klein, R., and Geratz, K.J.The extension of incompressible flow solvers to the weakly compressible regime. Comput. Fluids, 32(2): 173196, 2003.CrossRefGoogle Scholar
[36]Park, J.H. and Munz, C.-D.Multiple pressure variables methods for fluid flow at all Mach numbers. Int. J. Numer. Methods Fluids, 49(8): 905931, 2005.CrossRefGoogle Scholar
[37]Prasad, P.Ray theories for hyperbolic waves, kinematical conservation laws (KCL) and applications. Indian J. Pure Appl. Math., 38(5): 467490, 2007.Google Scholar
[38]Restelli, M.Semi-Lagrangian and Semi-Implicit Discontinuous Galerkin Methods for Atmospheric Modeling Applications. PhD thesis, Politecnico di Milano, 2007.Google Scholar
[39]Restelli, M. and Giraldo, F.X.A conservative discontinuous Galerkin semi-implicit formulation for the Navier-Stokes equations in nonhydrostatic mesoscale modeling. SIAM J. Sci. Comput., 31(3): 22312257, 2009.CrossRefGoogle Scholar
[40]Ricchiuto, M. and Bollermann, A.Stabilized residual distribution for shallow water simulations. J. Comput. Phys., 228(4): 10711115, 2009.CrossRefGoogle Scholar
[41]Rogers, B., Borthwick, A.G.L., and Taylor, P.H.Mathematical balancing of flux gradient and source terms prior to using Roe’s approximate Riemann solver. J. Comput. Phys., 192(2):422451, 2003.CrossRefGoogle Scholar
[42]Rogers, B., Fujihara, M., and Borthwick, A.G.L.Adaptive Q-tree Godunov-type scheme for shallow water equations. Int. J. Numer. Methods Fluids, 35(3): 247280, 2001.3.0.CO;2-E>CrossRefGoogle Scholar
[43]Sun, Y. and Ren, Y.The finite volume local evolution Galerkin method for solving the hyperbolic conservation laws. J. Comput. Phys., 228(13): 49454960, 2009.CrossRefGoogle Scholar
[44]Stecca, G., Siviglia, A. and Toro, E. F.A finite volume upwind-biased centred scheme for hyperbolic systems of conservation laws. Application to shallow water equations. Commun. Comput. Phys., 12(4): 11831214, 2012.Google Scholar
[45]Vallis, G.K.Atmospheric and Oceanic Fluid Dynamics. Cambridge University Press, 2006.Google Scholar
Cited by

Send article to Kindle

To send this article to your Kindle, first ensure is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

Note you can select to send to either the or variations. ‘’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

Find out more about the Kindle Personal Document Service.

IMEX Large Time Step Finite Volume Methods for Low Froude Number Shallow Water Flows
Available formats

Send article to Dropbox

To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

IMEX Large Time Step Finite Volume Methods for Low Froude Number Shallow Water Flows
Available formats

Send article to Google Drive

To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

IMEX Large Time Step Finite Volume Methods for Low Froude Number Shallow Water Flows
Available formats

Reply to: Submit a response

Please enter your response.

Your details

Please enter a valid email address.

Conflicting interests

Do you have any conflicting interests? *