Hi,
I solved for a Poisson equation with finite elements, using piecewise linear basis functions on 2d triangles.  

Now, I want to evaluate the following expressions:  

$$ \int\_\Omega \Delta u ~dx$$
and 
$$ \int\_\Omega (\Delta u)^2 ~dx$$
I want to evaluate these expressions using my approximated solution $u$ which has been computed on piecewise linear basis functions.


For the first one, I thought of using the identity $\int\_\Omega div \nabla u = \int\_{\partial\Omega} \nabla u . \vec{n}~ds$ and summing this expression over each triangle.However, as expected, the result is strictly 0 since the basis functions are linear.  

I also tried to use a kind of jump formula (like $f'(x)=\tilde{f}'(x) + f^+-f^-$ where $\tilde{f}'$ is the derivative of the smooth part of f) but I'm stuck on how to do that for each triangle in 2D (the outer normal is likely to cancel out when computing the same formula for two adjacent triangles sharing an edge) - and I'm wondering if it is supposed to work.


For the second one, I just have no clue.


Am I forced to use higher order elements ? Any idea ? Thanks!



integration of a laplacian
Since your approximate solution is piecewise linear, it is $H^1$ but not $H^2$. Therefore your calculation is impossible. You can do to things to overcome the difficulty:


* Either use higher-order elements,
* or post-process your approximate solution $u$. This means constructing a smoother $\bar u$ using some convolution by an appropriate $\phi(x/h)$, where $h$ is the typical mesh size. In one space dimension, this amounts to using splines. Then $\bar u\in H^2$ and your calculation is meaningful.


