In classical PDE courses it is common to learn to perform a change of variables without really learning how to find the adequate equations of the change (polar, cylindrical or spherical coordinates are just plain easy to detect).
Now, is it always possible to find a change of variables that transforms a PDE into a product of independent functions $f_i(x_i)$, with $$\frac{\partial f_i}{\partial x_j} = \delta_{ij}$$ How is this change of variables found? And, in case it doesn't always exist, what are the symmetries the expression of the PDE must have so that we can prove its existence?
Bonus question: what relationship does it have with the Hamilton-Jacobi equations?