Sebastian Noelle

Yulong Xing, Chi-Wang Shu, Sebastian Noelle

This note aims at demonstrating the advantage of moving-water well-balanced schemes over still-water well-balanced schemes for the shallow water equations. We concentrate on numerical examples with solutions near a moving-water equilibrium. For such examples, still-water well-balanced methods are not capable of capturing the small perturbations of the moving-water equilibrium and may generate significant spurious oscillations, unless an extremely refined mesh is used. On the other hand, moving- water well-balanced methods perform well in these tests. The numerical examples in this note clearly demonstrate the importance of utilizing moving-water well-balanced methods for solutions near a moving-water equilibrium.

Sebastian Noelle, Michael Westdickenberg

In this paper we consider convergence of approximate solutions of conservation laws. We start with an overview over the historical developments since the 1950s, and the analytical tools used in this context. Then we present some of our own results on the convergence of numerical approximations, discuss recent related work and open problems.

Andreas Bollermann, Sebastian Noelle, Maria Lukáčová - Medvidová

We present a new Finite Volume Evolution Galerkin (FVEG) scheme for the solution of the shallow water equations (SWE) with the bottom topography as a source term. Our new scheme will be based on the FVEG methods presented in (Lukáčová, Noelle and Kraft, J. Comp. Phys. 221, 2007), but adds the possibility to handle dry boundaries. The most important aspect is to preserve the positivity of the water height. We present a general approach to ensure this for arbitrary finite volume schemes. The main idea is to limit the outgoing fluxes of a cell whenever they would create negative water height. Physically, this corresponds to the absence of fluxes in the presence of vacuum. Well-balancing is then re-established by splitting gravitational and gravity driven parts of the flux. Moreover, a new entropy fix is introduced that improves the reproduction of sonic rarefaction waves.

Maria Lukáčová - Medvidová, Sebastian Noelle, Marcus Kraft

We present a new well-balanced finite volume method within the framework of the finite volume evolution Galerkin (FVEG) schemes. The methodology will be illustrated for the shallow water equations with source terms modelling the bottom topography and Coriolis forces. Results can be generalized to more complex systems of balance laws. The FVEG methods couple a finite volume formulation with approximate evolution operators. The latter are constructed using the bicharacteristics of multidimensional hyperbolic systems, such that all of the infinitely many directions of wave propagation are taken into account explicitly. We derive a well- balanced approximation of the integral equations and prove that the FVEG scheme is well-balanced for the stationary steady states as well as for the steady jets in the rotational frame. Several numerical experiments for stationary and quasi-stationary states as well as for steady jets confirm the reliability of the well-balanced FVEG scheme.

Sebastian Noelle, Wolfram Rosenbaum, Martin Rumpf

Central schemes are frequently used for incompressible and compressible flow calculations. The present paper is the first in a forthcoming series where a new approach to a 2nd order accurate Finite Volume scheme operating on cartesian grids is discussed. Here we start with an adaptively refined cartesian primal grid in 3D and present a construction technique for the staggered dual grid based on $L^{\infty}$-Voronoi cells. The local refinement constellation on the primal grid leads to a finite number of uniquely defined local patterns on a primal cell. Assembling adjacent local patterns forms the dual grid. All local patterns can be analysed in advance. Later, running the numerical scheme on staggered grids, all necessary geometric information can instantly be retrieved from lookup-tables. The new scheme is compared to established ones in terms of algorithmical complexity and computational effort.

Andreas Bollermann, Guoxian Chen, Alexander Kurganov et al.

In this paper, we construct a well-balanced, positivity preserving finite volume scheme for the shallow water equations based on a continuous, piecewise linear discretization of the bottom topography. The main new technique is a special reconstruction of the flow variables in wet-dry cells, which is presented in this paper for the one dimensional case. We realize the new reconstruction in the framework of the second-order semi-discrete central-upwind scheme from (A. Kurganov and G. Petrova, Commun. Math. Sci., 2007). The positivity of the computed water height is ensured following (A. Bollermann, S. Noelle and M. Lukáčová-Medviďová, Commun. Comput. Phys., 2010): The outgoing fluxes are limited in case of draining cells.

Normann Pankratz, Jostein R. Natvig, Bjørn Gjevik et al.

In this paper we compare a classical finite-difference and a high order finite- volume scheme for barotropic ocean flows. We compare the schemes with respect to their accuracy, stability, and study various outflow and inflow boundary conditions. We apply the schemes to the problem of eddy formation in shelf slope jets along the Ormen Lange section of the Norwegian shelf. Our results strongly confirm the development of mesoscale eddies caused by instability of the flows.

Jochen Schütz, Sebastian Noelle

For low Mach number flows, there is a strong recent interest in the development and analysis of IMEX (implicit/explicit) schemes, which rely on a splitting of the convective flux into stiff and nonstiff parts. A key ingredient of the analysis is the so-called Asymptotic Preserving (AP) property, which guarantees uniform consistency and stability as the Mach number goes to zero. While many authors have focussed on asymptotic consistency, we study asymptotic stability in this paper: does an IMEX scheme allow for a CFL number which is independent of the Mach number? We derive a stability criterion for a general linear hyperbolic system. In the decisive eigenvalue analysis, the advective term, the upwind diffusion and a quadratic term stemming from the truncation in time all interact in a subtle way. As an application, we show that a new class of splittings based on characteristic decomposition, for which the commutator vanishes, avoids the deterioration of the time step which has sometimes been observed in the literature.

Sebastian Noelle, Georgij Binev, K. R. Arun et al.

We propose a low Mach number, Godunov-type finite volume scheme for the numerical solution of the compressible Euler equations of gas dynamics. The scheme combines Klein's non-stiff/stiff decomposition of the fluxes (J. Comput. Phys. 121:213-237, 1995) with an explicit/implicit time discretization (Cordier et al., J. Comput. Phys. 231:5685- 5704, 2012) for the split fluxes. This results in a scalar second order partial differential equation (PDE) for the pressure, which we solve by an iterative approximation. Due to our choice of a crucial reference pressure, the stiff subsystem is hyperbolic, and the second order PDE for the pressure is elliptic. The scheme is also uniformly asymptotically consistent. Numerical experiments show that the scheme needs to be stabilized for low Mach numbers. Unfortunately, this affects the asymptotic consistency, which becomes non-uniform in the Mach number, and requires an unduly fine grid in the small Mach number limit. On the other hand, the CFL number is only related to the non-stiff characteristic speeds, independently of the Mach number. Our analytical and numerical results stress the importance of further studies of asymptotic stability in the development of AP (asymptotic preserving) schemes.