Higher‐order solutions of square driven cavity flow using a variable‐order multi‐grid method