Source author record

Philipp Öffner

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

math.NA Numerical Analysis

Catalog footprint

Research overview

13works

2topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2026arXiv

A Structure Preserving Finite Volume Scheme for the Navier-Stokes-Korteweg Equations

We present a semi-discrete finite volume scheme for the local NavierStokes-Korteweg and Euler-Korteweg systems. Our scheme is applicable for equidistant Cartesian meshes in one and two space dimensions. In contrast to other works, which employ, for example, hyperbolic approximations of the equations or auxiliary-variable approaches leading to extended systems, our scheme operates directly on the original system. We prove that it conserves mass and momentum and is energy stable. Numerical experiments complement our theoretical findings, showing that the scheme is convergent of order one if employed with explicit or implicit time discretisation.

preprint2022arXiv

Comparison to control oscillations in high-order Finite Volume schemes via physical constraint limiters, neural networks and polynomial annihilation

The construction of high-order structure-preserving numerical schemes to solve hyperbolic conservation laws has attracted a lot of attention in the last decades and various different ansatzes exist. In this paper, we compare three completely different approaches, i.e. physical constraint limiting, deep neural networks and the application of polynomial annihilation to construct high-order oscillation free Finite Volume (FV) blending schemes. We further analyze their analytical and numerical properties. We demonstrate that all techniques can be used and yield highly efficient FV methods but also come with some additional drawbacks which we point out. Our investigation of the different blending strategies should lead to a better understanding of those techniques and can be transferred to other numerical methods as well which use similar ideas.

preprint2022arXiv

Convergence of Discontinuous Galerkin Schemes for the Euler Equations via Dissipative Weak Solutions

In this paper, we present convergence analysis of high-order finite element based methods, in particular, we focus on a discontinuous Galerkin scheme using summation-by-parts operators. To this end, it is crucial that structure preserving properties, such as positivity preservation and entropy inequality hold. We demonstrate how to ensure them and prove the convergence of our multidimensional high-order DG scheme via dissipative weak solutions. In numerical simulations, we verify our theoretical results.

preprint2022arXiv

Energy-stable global radial basis function methods on summation-by-parts form

Radial basis function methods are powerful tools in numerical analysis and have demonstrated good properties in many different simulations. However, for time-dependent partial differential equations, only a few stability results are known. In particular, if boundary conditions are included, stability issues frequently occur. The question we address in this paper is how provable stability for RBF methods can be obtained. We develop and construct energy-stable radial basis function methods using the general framework of summation-by-parts operators often used in the Finite Difference and Finite Element communities.

preprint2022arXiv

Reinterpretation and Extension of Entropy Correction Terms for Residual Distribution and Discontinuous Galerkin Schemes: Application to Structure Preserving Discretization

For the general class of residual distribution (RD) schemes, including many finite element (such as continuous/discontinuous Galerkin) and flux reconstruction methods, an approach to construct entropy conservative/ dissipative semidiscretizations by adding suitable correction terms has been proposed by Abgrall (J.~Comp.~Phys. 372: pp. 640--666, 2018). In this work, the correction terms are characterized as solutions of certain optimization problems and are adapted to the SBP-SAT framework, focusing on discontinuous Galerkin methods. Novel generalizations to entropy inequalities, multiple constraints, and kinetic energy preservation for the Euler equations are developed and tested in numerical experiments. For all of these optimization problems, explicit solutions are provided. Additionally, the correction approach is applied for the first time to obtain a fully discrete entropy conservative/dissipative RD scheme. Here, the application of the deferred correction (DeC) method for the time integration is essential. This paper can be seen as describing a systematic method to construct structure preserving discretization, at least for the considered example.

preprint2021arXiv

Towards Stable Radial Basis Function Methods for Linear Advection Problems

In this work, we investigate (energy) stability of global radial basis function (RBF) methods for linear advection problems. Classically, boundary conditions (BC) are enforced strongly in RBF methods. By now it is well-known that this can lead to stability problems, however. Here, we follow a different path and propose two novel RBF approaches which are based on a weak enforcement of BCs. By using the concept of flux reconstruction and simultaneous approximation terms (SATs), respectively, we are able to prove that both new RBF schemes are strongly (energy) stable. Numerical results in one and two spatial dimensions for both scalar equations and systems are presented, supporting our theoretical analysis.

preprint2020arXiv

General polytopal H(div) conformal finite elements and their discretisation spaces

