Hi, let $\Omega \subset R^3$ be a bounded smooth domain, consider the elliptic equation
\begin{equation}
-\triangle u + u^2div u = f, \quad x \in \Omega, \quad u\big|\_{\partial \Omega} = 0.
\end{equation}
Here the force $f \in L^2$, and $div u = \sum\_j \partial u/\partial x\_j $. Is there any solutions of the equation in $H^1$? And how to prove it. Thanks for any hints and references.


\textbf{Edit}: In fact, the $1-d$ equation we considered is the stationary equation of the BBM equation 
$$
u\_t - u\_{txx} - u\_{xx} + u^2u\_x = f.
$$ 



Solvability of a nonlinear elliptic equation
Firstly, turning to the $\mathbb{R}^1$ case, we note that
$$
 -u\_{xx}+u^2u\_x=-\partial\_x\left(u\_x-\frac{1}{3}u^3\right)=f(x)
$$
and so we get a Chini equation (see [Inhomogeneous Bernoulli Equation](https://mathoverflow.net/questions/104006/inhomogeneous-bernoulli-equation)):
$$
 -u\_x+\frac{1}{3}u^3=\int\_{x\_0}^xf(x')dx'+C
$$
and a general solution to this equation is not known. Turning to $\mathbb{R}^3$, you can see that your equation can be written down in the form
$$ 
 \nabla\cdot {\bf v}=-f
$$
being
$$
 {\bf v}=\left(u\_x-\frac{1}{3}u^3,u\_y-\frac{1}{3}u^3,u\_z-\frac{1}{3}u^3\right)
$$
and you have to satisfy the Gauss identity
$$
 \int\_{\partial\Omega}{\bf v}\cdot d{\bf S}=-\int\_{\Omega}dV f.
$$
This means that one has to find three functions $v\_1,\ v\_2,\ v\_3$ that satisfy three different Chini equations. Again, a general solution misses in this case. The conclusion to draw is that you need to specialize the problem to see if this becomes manageable in some case or resort to some perturbation technique.


