Philipp Öffner

Thomas Eiter, Jan Giesselmann, Robert Lasarzik et al.

We consider a structure-preserving finite-volume scheme for the Euler-Korteweg (EK) and Navier-Stokes-Korteweg (NSK) equations. We prove that its numerical solutions converge to energy-variational solutions of EK or NSK under mesh refinement. Energy-variational solutions constitute a novel solution concept that has recently been introduced for hyperbolic conservation laws, including the EK system, and which we extend to the NSK model. Our proof is based on establishing uniform estimates following from the properties of the structure-preserving scheme, and using the stability of the energy-variational formulation under weak convergence in the natural energy spaces.

Philipp Öffner, Per Pettersson, Andrew R. Winters

Recently, two independent research efforts have been made to study the stochastic Galerkin formulation of the shallow water equations. %In particular, Bender and Öffner developed entropy-conservative discontinuous Galerkin (DG) methods to solve the stochastic shallow water equations in an stochastic Galerkin framework using Roe variable transformation, while Dai, Epshteyn and collaborators proposed second-order, energy-stable and well-balanced schemes for the same class of problems with a specific projection step used inside the Galerkin projection together with high-order quadrature rules and a time-step restriction. In this paper, we provide a comprehensive comparison of the two methodologies mentioned, focusing on their theoretical properties and practical implementation aspects. We highlight shared foundational concepts and key differences of both approaches, with a particular focus on the selection of basis functions in the stochastic domain. As a highlight, we show that under specific conditions, the two formulations align, offering a unified framework that connects these distinct approaches. From our theoretical findings, we extend the development of high-order entropy conservative DG methods for the one-dimensional stochastic Galerkin shallow equations to two space dimensions; constructing entropy conservative two-point fluxes via primitive variables instead of entropy variables and applying it in our high-order DG setting. In numerical simulations, we verify and support our theoretical findings of a well-balanced and entropy-stable DG scheme which can be used to solve geophyiscal fluid flows with uncertainty.