1 Introduction
In this paper, we are interested in the following biharmonic equation:
(Dτ)
{
L
τ
u
:=
Δ
2
u
-
τ
Δ
u
=
f
(
u
)
in
Ω
,
u
=
∇
u
=
0
on
∂
Ω
,
where Ω is a smooth bounded domain in ℝN, f is a smooth function and τ≥0.
The difficulty to handle existence or nonexistence of solutions with polyharmonic operators is often due to the fact that several fundamental tools like the maximum principle, the Harnack type inequality or the unique continuation argument are much less available in general, comparing with the second-order elliptic equations. When Ω is a ball, under the Dirichlet boundary conditions, we already know by Boggio [3, 4] that the Green’s function for Δ2 is positive, hence the positivity preserving property holds true, that is, the solution of Δ2u=f is positive if f is positive. On the other hand, the lack of the maximum principle in the Dirichlet boundary case generated many research works, the readers can find many interesting discussions in the book [9] and the references therein.
Our interest in problem ((Dτ)) comes also from its relation with microelectromechanical systems (MEMS), modeling the electrostatic actuation for membranes deflecting on thin plates in the field of nanotechnology detection systems. The Dirichlet boundary conditions are also called clamped boundary conditions, where the actuator at the boundary is clamped, giving rise to zero vertical displacement and zero slope. Actually, small deformations of a membrane clamped at its boundary are governed by
B
Δ
2
u
-
T
Δ
u
=
f
in
Ω
,
u
=
∇
u
=
0
on
∂
Ω
,
where BΔ2u with B>0 accounts for bending and -TΔu with T>0 for stretching, we refer to [7, 13, 19] for more detailed discussions.
On the contrary, for the Navier boundary conditions, which are also called pinned boundary conditions (u=Δu=0 on ∂Ω), we can easily apply the maximum principle.
For problem (D0), consider f(t)=|t|p*-1t, with p*:=N+4N-4, the critical exponent for the Sobolev embedding H2(Ω)⊂Lp(Ω). Mitidieri [15] and Oswald [16] proved that no positive solution exists when Ω is a star-shaped or starlike domain, but the nonexistence of any nontrivial solution remains an open problem, even for general smooth bounded convex domains. Reichel generalized in [22] this nonexistence result for conformally contractible domains which are not always star-shaped. Gazzola, Grunau and Squassina [8] proved the existence of a nontrivial solution for some contractible domain with a special geometry. To conclude, similar to the classical Yamabe problem (see [1] and the references therein), both topology and geometry of the domain Ω play important roles for the existence of solutions to problem (D0).
In the supercritical case, i.e., when f(t)=|t|p-1t, where p>p*, results due to Van Der Vorst [25] and Pucci and Serrin [21] show that (D0) has no nontrivial solutions when Ω is star-shaped. In [11], Khenissy and Ye considered some questions about the effect of topology and geometry of Ω on the existence of solutions to Δ2u=f(u) under the Dirichlet or Navier boundary conditions. In particular, they proved the following result.
Theorem 1.1
Theorem 1.1 ([11])
Let N≥7 and 1≤k≤N-6. Then there exist smooth bounded domains Ω⊂ℝN, homotopic to 𝕊k, such that for any continuous function f verifying tf(t)≥pF(t) in ℝ, with
p
>
2
(
N
-
k
-
1
)
N
-
k
-
5
𝑎𝑛𝑑
F
(
t
)
=
∫
o
t
f
(
s
)
𝑑
s
,
the only solution of (D0) is u≡0.
For the Laplacian case, McGough and Mortensen [14] introduced the notion of h-starlikeness of a domain Ω, and Schaaf [24] introduced a generalized critical exponent through some classifying number M(Ω). In [10], Khenissy adapted these definitions to the bilaplacian case under the Navier boundary conditions and proved some nonexistence results for h-starlike domains Ω under suitable growth condition on f.
In this paper, we are again interested in the effect of topology and geometry of the domain Ω for the existence of solutions to problem ((Dτ)). We will define the corresponding class of h-starlike domains for which equations ((Dτ)) do not have nontrivial solutions, under suitable growth assumption on f. The main tool is a Pohozaev identity computed with appropriate vector fields h.
In Section 2, we define a classifying number M*(Ω) for h-starlike domains with respect to the Dirichlet boundary conditions. This will allow us to get nonexistence results for ((Dτ)). The nonexistence theorem we obtain for (D0) generalizes Theorem 1.1, and gives the analogues to the results in [10] under the Navier boundary conditions. Our results show again that both geometrical and topological properties affect the existence of nontrivial solutions. In Section 3, some examples of such domains with rich geometry or topology are studied. Section 4 is devoted to prove a uniqueness result, using again the Pohozaev identity. Precisely, we consider the equation Lτu=λf(u) with the Dirichlet boundary conditions. This generalizes, for example, the result of Esposito and Ghoussoub for the second-order MEMS modelling (see [6]) and gives an analogous result to the Navier boundary one in [10].
2 Generalized Critical Exponent Under the Dirichlet Boundary Conditions
First, we show a Pohozaev identity for solutions to Lτu=f(u,x). We multiply the equation by h⋅∇u+au where h:Ω¯→ℝN is a smooth vector field, a∈ℝ, and integrate over Ω.
We will use the notation A⋅B:=∑aijbij for any matrices A=(aij), B=(bij)∈Mn(ℝ). Define also
F
(
s
,
x
)
:=
∫
0
s
f
(
t
,
x
)
𝑑
t
,
F
x
(
s
,
x
)
:=
(
∂
F
∂
x
1
(
s
,
x
)
,
…
,
∂
F
∂
x
N
(
s
,
x
)
)
=
∇
x
F
(
s
,
x
)
.
We have then the following result.
Lemma 2.1
Assume that Lτu=f(u,x) in Ω and u=∇u=0 on ∂Ω. For any smooth vector field h=(hi) from Ω¯ into ℝN and a∈ℝ,
∫
Ω
[
F
(
u
,
x
)
div
(
h
)
-
a
u
f
(
u
,
x
)
+
F
x
(
u
,
x
)
⋅
h
]
𝑑
x
=
∫
Ω
[
div
(
h
)
2
-
a
]
(
Δ
u
)
2
𝑑
x
-
2
∫
Ω
Δ
u
(
D
h
⋅
∇
2
u
)
𝑑
x
-
∫
Ω
Δ
u
(
Δ
h
⋅
∇
u
)
𝑑
x
+
∫
∂
Ω
Δ
u
∂
(
h
⋅
∇
u
)
∂
n
𝑑
σ
-
∫
∂
Ω
(
Δ
u
)
2
2
(
h
⋅
n
)
𝑑
σ
(2.1)
+
τ
∫
Ω
{
[
div
(
h
)
2
-
a
]
|
∇
u
|
2
-
∇
u
⋅
D
h
∇
u
}
𝑑
x
.
Here n denotes the unit external normal vector on ∂Ω and Dh is the differential map of h.
In order to have more precise information, we choose the vector field h satisfying
(2.2)
Δ
h
≡
0
ℝ
N
,
div
(
h
)
≡
1
and
h
⋅
n
≥
0
on
∂
Ω
.
Moreover, thanks to the Dirichlet boundary conditions, we know that (see for example [21])
∂
2
u
∂
x
i
∂
x
j
=
∂
2
u
∂
n
2
n
i
n
j
on
∂
Ω
,
hence
∫
∂
Ω
Δ
u
∂
(
h
⋅
∇
u
)
∂
n
𝑑
σ
=
∫
∂
Ω
Δ
u
∑
i
,
j
=
1
N
h
i
∂
2
u
∂
x
i
∂
x
j
n
j
d
σ
=
∫
∂
Ω
Δ
u
∂
2
u
∂
n
2
(
h
⋅
n
)
𝑑
σ
=
∫
∂
Ω
(
Δ
u
)
2
(
h
⋅
n
)
𝑑
σ
≥
0
,
where ni denotes the ith component of the normal vector n. Finally, with (2.2), we get
∫
Ω
[
F
(
u
,
x
)
-
a
u
f
(
u
,
x
)
+
F
x
(
u
,
x
)
⋅
h
]
𝑑
x
(2.3)
≥
1
-
2
a
2
∫
Ω
(
Δ
u
)
2
𝑑
x
-
2
∫
Ω
Δ
u
(
D
h
⋅
∇
2
u
)
𝑑
x
+
τ
∫
Ω
[
1
-
2
a
2
|
∇
u
|
2
-
∇
u
⋅
D
h
∇
u
]
𝑑
x
.
Given a bounded smooth domain Ω in ℝN, consider the set
ℋ
(
Ω
)
:=
{
h
∈
C
2
(
Ω
¯
,
ℝ
N
)
:
h
satisfies (2.2)
}
.
When ℋ(Ω)≠∅, we call the domain Ω to be a Navier h-starlike domain as in [10]. In order to define h-starlike domains with respect to the Dirichlet boundary conditions, we need more considerations on the vector field h. As ∥Δu∥2 and ∥u∥H02 are equivalent norms in H02(Ω), there exists Ch>0 such that
(2.4)
∫
Ω
Δ
v
(
D
h
⋅
∇
2
v
)
𝑑
x
≤
C
h
∫
Ω
(
Δ
v
)
2
𝑑
x
for all
v
∈
H
0
2
(
Ω
)
.
We denote by μ¯(h) the best constants in (2.4), that is,
μ
¯
(
h
)
:=
min
{
C
h
:
C
h
satisfies (2.4)
}
.
First, consider the case τ=0. If u is a solution of (D0), inequality (2.3) becomes
(2.5)
∫
Ω
[
F
(
u
)
-
a
u
f
(
u
)
]
𝑑
x
≥
∫
Ω
(
1
-
4
μ
¯
(
h
)
2
-
a
)
(
Δ
u
)
2
𝑑
x
.
In this case, we define just ℋ*(Ω)=ℋ(Ω) and
(2.6)
M
*
(
Ω
)
:=
sup
h
∈
ℋ
*
(
Ω
)
1
-
4
μ
¯
(
h
)
2
.
For τ>0, we choose h∈ℋ(Ω) satisfying furthermore
(2.7)
1
2
|
y
|
2
-
y
⋅
D
h
(
x
)
y
≥
c
h
|
y
|
2
for all
(
x
,
y
)
∈
Ω
×
ℝ
N
,
with ch a positive constant. Condition (2.7) was introduced first by McGough and Mortensen in [14] to define h-starlikeness for the second-order Dirichlet problem. Let, for τ>0,
ℋ
*
(
Ω
)
:=
{
h
∈
ℋ
(
Ω
)
:
h
satisfies (2.7)
}
.
Definition 2.2
A bounded smooth domain Ω is said to be Dirichlet h-starlike if ℋ*(Ω)≠∅.
For any τ≥0, we see that ℋ*(Ω)⊂ℋ(Ω), hence a Dirichlet h-starlike domain is always a Navier h-starlike domain. We denote by μ(h) the best constants in (2.7), that is, for h∈ℋ*(Ω),
μ
(
h
)
:=
max
{
c
h
:
c
h
satisfies (2.7)
}
.
Finally, for τ>0, we define the classifying number of a Dirichlet h-starlike domain Ω by
(2.8)
M
*
(
Ω
)
:=
sup
h
∈
ℋ
*
(
Ω
)
min
(
1
-
4
μ
¯
(
h
)
2
,
μ
(
h
)
)
.
If u is a solution of ((Dτ)), then
(2.9)
∫
Ω
[
F
(
u
)
-
a
u
f
(
u
)
]
𝑑
x
≥
∫
Ω
(
1
-
4
μ
¯
(
h
)
2
-
a
)
(
Δ
u
)
2
𝑑
x
+
τ
∫
Ω
(
μ
(
h
)
-
a
)
|
∇
u
|
2
𝑑
x
.
Assuming now M*(Ω)>0, 0<a<M*(Ω) and the growth condition F(s)≤asf(s) in ℝ, we easily conclude from (2.5) or (2.9) that no nontrivial solution to problem ((Dτ)) can exist.
Theorem 2.3
Let τ>0 (resp. τ=0) and Ω be a Dirichlet h-starlike domain in ℝN satisfying that M*(Ω) given by (2.8) (resp. (2.6)) is positive. Assume that f=f(u) verifies F(s)≤asf(s) in ℝ with 0<a<M*(Ω). Then there exists no nontrivial solution to problem ((Dτ)).
Remark 2.4
Consider a star-shaped domain Ω with respect to the origin up to a translation. Let h(x)=xN. Clearly h verifies (2.2), μ(h)=12-1N and μ¯(h)=1N. Therefore M*(Ω)>0 for N≥5 and we get some corresponding nonexistence results for ((Dτ)). In particular, M*(Ω)≥N-42N for (D0), hence if f(u)=|u|p-1u, no nonzero solution could exist to (D0) on a star-shaped domain Ω⊂ℝN with N≥5, if p>N+4N-4, which is just the classical result in [25, 21].
3 Examples
We exhibit here some examples of Dirichlet h-starlike domains with M*(Ω)>0. The first example gives some topologically nontrivial domains, it was constructed by Passaseo in [17, 18] to show nonexistence results for the Laplacian case, and adapted in [10, 11] for the biharmonic problem. In the second example, we revisit a domain with rich geometry mentioned in [24, 10].
Example 3.1
Let N≥3. We consider domains Ω in ℝN homotopic to spheres 𝕊k. Precisely, let η>0, 1≤k<N and consider the domain
Ω
=
𝕋
k
(
η
)
:=
{
x
∈
ℝ
N
:
(
∥
P
k
(
x
)
∥
-
η
-
1
)
2
+
∥
Q
k
(
x
)
∥
2
≤
η
2
}
,
where, for x=(xi)∈ℝN,
P
k
(
x
)
=
(
x
1
,
…
,
x
k
+
1
,
0
,
…
,
0
)
∈
ℝ
N
,
Q
k
(
x
)
=
x
-
P
k
(
x
)
.
Set Bk={x∈ℝN:Pk(x)=0} and define the vector field h as follows:
h
(
x
)
=
1
N
-
k
[
φ
(
r
k
)
P
k
(
x
)
+
Q
k
(
x
)
]
for any
x
∈
ℝ
N
∖
B
k
,
with rk(x)=∥Pk(x)∥ where ∥⋅∥ denotes the Euclidean norm and φ:(0,+∞)→ℝ is the function
φ
(
r
)
=
1
k
+
1
(
1
-
1
r
k
+
1
)
.
We have Ω=𝕋k(η)⊂ℝℕ∖Bk. Clearly, Ω is homotopic to 𝕊k, since
𝕋
k
(
η
)
=
{
x
∈
ℝ
N
:
dist
(
x
,
S
k
(
η
)
)
≤
η
}
,
where Sk(η)={x∈ℝN:∥x∥=η+1,Qk(x)=0}. Similarly to [11], we see that h verifies (2.2). Let us prove that h satisfies also (2.7). Let y∈ℝN and x∈Ω. The definition of h and a direct calculus give
1
2
|
y
|
2
-
y
⋅
D
h
y
=
(
1
2
-
φ
(
r
k
)
N
-
k
)
∑
i
=
1
k
+
1
y
i
2
+
(
1
2
-
1
N
-
k
)
∑
i
=
k
+
2
N
y
i
2
-
1
(
N
-
k
)
r
k
k
+
3
∑
i
,
j
=
1
k
+
1
x
i
x
j
y
i
y
j
.
Since
∑
i
,
j
=
1
k
+
1
x
i
x
j
y
i
y
j
≤
r
k
2
∑
i
=
1
k
+
1
y
i
2
,
we obtain
1
2
|
y
|
2
-
y
⋅
D
h
y
≥
(
1
2
-
φ
(
r
k
)
+
r
k
-
k
-
1
N
-
k
)
∑
i
=
1
k
+
1
y
i
2
+
(
1
2
-
1
N
-
k
)
∑
i
=
k
+
2
N
y
i
2
.
As 1≤rk≤2η+1 for any x∈Ω=𝕋k(η), we have
1
2
-
φ
(
r
k
)
+
r
k
-
k
-
1
N
-
k
≥
1
2
-
1
N
-
k
>
0
,
for all 1≤k<N-2. Then, we get (2.7) with ch=12-1N-k. We conclude that for k<N-2, the domain Ω is a Dirichlet h-starlike domain, non star-shaped and μ(h)≥N-k-22(N-k).
In order to estimate M*(Ω), we now focus on inequality (2.4). For any v∈H02(Ω),
D
h
⋅
∇
2
v
=
∑
i
,
j
=
1
N
∂
h
i
∂
x
j
∂
2
v
∂
x
i
∂
x
j
=
1
N
-
k
{
Δ
v
-
k
k
+
1
∑
i
=
1
k
+
1
∂
2
v
∂
x
i
2
}
+
∑
i
,
j
=
1
k
+
1
a
i
j
∂
2
v
∂
x
i
∂
x
j
with
a
i
j
=
1
N
-
k
(
x
i
x
j
r
k
k
+
3
-
δ
i
j
(
k
+
1
)
r
k
k
+
1
)
for all
1
≤
i
,
j
≤
k
+
1
.
Since rk=∥Pk(x)∥≥1 in Ω, there holds
∑
i
,
j
=
1
k
+
1
a
i
j
2
=
k
(
N
-
k
)
2
(
k
+
1
)
r
k
2
k
+
2
≤
k
(
N
-
k
)
2
(
k
+
1
)
in
Ω
.
By the Cauchy–Schwarz inequality, we have
∫
Ω
(
∑
i
,
j
=
1
k
+
1
a
i
j
∂
2
v
∂
x
i
∂
x
j
)
2
𝑑
x
≤
k
(
N
-
k
)
2
(
k
+
1
)
∫
Ω
∑
i
,
j
=
1
k
+
1
(
∂
2
v
∂
x
i
∂
x
j
)
2
d
x
:=
k
(
N
-
k
)
2
(
k
+
1
)
A
(
v
)
.
Therefore
(3.1)
∫
Ω
Δ
v
∑
i
,
j
=
1
k
+
1
a
i
j
∂
2
v
∂
x
i
∂
x
j
d
x
≤
∥
Δ
v
∥
2
×
∥
∑
i
,
j
=
1
k
+
1
a
i
j
∂
2
v
∂
x
i
∂
x
j
∥
2
≤
1
N
-
k
k
k
+
1
∥
Δ
v
∥
2
A
(
v
)
.
On the other hand, using ∇v=0 on ∂Ω and integrating by parts,
∫
Ω
∂
2
v
∂
x
i
2
∂
2
v
∂
x
j
2
𝑑
x
=
∫
Ω
(
∂
2
v
∂
x
i
∂
x
j
)
2
𝑑
x
≥
0
for all
1
≤
i
,
j
≤
N
.
Hence
(3.2)
B
(
v
)
:=
∫
Ω
Δ
v
∑
j
=
1
k
+
1
∂
2
v
∂
x
j
2
d
x
=
∫
Ω
∑
i
=
1
N
∑
j
=
1
k
+
1
(
∂
2
v
∂
x
i
∂
x
j
)
2
d
x
≥
A
(
v
)
.
Combining (3.1) and (3.2), we get
∫
Ω
Δ
v
(
D
h
⋅
∇
2
v
)
𝑑
x
≤
1
N
-
k
[
∫
Ω
(
Δ
v
)
2
𝑑
x
-
k
k
+
1
B
(
v
)
]
+
1
N
-
k
k
k
+
1
∥
Δ
v
∥
2
B
(
v
)
.
As, for any C>0 and D≥0,
max
ℝ
+
(
-
C
x
+
D
x
)
=
D
2
4
C
,
we obtain
∫
Ω
Δ
v
(
D
h
⋅
∇
2
v
)
𝑑
x
≤
5
4
(
N
-
k
)
∫
Ω
(
Δ
v
)
2
𝑑
x
for all
v
∈
H
0
2
(
Ω
)
.
So (2.4) holds with Ch=54(N-k), hence μ¯(h)≤54(N-k) for any k<N.
Finally, for k<N-2,
M
*
(
Ω
)
≥
min
(
1
-
4
μ
¯
(
h
)
2
,
μ
(
h
)
)
≥
min
(
1
2
-
5
2
(
N
-
k
)
,
N
-
k
-
2
2
(
N
-
k
)
)
=
N
-
k
-
5
2
(
N
-
k
)
.
For N≥7, we can take k≥1 such that M*(Ω)>0. Hence we can generalize [11, Theorem 1.1] as follows.
Theorem 3.2
Let τ>0. There exist smooth bounded domains in ℝN (N≥7), homotopic to 𝕊k with 1≤k≤N-6 which are Dirichlet h-starlike, satisfying M*(Ω)>0. Assume that f(u) verifies the growth condition F(s)≤asf(s) in ℝ with 0<a<M*(Ω). Then there exists no nontrivial solution to problem ((Dτ)).
Remark 3.3
More precisely, if f(t)=|t|p-1t and k≤N-6, then we get a nonexistence result on 𝕋k(η) for p>N-k+5N-k-5.
Example 3.4
Let b>0 and φ∈C2,α(-b,b) such that φ(t)>0 in (-b,b), φ(±b)=0 and φ<0 outside (-b,b). For x=(x1,x¯)∈ℝ×ℝN-1, we set Φ(x):=|x¯|2-φ(x1). Consider the domain
Ω
=
{
x
=
(
x
1
,
x
¯
)
∈
ℝ
×
ℝ
N
-
1
:
-
b
<
x
1
<
b
,
Φ
(
x
)
<
0
}
.
Moreover, suppose φ(t)=eξ(t) in (-b,b) with limt→±∞tξ′(t)=-∞ and assume that Ω is a regular domain. As ∇Φ(x)=(-φ′(x1),2x¯)≠0 on ∂Ω, the exterior normal to ∂Ω is given by
n
(
x
)
=
∇
Φ
|
∇
Φ
|
=
1
|
∇
Φ
|
(
-
φ
′
(
x
1
)
,
2
x
¯
)
.
For δ≥0, we consider the vector field
h
δ
(
x
)
=
(
δ
x
1
,
1
-
δ
N
-
1
x
¯
)
.
We have divhδ≡1 and Δhδ=0 in Ω for any δ. The condition hδ⋅n≥0 holds on ∂Ω if and only if
2
1
-
δ
N
-
1
φ
(
x
1
)
-
δ
x
1
φ
′
(
x
1
)
≥
0
for all
x
∈
∂
Ω
,
which is equivalent to
δ
≤
1
φ
¯
,
where
φ
¯
:=
max
-
b
<
t
<
b
(
1
+
N
-
1
2
t
φ
′
(
t
)
φ
(
t
)
)
.
Note that the maximum exists and there holds always 1φ¯≤1. Hence for 0≤δ≤1φ¯, we have hδ∈ℋ(Ω).
Now, for any x∈Ω and y=(y1,…,yN)∈ℝN, we have
1
2
|
y
|
2
-
y
⋅
D
h
δ
(
x
)
y
=
(
1
2
-
1
-
δ
N
-
1
)
|
y
|
2
+
1
-
N
δ
N
-
1
y
1
2
.
So, for N≥3, 0<δ≤δ*:=min(1φ¯,1N), we have
1
2
|
y
|
2
-
y
⋅
D
h
δ
(
x
)
y
≥
c
h
δ
|
y
|
2
with
c
h
:=
N
-
3
+
2
δ
2
(
N
-
1
)
>
0
.
Hence the domain Ω is Dirichlet h-starlike, i.e., ℋ*(Ω)≠∅ and μ(hδ)≥N-3+2δ2(N-1).
If v∈H02(Ω), then
∫
Ω
Δ
v
(
∇
h
δ
⋅
∇
2
v
)
𝑑
x
=
∫
Ω
Δ
v
[
δ
v
x
1
x
1
+
1
-
δ
N
-
1
(
Δ
v
-
v
x
1
x
1
)
]
𝑑
x
=
1
-
δ
N
-
1
∫
Ω
(
Δ
v
)
2
𝑑
x
+
(
δ
-
1
-
δ
N
-
1
)
∫
Ω
v
x
1
x
1
Δ
v
𝑑
x
.
As δ≤1N, it holds δ-1-δN-1≤0. And since ∇v=0 on ∂Ω, we have
∫
Ω
v
x
1
x
1
Δ
v
𝑑
x
=
∑
i
=
1
N
∫
Ω
v
x
1
x
i
2
𝑑
x
≥
0
.
Therefore, we are led to
∫
Ω
Δ
v
(
∇
h
δ
⋅
∇
2
v
)
𝑑
x
≤
1
-
δ
N
-
1
∫
Ω
(
Δ
v
)
2
𝑑
x
,
so μ¯(hδ)≤1-δN-1. Thus,
M
*
(
Ω
)
≥
sup
0
<
δ
≤
δ
*
min
(
N
-
5
+
4
δ
2
(
N
-
1
)
,
N
-
3
+
2
δ
2
(
N
-
1
)
)
=
N
-
5
+
4
δ
*
2
(
N
-
1
)
.
For N≥5, we have M*(Ω)>0. As φ¯ is an arbitrary function satisfying our assumptions, we can obtain domains with rich geometry as dumbbell shaped or multi-dumbbell shaped domains.
4 A Uniqueness Result in Dirichlet h-Starlike Domains
Our aim in this section is to study the following biharmonic equation with a bifurcation parameter
(Pλ)
{
Δ
2
u
-
τ
Δ
u
=
λ
f
(
u
)
in
Ω
,
u
=
∇
u
=
0
on
∂
Ω
,
where τ≥0, λ>0, and Ω⊂ℝN is a Dirichlet h-starlike domain with M*(Ω)>0, and f:ℝ→(0,∞) is assumed to be nondecreasing and C2 such that f(s)→+∞ as s→+∞ and
(4.1)
lim sup
s
→
+
∞
F
(
s
)
s
f
(
s
)
=
η
<
M
*
(
Ω
)
,
F
(
s
)
:=
∫
0
s
f
(
t
)
𝑑
t
.
Our approach is again based on the Pohozaev identity. We consider only classical solutions. We use ideas of Dolbeault and Stanczy [5] for the second-order problem with -Δ. The corresponding problem with τ=0 and Navier boundary conditions was considered in [10].
Let Λ1τ denote the first eigenvalue of the operator Lτ=Δ2-τΔ in H02(Ω), given by
Λ
1
τ
=
inf
u
∈
H
0
2
(
Ω
)
∖
{
0
}
∥
Δ
u
∥
2
2
+
τ
∥
∇
u
∥
2
2
∥
u
∥
2
2
>
0
and Φ1τ a corresponding eigenfunction. A result in [2] shows that if we assume that Φ1τ>0 in Ω, then the operator Lτ-1 is positivity preserving in H02(Ω), i.e.,
(4.2)
L
τ
u
=
f
,
u
∈
H
0
2
(
Ω
)
,
f
≥
0
⟹
u
≥
0
in
Ω
.
A very simple case of a Dirichlet h-starlike domain satisfying M*(Ω)>0 and Φ1τ>0 is a ball in ℝN, see [12].
Assume now that the positivity preserving property (4.2) holds true and that we have a nonnegative supersolution u¯ for ((Pλ)). Then we are able to use the monotone iterative scheme: let w0=0, and for n≥1,
Δ
2
w
n
-
τ
Δ
w
n
=
λ
f
(
w
n
-
1
)
in
Ω
,
w
n
=
∇
w
n
=
0
on
∂
Ω
to get a monotone sequence {wn} satisfying w0≤w1≤⋯≤wn≤u¯ by the monotonicity of f. In particular, the limit uλ=limwn is a classical solution of equation ((Pλ)) and uλ≤u¯. Therefore, uλ is the minimal solution of equation ((Pλ)), in other words, any solution v to ((Pλ)) verifies uλ≤v in Ω, since uλ does not depend on u¯ and we can take just u¯=v. As any solution of ((Pλ)) is a supersolution to (Pλ′) with λ′∈(0,λ), the set of λ≥0 such that ((Pλ)) admits a classical nonnegative solution is an interval I. Moreover, ξ=Lτ-1(1) is a supersolution to ((Pλ)) for λ>0 but small enough. We conclude then I≠∅. We can also see that the application λ↦uλ is nondecreasing and
(4.3)
lim
λ
→
0
∥
u
λ
∥
C
1
(
Ω
¯
)
=
0
.
Our main result is to claim that uλ is the unique solution for small λ>0.
Theorem 4.1
Let τ≥0. Suppose that the domain Ω is Dirichlet h-starlike with M*(Ω)>0 and Lτ satisfies (4.2). Assume that the nonlinearity f(u) satisfies (4.1). Then, there exists λ¯>0 such that the minimal solution uλ is the unique solution of ((Pλ)) for λ∈[0,λ¯].
Proof.
Assume that ((Pλ)) has two solutions, uλ and uλ+v, with v≥0, v≢0. Then v solves
(4.4)
{
Δ
2
v
-
τ
Δ
v
=
λ
[
f
(
v
+
u
λ
(
x
)
)
-
f
(
u
λ
(
x
)
)
]
in
Ω
,
v
=
∇
v
=
0
on
∂
Ω
.
Multiplying (4.4) by v, we obtain
(4.5)
τ
∫
Ω
|
∇
v
|
2
𝑑
x
+
∫
Ω
(
Δ
v
)
2
𝑑
x
=
λ
∫
Ω
v
[
f
(
v
+
u
λ
)
-
f
(
u
λ
)
]
𝑑
x
.
On the other hand, the Pohozaev identity (2.1) with a=0 and h∈ℋ*(Ω) yields
1
2
∫
Ω
(
Δ
v
)
2
𝑑
x
-
2
∫
Ω
Δ
v
(
D
h
⋅
∇
2
v
)
𝑑
x
+
1
2
∫
∂
Ω
(
Δ
v
)
2
(
h
⋅
n
)
𝑑
σ
+
τ
∫
Ω
[
1
2
|
∇
v
|
2
-
∇
v
⋅
D
h
∇
v
]
𝑑
x
(4.6)
=
λ
∫
Ω
[
F
(
v
+
u
λ
)
-
F
(
u
λ
)
-
f
(
u
λ
)
v
+
h
⋅
∇
u
λ
[
f
(
u
λ
+
v
)
-
f
(
u
λ
)
-
f
′
(
u
λ
)
v
]
]
𝑑
x
.
Recall that F is a primitive of f such that F(0)=0. Set
μ
*
(
h
)
=
min
(
μ
(
h
)
,
1
2
-
2
μ
¯
(
h
)
)
with μ(h), μ¯(h) defined as in Section 2, i.e. the best constants for (2.4), (2.7). By (4.1), we have M*(Ω)>η. We may choose a vector field h∈ℋ*(Ω) such that μ*(h)>η.
Using (4.3), there holds ∥h⋅∇uλ∥∞≤ε:=μ*(h)-μ, provided λ>0 is small enough. Define Hε by
H
ε
(
u
λ
,
v
)
:=
F
(
u
λ
+
v
)
-
F
(
u
λ
)
-
F
′
(
u
λ
)
v
+
ε
|
f
(
u
λ
+
v
)
-
f
(
u
λ
)
-
f
′
(
u
λ
)
v
|
-
η
1
v
[
f
(
u
λ
+
v
)
-
f
(
u
λ
)
]
,
with η1∈(η,μ*(h)). Because of the smoothness of f and (4.1), the function v-2Hε(uλ,v) is bounded from above by some constant Σ, uniformly in λ>0 small enough. From (4.5), (4.6), and the definitions of μ(h) and μ¯(h), it follows that any solution v of (4.4) satisfies
[
1
2
-
2
μ
¯
(
h
)
]
∫
Ω
(
Δ
v
)
2
𝑑
x
+
τ
μ
(
h
)
∫
Ω
|
∇
v
|
2
𝑑
x
≤
∫
Ω
1
2
(
Δ
v
)
2
𝑑
x
-
2
∫
Ω
Δ
v
(
∇
h
⋅
∇
2
v
)
𝑑
x
+
1
2
∫
∂
Ω
(
Δ
v
)
2
(
h
⋅
n
)
𝑑
σ
+
τ
∫
Ω
[
1
2
|
∇
v
|
2
-
∇
v
⋅
D
h
∇
v
]
𝑑
x
≤
λ
∫
Ω
{
H
ε
(
u
λ
,
v
)
+
η
1
v
[
f
(
u
λ
+
v
)
-
f
(
u
λ
)
]
}
𝑑
x
≤
λ
Σ
∫
Ω
v
2
𝑑
x
+
η
1
∫
Ω
[
(
Δ
v
)
2
+
τ
|
∇
v
|
2
]
𝑑
x
.
Applying Poincaré’s inequality, there exists C0>0 such that ∥v∥2≤C0∥∇v∥2, for any v∈H01(Ω). So we get
(4.7)
[
1
2
-
2
μ
¯
(
h
)
-
η
1
]
∫
Ω
(
Δ
v
)
2
𝑑
x
+
[
τ
μ
(
h
)
-
τ
η
1
-
C
0
2
λ
Σ
]
∫
Ω
|
∇
v
|
2
𝑑
x
≤
0
.
By the definition of μ*(h) and η1<μ*(h), for λ>0 sufficiently small, the left-hand side of (4.7) is positive as v≢0, which is absurd. To conclude, for λ>0 small enough, uλ is the unique solution of ((Pλ)). ∎
In the MEMS modelling, we consider Lτu=λf(u) with u∈H02(Ω) and f(u)=(1-u)-p with p>0. There holds
lim
s
→
1
F
(
s
)
s
f
(
s
)
=
0
,
F
(
s
)
:=
∫
0
s
f
(
t
)
𝑑
t
.
We can carry on exactly the same study as above to claim the uniqueness of solution for λ>0 small enough.
Part of this work was carried out while the first and third author were visiting the Institut Elie Cartan de Lorraine in Metz, France. They wish to express their gratitude to IECL for its warm hospitality and to Professor Dong Ye for the invitation.
References
[1] Bahri A. and Coron J. M., On a nonlinear elliptic equation involving the critical Sobolev exponent: The effect of the topology of the domain, Comm. Pure Appl. Math. 41 (1988), 253–294. 10.1002/cpa.3160410302Search in Google Scholar
[2] Ben Omrane H. and Khenissy S., A positivity preserving result for a biharmonic equation under Dirichlet boundary conditions, Opuscula Math. 34 (2014), no. 3, 601–608. 10.7494/OpMath.2014.34.3.601Search in Google Scholar
[3] Boggio T., Sull’equilibrio delle piastre elastiche incastrate, Rom. Acc. L. Rend. (5) 10 (1901), 197–205. Search in Google Scholar
[4] Boggio T., Sulle funzioni di Green dordine m, Rend. Circ. Mat. Palermo 20 (1905), 97–135. 10.1007/BF03014033Search in Google Scholar
[5] Dolbeault J. and Stańczy R., Non-existence and uniqueness results for supercritical semilinear elliptic equations, Ann. Henri Poincaré 10 (2010), 1311–1333. 10.1007/s00023-009-0016-9Search in Google Scholar
[6] Esposito P. and Ghoussoub G., Uniqueness of solutions for an elliptic equation modeling MEMS, Methods Appl. Anal. 15 (2008), 341–354. 10.4310/MAA.2008.v15.n3.a6Search in Google Scholar
[7] Esposito P., Ghoussoub G. and Guo Y., Mathematical Analysis of Partial Differential Equations Modeling Electrostatic MEMS, Courant Lect. Notes Math. 20, Courant Institute of Mathematical Sciences, New York, 2010. 10.1090/cln/020Search in Google Scholar
[8] Gazzola F., Grunau H.-C. and Squassina M., Existence and nonexistence results for critical growth biharmonic elliptic equations, Calc. Var. Partial Differential Equations 18 (2003), 117–143. 10.1007/s00526-002-0182-9Search in Google Scholar
[9] Gazzola F., Grunau H.-C. and Sweers G., Polyharmonic Boundary Value Problems, Positivity Preserving and Nonlinear Higher Order Elliptic Equations in Bounded Domains, Lecture Notes in Math. 1991, Springer, Berlin, 2010. 10.1007/978-3-642-12245-3Search in Google Scholar
[10] Khenissy S., Nonexistence and uniqueness for biharmonic problems with supercritical growth and domain geometry, Differential Integral Equations 24 (2011), 1093–1106. 10.57262/die/1356012878Search in Google Scholar
[11] Khenissy S. and Ye D., Nonexistence results for biharmonic boundary value problems with supercritical growth, Commun. Contemp. Math. 10 (2008), 195–204. 10.1142/S0219199708002764Search in Google Scholar
[12] Laurencot P. and Walker C., Sign-preserving property for some fourth-order elliptic operators in one dimension and radial symmetry, preprint 2013, http://arxiv.org/abs/1303.2237. 10.1007/s11854-015-0024-2Search in Google Scholar
[13] Lin F. H. and Yang Y. S., Nonlinear non-local elliptic equation modelling electrostatic actuation, Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 463 (2007), 1323–1337. 10.1098/rspa.2007.1816Search in Google Scholar
[14] McGough J. and Mortensen J., Pohozaev obstructions on non-starlike domains, Calc. Var. Partial Differential Equations 18 (2003), 189–205. 10.1007/s00526-002-0188-3Search in Google Scholar
[15] Mitidieri E., A Rellich type identity and applications, Comm. Partial Differential Equations 18 (1993), 125–151. 10.1080/03605309308820923Search in Google Scholar
[16] Oswald P., On a priori estimates for positive solutions of a semilinear biharmonic equation in a ball, Comment. Math. Univ. Carolin. 26 (1985), 565–577. Search in Google Scholar
[17] Passaseo D., Nonexistence results for elliptic problems with supercritical nonlinearity in non-trivial domains, J. Funct. Anal. 114 (1993), 97–105. 10.1006/jfan.1993.1064Search in Google Scholar
[18] Passaseo D., New nonexistence results for elliptic problems with supercritical nonlinearity, Differential Integral Equations 8 (1995), 577–586. 10.57262/die/1369316507Search in Google Scholar
[19] Pelesko J. A. and Bernstein D. H., Modeling MEMS and NEMS, Chapman & Hall/CRC, Boca Raton, 2003. 10.1201/9781420035292Search in Google Scholar
[20] Pohozaev S. I., Eigenfunctions of the equation Δu+λf(u)=0, Soviet, Math. Dokl 6 (1986), 1408–1411. Search in Google Scholar
[21] Pucci P. and Serrin J., A general variational identity, Indiana Univ. Math. J. 35 (1986), 681–703. 10.1512/iumj.1986.35.35036Search in Google Scholar
[22] Reichel W., Uniqueness results for semilinear polyharmonic boundary value problems on conformally contractible domains. II, J. Math. Anal. Appl. 287 (2003), 75–89. 10.1016/S0022-247X(03)00492-XSearch in Google Scholar
[23] Rellich F., Darstellung der Eigenwerte von Δu+λu=0 durch ein Randintegral, Math. Z. 46 (1940), 635–636. 10.1007/BF01181459Search in Google Scholar
[24] Schaaf R., Uniqueness for semilinear elliptic problems: Supercritical growth an domain geometry, Adv. Differential Equations 5 (2000), 1201–1220. 10.57262/ade/1356651221Search in Google Scholar
[25] Van Der Vorst R. C. A. M., Fourth order elliptic equations with critical growth, C. R. Math. Acad. Sci. Paris 320 (1995), 295–299. Search in Google Scholar