On the construction of suitable weak solutions to the 3D Navier–Stokes equations in a bounded domain by an artificial compressibility method