3
$\begingroup$

I have a question regarding the weak formulation of the poisson equaiton:

$-\nabla^2 u = f$ on $\Omega \subset \mathbb{R}^n$, $u_{\delta\omega}=0$

In my notes as well as on wiki, it says that the weak formulation can be derived using integration by parts and greens identity. Can someone do this explicit? Where does the use of greens identity come in?

Regards

$\endgroup$

1 Answer 1

9
$\begingroup$

Take a test function $v$ such that $v=0$ on $\partial \Omega$, multiply both sides by this $v$ and integrate over $\Omega$, we get (I implicitly suppose that all fields and functions are sufficiently regular to apply those operators involved):

$$-\int_{\Omega} (\Delta u) v = \int_{\Omega} fv $$

Now, remember the identity $ \text{div} (v \nabla u )= (\nabla u )(\nabla v) + v \Delta u $ to get:

$$\int_{\Omega} \nabla u \nabla v - \int_{\Omega} \text{div} (v \nabla u ) = \int_{\Omega} fv $$

By Green's theorem we can write (if $\partial \Omega$ is sufficiently regular):

$$ \int_{\Omega} \nabla u \nabla v - \int_{\partial \Omega} (v \nabla u) \cdot \mathbf{n} \, dS = \int_{\Omega} fv $$

Finally, because $v=0$ on $\partial \Omega$ we conclude:

$$ \int_{\Omega} \nabla u \nabla v = \int_{\Omega} fv $$

Note that the identity $\text{div} (f \mathbf{F})= \nabla f \cdot \mathbf{F} + f \text{div} \mathbf{F}$ plays the role of "integration by parts", as Wikipedia states.

$\endgroup$
7
  • $\begingroup$ Thanks! Yeah I have one further question wrt. integration by parts. If you just do it clasically isn't that enough then. You get the surface integral as well right? $\endgroup$
    – 1233023
    Commented Feb 21, 2017 at 16:28
  • $\begingroup$ Edited it:) I think this is good as a general approach, if you want something more similar to the integration by parts in one dimension, we should work in a Sobolev space and test on compactly supported functions. $\endgroup$
    – GaC
    Commented Feb 21, 2017 at 16:31
  • $\begingroup$ Well, in one dimension the gradient of a scalar function is just its derivative and the Laplacian is just the second derivative! $\endgroup$
    – GaC
    Commented Feb 21, 2017 at 16:32
  • $\begingroup$ I can add the example in one dimension if you need it ! $\endgroup$
    – GaC
    Commented Feb 21, 2017 at 16:33
  • $\begingroup$ The problem with the one dimensional example, is that then classical integration by parts just does it, you don't get the insight with the boundary integral, since it's just two points you have to substract. $\endgroup$
    – 1233023
    Commented Feb 21, 2017 at 16:34

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .