Accuracy of operator splitting for advection‐dispersion‐reaction problems