In this lecture “PDEs and applications”, Prof. Hiptmair mentioned that the blow up of Galerkin matrix condition number is related to the intrinsic instability of the boundary element basis in the energy norm. On the other hand, the discretization of the operator in the \(V\)-norm is perfectly stable.

In (Brenner, page 4), energy inner product is defined as the symmetric bilinear form \(a(\cdot ,\cdot )\). For a general second order PDE, the bilinear form is $$ \begin{equation} \label {eq:bilinear-form-2nd-order-pde} a(u,v):= \sum _{i,j=1}^d\int _{\Omega }a_{ji}(x)\frac {\pdiff u}{\pdiff x_{i}}\frac {\pdiff v}{\pdiff x_{j}}\intd x + \int _{\Omega }a_0uv \intd x. \end{equation} $$

The symmetry of the bilinear form is due to the self-adjointness of the differential operator \(L\) $$ \begin{equation} \label {eq:diff-operator-2nd-order} (Lu)(x):= -\sum _{i,j=1}^d \frac {\pdiff }{\pdiff x_j} \left [ a_{ji}(x)\frac {\pdiff u}{\pdiff x_{i}} \right ] + a_0(x)u(x). \end{equation} $$

Then in (Brenner, page 5), energy norm is defined via the bilinear form: $$ \begin{equation} \label {eq:energy-norm} \lVert v \rVert _E = \sqrt {a(v,v)} \quad \forall v\in X. \end{equation} $$ Meanwhile, \(a(\cdot ,\cdot )\) can also be considered as an inner product on \(X\). Cauchy-Schwartz inequality also holds for such inner product and energy norm: $$ \begin{equation} \lvert a(v,w) \rvert \leq \lVert v \rVert _E \lVert w \rVert _E. \end{equation} $$ On the other hand, the original norm and inner product on the space \(X\) can still be adopted, which has nothing to do with the differential operator \(L\).

Take the \(H^1(\Omega )\) Sobolev space as example, its standard norm is $$ \begin{equation} \| u \|_{H^1(\Omega )} = \left ( \| u \|_{L^2(\Omega )}^2 + \| \nabla u \|_{(L^2(\Omega ))^n}^2 \right )^{1/2}. \end{equation} $$ If the differential operator is \(L=-\Delta + I\), its corresponding bilinear form is \(a(u,v) = \left \langle \nabla u, \nabla v \right \rangle + \left \langle u,v \right \rangle \). Then the norm \(\lVert u \rVert _{H^1(\Omega )}\) is equivalent to \(\sqrt {a(u,u)}\), so it is an energy norm.

If the differential operator is \(L=-\Delta \), the bilinear form is \(a(u,v) = \left \langle \nabla u, \nabla v \right \rangle \). Then \(a(u,u) = \left \langle \nabla u, \nabla u \right \rangle \) and the energy norm \(\sqrt {a(u,u)}\) is actually the semi-norm \(\lvert u \rvert _{H^1(\Omega )}\).

In addition, because \(H^1(\Omega )\) is a subspace of \(L_2(\Omega )\), we can still assign the \(L_2\)-norm to \(H^1(\Omega )\). Of course, this norm is not the energy norm.

According to the definition of the standard \(H^1\)-norm, when the differential operator is \(L=-\Delta + I\), both \(L_2\)-norm and \(H^1\)-semi-norm is smaller than the energy norm. When the differential operator is \(L=-\Delta \), while the \(H^1\)-semi-norm is the same as the energy norm, whether \(L_2\)-norm is smaller than the energy norm seems unknown. However, in (Brenner, page 6), it is proved generally that the \(L_2\)-norm is much smaller than the energy norm.

N.B. The Ritz-Galerkin method ensures that the solution error is optimal in the energy norm. However, if we use \(L_2\)-norm to measure the error, since for a same error function, its value is smaller than the energy norm, the method is actually more stable. This explains the statement of Prof. Hiptmair at the beginning.

Finally, we can also see that the usually adopted function space norm and inner product on \(X\) only focus on the continuity and smoothness of the function itself. The energy norm and inner product defined via the bilinear form on \(X\) is related to the differential operator, or rather the physical property of the system.

References

   Susanne C. Brenner. The Mathematical Theory of Finite Element Methods. Springer New York, 3 edition. ISBN 978-1-4419-2611-1. URL https://book.douban.com/subject/4219448/#.

Backlinks: 《Spectral equivalence》