  We introduce a new class of Discontinuous Galerkin (DG) methods for solving
nonlinear conservation laws on unstructured Voronoi meshes that use a
nonconforming Virtual Element basis defined within each polygonal control
volume. The basis functions are evaluated as an L2 projection of the virtual
basis which remains unknown, along the lines of the Virtual Element Method
(VEM). Contrarily to the VEM approach, the new basis functions lead to a
nonconforming representation of the solution with discontinuous data across the
element boundaries, as typically employed in DG discretizations. To improve the
condition number of the resulting mass matrix, an orthogonalization of the full
basis is proposed. The discretization in time is carried out following the ADER
(Arbitrary order DERivative Riemann problem) methodology, which yields one-step
fully discrete schemes that make use of a coupled space-time representation of
the numerical solution. The space-time basis functions are constructed as a
tensor product of the virtual basis in space and a one-dimensional Lagrange
nodal basis in time. The resulting space-time stiffness matrix is stabilized by
an extension of the dof-dof stabilization technique adopted in the VEM
framework, hence allowing an element-local space-time Galerkin finite element
predictor to be evaluated. The novel methods are referred to as VEM-DG schemes,
and they are arbitrarily high order accurate in space and time. The new VEM-DG
algorithms are rigorously validated against a series of benchmarks in the
context of compressible Euler and Navier-Stokes equations. Numerical results
are verified with respect to literature reference solutions and compared in
terms of accuracy and computational efficiency to those obtained using a
standard modal DG scheme with Taylor basis functions. An analysis of the
condition number of the mass and space-time stiffness matrix is also forwarded.
