• No results found

Blow-up of solutions for the sixth-order thin film equation with positive initial energy

N/A
N/A
Protected

Academic year: 2022

Share "Blow-up of solutions for the sixth-order thin film equation with positive initial energy"

Copied!
6
0
0

Loading.... (view fulltext now)

Full text

(1)

P

RAMANA c Indian Academy of Sciences Vol. 85, No. 4

— journal of October 2015

physics pp. 577–582

Blow-up of solutions for the sixth-order thin film equation with positive initial energy

WENJUN LIUand KEWANG CHEN

College of Mathematics and Statistics, Nanjing University of Information Science and Technology, Nanjing 210044, China

Corresponding author. E-mail: wjliu@nuist.edu.cn

MS received 29 October 2013; revised 06 August 2014; accepted 21 August 2014 DOI:10.1007/s12043-014-0907-2; ePublication:20 May 2015

Abstract. In this paper, a sixth-order parabolic thin film equation with the initial boundary condi- tion is considered. By using the improved energy estimate method and by constructing second-order elliptic problem, a blow-up result for certain solution with positive initial energy is established, which is an improve over the previous result of Li and Liu.

Keywords.Blow-up; sixth-order thin film equation; positive initial energy.

PACS Nos 02.30.Jr; 02.30.Sa

1. Introduction and main result

In the last 20 years, higher-order nonlinear parabolic partial differential equations (PDEs), as models for applications in mechanics and physics, have become more common in the study on pure and applied PDEs. For instance, the pure sixth-order parabolic thin film equation (TFE) was first introduced in [1,2] to describe the spreading of a thin viscous fluid (with possible slip at the solid interface) under the driving force of an elastica (or light plate). There are many related works on blow-up solutions to these kinds of parabolic equations and systems (see [3–9] and the references therein).

In this paper, we consider the following initial boundary problem of the sixth-order TFE:

ut−(2u− |u|p−1u)=0, x∈, t∈(0, T ), u=u=2u=0, x∈∂, t∈ [0, T ),

u=u0(x), x∈, t=0,

(1) where⊂Rn, n≥ 1 is a bounded smooth domain,p >1. Recently, Li and Liu [10]

defined the energy Lyapunov functional E(t )= 1

2

|u(x, t )|2dx− 1 p+1

|u|p+1dx, t≥0, (2)

(2)

and proved that if the initial datumu0∈C6+α()¯ with the initial energy E(0)= 1

2

|u0|2dx− 1 p+1

|u0|p+1dx ≤0, (3) then the solution to problem (1) should blow up in finite time. In this paper, the above result is improved to show that certain solutions with positive initial energy can also blow up in finite time. For a Banach spaceX,·Xdenotes the norm ofX. For simplicity, we denote · L2()by · 2.

For our purpose, it is assumed that B is the optimal constant of the embedding inequality

up+1≤Bu2, u∈H02() (4) i.e.,

B−1= inf

u∈H2 0() u =0

u2

up+1

,

and set

α1=B−(p+1)/(p−1)

, E1= 1

2 − 1

p+1

B−(2(p+1))/(p−1)

. (5)

Then, from [10] we get (with some minor corrections) E(t )= −

|∇(2u− |u|p−1u)|2dx≤0, fort ≥0, (6) which implies that the equilibrium is stable.

Our main result reads as follows.

Theorem 1. Assume that the initial datumu0∈C6+α()¯ satisfy

E(0) < E1 (7)

and

u02> α1. (8)

Then the solutionu(x, t )of problem(1)blows up in a finite time.

2. Proof of Theorem 1

We shall use the improved energy estimate method, which has been successfully applied in [11] to deal with the second-orderp-Laplacian equation (see also [12–14] for further applications of this method). To extend the method to the sixth-order thin film equation herein, we should combine it with the construction of a second-order elliptic problem (see (18)).

We first prove the following lemmas by applying the idea of Vitillaro in [15] where a different type of equation was discussed. The first lemma gives a lower bound estimate ofu2.

(3)

Lemma2. Suppose u is a solution of the system (1). Assume that E(0) < E1 and u02> α1. Then there exists a positive constantα2> α1,such that

u2≥α2, ∀t ≥0 (9)

and

up+1≥Bα2, ∀t ≥0. (10) Proof. We first note that, by (2) and (4),

