The answer to the general question you asked
What are the conditions[General Criteria] for the existence or non existence of the solutions to a PDE[Elliptic type] subject to given boundary conditions?
Is "there are none". To be a bit more verbose: the answer is that "it depends on what form of equation" you are allowed to consider and "what you admit as solutions".
Let me give one extremely general example. Consider the family of partial differential equations $F(\partial^2 u) = 0$, where $F$ is some smooth function defined on the space of $N\times N$ matrices. The equation is said to be elliptic if $\frac{\partial}{\partial X_{ij}} F(X_{ij})$ is elliptic (as a matrix, evaluated at all $X$). For these equations,
- In dimension $N = 2$, Nirenberg showed that the Dirichlet problem admits classical solutions.
- In arbitrary dimensions, viscosity solutions exist (using the Perron method).
- If $F$ is convex, viscosity solutions are in fact classical (Krylov-Evans regularity).
- If $F$ is not convex, there can be cases where a viscosity solution exists but is not classical, which by the uniqueness properties of viscosity solutions, imply that such equations do not admit classical solutions.
If you restrict to only quasilinear instead of fully nonlinear equations, there are a wealth of existence results for the Dirichlet problem. In fact, in section 15.5 of their book, Gilbarg and Trudinger wrote that
It is not feasible to present here a comprehensive account of existence theorems for the classical Dirichlet problem that follow by combination of the results of Chapters 10, 13, 14, and 15.
... and they proceed to list 10 different existence theorems applicable in different contexts. What's important to note, however, is that there is a fairly strong non-existence theorem that applies even in the case of the semilinear problem. I will state here a special case
Theorem (c.f. Theorem 14.11 in Gilbarg-Trudinger) Let $\Omega\subset\mathbb{R}^n$ be a bounded domain such that the largest ball contained in $\Omega$ has radius $R$, and consider the Dirichlet problem for $ \triangle u(x) + b(x,u, \nabla u) = 0~. $ If we assume that there is a positive constant $\theta$ such that for every $p\in \mathbb{R}^n$, $z \in \mathbb{R}$, and $x\in\mathbb{R}^n$ you have that $ |b(x,z,p)| \geq \frac{n |p|}{R} + |p|^{2+\theta} $ then one can find a smooth boundary data $\phi$ for which the Dirichlet problem does not admit a solution.
Note, however, that this theorem does not apply to your special case, since your $f(u)$ does not depend on the first derivative.
Now, to address your specific problem: you outlined an approach that requires expanding $u$ in terms of the eigenfunctions of the Laplacian. This is roughly equivalent to taking the Fourier transform of the problem, and so you see that to directly apply the solutions to the reduced wave equation you need that $f(u)$ to be linear ($f(au+bv) = af(u) + bf(v)$), which puts strong restrictions on $f$.
Now, for a semilinear problem, one can always convert a problem with nonzero boundary condition to one with zero boundary conditions: say you want to solve
$ \triangle u = f(x,u) $
with boundary condition $\phi$. First solve the homogeneous problem $\triangle w = 0$ with boundary condition $\phi$. Then write the ansatz that $u = v + w$ where $v$ vanishes on the boundary. Then $v$ solves
$ \triangle v = \tilde{f}(x,v) := f(x,v+w) $
In general, however, for a problem like the one you proposed $\triangle u = f(u)$, you are often better off demonstrating existence of a solution using a variational method (if possible) or a fixed-point/iteration method. (If you are really lucky you can use Perron method, but for general nonlinearities $f(u)$ the comparison results for super and sub solutions may not be available.)