SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research

Energy decay and nonexistence of solution for a reaction-diffusion equation with exponential nonlinearity

Huiya Dai12* and Hongwei Zhang2

Author Affiliations

1 School of Mathematics and Statistics, Xi’an Jiaotong University, Xianning Road, Xi’an, P.R. China

2 Department of Mathematics, Henan University of Technology, Lianhua Street, Zhengzhou, P.R. China

For all author emails, please log on.

Boundary Value Problems 2014, 2014:70  doi:10.1186/1687-2770-2014-70

The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2014/1/70


Received:2 December 2013
Accepted:13 March 2014
Published:25 March 2014

© 2014 Dai and Zhang; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Abstract

In this work we consider the energy decay result and nonexistence of global solution for a reaction-diffusion equation with generalized Lewis function and nonlinear exponential growth. There are very few works on the reaction-diffusion equation with exponential growth f as a reaction term by potential well theory. The ingredients used are essentially the Trudinger-Moser inequality.

Keywords:
reaction-diffusion equation; stable and unstable set; exponential reaction term; decay rate; global nonexistence

1 Introduction

In this paper, we study the following initial boundary value problem with generalized Lewis function a ( x , t ) which depends on both spacial variable and time:

a ( x , t ) u t Δ u = f ( u ) , x Ω , t > 0 , (1)

u ( x , t ) = 0 , x Ω , t > 0 , (2)

u ( x , 0 ) = u 0 ( x ) , x Ω , (3)

here f ( s ) is a reaction term with exponential growth at infinity to be specified later, Ω is a bounded domain with smooth boundary Ω in R 2 .

For the reaction-diffusion equation with polynomial growth reaction terms (that is, equation (1) with a ( x , t ) = 1 and f ( u ) = | u | p 1 u ), there have been many works in the literature; one can find a review of previous results in [1,2] and references therein, which are not listed in this paper just for concision. Problem (1)-(3) with a ( x , t ) > 0 describes the chemical reaction processes accompanied by diffusion [2]. The author of work [1] proved the existence and asymptotic estimates of global solutions and finite time blow-up of problem (1)-(3) with a ( x , t ) > 0 and the critical Sobolev exponent p = n + 2 n 2 for f ( u ) = u p .

In this paper we assume that f ( s ) is a reaction term with exponential growth like e s 2 at infinity. When a ( x , t ) = 1 , f ( u ) = e u , model (1)-(3) was proposed by [3] and [4]. In this case, Fujita [5] studied the asymptotic stability of the solution. Peral and Vazquez [6] and Pulkkinen [7] considered the stability and blow-up of the solution. Tello [8] and Ioku [9] considered the Cauchy problem of heat equation with f ( u ) e u 2 for | u | 1 .

Recently, Alves and Cavalcanti [10] were concerned with the nonlinear damped wave equation with exponential source. They proved global existence as well as blow-up of solutions in finite time by taking the initial data inside the potential well [11]. Moreover, they also got the optimal and uniform decay rates of the energy for global solutions.

Motivated by the ideas of [1,10], we concentrate on studying the uniform decay estimate of the energy and finite time blow-up property of problem (1)-(3) with generalized Lewis function a ( x , t ) and exponential growth f as a reaction term. To the authors’ best knowledge, there are very few works in the literature that take into account the reaction-diffusion equation with exponential growth f as a reaction term by potential well theory. The majority of works in the literature make use of the potential well theory when f possesses polynomial growth. See, for instance, the works [12-16] and a long list of references therein. The ingredients used in our proof are essentially the Trudinger-Moser inequality (see [17,18]). We establish decay rates of the energy by considering ideas from the work of Messaoudi [15]. The case of nonexistence results is also treated, where a finite time blow-up phenomenon is exhibited for finite energy solutions by the standard concavity method adapted for our context.

The remainder of our paper is organized as follows. In Section 2 we present the main assumptions and results, Section 3 and Section 4 are devoted to the proof of the main results.

Throughout this study, we denote by , p , H 0 1 the usual norms in spaces L 2 ( Ω ) , L p ( Ω ) and H 0 1 ( Ω ) , respectively.

2 Assumptions and preliminaries

In this section, we present the main assumptions and results. We always assume that:

(A1) a ( x , t ) is a positive differentiable function and is bounded for t [ 0 , + ) , x Ω .

(A2) f : R R is a C 1 function. The function f ( t ) / t is increasing in ( 0 , ) , and for each β > 0 , there exists a positive constant C β such that

