A four-way coupled Euler—Lagrange approach using a variational multiscale method for simulating cavitation