A novel hybrid grid-based and grid-free computational method which is called VLES is proposed and validated over several benchmark cases. VLES splits the flow structures into large scale ones, resolved on the grid (Eulerian approach), and small scale ones, represented by vortex particles (Lagrangian approach). Two transport equations for grid and particle solution are derived which are dynamically coupled through the existence of coupling terms in each of them. The method resembles LES with an effort to directly reproduce the subgrid motion at least in the statistical sense.<eng>