| f ( t ) | C β e β t 2 , | f ( t ) | C β e β t 2 . (4)

(A3) For each ε > 0 , β > 0 and p > 1 fixed, there exists a positive constant C ( ε , β ) such that

| f ( t ) | ε | t | + C ( ε , β ) | t | p 1 e β t 2 , (5)

| F ( t ) | ε | t | 2 + C ( ε , β ) | t | p e β t 2 , (6)

where F ( t ) = 0 t f ( s ) d s .

(A4) There exists a positive constant θ > 2 such that

0 < θ F ( t ) < f ( t ) t , t R / { 0 } . (7)

A typical example of functions satisfies (A2)-(A4) is f ( t ) = C | t | p 1 t e M t α , with given p > 1 , M > 0 , C > 0 , and α ( 1 , 2 ) .

Now we define some functional as follows:

E ( t ) = E ( u ) = 1 2 u 2 Ω F ( u ) d x , (8)

I ( t ) = I ( u ) = u 2 Ω u f ( u ) d x , (9)

then the ‘potential depth’ given by

d = inf { sup λ R E ( λ u ) , u H 0 1 / { 0 } }

is a positive constant [10]. Hence, we are able to define stable and unstable sets respectively as follows:

W 1 = { u H 0 1 , E ( u ) < d , I ( u ) > 0 } , W 2 = { u H 0 1 , E ( u ) < d , I ( u ) < 0 } .

We also need the following lemmas.

Lemma 2.1[17,18]

Let Ω be a bounded domain in R 2 . For all u H 0 1 ( Ω ) ,

e α | u | 2 L 1 ( Ω ) for all  α > 0 , (10)

and there exist positive constants m 2 such that

sup u H 0 1 ( Ω ) , u 1 Ω e α | u | 2 d x = m 2 < for all  α 4 π . (11)

Lemma 2.2[19]

Let ϕ ( t ) be a nonincreasing and nonnegative function on [ 0 , ) , such that

sup s [ t , t + 1 ] ϕ ( s ) C ( ϕ ( t ) ϕ ( t + 1 ) ) , t > 0 , (12)

then

ϕ ( t ) C e ω t ,

whereC, ωare positive constants depending on ϕ ( 0 ) and other known qualities.

Lemma 2.3[20]

Suppose that a positive, twice-differentiable function H ( t ) satisfies on t 0 the inequality

H ( t ) H ( t ) ( δ + 1 ) ( H ( t ) ) 2 0 , (13)

where δ > 0 , then there is t 1 < t 2 = H ( 0 ) δ H ( 0 ) such that H ( t ) as t t 1 .

In order to state and prove our main results, we remind that by the embedding theorem there exists a constant C 0 depending on p and Ω only such that

u p C 0 u . (14)

By multiplying equation (1) by u t , integrating over Ω, using integration by parts and a ( x , t ) > 0 , we get

E ( t ) = Ω a ( x , t ) u t 2 ( x , t ) d x 0 . (15)

Our main results read as follows.

Theorem 2.1Let (A1)-(A4) hold. Assume further that u 0 W 1 satisfies

ε 0 C 0 2 + C ε 0 C 0 p m 2 1 2 ( 2 θ θ 2 E ( 0 ) ) p 2 < 1 (16)

for some sufficiently small ε 0 > 0 and C ε 0 > 0 . Then there exist positive constantsKandksuch that the energy E ( t ) satisfies the decay estimates for larget

E ( t ) K e k t . (17)

Theorem 2.2Let (A1)-(A4) hold. Assume further that a t ( x , t ) 0 , u 0 W 2 and E ( 0 ) < ( θ 2 ) d θ < d , then the solutions of (1)-(3) blow up in finite time.

3 Proof of decay of the energy

In this section we prove Theorem 2.1. We divide the proof into two lemmas.

Lemma 3.1Under the assumptions of Theorem 2.1, we have, for all t 0 , u ( t ) W 1 .

Proof Since I ( u 0 ) 0 , then there exists (by continuity) T m < T such that

I ( u ( t ) ) 0 , t [ 0 , T m ] .

This and (A4) give

E ( t ) = ( 1 2 1 θ ) u 2 + 1 θ I ( u ) + Ω ( 1 θ u f ( u ) F ( u ) ) d x ( 1 2 1 θ ) u 2 , t [ 0 , T m ] . (18)

So, by (15) we have

u 2 2 θ θ 2 E ( t ) 2 θ θ 2 E ( 0 ) 2 θ d θ 2 , t [ 0 , T m ] . (19)

