A multigrid-based multiscale numerical scheme for shallow water flows at low froude number

Abstract

A new multiscale semi-implicit scheme for the computation of low Froude number shallow water flows is presented. Motivated by the needs of atmospheric flow applications, it aims to minimize dispersion and amplitude errors in the computation of long-wave gravity waves. While it correctly balances "slaved" dynamics of short-wave solution components induced by slow forcing, the method eliminates freely propagating compressible short-wave modes, which are under- resolved in time. This is achieved through a multilevel approach borrowing ideas from multigrid schemes for elliptic equations. The scheme is second- order accurate and admits time steps depending essentially on the flow velocity. First, a multilevel method is derived for the one-dimensional linearized shallow water equations. Scale-wise decomposition of the data enables a scale-dependent blending of time integrators with different principal features. To guide the selection of these integrators, the discrete- dispersion relations of some standard second-order schemes are analyzed, and their response to high-wave-number low-frequency source terms is discussed. The resulting method essentially consists of the solution of a Helmholtz problem on the original fine grid, where the differencing operator and the right hand side incorporate the multiscale information of the discretization. The performance of the new multilevel method is illustrated on a test case with "multiscale" initial data and a problem with a slowly varying high-wave- number source term. The scheme for simulating fully nonlinear shallow water flows is a generalization of a projection method for the zero Froude number equations. Therefore, the method described in Vater and Klein (Numer. Math. 113, pp. 123-161, 2009) is extended to account for time dependent bottom topography. Numerical simulations show that the method is well-balanced in that it can reproduce the steady state of a lake at rest with non-trivial bottom topography. Further results verify the correct representation of time dependent bottom topography. Finally, the multiscale semi-implicit method for the nonlinear shallow water equations at low Froude numbers is derived. This incorporates two different extensions of the aforementioned projection method by incorporating the local time derivatives of the height. The two variants are combined in the multiscale method, which incorporates the solution of a multilevel Helmholtz problem, similar to the linear case. The method is implemented in one space dimension. The convergence of the method is analyzed by means of different test cases. Moreover, the method's balancing properties are addressed for the lake at rest, and the asymptotic regime of fast gravity waves traveling over short-range topography is considered. The numerical results of the proposed method suggest that the scheme correctly reproduces this regime, and can be therefore considered as a so-called asymptotically adaptive numerical method.

In dieser Arbeit wird ein neues semi-implizites Mehrskalenverfahren zur Berechnung von Flachwasserströmungen bei kleinen Froudezahlen vorgestellt. Motiviert durch meteorologische Anwendungen zielt es darauf ab, Dispersions- und Amplitudenfehler bei der Berechnung von langwelligen Schwerewellen zu minimieren. Während die durch langsame Anregung aufgezwungene Dynamik in kurzwelligen Lösungskomponenten korrekt balanciert wird, unterdrückt das Verfahren sich frei ausbreitende kurzwellige Moden, die in der Zeit nicht aufgelöst sind. Dies wird durch einen Multilevelansatz erreicht, der sich Konzepten bedient, die von Mehrgitterverfahren zur Lösung von elliptischen Gleichungen bekannt sind. Das Verfahren ist von zweiter Ordnung und erlaubt Zeitschritte, die nur von der Strömungsgeschwindigkeit abhängen. In einem ersten Schritt wird ein Multilevelverfahren zur Lösung der eindimensionalen linearisierten Flachwassergleichungen hergeleitet. Die Zerlegung der Daten nach verschiedenen Skalen ermöglicht eine skalenabhängige Verschneidung von Zeitintegratoren mit verschiedenen wesentlichen Merkmalen. Zur Auswahl der Integratoren werden die diskreten Dispersionsrelationen einiger klassischer Zweite-Ordnung-Verfahren berechnet und deren Verhalten im Falle von niederfrequenten räumlich stark variierenden Quelltermen diskutiert. Das resultierende Verfahren besteht im Wesentlichen aus der Lösung eines Helmholtzproblems auf dem ursprünglichen feinen Gitter, wobei der Differenzenoperator und die rechte Seite die Mehrskaleninformation der Diskretisierung beinhalten. Die Güte des Verfahrens wird in einem Testfall mit "mehrskaligen" Anfangsdaten analysiert. Ein weiterer Testfall gibt Aufschluss über das Verhalten in Anwesenheit eines sich langsam in der Zeit ändernden Quellterms, der räumlich stark variiert. Das Verfahren zur Simulation der vollen nichtlinearen Flachwassergleichungen baut auf ein Projektionsverfahren für die Grenzgleichungen im Limes Froudezahl gegen Null auf. Dafür wird die in Vater und Klein (Numer. Math. 113, S. 123-161, 2009) beschriebene Methode auf Probleme mit zeitlich variierender Bodentopographie erweitert. Numerische Simulationen zeigen, dass das Verfahren gut-balanciert ist und einen ruhenden See mit nichttrivialer Bodentopographie berechnen kann. Weitere Ergebnisse bestätigen die korrekte Darstellung von zeitabhängigen Bodentopographie. In einem letzten Schritt wird das semi-implizite Mehrskalenverfahren für die numerische Lösung der Flachwassergleichungen bei kleinen Froudezahlen hergeleitet. Dies beinhaltet zwei verschiedene Erweiterungen des oben genannten Projektionsverfahrens durch Berücksichtigung der lokalen Zeitableitungen des Höhenfeldes. Im Mehrskalenverfahren werden diese beiden Varianten kombiniert. Dies führt zur Lösung eines Multilevel-Helmholtzproblems analog zum linearen Fall. Das Verfahren wird für den eindimensionalen Fall implementiert, und seine Konvergenzeigenschaften werden anhand verschiedener Testfälle untersucht. Außerdem werden die Balancierungseigenschaften bzgl. eines ruhenden Sees getestet und das asymptotische Regime schneller Schwerewellen, die über kurzwellige Bodentopographie wandern, betrachtet. Die numerischen Ergebnisse des Verfahrens lassen darauf schließen, dass dieses Regime korrekt reproduziert wird, und es sich somit um ein sogenanntes asymptotisch adaptives numerisches Verfahren handelt.