Abstract:We introduce a new family of high order accurate semi-implicit schemes for the solution of non-linear hyperbolic partial differential equations on unstructured polygonal meshes. The time discretization is based on a splitting between explicit and implicit terms that may arise either from the multi-scale nature of the governing equations, which involve both slow and fast scales, or in the context of projection methods, where the numerical solution is projected onto the physically meaningful solution manifold. We propose to use a high order finite volume (FV) scheme for the explicit terms, ensuring conservation property and robustness across shock waves, while the virtual element method (VEM) is employed to deal with the discretization of the implicit terms, which typically requires an elliptic problem to be solved. The numerical solution is then transferred via suitable L2 projection operators from the FV to the VEM solution space and vice-versa. High order time accuracy is achieved using the semi-implicit IMEX Runge-Kutta schemes, and the novel schemes are proven to be asymptotic preserving and well-balanced. As representative models, we choose the shallow water equations (SWE), thus handling multiple time scales characterized by a different Froude number, and the incompressible Navier-Stokes equations (INS), which are solved at the aid of a projection method to satisfy the solenoidal constraint of the velocity field. Furthermore, an implicit discretization for the viscous terms is devised for the INS model, which is based on the VEM technique. Consequently, the CFL-type stability condition on the maximum admissible time step is based only on the fluid velocity and not on the celerity nor on the viscous eigenvalues. A large suite of test cases demonstrates the accuracy and the capabilities of the new family of schemes to solve relevant benchmarks in the field of incompressible fluids.

GEA: a new finite volume-based open source code for the numerical simulation of atmospheric and ocean flows

ECOGEN, an open-source tool dedicated to multiphase compressible multiphysics flows

A comparative computational study of different formulations of the compressible Euler equations for mesoscale atmospheric flows in a finite volume framework

Computational study of numerical flux schemes for mesoscale atmospheric flows in a Finite Volume framework

An Open-Source Code for High-Speed Non-Equilibrium Gas–Solid Flows in OpenFOAM

The GeoClaw software for depth-averaged flows with adaptive refinement

A new family of semi-implicit Finite Volume / Virtual Element methods for incompressible flows on unstructured meshes

OpenGeoSys: an open-source initiative for numerical simulation of thermo-hydro-mechanical/chemical (THM/C) processes in porous media

Efficient implementation of complex equations of state in a high-order framework

voFoam - A geometrical Volume of Fluid algorithm on arbitrary unstructured meshes with local dynamic adaptive mesh refinement using OpenFOAM

Numerical Wave Simulation Using Geometrical VOF Method Based on OpenFOAM

A finite element implementation of the incompressible Schrödinger flow method

An open-source finite volume toolbox for solid mechanics and fluid-solid interaction simulations

Ocellaris: a discontinuous Galerkin finite element solver for two-phase flows with high density differences

Numerical Simulations of Highly Nonlinear Steady and Unsteady Free Surface Flows

A potential energy conserving finite element method for turbulent variable density flow: application to glacier-fjord circulation

incompressibleFoam: a new time consistent framework with BDF and DIRK integration schemes

A finite element level set method based on adaptive octree meshes for thermal free‐surface flows

An open source massively parallel solver for Richards equation: Mechanistic modelling of water fluxes at the watershed scale

Enhanced Level-Set Method for free surface flow applications

An Open-Source Python-Based Boundary-Element Method Code for the Three-Dimensional, Zero-Froude, Infinite-Depth, Water-Wave Diffraction-Radiation Problem