INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES
WARSZAWA 1992
COEFFICIENTS OF THE SINGULARITIES ON DOMAINS WITH CONICAL POINTS
M O N I Q U E D A U G E
D´ epartement de Math´ ematiques, Universit´ e de Nantes, U.R.A. C.N.R.S. 758 2, rue de la Houssini` ere, F-44072 Nantes Cedex 03, France
S E R G E N I C A I S E
U.F.R. de Math´ ematiques Pures et Appliqu´ ees, Universit´ e des Sciences et Techniques de Lille Flandres Artois, U.R.A. C.N.R.S. 751, F-59655 Villeneuve D’Ascq Cedex, France
Abstract. As a model for elliptic boundary value problems, we consider the Dirichlet prob- lem for an elliptic operator. Solutions have singular expansions near the conical points of the domain. We give formulas for the coefficients in these expansions.
1. Introduction. We consider bounded n-dimensional domains with conical points, as Kondrat’ev in [4]. For simplicity, we suppose that there is only one conical point and that it is located at 0. We denote by Ω our domain and we assume that its boundary is C
∞outside 0 and that it coincides with a cone Γ in a neighborhood of 0. We denote by x the cartesian coordinates in R
nand by (r, θ) the spherical coordinates. The spherical section of Γ is denoted by G:
Γ ∩ S
n−1=: G .
We are interested in the Dirichlet boundary value problem for an elliptic operator P (x; D
x) of order 2m. We assume that the coefficients of this operator are C
∞(Ω \ 0). We have to sharpen this assumption. We will consider three cases (C1), (C2) and (C3), each of them being more general than the previous one:
• (C1): P is homogeneous with constant coefficients; then there exists an operator L with C
∞(G) coefficients such that
P (D
x) = r
−2mL(θ; r∂
r, ∂
θ) .
• (C2): P has C
∞(Ω) coefficients; then, if L denotes the principal part of
[91]
P (0; D
x), then L satisfies the assumption of (C1) and the difference R(x; D
x) := P (x; D
x) − L(D
x)
is a remainder.
• (C3): there exists an operator L with C
∞(G) coefficients such that the difference R(x; D
x) := P (x; D
x) − r
−2mL(θ; r∂
r, ∂
θ)
is a remainder in a sense we are going to explain.
The Coulomb operator −∆ + 1/r satisfies the assumptions of (C3). To explain what we mean by remainder, we need some weighted Sobolev spaces.
As usual, H
om(Ω) denotes the closure of D(Ω) in H
m(Ω) and H
−m(Ω) is its dual space. For any positive integer k and any real β, H
βk(Ω) is defined as
H
βk(Ω) = {u ∈ D
0(Ω) | r
β−k+|α|D
αxu ∈ L
2(Ω) ∀α, |α| ≤ k} .
We also define H
βs(Ω) for any positive real s in a natural way (cf. for instance appendix A in [1]), and for any negative s by duality.
For any s > 0 and any β, the operators P and L in the case (C2), and r
−2mL in the cases (C1) and (C3) are continuous H
βs+m(Ω) → H
βs−m(Ω). Moreover, in the case (C2), the remainder R is continuous H
βs+m(Ω) → H
β−s−m(Ω).
Now, the assumption in the case (C3) is that there exists δ ∈ ]0, 1] such that for any s ≥ 0 and β the remainder R is continuous
(1.1) H
βs+m(Ω) → H
β−δs−m(Ω) . If a
αdenote the coefficients of R:
R(x; D
x) = X
|α|≤2m
a
αD
xα, the assumption (1.1) holds if
∀γ ∈ N
n, D
xγa
αis O(r
|α|−2m+δ−|γ|) .
With the above assumptions, we are interested in the structure of any solution u of the Dirichlet problem
(1.2) u ∈ H
om(Ω), P u ∈ H
βs−m(Ω), with s > 0 and s − β > 0 .
Since s > 0 and s − β > 0, H
βs−m(Ω) is (compactly) embedded in H
−m(Ω). So P u has in a sense more regularity than u. Of course, it is possible to consider more general situations than (1.2), i.e. to assume that u belongs to some weighted space. The solution of this problem would be essentially the same as for (1.2). We chose (1.2), because it is the natural framework when P is strongly elliptic.
The assumptions of the case (C2) are those of Kondrat’ev [4]. We also took these assumptions in our earlier works [2] and [3]. The assumptions of the case (C3) were introduced by Maz’ya and Plamenevski˘ı [7]. Kondrat’ev proved the existence of an expansion of the solutions of (1.2) in the form of a sum P
i
c
iS
iwhere the S
ionly depend on Ω and P and the c
iare some coefficients. Maz’ya and
Plamenevski˘ı in [6] and [7] gave formulas for these coefficients; we also studied them in [2] and [3] in a different framework. What we give here extends in a certain sense [7] and [2].
To end this section, let us state a Fredholm theorem. Such a result is re- lated to asymptotics of solutions: a solution of (1.2) can be split into a singular part (asymptotics) and a regular part (remainder) when the assumptions of the following theorem hold.
We need some notations. We denote by P
βsthe operator P
βs: H
oβ−sm(Ω) ∩ H
βs+m(Ω) → H
βs−m(Ω), u 7→ P u ,
where H
oγm(Ω) is the closure of D(Ω) in H
γm(Ω). We will simply denote by P
βthe operator P
β0.
P
−β∗denotes the adjoint of P
β. It acts
P
−β∗: H
o−βm(Ω) → H
−β−m(Ω) . For any λ ∈ C, L(λ) is the operator
L(θ; λ, ∂
θ) : H
om(G) → H
−m(G) .
L(λ) is one-to-one except when λ belongs to a countable set in C, which can be called the spectrum of L and is denoted by Sp(L).
Theorem 1.1. In the case (C3), assume that s ≥ 0, β ∈ R, s − β ≥ 0 and that
∀λ such that Re λ = s + m − β − n/2, λ 6∈ Sp(L) . Then P
βsis a Fredholm operator.
2. The model problem. In this section, we will only study the case (C1), when the operator is homogeneous with constant coefficients. We recall that then P = L.
For each λ ∈ Sp(L), the space Z
λ:= n
u
u = r
λX
q
Log
qr u
q(θ), u
q∈ H
om(G), Lu = 0 o
does not reduce to 0: all functions of the form r
λu
0where u
0∈ Ker L(λ) are in Z
λ. Let σ
νλ, for ν = 1, . . . , N
λ, denote a basis of this space.
Theorem 2.1. In the case (C1), assume the same hypotheses about s and β as
in Theorem 1.1. Let η be a cut-off function which is equal to 1 in a neighborhood
of 0 and has its support in another neighborhood of 0 where Ω coincides with the
cone Γ . Assume that u ∈ H
om(Ω) is such that P u ∈ H
βs−m(Ω). Then there exist
coefficients c
λνsuch that
u − X
λ∈Sp(L)
m−n/2<Re λ<s−β+m−n/2 Nλ
X
ν=1
c
λνησ
νλ∈ H
βs+m(Ω) .
If Ker P
0⊆ H
βs+m(Ω), then the coefficients c
λνonly depend on P u. This is the reason for the introduction of the following assumption.
(2.1) If u ∈ H
om(Ω) is such that P u = 0 then u ∈ H
β−sm(Ω).
If (2.1) holds, as a consequence of a well-known regularity result for corner prob- lems, such an element of the kernel belongs to any space H
β−s+tm+t(Ω) for any t ≥ 0 (see Theorem 2.9 of [2] for instance).
Now, we are going to construct dual singular functions. We start with a result from [6].
Lemma 2.2. For all λ ∈ Sp(L), there exists a basis τ
νλ, for ν = 1, . . . , N
λ, of the space
Z e
λ:= n u
u = r
−¯λ+2m−nX
q
Log
qrv
q(θ), v
q∈ H
om(G), L
∗u = 0 o such that ∀µ, µ
0∈ Sp(L), ∀ν, ν
0R
Ω
L(ησ
µν)τ
νµ00= δ
µµ0δ
νν0. See [6] and [2] for more details.
We set T
νλ:= ητ
νλ. We have
∀λ, m − n/2 < Re λ < s − β + m − n/2 ,
T
νλ6∈ H
0m(Ω) and T
νλ∈ H
os−βm(Ω).
Due to assumption (2.1), there exists Y
νλ∈ H
om(Ω) such that P
∗Y
νλ= P
∗T
νλ. We set
K
νλ:= T
νλ− Y
νλ.
Theorem 2.3. In the case (C1) and with the hypothesis (2.1), assume the same hypotheses about s and β as in Theorems 1.1 and 2.1. Then
c
λν= R
Ω
P uK
νλdx .
When P = ∆ and when Γ is a plane sector with opening ω, the σ
λνare the functions r
kπ/ωsin(kπθ/ω) for k ∈ N
∗and the τ
νλare the functions −(kπ)
−1× r
−kπ/ωsin(kπθ/ω).
3. The general problem. Now, we will work in the framework of the general
case (C3). All the above results can be extended in a certain sense to the case
(C3). We are going to introduce auxiliary functions. In the case (C2), the structure of these functions is more precisely known.
First, we construct elements of the kernel of P
∗in the same way by subtracting a corrective function Y
νλfrom T
νλ. The difference lies in the construction of Y
νλ. They cannot be found in H
om(Ω) in general but in a larger space.
Proposition 3.1. In the case (C3) and with the hypothesis (2.1), let λ ∈ Sp(L) such that m − n/2 < Re λ. Set γ(λ) := Re λ − m + n/2. Then ∀ε > 0, we have (3.1) T
νλ∈ H
γ(λ)+εm(Ω) and T
νλ6∈ H
γ(λ)m(Ω) .
Let δ
0= min{δ, γ(λ)}, where δ was introduced in (1.1). Then there exists Y
νλwhich satisfies the homogeneous Dirichlet conditions and such that
P
∗T
νλ= P
∗Y
νλand ∀ε > 0, Y
νλ∈ H
γ(λ)−δm 0+ε(Ω) .
Set K
νλ:= T
νλ− Y
νλ. The K
νλfor λ ∈ Sp(L), m − n/2 < Re λ < s − β + m − n/2 and for ν = 1, . . . , N
λform a basis of Ker P
s−β∗/ Ker P
0∗.
R e m a r k 3.2. In the case (C2), the Y
νλcan be constructed as a sum of terms T
ν,jλ= ηr
−λ+2m−n+jX
Log
qrv
ν,j,qλ(θ)
with 1 ≤ j ≤ Re λ − m + n/2 and of an element X
νλ∈ H
om(Ω). In the case (C1), the T
ν,jλare zero and Y
νλ= X
νλ(see §4 of [2]).
P r o o f. First step. Let us prove the existence of Y
νλ. By construction, L
∗T
νλ= 0 in a neighborhood of 0; as a consequence of the assumption (1.1), P
∗T
νλ∈ H
γ(λ)+ε−δ−m 0(Ω). We want to prove that
(3.2) P
∗T
νλ∈ Rg P
γ(λ)+ε−δ∗ 0.
But the regularity of T
νλyields P
∗T
νλ∈ Rg P
γ(λ)+ε∗. We choose ε small enough such that the ranges of P
γ(λ)+ε∗and of P
γ(λ)+ε−δ∗ 0are closed. We have
Rg P
γ(λ)+ε∗= (Ker P
−γ(λ)−ε)
⊥, Rg P
γ(λ)+ε−δ∗ 0= (Ker P
−γ(λ)−ε+δ0)
⊥. The hypothesis (2.1) yields that Ker P
−γ(λ)−ε= Ker P
−γ(λ)−ε+δ0. So, we have obtained (3.2).
Second step. Let us prove that the K
νλare linearly independent modulo H
om(Ω). Suppose that there exist nonzero coefficients c
λνsuch that P c
λνK
νλ∈
o
H
m(Ω). Let ξ be the largest real part of the λ which are associated with a nonzero coefficient. Since the whole sum belongs to H
om(Ω), we deduce by construction of the K
νλthat
∃% > 0, X
Re λ=ξ
c
λνT
νλ∈ H
ξ−m+n/2−%m(Ω) .
The form of the T
νλ(cf. (3.1)) allows us to show that the coefficients in the above
sum are all zero. We have obtained a contradiction.
Third step. Let γ be s − β and let n
γbe the cardinal of the set {K
νλ| m − n/2 < Re λ < s − β + m − n/2 and ν = 1, . . . , N
λ} .
We have to show that the dim Ker P
γ∗= dim Ker P
0∗+ n
γ. We rely on an index calculus. Choose γ
0, . . . , γ
Jsuch that
(3.3)
0 ≤ γ
0≤ . . . ≤ γ
J= γ ,
∀j = 1, . . . , J, γ
j− γ
j−1≤ δ ,
Sp(L) ∩ {λ ∈ C | m − n/2 < Re λ < m − n/2 + γ
0} = ∅ ,
∀j = 1, . . . , J, Sp(L) ∩ {λ ∈ C | Re λ = m − n/2 + γ
j} = ∅.
For each j = 1, . . . , J , the functions u ∈ H
o−γmj−1(Ω) such that Lu ∈ H
−γ−mj(Ω) can be written as a sum of a regular part in H
−γmj(Ω) and a singular part which is a combination of the ησ
νλwith λ ∈ Sp(L) and m−n/2+γ
j−1< Re λ < m−n/2+γ
j. Due to (1.1), the same holds for the operator P . Applying the result of appendix B of [1] for each pair (P
−γj−1, P
−γj) and summing over j = 1, . . . , J , we get
Ind P
−γ0− Ind P
−γ= n
γ.
As a consequence of the assumption (2.1), Ker P
−γ0= Ker P
−γ. Then Codim Rg P
−γ− Codim Rg P
−γ0= n
γ.
So for the adjoints, we get
dim Ker P
γ∗− dim Ker P
γ∗0= n
γ.
We end the proof by noting that the construction of γ
0implies Ker P
γ∗0= Ker P
0∗. We are now going to construct the singularities, i.e. a basis of functions be- longing to H
om(Ω), which are not in H
βs+m(Ω) and such that P u ∈ H
βs−m(Ω). In the case (C1), such a basis is formed by the ησ
νλ(cf. Theorem 2.1). Such a result extends to the case (C3) only if s − β ≤ δ. Let us state that with s = 0:
Lemma 3.3. In the case (C3), let τ and τ
0be such that 0 < τ
0− τ ≤ δ. Assume that
Sp(L) ∩ {λ ∈ C | Re λ = m − n/2 + τ } = ∅
and that u ∈ H
oτm0(Ω) is such that P u ∈ H
τ−m(Ω). Then there exist coefficients c
λνsuch that
u − X
λ∈Sp(L)
m−n/2−τ0<Re λ<m−n/2−τ Nλ
X
ν=1
c
λνησ
νλ∈ H
τm(Ω) .
As we already explained in the above proof, this is a simple consequence of the assumptions of (1.1) and of the corresponding result for L which is known [4].
In the case (C2), when s − β > 1, P (ησ
νλ) does not belong to H
βs−m(Ω) in general but there exist
σ
ν,jλ= r
λ+jX
Log
qr u
λν,j,q(θ)
where u
λν,j,q∈ H
om(G) and such that P h
η
σ
λν+ X
1≤j≤s−β+m−n/2
σ
λν,ji
∈ H
βs−m(Ω) (see §4.B of [2]).
In the general case (C3), we have another construction, which is less explicit, as in the previous Proposition 3.1.
Proposition 3.4. In the case (C3) and with the hypothesis (2.1), let λ ∈ Sp(L) such that m − n/2 < Re λ. With the notation of Proposition 3.1 for all ε > 0, we have
(3.4) ησ
νλ∈ H
−γ(λ)+εm(Ω) and ησ
νλ6∈ H
−γ(λ)m(Ω) .
Then there exists Z
νλwhich satisfies the homogeneous Dirichlet conditions and such that
P (ησ
λν− Z
νλ) ∈ C
0∞(Ω \ 0) and ∀ε > 0, Z
νλ∈ H
−γ(λ)−δ+εm(Ω) . Set S
νλ:= η σ
νλ− Z
νλand F
νλ:= P S
νλ. The F
νλfor λ, ν satisfying
(3.5) λ ∈ Sp(L), m − n/2 < Re λ < s − β + m − n/2 and ν = 1, . . . , N
λform a basis of (H
βs−m(Ω) ∩ Rg P
0)/ Rg P
βs.
P r o o f. First step. As a consequence of Proposition 3.1, for any u ∈ H
om(Ω) such that P u ∈ H
β−s−m(Ω), the following equivalence holds:
u ∈ H
β−sm(Ω) ⇔ ∀λ, ν as in (3.5), hP u, K
νλi = 0 .
Due to a classical regularity result for corner problems (see for instance the statement given in [2], p. 33), if u ∈ H
oβ−sm(Ω) satisfies P u ∈ H
βs−m(Ω), then u ∈ H
βs+m(Ω). Thus we only have to consider the above equivalence.
Since the K
νλare functions as well as all elements of the kernel of any operator P
τ∗, there exist e F
νλ∈ C
0∞(Ω \ 0) ∩ Rg P
0such that
∀λ, ν and λ
0, ν
0as in (3.5), h e F
νλ, K
νλ00i = δ
λ,λ0δ
ν,ν0.
Let e S
νλ∈ H
om(Ω) be such that P e S
νλ= e F
νλ. We now have to construct the S
νλsatisfying the assertions of Proposition 3.4 as linear combinations of the e S
νλ.
Second step. We use again the γ
jsatisfying (3.3) we have introduced in the previous proof. Applying Lemma 3.3 for τ = −γ
1and τ
0= −γ
0, we deduce that the e S
νλfor m − n/2 + γ
0< Re λ < m − n/2 + γ
1generate
{η σ
νλ| m − n/2 + γ
0< Re λ < m − n/2 + γ
1}
modulo H
o−γm1(Ω). Thus, for any such λ, there exist Z
νλ∈ H
o−γm1(Ω) such that
ησ
νλ− Z
νλis a linear combination of the e S
νλ00. So, the functions S
νλare constructed
for m − n/2 + γ
0< Re λ < m − n/2 + γ
1.
For the next step, corresponding to the weights −γ
1and −γ
2, we use the same arguments where we replace the e S
νλfor m − n/2 + γ
1< Re λ < m − n/2 + γ
2by the functions
S e
νλ− X
m−n/2+γ0<Re λ0<m−n/2+γ1
d
λν00S
νλ00where, according to Lemma 3.3, the coefficients d
λν00are chosen such that all the above functions belong to H
−γm1(Ω).
Step by step, we reach γ
J= γ and our S
νλare independent and their number is n
γ, which is what we need.
Now it is not too difficult to deduce from the two previous propositions and from the Green formula the three following statements.
With the functions S
νλwe have just constructed, we have the extension of Theorem 2.1 to the case (C3).
Theorem 3.5. In the case (C3) and with the hypothesis (2.1), assume the same hypotheses about s and β as in Theorems 1.1 and 2.1. Assume that u ∈ H
om(Ω) is such that P u ∈ H
βs−m(Ω). Then there exist coefficients c
λνsuch that
u − X
λ∈Sp(L)
m−n/2<Re λ<s−β+m−n/2 Nλ
X
ν=1
c
λνS
νλ∈ H
βs+m(Ω) .
As a result of the previous constructions, we have some independent functions S e
νλsuch that
hP e S
νλ, K
νλ00i = δ
λ,λ0δ
ν,ν0,
and the singularities S
λνare a basis of the space generated by the e S
νλ. Thus, we can show
Lemma 3.6. Under the assumptions of Theorem 3.5, there exists a basis e K
νλof the space generated by the K
νλfor m − n/2 < Re λ < s − β + m − n/2 such that
hP S
νλ, e K
νλ00i = δ
λ,λ0δ
ν,ν0. The e K
νλhave the form
K e
νλ= X
Re λ0≥Re λ
d
λν00K
νλ00. With these new elements of the kernel of P
∗we have Theorem 3.7. Under the assumptions of Theorem 3.5
c
λν= R
Ω
P u e K
νλdx .
The above results have to be compared with the following statements of [7]:
Corollaries 3.1 and 3.2, Theorems 3.3 and 3.4. Our hypothesis (2.1) is more general
than the hypothesis of [7], which in our framework would correspond to P is one-to-one H
om(Ω) → H
−m(Ω) .
The paper [5] gives similar expressions for the coefficients of the singularities in a different framework.
In the case (C2), under the extra assumption
(3.6) ∀λ, λ
0∈ Sp(L) such that Re λ, Re λ
0∈ ]m − n/2, s − β + m − n/2[ , λ − λ
06∈ N \ 0 ,
the e K
νλand the K
νλcoincide and the formula for the coefficients is the same as in Theorem 2.3:
c
λν= R
Ω