Multi-component flow problems are typical of many technological and engineering applications. In this work, we propose an implicit high-order discontinuous Galerkin discretization of the variable density incompressible (VDI) flow model for the simulation of multi-component problems. Indeed, the peculiarity of the VDI model is that the density is treated as an advected property, which can be used to possibly track multiple (more than two) components. The interface between fluids is described by a smooth, but sharp, variation in the density field, thus not requiring any geometrical reconstruction. Godunov numerical fluxes, density positivity, mass conservation, and Gibbs-type phenomena at material interfaces are challenges that are considered during the numerical approach development. To avoid Courant-related time step restrictions, high-order single-step multi-stage implicit schemes are applied for the temporal integration. Several test cases with known analytical solutions are used to assess the current approach in terms of space, time, and mass conservation accuracy. As a challenging application, the simulation of a 2D droplet impinging on a thin liquid film is performed and shows the capabilities of the proposed DG approach when dealing with high-density (water–air) multi-component problems.