E(t ) ≥ 1

2u22− 1

p+1Bp+1up+12

= 1

2− 1

p+1Bp+1αp+1=:g(α), (11)

whereα= u2. The method used here is also related to the so-called Fibering method introduced by Pohozev in the 1970s [16–19]. It is easy to verify thatgis increasing for 0< α < α1, decreasing forα > α1;g(α)→ −∞asα→ +∞andg(α1)=E1, where α1 is given in (5). AsE(0) < E1, there existsα2 > α1 such thatg(α2) = E(0). Let α0 = u02, then by (11) we haveg(α0)≤E(0)=g(α2), which implies thatα0≥α2. To establish (9), we suppose by contradiction thatu(·, t0)2 < α2for somet0 >0.

By the continuity ofu(·, t )2we can chooset0such thatu(·, t0)2 > α1. It follows from (11) that

E(t0)≥g(u(·, t0)2) > g(α2)=E(0).

This is impossible becauseE(t )≤E(0)for allt ≥0. Hence (9) is established.

To prove (10), we exploit (2) and (6) to obtain that 1

p+1up+1p+1≥ 1

2u22−E(0). (12)

Consequently, by (9), we have 1

p+1up+1p+1≥ 1

22−E(0)≥ 1

22−g(α2)= 1

p+1Bp+1αp+12 . (13) Therefore (10) is concluded.

In the remainder of this section we consider the case thatE(0) < E1andu02> α1. We set

H (t )=E1−E(t ), t ≥0. (14)

Then we have the following lemma.

Lemma3. For allt ≥0,

0< H (0)≤H (t )≤ 1

p+1up+1p+1. (15)

Proof. By (6) we see thatH(t )≥0. Thus

H (t )≥H (0)=E1−E(0) >0, t≥0. (16)

(4)

From (2), we obtain H (t )=E1−1

2u22+ 1

p+1up+1p+1, and exploiting (9) and (5) we get

E1−1

2u22≤E1−1

12= − 1

p+1Bp+1α1p+1<0, ∀t≥0.

Hence

H (t )≤ 1

p+1up+1p+1, ∀t ≥0. (17)

Then (15) follows from (16) and (17).

Letφbe the unique solution to −φ=u, x∈,

φ=0, x ∈∂. (18)

Due to the ellipticL2-theory, we have

∇φ22 ≤Cu22. (19)

Completion of the proof of Theorem1. We define G(t ):=1

2

|∇φ (x, t )|2dx. (20) DifferentiatingG(t )and exploiting (18), (1), Green’s first formula

uvdx=

u∂v

∂nds−

∇u· ∇vdx and Green’s second formula

(uv−vu)dx=

u∂v

∂n−v∂u

∂n

ds, we get

G(t ) =

∇φ· ∇φtdx=

φ∂φt

∂nds−

φφtdx

=

φutdx=

φ(2u− |u|p−1u)dx

=

φ (2u− |u|p−1u)dx−

(2u− |u|p−1u)∂φ

∂n

−φ∂(2u− |u|p−1u)

∂n

ds

= −

u(2u− |u|p−1u)dx = −

|u|2dx+

|u|p+1dx +

u∂u

∂n−u∂(u)

∂n

ds

= −

|u|2dx+

|u|p+1dx. (21)

(5)

From (2) and (14), we get

|u|2dx = −2E1− 2 p+1

|u|p+1dx+2H (t ). (22) It follows from (21) and (22) that

G(t )= −2E1+

1− 2 p+1

up+1p+1+2H (t ). (23)

By using (5) and (10), we have 2E1 = 2(p−1)

p+1 α21= 2(p−1)

p+1 Bp+1α1p+1

= αp+11 αp+12

2(p−1)

p+1 Bp+1α2p+1

≤ αp+11 αp+12

2(p−1)

p+1 up+1p+1. (24)

It follows from (23), (24) and (15) that G(t ) ≥ 2(p−1)

p+1

1−α1p+1 α2p+1

up+1p+1+2H (t )

= C0up+1p+1+2H (t )≥C0up+1p+1≥0, (25) where

C0= 2(p−1) p+1

1−α1p+1 α2p+1

>0 due to the fact thatp >1 andα2 > α1.

