BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Sebastian Noelle

Sebastian Noelle appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysis

Catalog footprint

What is connected

14works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

A class of boundary conditions for time-discrete Green-Naghdi equations with bathymetry

This work is devoted to the structure of the time-discrete Green-Naghdi equations including bathymetry. We use the projection structure of the equations to characterize homogeneous and inhomogeneous boundary conditions for which the semi-discrete equations are well-posed. This structure allows us to propose efficient and robust numerical treatment of the boundary conditions that ensures entropy stability of the scheme by construction. Numerical evidence is provided to illustrate that our approach is suitable for situations of practical interest that are not covered by existing theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

3D Adaptive Central Schemes: part I Algorithms for Assembling the Dual Mesh

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite Volume Evolution Galerkin Methods for the Shallow Water Equations with Dry Beds

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the advantage of well-balanced schemes for moving-water equilibria of the shallow water equations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Well-balanced finite volume evolution Galerkin methods for the shallow water equations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A note on adjoint error estimation for one-dimensional stationary balance laws with shocks

We consider one-dimensional steady-state balance laws with discontinuous solutions. Giles and Pierce realized that a shock leads to a new term in the adjoint error representation for target functionals.This term disappears if and only if the adjoint solution satisfies an internal boundary condition. Curiously, most computer codes implementing adjoint error estimation ignore the new term in the functional, as well as the internal adjoint boundary condition. The purpose of this note is to justify this omission as follows: if one represents the exact forward and adjoint solutions as vanishing viscosity limits of the corresponding viscous problems, then the internal boundary condition is naturally satisfied in the limit.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Weakly Asymptotic Preserving Low Mach Number Scheme for the Euler Equations of Gas Dynamics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A well-balanced reconstruction for wetting/drying fronts

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Adaptive timestep control for nonstationary solutions of the Euler equations

In this paper we continue our work on adaptive timestep control for weakly non- stationary problems. The core of the method is a space-time splitting of adjoint error representations for target functionals due to Süli and Hartmann. The main new ingredients are (i) the extension from scalar, 1D, conservation laws to the 2D Euler equations of gas dynamics, (ii) the derivation of boundary conditions for a new formulation of the adjoint problem and (iii) the coupling of the adaptive time-stepping with spatial adaptation. For the spatial adaptation, we use a multiresolution-based strategy developed by Müller, and we combine this with an implicit time discretization. The combined space-time adaptive method provides an efficient choice of timesteps for implicit computations of weakly nonstationary flows. The timestep will be very large in time intervalls of stationary flow, and becomes small when a perturbation enters the flow field. The efficiency of the solver is investigated by means of an unsteady inviscid 2D flow over a bump.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Flux Splitting for stiff equations: A notion on stability

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-order well-balanced finite-volume schemes for barotropic flows. Development and numerical comparisons

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On adaptive timestepping for weakly instationary solutions of hyperbolic conservation laws via adjoint error control

We study a recent timestep adaptation technique for hyperbolic conservation laws. The key tool is a space-time splitting of adjoint error representations for target functionals due to Süli and Hartmann. It provides an efficient choice of timesteps for implicit computations of weakly instationary flows. The timestep will be very large in regions of stationary flow, and become small when a perturbation enters the flow field. Besides using adjoint techniques which are already well-established, we also add a new ingredient which simplifies the computation of the dual problem. Due to Galerkin orthogonality, the dual solution ϕ does not enter the error representation as such. Instead, the relevant term is the difference of the dual solution and its projection to the finite element space, ϕ-ϕh . We can show that it is therefore sufficient to compute the spatial gradient of the dual solution, $w = {\nabla} ϕ$. This gradient satisfies a conservation law instead of a transport equation, and it can therefore be computed with the same algorithm as the forward problem, and in the same finite element space. We demonstrate the capabilities of the approach for a weakly instationary test problem for scalar conservation laws.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Timestep control for weakly instationary flows

We report on recent work on adaptive timestep control for weakly instationary gas flows [16, 18, 17] carried out within SFB 401, TPA3. The method which we implement and extend is a space-time splitting of adjoint error representations for target functionals due to Süli [19] and Hartmann [10]. In this paper, we first review the method for scalar, 1D, conservation laws. We design a test problem for weakly instationary solutions and show numerical experiments which clearly show the possible benefits of the method. Then we extend the approach to the 2D Euler equations of gas dynamics. New ingredients are (i) a conservative formulation of the adjoint problem which makes its solution robust and efficient, (ii) the derivation of boundary conditions for this new formulation of the adjoint problem and (iii) the coupling of the adaptive time-stepping with the multiscale spatial adaptation due to M üller [12, 3], also developed within SFB 401. The combined space-time adaptive method provides an efficient choice of timesteps for implicit computations of weakly instationary flows. The timestep will be very large in regions of stationary flow, and becomes small when a perturbation enters the flow field. The efficiency of the Euler solver is investigated by means of an unsteady inviscid 2D flow over a bump.

0 citations0 reviewsNo code linkedOpen access