Philipp Öffner

Hendrik Ranocha, Philipp Öffner, Thomas Sonar

The correction procedure via reconstruction (CPR, formerly known as flux reconstruction) is a framework of high order methods for conservation laws, unifying some discontinuous Galerkin, spectral difference and spectral volume methods. Linearly stable schemes were presented by Vincent et al. (2011, 2015), but proofs of non-linear (entropy) stability in this framework have not been published yet (to the knowledge of the authors). We reformulate CPR methods using summation-by-parts (SBP) operators with simultaneous approximation terms (SATs), a framework popular for finite difference methods, extending the results obtained by Gassner (2013) for a special discontinuous Galerkin spectral element method. This reformulation leads to proofs of conservation and stability in discrete norms associated with the method, recovering the linearly stable CPR schemes of Vincent et al. (2011, 2015). Additionally, extending the skew-symmetric formulation of conservation laws by additional correction terms, entropy stability for Burgers' equation is proved for general SBP CPR methods not including boundary nodes.

Hendrik Ranocha, Philipp Öffner, Thomas Sonar

A generalised analytical notion of summation-by-parts (SBP) methods is proposed, extending the concept of SBP operators in the correction procedure via reconstruction (CPR), a framework of high-order methods for conservation laws. For the first time, SBP operators with dense norms and not including boundary points are used to get an entropy stable split-form of Burgers' equation. Moreover, overcoming limitations of the finite difference framework, stability for curvilinear grids and dense norms is obtained for SBP CPR methods by using a suitable way to compute the Jacobian.

Philipp Öffner, Jan Glaubitz, Hendrik Ranocha

In this paper, we consider Burgers' equation with uncertain boundary and initial conditions. The polynomial chaos (PC) approach yields a hyperbolic system of deterministic equations, which can be solved by several numerical methods. Here, we apply the correction procedure via reconstruction (CPR) using summation-by-parts operators. We focus especially on stability, which is proven for CPR methods and the systems arising from the PC approach. Due to the usage of split-forms, the major challenge is to construct entropy stable numerical fluxes. For the first time, such numerical fluxes are constructed for all systems resulting from the PC approach for Burgers' equation. In numerical tests, we verify our results and show also the advantage of the given ansatz using CPR methods. Moreover, one of the simulations, i.e. Burgers' equation equipped with an initial shock, demonstrates quite fascinating observations. The behaviour of the numerical solutions from several methods (finite volume, finite difference, CPR) differ significantly from each other. Through careful investigations, we conclude that the reason for this is the high sensitivity of the system to varying dissipation. Furthermore, it should be stressed that the system is not strictly hyperbolic with genuinely nonlinear or linearly degenerate fields.

Philipp Öffner, Hendrik Ranocha

For practical applications, the long time behaviour of the error of numerical solutions to time-dependent partial differential equations is very important. Here, we investigate this topic in the context of hyperbolic conservation laws and flux reconstruction schemes, focusing on the schemes in the discontinuous Galerkin spectral element framework. For linear problems with constant coefficients, it is well-known in the literature that the choice of the numerical flux (e.g. central or upwind) and the selection of the polynomial basis (e.g. Gauß-Legendre or Gauß-Lobatto-Legendre) affects both the growth rate and the asymptotic value of the error. Here, we extend these investigations of the long time error to variable coefficients using both Gauß-Lobatto-Legendre and Gauß-Legendre nodes as well as several numerical fluxes. We derive conditions guaranteeing that the errors are still bounded in time. Furthermore, we analyse the error behaviour under these conditions and demonstrate in several numerical tests similarities to the case of constant coefficients. However, if these conditions are violated, the error shows a completely different behaviour. Indeed, by applying central numerical fluxes, the error increases without upper bound while upwind numerical fluxes can still result in uniformly bounded numerical errors. An explanation for this phenomenon is given, confirming our analytical investigations.

Jan Glaubitz, Philipp Öffner, Thomas Sonar

We adapt the spectral viscosity (SV) formulation implemented as a modal filter to a Spectral Difference Method (SD) solving hyperbolic conservation laws. In the SD Method we use selections of different orthogonal polynomials (APK polynomials). Furthermore we obtain new error bounds for filtered APK extensions of smooth functions. We demonstrate that the modal filter also depends on the chosen polynomial basis in the SD Method. Spectral filtering stabilizes the scheme and leaves weaker oscillations. Hence, the selection of the family of orthogonal polynomials on triangles and their specific modal filter possesses a positive influence on the stability and accuracy of the SD Method. In the second part, we initiate a stability analysis for a linear scalar test case with periodic initial condition to find the best selection of APK polynomials and their specific modal filter. To the best of our knowledge, this work is the first that gives a stability analysis for a scheme with spectral filtering. Finally, we demonstrate the influence of the underlying basis of APK polynomials in a well-known test case.

Hendrik Ranocha, Jan Glaubitz, Philipp Öffner et al.

The correction procedure via reconstruction (CPR, also known as flux reconstruction) is a framework of high order semidiscretisations used for the numerical solution of hyperbolic conservation laws. Using a reformulation of these schemes relying on summation-by-parts (SBP) operators and simultaneous approximation terms (SATs), artificial dissipation / spectral viscosity operators are investigated in this first part of a series. Semidiscrete stability results for linear advection and Burgers' equation as model problems are extended to fully discrete stability by an explicit Euler method. As second part of this series, Glaubitz, Ranocha, Öffner, and Sonar (Enhancing stability of correction procedure via reconstruction using summation-by-parts operators II: Modal filtering, 2016) investigate connections to modal filters and their application instead of artificial dissipation.

