Las soluciones débiles de la ecuación son puntos críticos del funcional energético $$ E(u)=\int_\Omega\Bigl(\frac12|\nabla u|^2+\frac1{3}|u|^3-f\,u\bigr)dx. $$ Observe que $E(u)$ está bien definida en el espacio de Banach $X=H^1_0(\Omega)\cap L^3(\Omega)$ .
Para demostrar la existencia de una solución basta con demostrar que $E$ alcanza su mínimo en $X$ . Esto se hace demostrando que $E$ es estrictamente convexo y coercitivo (es decir $E(u)\to\infty$ si $\|u\|_X\to\infty$ ).
La unicidad se deduce (por un método diferente) porque la función $g(u)=u\,|u|$ está aumentando.
Hay muchos libros sobre métodos variacionales. Uno de ellos es Introducción a la teoría del punto crítico por O. Kavian (en francés.) Un caso más general de su pregunta se resuelve en la página 139.
Lo que sigue es una explicación formal; no entraré en los detalles de los espacios correctos,...
La formulación débil del problema es $$ \int_\Omega\bigl(\nabla u\cdot\nabla u+u\,|u|\,v-f\,v\bigr)dx=0\quad v\text{ a test function.} $$ El funcional de energía es un funcional $E:X\to\mathbb{R}$ cuyo derivada con respecto a $u$ viene dada precisamente por la formulación débil. ¿Qué significa esto? El diferencial en $u$ se define como un elemento $E'(u)$ del dual de $X$ tal que $$ E(u+v)=E(u)+\langle E'(u),v\rangle+o(\|v\|). $$ Una forma de calcularlo es $$ \langle E'(u),v\rangle=\lim_{t\to0}\frac{E(u+t\,v)-E(v)}{t}. $$ Utilizando esto es fácil demostrar que $E(u)$ es de hecho el funcional de energía asociado a la ecuación. Obsérvese que $|u|^3/3$ es una primitiva de $u\,|u|$ .