We obtain a new maximum principle of the periodic solutions when the corresponding impulsive equation is linear. If the nonlinear is quasi-monotonicity, we study the existence of the minimal and maximal periodic mild solutions for impulsive partial differential equations by using the perturbation method, the monotone iterative technique and the method of upper and lower solution. We give an example in last part to illustrate the main theorem.
Citation: Huanhuan Zhang, Jia Mu. Periodic problem for non-instantaneous impulsive partial differential equations[J]. AIMS Mathematics, 2022, 7(3): 3345-3359. doi: 10.3934/math.2022186
[1] | Xiongmei Fan, Sen Ming, Wei Han, Zikun Liang . Lifespan estimate of solution to the semilinear wave equation with damping term and mass term. AIMS Mathematics, 2023, 8(8): 17860-17889. doi: 10.3934/math.2023910 |
[2] | Jiangyan Yao, Sen Ming, Wei Han, Xiuqing Zhang . Blow-up of solutions to the coupled Tricomi equations with derivative type nonlinearities. AIMS Mathematics, 2022, 7(7): 12514-12535. doi: 10.3934/math.2022694 |
[3] | Salah Boulaaras, Abdelbaki Choucha, Bahri Cherif, Asma Alharbi, Mohamed Abdalla . Blow up of solutions for a system of two singular nonlocal viscoelastic equations with damping, general source terms and a wide class of relaxation functions. AIMS Mathematics, 2021, 6(5): 4664-4676. doi: 10.3934/math.2021274 |
[4] | Mohammad Kafini, Shadi Al-Omari . Local existence and lower bound of blow-up time to a Cauchy problem of a coupled nonlinear wave equations. AIMS Mathematics, 2021, 6(8): 9059-9074. doi: 10.3934/math.2021526 |
[5] | Hatice Taskesen . Qualitative results for a relativistic wave equation with multiplicative noise and damping terms. AIMS Mathematics, 2023, 8(7): 15232-15254. doi: 10.3934/math.2023778 |
[6] | Takiko Sasaki, Shu Takamatsu, Hiroyuki Takamura . The lifespan of classical solutions of one dimensional wave equations with semilinear terms of the spatial derivative. AIMS Mathematics, 2023, 8(11): 25477-25486. doi: 10.3934/math.20231300 |
[7] | Xiaoming Peng, Xiaoxiao Zheng, Yadong Shang . Lower bounds for the blow-up time to a nonlinear viscoelastic wave equation with strong damping. AIMS Mathematics, 2018, 3(4): 514-523. doi: 10.3934/Math.2018.4.514 |
[8] | Salim A. Messaoudi, Mohammad M. Al-Gharabli, Adel M. Al-Mahdi, Mohammed A. Al-Osta . A coupled system of Laplacian and bi-Laplacian equations with nonlinear dampings and source terms of variable-exponents nonlinearities: Existence, uniqueness, blow-up and a large-time asymptotic behavior. AIMS Mathematics, 2023, 8(4): 7933-7966. doi: 10.3934/math.2023400 |
[9] | Noufe H. Aljahdaly . Study tsunamis through approximate solution of damped geophysical Korteweg-de Vries equation. AIMS Mathematics, 2024, 9(5): 10926-10934. doi: 10.3934/math.2024534 |
[10] | Sen Ming, Jiayi Du, Yaxian Ma . The Cauchy problem for coupled system of the generalized Camassa-Holm equations. AIMS Mathematics, 2022, 7(8): 14738-14755. doi: 10.3934/math.2022810 |
We obtain a new maximum principle of the periodic solutions when the corresponding impulsive equation is linear. If the nonlinear is quasi-monotonicity, we study the existence of the minimal and maximal periodic mild solutions for impulsive partial differential equations by using the perturbation method, the monotone iterative technique and the method of upper and lower solution. We give an example in last part to illustrate the main theorem.
Nonlinear is one of the important scientific phenomena which widely exists in different physical systems. The nonlinear wave equation is a typical hyperbolic equation which possesses important physical background and theoretical significance. It can be used to explain and simulate many physical phenomena, such as the propagation of sound waves, small vibrations of elastic rods, and so on. The stability and fracture behavior of nonlinear problems are extremely important research topics. In the nonlinear systems, different initial exponential components can lead to different characteristics of their mathematical solutions. Due to the influence of nonlinear factors, waves will become steep during propagation until they break. Due to the existence of friction phenomena in reality, factors such as damping are inevitable, which could have a certain impact on the energy of physical systems. The coefficient in the damping term depending on time variable could affect the energy decay rate of the solution to the wave equation. Nonlinear external forces could affect the fracture behavior and lifespan estimation of the solution. Recently, many mathematicians pay attention to the blow-up and lifespan estimate of the solution to the nonlinear wave equation together with its generalized forms, which are related to the famous Strauss conjecture (see detailed illustrations in [1,2,3]).
In the present work, we study the Cauchy problem for a coupled system of semi-linear wave equations with damping terms and derivative type nonlinearities in Einstein-de Sitter space-time, which is shown as follows:
{utt−t−2k1Δu+μ1tut=|vt|p,(x,t)∈Rn×[1,∞),vtt−t−2k2Δv+μ2tvt=|ut|q,(x,t)∈Rn×[1,∞),u(x,1)=εf1(x),ut(x,1)=εg1(x),x∈Rn,v(x,1)=εf2(x),vt(x,1)=εg2(x),x∈Rn. | (1.1) |
Here, n≥1, p,q>1, ki∈[0,1), μi≥0, i=1,2. fi(x), gi(x)(i=1,2) are nonnegative functions which satisfy supp(fi,gi)⊂BR(0) for i=1,2, where BR(0)={x||x|≤R}, R≥2. ε is a small positive parameter describing the size of initial data.
To begin, we recall some known results for the classical wave equation
{utt−Δu=f(u,ut),x∈Rn,t>0,u(x,0)=εf(x),ut(x,0)=εg(x),x∈Rn. | (1.2) |
It is well-known that problem (1.2) with power nonlinearity f(u,ut)=|u|p possesses the Strauss exponent pS(n) (see [2]). In addition, pS(1)=∞. When n≥2, pS(n) is the positive solution to quadratic equation
−(n−1)p2+(n+1)p+2=0. |
More precisely, a solution to the Cauchy problem of the wave equation with small initial values blows up in finite time when 1<p≤pS(n). If p>pS(n), there exists a unique global (in time) solution. Lifespan estimate of solution to the initial boundary value problem of the variable coefficient wave equation with f(u,ut)=|u|p on exterior domain in two dimensions is obtained by employing the Kato lemma (see [4]). Both the blow-up result and existence of the global solution to problem (1.2) with f(u,ut)=|ut|p are discussed in [5,6]. Problem (1.2) with f(u,ut)=|ut|p possesses the Glassey index pG(n)=1+2n−1. When 1<p≤pG(n), a solution of the equation blows up in finite time. When p>pG(n), the problem admits global solution. Zhou et al. [7] established blow-up of the solution to the initial boundary value problem of the variable coefficient wave equation with derivative nonlinearity by solving inequalities of ordinary differential equations. Nonexistence of the global solution to problem (1.2) with f(u,ut)=|ut|p+|u|q is derived (see [8,9]). The proof is based on the Kato lemma and test function method. Existence of global solution to the wave equations is considered in [10,11,12]. Kitamura et al. [13] obtained the lifespan estimate of classical solution of nonlinear wave equation with spatial weight in one space dimension. Lai et al. [14] derived lifespan estimate of solution for 2-dimensional semi-linear wave equation in asymptotically Euclidean exterior domain. Formation of singularities of solution to the semi-linear wave equation in general dimensions are investigated (see [15,16,17,18,19]). Ming et al. [20] considered blow-up of solution to the semilinear Moore-Gibson-Thompson equation, which is structurally similar to the wave equation. Concerning the related study of hyperbolic type equations and other models, we refer readers to the related references [21,22,23,24] for more details.
Taking k1=k2=μ1=μ2=0 in problem (1.1), we obtain the following small initial value problem for the coupled system with classical wave equations:
{utt−Δu=|vt|p,(x,t)∈Rn×[0,∞),vtt−Δv=|ut|q,(x,t)∈Rn×[0,∞),(u,ut,v,vt)(x,0)=(εu0,εu1,εv0,εv1)(x),x∈Rn. | (1.3) |
We bear in mind that problem (1.2) is a special case of problem (1.3) when p=q and u=v. Now, we make brief review of the existing literature on the dynamic properties of solutions to the coupled system of wave equations. There are several research results on the blow-up of the solution to the problem (1.3). Nonexistence of the global solution to the coupled system is shown in [25,26]. Ikeda et al. [27] studied the behavior of solutions to the problem in different coupling cases with a variety of nonlinear terms. Series blow-up results and upper bound lifespan estimates of solutions to the problem are obtained. Based on the above results, it can be seen that there is a (p,q) critical curve, which describes the blow-up and existence of global solutions, namely,
Ω(n,p,q):=max(Λ(n,p,q),Λ(n,q,p))=0, |
where
Λ(n,p,q)=p+1pq−1−n−12, | (1.4) |
and p,q>1. When Ω(n,p,q)<0, there are global solutions to the Cauchy problem (1.3). When Ω(n,p,q)≥0, no matter how small the initial values are, the solutions always blow up in finite time. Upper bound estimate of the lifespan of solutions satisfies:
T(ε)≤{Cε−Ω(n,p,q)−1,Ω(n,p,q)>0,exp(Cε−(pq−1)),Ω(n,p,q)=0,p≠q,exp(Cε−(p−1)),Ω(n,p,q)=0,p=q. | (1.5) |
Kubo et al. [28] proved the existence of the global solution to the problem in the three-dimensional and radial symmetry case.
In problem (1.1), let k1=k2=−m, m>0, μ1=μ2=0. We arrive at the following small initial value problem for coupled Tricomi equations:
{utt−t2mΔu=|vt|p,(x,t)∈Rn×[0,∞),vtt−t2mΔv=|ut|q,(x,t)∈Rn×[0,∞),(u,ut,v,vt)(x,0)=(εu0,εu1,εv0,εv1)(x),x∈Rn. | (1.6) |
When p=q, and u=v in problem (1.6), we obtain the Cauchy problem for the single Tricomi equation, namely,
{utt−t2mΔu=f(u,ut),x∈Rn,t>0,u(x,0)=εf(x),ut(x,0)=εg(x),x∈Rn. | (1.7) |
The Tricomi equation arises in gas dynamic problems, which is associated with the gas flows with nearly sonic speed describing the transition from subsonic flow to supersonic flow. Let us recall some related investigation of the Cauchy problem for the Tricomi equation. Problem (1.7) with f(u,ut)=|u|p possesses the Strauss critical exponent pS(m,n), which is the biggest root of quadratic equation
−[n−1+(1−2m+2)]p2+[n+1−3(1−2m+2)]p+2=0. |
Applying the iteration argument, Lin et al. [29] derived blow-up dynamics and lifespan estimates of solutions to problem (1.7) with f(u,ut)=|u|p in the subcritical and critical cases. Lifespan estimate of solution to problem (1.7) with f(u,ut)=|ut|p is deduced by applying the test function technique (Φ(x,t)=−t−2m∂t(η2p′M(t)l(t))ϕ(x)) (see [30]). Chen et al. [31] investigated the blow-up result and upper bound lifespan estimate of solution to problem (1.7) with f(u,ut)=|ut|p+|u|q, where the iteration method is used. It is worth it to mention that pS(0,n) coincides with the Strauss critical exponent pS(n) of the classical wave equation. We refer the readers to [31,32,33,34] for the relevant illustrations.
Recently, the study of the Cauchy problem of the semi-linear Tricomi equation with damping term and mass term
{utt−t2mΔu+g(ut)+h(u)=f(u,ut),x∈Rn,t>0,u(x,0)=εu0(x),ut(x,0)=εu1(x),x∈Rn | (1.8) |
attracts more attention, where g(ut) is the damping term and h(u) is the mass term. For the case of m=0 in problem (1.8), using the test function technique, Ikeda et al. [35] studied formation of singularity for the solution to the semi-linear wave equation with damping term and mass term in the subcritical and critical cases. Lai et al. [36] showed the blow-up result and lifespan estimate of solutions to problem (1.8) with space dependent damping term, potential term, and power nonlinearity by employing the test function technique. Ming et al. [37] discussed the Cauchy problem of the semi-linear wave equation with scattering time dependent damping term and divergence form nonlinearities. Upper bound lifespan estimate of solution in the subcritical and critical cases is deduced by making use of the rescaled test function approach and iteration method. Hamouda et al. [38,39] verified the blow-up result of the solution to problem (1.8) with −1<m<0, g(ut)=μt−1ut, h(u)=0, and derivative nonlinearity in the case of t>1. That is,
{utt−t−2kΔu+μtut=|ut|p,x∈Rn,t>1,u(x,0)=εu0(x),ut(x,0)=εu1(x),x∈Rn, | (1.9) |
where k∈(0,1). Under certain assumptions, Hamouda et al. [39] applied the test function method to derive the functional differential inequality of solution in the generalized Einstein-de Sitter space-time. Blow-up index p=pE(n,k,μ)=1+2(1−k)(n−1)+k+μ is presented. When 1<p≤pE(n,k,μ), a solution of the equation blows up in finite time. In the special case μ=2, Hamouda et al. [38] derived the blow-up result for the semi-linear wave equation in generalized Einstein-de Sitter space-time with derivative type nonlinearity by taking advantage of the integral representation formula for the solution and Yagdjian integral transform technique. Palmieri [40] deduces formation of singularity for solution to problem (1.8) with m>−1, g(ut)=μt−1ut, and h(u)=ν2t−2u for t>1 by applying the iteration method, where the nonlinear term is power type. Similar to the approach utilized in [39], Hamouda et al. [41] considered formation of singularity of the solution to the corresponding small initial boundary value problem (1.8) with m>0, g(ut)=μt−1ut, and h(u)=0 on exterior domain. Blow-up dynamics of the small initial value problem of the Tricomi equation with the general damping term, and mass term are presented (see [42]). The related results of generalized variable coefficients wave equations are obtained in [43,44,45,46]. We also refer the interested readers to the references [47,48,49,50,51] for more details.
Now, we turn to problem (1.6). Under certain hypothetical conditions, there is a (p,q) curve (which has not been proved to be critical or not) regarding the formation of singularity of solutions to system (1.6), that is,
Ω(n,m,p,q):=max(Λ(n,m,p,q),Λ(n,m,q,p))=0, |
where Λ(n,m,p,q)=p+1pq−1−n(m+1)−2m−12. When Ω(n,m,p,q)≥0, the solutions (u,v) of the problem blow up in finite time T(ε), which satisfies:
T(ε)≤{Cε−Ω(n,m,p,q),Ω(n,m,p,q)>0,exp(Cε−(pq−1)),Ω(n,m,p,q)=0,p≠q,exp(Cε−(p−1)),Ω(n,m,p,q)=0,p=q. |
For the mixed form of nonlinear terms in problem (1.6), relevant research of upper bound lifespan estimates of solutions can be referred to [52,53]. The related study of the blow-up of solutions to problem (1.6) containing scale invariant damping terms and mass terms is shown in [54].
In problem (1.1), letting k1=k2=0 and modifying the time coefficient of damping terms, we obtain the small initial value problem of coupled system with scale invariant damping terms as follows:
{utt−Δu+μ11+tut=|vt|p,(x,t)∈Rn×[0,∞),vtt−Δv+μ21+tvt=|ut|q,(x,t)∈Rn×[0,∞),(u,ut,v,vt)(x,0)=(εu0,εu1,εv0,εv1)(x),x∈Rn. | (1.10) |
Chen and Palmieri [55] first show the critical curve of system (1.10) in the effective case by applying the contraction mapping technique and test function method. By constructing appropriate functions, Hamouda et al. [56,57] obtained the coupled integral inequalities. Blow-up of solutions to problem (1.10) is established by taking advantage of the iterative method. In addition, lifespan estimation T(ε) of solutions satisfies:
T(ε)≤{Cε−(Ω(n,μ1,μ2,p,q))−1,Ω(n,μ1,μ2,p,q)>0,exp(Cε−(pq−1)),Ω(n,μ1,μ2,p,q)=0, p≠q,exp(Cε−min(pq−1p+1,pq−1q+1)),Ω(n,μ1,μ2,p,q)=0,p=q, |
where Ω(n,μ1,μ2,p,q)=max(Λ(n+μ1,p,q),Λ(n+μ2,q,p)) (here, Λ(n,p,q) is defined by (1.4)). Blow-up region of solutions to the problem is presented. In problem (1.10), by replacing the coefficient of damping terms μi1+t(i=1,2) with general nonnegative continuous and integrable functions bi(t)(i=1,2), we obtain a general weakly coupled system of semi-linear damped wave equations. By applying the multiplier technique, Palmieri and Takamura [58] proved that there is a (p,q) critical curve for solutions to the problem. Compared with that of problem (1.3), the (p,q) critical curve does not change. Upper bound estimate for the lifespan of solution is the same as (1.5). The studies of the blow-up and existence of global solutions to coupled systems with other types of damping terms and nonlinear terms are presented in [27,59,60,61].
Since the wave speeds change with time in problem (1.1), we transform the damping terms to obtain scale invariant dampings, namely, u1(x,τ1)=u(x,t), v1(x,τ2)=v(x,t), where
τi=ϕki(t):=t1−ki1−ki,i=1,2. | (1.11) |
Then, u1(x,τ1) and v1(x,τ2) satisfy the following equations, respectively:
(u1)τ1τ1−Δu1+μ1−k1(1−k1)τ1(u1)τ1=(1−k1)2k11−k1(1−k2)k2pk2−1τ2k11−k11τk2pk2−12|(v1)τ2|p, |
(v1)τ2τ2−Δv1+μ2−k2(1−k2)τ2(v1)τ2=(1−k2)2k21−k2(1−k1)k1qk1−1τ2k21−k22τk1qk1−11|(u1)τ1|q. |
In fact, we make the transformation of time variable t in (1.11), so that the ranges of compact support of u(x,t) and v(x,t) are
supp(u,v)⊂{(x,t)∈Rn×[1,∞)||x|≤ϕki(t)+R},i=1,2. |
Inspired by the works in [38,39,41,42,54,56,62], our target in this paper is to investigate blow-up and lifespan estimates of solutions to the weakly coupled system of the semi-linear wave equation with scale invariant dampings and derivative nonlinearities in generalized Einstein-de Sitter space-time. We observe that nonexistence of the global solution to the Cauchy problem of the semi-linear wave equation with derivative nonlinearity in the generalized Einstein-de Sitter space-time is derived in [38]. The proof is based on the integral representation formula for the solution to the corresponding linear problem in the one-dimensional case, which is determined through the Yagdjian integral transform approach. It is worth noticing that Hassen et al. [42] verified the nonexistence of global solution to the generalized Tricomi equation with scale invariant damping, mass term, and derivative nonlinearity in the subcritical and critical cases. Upper bound lifespan estimates of solution to the problem are obtained by making use of the method of solving inequalities of ordinary differential equations. Employing the similar method in [42], Hamouda et al. [41] established the blow-up dynamic and lifespan estimate of solutions to the initial boundary problem for the generalized Tricomi equation with scale invariant damping and derivative nonlinearity on exterior domain. For problem (1.9), Hamouda et al. [39] exploited the test function method to establish functional differential inequality of the solution in the generalized Einstein-de Sitter space-time. Blow-up region and upper bound lifespan estimates of solutions to problem (1.9) are established. Hamouda et al. [56] obtained the blow-up of solutions to the Cauchy problem of the semi-linear wave equation with scale invariant damping terms (namely, problem (1.10)) by deriving the coupled integral inequalities and using the iterative approach. Improvement on the blow-up results of solutions to the weakly coupled system of wave equations with scale invariant damping terms, mass terms, and derivative nonlinearities are considered (see [62]). Taking advantage of the test function technique, Hassen et al. [54] investigated the blow-up phenomenons of solutions to the weakly coupled system of the generalized Tricomi equation with scale invariant dampings, mass terms, and derivative type nonlinear terms. From our observation, there is no related results for blow-up and lifespan estimate of solutions to problem (1.1) with scale invariant damping terms and derivative nonlinearities in the generalized Einstein-de Sitter space-time. The main difficulties in this work are to establish appropriate iterative frameworks in the coupled system case and lower bounds for the constructed functionals, which are used to derive the finite time blow-up and upper bound lifespan estimate of solutions. We extend the problem studied in [39] to problem (1.1) by utilizing the test function method, which is different from the constructed functions in [54,62] (see Theorem 1.1). The results in Theorem 1.1 coincide with the results in [39] when p=q and k1=k2. As a consequence, the results in [39] are a special case of the results in Theorem 1.1. Compared with blow-up results in the work [56], when ki=0,i=1,2 in problem (1.1) (that is, σi=μi), blow-up region Ω(n,σ1,σ2,p,q) of problem (1.1) is same as that in [56]. The blow-up region of solutions obtained in this paper varies due to the influence of the coefficients k1,k2. To the best of our knowledge, the results in Theorem 1.1 are new.
For the problem (1.1), first of all, we present the definition of its energy solutions.
Definition 1.1. Let (fi,gi)∈H1(Rn)×L2(Rn), i=1,2. The solutions (u,v) are called energy solutions of problem (1.1) on [1,T), if
{(u,v)∈C([1,T),H1(Rn))∩C1([1,T),L2(Rn)),ut∈Lqloc((1,T)×Rn), vt∈Lploc((1,T)×Rn), |
and the following two equalities are satisfied:
∫Rnut(x,t)Φ(x,t)dx−ε∫Rng1(x)Φ(x,1)dx−∫t1∫Rnut(x,s)Φt(x,s)dxds+∫t1s−2k1∫Rn∇u(x,s)∇Φ(x,s)dxds+∫t1∫Rnμ1sut(x,s)Φ(x,s)dxds=∫t1∫Rn|vt(x,s)|pΦ(x,s)dxds, |
∫Rnvt(x,t)˜Φ(x,t)dx−ε∫Rng2(x)˜Φ(x,1)dx−∫t1∫Rnvt(x,s)˜Φt(x,s)dxds+∫t1s−2k2∫Rn∇v(x,s)∇˜Φ(x,s)dxds+∫t1∫Rnμ2svt(x,s)˜Φ(x,s)dxds=∫t1∫Rn|ut(x,s)|q˜Φ(x,s)dxds, |
where Φ,˜Φ∈C∞0(Rn×[1,T)).
The main result of this paper is stated below.
Theorem 1.1. Let Ω(n,σ1,σ2,p,q)=max(Λ(n+σ1,p,q),Λ(n+σ2,q,p))≥0, where Λ(n,p,q) is defined by (1.4), σi=ki+μi, i=1,2. Suppose that fi∈H1(Rn), gi∈L2(Rn) are non-negative functions which are not identically zero. If (u,v) are the energy solutions of problem (1.1) which satisfy supp(u,v)⊂{(x,t)∈Rn×[1,∞)||x|≤ϕki(t)+R}, then there exists a positive constant ε0(n,k1,k2,μ1,μ2,f1,f2,g1,g2,T3) (T3 is a positive constant), such that for all ε∈(0,ε0], the upper bound estimate of lifespan T(ε) of (u,v) satisfies:
T(ε)≤{Cεpq−1A1−B1,Ω(n,σ1,σ2,p,q)=Λ(n+σ1,p,q)≥0,σ1≥σ2,Cεpq−1A2−B2,Ω(n,σ1,σ2,p,q)=Λ(n+σ2,q,p)≥0,σ1≤σ2,exp(Cε−min(pq−1p+1,pq−1q+1)),Λ(n+σ1,p,q)=Λ(n+σ2,q,p)=0, | (1.12) |
where the constant C is independent of ε, A1=n−12(pq−1)+σ22pq+σ12p, B1=σ22p+σ12+p+1, A2=n−12(pq−1)+σ12pq+σ22q, B2=σ12q+σ22+q+1.
Remark 1.1. It is worth it to mention that the upper bound lifespan estimate (1.12) of solutions to problem (1.1) in Theorem 1.1 coincides with the results in [39] when p=q and k1=k2. In addition, when ki=0,i=1,2 in problem (1.1) (that is, σi=μi), blow-up region Ω(n,σ1,σ2,p,q) of problem (1.1) is the same as that in [56]. Therefore, the results in [39,56] could be viewed as special cases of the results in Theorem 1.1. Moreover, the results in Theorem 1.1 complement the case that propagation speeds of generalized coupled wave equations with dampings and derivative type nonlinearities contain power functions with negative exponents.
Suppose that ρi(t)(i=1,2) satisfies the following differential equation:
d2dt2ρi(t)−t−2kiρi(t)−ddt(μitρi(t))=0,t≥1,i=1,2. |
According to [34,39], ρi(t) can be expressed as
ρi(t)=t1+μi2Kμi−12(1−ki)(ϕki(t)),t≥1,i=1,2, | (2.1) |
where ϕki(t) are defined by (1.11),
Kν(t)=∫∞0exp(−tcoshξ)cosh(νξ)dξ,ν∈R. |
It holds that
limt→∞tkiρ′i(t)ρi(t)=−1,i=1,2, | (2.2) |
ρ1(t)eϕk1(t)≥C−10tk1+μ12,ρ−p2(t)e−pϕk1(t)≥C−p0t−p(k1+μ1)2, | (2.3) |
where C0 is a positive constant. Similarly, we acquire
ρ2(t)eϕk2(t)≥˜C−10tk2+μ22,ρ−q1(t)e−qϕk2(t)≥˜C−q0t−q(k2+μ2)2, |
where ˜C0 is a positive constant.
We define the function φi(x)(i=1,2) as follows:
φi(x)={ex+e−x,n=1,∫Sn−1ex⋅ωdω,n≥2. |
Therefore, we obtain Δφi(x)=φi(x). Let Ψi(x,t)=φi(x)ρi(t)(i=1,2). It holds that ΔΨi(t)=Ψi(x,t) and
d2dt2Ψi(x,t)−t−2kiΔΨi(x,t)−ddt(μitΨi(x,t))=0. |
Lemma 2.1. [17] Let r>1. It holds that
∫{|x|≤ϕki(t)+R}Ψri(x,t)dx≤C1ρri(t)erϕki(t)(1+ϕki(t))(2−r)(n−1)2,t≥1,i=1,2, | (2.4) |
where C1 is a positive number C1=C(n,μ1,μ2,k1,k2,p,q,R,r).
Now, we define the following functions related to the solutions u(x,t) and v(x,t) of problem (1.1):
G1(t)=∫Rnu(x,t)Ψ1(x,t)dx,˜G1(t)=∫Rnut(x,t)Ψ1(x,t)dx, |
G2(t)=∫Rnv(x,t)Ψ2(x,t)dx,˜G2(t)=∫Rnvt(x,t)Ψ2(x,t)dx. |
In the lemma below, we present lower bound estimates for Gi(t) and ˜Gi(t).
Lemma 2.2. [39] Suppose that (u,v) are solutions to problem (1.1). The initial value functions fi,gi(i=1,2) satisfy the assumptions of Theorem 1.1. Then, there exists a constant T0=T0(k1,k2,μ1,μ2)>2 such that
Gi(t)≥CGiεtki,t≥T0,i=1,2, |
where CGi=CGi(f1,g1,f2,g2,n,k1,k2,μ1,μ2,R).
Lemma 2.3. [39] Suppose that (u,v) are solutions to problem (1.1). The initial value functions fi,gi(i=1,2) satisfy the assumptions of Theorem 1.1. Then, there exists T1=T1(k1,k2,μ1,μ2)>T0 such that
˜Gi(t)≥C˜Giε,t≥T1,i=1,2, |
where C˜Gi=C˜Gi(f1,g1,f2,g2,n,k1,k2,μ1,μ2,R). It holds that
˜G1(t)+[μ1t−ρ′1(t)ρ1(t)]G1(t)=∫t1∫Rn|vt(x,s)|pΨ1(x,s)dxds+εC(f1,g1), | (2.5) |
˜G′1(t)+[μ1t−ρ′1(t)ρ1(t)]˜G1(t)=t−2k1G1(t)+∫Rn|vt(x,t)|pΨ1(x,t)dx, | (2.6) |
where ρ1(t) is defined by (2.1).
Remark 2.1. For ˜G2(t) and ˜G′2(t), we acquire similar forms of (2.5) and (2.6) with suitable modification, which are omitted here for simplification.
Let α be a positive constant. Multiplying both sides of (2.5) by αρ′1(t)ρ1(t) and using (2.6), we arrive at
˜G′1(t)+(μ1t−(1−α)ρ′1(t)ρ1(t))˜G1(t)=−εαρ′1(t)ρ1(t)C(f1,g1)+(t−2k1+αρ′1(t)ρ1(t)(μ1t−ρ′1(t)ρ1(t)))G1(t)+∫Rn|vt(x,t)|pΨ1(x,t)dx−αρ′1(t)ρ1(t)∫t1∫Rn|vt(x,s)|pΨ1(x,s)dxds,t≥1. |
Noticing the condition (2.2), we choose T2≥T1 such that
˜G′1(t)+(μ1t−(1−α)ρ′1(t)ρ1(t))˜G1(t)≥εαt−k12C(f1,g1)+(1−4α)t−2k1G1(t)+∫Rn|vt(x,t)|pΨ1(x,t)dx+αt−k12∫t1∫Rn|vt(x,s)|pΨ1(x,s)dxds,t≥T2. | (3.1) |
We set α∈(17,14). Using Lemma 2.2, formula (3.1) can be rewritten as
˜G′1(t)+(μ1t−(1−α)ρ′1(t)ρ1(t))˜G1(t)≥εαt−k12C(f1,g1)+∫Rn|vt(x,t)|pΨ1(x,t)dx+αt−k12∫t1∫Rn|vt(x,s)|pΨ1(x,s)dxds,t≥T2. | (3.2) |
In the following, we set
L1(t)=116∫tT3∫Rn|vt(x,s)|pΨ1(x,s)dxds+εC2,L2(t)=116∫tT3∫Rn|ut(x,s)|qΨ2(x,s)dxds+εC2, | (3.3) |
where C2=min(αC(f1,g1)4(1+α),αC(f2,g2)4(1+α),C˜G1,C˜G2).
Let Fi(t)=˜G1(t)−L1(t)(i=1,2). Making use of (3.2), we know
F′1(t)+(μ1t−(1−α)ρ′1(t)ρ1(t))F1(t)≥1516∫Rn|vt(x,t)|pΨ1(x,t)dx+(α2−116(μ1t1−k1−(1+α)tk1ρ′1(t)ρ1(t)))+(α2C(f1,g1)−C2(μ1t1−k1−(1+α)tk1ρ′1(t)ρ1(t)))εt−k1,t≥T3, |
where T3≥T2. When t≥T3, we arrive at the following inequalities:
α2−116tk1(μ1t−(1+α)ρ′1(t)ρ1(t))≥0, |
α2C(f1,g1)−C2tk1(μ1t−(1+α)ρ′1(t)ρ1(t))≥0. |
Thus, we obtain
F′1(t)+(μ1t−(1−α)ρ′1(t)ρ1(t))F1(t)≥0,t≥T3. | (3.4) |
Multiplying both sides of (3.4) by tμ1ρ1+α1(t) and integrating over (T3,t), we acquire
F1(t)≥F1(T3)Tμ13ρ1+α1(T3)ρ1+α1(t)tμ1,t≥T3. |
Combining the definition of C2 and applying
F1(T3)=˜G1(T3)−L1(T3)≥˜G1(T3)−C˜G1ε≥0, |
we conclude F1(t)≥0. It follows that
˜G1(t)≥L1(t),t≥T3. |
Similarly, we have
˜G2(t)≥L2(t),t≥T3. | (3.5) |
Making use of the H¨older inequality and (2.4) yields
∫Rn|vt(x,t)|pΨ1(x,t)dx≥(∫{|x|≤ϕk1(t)+R}Ψpp−12(x,t)Ψ−1p−11(x,t)dx)−(p−1)(∫Rnvt(x,t)Ψ2(x,t)dx)p≥C1−p1ρ1(t)ρ−p2(t)e−(p−1)ϕk1(t)(1+ϕk1(t))−(p−1)(n−1)2˜Gp2(t). | (3.6) |
Substituting (2.4) into (3.6) leads to
∫Rn|vt(x,t)|pΨ1(x,t)dx≥C−p−10C1−p1t−(k1+μ1)(p−1)2(1+ϕk1(t))−(p−1)(n−1)2˜Gp2(t). |
Combining (3.3) and (3.5), we have
L′1(t)≥C3t−(p−1)(k1+μ1+n−1)2Lp2(t),t≥T3, | (3.7) |
where C3=C3(p,q). Similarly, we obtain
L′2(t)≥C3t−(q−1)(k2+μ2+n−1)2Lq1(t),t≥T3. | (3.8) |
Integrating (3.7) and (3.8) over (T3,t) gives rise to
L1(t)≥C4∫tT3(s+T3)−(p−1)(k1+μ1+n−1)2Lp2(s)ds+εC3,t≥T3, | (3.9) |
L2(t)≥C4∫tT3(s+T3)−(q−1)(k2+μ2+n−1)2Lq1(s)ds+εC3,t≥T3. | (3.10) |
In the following, we define
σi=ki+μi,(i=1,2). | (3.11) |
We note the constant T3>2. Using the definition of σi(i=1,2), (3.9) and (3.10) can be rewritten as follows:
L1(t)≥C4∫tT3(s+T3)−n−12(p−1)+σ12−σ12pLp2(s)ds+C3ε,t≥T3, | (3.12) |
L2(t)≥C4∫tT3(s+T3)−n−12(q−1)+σ22−σ22qLq1(s)ds+C3ε,t≥T3. | (3.13) |
Suppose the condition Ω(n,σ1,σ2,p,q)=Λ(n+σ1,p,q)>0. In the following, we assume
L1(t)≥αj(t+T3)−βj(t−T3)γj,t≥T3, | (3.14) |
where {αj}j∈N, {βj}j∈N, {γj}j∈N are nonnegative sequences of real numbers.
When j=0, we set α0=C2ε, β0=0, γ0=0. We intend to prove that (3.14) is true for all j≥0. Substituting (3.14) into (3.13), we get
L2(t)≥C4∫tT3(s+T3)−n−12(q−1)+σ22−σ22qLq1(s)ds≥C4αqj∫tT3(s+T3)−n−12(q−1)+σ22−σ22q−qβj(s−T3)qγjds≥C4αqj(t+T3)−n−12(q−1)−σ22q−qβj∫tT3(s−T3)qγj+σ22ds=C4αqj(qγj+σ22+1)−1(t+T3)−n−12(q−1)−σ22q−qβj(t−T3)qγj+σ22+1, |
where t≥T3. Substituting the above lower bound estimate for L2(t) into (3.12) yields
L1(t)≥C4∫tT3(s+T3)−n−12(p−1)+σ12−σ12pLp2(s)ds≥Cp+14αpqj(qγj+σ22+1)p∫tT3(s+T3)−n−12(pq−1)−σ22pq−σ12p−pqβj×(s−T3)pqγj+σ22p+σ12+pds≥Cp+14αpqj(qγj+σ22+1)p(t+T3)−n−12(pq−1)−σ22pq−σ12p−pqβj×∫tT3(s−T3)pqγj+σ22p+σ12+pds=Cp+14αpqj(qγj+σ22+1)p(pqγj+σ22p+σ12+p+1)×(t+T3)−n−12(pq−1)−σ22pq−σ12p−pqβj×(t−T3)pqγj+σ22p+σ12+p+1, t≥T3. |
Making use of (3.14), we have
αj+1=Cp+14(σ22+1+qγj)−p(σ22p+σ12+p+1+pqγj)αpqj, | (3.15) |
βj+1=n−12(pq−1)+σ22pq+σ12p+pqβj,γj+1=σ22p+σ12+p+1+pqγj. | (3.16) |
Let A1=n−12(pq−1)+σ22pq+σ12p, β0=0. Employing (3.16) leads to
βj=A1+pqβj−1=A1j−1∑k=0(pq)k+(pq)jβ0. | (3.17) |
In a similar way, we have
γj=B1(pq)j−1pq, | (3.18) |
where B1=σ22p+σ12+p+1.
According to αj, γj≤B1pq−1(pq)j, and (3.15), we acquire
αj=Cp+14(σ22+1+qγj−1)−p(σ22p+σ12+p+1+pqγj−1)−1αpqj−1≥Cp+14(σ22p+σ12+p+1+pqγj−1)−(p+1)αpqj−1≥Cp+14(B1pq−1)−(p+1)⏟C5(pq)−(p+1)jαpqj−1=C5(pq)−(p+1)jαpqj−1. | (3.19) |
Taking the logarithm of both sides of (3.19), we derive
logαj≥pqlogαj−1−jlog(pq)p+1+logC5≥(pq)2logαj−2−(j+(j−1)pq)log(pq)p+1+(1+pq)logC5≥(pq)jlogα0−(j−1∑k=0(j−k)(pq)k)log(pq)p+1+(j−1∑k=0(pq)k)logC5=(pq)j(logα0−pq(pq−1)2log(pq)p+1+logC5pq−1)+(j+1)log(pq)p+1pq−1+log(pq)p+1(pq−1)2−logC5pq−1⏟≥0. |
Let j≥˜j=max{0,logC4log(pq)p+1−pqpq−1}. It follows that
logαj≥(pq)jlog(C6ε), | (3.20) |
where C6=C3(pq)−pq(p+1)(pq−1)2C1pq−15. Finally, combining the formulas (3.14), (3.17), (3.18), and (3.20), we arrive at
L1(t)≥αj(t+T3)−A1(pq)j−1pq−1(t−T3)B1(pq)j−1pq−1≥exp((pq)jlog(C6ε))(t+T3)−A1(pq)j−1pq−1(t−T3)B1(pq)j−1pq−1=exp((pq)j(log(C6ε)−A1pq−1log(t+T3)+B1pq−1log(t−T3)))×(t+T3)A1pq−1(t−T3)−B1pq−1,j≥˜j,t≥T3. | (3.21) |
Let t≥3T3. Hence, 2(t−T3)≥t+T3. According to (3.21), we have
L1(t)≥exp((pq)jlog(2−B1pq−1C6ε(t+T3)B1−A1pq−1))×(t+T3)A1pq−1(t−T3)−B1pq−1,j≥˜j,t≥T3. | (3.22) |
Combining the definition of Λ(n+σ1,p,q) and its nonnegativity, we know
B1−A1pq−1≥−σ2p(q−1)+σ1(p−1)2(pq−1)+p+1pq−1−n−12≥σ2p(1−q)+σ1(1−p)2(pq−1)+σ12=p(q−1)(σ1−σ2)2(pq−1). | (3.23) |
Suppose σ1≥σ2. Then, we have B1−A1pq−1≥0. Because of t+T3≥4T3, we choose ε0=ε0(n,k1,k2,μ1,μ2,f1,f2,g1,g2,T3)>0 which is small enough such that t+T3≥4T3. Thus, we have
(2−B1pq−1C6ε0)−(B1−A1pq−1)−1≥4T3. |
Hence, for all ε∈(0,ε0] and t>(2−B1pq−1C6ε)−(B1−A1pq−1)−1≥4T3, we obtain 2−B1pq−1C6εtB1−A1pq−1>1. Therefore, sending j→∞ in (3.22), we arrive at the blow-up result of L1(t) in finite time. Blow-up of L2(t) and L1(t) occurs simultaneously. Let 2−B1pq−1C6εtB1−A1pq−1≤1. Upper bound estimate for the lifespan of L1(t) and L2(t) is
T(ε)≤Cεpq−1A1−B1, | (3.24) |
where A1=n−12(pq−1)+σ22pq+σ12p, B1=σ22p+σ12+p+1, σ1, and σ2 are defined by (3.11), and A1−B1<0 if σ1≥σ2.
On the other hand, let Ω(n,σ1,σ2,p,q)=Λ(n+σ2,q,p)>0, and using the same method, we see that L1(t) and L2(t) blow up in finite time. The upper bound estimate of their lifespan satisfies:
T(ε)≤Cεpq−1A2−B2, | (3.25) |
where A2=n−12(pq−1)+σ12pq+σ22q, B2=σ12q+σ22+q+1, σ1, and σ2 are defined by (3.11), and A2−B2<0 if σ2≥σ1.
Without loss of generality, we may assume Ω(n,σ1,σ2,p,q)=Λ(n+σ1,p,q)=0>Λ(n+σ2,q,p). Then, we deduce from (3.23) that
B1−A1pq−1=σ2p(1−q)+σ1(1−p)2(pq−1)+σ12=p(q−1)(σ1−σ2)2(pq−1). |
Because of pq−p>0, we get B1−A1pq−1≥0 if σ1≥σ2. Similar to the derivations in Subsection 3.1, we obtain the same estimate as in (3.24), namely,
T(ε)≤Cεpq−1A1−B1. |
On the other hand, letting Ω(n,σ1,σ2,p,q)=Λ(n+σ2,q,p)=0>Λ(n+σ1,p,q) and using the similar method, we see that the upper bound estimate for the lifespan of L1(t) and L2(t) is same as (3.25) if σ2≥σ1, that is,
T(ε)≤Cεpq−1A2−B2. |
Using the condition Λ(n+σ1,p,q)=Λ(n+σ2,q,p)=0, we rewrite (3.12) and (3.13) as
L1(t)≥C4∫tT3(s+T3)−p2−1pq−1Lp2(s)ds+C3ε≥C4∫tT3(s+T3)−(p+1)Lp2(s)ds+C3ε,t≥T3, | (3.26) |
L2(t)≥C4∫tT3(s+T3)−q2−1pq−1Lq1(s)ds+C3ε≥C4∫tT3(s+T3)−(q+1)Lq1(s)ds+C3ε,t≥T3. | (3.27) |
In the following, we estimate the lower bound of integrals in (3.26) and (3.27).
Suppose
L1(t)≥Ej(t+T3)−pαj(logtT3)βj,t≥T3, | (3.28) |
where {Ej}j∈N, {αj}j∈N, {βj}j∈N are nonnegative real numbers sequences. When j=0, we derive E0=C3ε, α0=β0=0. Noticing t+T3≤2t, for all t≥T3, and substituting (3.28) into (3.27), we conclude
L2(t)≥C4(t+T3)−q∫tT3(s+T3)−1Lq1(s)ds≥C4Eqj(t+T3)−q∫tT3(s+T3)−1−pqαj(logsT3)qβjds≥2−1C4Eqj(t+T3)−q−pqαj∫tT31s(logsT3)qβjds≥2−1C4T−13Eqj(qβj+1)−1(t+T3)−q−pqαj(logtT3)qβj+1, | (3.29) |
where t≥T3. From (3.29) and (3.26), we deduce
L1(t)≥C4(t+T3)−p∫tT3(s+T3)−1Lp2(s)ds≥2−pC1+p4T−p3(qβj+1)−pEpqj(t+T3)−p×∫T3t(s+T3)−1−pq−p2qαj(logsT3)pqβj+pds≥2−1−pC1+p4T−p3(qβj+1)−pEpqj(t+T3)−p−pq−p2qαj×∫T3t1s(logsT3)pqβj+pds≥2−1−pC1+p4T−p3(qβj+1)−p(pqβj+p+1)−1Epqj×(t+T3)−p−pq−p2qαj(logtT3)pqβj+p+1,t≥T3. |
Therefore, the following derivation holds for j+1, namely,
αj+1=−p−pq−p2qαj, |
βj+1=p+1+pqβj, |
Ej+1=2−1−pC1+p4(qβj+1)−p(pqβj+p+1)−1Epqj. |
As a consequence, we have
αj=−p−pq−p2qαj−1=−p−pq+p3q+p3q2+p4q2αj−2=−(p+pq)j−1∑k=0(−p2q)k=p+pqp2q+1((−p2q)j−1−1), | (3.30) |
βj=p+1+pqβj−1=(p+1)(1+pq)+(pq)2βj−2=(p+1)j−1∑k=0(pq)k+(pq)jβ0=p+1pq−1((pq)j−1), | (3.31) |
Ej=2−1−pC1+p4(qβj−1+1)−p(pqβj−1+p+1)−1Epqj≥2−1−pC1+p4(pqβj−1+p+1)−1−pEpqj−1=2−1−pC1+p4β−1−pjEpqj−1≥˜E(pq)−(p+1)jEpqj−1, | (3.32) |
where ˜E=2−1−pC1+p4(p+1pq−1)−1−p. Taking j≥˜j=max{0,log˜Elog(pq)p+1−pqpq−1}, we obtain
Ej≥exp((pq)jlog(˜˜Eε)), |
where ˜˜E=C3(pq)−pq(p+1)(pq−1)2˜E1pq−1.
Combining (3.28) with (3.30)–(3.32), we have
L1(t)≥Ej(t+T3)−pαj(logtT3)βj≥exp((pq)jlog(˜˜Eε)(logtT3)p+1pq−1((pq)j−1)(t+T3)−p(p+pq)p2q+1((−p2q)j−1−1)). |
Let j=2k−1,k∈N∗. For all t≥T3, j≥˜j, we acquire
L1(t)≥exp((pq)jlog(˜˜Eε))(logtT3)p+1pq−1((pq)j−1)=exp((pq)jlog(˜˜Eε(logtT3)p+1pq−1))(logtT3)−p−1pq−1. | (3.33) |
Choose ε0=ε0(n,k1,k2,μ1,μ2,f1,f2,g1,g2,T3) small enough such that
exp((˜˜Eε0)−pq−1p+1)≥1. |
Thus, we have
˜˜Eε(logtT3)p+1pq−1>1 |
for all ε∈(0,ε0] and t>T3exp((˜˜Eε0)−pq−1p+1).
Therefore, from (3.33), we know that L1(t) and L2(t) blow up in finite time when j→∞. On the other hand, when t≤T3exp((˜˜Eε0)−pq−1p+1), we derive that the upper bound estimate of the lifespan satisfies:
T(ε)≤exp(Cε−pq−1p+1), |
where C is a positive constant.
Similar to (3.28), we make the following assumption for L2(t), that is,
L2(t)≥Ej(t+T3)−qαj(logtT3)βj,t≥T3. |
Similarly, we acquire that L1(t) and L2(t) blow up in finite time. The upper bound estimate of their lifespan satisfies
T(ε)≤exp(Cε−pq−1q+1), |
where C is a positive constant. The proof of Theorem 1.1 is finished.
This work is dedicated to investigating the weakly coupled system of semi-linear wave equations with time dependent speed of propagation, damping terms, and derivative nonlinear terms in generalized Einstein-de Sitter space-time on Rn. Taking advantage of the iteration method, we deduce formation of singularity of solutions to the coupled system (1.1). Blow-up region and upper bound lifespan estimate of solutions to the problem are established. The blow-up region of solutions in this work varies due to the influence of coefficients k1,k2. It is worth it to mention that upper bound lifespan estimate (1.12) of solutions to problem (1.1) in Theorem 1.1 coincides with the results in [39] when p=q and k1=k2. When ki=0,i=1,2 in problem (1.1) (namely, σi=μi), blow-up region Ω(n,σ1,σ2,p,q) of problem (1.1) is the same as that in [56]. In addition, the results in Theorem 1.1 complement the case that propagation speeds of generalized coupled wave equations with dampings and derivative type nonlinearities contain power functions with negative exponents.
All authors contributed equally to this work. All authors have read and approved the final version of the manuscript for publication.
The author Sen Ming would like to express his sincere thank to Professor Han Yang for his helpful suggestions and discussions. The project is supported by Natural Science Foundation of Shanxi Province of China (No. 201901D211276), Fundamental Research Program of Shanxi Province (No. 20210302123045, No. 202103021223182), National Natural Science Foundation of China (No. 11601446).
The authors declare that they have no conflict of interest.
[1] |
N. Abada, M. Benchohra, H. Hammouche, Existence and controllability results for nondensely defined impulsive semilinear functional differential inclusions, J. Differ. Equations, 246 (2009), 3834–3863. doi: 10.1016/j.jde.2009.03.004. doi: 10.1016/j.jde.2009.03.004
![]() |
[2] |
S. Abbas, M. Benchohra, Upper and lower solutions method for impulsive partial hyperbolic differential equations with fractional order, Nonlinear Anal. Hybrid Syst., 4 (2010), 406–413. doi: 10.1016/j.nahs.2009.10.004. doi: 10.1016/j.nahs.2009.10.004
![]() |
[3] |
S. M. Afonso, E. M. Bonotto, M. Federson, ˇS. Schwabik, Discontinuous local semiflows for Kurzweil equations leading to LaSalle's invariance principle for differential systems with impulses at variable times, J. Differ. Equations, 250 (2011), 2969–3001. doi: 10.1016/j.jde.2011.01.019. doi: 10.1016/j.jde.2011.01.019
![]() |
[4] |
N. U. Ahmed, K. L. Teo, S. H. Hou, Nonlinear impulsive systems on infinite dimensional spaces, Nonlinear Anal., 54 (2003), 907–925. doi: 10.1016/S0362-546X(03)00117-2. doi: 10.1016/S0362-546X(03)00117-2
![]() |
[5] | J. Banasiak, L. Arlotti, Perturbations of Positive Semigroups with Applications, London: Springer Verlag, 2006. |
[6] |
P. Y. Chen, Y. X. Li, Mixed monotone iterative technique for a class of semilinear impulsive evolution equations in Banach spaces, Nonlinear Anal., 74 (2011), 3578–3588. doi: 10.1016/j.na.2011.02.041. doi: 10.1016/j.na.2011.02.041
![]() |
[7] |
P. Y. Chen, Y. X. Li, H. Yang, Perturbation method for nonlocal impulsive evolution equations, Nonlinear Anal. Hybrid Syst., 8 (2013), 22–30. doi: 10.1016/j.nahs.2012.08.002. doi: 10.1016/j.nahs.2012.08.002
![]() |
[8] |
V. Colao, L. Muglia, H. K. Xu, An existence result for a new class of impulsive functional differential equations with delay, J. Math. Anal. Appl., 441 (2016), 668–683. doi: 10.1016/j.jmaa.2016.04.024. doi: 10.1016/j.jmaa.2016.04.024
![]() |
[9] | K. Deimling, Nonlinear Functional Analysis, New York: Springer Verlag, 1985. |
[10] |
Z. Fan, G. Li, Existence results for semilinear differential equations with nonlocal and impulsive conditions, J. Funct. Anal., 258 (2010), 1709–1727. doi: 10.1016/j.jfa.2009.10.023. doi: 10.1016/j.jfa.2009.10.023
![]() |
[11] |
M. Feˇckan, J. R. Wang, Y. Zhou, Periodic solutions for nonlinear evolution equations with non-instantaneous impulses, Nonauton. Dyn. Syst., 1 (2014), 93–101. doi: 10.2478/msds-2014-0004. doi: 10.2478/msds-2014-0004
![]() |
[12] |
M. Frigon, D. O'Regan, Existence results for first-order impulsive differential equations, J. Math. Anal. Appl., 193 (1995), 96–113. doi: 10.1006/jmaa.1995.1224. doi: 10.1006/jmaa.1995.1224
![]() |
[13] |
M. Frigon, D. O'Regan, First order impulsive initial and periodic problems with variable moments, J. Math. Anal. Appl., 233 (1999), 730–739. doi: 10.1006/jmaa.1999.6336. doi: 10.1006/jmaa.1999.6336
![]() |
[14] |
G. R. Gautam, J. Dabas, Mild solutions for class of neutral fractional functional differential equations with not instantaneous impulses, Appl. Math. Comput., 259 (2015), 480–489. doi: 10.1016/j.amc.2015.02.069. doi: 10.1016/j.amc.2015.02.069
![]() |
[15] | D. Guo, V. Lakshmikantham, Nonlinear Problems in Abstract Cones, New York: Academic Press, 1988. |
[16] |
H. D. Gou, B. L. Li, Local and global existence of mild solution to impulsive fractional semilinear integro-differential equation with noncompact semigroup, Commun. Nonlinear Sci. Numer. Simul., 42 (2017), 204–214. doi: 10.1016/j.cnsns.2016.05.021. doi: 10.1016/j.cnsns.2016.05.021
![]() |
[17] |
D. Guo, X. Liu, Extremal solutions of nonlinear impulsive integro differential equations in Banach spaces, J. Math. Anal. Appl., 177 (1993), 538–552. doi: 10.1006/jmaa.1993.1276. doi: 10.1006/jmaa.1993.1276
![]() |
[18] |
E. Hernandez, D. O'Regan, On a new class of abstract impulsive differential equations, Proc. Am. Math. Soc., 141 (2013), 1641–1649. doi: 10.1090/S0002-9939-2012-11613-2. doi: 10.1090/S0002-9939-2012-11613-2
![]() |
[19] |
H. M. Eduardo, S. M. Tanaka Aki, Global solutions for abstract impulsive differential equations, Nonlinear Anal., 72 (2010), 1280–1290. doi: 10.1016/j.na.2009.08.020. doi: 10.1016/j.na.2009.08.020
![]() |
[20] |
T. Jankowski, Monotone iterative method for first-order differential equations at resonance, Appl. Math. Comput., 233 (2014), 20–28. doi: 10.1016/j.amc.2014.01.123. doi: 10.1016/j.amc.2014.01.123
![]() |
[21] |
H. Jian, B. Liu, S. F. Xie, Monotone iterative solutions for nonlinear fractional differential systems with deviating arguments, Appl. Math. Comput., 262 (2015), 1–14. doi: 10.1016/j.amc.2015.03.127. doi: 10.1016/j.amc.2015.03.127
![]() |
[22] |
V. Lakshmikanthama, A. S. Vatsala, General uniqueness and monotone iterative technique for fractional differential equations, Appl. Math. Lett., 21 (2008), 828–834, doi: 10.1016/j.aml.2007.09.006. doi: 10.1016/j.aml.2007.09.006
![]() |
[23] |
Q. Li, Y. X. Li, Monotone iterative technique for second order delayed periodic problem in Banach spaces, Appl. Math. Comput., 270 (2015), 654–664. doi: 10.1016/j.amc.2015.08.070. doi: 10.1016/j.amc.2015.08.070
![]() |
[24] |
Y. X. Li, Maximum principles and the method of upper and lower solutions for time-periodic problems of the telegraph equations, J. Math. Anal. Appl., 327 (2007), 997–1009. doi: 10.1016/j.jmaa.2006.04.066. doi: 10.1016/j.jmaa.2006.04.066
![]() |
[25] |
Y. X. Li, A monotone iterative technique for solving the bending elastic beam equations, Appl. Math. Comput., 217 (2010), 2200–2208. doi: 10.1016/j.amc.2010.07.020. doi: 10.1016/j.amc.2010.07.020
![]() |
[26] |
Y. X. Li, Z. Liu, Monotone iterative technique for addressing impulsive integro-differential equations in Banach spaces, Nonlinear Anal., 66 (2007), 83–92. doi: 10.1016/j.na.2005.11.013. doi: 10.1016/j.na.2005.11.013
![]() |
[27] |
J. Liang, J. H. Liu, T. J. Xiao, Nonlocal impulsive problems for nonlinear differential equations in Banach spaces, Math. Comput. Modelling, 49 (2009), 798–804. doi: 10.1016/j.mcm.2008.05.046. doi: 10.1016/j.mcm.2008.05.046
![]() |
[28] | A. Pazy, Semigroup of linear operators and applications to partial differential equations, Berlin: Springer-Verlag, 1983. |
[29] |
M. Pierri, D. O'Regan, V. Rolnik, Existence of solutions for semi-linear abstract differential equations with not instantaneous impulses, Appl. Math. Comput., 219 (2013), 6743–6749. doi: 10.1016/j.amc.2012.12.084. doi: 10.1016/j.amc.2012.12.084
![]() |
[30] |
J. R. Wang, X. Z. Li, Periodic BVP for integer/fractional order nonlinear differential equations with non-instantaneous impulses, J. Appl. Math. Comput., 46 (2014), 321–334. doi: 10.1007/s12190-013-0751-4. doi: 10.1007/s12190-013-0751-4
![]() |
[31] |
J. R. Wang, Y. Zhou, M. Fe˘ckan, Nonlinear impulsive problems for fractional differential equations and Ulam stability, Comput. Math. Appl., 64 (2012), 3389–3405. doi: 10.1016/j.camwa.2012.02.021. doi: 10.1016/j.camwa.2012.02.021
![]() |
[32] |
J. R. Wang, Y. Zhou, Z. Lin, On a new class of impulsive fractional differential equations, Appl. Math. Comput., 242 (2014), 649–657. doi: 10.1016/j.amc.2014.06.002. doi: 10.1016/j.amc.2014.06.002
![]() |
[33] |
X. L. Yu, J. R. Wang, Periodic boundary value problems for nonlinear impulsive evolution equations on Banach spaces, Commun. Nonlinear Sci. Numer. Simul., 22 (2015), 980–989. doi: 10.1016/j.cnsns.2014.10.010. doi: 10.1016/j.cnsns.2014.10.010
![]() |
[34] |
V. I. Slyn'ko, C. Tunç, Stability of abstract linear switched impulsive differential equations, Automatica, 107 (2019), 433–441. doi: 10.1016/j.automatica.2019.06.001. doi: 10.1016/j.automatica.2019.06.001
![]() |
[35] |
V. I. Slyn'ko, C. Tunç, Instability of set differential equations, J. Math. Anal. Appl., 467 (2018), 935–947. doi: 10.1016/j.jmaa.2018.07.048. doi: 10.1016/j.jmaa.2018.07.048
![]() |