Next, we use (19) and Hölder’s inequality to estimateG(p+1)/2(t )as

G(p+1)/2(t )≤Cup+12 ≤C||(p−1)/2up+1p+1. (26) By combining (25) and (26) we have

G(t )≥γ G(p+1)/2(t ), (27)

whereγ =C0/

C||(p−1)/2

. A direct integration of (27) then yields G(p−1)/2(t )≥ 1

G(1−p)/2(0)−(p−1)γ t /2.

Therefore G(t ) blows up in a time t ≤ (2G(1−p)/2(0))/((p−1)γ ). So does u22, according to (19).

Acknowledgements

This work was partly supported by the National Natural Science Foundation of China (Grant No. 11301277), the Qing Lan Project of Jiangsu Province, the Overseas Scholar- ship of Jiangsu Provincial Government, and the Training Abroad Project of Outstanding

(6)

Young and Middle-Aged University Teachers and Presidents. The authors wish to thank the anonymous referees and the editor for their valuable comments.

References

[1] J R King,Mathematical aspects of semiconductor process modelling, Ph.D. Thesis (University of Oxford, Oxford, 1986)

[2] J R King,SIAM J. Appl. Math.49, 1064 (1989)

[3] A L Bertozzi and M C Pugh,Comm. Pure Appl. Math.51, 625 (1998) [4] J D Evans, V A Galaktionov and J R King,Eur. J. Appl. Math.18, 195 (2007) [5] J D Evans, V A Galaktionov and J R King,Nonlinearity20, 1799 (2007) [6] P Souplet,Math. Methods Appl. Sci.19, 1317 (1996)

[7] F Q Sun,Electron. J. Differ. Eq.2010, 9 (2010) [8] S Q Yu,Appl. Anal.91, 1321 (2012)

[9] J P Zhao,Z. Naturforsch.67a, 479 (2012)

[10] Z B Li and C C Liu,Miskolc Math. Notes13, 429 (2012) [11] W J Liu and M X Wang,Acta Appl. Math.103, 141 (2008) [12] Z B Fang, L Sun and C Li,Bound. Value Probl.2013, 181 (2013) [13] W J Gao and Y Z Han,Appl. Math. Lett.24, 784 (2011)

[14] X L Wu and W J Gao,Acta Math. Sci. Ser. B Engl. Ed.33, 1024 (2013) [15] E Vitillaro,Arch. Ration. Mech. Anal.149, 155 (1999)

[16] P Álvarez-Caudevilla and V A Galaktionov,Adv. Nonlinear Stud.12, 315 (2012) [17] K J Brown and T-F Wu,Differential Integral Equations22, 1097 (2009)

[18] S I Pohožaev, An approach to nonlinear equations,Dokl. Akad. Nauk SSSR247, 1327 (1979), English trans1.:Sov. Math. Dokl.20, 912 (1979)

[19] S I Pohožaev, The fibering method in nonlinear variational problems, inTopological and variational methods for nonlinear boundary value problems (Cholín, 1995), pp. 35–88, Pitman Res. Notes Math. Ser., 365, Longman, Harlow

References

Related documents

In order to identify priorities for the implementation of renewable energy services and programs and to provide strategic sustainable energy solutions, four management units,

l Scale up to Access of Clean Energy Scheme (ACE): The Ministry of New and Renewable Energy (MNRE), is implementing a project titled “Scale Up of Access to Clean Energy (ACE)

wightii was treated in cadmium and lead solutions having different initial concentrations, sorption capacity of the biomass increased with increase in initial metal concentration

Introduction to advanced blow-moulding techniques: Blow moulding with multi- parison system, injection stretch blow moulding process, extrusion stretch blow moulding process,

only up to 5%, if there is no load. This is the result of Second Way of energy saving. However, if there is load F with constant direction, and if it can be used for support of

Abstract ; A method foi the determination of the activation energy E of a thermoluminescence (TL) peak obeying mixed order kinetics by using a set of expressions, is

guuld.aae and wfprrvle tam. The thesia his ieaohed the standa rd roltL l ltlg the requlrcmenfs of the r gulnUm r relatlm to the degree, the results a tstwd to the beats hive

the Airy’s differential equation is a second order differen- tial equation, it will have two linearly independent simultaneous solutions which are denoted by Ai(x) and Bi(x)..