We present a class of discretisation spaces and H(div)-conformal elements that can be built on any polytope. Bridging the flexibility of the Virtual Element spaces towards the element's shape with the divergence properties of the Raviart-Thomas elements on the boundaries, the designed frameworks offer a wide range of H(div)-conformal discretisations. As those elements are set up through degrees of freedom, their definitions are easily amenable to the properties the approximated quantities are wished to fulfill. Furthermore, we show that one straightforward restriction of this general setting share its properties with the classical Raviart-Thomas elements at each interface, for any order and any polytopial shape. Then, to close the introduction of those new elements by an example, we investigate the shape of the basis functions corresponding to particular elements in the two dimensional case.

preprint2019arXiv

Analysis of Artificial Dissipation of Explicit and Implicit Time-Integration Methods

Stability is an important aspect of numerical methods for hyperbolic conservation laws and has received much interest. However, continuity in time is often assumed and only semidiscrete stability is studied. Thus, it is interesting to investigate the influence of explicit and implicit time integration methods on the stability of numerical schemes. If an explicit time integration method is applied, spacially stable numerical schemes for hyperbolic conservation laws can result in unstable fully discrete schemes. Focusing on the explicit Euler method (and convex combinations thereof), undesired terms in the energy balance trigger this phenomenon and introduce an erroneous growth of the energy over time. In this work, we study the influence of artificial dissipation and modal filtering in the context of discontinuous spectral element methods to remedy these issues. In particular, lower bounds on the strength of both artificial dissipation and modal filtering operators are given and an adaptive procedure to conserve the (discrete) $\mathbf{L}_2$ norm of the numerical solution in time is derived. This might be beneficial in regions where the solution is smooth and for long time simulations. Moreover, this approach is used to study the connections between explicit and implicit time integration methods and the associated energy production. By adjusting the adaptive procedure, we demonstrate that filtering in explicit time integration methods is able to mimic the dissipative behavior inherent in implicit time integration methods. This contribution leads to a better understanding of existing algorithms and numerical techniques, in particular the application of artificial dissipation as well as modal filtering in the context of numerical methods for hyperbolic conservation laws together with the selection of explicit or implicit time integration methods.

preprint2019arXiv

Arbitrary high-order, conservative and positive preserving Patankar-type deferred correction schemes

Production-destruction systems (PDS) of ordinary differential equations (ODEs) are used to describe physical and biological reactions in nature. The considered quantities are subject to natural laws. Therefore, they preserve positivity and conservation of mass at the analytical level. In order to maintain these properties at the discrete level, the so-called modified Patankar-Runge-Kutta (MPRK) schemes are often used in this context. However, up to our knowledge, the family of MPRK has been only developed up to third order of accuracy. In this work, we propose a method to solve PDS problems, but using the Deferred Correction (DeC) process as a time integration method. Applying the modified Patankar approach to the DeC scheme results in provable conservative and positivity preserving methods. Furthermore, we demonstrate that these modified Patankar DeC schemes can be constructed up to arbitrarily high order. Finally, we validate our theoretical analysis through numerical simulations.

preprint2016arXiv

Enhancing stability of correction procedure via reconstruction using summation-by-parts operators I: Artificial dissipation

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.

preprint2016arXiv

Enhancing stability of correction procedure via reconstruction using summation-by-parts operators II: Modal filtering

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.

preprint2016arXiv

On Hahn polynomial expansion of a continuous function of bounded variation

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]$ .

preprint2016arXiv

Spectral accuracy for the Hahn polynomials

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.

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.

Philipp Öffner

Research overview

Is this your work?

13 publications

A Structure Preserving Finite Volume Scheme for the Navier-Stokes-Korteweg Equations

Comparison to control oscillations in high-order Finite Volume schemes via physical constraint limiters, neural networks and polynomial annihilation

Convergence of Discontinuous Galerkin Schemes for the Euler Equations via Dissipative Weak Solutions

Energy-stable global radial basis function methods on summation-by-parts form

Reinterpretation and Extension of Entropy Correction Terms for Residual Distribution and Discontinuous Galerkin Schemes: Application to Structure Preserving Discretization

Towards Stable Radial Basis Function Methods for Linear Advection Problems

General polytopal H(div) conformal finite elements and their discretisation spaces

Analysis of Artificial Dissipation of Explicit and Implicit Time-Integration Methods

Arbitrary high-order, conservative and positive preserving Patankar-type deferred correction schemes

Enhancing stability of correction procedure via reconstruction using summation-by-parts operators I: Artificial dissipation

Enhancing stability of correction procedure via reconstruction using summation-by-parts operators II: Modal filtering

On Hahn polynomial expansion of a continuous function of bounded variation

Spectral accuracy for the Hahn polynomials

See the researcher in context

Building this map preview