I'll give you an example that shows that that one non-separable solution of a PDE, could have been obtained a sum of separable solutions.
$\beta x \frac{\partial \psi}{\partial x} - y \frac{\partial \psi}{\partial y} = 0$
Using separation of variables you find a set of solutions parametrized by $\alpha$
$X_\alpha(x)Y_\alpha(y) = (xy^\beta)^{\alpha/2}$
However there exists an inseparable set of solutions parametrized by $\gamma$
$\text{exp}(\gamma xy^\beta)$
However in this particular example the inseparable solution can be cast as an infinite sum of separable ones (owing to the linearity of the PDE) and the specific form of the answers
$\text{exp}(\gamma xy^\beta) = \sum_{\alpha\in2\mathbb{N}}C_\alpha X_\alpha(x)Y_\alpha(y) =\sum_{\alpha\in2\mathbb{N}} \frac{\gamma^{\alpha/2}}{(\alpha/2)!}(xy^\beta)^{\alpha/2}$
so that $C_\alpha \equiv \frac{\gamma^{\alpha/2}}{(\alpha/2)!}$
This does not say in generality whether we can do the same thing in all situations. But as has been mentioned earlier, it does if the basis of solutions is a complete basis