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

被引:3
|
作者
Boscheri, Walter [1 ]
Chiozzi, Andrea [2 ]
Carlino, Michele Giuliano [1 ]
Bertaglia, Giulia [2 ]
机构
[1] Univ Ferrara, Dept Math & Comp Sci, Via Machiavelli 30, I-44121 Ferrara, Italy
[2] Univ Ferrara, Dept Environm & Prevent Sci, Cso Ercole 1dEste 32, I-44121 Ferrara, Italy
关键词
Semi-implicit schemes; Finite volume; Virtual Element Method; High order in space and time; Asymptotic preserving; Incompressible; flows; NAVIER-STOKES EQUATIONS; DISCONTINUOUS GALERKIN METHOD; SHALLOW-WATER EQUATIONS; CENTRAL WENO SCHEMES; RUNGE-KUTTA SCHEMES; FREE-SURFACE FLOWS; VELOCITY RECONSTRUCTION; HYPERBOLIC SYSTEMS; DIFFERENCE METHODS; VIRTUAL ELEMENTS;
D O I
10.1016/j.cma.2023.116140
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
We introduce a new family of high order accurate semi-implicit schemes for the solution of nonlinear time-dependent systems of partial differential equations (PDE) 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, hence 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 then achieved using the semi-implicit IMEX Runge-Kutta schemes, and the novel schemes are proven to be asymptotic preserving (AP) and well-balanced (WB). 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 also 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.(c) 2023 Elsevier B.V. All rights reserved.
引用
收藏
页数:39
相关论文
共 50 条