The complexity of the problem in three dimensional domains can only be captured in higher-order polynomial equations, and FEM is not efficient to integrate higher-order polynomial equations. Spectral-element method (SEM) has excellent error properties (exponential convergence) being the fastest possible: 1) geometrical flexibility, 2) high computational efficiency and 3) reliable spectral accuracy. SEM combines the flexibility of FEM with the accuracy...