TY - JOUR
T1 - The Bramble-Hilbert lemma for convex domains
AU - Dekel, S.
AU - Leviatan, D.
PY - 2004
Y1 - 2004
N2 - The Bramble-Hilbert lemma is a fundamental result on multivariate polynomial approximation. It is frequently applied in the analysis of finite elements methods (FEM) used for numerical solutions of PDEs. However, this classical estimate depends on the geometry of the domain and may "blow up" for simple examples such as a sequence of triangles of equivalent diameter that become thinner and thinner. Thus, in FEM applications one usually requires that the mesh has "quasi-uniform" geometry. This assumption is perhaps too restrictive when one tries to obtain estimates of nonlinear approximation methods that use piecewise polynomials. Our main result that improves upon this point is the following. Let Ω ⊂ ℝn be a bounded convex domain and let g ∈ Wpm(Ω), m ∈ ℕ, 1 ≤ p ≤ ∞, where Wpm(Ω) is the Sobolev space. Then there exists a polynomial P of total degree m - 1 for which |g - P|k,P ≤ C(n, m)(diam Ω)m-k|g| m,p, k = 0, 1,..., m, where|·|k,p := Σ |α|=k ∥Dα ∥Lp(Ω) is the Sobolev seminorm of order k. As a consequence we get that for f ∈ Lp(Ω), Em-1(f, Ω)p ≈ K m(f, (diam Ω)m)p, where Em-1 (f, Ω)p := infp∈m-1∥ f-P∥ Lp(Ω) is the error of polynomial approximation of degree m - 1 and Km( , )p is the K-functional associated with the pair (Lp(Ω), Wpm(Ω)), and where the constants of equivalence depend only on m and n. For the case of convex domains (elements) this extends a recent result for p = 2, and for m = 1 and 2 < p ≤ ∞. This also improves previous results where the constant in the estimate further depends on the geometry of the domain, or where there is a constraint p > n(≥ 2).
AB - The Bramble-Hilbert lemma is a fundamental result on multivariate polynomial approximation. It is frequently applied in the analysis of finite elements methods (FEM) used for numerical solutions of PDEs. However, this classical estimate depends on the geometry of the domain and may "blow up" for simple examples such as a sequence of triangles of equivalent diameter that become thinner and thinner. Thus, in FEM applications one usually requires that the mesh has "quasi-uniform" geometry. This assumption is perhaps too restrictive when one tries to obtain estimates of nonlinear approximation methods that use piecewise polynomials. Our main result that improves upon this point is the following. Let Ω ⊂ ℝn be a bounded convex domain and let g ∈ Wpm(Ω), m ∈ ℕ, 1 ≤ p ≤ ∞, where Wpm(Ω) is the Sobolev space. Then there exists a polynomial P of total degree m - 1 for which |g - P|k,P ≤ C(n, m)(diam Ω)m-k|g| m,p, k = 0, 1,..., m, where|·|k,p := Σ |α|=k ∥Dα ∥Lp(Ω) is the Sobolev seminorm of order k. As a consequence we get that for f ∈ Lp(Ω), Em-1(f, Ω)p ≈ K m(f, (diam Ω)m)p, where Em-1 (f, Ω)p := infp∈m-1∥ f-P∥ Lp(Ω) is the error of polynomial approximation of degree m - 1 and Km( , )p is the K-functional associated with the pair (Lp(Ω), Wpm(Ω)), and where the constants of equivalence depend only on m and n. For the case of convex domains (elements) this extends a recent result for p = 2, and for m = 1 and 2 < p ≤ ∞. This also improves previous results where the constant in the estimate further depends on the geometry of the domain, or where there is a constraint p > n(≥ 2).
KW - Bramble-Hilbert lemma
KW - Finite element methods
KW - Multivariate nonlinear approximation
UR - http://www.scopus.com/inward/record.url?scp=4644280544&partnerID=8YFLogxK
U2 - 10.1137/S0036141002417589
DO - 10.1137/S0036141002417589
M3 - ???researchoutput.researchoutputtypes.contributiontojournal.article???
AN - SCOPUS:4644280544
SN - 0036-1410
VL - 35
SP - 1203
EP - 1212
JO - SIAM Journal on Mathematical Analysis
JF - SIAM Journal on Mathematical Analysis
IS - 5
ER -