3Function spaces
III Analysis of Partial Differential Equations
3.3 Approximation of functions in Sobolev spaces
It would be nice if we could approximate functions in
W
k,p
(
U
) with something
more tractable. For example, it would be nice if we could approximate them by
smooth functions, so that the weak derivatives are genuine derivatives. A useful
trick to improve regularity of a function is to convolve with a smooth mollifier.
Definition (Standard mollifier). Let
η(x) =
(
Ce
1/(|x|
2
−1)
|x| < 1
0 |x| ≥ 1
,
where C is chosen so that
R
R
n
η(x) dx = 1.
One checks that this is a smooth function on R
n
, peaked at x = 0.
For each ε > 0, we set
η
ε
(x) =
1
ε
n
η
x
ε
.
Of course, the pre-factor of
1
ε
n
is chosen so that
η
ε
is appropriately normalized.
We call η
ε
the standard mollifier, and it satisfies supp η
ε
⊆ B
ε
(0).
We think of these η
ε
as approximations of the δ-function.
Now suppose U ⊆ R
n
is open, and let
U
ε
= {x ∈ U : dist(x, ∂U) > ε}.
Definition
(Mollification)
.
If
f ∈ L
1
loc
(
U
), we define the mollification
f
ε
:
U
ε
→
R by the convolution
f
ε
= η
ε
∗ f.
In other words,
f
ε
(x) =
Z
U
η
ε
(x − y)f(y) =
Z
B
ε
(x)
η
ε
(x − y)f(y) dy.
Thus,
f
ε
is the “local average” of
f
around each point, with the weighting
given by
η
ε
. The hope is that
f
ε
will have much better regularity properties
than f.
Theorem. Let f ∈ L
1
loc
(U). Then
(i) f
ε
∈ C
∞
(U
ε
).
(ii) f
ε
→ f almost everywhere as ε → 0.
(iii) If in fact f ∈ C(U), then f
ε
→ f uniformly on compact subsets.
(iv)
If 1
≤ p < ∞
and
f ∈ L
p
loc
(
U
), then
f
ε
→ f
in
L
p
loc
(
U
), i.e. we have
convergence in L
p
on any V b U.
In general, the difficulty of proving these approximation theorems lie in what
happens at the boundary
Lemma. Assume u ∈ W
k,p
(U) for some 1 ≤ p < ∞, and set
u
ε
= η
ε
∗ u on U
ε
.
Then
(i) u
ε
∈ C
∞
(U
ε
) for each ε > 0
(ii) If V b U, then u
ε
→ u in W
k,p
(V ).
Proof.
(i) As above.
(ii) We claim that
D
α
u
ε
= η
ε
∗ D
α
u
for |α| ≤ k in U
ε
.
To see this, we have
D
α
u
ε
(x) = D
α
Z
U
η
ε
(x − y)u(y) dy
=
Z
U
D
α
x
η
ε
(x − y)u(y) dy
=
Z
U
(−1)
|α|
D
α
y
η
ε
(x − y)u(y) dy
For a fixed
x ∈ U
ε
,
η
ε
(
x − ·
)
∈ C
∞
c
(
U
), so by the definition of a weak
derivative, this is equal to
=
Z
U
η
ε
(x − y)D
α
u(y) dy
= η
ε
∗ D
α
u.
It is an exercise to verify that we can indeed move the derivative past the
integral.
Thus, if we fix
V b U
. Then by the previous parts, we see that D
α
u
ε
→
D
α
u in L
p
(V ) as ε → 0 for |α| ≤ k. So
ku
ε
− uk
p
W
k.p
(V )
=
X
|α|≤k
kD
α
u
ε
− D
α
uk
p
L
p
(V )
→ 0
as ε → 0.
Theorem
(Global approximation)
.
Let 1
≤ p < ∞
, and
U ⊆ R
n
be open and
bounded. Then C
∞
(U) ∩ W
k,p
(U) is dense in W
k,p
(U).
Our main obstacle to overcome is the fact that the mollifications are only
defined on U
ε
, and not U.
Proof. For i ≥ 1, define
U
i
=
x ∈ U | dist(x, ∂U) >
1
i
V
i
= U
i+3
−
¯
U
i+1
W
i
= U
i+4
−
¯
U
i
.
We clearly have
U
=
S
∞
i=1
U
i
, and we can choose
V
0
b U
such that
U
=
S
∞
i=0
V
i
.
Let
{ζ
i
}
∞
i=0
be a partition of unity subordinate to
{V
i
}
. Thus, we have
0 ≤ ζ
i
≤ 1, ζ
i
∈ C
∞
c
(V
i
) and
P
∞
i=0
ζ
i
= 1 on U.
Fix δ > 0. Then for each i, we can choose ε
i
sufficiently small such that
u
i
= η
ε
i
∗ ζ
i
u
satisfies supp u
i
⊆ W
i
and
ku
i
− ζ
i
uk
W
k.p
(U)
= ku
i
− ζ
i
uk
W
k.p
(W
i
)
≤
δ
2
i+1
.
Now set
v =
∞
X
i=0
u
i
∈ C
∞
(U).
Note that we do not know (yet) that
v ∈ W
k.p
(
U
). But it certainly is when we
restrict to some V b U.
In any such subset, the sum is finite, and since u =
P
∞
i=0
ζ
i
u, we have
kv −uk
W
k,p
(V )
≤
∞
X
i=0
ku
i
− ζ
i
uk
W
k.p
(V )
≤ δ
∞
X
i=0
2
−(i+1)
= δ.
Since the bound
δ
does not depend on
V
, by taking the supremum over all
V
,
we have
kv −uk
W
k.p
(U)
≤ δ.
So we are done.
It would be nice for
C
∞
(
¯
U
) to be dense, instead of just
C
∞
(
U
). It turns out
this is possible, as long as we have a sensible boundary.
Definition
(
C
k,δ
boundary)
.
Let
U ⊆ R
n
be open and bounded. We say
∂U
is
C
k,δ
if for any point in the boundary
p ∈ ∂U
, there exists
r >
0 and a function
γ ∈ C
k,δ
(
R
n−1
) such that (possibly after relabelling and rotating axes) we have
U ∩ B
r
(p) = {(x
0
, x
n
) ∈ B
r
(p) : x
n
> γ(x
0
)}.
Thus, this says our boundary is locally the graph of a C
k,δ
function.
Theorem
(Smooth approximation up to boundary)
.
Let 1
≤ p < ∞
, and
U ⊆ R
n
be open and bounded. Suppose
∂U
is
C
0,1
. Then
C
∞
(
¯
U
)
∩ W
k,p
(
U
) is
dense in W
k,p
(U).
Proof.
Previously, the reason we didn’t get something in
C
∞
(
¯
U
) was that we
had to glue together infinitely many mollifications whose domain collectively
exhaust
U
, and there is no hope that the resulting function is in
C
∞
(
¯
U
). In the
current scenario, we know that U locally looks like
x
0
The idea is that given a
u
defined on
U
, we can shift it downwards by some
ε
.
It is a known result that translation is continuous, so this only changes
u
by a
tiny bit. We can then mollify with a
¯ε < ε
, which would then give a function
defined on U (at least locally near x
0
).
So fix some
x
0
∈ ∂U
. Since
∂U
is
C
0,1
, there exists
r >
0 such that
γ ∈ C
0,1
(R
n−1
) such that
U ∩ B
r
(x
0
) = {(x
0
, x
n
) ∈ B
r
(x
0
) | x
n
> γ(x
0
)}.
Set
V = U ∩ B
r/2
(x
0
).
Define the shifted function u
ε
to be
u
ε
(x) = u(x + εe
n
).
Now pick ¯ε sufficiently small such that
v
ε,¯ε
= η
¯ε
∗ u
ε
is well-defined. Note that here we need to use the fact that
∂U
is
C
0,1
. Indeed,
we can see that if the slope of ∂U is very steep near a point x:
ε
then we need to choose a
¯ε
much smaller than
ε
. By requiring that
γ
is 1-H¨older
continuous, we can ensure there is a single choice of
¯ε
that works throughout
V
.
As long as ¯ε is small enough, we know that v
ε,¯ε
∈ C
∞
(
¯
V ).
Fix δ > 0. We can now estimate
kv
ε,˜ε
− uk
W
k.p
(V )
= kv
ε,˜ε
− u
ε
+ u
ε
− uk
W
k,p
(V )
≤ kv
ε,˜ε
− u
ε
k
W
k,p
(V )
+ ku
ε
− uk
W
k.p
(V )
.
Since translation is continuous in the
L
p
norm for
p < ∞
, we can pick
ε >
0
such that
ku
ε
− uk
W
k.p
(V )
<
δ
2
. Having fixed such an
ε
, we can pick
˜ε
so small
that we also have kv
ε,˜ε
− u
ε
k
W
k.p
(V )
<
δ
2
.
The conclusion of this is that for any
x
0
∈ ∂U
, we can find a neighbourhood
V ⊆ U
of
x
0
in
U
such that for any
u ∈ W
k,p
(
U
) and
δ >
0, there exists
v ∈ C
∞
(
¯
V ) such that ku − vk
W
k,p
(V )
≤ δ.
It remains to patch all of these together using a partition of unity. By the
compactness of
∂U
, we can cover
∂U
by finitely many of these
V
, say
V
1
, . . . , V
N
.
We further pick a V
0
such that V
0
b U and
U =
N
[
i=0
V
i
.
We can pick approximations
v
i
∈ C
∞
(
¯
V
i
) for
i
= 0
, . . . , N
(the
i
= 0 case is given
by the previous global approximation theorem), satisfying
kv
i
− uk
W
k,p
(V
i
)
≤ δ
.
Pick a partition of unity {ζ
i
}
N
i=0
of
¯
U subordinate to {V
i
}. Define
v =
N
X
i=0
ζ
i
v
i
.
Clearly v ∈ C
∞
(
¯
U), and we can bound
kD
α
v −D
α
uk
L
p
(U)
=
D
α
N
X
i=0
ζ
i
v
i
− D
α
N
X
i=0
ζ
i
u
L
p
(U)
≤ C
k
N
X
i=0
kv
i
− uk
W
k.p
(V
i
)
≤ C
k
(1 + N)δ,
where
C
k
is a constant that solely depends on the derivatives of the partition of
unity, which are fixed. So we are done.