4Elliptic boundary value problems

III Analysis of Partial Differential Equations



4.2 The Fredholm alternative
To understand the second problem, we shall seek to prove the following theorem:
Theorem (Fredholm alternative). Consider the problem
Lu = f, u|
U
= 0. ()
For
L
a uniformly elliptic operator on an open bounded set
U
with
C
1
boundary,
either
(i)
For each
f L
2
(
U
), there is a unique weak solution
u H
1
0
(
U
) to (
); or
(ii)
There exists a non-zero weak solution
u H
1
0
(
U
) to the homogeneous
problem, i.e. () with f = 0.
This is similar to what we know about solving matrix equations
Ax
=
b
either there is a solution for all
b
, or there are infinitely many solutions to the
homogneous problem.
Similar to the previous theorem, this follows from some general functional
analytic result. Recall the definition of a compact operator:
Definition
(Compact operator)
.
A bounded operator
K
:
H H
0
is compact
if every bounded sequence (
u
m
)
m=1
has a subsequence
u
m
j
such that (
Ku
m
j
)
j=1
converges strongly in H.
Recall (or prove as an exercise) the following theorem regarding compact
operators.
Theorem
(Fredholm alternative)
.
Let
H
be a Hilbert space and
K
:
H H
be a compact operator. Then
(i) ker(I K) is finite-dimensional.
(ii) im(I K) is finite-dimensional.
(iii) im(I K) = ker(I K
)
.
(iv) ker(I K) = {0} iff im(I K) = H.
(v) dim ker(I K) = dim ker(I K
) = dim coker(I K).
How do we apply this to our situation? Our previous theorem told us that
L
+
γ
is invertible for large
γ
, and we claim that (
L
+
γ
)
1
is compact. We can
then deduce the previous result by applying (iv) of the Fredholm alternative
with K a (scalar multiple of) (L + γ)
1
(plus some bookkeeping).
So let us show that (
L
+
γ
)
1
is compact. Note that this maps sends
f L
2
(
U
)
to
u H
1
0
(
U
). To make it an endomorphism, we have to compose this with the
inclusion
H
1
0
(
U
)
L
2
(
U
). The proof that (
L
+
γ
)
1
is compact will not involve
(
L
+
γ
)
1
in any way we shall show that the inclusion
H
1
0
(
U
)
L
2
(
U
) is
compact!
We shall prove this in two steps. First, we need the notion of weak conver-
gence.
Definition
(Weak convergence)
.
Suppose (
u
n
)
n=1
is a sequence in a Hilbert
space H. We say u
n
converges weakly to u H if
(u
n
, w) (u, w)
for all w H. We write u
n
* u.
Of course, we have
Lemma. Weak limits are unique.
Lemma. Strong convergence implies weak convergence.
We shall show that given any bounded sequence in
H
1
0
(
U
), we can find
a subsequence that is weakly convergent. We then show that every weakly
convergent sequence in H
1
0
(U) is strongly convergent in L
2
(U).
In fact, the first result is completely general:
Theorem
(Weak compactness)
.
Let
H
be a separable Hilbert space, and suppose
(
u
m
)
m=1
is a bounded sequence in
H
with
ku
m
k K
for all
m
. Then
u
m
admits
a subsequence (u
m
j
)
j=1
such that u
m
j
* u for some u H with kuk K.
One can prove this theorem without assuming
H
is separable, but it is slightly
messier.
Proof.
Let (
e
i
)
i=1
be an orthonormal basis for
H
. Consider (
e
1
, u
m
). By Cauchy–
Schwarz, we have
|(e
1
, u
m
)| ke
1
kke
m
k K.
So by Bolzano–Weierstrass, there exists a subsequence (
u
m
j
) such that (
e
1
, u
m
j
)
converges.
Doing this iteratively, we can find a subsequence (
v
`
) such that for each
i
,
there is some c
i
such that (e
i
, v
`
) c
i
as ` .
We would expect the weak limit to be
P
c
i
e
i
. To prove this, we need to first
show it converges. We have
p
X
j=1
|c
j
|
2
= lim
k→∞
p
X
j=1
|(e
j
, v
`
)|
2
sup
p
X
j=1
|(e
j
, v
`
)|
2
sup kv
k
k
2
K
2
,
using Bessel’s inequality. So
u =
X
j=1
c
j
e
j
converges in H, and kuk K. We already have
(e
j
, v
`
) (e
j
, u)
for all
j
. Since
kv
`
uk
is bounded by 2
K
, it follows that the set of all
w
such
that
(w, v
`
) (v, u) ()
is closed under finite linear combinations and taking limits, hence is all of
H
.
To see that it is closed under limits, suppose
w
k
w
, and
w
k
satisfy (
). Then
|(w, v
`
)(w, u)| |(w w
k
, v
`
u)|+|(w
k
, v
`
u)| 2Kkww
k
k+|(w
k
, v
`
u)|
So we can first find
k
large enough such that the first term is small, then pick
`
such that the second is small.
We next want to show that if
u
m
* u
in
H
1
(
U
), then
u
m
u
in
L
1
. We
may as well assume that
U
is some large cube of length
L
by extension. Notice
that since
U
is bounded, the constant function 1 is in
H
1
0
(
U
). So
u
m
* u
in
particular implies
R
U
(u
m
u) dx 0.
Recall that the Poincar´e inequality tells us if
u H
1
0
(
U
), then we can bound
kuk
L
2
(Q)
by some multiple of
k
D
uk
L
2
(U)
. If we try to prove this without the
assumption that
u
vanishes on the boundary, then we find that we need a
correction term. The resulting lemma is as follows:
Lemma
(Poincar´e revisited)
.
Suppose
u H
1
(
R
n
). Let
Q
= [
ξ
1
, ξ
1
+
L
]
×···×
[ξ
n
, ξ
n
+ L] be a cube of length L. Then we have
kuk
2
L
2
(Q)
1
|Q|
Z
Q
u(x) dx
2
+
nL
2
2
kDuk
2
L
2
(Q)
.
We can improve this to obtain better bounds by subdividing
Q
into smaller
cubes, and then applying this to each of the cubes individually. By subdividing
enough, this leads to a proof that u
m
* u in H
1
implies u
m
u in H
0
.
Proof. By approximation, we can assume u C
(
¯
Q). For x, y Q, we write
u(x) u(y) =
Z
x
1
y
1
d
dt
u(t, x
1
, . . . , x
n
) dt
+
Z
x
2
y
2
d
dt
u(y
1
, t, x
3
, . . . , x
n
) dt
+ ···
+
Z
x
n
y
n
d
dt
u(y
1
, . . . , y
n1
, t) dt.
Squaring, and using 2ab a
2
+ b
2
, we have
u(x)
2
+ u(y)
2
2u(x)u(y) n
Z
x
1
y
1
d
dt
u(t, x
1
, . . . , x
n
) dt
2
+ ···
+ n
Z
x
n
y
n
d
dt
u(y
1
, . . . , y
n1
, t) dt
2
.
Now integrate over x and y. On the left, we get
ZZ
Q×Q
dx dy (u(x)
2
+ u(y)
2
2u(x)u(y)) = 2|Q|kuk
2
L
2
(Q)
2
Z
Q
u(x) dx
2
.
On the right we have
I
1
=
Z
x
1
y
1
d
dt
u(t, x
2
, . . . , x
n
) dt
2
Z
x
1
y
1
dt
Z
x
1
y
1
d
dt
u(t, x
2
, . . . , x
n
)
2
dt (Cauchy–Schwarz)
L
Z
ξ
1
+L
ξ
1
d
dt
u(t, x
2
, . . . , x
n
)
2
dt.
Integrating over all x, y Q, we get
ZZ
Q×Q
dx dy I
1
L
2
|Q|kD
1
uk
2
L
2
(Q)
.
Similarly estimating the terms on the right-hand side, we find that
2|Q|kuk
L
2
(Q)
2
Z
Q
u(x) dx
2
n|Q|
n
X
i=1
kD
i
uk
2
L
2
(Q)
= n|Q|L
2
kDuk
2
L
2
(Q)
.
It now follows that
Theorem
(Rellich–Kondrachov)
.
Let
U R
n
be open, bounded with
C
1
boundary. Then if (
u
m
)
m=1
is a sequence in
H
1
(
U
) with
u
m
* u
, then
u
m
u
in L
2
.
In particular, by weak compactness any sequence in
H
1
(
U
) has a subsequence
that is convergent in L
2
(U).
Note that to obtain the “in particular” part, we need to know that
H
1
(
U
) is
separable. This is an exercise on the example sheet. Alternatively, we can appeal
to a stronger version of weak compactness that does not assume separability.
Proof.
By the extension theorem, we may assume
U
=
Q
for some large cube
Q
with U b Q.
We subdivide
Q
into
N
many cubes of side length
δ
, such that the cubes
only intersect at their faces. Call these {Q
a
}
N
a=1
.
We apply Poincar´e separately to each of these to obtain
ku
j
uk
2
L
2
(Q)
=
N
X
a=1
ku
j
uk
2
L
2
(Q
a
)
N
X
a=1
"
1
|Q
a
|
Z
Q
a
(u
i
u) dx
2
+
2
2
kDu
i
Duk
2
L
2
(Q
a
)
#
=
N
X
a=1
1
|Q
a
|
Z
Q
a
(u
i
u) dx
2
+
2
2
kDu
i
Duk
2
L
2
(Q)
.
Now since
k
D
u
i
D
uk
2
L
2
(Q)
is fixed, for
δ
small enough, the second term is
<
ε
2
.
Then since u
i
* u, we in particular have
Z
Q
1
(u
i
u) dx 0 as i
for all
a
, since this is just the inner product with the constant function 1. So for
i large enough, the first term is also <
ε
2
.
The same result holds with
H
1
(
U
) replaced by
H
1
0
(
U
). The proof is in fact
simpler, and we wouldn’t need the assumption that the boundary is C
1
.
Corollary.
Suppose
K
:
L
2
(
U
)
H
1
(
U
) is a bounded linear operator. Then
the composition
L
2
(U) H
1
(U) L
2
(U)
K
is compact.
The slogan is that we get compactness whenever we improve regularity, which
is something that happens in much more generality.
Proof.
Indeed, if
u
m
L
2
(
U
) is bounded, then
Ku
m
is also bounded. So by
Rellich–Kondrachov, there exists a subsequence u
m
j
u in L
2
(U).
We are now ready to prove the Fredholm alternative for elliptic boundary
value problems. Recall that in our description of the Fredholm alternative, we
had the direct characterizations
im
(
I K
) =
ker
(
I K
)
. We can make the
analogous statement here. To do so, we need to talk about the adjoint of
L
.
Since
L
is not an operator defined on
L
2
(
U
), trying to write down what it means
to be an adjoint is slightly messy. Instead, we shall be content with talking
about “formal adjoints”.
It’s been a while since we’ve met a PDE, so let’s recall the setting we had.
We have a uniformly elliptic operator
Lu =
n
X
i,j=1
(a
ij
(x)u
x
j
)
x
i
+
n
X
i=1
b
i
(x)u
x
i
+ c(x)u
on an open bounded set U with C
1
boundary. The associated bilinear form is
B[u, v] =
Z
U
u
X
i,j
a
ij
(x)u
x
i
v
x
j
+
n
X
i=1
b
i
(x)u
x
i
v + c(x)uv
dx.
We are interested solving in the boundary value problem
Lu = f, u|
u
= 0
with f L
2
(U).
The formal adjoint of L is defined by the relation
(Lφ, ψ)
L
2
(U)
= (φ, L
ψ)
L
2
(U)
for all φ, ψ C
c
(U). By integration by parts, we know L
should be given by
L
v =
n
X
i,j=1
(a
ij
v
x
j
)
x
i
n
X
i=1
b
i
(x)v
x
j
+
c
n
X
i=1
b
i
x
i
!
v.
Note that here we have to assume that
b
i
C
1
(
¯
U
). However, what really
interests us is the adjoint bilinear form, which is simply given by
B
[v, u] = B[u, v].
We are actually just interested in
B
, and not
L
, and we can sensibly talk about
B
even if b
i
is not differentiable.
As usual, we say
v H
1
0
(
U
) is a weak solution of the adjoint problem
L
v = f, v|
U
= 0 if
B
[v, u] = (f, u)
L
2
(U)
for all u H
1
0
(U).
Given this set up, we can now state and prove the Fredholm alternative.
Theorem
(Fredholm alternative for elliptic BVP)
.
Let
L
be a uniformly elliptic
operator on an open bounded set U with C
1
boundary. Consider the problem
Lu = f, u|
U
= 0. ()
Then exactly one of the following are true:
(i) For each f L
2
(U), there is a unique weak solution u H
1
0
(U) to ()
(ii)
There exists a non-zero weak solution
u H
1
0
(
U
) to the homogeneous
problem, i.e. () with f = 0.
If this holds, then the dimension of
N
=
ker L H
1
0
(
U
) is equal to the
dimension of N
= ker L
H
1
0
(U).
Finally, () has a solution if and only if (f, v)
L
2
(U)
= 0 for all v N
Proof.
We know that there exists
γ >
0 such that for any
f L
2
(
U
), there is a
unique weak solution u H
1
0
(U) to
L
γ
u = Lu + γu = f, u|
U
= 0.
Moreover, we have the bound
kuk
H
1
(U)
Ckf k
L
2
(U)
(which gives uniqueness).
Thus, we can set
L
1
γ
f
to be this
u
, and then
L
1
γ
:
L
2
(
U
)
H
1
0
(
U
) is a
bounded linear map. Composing with the inclusion
L
2
(
U
), we get a compact
endomorphism of L
2
(U).
Now suppose u H
1
0
is a weak solution to (). Then
B[u, v] = (f, v)
L
2
(U)
for all v H
1
0
(U)
is true if and only if
B
γ
[u, v] B[u, v] + γ(u, v) = (f + γu, v) for all v H
1
0
(U).
Hence, u is a weak solution of () if and only if
u = L
1
γ
(f + γu) = γL
1
γ
u + L
1
γ
f.
In other words, u solves () iff
u Ku = h,
for
K = γL
1
γ
, h = L
1
γ
f.
Since we know that
K
:
L
2
(
U
)
L
2
(
U
) is compact, by the Fredholm alternative
for compact operators, either
(i) u Ku = h admits a solution u L
2
(U) for all h L
2
(U); or
(ii)
There exists a non-zero
u L
2
(
U
) such that
u Ku
= 0. Moreover,
im(I K) = ker(I K
)
and dim ker(I K) = dim im(I K)
.
There is a bit of bookkeeping to show that this corresponds to the two alternatives
in the theorem.
(i) We need to show that u H
1
0
(U). But this is trivial, since we have
u = γL
1
γ
u + L
1
γ
f,
and we know that L
1
γ
maps L
2
(U) into H
1
0
(U).
(ii)
As above, we know that the non-zero solution
u
. There are two things to
show. First, we have to show that
v K
v
= 0 iff
v
is a weak solution to
L
v = 0, v|
U
= 0.
Next, we need to show that h = L
1
γ
f (N
)
iff f (N
)
.
For the first part, we want to show that
v ker
(
I K
) iff
B
[
v, u
] =
B[u, v] = 0 for all u H
1
0
(U).
We are good at evaluating
B
[
u, v
] when
u
is of the form
L
1
γ
w
, by definition
of a weak solution. Fortunately,
im L
1
γ
contains
C
c
(
U
), since
L
1
γ
L
γ
φ
=
φ
for all
φ C
c
(
U
). In particular,
im L
1
γ
is dense in
H
1
0
(
U
). So it
suffices to show that
v ker
(
I K
) iff
B
[
L
1
γ
w, v
] = 0 for
w L
2
(
U
).
This is immediate from the computation
B[L
1
γ
w, v] = B
γ
[L
1
γ
w, v]γ(L
1
γ
w, v) = (w, v)(Kw, v) = (w, vK
v).
The second is also easy if v N
= ker(I K
), then
(L
1
γ
f, v) =
1
γ
(Kf, v) =
1
γ
(f, K
v) =
1
γ
(f, v).