Jan Glaubitz, Hendrik Ranocha, Philipp Öffner et al.

A recently introduced framework of semidiscretisations for hyperbolic conservation laws known as correction procedure via reconstruction (CPR, also known as flux reconstruction) is considered in the extended setting of summation-by-parts (SBP) operators using simultaneous approximation terms (SATs). This reformulation can yield stable semidiscretisations for linear advection and Burgers' equation as model problems. In order to enhance these properties, modal filters are introduced to this framework. As a second part of a series, the results of Ranocha, Glaubitz, Öffner, and Sonar ("Enhancing stability of correction procedure via reconstruction using summation-by-parts operators I: Artificial dissipation", 2016) concerning artificial dissipation / spectral viscosity are extended, yielding fully discrete stable schemes. Additionally, a new adaptive strategy to compute the filter strength is introduced and different possible applications of modal filters are compared both theoretically and numerically.

Rémi Abgrall, Philipp Öffner, Yongle Liu

In this paper, we provide a few new properties of Active Flux (AF)/Point-Average-Moment PolynomiAl-interpreted (\pampa) schemes. First, we show, in full generality, that the AF/pampa schemes can be interpreted in such a way that the discontinuous Galerkin (dG) scheme is one of their building blocks. Secondly we provide intrinsic bound preserving properties of the current variant of pampa. This is also illustrated numerically. Last, we show, at least in one dimension, that the pampa scheme has the summation by part (SBP) property.

Jaya Agnihotri, Philipp Öffner

We present a convergence analysis of a finite volume (FV) scheme for the multicomponent compressible Euler system in the framework of dissipative weak (DW) solutions. DW solutions were introduced as a generalized solution framework in computational fluid dynamics and have recently gained considerable attention. They extend the well-known Lax Equivalence Theorem to nonlinear settings, meaning that if a numerical scheme is both consistent and stable, it will also converge. The FV scheme under consideration preserves key physical properties of the fluid mixture, in particular, positivity of partial densities, pressure, and temperature. Using uniform stability bounds and consistency estimates, we prove that the numerical solutions converge in the framework of DW solutions of the multicomponent Euler system. Applying the relative entropy method and the weak-strong uniqueness principle, we further show that the approximate solutions converge strongly to the classical solution as long as it exists. Numerical experiments confirm the theoretical results, not only for low-order FV methods but also through extended numerical investigations of a higher-order, structure-preserving discontinuous Galerkin scheme.

Jan Giesselmann, Philipp Öffner, Robert Sauerborn

We propose a concept of dissipative weak (DW) solutions for the Navier-Stokes-Korteweg (NSK) system and prove conditional convergence of a structure-preserving finite volume scheme towards such a solution. DW solutions provide a generalized solution concept in computational fluid dynamics and have recently attracted significant attention. They provide an extension of the famous Lax Equivalence Theorem to nonlinear problems, i.e. consistency and stability of a numerical scheme imply convergence. Our work builds on recent advances where convergence towards DW solutions of structure-preserving schemes has been established for the Euler and Navier-Stokes equations. Indeed, we prove convergence of a recently proposed FV scheme by leveraging its conservation and dissipation properties as well as its consistency.

René Goertz, Philipp Öffner

We consider the well-known method of least squares on an equidistant grid with $N+1$ nodes on the interval $[-1,1]$. We investigate the following problem: For which ratio $N/n$ and which functions, do we have pointwise convergence of the least square operator ${LS}_n^N:\mathcal{C}\left[-1,1\right]\rightarrow\mathcal{P}_n$? To solve this problem we investigate the relation between the Jacobi polynomials $P_k^{α,β}$ and the Hahn polynomials $Q_k\left(\cdot;α,β,N\right)$. Thereby we describe the least square operator ${LS}_n^N$ by the expansion of a function by Hahn polynomials. In particular we present the following result: The series expansion $\sum_{k=0}^n{\hat{f} Q_k}$ of a function $f$ by Hahn polynomials $Q_k$ converges pointwise, if the series expansion $\sum_{k=0}^n{\hat{f} P_k}$ of the function $f$ by Jacobi polynomials $P_k$ converges pointwise and if ${n^4}/N\rightarrow 0$ for $n,N\rightarrow\infty$. Furthermore we obtain the following result: Let $f\in\left\{g\in\mathcal{C}^1\left[-1,1\right]:g^\prime\in\mathcal{BV}\left[-1,1\right]\right\}$ and let $(N_n)_{n}$ be a sequence of natural numbers with ${n^4}/{N_n}\rightarrow 0$. Then the least square method ${LS}_n^{N_n}[f]$ converges for each $x\in[-1,1]$.

René Goertz, Philipp Öffner

We consider in this paper the Hahn polynomials and their application in numerical methods. The Hahn polynomials are classical discrete orthogonal polynomials. We analyse the behaviour of these polynomials in the context of spectral approximation of partial differential equations. We study series expansions $u=\sum_{n=0}^\infty \hat{u}_n ϕ_n$, where the $ϕ_n$ are the Hahn polynomials. We examine the Hahn coefficients and proof spectral accuracy in some sense. We substantiate our results by numericals tests. Furthermore we discuss a problem which arise by using the Hahn polynomials in the approximation of a function $u$, which is linked to the Runge phenomenon. We suggest two approaches to avoid this problem. These will also be the motivation and the outlook of further research in the application of discrete orthogonal polynomials in a spectral method for the numerical solution of hyperbolic conservation laws.