Coupling OpenFOAM(R) with BoSSS, a discontinuous Galerkin solver written in C#

TL;DR

This paper introduces foam-dg, a flexible coupling framework that integrates OpenFOAM with the high-order DG solver BoSSS, enabling mixed simulations of complex PDE systems.

Contribution

It presents a novel coupling approach between OpenFOAM and BoSSS, allowing for mixed high-order and finite volume calculations within a single framework.

Findings

Successful coupling demonstrated with convective Cahn-Hilliard equation

Mixed calculations are feasible and produce reasonable results

Framework enables flexible multi-method PDE solving

Abstract

In this article, we present the foam-dg project, which provides a bridge between OpenFOAM(R) and the high-order DG (discontinuous Galerkin) framework BoSSS. Thanks to the flexibility of the coupling approach, mixed calculations where some parts of the equation system are solved in OpenFOAM(R) and others are solved in BoSSS are easily possible. This is showcased using the convective Cahn-Hilliard equation, where the Cahn-Hilliard part is solved in BoSSS and the Navier-Stokes part is solved in OpenFOAM(R). The obtained results appear reasonable, though the main focus of this paper is to present and document the foam-dg project rather than on quantitative results.