One can solve Poisson’s problem $-\Delta u = f$ in $d$ dimensions with homogeneous Dirichlet boundary conditions using a mixed formulation as explained below:

Let $\sigma = \nabla u$, then for a sufficiently smooth function $\tau$, by Green’s theorem \begin{align*} (\sigma, \tau) &= (\nabla u, \tau) &= -(u, \textrm{div } \tau). \end{align*} Again, choosing $v$ a function sufficiently smooth, we have \begin{align*} f = -\textrm{div } \sigma \implies (f, v) = (-\textrm{div } \sigma, v). \end{align*} This gives the saddle-point problem: find $(u, \sigma) \in V \times M$ such that \begin{align*} (\sigma, \tau) + (u, \textrm{div } \tau) &= 0 (\textrm{div } \sigma, v) &= -(f, v) \end{align*} hold for all $(\tau, v) \in V \times M$. Note that we don’t have to take a derivative of $u$, hence it’s natural to try $M = L^2$, but what about the space $V$?

One very easy choice to guess is $V = [H^1(\Omega)]^d$ as we want the divergence to be all defined, but unfortunately this doesn’t work as the gradient of the solution to Poisson’s problem can easily not be in $[H^1(\Omega)]^d$.

In order to illustrate this, consider $u =\left(r^{2/3}-r^{5/3}\right) \sin \left(\frac{2 \theta }{3}\right)$ on the domain of the unit circle with bottom left quarter taken out. It’s not hard to see that $u = 0$ on the boundary of the domain, and we can easily find the $f$ such that it satisfies Poisson’s equation. Now, we can either calculate the gradient exactly or argue as follows.

First, recall how to take a gradient in polar coordinates. Note that $tial_r u \approx r^{-1/3}$ plus higher order terms and that $\frac{1}{r}tial_\theta u \approx r^{-1/3}$ plus higher order terms also. Now, one can easily calculate the $H^1$ seminorm to see that the derivative is unbounded as we’re integrating over $[0,1]$ with $(r^{-4/3})^2r$ terms (the extra $r$ comes from the change of basis from polar integration).

The above is an example of why the space $H(\textrm{div})$ is needed.