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.

Implicit like time discretization for the one-phase Hele-Shaw problem with surface tension

Implicit like time discretization for the one-phase Hele-Shaw problem with surface tension

A Volume of Fluid Method Algorithm for Simulation of Surface Tension Dominant Two-Phase Flows

Removing the stiffness from interfacial flows with surface tension

Local Smoothing Effect and Existence for the One-Phase Hele–Shaw Problem with Zero Surface Tension

Diffuse Interface Model for Two-Phase Flows on Evolving Surfaces with Different Densities: Global Well-Posedness

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

A fully-coupled algorithm with implicit surface tension treatment for interfacial flows with large density ratios

Diffuse Interface Model for Two-Phase Flows on Evolving Surfaces with Different Densities: Local Well-Posedness

Two-phase flows with bulk-surface interaction: thermodynamically consistent Navier-Stokes-Cahn-Hilliard models with dynamic boundary conditions

Simplified conservative discretization of the Cahn-Hilliard-Navier-Stokes equations

Efficient solver of relativistic hydrodynamics with implicit Runge-Kutta method

Phase-field model for Hele-Shaw flows with arbitrary viscosity contrast. II. Numerical study

A stable and linear time discretization for a thermodynamically consistent model for two-phase incompressible flow

A cut-cell method for the numerical simulation of 3D multiphase flows with strong interfacial effects

Quinpi: Integrating stiff hyperbolic systems with implicit high order finite volume schemes

Implicit-Explicit Schemes for Compressible Cahn–Hilliard–Navier–Stokes Equations

Implicit and semi-implicit second-order time stepping methods for the Richards equation

On the two-phase Navier-Stokes equations with surface tension

An efficient numerical approach for singularly perturbed time delayed parabolic problems with two-parameters