Simulating thermal density operators with cluster expansions and tensor networks