Understanding about energy norm used in Galerkin method
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(⋅,⋅). For a general second order PDE, the bilinear form is a(u,v):=d∑i,j=1∫Ωaji(x)∂u∂xi∂v∂xjdx+∫Ωa0uvdx.
The symmetry of the bilinear form is due to the self-adjointness of the differential operator L (Lu)(x):=−d∑i,j=1∂∂xj[aji(x)∂u∂xi]+a0(x)u(x).
Then in (Brenner, page 5), energy norm is defined via the bilinear form: ‖v‖E=√a(v,v)∀v∈X. Meanwhile, a(⋅,⋅) can also be considered as an inner product on X. Cauchy-Schwartz inequality also holds for such inner product and energy norm: |a(v,w)|≤‖v‖E‖w‖E. 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 H1(Ω) Sobolev space as example, its standard norm is ‖u‖H1(Ω)=(‖u‖2L2(Ω)+‖∇u‖2(L2(Ω))n)1/2. If the differential operator is L=−Δ+I, its corresponding bilinear form is a(u,v)=⟨∇u,∇v⟩+⟨u,v⟩. Then the norm ‖u‖H1(Ω) is equivalent to √a(u,u), so it is an energy norm.
If the differential operator is L=−Δ, the bilinear form is a(u,v)=⟨∇u,∇v⟩. Then a(u,u)=⟨∇u,∇u⟩ and the energy norm √a(u,u) is actually the semi-norm |u|H1(Ω).
In addition, because H1(Ω) is a subspace of L2(Ω), we can still assign the L2-norm to H1(Ω). Of course, this norm is not the energy norm.
According to the definition of the standard H1-norm, when the differential operator is L=−Δ+I, both L2-norm and H1-semi-norm is smaller than the energy norm. When the differential operator is L=−Δ, while the H1-semi-norm is the same as the energy norm, whether L2-norm is smaller than the energy norm seems unknown. However, in (Brenner, page 6), it is proved generally that the L2-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 L2-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》