  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.
