“…It is well known that for problems with smooth solutions, the approximation obtained with high-order methods converges exponentially with the order of the approximating polynomial. More generally, high-order methods have been shown to deliver higher accuracy with a lower computational cost than low-order methods in many practical applications [6,7,8,9,10,11,12,13,14,15]. In addition, the accurate approximation of the domain geometry eliminates the spurious effects in the solution that can arise from a piecewise linear representation of the curved domain boundaries [16,17,18,19,20].…”