We then use (5), the Holder inequality and the embedding theorem to obtain, for each t [ 0 , T m ] ,

Ω u f ( u ) d x Ω [ ε | u | 2 + C ( ε , β ) | u | p e β u 2 ] d x ε C 0 2 u 2 + C ( ε , β ) ( Ω | u | 2 p d x ) 1 2 ( Ω e 2 β u 2 d x ) 1 2 . (20)

Once u 2 2 θ d θ 2 , we choose β such that θ β d θ 2 < π , then, from Trudinger-Moser inequality (11),

Ω e 2 β | u | 2 d x Ω e 2 β u 2 ( u u ) 2 d x m 2 ,

and therefore, by (16) for ε 0 > 0 and C ε 0 > 0 , we have

Ω u f ( u ) d x ε 0 C 0 2 u 2 + C ε 0 C 0 p m 2 1 2 u p = ε 0 C 0 2 u 2 + C ε 0 C 0 p m 2 1 2 ( 2 θ θ 2 E ( t ) ) p 2 u 2 [ ε 0 C 0 2 + C ε 0 C 0 p m 2 1 2 ( 2 θ θ 2 E ( 0 ) ) p 2 ] u 2 < u 2 . (21)

By virtue of (21) and the definition of I ( t ) , we have

I ( t ) = u 2 Ω u f ( u ) d x > 0 .

This shows that u ( t ) W 1 for all t [ 0 , T m ] . By repeating this procedure and the fact that E ( t ) E ( 0 ) , we obtain

lim t T m ( ε 0 C 0 2 + C ε 0 C 0 p m 2 1 2 ( 2 θ θ 2 E ( t ) ) p 2 ) ρ < 1 .

This is extended to T. □

Lemma 3.2Under the assumptions of Theorem 2.1, we have, for η = 1 [ ε 0 C 0 2 + C ε 0 C 0 p m 2 1 2 ( 2 θ θ 2 E ( 0 ) ) p 2 ] ,

η u 2 < I ( t ) . (22)

Proof It suffices to rewrite (21) as

Ω u f ( u ) d x [ ε 0 C 0 2 + C ε 0 C 0 p m 2 1 2 ( 2 θ θ 2 E ( 0 ) ) p 2 ] u 2 = ( 1 η ) u 2 = u 2 η u 2 . (23)

Thus (22) follows from (23). □

Proof of Theorem 2.1 We integrate (15) over [ t , t + 1 ] to obtain

E ( t ) E ( t + 1 ) = t t + 1 Ω a ( x , s ) u t 2 ( x , s ) d x d s = D 2 ( t ) . (24)

Now we multiply (1) by u and integrate over Ω × [ t , t + 1 ] to arrive at

t t + 1 I ( s ) d s = t t + 1 [ u 2 Ω u f ( u ) d x ] d s = t t + 1 Ω a ( x , t ) u t ( x , t ) u ( x , t ) d x d s A t t + 1 a 1 2 u t ( s ) u ( s ) d s , (25)

where A 2 = sup ( x , t ) Ω × [ 0 , + ) | a ( x , t ) | . Exploiting (14) and (19), we obtain

t t + 1 I ( s ) d s A C 0 ( 2 θ θ 2 ) 1 2 ( sup s [ t , t + 1 ] E 1 2 ( t ) ) t t + 1 a 1 2 u t d s A C 0 ( 2 θ θ 2 ) 1 2 ( sup s [ t , t + 1 ] E 1 2 ( t ) ) D ( t ) . (26)

Using (7), (23) and (22), we have

E ( t ) = θ 2 2 θ u 2 + 1 θ I ( t ) + Ω ( 1 θ u f ( u ) F ( u ) ) d x θ 2 2 θ u 2 + 1 θ I ( t ) + Ω 2 θ u f ( u ) d x θ 2 2 θ u 2 + 1 θ I ( t ) + 2 θ ( 1 η ) u 2 [ θ 2 2 θ η + 1 θ + 2 θ ( 1 η ) ] I ( t ) . (27)

Integrating both sides of (27) over [ t , t + 1 ] and using (26), one can write

t t + 1 E ( s ) d s = [ θ 2 2 θ η + 1 θ + 2 θ ( 1 η ) ] A C 0 ( 2 θ θ 2 ) 1 2 ( sup s [ t , t + 1 ] E 1 2 ( t ) ) D ( t ) . (28)

By using (15) again, we have E ( s ) E ( t + 1 ) , s t + 1 , hence

t t + 1 E ( s ) d s E ( t + 1 ) . (29)

