This paper deals with blow-up and global solutions of the following nonlinear reaction-diffusion equations under nonlinear boundary conditions: g(u)t=∇·au∇u+fu in Ω×0,T,∂u/∂n=bx,u,t on ∂Ω×(0,T),u(x,0)=u0(x)>0, in Ω¯, where Ω⊂RN(N≥2) is a bounded domain with smooth boundary ∂Ω. We obtain the conditions under which the solutions either exist globally or blow up in a finite time by constructing auxiliary functions and using maximum principles. Moreover, the upper estimates of the “blow-up time,” the “blow-up rate,” and the global solutions are also given.
1. Introduction
During the past few decades, the blow-up phenomena for the nonlinear reaction-diffusion equations have been studied by a large number of authors, and the reader is referred to [1–8] and the references therein. In this paper, we consider the following nonlinear reaction-diffusion problem under nonlinear conditions:(1)gut=∇·au∇u+fuin Ω×0,T,∂u∂n=bx,u,ton ∂Ω×0,T,ux,0=u0x>0in Ω¯,where Ω⊂RN(N≥2) is a bounded domain with smooth boundary ∂Ω, ∂u/∂n represents the outward normal derivative on ∂Ω, and T is the maximal existence time of u. Set R+≔(0,+∞). We assume, throughout the paper, that a(s) is a positive C2(R+) function, g(s) is a C2(R+) function, g′(s)>0 for any s∈R+, f(s) is a nonnegative C1(R+) function, b(x,s,t) is a nonnegative C1(∂Ω×R+×R+) function, and u0(x) is a positive C2(R+¯) function and satisfies the compatibility conditions. Under the above assumptions, the local existence and uniqueness of classical solution of problem (1) were established by Amann [9]. Furthermore, it follows from maximum principle [10] and regularity theorem [11] that the solution u(x,t) is positive and u(x,t)∈C3(Ω×0,T)∩C2(Ω¯×(0,T)).
Many authors have investigated blow-up and global solutions of nonlinear reaction-diffusion equations under nonlinear boundary conditions and have obtained a lot of interesting results (see, e.g., [12–20]). To my knowledge, some special cases of (1) have been studied. Zhang [21] considered the following problem: (2)gut=Δu+fuinΩ×0,T,∂u∂n=buon ∂Ω×0,T,ux,0=u0x>0in Ω¯,where Ω⊂RN(N≥2) is a bounded domain with smooth boundary ∂Ω. By constructing auxiliary functions and using maximum principles, the existence of blow-up and global solutions were obtained under appropriate assumptions on the functions b,f,g, and u0. Zhang et al. [22] dealt with the following problem: (3)gut=∇·au∇u+fuin Ω×0,T,∂u∂n=buon ∂Ω×0,T,ux,0=u0x>0in Ω¯,where Ω⊂RN(N≥2) is a bounded domain with smooth boundary ∂Ω. Some conditions on nonlinearities and the initial data were given to ensure that u(x,t) exists globally or blows up at some finite time T. In addition, the upper estimates of the global solution, the “blow-up time,” and the “blow-up rate” were also established.
In this paper, we study reaction-diffusion problem (1). It is well known that f(u), g(u), a(u), and b(x,u,t) are nonlinear reaction, nonlinear diffusion, nonlinear convection, and nonlinear boundary flux, respectively. What interactions among the four nonlinear mechanisms result in the blow-up and global solutions of (1) is investigated in this work. We note that the boundary flux function b(x,u,t) depends not only on the concentration variable u but also on the space variable x and the time variable t. Hence, it seems that the methods of [21, 22] are not applicable for problem (1). In this paper, by constructing completely different auxiliary functions from those in [21, 22] and technically using maximum principles, we obtain the existence theorems of the blow-up and global solution. Moreover, the upper estimates of “blow-up time,” “blow-up rate,” and global solution are also given. Our results can be seen as the extension and supplement of those obtained in [21, 22].
The present work is organized as follows. In Section 2, we deal with the blow-up solution of (1). Section 3 is devoted to the global solution of (1). As applications of the obtained results, some examples are presented in Section 4.
2. Blow-Up Solution
In this section, we discuss what interactions among the four nonlinear mechanisms of (1) result in the blow-up solution. Our main result in this section is the following theorem.
Theorem 1.
Let u(x,t) be a solution of problem (1). Assume that the following conditions (i)–(iv) are satisfied.
For s∈R+,(4)as-a′s≥0,f′s-fs≥0,asg′s′≥0.
For (x,s,t)∈∂Ω×R+×R+, (5)bsx,s,t-bx,s,t≤0,asbx,s,t-asbx,s,ts≤0,btx,s,t≥0.
Consider the following:(6)α≔minD¯au0∇·au0∇u0+fu0eu0g′u0>0.
Consider the following:(7)∫M0+∞asesds<+∞,M0≔maxD¯u0x.
Then u(x,t) blows up in a finite time T and (8)T≤1α∫M0+∞asesds,ux,t≤H-1αT-t,where (9)Hy≔∫y+∞asesds,y>0,and H-1 is the inverse function of H.
Proof.
Introduce an auxiliary function (10)Px,t≔-e-uut+α1au,and then we have (11)∇P=e-uut∇u-e-u∇ut-αa′a2∇u,(12)ΔP=-e-uut∇u2+2e-u∇u·∇ut+e-uutΔu-e-uΔut+2αa′2a3-αa′′a2∇u2-αa′a2Δu,(13)Pt=e-uut2-e-uutt-αa′a2ut=e-uut2-e-uag′Δu+a′g′∇u2+fg′t-αa′a2ut=e-uut2+ag′′g′2-a′g′e-uutΔu-ag′e-uΔut+a′g′′g′2-a′′g′e-uut∇u2-2a′g′e-u∇u·∇ut+fg′′g′2-f′g′e-uut-αa′a2ut.It follows from (12) and (13) that (14)ag′ΔP-Pt=a′′g′-ag′-a′g′′g′2e-uut∇u2+2ag′+2a′g′e-u∇u·∇ut+a′g′+ag′-ag′′g′2e-uutΔu+2αa′2a2g′-αa′′ag′∇u2-αa′ag′Δu-e-uut2+f′g′-fg′′g′2e-uut+αa′a2ut.The first equation of (1) implies (15)Δu=g′aut-a′a∇u2-fa.Inserting (15) into (14), we obtain (16)ag′ΔP-Pt=a′′g′-ag′-a′2ag′-a′g′e-uut∇u2+2ag′+2a′g′e-u∇u·∇ut+g′aag′′e-uut2+f′g′-a′fag′-fg′e-uut+3αa′2a2g′-αa′′ag′∇u2+αa′fa2g′.It follows from (11) that (17)∇ut=-eu∇P+ut∇u-αa′a2eu∇u.Substituting (17) into (16), we get (18)ag′ΔP+2a+a′g′∇u·∇P-Pt=a′′g′+ag′-a′2ag′+a′g′e-uut∇u2+αa′2a2g′-2αa′ag′-αa′′ag′∇u2+g′aag′′e-uut2+f′g′-a′fag′-fg′e-uut+αa′fa2g′.By (10), we have (19)ut=-euP+αeu1a.Now, we insert (19) into (18) to deduce (20)ag′ΔP+2a′+ag′∇u·∇P+ag′1+a′a+a′a′∇u2+fa′-faP-Pt=α1ag′a-a′∇u2+g′aag′′e-uut2+α1ag′f′-f.Assumptions (4) ensure that the right side in equality (20) is nonnegative; that is, (21)ag′ΔP+2a′+ag′∇u·∇P+ag′1+a′a+a′a′∇u2+fa′-faP-Pt≥0in Ω×0,T.Next, it follows from (1) and (10) that (22)∂P∂n=e-uut∂u∂n-e-u∂ut∂n-αa′a2∂u∂n=e-uutb-e-u∂u∂nt-αa′a2b=e-uutb-e-ubx,u,tt-αa′a2b=b-bue-uut-e-ubt-αa′a2b=b-bu-P+α1a-e-ubt-αa′a2b=bu-bP+αab-abua2-e-ubton ∂Ω×0,T.We note that (6) implies (23)maxΩ¯Px,0=maxΩ¯-e-u0u0t+α1au0=maxΩ¯-∇·au0∇u0+fu0eu0g′u0+α1au0=maxΩ¯1au0α-au0∇·au0∇u0+fu0eu0g′u0=0.There, (21)–(23), assumption (5) and the maximum principle [10] imply that the maximum of the function P in Ω¯×[0,T) is zero. In fact, if the function P takes a positive maximum at point (x0,t0)∈∂Ω×(0,T), then we have (24)Px0,t0>0,∂P∂nx0,t0>0.Using assumption (5) and the fact that P(x0,t0)>0, it follows from (22) that (25)∂P∂nx0,t0≤0,which contradicts the second inequality in (24). Hence, the maximum of the function P in Ω¯×[0,T) is zero. Now, we have (26)P≤0in Ω¯×0,T;that is, (27)aueuut≥α.At the point x0∈Ω¯, where u0(x0)=M0, integrating inequality (27) from 0 to t, we arrive at(28)∫0taueuutdt=∫M0ux0,tasesds≥αt.Inequality (28) and assumption (7) imply that u(x,t) blows up in finite time t=T. Now, we let t→T in (28) to deduce (29)T≤1α∫M0+∞asesds.For each fixed x∈Ω¯, integrating inequality (27) from t to s(0<t<s<T), we get (30)Hux,t≥Hux,t-Hux,s=∫ux,t+∞asesds-∫ux,s+∞asesds=∫ux,tux,sasesds=∫tsaueuutdt≥αs-t.In the above inequality, letting s→T, we obtain (31)Hux,t≥αT-t.We note that H is a strictly decreasing function. Hence, (32)ux,t≤H-1αT-t.The proof is complete.
3. Global Solution
In this section, we study what interactions among the four nonlinear mechanisms of (1) result in the global solution of (1). The main results of this section are formulated in the following theorem.
Theorem 2.
Let u(x,t) be a solution of problem (1). Assume that the following conditions (i)–(iv) are fulfilled.
For s∈R+,(33)as+a′s≤0,fs+f′s≤0,asg′s′≤0.
For (x,s,t)∈∂Ω×R+×R+,(34)bx,s,t+bsx,s,t≤0,asbx,s,t+asbx,s,ts≤0,btx,s,t≤0.
Consider the following:(35)β≔maxΩ¯au0∇·au0∇u0+fu0e-u0g′u0>0.
Consider the following:(36)∫m0+∞ase-sds=+∞,m0≔minΩ¯u0x.
Then u(x,t) must be a global solution and (37)ux,t≤K-1βt+Ku0x,where(38)Ky≔∫m0yase-sds,y≥m0,and K is the inverse function of K.
Proof.
We consider an auxiliary function (39)Qx,t≔-euut+β1au.Using the same reasoning process as that of (11)–(20), we obtain (40)ag′ΔQ+2a′-ag′∇u·∇Q+ag′1-a′a+a′a′∇u2+fa′+faQ-Qt=β1ag′a+a′∇u2+g′aag′′euut2+β1ag′f′+f.Assumptions (33) imply that the right side of (40) is nonpositive; that is, (41)ag′ΔQ+2a′-ag′∇u·∇Q+ag′1-a′a+a′a′∇u2+fa′+faQ-Qt≤0in Ω×0,T.By (1) and (39), we get (42)∂Q∂n=-euut∂u∂n-eu∂ut∂n-βa′a2∂u∂n=-euutb-eu∂u∂nt-βa′a2b=-euutb-eubx,u,tt-βa′a2b=-b+bueuut-eubt-βa′a2b=b+buQ-β1a-eubt-βa′a2b=b+buQ-βab+abua2-eubton ∂Ω×0,T.Assumption (35) implies (43)minΩ¯Qx,0=minΩ¯-eu0u0t+β1au0=minΩ¯-∇·au0∇u0+fu0e-u0g′u0+β1au0=minΩ¯1au0β-au0∇·au0∇u0+fu0e-u0g′u0=0.It follows from (41)–(43), (34), and the maximum principle that the minimum of Q in Ω¯×[0,T) is zero. Hence, we have the following inequality: (44)Q≥0in Ω¯×0,T;that is (45)aue-uut≤β.For each fixed x∈Ω¯, we integrate (45) from 0 to t to deduce (46)1β∫0taue-uutdt=1β∫u0xux,tase-sds≤t.Inequality (46) and assumption (36) imply that u must be a global solution. Furthermore, it follows from (45) that (47)Kux,t-Ku0x=∫m0ux,tase-sds-∫m0u0xase-sds=∫u0xux,tase-sds=∫0taue-uutdt≤βt.Since K is a strictly increasing function, we obtain (48)ux,t≤K-1βt+Ku0x.The proof is complete.
4. Applications
When a(u)=1 and b(x,u,t)=b(u), the results of Theorems 1–2 still hold. In this sense, our results extend and supplement those obtained in [21, 22].
In the following, we give a few examples to demonstrate the applications of Theorems 1–2.
Example 3.
Let u be a solution of the following problem:(49)u2+eu/2t=∇·eu/2∇u+euin Ω×0,T,∂u∂n=e3/4u-2+eu-2+tx2on ∂Ω×0,T,ux,0=1+x2in Ω¯,where Ω=x=(x1,x2,x3)∣|x|2=x12+x22+x32<1. We note (50)au=e2/u,gu=u2+eu/2,fu=eu,bx,u,t=e3/4u-2+eu-2+tx2.In order to calculate the constant α, we set (51)ω=x2,and then 0≤ω≤1 and (52)α=minΩ¯au0∇·au0∇u0+fu0eu0g′u0=minΩ¯26+2x2+e1/21+x2e1/21+x2+1=min0≤ω≤126+2ω+e1/21+ωe1/21+ω+1=28+e1+e.We can check that (4), (5), and (7) hold. It follows from Theorem 1 that u(x,t) blows up in a finite time T and (53)T≤1α∫M0+∞asesds=1+e28+e∫2+∞e-s/2ds=1+ee8+e,ux,t≤H-1αT-t=2ln1+e8+eT-t.
Example 4.
Let u be a solution of the following problem: (54)ueut=∇·e-u1+u∇u+e-uin Ω×0,T,∂u∂n=e-2u+e-u-tx2on ∂Ω×0,T,ux,0=1+x2in Ω¯,where Ω=x=(x1,x2,x3)∣|x|2=x12+x22+x32<1. Now we have (55)au=e-u1+u,gu=ueu,fu=e-u,bx,u,t=e-2u+e-u-tx2.Setting (56)ω=x2,we get 0≤ω≤1 and (57)β=maxΩ¯au0∇·au0∇u0+fu0e-u0g′u0=maxΩ¯16-2x2-3x4e21+x22+x24=max0≤ω≤116-2ω-3ω2e21+ω2+ω4=e-2.We can also check that (33), (34), and (36) hold. Hence, Theorem 2 implies u must be a global solution and (58)ux,t≤K-1βt+Ku0x=2+x2ete-2-1.
Remark 5.
When the functions a,b,f, and g are all exponential functions, Theorems 1 and 2 can be used. When not all of them are exponential functions, Theorems 1 and 2 can be used in some special cases.
Competing Interests
The author declares that there are no competing interests regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (no. 61473180).
QuittnerP.SoupletP.2007Basel, SwitzerlandBirkhäuserBirkhäuser Advanced TextsMarrasM.Vernier PiroS.On global existence and bounds for blow-up time in nonlinear parabolic problems with time dependent coefficients20132013535544PayneL. E.PhilippinG. A.SchaeferP. W.Blow-up phenomena for some nonlinear parabolic problems200869103495350210.1016/j.na.2007.09.0352-s2.0-51349123494ChenS. H.Global existence and nonexistence for some degenerate and quasilinear parabolic systems200824541112113610.1016/j.jde.2007.11.008MR24274072-s2.0-37549016066ChenS. H.MacDonaldK.Global and blowup solutions for general quasilinear parabolic systems201314142343310.1016/j.nonrwa.2012.07.006MR29698452-s2.0-84866357842DingJ. T.HuH. J.Blow-up and global solutions for a class of nonlinear reaction diffusion equations under Dirichlet boundary conditions20164332171817351974310.1016/j.jmaa.2015.08.046MR33987872-s2.0-84942200091DingJ. T.WangM.Blow-up solutions, global existence, and exponential decay estimates for second order parabolic problems201520151681610.1186/s13661-015-0429-yMR33956152-s2.0-84941906711ZhangL. L.Blow-up of solutions for a class of nonlinear parabolic equations200625447948610.4171/zaa/1303MR2285097AmannH.Quasilinear parabolic systems under nonlinear boundary conditions198692215319210.1007/BF00251255MR816618ZBL0596.350612-s2.0-0000324117SperbR.1981New York, NY, USAAcademic PressFriedmanA.1964Englewood Cliffs, NJ, USAPrentice-HallMR0181836HaradaJ.Non self-similar blow-up solutions to the heat equation with nonlinear boundary conditions20141022368310.1016/j.na.2014.01.0282-s2.0-84894594569HaradaJ.MiharaK.Blow-up rate for radially symmetric solutions of some parabolic equations with nonlinear boundary conditions201225351647166310.1016/j.jde.2012.05.008MR29273942-s2.0-84861951363LiuY.Lower bounds for the blow-up time in a non-local reaction diffusion problem under nonlinear boundary conditions2013573-492693110.1016/j.mcm.2012.10.002MR30112122-s2.0-84870480649LiuY.LuoS. G.YeY. H.Blow-up phenomena for a parabolic problem with a gradient nonlinearity under nonlinear boundary conditions20136581194119910.1016/j.camwa.2013.02.0142-s2.0-84876140058LiF. S.LiJ. L.Global existence and blow-up phenomena for nonlinear divergence form parabolic equations with inhomogeneous Neumann boundary conditions201238521005101410.1016/j.jmaa.2011.07.018MR28349062-s2.0-80052151177LiangF.Blow-up and global solutions for nonlinear reaction-diffusion equations with nonlinear boundary condition201121883993399910.1016/j.amc.2011.10.021MR28620722-s2.0-81855170054PayneL. E.PhilippinG. A.Vernier PiroS.Blow-up phenomena for a semilinear heat equation with nonlinear boundary condition, II201073497197810.1016/j.na.2010.04.0232-s2.0-77955317207DingJ.GaoX.LiS.Global existence and blow-up problems for reaction diffusion model with multiple nonlinearities2008343115916910.1016/j.jmaa.2008.01.038MR2409465ZBL1140.350102-s2.0-41449084963QiY. W.WangM. X.WangZ. J.Existence and non-existence of global solutions of diffusion systems with nonlinear boundary conditions200413461199121710.1017/s030821050000370xMR2107490ZhangH. L.Blow-up solutions and global solutions for nonlinear parabolic problems200869124567457410.1016/j.na.2007.11.0132-s2.0-55049090688ZhangL. L.ZhangN.LiL. X.Blow-up solutions and global existence for a kind of quasilinear reaction-diffusion equations201433324725810.4171/zaa/1509MR32295852-s2.0-84904915362