I want to avoid the use of numerical method in finding the roots of the inherent eigenvalues in the characteristics equation of the solution. This will help in increasing the accuracy of the solution and in providing a total analytical method.
Sure, you can use Fourier transforms of varying dimensions to solve multi-dimensional PDE's. The complexity depends much the same on lower dimensions, i.e. whether or not the space is simple and regular.