Inserting (29) in (24) and using (27), we easily have

E ( t ) t t + 1 E ( s ) d s + t t + 1 Ω a ( x , s ) u t 2 ( x , s ) d x d s [ θ 2 2 θ η + 1 θ + 2 θ ( 1 η ) ] A C 0 ( 2 θ θ 2 ) 1 2 ( sup s [ t , t + 1 ] E 1 2 ( t ) ) D ( t ) + D 2 ( t ) C 1 [ E 1 2 ( t ) D ( t ) + D 2 ( t ) ] (30)

for C 1 a constant depending on C 0 , A, θ, η only. We then use Young’s inequality to get from (30) and (24)

sup s [ t , t + 1 ] E ( t ) C 2 D 2 ( t ) C 2 ( E ( t ) E ( t + 1 ) ) . (31)

By (12) in Lemma 2.2 we then get the results. □

4 Proof of the blow-up result

In this section, we shall prove Theorem 2.2 by adapting the concavity method (see Levine [20]). We recall the following lemma in [10].

Lemma 4.1[10]

Assume that u 0 W 2 and E ( 0 ) < d , then it holds that

u ( t ) W 2 for  t [ 0 , T max ) , (32)

u 2 2 d for  t [ 0 , T max ) . (33)

Proof of Theorem 2.2 Assume by contradiction that the solution is global. Then, for any T > 0 , we consider the function H ( ) : [ 0 , T ] R + defined by

H ( t ) = 0 t Ω a ( x , s ) u 2 ( x , s ) d x d s + 0 t Ω ( s t ) a t ( x , s ) u 2 ( x , s ) d x d s + ( T t ) Ω a ( x , 0 ) u 0 2 ( x ) d x + ρ ( t + t 0 ) 2 , (34)

where t 0 , T, ρ are positive constants which will be fixed later. Direct computations show that

H ( t ) = Ω a ( x , t ) u 2 ( x , t ) d x 0 t Ω a t ( x , s ) u 2 ( x , s ) d x d s Ω a ( x , 0 ) u 0 2 d x + 2 ρ ( t + t 0 ) H ( t ) = 2 0 t Ω a ( x , s ) u t ( x , s ) u ( x , s ) d x d s + 2 ρ ( t + t 0 ) , (35)

H ( t ) = 2 Ω a ( x , t ) u ( x , t ) u t ( x , s ) d x + 2 ρ . (36)

Then, due to equations (1), (7) and (33), we have

H ( t ) = 2 u 2 + 2 Ω u f ( u ) d x + 2 ρ 2 u 2 + 2 θ Ω F ( u ) d x + 2 ρ = ( θ 2 ) u 2 2 θ E ( t ) + 2 ρ = ( θ 2 ) u 2 2 θ E ( 0 ) + 2 θ 0 t Ω a ( x , s ) u t 2 ( x , s ) d x d s + 2 ρ 2 ( θ 2 ) d 2 θ E ( 0 ) + 2 θ 0 t Ω a ( x , s ) u t 2 ( x , s ) d x d s + 2 ρ . (37)

Now we take 0 < ρ < ( θ 2 ) d θ E ( 0 ) θ 1 such that 2 ( θ 2 ) d 2 θ E ( 0 ) + 2 ρ > 2 θ ρ (this ρ can be chosen since E ( 0 ) < ( θ 2 ) d θ ), and then

H ( t ) 2 θ ρ + 2 θ 0 t Ω a ( x , s ) u t 2 ( x , s ) d x d s . (38)

We also note that

H ( 0 ) = T Ω a ( x , 0 ) u 0 2 ( x ) d x + ρ t 0 2 > 0 , H ( 0 ) = 2 ρ t 0 > 0 , H ( t ) 2 θ ρ > 0 , t 0 .

Therefore H ( t ) and H ( t ) are both positive. Since a t ( x , t ) 0 for all x Ω and t 0 , by the construction of H ( t ) , it is clear that

H ( t ) 0 t Ω a ( x , s ) u 2 ( x , s ) d x d s + ρ ( t + t 0 ) 2 . (39)

Thus, for all ( ξ , η ) R 2 , from (35), (38) and (39) it follows that

H ( t ) ξ 2 + H ( t ) ξ η + 1 2 θ H ( t ) η 2 ( 0 t Ω a ( x , s ) u ( x , s ) 2 d x d s + ρ ( t + t 0 ) 2 ) ξ 2 + 2 ξ η 0 t Ω a ( x , s ) u ( x , s ) u t ( x , s ) d x d s + 2 ρ ( t + t 0 ) ξ η + ρ η 2 + η 2 0 t Ω a ( x , s ) u t 2 ( x , s ) d x d s 0 ,

