<p>This study evaluates the capability of a Remeshed Vortex Method (RVM) to accurately simulate complex wall-bounded turbulent flows by performing both Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) of the flow over periodic hills. The RVM stems from Lagrangian Vortex methods, discretizing the flow fields on numerical particles, where a remeshing procedure is applied to redistribute the particles onto a Cartesian grid in order to control the interparticle distance. In DNS, different algorithmic configurations inherent to the RVM’s remeshing procedure are explored. The RVM shows good accuracy and computational efficiency in terms of mesh resolution and time steps compared to the DNS results reported in literature, despite the limitation of using a uniform Cartesian grid. In the LES framework, first a systematic analysis of the directional splitting sequence in the RVM algorithm for anisotropic flows allows us to select a priori the optimal direction ordering in the operators splitting procedure. Secondly, different subgrid-scale models are selected and assessed, including some models previously evaluated and calibrated in a HIT configuration. Among these, the Variational Multiscale (VMS) eddy-viscosity models and the Spectral Vanishing Viscosity (SVV) approaches, which introduce dissipation only at the smallest resolved vorticity scales, are found to be the most suitable for use with the RVM. In particular, the calibrated VMS Smagorinsky model shows the best performance, confirming our previous findings on the robustness of the model calibration.</p>