Convergence rates to stationary solutions of a gas–liquid model with external forces