which implies

( H ( t ) ) 2 2 θ H ( t ) H ( t ) 0 .

That is,

H ( t ) H ( t ) θ 2 ( H ( t ) ) 2 0 .

Then we complete the proof by the standard concavity method (Lemma 2.3) since θ > 2 . □

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The authors declare that the study was realized in collaboration with the same responsibility. All authors read and approved the final manuscript.

Acknowledgements

We thank the referees for their valuable suggestions which helped us improve the paper so much. This work was supported by the National Natural Science Foundation of China (11171311).

References

  1. Tan, Z: The reaction-diffusion equation with Lewis function and critical Sobolev exponent. J. Math. Anal. Appl.. 2, 480–495 (2002)

  2. Pao, CV: Nonlinear Parabolic and Elliptic Equations, Plenum, New York (1992)

  3. Frank-Kamenetskii, DA: Diffusion and Heat Transfer in Chemical Kinetics, Plenum, New York (1969)

  4. Gel’fand, IM: Some problems in the theory of quasilinear equations. Usp. Mat. Nauk. 2, 87–158 (1959) English transl.: Am. Math. Soc. Transl. (2) 29, 295-381 (1963)

  5. Fujuita, H: On the nonlinear equations δ u + e u = 0 and v / t = δ v + e v . Bull. Am. Math. Soc.. 75, 132–135 (1969). Publisher Full Text OpenURL

  6. Peral, I, Vazquez, JL: On the stability or instability of the singular solution of the semilinear heat equation with exponential reaction term. Arch. Ration. Mech. Anal.. 129, 201–224 (1995). Publisher Full Text OpenURL

  7. Pulkkinen, A: Blow-up profiles of solutions for the exponential reaction-diffusion equation. Math. Methods Appl. Sci.. 34, 2011–2030 (2011). Publisher Full Text OpenURL

  8. Tello, J: Stability of steady states of the Cauchy problem for the exponential reaction-diffusion equation. J. Math. Anal. Appl.. 324, 381–396 (2006). Publisher Full Text OpenURL

  9. Ioku, N: The Cauchy problem for heat equations with exponential nonlinearity. J. Differ. Equ.. 251, 1172–1194 (2011). Publisher Full Text OpenURL

  10. Alves, CO, Cavalcanti, MM: On existence, uniform decay rates and blow up for solutions of the 2-d wave equation with exponential source. Calc. Var.. 34, 377–411 (2009). Publisher Full Text OpenURL

  11. Sattiger, DH: On global solutions of nonlinear hyperbolic equations. Arch. Ration. Mech. Anal.. 30, 148–172 (1968)

  12. Tsutsumi, M: Existence and nonexistence of global solutions for nonlinear parabolic equations. Publ. Res. Inst. Math. Sci.. 8, 211–229 (1972/73). Publisher Full Text OpenURL

  13. Iekhata, R, Suzuki, T: Stable and unstable sets for evolution equations of parabolic and hyperbolic type. Hiroshima Math. J.. 26, 475–491 (1996)

  14. Levine, HA, Park, SR, Serrin, J: Global existence and nonexistence theorems for quasilinear evolution equations of formally parabolic type. J. Differ. Equ.. 142, 212–229 (1998). Publisher Full Text OpenURL

  15. Messaoudi, SA: A decay result for a quasilinear parabolic system type. Elliptic and Parabolic Problems. 43–50 (2005)

  16. Chen, H, Liu, GW: Global existence and nonexistence theorems for semilinear parabolic equations with conical degeneration. J. Pseud.-Differ. Oper. Appl.. 3(3), 329–349 (2012). Publisher Full Text OpenURL

  17. Moser, J: A sharp form of an inequality by N. Trudinger. Indiana Univ. Math. J.. 20, 1077–1092 (1971). Publisher Full Text OpenURL

  18. Trudinger, NS: On imbeddings into Orlicz spaces and some applications. J. Math. Mech.. 17, 473–483 (1967)

  19. Nakao, M, Ono, K: Global existence to the Cauchy problem of the semi-linear evolution equations with a nonlinear dissipation. Funkc. Ekvacioj. 38, 417–431 (1995)

  20. Levine, HA: Some nonexistence and instability theorems for solutions of formally parabolic equations of the form p u t = a u + f ( u ) . Arch. Ration. Mech. Anal.. 51, 371–386 (1973)