A massively-parallel GPU-accelerated high-order solver for direct numerical simulations of non-ideal wall-bounded flows