Kirchhoff-Carrier Strings in Noncylindrical Domains
Kirchhoff-Carrier Strings in Noncylindrical Domains
Abstract: The existence and uniqueness of local and global solutions for the
Kirchhoff–Carrier nonlinear model for the vibrations of elastic strings in noncylindrical
domains are investigated by means of the Galerkin Method. The asymptotic behaviour
of the energy is also studied.
1 – Introduction
b of Q
The lateral boundary Σ b is given by
[
b =
Σ {α(t), β(t)} × {t} .
0<t<T
In this work we investigate the existence, uniqueness and the asymptotic be-
haviour of solutions of the following noncylindrical mixed problem, for Kirchhoff–
Carrier elastic strings:
Received : January 12, 1998; Revised : February 18, 1998.
AMS Subject Classification: 35F30, 35K35.
Keywords: Kirchhoff–Carrier, Strings, Noncylindrical, Global solutions, Asymptotic
behaviour.
466 J. LIMACO FERREL and L.A. MEDEIROS
µZ β(t) µ ¶ ¶ 2
∂2u ∂u 2 ∂ u b
−M dx = f (x, t) in Q
∂t 2 ∂x ∂x2
α(t)
(1.1) b
u = 0 on Σ
∂u
u(x, 0) = u0 (x), (x, 0) = u1 (x) in α(0) < x < β(0) .
∂t
In the bibliography at the end of this article a complete list of papers dealing
with the Kirchhoff–Carrier operator
µZ ¶
∂2w 2
Lw = −M |∇w| dx ∆w
∂t2 Ω
We consider real functions α(t), β(t) and M (λ) satisfying the following con-
ditions:
(H1) α, β ∈ C 2 ([0, +∞]; R) with α(t) < β(t), α0 (t) ≤ 0, β 0 (t) ≥ 0 and
n o µ ¶1/2
0 0 m0
max |α (t)|, |β (t)| ≤
2
for all 0 ≤ t < ∞.
(H2) M ∈ C 1 ([0, ∞[; R) such that M (λ) ≥ m0 > 0 for all λ ≥ 0.
Remark 2.1. Note that the assumption α0 (t) ≤ 0 and β 0 (t) ≥ 0 means that
b is increasing in the sense that γ(t) = β(t) − α(t) is increasing.
Q
x−α
b the point (y, t), with y =
Observe that when (x, t) varies in Q varies in the
γ
b → Q, given by τ : (x, t) → (y, t) is
cylinder Q = ]0, 1] × ]0, T [ . The mapping τ : Q
a diffeomorphism. We transform system (1.1) by means of the change of variables
x−α
u(x, t) = v(y, t) with y = ,
γ
Here we have
• dx = γ dy ,
m0 m0
• M̌ (λ) = M (λ) − ≥ >0,
2 2
µ ¶2
m0 α0 + γ 0 y
• a(y, t) = − >0,
2 γ2 γ
µ ¶
α0 + γ 0 y
• b(y, t) = −2 ,
γ
µ ¶
α00 + γ 00 y
• c(y, t) = − .
γ
Theorem 2.1. Let Ωt be the interval ]α(t), β(t)[ , 0 < t < T , and suppose
that (H1) and (H2) hold. Then, given
u ∈ L∞ (0, T0 ; H01 (Ωt )∩H 2 (Ωt )), u0 ∈ L∞ (0, T0 ; H01 (Ωt )), u00 ∈ L∞ (0, T0 ; L2 (Ωt )) ,
b0.
solution of (1.1) in Q
v0 ∈ H01 (0, 1) ∩ H 2 (0, 1) , v1 ∈ H01 (0, 1) and g ∈ L∞ ([0, T ]; H01 (0, 1)) ,
v : Q0 → R
giving a positive symmetric bilinear continuous form a(t, v, w), which can be
handled easily when getting a priori estimates. Thus the operator Ľv is well
adapted to apply the energy method.
Remark 2.5. When getting estimates for v 0 in H01 (0, 1) and v in H 2 (0, 1)
terms of the form · ¸ · ¸
β 0 ∂ v 0 (1, t) 2 (−α0 ) ∂ v 0 (0, t) 2
+
γ ∂y γ ∂y
appear. In order to guarantee its positivity we need that α0 ≤ 0 and β 0 ≥ 0 or, in
b is increasing.
other words, that Q
as m → ∞.
0 1 d 1
(3.4) a(t, vm , vm )= a(t, vm , vm ) − a0 (t, vm , vm ) ,
2 dt 2
where Z 1 ∂v ∂w
a0 (t, v, w) = a0 (y, t) dy .
0 ∂y ∂y
We also have:
µ ¶ ¯ Z
0
∂vm 1 ¯y=1 1 1 ∂b
b(y, t) 0
, vm 0
= b(y, t) vm (y, t)2 ¯¯ − 0
(vm (t))2 dy
∂y 2 y=0 2 0 ∂y
(3.5) Z 1
1 ∂b 0
=− (v (t))2 dy .
2 0 ∂y m
472 J. LIMACO FERREL and L.A. MEDEIROS
∂ 2 vm
0
Estimate II. In the approximate system (3.1) we take v = − . This
∂y 2
gives:
µ ¶ ¯ ¯
1 d 0 1 1 d ¯¯ ∂ 2 vm ¯¯2
(3.12) kvm (t)k2 + M̌ kv (t)k 2
+
dt ¯ ∂y 2 ¯
m
2 dt 2 γ2 γ
µ 0 ¶ µ 0 ¶ µ 0 ¶
∂ 2 vm ∂vm0 ∂ 2 vm ∂vm ∂ 2 vm
+ a t, vm , − + b(y, t) , − + c(y, t) , − =
∂y 2 ∂y ∂y 2 ∂y ∂y 2
µ ¶
∂ 2 vm
= g, − .
∂y 2
KIRCHHOFF–CARRIER ELASTIC STRINGS 473
We have:
µ ¶ Z µ ¶
∂ 2 vm
0 1 d 1 ∂ 2 vm 2
a t, vm , − = a(y, t) dy
∂y 2 2 dt 0 ∂y 2
Z µ 2 ¶
1 1 0 ∂ vm 2
(3.13) − a (y, t) dy
2 0 ∂y 2
Z 1 ·
∂ ∂a ∂vm ∂vm
¸ 0 0 ¯¯y=1
∂a ∂vm ∂vm
− dy + ¯ .
0 ∂y ∂y ∂y ∂y ∂y ∂y ∂y ¯y=0
Note that:
µ ¶ Z
∂v 0 ∂ 2 vm
0 0 ∂ 2v0
∂vm 1
b(y, t) m , − =− m
b(y, t)
dy
∂y ∂y 2 0 ∂y ∂y 2
Z 1 µ 0 ¶2
1 ∂ ∂vm
= − b(y, t) dy .
0 2 ∂y ∂y
α0 + γ 0 y
Remark 3.1. Since b(y, t) = −2 ,
γ
Z µ ¶2 Z µ ¶2
1 d 1 ∂ 2 vm 1 1 ∂ 2 vm
+ a(y, t) dy − a0 (y, t) dy
2 dt 0 ∂y 2 2 0 ∂y 2
Z 1 · 0 ¸
0 ¯y=1 ¯
∂ ∂a ∂vm ∂vm ∂a ∂vm ∂vm ¯
− dy +
0 ∂y ∂y ∂y ∂y ∂y ∂y ∂y ¯y=0
(3.16) ¶2 ¯y=1
1
Z 1 ∂b µ ∂v 0 ¶2 1 0
∂vm
µ
¯
m ¯
+ dy − b(y, t) ¯
2 0 ∂y ∂y 2 ∂y y=0
Z 1 · ¸ 0 0 ¯y=1 ¯
∂ ∂vm ∂vm ∂vm ∂vm ¯
+ c(y, t) dy − c(y, t) =
0 ∂y ∂y ∂y ∂y ∂y ¯y=0
µ 0 ¶
∂g ∂vm
= , .
∂y ∂y
µ ¶
1 1
Remark 3.2. Denoting by µ(t) = M̌ kvk2 , we have:
γ2 γ
µ ¶ µ ¶ µ ¶
2 γ0 1 2 1 γ0 1
0
µ (t) = − 3 M̌ kvk2 + 3 M̌ 0 kvk2 ((v, v 0 )) − 4 M̌ 0 kvk2 kvk2
γ γ γ γ γ γ
and, by Estimate I and the fact that M ∈ C 1 ([0, ∞[, R) we obtain, by the hypo-
thesis (H1) on α, β:
|µ0 (t)| ≤ C + kvm
0
k.
or ¯ ¯ ¯ ¯ ¯ ¯
¯ ∂vm ¯ ¯ ∂ 2 vm ¯ ¯ ∂ 2 vm ¯
¯ (0, t) ¯ ≤ ¯¯ ¯ + kv k ≤ C + ¯ ¯ ,
¯ ∂y ¯ ∂y 2 ¯L2 (0,1)
m ¯ ∂y 2 ¯ 2
R L (0,1)
∂vm
A similar estimate is true for (1, t).
∂y
Remark 3.4. We have:
µ ¶2
m0 α0 + γ 0 y ∂a α0 + γ 0 y γ 0
a(y, t) = − >0 and = −2 ,
2 γ2 γ ∂y γ γ
α00 + γ 00 y α0 + γ 0 y γ 0
c(y, t) = − − .
γ γ γ
1 β0 1 β 0 ³ ∂vm0 ´2
If = , then (1, t) goes to the left side of (3.17) and it is compen-
λ γ 4 γ ∂y
β 0 ³ ∂v 0 ´2
sated with (1, t) which gives a positive contribution in the first member.
γ ∂y
¯ β 0 γ 0 ¯ ¯ ∂v ¯2
¯ ¯¯ m ¯
The term ¯ 2 ¯ ¯ (1, t)¯ can be estimated as in Remark 3.3. We get
γ ∂y
¯ 0 0 ¯2 ¯ ¯2 ¯ 2 ¯
¯β γ ¯ ¯ ∂vm ¯ ¯ ∂ vm ¯ 2
¯ ¯ ¯ ¯ ¯ ¯
¯ γ2 ¯ ¯ ∂y (1, t)¯ ≤ C + ¯ ∂y ¯ ,
The same argument is true for all the other terms above.
476 J. LIMACO FERREL and L.A. MEDEIROS
Remark 3.5. From the hypothesis on α and β, we can estimate all the other
¯ ∂ 2 v ¯2
¯ m¯ 0 (t)k2 .
terms in the left side of (3.17) by C ¯ 2
¯ and Ckvm
∂y
Then by the Remarks above, we modify (3.17) obtaining:
¯ ¯ Z µ ¶2
d 0 2 d ¯ ∂ 2 vm ¯¯2 d 1 ∂ 2 vm
(3.18) kvm k + µ(t) ¯¯ + a(y, t) dy ≤
dt dt ∂y 2 ¯ dt 0 ∂y
¯ 2 ¯
¯ ∂ vm ¯ 2
0
≤ C0 + C1 kvm (t)k2 + ¯¯ ¯ .
∂y 2 ¯
· ¸
d ¯¯ ∂ 2 vm ¯¯2 d ¯ ∂2v ¯
¯ m¯
¯ 2 ¯2
0 (t)¯ ∂ vm ¯ and by
Substituting µ ¯ ¯ in (3.18) by µ(t) ¯ ¯ − µ ¯ ¯
dt ∂y 2 dt ∂y 2 ∂y 2
Remark 3.2, we obtain:
" ¯ ¯ Z µ ¶2 #
d ¯ ∂ 2 vm ¯ 2 1 ∂vm
(3.19) 0
kvm (t)k2 + µ(t) ¯¯ (t)¯¯ + a(y, t) dy ≤
dt ∂y 2 0 ∂y
¯ 2 ¯
¯ ∂ vm ¯ 2
0
≤ C0 + C1 kvm (t)k2 + kvm
0
(t)k ¯¯ ¯ in [0, T ] .
∂y 2 ¯
If ¯ 2 ¯2 Z µ ¶2
¯ ∂ vm ¯ 1 ∂vm
0
hm (t) = kvm (t)k2 + µ(t) ¯¯ (t)¯¯ + a(y, t) dy ,
∂y 2 0 ∂y
we have from (3.19):
dhm
≤ C0 + C1 hm + h3/2
m in [0, T ] .
dt
From this inequality, we get a number 0 < T0 < T , such that hm (t) is bounded
in [0, T0 ] independently of m. This gives the second estimate
¯ 2 ¯
¯ ∂ vm ¯ 2
(3.20) 0
kvm (t)k2 + ¯¯ ¯ <C on [0, T0 ] .
∂y 2 ¯
Note that µ(t) is strictly positive on [0, T ].
(3.21) 00
|vm (t)|2 < C on [0, T0 ] .
KIRCHHOFF–CARRIER ELASTIC STRINGS 477
Proof of Theorem 2.2: In this step we prove that the estimates obtained
above are sufficient to take limits in the approximate equation (3.1). In view of
(3.11), (3.20) and (3.21) a subsequence represented by (vk ) can be extracted from
(vm ) such that:
(4.1) vk * v weak star in L∞ (0, T0 ; H01 (0, 1) ∩ H 2 (0, 1)) ,
(4.2) vk0 * v 0 weak star in L∞ (0, T0 ; H01 (Ω)) .
By the classical compactness argument of Aubin–Lions, cf. Lions [13], it fol-
lows that:
(4.3) vk → v strongly L2 (0, T0 ; H01 (0, 1)) .
Because of the estimate (3.21) the subsequence satisfies also:
(4.4) vk00 * v 00 weak star in L∞ (0, T0 ; L2 (Ω)) .
From hypothesis (H1), (H2) and the estimates (3.11) and (3.20), we obtain:
¯ µ ¶2 2 ¯
¯1 1 2 ∂ v ¯¯
(4.5) ¯
¯ γ 2 M̌ γ kvm (t)k ∂y 2 ¯
<C on [0, T0 [ .
By (4.6) the first right side integral of (4.7) goes to zero when k → ∞ and the
second one, by (3.20), also goes to zero. To analyse the third member of the right
side of (4.7) we employ hypothesis (H2) on M (λ). Then, we have:
¯ ¯ ³ ´
¯ ¯
|µk (t) − µ(t)| ≤ c ¯ kvk (t)k2 − kv(t)k2 ¯ ≤ c kvk (t) − v(t)k kvk (t)k + kv(t)k .
It follows from (4.3), estimates (3.11), (3.20) and Lebesgue convergence the-
orem, that the last term of (4.7) goes to zero when k → ∞.
or, equivalently,
Uniqueness
If v and v̂ are two solutions in the conditions of Theorem 2.2, then w = v − v̂
satisfies:
µ 2 ¶ 2µ ¶ µ ¶
1 1 2 ∂ v 1 1 2 ∂ v̂ ∂ ∂w
(4.13) w 00 − 2
M̌ kvk 2
+ 2
M̌ kv̂k 2
− a(y, t) +
γ γ ∂y γ γ ∂y ∂y ∂y
∂w0 ∂w
+ b(y, t) + c(y, t) = 0 in L2 (0, T ; L2 (0, 1)) ,
∂y ∂y
w(0) = w 0 (0) = 0 in Ω and w = 0 on ]0, 1[ × ]0, T0 [ .
KIRCHHOFF–CARRIER ELASTIC STRINGS 479
We have:
µ ¶ µ µ ¶ ¶
1 1 d 1 d 1 1
2
M̌ kvk2 kwk2 = 2
M̌ kvk2 kwk2
2γ γ dt 2 dt γ γ
µ ¶
γ0 1
(4.15) + 3 M̌ kvk2
γ γ
µ ¶ µ ¶
1 1 d 1
− 2 M̌ 0 kvk2 kvk2 kwk2
2γ γ dt γ
Z Z
1 ∂w0 0 1 1 ∂b
(4.16) b(y, t) w dy = − (w0 )2 dy .
0 ∂y 2 0 ∂y
We also set
³ ´
(4.20) v0 (y) = u(x, 0) = u0 α(0) + γ(0) y ,
³ ´
v1 (y) = u0 (x, 0) = u1 α(0) + γ(0) y
(4.21) ³ ´ ³ ´
+ α0 (0) + γ 0 (0) y u00 α(0) + γ(0) y .
The function u(x, t) defined by (4.18) is the solution of Theorem 2.1. To see
this it is sufficient to observe that the application
µ ¶
x−α
(x, t) → ,t
γ
∂2u 1 ∂2v
(4.22) (x, t) = (y, t) ,
∂x2 γ 2 ∂y
µ ¶
00 00 ∂ ∂v ∂v 0 ∂v
(4.23) u (x, t) = v (y, t)− a(y, t) +b(y, t) (y, t)+c(y, t) (y, t) ,
∂y ∂y ∂y ∂y
x−α
with y = ,
γ
Z β(t) µ ∂u ¶2 1
Z 1µ
∂v
¶2
(4.24) dx = dy .
α(t) ∂x γ 0 ∂y
x−α
with y = . Then u solves (1.1), with initial conditions u0 and u1 .
γ
The regularity of v given by Theorem 2.2, implies the regularity of u claimed
in Theorem 2.1.
To prove the uniqueness observe that from (4.22)–(4.25) we have the equiv-
alence between the mixed problems (1.1) and (2.1). Then, if u and û are two
solutions of (1.1) given by Theorem 2.1, then v and v̂ obtained by (4.18) are
solutions in the conditions of Theorem 2.2. Since we have uniqueness for v, i.e.,
v = v̂ it implies u = û.
KIRCHHOFF–CARRIER ELASTIC STRINGS 481
5 – Applications
We will give examples of functions α(t) and β(t) in C 1 ([0, ∞[ , R), such that
α(t) < β(t), α0 (t) ≤ 0 and β 0 (t) ≥ 0. If γ(t) = β(t) − α(t), we have γ 0 (t) ≥ 0 and,
by hypothesis, (H1), Remark 2.2, we must have:
¯ ¯ µ ¶
¯ 0 0 ¯ m0 1/2
¯α (t) + γ (t) y ¯ ≤ ,
2
for all 0 ≤ t ≤ T and 0 ≤ y ≤ 1.
Let us consider the family of straight lines:
x = α0 (t) + γ 0 (t) y
depending of the parameter t ≥ 0.
t0 = (2k)
2k
1−2k
( m20 ) 1−2k .
6 – Global solutions
Then
γ(t) = 2β(t) for all t ≥ 0
and by Remark 2.2, we must have:
µ ¶1/2
0 m0
(6.1) 0 < β (t) < C .
2
with µ ¶1/2 µ ¶
1 2 π
C= .
6 5 π+1
We suppose, in the present section, that
x+β
(6.3) y= or x = (2 y − 1) β ,
2β
with 0 ≤ y ≤ 1.
We consider the perturbated system (6.4). The modified Kirchhoff–Carrier
model with damping is, for δ > 0 fixed, of the type:
µZ β(t) µ ¶ ¶ 2 µ 0 ¶
00 ∂u 2 ∂ u β ∂u 0 b ,
u −M dx +δ x +u = 0 in Q
−β(t) ∂x ∂x2 β ∂x
(6.4) b ,
u = 0 on Σ
u(x, 0) = u (x), u0 (x, 0) = u (x) on − β0 < x < β0 ,
0 1
where β0 = β(0).
KIRCHHOFF–CARRIER ELASTIC STRINGS 483
We have:
· ¸
1 β0 m0 ³ 0 ´2
− a0 (y, t) = 3
β 00 β (1 − 2 y)2 + − β (1 − 2 y) .
2 4β 2
1 m0 β 0
(6.6) − a0 (y, t) ≥ .
2 16 β 3
484 J. LIMACO FERREL and L.A. MEDEIROS
Note that
β 0 ∂u
v0 = x + u0
β ∂x
which gives in (6.5) the damping δv 0 .
Now we determine the condition on the initial data v0 , v1 which implies the
global existence of solution for the mixed problem (6.5).
Let us fix the number δ such that
µ ¶1/2 µ ¶
δ 5 m0 4 + π2
(6.7) > .
2 2 β0 2 π2
δ δ2
H(t) = kv 0 (t)k2 + (∇v 0 (t), ∇v(t)) + kv(t)k2
(6.8) 4 8
µ ¶ Z 1
1 1 2 2
+ M̌ kv(t)k |∆v(t)| + a(y, t) (∆v)2 dy ,
8 β2 2β 0
∂2v ∂v
where ∆v = 2
, ∇v = and v is a solution of the approximate problem for
∂y ∂y
(6.5).
Proof: In fact,
¯ ¯ ¯ ¯ 2
¯δ ¯
¯ (∇v, ∇v 0 )¯ ≤ δ ¯¯(∇v, ∇v 0 )¯¯ ≤ δ kvk2 + |v 0 |2 .
¯4 ¯ 4 32
m0 m0
We know that M̌ (λ) = + m1 λ ≥ , then
2 2
1 m0
2
M̌ (λ) ≥
8β 16 β 2
1 0 2 3 δ2 m0
H(t) ≥ kv k + kvk2 + |∆v|2 .
2 32 16 β 2
Theorem 6.1. Suppose (H1) and (H3) are true and α(t) = −β(t). Suppose
also that M = M (λ) is of the form (6.2). Then given
Remark 6.1. Before going into the proof of Theorem 6.1, we analyse the
restriction (6.10) on the initial data u0 , u1 of (6.4).
We have
³ ´
(i) v0 (y) = u0 (2 y − 1) β0 , with β0 = β(0) ,
³ ´ d ³ ´
(ii) v1 (y) = u1 (2 y − 1) β0 + (2 y − 1) β 0 (0) u0 (2 y − 1) β0 .
dx
Let us represent (−β0 , β0 ) by Ω0 . Then we have:
(iii) kv0 k2 = 2 β0 ku0 k2H 1 (Ω0 ) ,
m0 m0
(iv) a(t, v0 , v0 ) ≤ 2 kv0 k2 = ku0 k2H 1 (Ω0 ) ,
8 β0 4 β0
(v) |4 v0 |2 = 8 β02 |4 u0 |2L2 (Ω0 ) ,
1 (β 0 (0))2
(vi) |v1 |2 ≤ |u1 |2L2 (Ω0 ) + ku0 k2H 1 (Ω0 ) ,
β0 β0 0
(vii) kv1 k2 ≤ 4 β0 ku1 k2H 1 (Ω0 ) + 8 β0 β 02 (0) |∆u0 |2L2 (Ω0 ) + 8 β0 ku0 k2H 1 (Ω0 ) ,
Z 1
0 0
2
(viii) a(y, 0, (∆v0 ) ) dy ≤ β0 m0 |∆u0 |2L2 (Ω0 ) .
0
and
(x) H(0) ≤ G0 = 12 β0 ku1 k2H 1 (Ω0 ) + 13 β0 m0 |∆u0 |2L2 (Ω0 )
µ 0
¶
582
+ 24 β0 + β0 ku0 k2H 1 (Ω0 )
³ 16 ´ 0
Therefore, if we take
Note that this basis can be obtained explicitly as in Section 3. Let us obtain the
a priori estimates.
(β 0 )2
Because of (H3), i.e., |β 00 | < , we obtain:
β
¯Z 1 ¯ Z 1 ¯ 00 ¯
¯ ¯ ¯ β (1−2y) ∂v 0 ¯
¯ c(y, t) ∂v v 0 dy ¯ ≤ ¯ v ¯ dy
¯ ∂y ¯ ¯ 2β ∂y ¯
0 0
Z ¯ ¯
1 (β 0 )2 ¯ ∂v ¯
¯ ¯ · |v 0 | dy
≤ ¯ ∂y ¯
0 2 β2
Z ¯ ¯2 Z
1 µ(β 0 )4 ¯ ∂v ¯
¯ ¯ dy +
1 1 02
≤ ¯ ∂y ¯ |v | dy .
0 4 β4 0 4µ
µ(β 0 )4 m0 β 0 1 (β 0 )3
If we consider µ such that = , we obtain = . Then:
4 β4 16 β 3 4µ m0 β
¯Z 1 ¯ µ ¶
¯ ¯ 0
¯ c(y, t) ∂v v 0 dy ¯ ≤ m0 β kvk2 + 1 m0 1/2 0 2
¯ ¯ |v | .
0 ∂y 16 β 3 16 β0 2
Whence,
Z √
1 ∂v 0 m0 β 0 2 m0
(6.13) c(y, 1) v dy ≥ − 3
kvk − √ |v 0 |2 .
0 ∂y 16 β 16 2 β0
From (6.12) and (6.13), we have:
Z 1
1 0 ∂v 0
(6.14) − a (t, v, v) + c(y, t) v dy + δ |v 0 |2 ≥
2 0 ∂y µ µ ¶ ¶
1 m0 1/2
≥ δ− |v 0 |2 .
16 β0 2
By (6.14) we modify (6.11), obtaining:
µ µ ¶¶
1 d 0 2 d 1 c 1
(6.15) |v (t)| + M kv(t)k2 +
2 dt dt 4 β 2β
µ µ ¶ ¶
1 d 1 m0 1/2
+ a(t, v, v) + δ − |v 0 |2 ≤ 0 .
2 dt 8 β0 2
The parameter δ was fixed satisfying the condition (6.7), what implies:
µ ¶1/2
δ 1 m0
> .
2 8 β0 2
Then, from (6.15) we get:
· µ ¶¸
1 d 0 2 d 1 c 1 1 d δ
|v (t)| + M kv(t)k2 + a(t, v, v) + |v 0 (t)|2 ≤ 0 .
2 dt dt 4 β 2β 2 dt 2
488 J. LIMACO FERREL and L.A. MEDEIROS
for all t ≥ 0.
If we set
µ ¶
1 1 c 1
C1 = |v1 |2 + M kv0 k2 + a(0, v0 , v0 )
2 4 β0 2 β0
for all t > 0. Note that C1 does not depend on t, but depends on v0 , v1 , β0 ,
a(y, 0).
∂ 2v0
Estimate II. We multiply both sides of the equation (6.5)1 by −∆v 0 = −
∂y 2
and integrate on (0, 1).
We obtain:
µ ¶
1 d 0 1 1 d
(6.17) kv (t)k2 + 2
M̌ kv(t)k2 |∆v(t)|2 +
2 dt 8β 2β dt
Z 1 Z 1
1 d 1 2
+ a(y, t) (∆v) dy − a0 (y, t) (∆v)2 dy
2 dt 0 2 0
Z µ ¶2 µ ¶¯
1 1 ∂b ∂v 0 1 ∂v 0 ¯¯y=1
+ dy − b(y, t)
2 0 ∂y ∂y 2 ∂y ¯y=0
Z µ ¶ ¯
1 ∂ ∂a ∂v ∂v 0 ∂a ∂v ∂v 0 ¯¯y=1
− dy +
0 ∂y ∂y ∂y ∂y ∂y ∂y ∂y ¯y=0
Z · ¸ ¯
1 ∂ ∂v ∂v 0 ∂v ∂v 0 ¯¯y=1
+ c(y, t) dy − c(y, t) + δ kv 0 (t)k2 = 0 .
0 ∂y ∂y ∂y ∂y ∂y ¯y=0
The inequalities below for the various terms of this identity are true. In fact,
we have:
Z
1 1 m0 β 0
(6.18) − a0 (y, t) (∆v)2 dy ≥ |∆v|2 ,
2 0 32 β 3
KIRCHHOFF–CARRIER ELASTIC STRINGS 489
¯Z · ¸ ¯ µ ¶1/2 µ ¶
¯ 1 ∂ ∂a ∂v ∂v 0 ¯¯ m0 β 0 2 m0 4 + π2
(6.19) ¯
¯ dy ≤ |∆v|2 + kv 0 k2 ,
0 ∂y ∂y ∂y ∂y ¯ 64 β 3 β0 2 π2
¯ Z 1 ∂b µ ∂v 0 ¶2
¯ µ ¶
¯1
¯
¯
¯ 1 m0 1/2 0 2
(6.20) ¯2 dy ¯ ≤ kv k ,
0 ∂y ∂y 2 β0 2
¯Z 1 · ¸ ¯ µ ¶ µ ¶
¯ ∂ ∂v ∂v 0 ¯¯ m0 β 0 2 1 m0 1/2 4+π 2
(6.21) ¯
¯ ∂y c(y, t) ∂y ∂y dy ¯ ≤ 64 β 2 |∆v| + 2β kv 0 k2 ,
0 0 2 π2
µ ¶2 ¯y=1 µ ¶ µ ¶
∂v 0 ¯
¯ β 0 ∂ v 0 (1, t) 2 β 0 ∂ v 0 (0, t) 2
(6.22) −b(y, t) ¯ = + ,
∂y y=0 2β ∂y 2β ∂y
¯ ¯
¯ ∂a ∂v ∂v 0 ¯¯y=1 ¯ µ
β 0 ∂v 0
¶2 µ
β 0 ∂v 0
¶2
¯ ¯ ¯
¯ ¯ ≤ (1, t) + (0, t)
¯ ∂y ∂y ∂y ¯y=0 ¯ 8 β ∂y 8 β ∂y
(6.23) µ 0 ¶3 µ ¶2
β π+1
+4 |∆v|2 ,
β π
¯ ¯ ¯ µ ¶2 µ ¶2
¯ ∂v ∂v 0 ¯¯1 ¯¯ β 0 ∂v β 0 ∂v 0
¯
¯ c(y, t) ¯ ≤ (1, t) + (0, t)
¯ ∂y ∂y ¯0 ¯ 8 β ∂y 8 β ∂y
(6.24) µ 0 ¶3 µ ¶2
β π+1
+ |∆v|2 .
β π
Using (6.18)–(6.24) the following inequality can be obtained from (6.17):
µ ¶
1 d 0 1 1 d
(6.25) kv (t)k2 + M̌ kv(t)k2 |∆v(t)|2 +
2 dt 8 β2 2β dt
Z µ ¶3 µ ¶2
1 d 1
2 β0 π+1 m0 β 0
+ a(y, t) (∆v) dy − |∆v|2 + |∆v|2
2 dt 0 β π 32 β 2
µ ¶2 µ ¶2
β0
∂v 0 β 0 ∂v 0
+ (1, t) + (0, t)
4 β ∂y 4 β ∂y
· µ ¶1/2 µ ¶¸
5 m0 4+π 2
+ δ− kv 0 (t)k2 ≤ 0 .
2 β0 2 π2
Whence,
µ ¶
d 1 1
(6.27) (∇v 0 , ∇v) − kv 0 k2 + M̌ kvk2 |∆v|2 −
dt 4 β2 2β
µ ¶ µ ¶
δ d 2 ∂v 0 ∂v
− a(t, v, ∆v) − kvk + b , −∆v + c , −∆v = 0 .
2 dt ∂y ∂y
d δ d m0
(6.31) (∇v 0 , ∇v) + kvk2 + |∆v|2 − 2 kv 0 k2 ≤ 0 ,
dt 2 dt 8 β2
m0
since M̌ (λ) > > 0 and the coefficient of |∆v|2 is also positive. Note that δ is
2
δ
a fixed parameter. Multiply (6.31) by and adding to (6.26) we cancel −2 kv 0 k2 ,
4
KIRCHHOFF–CARRIER ELASTIC STRINGS 491
obtaining:
" µ ¶
1 d δ δ2 0 2 1 1
(6.32) kv 0 k2 + (∇v 0 , ∇v) + kv k + 2
M̌ kvk2 ) |∆v|2 +
2 dt 4 8 8β 2β
Z #
1 1 m0 δ
2
+ a(y, t) (∆v) dy + |∆v|2 ≤
2 0 32 β 2
· µ ¶¸0
1 1
≤ 2
M̌ kv(t)k2 |∆v|2 .
8β 2β
m0
Note that M̌ (λ) = + m1 λ with m0 , m1 > 0. Then M̌ 0 (λ) = m1 is constant.
2
We obtain:
· µ ¶¸0 µ ¶
1 1 β0 1 m1
M̌ kv(t)k2 = − 2 M̌ kvk2 + (∇v, ∇v 0 )
8 β2 2β 4β 2β 8 β3
β0
− kvk2 , with − β 0 = α0 < 0 ,
16 β 4
or
· µ ¶¸0
1 1 m1 kvk 0
(6.33) 2
M̌ kvk2 ≤ kv k .
8β 2β 8 β2 β
with µ ¶
1 1 c 1
C1 = |v1 |2 + M kv0 k2 a(0, v0 , v0 ) .
2 4 β0 2 β0
Then, from (6.31) and (6.32) it follows:
· µ ¶¸0 µ ¶1/2
1 1 m 1 C1
(6.35) 2
M̌ kvk2 2
|∆v| ≤ kv 0 k |∆v|2 .
8β 2β 2 β m0
where µ ¶1/2
C1 m1 0
γ(t) = kv (t)k
m0 2
p
and kv 0 (t)k < H(t).
m0 δ
Suppose that γ(t) ≥ and let us see that this leads us to a contradiction.
32
In fact, we have:
1 m1 q
γ(t) < √ C1 H(t) ,
m0 2
then
1 m1 q
(6.38) γ(0) ≤ √ C1 H(0) .
m0 2
m0 δ
γ(0) < .
32
n mδo
Let us consider t∗ = min t > 0, γ(t) = . This minimum exists and it is
32
greater than zero, by continuity of γ(t). We have:
m0
on 0 ≤ t < t∗ ,
γ(t) < 32 δ
(6.39)
m
γ(t∗ ) = 0 δ .
32
We obtain:
Z t∗
∗
(6.40) H(t ) − H(0) = H 0 (s) ds .
0
By (6.34) we have H 0 (s) < 0 on (0, t∗ ), then H(t∗ ) < H(0) or H(t∗ )1/2 < H(0)1/2 .
Consequently,
µ ¶ µ ¶
∗ 1 m1 q ∗
1 m1 q m0
γ(t ) < √ C1 H(t ) < √ C1 H(0) < δ
m0 2 m0 2 32
Then, from (6.37) we have H 0 (t) ≤ 0 for all t ≥ 0 and, therefore, the global
estimate
The estimate for v 00 (t) and the uniqueness can be proved as in the local case.
7 – Asymptotic behaviour
which we still call the energy associated to the solution v of the mixed problem
(6.5).
m0 c(λ) = m0 λ + m1 λ2 and (7.1)
Remember that M̌ (λ) = + m1 λ, then M
2 2 2
takes the form
1 0 2 m0 m1 1
(7.2) E(v, t) = |v (t)| + kv(t)k2 + kv(t)k4 + a(t, v(t), v(t)) .
2 16 β 2 32 β 3 2
Theorem 7.1. If v is the global solution of the mixed problem (6.5), then
we have R t ds
−C2
E(v, t) ≤ C1 e 0 β(s)2 ,
1 + 2 k 16 (k + δ)
C0 = +
β02 m0
δ
with k such that 20 k = .
4
494 J. LIMACO FERREL and L.A. MEDEIROS
where k is a constant to be fixed; see Remark 7.1. This method of working with a
perturbed energy F (v, t) was employed by Komornik–Zuazua [11]. We also refer
to Nakao–Narazaki [15] and Hosya–Yamada [9].
Since
¯ ¯ k 0 2
¯ ¯
(7.4) ¯2 k (v 0 (t), v(t))¯ ≤ |v (t)| + k δ |v(t)|2 ,
δ
we obtain:
k
(7.5) 2 k (v 0 (t), v(t)) + k δ |v(t)|2 ≥ − |v 0 (t)|2 .
δ
Then, from (7.4) and (7.5) we modify (7.3), obtaining:
µ ¶
1 k m0 m1 1
F (v, t) ≥ − |v 0 (t)|2 + 2
kv(t)k2 + 3
kv(t)k4 + a(t, v, v) .
2 δ 16 β 32 β 2
k 1
If we take k > 0 such that < we obtain:
δ 4
1
(7.6) F (v, t) ≥ E(v, t) .
2
From (7.3) we also have
since
16 2
|v 0 (t)|2 ≤ 2 E(v, t) and |v(t)|2 ≤ β (t) E(v, t) .
m0
From (7.7) we obtain:
µ ¶
k+δ
F (v, t) ≤ (1 + 2 k) E(v, t) + 16 β 2 E(v, t)
m0
µ ¶
1 + 2k 2 k+δ
≤ β E(v, t) + 16 β 2 E(v, t) ,
β0 m0
with
1 + 2 k 16 (k+δ)
C0 = + .
β02 m0
Taking the time derivative of F (v, t) we have:
dF dE h i
(7.9) = + k 2 (v 00 (t), v(t)) + 2 |v 0 (t)|2 + 2 δ (v 0 (t), v(t)) .
dt dt
From the first estimate in Section 7, we have:
dE δ
(7.10) ≤ − |v 0 (t)|2 .
dt 2
Multiplying both sides of (6.5)1 by v and integrating on (0, 1), we obtain:
µ ¶
1 1
(7.11) (v 00 (t), v(t)) + M̌ kv(t)k2 kv(t)k2 + a(t, v(t), v(t)) +
4 β2 2β
µ ¶ µ ¶
∂v 0 ∂v
+ b(y, t) , v(t) + c(y, t) , v(t) + δ(v 0 (t), v(t)) = 0 .
∂y ∂y
We verify that:
¯µ 0 ¶¯ µ 0 ¶2
¯ ¯
(7.12) ¯ b(y, t) ∂v , v(t) ¯ ≤ 1 β kv(t)k2 + 9 |v 0 (t)|2 ,
¯ ∂y ¯ 4 β
¯µ ¶¯ µ 0 ¶2
¯ ¯
(7.13) ¯ c(y, t) ∂v , v(t) ¯ ≤ 1 β kv(t)k2 ,
¯ ∂y ¯ 2 β
µ ¶
1 m0
(7.14) a(t, v(t), v(t)) ≥ − β 02 kv(t)k2 .
4 β2 2
We have:
m0
0 < a(t, v(t), v(t)) ≤ kv(t)k2
8 β2
which substituted in (7.2) gives:
1 0 2 3 m0 m1
E(v, t) ≤ |v (t)| + 2
kv(t)k2 + kv(t)k4 ,
2 16 β 32 β 3
µ ¶
δ δ δ m0 m1
(7.18) E(v, t) ≤ |v 0 (t)|2 + 2
kv(t)k2 + kv(t)k4 .
20 4 20 16 β 16 β 3
dF δ
(7.19) + E(v, t) ≤ 0 for all t ≥ 0 .
dt 60
From (7.8) and (7.19) we obtain:
Rt
− 60δC ds
F (v, t) ≤ F (v0 , 0) e 0 0 β(s)2 .
Applications
We have: m−2
m C2
E(v, t) ≤ Cm e− m−2 (t+t0 )
m
β(t) = (t + t0 )1/2 .
We obtain:
E(v, t) ≤ C1 (t + t0 )−C2 , C2 > 0 ,
which gives an algebraic decay.
c(λ) = m0 λ + m1 λ2 . Then
We have M
2 2
" Z
1 1 +βµ β 0 x ∂u
¶2
1 c
µZ +βµ ∂u ¶2 ¶
E(u, t) = + u0 dx + M dx
2β 2 −β β ∂x 2 −β ∂x
(7.20) #
Z +βµ m µ
β0x
¶2 ¶ µ
∂u
¶2
0
+ − dx .
−β 2 β ∂x
where
1
Z +βµ β 0 x ∂u
¶2
1 c
µZ +βµ ∂u ¶2 ¶
b
E(u, t) = + u0 dx + M dx .
2 −β β ∂x 2 −β ∂x
b
The behaviour of E(u, t) depends essentially on β because
Rt
1 b C2 ds
0 β(s)2
E(u, t) ≤ β C1 e .
2
For example, if β(t) = (1+t)1/2 , then
1 b
E(u, t) ≤ C1 (1 + t)2−C2 ,
2
and we need C2 > 2 to have a decay of algebraic type. Note that C2 depends on
m0 , δ and β0 .
Observe that
µ ¶2 µ ¶2 µ ¶2
0 2 0 β 0 ∂u β 0 ∂u 1 β 0 ∂u 1 β 0 ∂u
(u ) = u + x − x ≤ u0 + x + x .
β ∂x β ∂x 2 β ∂x 2 β ∂x
Then, ÃZ
+β µZ +βµ ∂u ¶2 ¶!
1 0 2 1 c
E(u, t) ≥ (u ) dx + M dx .
2β −β 2 −β ∂x
If we represent by
Z β µZ β µ ¶2 ¶
b
b 1 c ∂u
E(u, t) = (u0 )2 dx + M dx
−β 2 −β ∂x
we have
1 bb
E(u, t) ≤ E(u, t), for all t ≥ 0 .
2β
Now suppose that β 0 is strictly positive and
δ 1
0 < β 0 (t) < for all t ≥ 0
120 C0 β
KIRCHHOFF–CARRIER ELASTIC STRINGS 499
with µ ¶1/2
δ m0
<C .
120 C0 β0 2
Note that δ, m0 , C0 , C are the constants fixed above.
We obtain Z t Z t 0
δ ds β
− 2
≤ −2 ds
60 C0 0 β (s) 0 β
or µ Z t ¶
δ ds β02
exp − < .
60 C 0 β 2 β(t)2
Substituting in the inequality of Theorem 7.1, we obtain:
C0 β 0
E(v, t) ≤ E(v, 0) .
β 2 (t)
bb
But E(u, t) < 2 β E(v, t). This gives
bb 2 C0 β04 E(v0 , 0)
E(u, t) ≤ .
β(t)
BIBLIOGRAPHY
[1] Arosio, A. and Spagnolo, S. – Global solutions to the Cauchy problem for
a nonlinear hyperbolic equation, in “Nonlinear Partial Differential Equations and
Their Applications”, Collège de France Seminar (H. Brezis and J.L. Lions, Eds.),
Pitman Research Notes in Mathematics Series, Vol. 16, London, 1984, pp. 1–24.
[2] Arosio, A. and Garivaldi, S. – On degenerate Kirchhoff strings, Math. Appl.
Science, 14 (1991), 177–195.
[3] Brezis, H. – Analyse Fonctionnelle (Théorie et Applications), Masson, Paris,
1983.
[4] Brito, E. – The damped elastic stretched string equation generalized: existence,
uniqueness, regularity and stability, Appl. Anal., 13 (1982), 219–233.
[5] Carrier, C.F. – On the vibrations problem of elastic string, Q. J. Appl. Math.,
3 (1945), 151–165.
[6] Bainov, D. and Minchev, E. – An estimate of the interval of existence of so-
lutions of the nonlinear Kirchhoff equation, Computers Math. Appl., 31(2) (1996),
63–65.
500 J. LIMACO FERREL and L.A. MEDEIROS
[7] Ebihara, Y., Medeiros, L.A. and Milla Miranda, M. – Local solutions for
a nonlinear degenerate hyperbolic equation, Nonlinear Analysis, 10 (1986), 27–40.
[8] Frota, C. and Larkin, N.A. – On global existence and uniqueness for unilateral
problem associated to the degenerate Kirchhoff equation, Nonlinear Analysis, 3
(1997), 443–452.
[9] Hosya, M. and Yamada, Y. – On some nonlinear wave equations. II – Global
existence and energy decay of solutions, Journal of the Faculty of Science, The
University of Tokyo, Sec IA, 38(2) (1991), 239–250.
[10] Kirchhoff, G. – Vorlesungen über Mechanik, ch. 9, §7, p. 144, Leipzig, Tauber,
1883.
[11] Komornik, V. and Zuazua, E. – A direct method for the boundary stabilization
of the wave equation, J. Math. Pure Appl., 69 (1980), 33–54.
[12] Kouemou Patcheu, S. – Existence globale et decroissance exponentielle de l’éner-
gie d’une équation quasilinéaire, C. R. Acad. Sci. Paris, t. 322, Série I (1996), 631–
634.
[13] Lions, J.L. – On some questions in boundary value problems of Mathematical
Physics, in “Contemporary Development in Continuum Mechanics and Partial Dif-
ferential Equations”, (G.M. de La Penha and L.A. Medeiros, Eds.), North Holland,
Amsterdam, (1978), pp. 284–329.
[14] Lions, J.L. – Quelques Méthodes de Résolution des Problèmes aux Limites Non
Linéaires, Dunod, Paris, 1969.
[15] Nakao, M. and Narazaki, T. – Existence and decay of solutions of some non-
linear wave equations in noncylindrical domains, Math. Rep., XI (1978), 117–125.
[16] Pohozhaev, S. – On a class of quasilinear hyperbolic equations, Math. Sbornik,
95 (1975), 152–166.
[17] Sidelnik, Ya.I. – Existence and uniqueness of a generalized solution of the mixed
problem for an equation of plate oscillation type in noncylindrical domain, Math.
Fisiko Mekhanicheskic, Polya, 28 (1988), 87–89.
[18] Yamada, Y. – Some nonlinear degenerate equations, Nonlinear Analysis, 11(10)
(1987), 1155–1168.
J. Limaco Ferrel,
Instituto de Matemática UFF,
Rua São Paulo s/n, CEP 24210-110 Niterói RJ – BRASIL
and
L.A. Medeiros,
Instituto de Matemática UFRJ,
C.P. 68530 – CEP 21944-970, Rio de Janeiro, RJ – BRASIL