Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element...
Transcript of Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element...
![Page 1: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/1.jpg)
Abstract Variational Problem
An abstract boundary value problem can be written in the form
Lu = f inD, Bu = 0 on ∂D
with a differential operator L and a boundary operator B.
Incorporating the boundary conditions in a Hilbert space H, the differentialequation usually admits a variational formulation
a(u, v) = λ(v), v ∈ H
with a bilinear form a and a linear functional λ.
This weak form of the boundary value problem is well suited for numericalapproximations, in particular because it requires less regularity. For adifferential operator of order 2m, the existence of weak derivatives up toorder m suffices.
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 1
![Page 2: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/2.jpg)
Abstract Variational Problem
An abstract boundary value problem can be written in the form
Lu = f inD, Bu = 0 on ∂D
with a differential operator L and a boundary operator B.Incorporating the boundary conditions in a Hilbert space H, the differentialequation usually admits a variational formulation
a(u, v) = λ(v), v ∈ H
with a bilinear form a and a linear functional λ.
This weak form of the boundary value problem is well suited for numericalapproximations, in particular because it requires less regularity. For adifferential operator of order 2m, the existence of weak derivatives up toorder m suffices.
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 1
![Page 3: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/3.jpg)
Abstract Variational Problem
An abstract boundary value problem can be written in the form
Lu = f inD, Bu = 0 on ∂D
with a differential operator L and a boundary operator B.Incorporating the boundary conditions in a Hilbert space H, the differentialequation usually admits a variational formulation
a(u, v) = λ(v), v ∈ H
with a bilinear form a and a linear functional λ.
This weak form of the boundary value problem is well suited for numericalapproximations, in particular because it requires less regularity. For adifferential operator of order 2m, the existence of weak derivatives up toorder m suffices.
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 1
![Page 4: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/4.jpg)
Ritz–Galerkin Approximation
The Ritz–Galerkin approximation uh =∑
i uiBi ∈ Bh ⊂ H of thevariational problem
a(u, v) = λ(v), v ∈ H,
is determined by the linear system∑i
a(Bi ,Bk) ui = λ(Bk),
which we abbreviate as GU = F .
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 2
![Page 5: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/5.jpg)
Example
model problem−∆u = f inD, u = 0 on ∂D
differential and boundary operators
L = −∆, Bu = u
bilinear form and linear functional
a(u, v) =
∫D
grad u grad v , λ(v) =
∫Dfv
Hilbert space: H = H10 (D)
simple finite element subspace Bh:piecewise linear functions on a triangulation of D
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 3
![Page 6: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/6.jpg)
Example
model problem−∆u = f inD, u = 0 on ∂D
differential and boundary operators
L = −∆, Bu = u
bilinear form and linear functional
a(u, v) =
∫D
grad u grad v , λ(v) =
∫Dfv
Hilbert space: H = H10 (D)
simple finite element subspace Bh:piecewise linear functions on a triangulation of D
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 3
![Page 7: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/7.jpg)
Example
model problem−∆u = f inD, u = 0 on ∂D
differential and boundary operators
L = −∆, Bu = u
bilinear form and linear functional
a(u, v) =
∫D
grad u grad v , λ(v) =
∫Dfv
Hilbert space: H = H10 (D)
simple finite element subspace Bh:piecewise linear functions on a triangulation of D
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 3
![Page 8: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/8.jpg)
Example
model problem−∆u = f inD, u = 0 on ∂D
differential and boundary operators
L = −∆, Bu = u
bilinear form and linear functional
a(u, v) =
∫D
grad u grad v , λ(v) =
∫Dfv
Hilbert space: H = H10 (D)
simple finite element subspace Bh:piecewise linear functions on a triangulation of D
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 3
![Page 9: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/9.jpg)
Example
model problem−∆u = f inD, u = 0 on ∂D
differential and boundary operators
L = −∆, Bu = u
bilinear form and linear functional
a(u, v) =
∫D
grad u grad v , λ(v) =
∫Dfv
Hilbert space: H = H10 (D)
simple finite element subspace Bh:piecewise linear functions on a triangulation of D
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 3
![Page 10: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/10.jpg)
Ellipticity
A bilinear form a on a Hilbert space H is elliptic if it is bounded andequivalent to the norm on H, i.e., if for all u, v ∈ H
|a(u, v)| ≤ cb‖u‖‖v‖, ce‖u‖2 ≤ a(u, u)
with positive constants cb and ce .
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 4
![Page 11: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/11.jpg)
Example
Poisson bilinear form
a(u, v) =
∫D
grad u grad v
Cauchy–Schwarz inequality =⇒
|a(u, v)| ≤ a(u, u)1/2 a(v , v)1/2 =
(∫D‖grad u‖2
)1/2(∫D‖grad v‖2
)1/2
≤ ‖u‖1‖v‖1
where
‖w‖1 =
(∫D|w |2 + ‖gradw‖2
)1/2
is the norm on H = H10 (D) ⊂ H1(D)
cb = 1
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 5
![Page 12: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/12.jpg)
Example
Poisson bilinear form
a(u, v) =
∫D
grad u grad v
Cauchy–Schwarz inequality =⇒
|a(u, v)| ≤ a(u, u)1/2 a(v , v)1/2 =
(∫D‖grad u‖2
)1/2(∫D‖grad v‖2
)1/2
≤ ‖u‖1‖v‖1
where
‖w‖1 =
(∫D|w |2 + ‖gradw‖2
)1/2
is the norm on H = H10 (D) ⊂ H1(D)
cb = 1
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 5
![Page 13: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/13.jpg)
Example
Poisson bilinear form
a(u, v) =
∫D
grad u grad v
Cauchy–Schwarz inequality =⇒
|a(u, v)| ≤ a(u, u)1/2 a(v , v)1/2 =
(∫D‖grad u‖2
)1/2(∫D‖grad v‖2
)1/2
≤ ‖u‖1‖v‖1
where
‖w‖1 =
(∫D|w |2 + ‖gradw‖2
)1/2
is the norm on H = H10 (D) ⊂ H1(D)
cb = 1
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 5
![Page 14: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/14.jpg)
Poincare–Friedrichs inequality =⇒∫D|u|2 ≤ const(D)
∫D‖grad u‖2, u ∈ H1
0 (D)
add∫D ‖grad u‖2 = a(u, u) to both sides
‖u‖21 ≤ (const(D) + 1)
∫D‖grad u‖2 ,
i.e., ce = (const(D) + 1)−1
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 6
![Page 15: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/15.jpg)
Poincare–Friedrichs inequality =⇒∫D|u|2 ≤ const(D)
∫D‖grad u‖2, u ∈ H1
0 (D)
add∫D ‖grad u‖2 = a(u, u) to both sides
‖u‖21 ≤ (const(D) + 1)
∫D‖grad u‖2 ,
i.e., ce = (const(D) + 1)−1
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 6
![Page 16: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/16.jpg)
Lax–Milgram Theorem
If a is an elliptic bilinear form and λ is a bounded linear functional on aHilbert space H, then the variational problem
a(u, v) = λ(v), v ∈ V ,
has a unique solution u ∈ V for any closed subspace V of H. Moreover, ifa is symmetric, the solution u can be characterized as the minimum of thequadratic form
Q(u) =1
2a(u, u)− λ(u)
on V .
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 7
![Page 17: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/17.jpg)
Special Case
finite element approximation: V = Bh
variational problem ⇔ Ritz–Galerkin system GU = Fellipticity of a =⇒ positive definiteness of G :
UGU =∑i ,k
uka(Bi ,Bk)ui = a(uh, uh) ≥ ce‖uh‖2 > 0
for uh 6= 0existence of G−1 =⇒ unique solvability of the Ritz–Galerkin system part of the Lax–Milgram theorem, relevant for numerical schemes
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 8
![Page 18: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/18.jpg)
Special Case
finite element approximation: V = Bh
variational problem ⇔ Ritz–Galerkin system GU = F
ellipticity of a =⇒ positive definiteness of G :
UGU =∑i ,k
uka(Bi ,Bk)ui = a(uh, uh) ≥ ce‖uh‖2 > 0
for uh 6= 0existence of G−1 =⇒ unique solvability of the Ritz–Galerkin system part of the Lax–Milgram theorem, relevant for numerical schemes
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 8
![Page 19: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/19.jpg)
Special Case
finite element approximation: V = Bh
variational problem ⇔ Ritz–Galerkin system GU = Fellipticity of a =⇒ positive definiteness of G :
UGU =∑i ,k
uka(Bi ,Bk)ui = a(uh, uh) ≥ ce‖uh‖2 > 0
for uh 6= 0
existence of G−1 =⇒ unique solvability of the Ritz–Galerkin system part of the Lax–Milgram theorem, relevant for numerical schemes
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 8
![Page 20: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/20.jpg)
Special Case
finite element approximation: V = Bh
variational problem ⇔ Ritz–Galerkin system GU = Fellipticity of a =⇒ positive definiteness of G :
UGU =∑i ,k
uka(Bi ,Bk)ui = a(uh, uh) ≥ ce‖uh‖2 > 0
for uh 6= 0existence of G−1 =⇒ unique solvability of the Ritz–Galerkin system
part of the Lax–Milgram theorem, relevant for numerical schemes
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 8
![Page 21: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/21.jpg)
Special Case
finite element approximation: V = Bh
variational problem ⇔ Ritz–Galerkin system GU = Fellipticity of a =⇒ positive definiteness of G :
UGU =∑i ,k
uka(Bi ,Bk)ui = a(uh, uh) ≥ ce‖uh‖2 > 0
for uh 6= 0existence of G−1 =⇒ unique solvability of the Ritz–Galerkin system part of the Lax–Milgram theorem, relevant for numerical schemes
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 8
![Page 22: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/22.jpg)
Proof
variational problem ⇔ identity between bounded linear functionals on V
left sideAu : v 7→ a(u, v)
bounded in view of the boundedness of a:
‖Au‖ = sup‖v‖=1
|(Au)(v)| = sup |a(u, v)| ≤ cb‖u‖
Riesz theorem representation for bounded linear functionals % : V → R:
%(v) = 〈R%, v〉
with 〈·, ·〉 the scalar product on V (identical to the scalar product on H)and R an isometry onto V
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 9
![Page 23: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/23.jpg)
Proof
variational problem ⇔ identity between bounded linear functionals on Vleft side
Au : v 7→ a(u, v)
bounded in view of the boundedness of a:
‖Au‖ = sup‖v‖=1
|(Au)(v)| = sup |a(u, v)| ≤ cb‖u‖
Riesz theorem representation for bounded linear functionals % : V → R:
%(v) = 〈R%, v〉
with 〈·, ·〉 the scalar product on V (identical to the scalar product on H)and R an isometry onto V
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 9
![Page 24: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/24.jpg)
Proof
variational problem ⇔ identity between bounded linear functionals on Vleft side
Au : v 7→ a(u, v)
bounded in view of the boundedness of a:
‖Au‖ = sup‖v‖=1
|(Au)(v)| = sup |a(u, v)| ≤ cb‖u‖
Riesz theorem representation for bounded linear functionals % : V → R:
%(v) = 〈R%, v〉
with 〈·, ·〉 the scalar product on V (identical to the scalar product on H)and R an isometry onto V
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 9
![Page 25: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/25.jpg)
Proof
variational problem ⇔ identity between bounded linear functionals on Vleft side
Au : v 7→ a(u, v)
bounded in view of the boundedness of a:
‖Au‖ = sup‖v‖=1
|(Au)(v)| = sup |a(u, v)| ≤ cb‖u‖
Riesz theorem representation for bounded linear functionals % : V → R:
%(v) = 〈R%, v〉
with 〈·, ·〉 the scalar product on V (identical to the scalar product on H)and R an isometry onto V
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 9
![Page 26: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/26.jpg)
equivalent form of the variational problem
RAu = Rλ
rewrite as fixed point equation
u = u − ωRAu︸ ︷︷ ︸Su
+ωRλ
with ω > 0show: S is contraction for small ω( Banach’s fixed point theorem completes proof of first part)estimate
‖S‖ = sup‖u‖=1
〈u − ωRAu, u − ωRAu〉1/2
using the ellipticity of afor ω = ce/c
2b
〈. . . , . . .〉 = ‖u‖2 − 2ωa(u, u) + ω2‖RAu‖2 ≤ 1− 2ωce + ω2c2b < 1
since 〈RAu, u〉 = (Au)(u) = a(u, u)
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 10
![Page 27: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/27.jpg)
equivalent form of the variational problem
RAu = Rλ
rewrite as fixed point equation
u = u − ωRAu︸ ︷︷ ︸Su
+ωRλ
with ω > 0
show: S is contraction for small ω( Banach’s fixed point theorem completes proof of first part)estimate
‖S‖ = sup‖u‖=1
〈u − ωRAu, u − ωRAu〉1/2
using the ellipticity of afor ω = ce/c
2b
〈. . . , . . .〉 = ‖u‖2 − 2ωa(u, u) + ω2‖RAu‖2 ≤ 1− 2ωce + ω2c2b < 1
since 〈RAu, u〉 = (Au)(u) = a(u, u)
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 10
![Page 28: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/28.jpg)
equivalent form of the variational problem
RAu = Rλ
rewrite as fixed point equation
u = u − ωRAu︸ ︷︷ ︸Su
+ωRλ
with ω > 0show: S is contraction for small ω( Banach’s fixed point theorem completes proof of first part)
estimate‖S‖ = sup
‖u‖=1〈u − ωRAu, u − ωRAu〉1/2
using the ellipticity of afor ω = ce/c
2b
〈. . . , . . .〉 = ‖u‖2 − 2ωa(u, u) + ω2‖RAu‖2 ≤ 1− 2ωce + ω2c2b < 1
since 〈RAu, u〉 = (Au)(u) = a(u, u)
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 10
![Page 29: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/29.jpg)
equivalent form of the variational problem
RAu = Rλ
rewrite as fixed point equation
u = u − ωRAu︸ ︷︷ ︸Su
+ωRλ
with ω > 0show: S is contraction for small ω( Banach’s fixed point theorem completes proof of first part)estimate
‖S‖ = sup‖u‖=1
〈u − ωRAu, u − ωRAu〉1/2
using the ellipticity of a
for ω = ce/c2b
〈. . . , . . .〉 = ‖u‖2 − 2ωa(u, u) + ω2‖RAu‖2 ≤ 1− 2ωce + ω2c2b < 1
since 〈RAu, u〉 = (Au)(u) = a(u, u)
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 10
![Page 30: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/30.jpg)
equivalent form of the variational problem
RAu = Rλ
rewrite as fixed point equation
u = u − ωRAu︸ ︷︷ ︸Su
+ωRλ
with ω > 0show: S is contraction for small ω( Banach’s fixed point theorem completes proof of first part)estimate
‖S‖ = sup‖u‖=1
〈u − ωRAu, u − ωRAu〉1/2
using the ellipticity of afor ω = ce/c
2b
〈. . . , . . .〉 = ‖u‖2 − 2ωa(u, u) + ω2‖RAu‖2 ≤ 1− 2ωce + ω2c2b < 1
since 〈RAu, u〉 = (Au)(u) = a(u, u)SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 10
![Page 31: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/31.jpg)
Proof (symmetric case)
symmetric elliptic bilinear form a scalar product and equivalent norm‖ · ‖a:
‖u‖2a = 〈u, u〉a = a(u, u) � ‖u‖2, u ∈ H
Riesz representation theorem =⇒
λ(v) = a(Rλ, v), v ∈ H
rewrite the quadratic form:
Q(u) =1
2a(u, u)− λ(u) =
1
2‖u − Rλ‖2a −
1
2‖Rλ‖2a
minimizing Q ⇔ best approximation to Rλ from Vcharacterization of the best approximation u,
a(u − Rλ, v) = 0, v ∈ V ,
⇔ variational equations
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 11
![Page 32: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/32.jpg)
Proof (symmetric case)
symmetric elliptic bilinear form a scalar product and equivalent norm‖ · ‖a:
‖u‖2a = 〈u, u〉a = a(u, u) � ‖u‖2, u ∈ H
Riesz representation theorem =⇒
λ(v) = a(Rλ, v), v ∈ H
rewrite the quadratic form:
Q(u) =1
2a(u, u)− λ(u) =
1
2‖u − Rλ‖2a −
1
2‖Rλ‖2a
minimizing Q ⇔ best approximation to Rλ from Vcharacterization of the best approximation u,
a(u − Rλ, v) = 0, v ∈ V ,
⇔ variational equations
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 11
![Page 33: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/33.jpg)
Proof (symmetric case)
symmetric elliptic bilinear form a scalar product and equivalent norm‖ · ‖a:
‖u‖2a = 〈u, u〉a = a(u, u) � ‖u‖2, u ∈ H
Riesz representation theorem =⇒
λ(v) = a(Rλ, v), v ∈ H
rewrite the quadratic form:
Q(u) =1
2a(u, u)− λ(u) =
1
2‖u − Rλ‖2a −
1
2‖Rλ‖2a
minimizing Q ⇔ best approximation to Rλ from Vcharacterization of the best approximation u,
a(u − Rλ, v) = 0, v ∈ V ,
⇔ variational equations
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 11
![Page 34: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/34.jpg)
Proof (symmetric case)
symmetric elliptic bilinear form a scalar product and equivalent norm‖ · ‖a:
‖u‖2a = 〈u, u〉a = a(u, u) � ‖u‖2, u ∈ H
Riesz representation theorem =⇒
λ(v) = a(Rλ, v), v ∈ H
rewrite the quadratic form:
Q(u) =1
2a(u, u)− λ(u) =
1
2‖u − Rλ‖2a −
1
2‖Rλ‖2a
minimizing Q ⇔ best approximation to Rλ from V
characterization of the best approximation u,
a(u − Rλ, v) = 0, v ∈ V ,
⇔ variational equations
SIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 11
![Page 35: Abstract Variational Problem · 2012-04-13 · SIAM FR26: FEM with B-Splines Basic Finite Element Concepts { Abstract Variational Problems 2-4, page 1. Abstract Variational Problem](https://reader034.fdocument.org/reader034/viewer/2022042804/5f59fdc573dc025b14298430/html5/thumbnails/35.jpg)
Proof (symmetric case)
symmetric elliptic bilinear form a scalar product and equivalent norm‖ · ‖a:
‖u‖2a = 〈u, u〉a = a(u, u) � ‖u‖2, u ∈ H
Riesz representation theorem =⇒
λ(v) = a(Rλ, v), v ∈ H
rewrite the quadratic form:
Q(u) =1
2a(u, u)− λ(u) =
1
2‖u − Rλ‖2a −
1
2‖Rλ‖2a
minimizing Q ⇔ best approximation to Rλ from Vcharacterization of the best approximation u,
a(u − Rλ, v) = 0, v ∈ V ,
⇔ variational equationsSIAM FR26: FEM with B-Splines Basic Finite Element Concepts – Abstract Variational Problems 2-4, page 11