We study a logarithmic Keller-Segel system proposed by Rodríguez for crime modeling as follows:
{ut=Δu−χ∇⋅(u∇lnv)−κuv+h1,vt=Δv−v+u+h2,
in a bounded and smooth spatial domain Ω⊂Rn with n≥3, with the parameters χ>0 and κ>0, and with the nonnegative functions h1 and h2. For the case that κ=0, h1≡0 and h2≡0, recent results showed that the corresponding initial-boundary value problem admits a global generalized solution provided that χ<χ0 with some χ0>0.
In the present work, our first result shows that for the case of κ>0 such problem possesses global generalized solutions provided that χ<χ1 with some χ1>χ0, which seems to confirm that the mixed-type damping −κuv has a regularization effect on solutions. Besides the existence result for generalized solutions, a statement on the large-time behavior of such solutions is derived as well.
Citation: Bin Li, Zhi Wang, Li Xie. Regularization effect of the mixed-type damping in a higher-dimensional logarithmic Keller-Segel system related to crime modeling[J]. Mathematical Biosciences and Engineering, 2023, 20(3): 4532-4559. doi: 10.3934/mbe.2023210
[1] | A. Q. Khan, Ibraheem M. Alsulami . Discrete Leslie's model with bifurcations and control. AIMS Mathematics, 2023, 8(10): 22483-22506. doi: 10.3934/math.20231146 |
[2] | Xiongxiong Du, Xiaoling Han, Ceyu Lei . Dynamics of a nonlinear discrete predator-prey system with fear effect. AIMS Mathematics, 2023, 8(10): 23953-23973. doi: 10.3934/math.20231221 |
[3] | Fatao Wang, Ruizhi Yang, Yining Xie, Jing Zhao . Hopf bifurcation in a delayed reaction diffusion predator-prey model with weak Allee effect on prey and fear effect on predator. AIMS Mathematics, 2023, 8(8): 17719-17743. doi: 10.3934/math.2023905 |
[4] | Abdul Qadeer Khan, Ayesha Yaqoob, Ateq Alsaadi . Neimark-Sacker bifurcation, chaos, and local stability of a discrete Hepatitis C virus model. AIMS Mathematics, 2024, 9(11): 31985-32013. doi: 10.3934/math.20241537 |
[5] | Xiaoming Su, Jiahui Wang, Adiya Bao . Stability analysis and chaos control in a discrete predator-prey system with Allee effect, fear effect, and refuge. AIMS Mathematics, 2024, 9(5): 13462-13491. doi: 10.3934/math.2024656 |
[6] | Ahmad Suleman, Rizwan Ahmed, Fehaid Salem Alshammari, Nehad Ali Shah . Dynamic complexity of a slow-fast predator-prey model with herd behavior. AIMS Mathematics, 2023, 8(10): 24446-24472. doi: 10.3934/math.20231247 |
[7] | Weili Kong, Yuanfu Shao . Bifurcations of a Leslie-Gower predator-prey model with fear, strong Allee effect and hunting cooperation. AIMS Mathematics, 2024, 9(11): 31607-31635. doi: 10.3934/math.20241520 |
[8] | Ali Al Khabyah, Rizwan Ahmed, Muhammad Saeed Akram, Shehraz Akhtar . Stability, bifurcation, and chaos control in a discrete predator-prey model with strong Allee effect. AIMS Mathematics, 2023, 8(4): 8060-8081. doi: 10.3934/math.2023408 |
[9] | Abdul Qadeer Khan, Ayesha Yaqoob, Ateq Alsaadi . Discrete Hepatitis C virus model with local dynamics, chaos and bifurcations. AIMS Mathematics, 2024, 9(10): 28643-28670. doi: 10.3934/math.20241390 |
[10] | Fethi Souna, Salih Djilali, Sultan Alyobi, Anwar Zeb, Nadia Gul, Suliman Alsaeed, Kottakkaran Sooppy Nisar . Spatiotemporal dynamics of a diffusive predator-prey system incorporating social behavior. AIMS Mathematics, 2023, 8(7): 15723-15748. doi: 10.3934/math.2023803 |
We study a logarithmic Keller-Segel system proposed by Rodríguez for crime modeling as follows:
{ut=Δu−χ∇⋅(u∇lnv)−κuv+h1,vt=Δv−v+u+h2,
in a bounded and smooth spatial domain Ω⊂Rn with n≥3, with the parameters χ>0 and κ>0, and with the nonnegative functions h1 and h2. For the case that κ=0, h1≡0 and h2≡0, recent results showed that the corresponding initial-boundary value problem admits a global generalized solution provided that χ<χ0 with some χ0>0.
In the present work, our first result shows that for the case of κ>0 such problem possesses global generalized solutions provided that χ<χ1 with some χ1>χ0, which seems to confirm that the mixed-type damping −κuv has a regularization effect on solutions. Besides the existence result for generalized solutions, a statement on the large-time behavior of such solutions is derived as well.
The topic of boundary value problems is an interesting area of research in view of its applications in applied and technical sciences. In the recent years, the class of nonlocal fractional order boundary value problems involving different fractional derivatives (such as Riemann-Liouville, Caputo, etc.) received an overwhelming interest from many researchers. For the details of a variety of nonlocal single-valued and multivalued boundary value problems involving different types of fractional order derivative operators, we refer the reader to the text [1], articles [2,3,4,5,6,7] and the references cited therein. There has been shown a great enthusiasm in developing the existence theory for Hilfer, ψ-Hilfer and (k,ψ) Hilfer type fractional differential equations equipped with different types of boundary conditions, for instance, see [8,9,10,11,12,13,14,15,16].
Nonlocal boundary conditions are found to be more plausible and practical in contrast to the classical boundary conditions in view of their applicability to describe the changes happening within the given domain. Closed boundary conditions are found to be of great help in describing the situation when there is no fluid flow along the boundary or through it. The free slip condition is also a type of the closed boundary conditions which describes the situation when there is a flow along the boundary, but there is no flow perpendicular to it. Such conditions are also useful in the study of sandpile model [17,18], honeycomb lattice [19], deblurring problems [20], closed-aperture wavefield decomposition in solid media [21], vibration analysis of magneto-electro-elastic cylindrical composite panel [22], etc.
Now we review some works on the boundary value problems with closed boundary conditions. In [23], the authors studied the single-valued and multivalued fractional boundary value problems with open and closed boundary conditions. A three-dimensional Neumann boundary value problem with a generalized boundary condition in a domain with a smooth closed boundary was discussed in [24]. For some interesting results on impulsive fractional differential equations with closed boundary conditions, see the articles [25,26].
The objective of the present work is to investigate a new class of mixed nonlinear boundary value problems involving a right Caputo fractional derivative, mixed Riemann-Liouville fractional integral operators, and multipoint variant of closed boundary conditions. In precise terms, we consider the following fractional order nonlocal and nonlinear problem:
CDαT−y(t)+λIρT−Iσ0+h(t,y(t))=f(t,y(t)),t∈J:=[0,T], | (1.1) |
y(T)=m∑ı=1(piy(ξi)+Tqiy′(ξi)),Ty′(T)=m∑ı=1(riy(ξi)+Tviy′(ξi)), | (1.2) |
where CDαT− denote the right Caputo fractional derivative of order α∈(1,2], IρT− and Iσ0+ represent the right and left Riemann-Liouville fractional integral operators of orders ρ,σ>0 respectively, f,h:[0,T]×R→R are given continuous functions and λ,pi,qi,ri,vi∈R,i∈{1,2,3,...,m}, and ξi∈(0,T). Notice that the integro-differential Eq (1.1) contains the usual and mixed Riemann-Liouville integrals type nonlinearities. The boundary conditions (1.2) can be interpreted as the values of the unknown function and its derivative at the right end-point T of the interval [0,T] are proportional to a linear combination of these values at arbitrary nonlocal positions ξi∈(0,T). Physically, the nonlocal multipoint closed boundary conditions provide a flexible mechanism to close the boundary at arbitrary positions in the given domain instead of the left end-point of the domain.
Here we emphasize that much of the literature on fractional differential equations contains the left-sided fractional derivatives and there are a few works dealing with the right-sided fractional derivatives. For instance, the authors in [27,28] studied the problems involving the right-handed Riemann–Liouville fractional derivative operators, while a problem containing the right-handed Caputo fractional derivative was considered in [29]. The problem studied in the present paper is novel in the sense that it solves an integro-differential equation with a right Caputo fractional derivative and mixed nonlinearities complemented with a new concept of nonlocal multipoint closed boundary conditions. The results accomplished for the problems (1.1) and (1.2) will enrich the literature on boundary value problems involving the right-sided fractional derivative operators. The present work is also significant as it produces several new results as special cases as indicated in the last section.
The rest of the paper is arranged as follows. In Section 2, we present an auxiliary lemma which is used to transform the given nonlinear problem into a fixed-point problem. Section 3 contains the main results and illustrative examples. Some interesting observations are presented in the last Section 4.
Let us begin this section with some definitions [30].
Definition 2.1. The left and right Riemann-Liouville fractional integrals of order β>0 for g∈L1[a,b], existing almost everywhere on [a,b], are respectively defined by
Iβa+g(t)=∫ta(t−s)β−1Γ(β)g(s)dsandIβb−g(t)=∫bt(s−t)β−1Γ(β)g(s)ds. |
Definition 2.2. For g∈ACn[a,b], the right Caputo fractional derivative of order β∈(n−1,n],n∈N, existing almost everywhere on [a,b], is defined by
CDβb−g(t)=(−1)n∫bt(s−t)n−β−1Γ(n−β)g(n)(s)ds. |
In the following lemma, we solve a linear variant of the fractional integro-differential equation (1.1) supplemented with multipoint closed boundary conditions (1.2).
Lemma 2.1. Let H,F∈C[0,T] and Δ≠0. Then the linear problem
{CDαT−y(t)+λIρT−Iσ0+H(t)=F(t),t∈J:=[0,T],y(T)=m∑ı=1(piy(ξi)+Tqiy′(ξi)),Ty′(T)=m∑ı=1(riy(ξi)+Tviy′(ξi)),0<ξi<T, | (2.1) |
is equivalent to the integral equation
y(t)=∫Tt(s−t)α−1Γ(α)[F(s)−λIρT−Iσ0+H(s)]ds+b1(t){m∑ı=1pi∫Tξi(s−ξi)α−1Γ(α)[F(s)−λIρT−Iσ0+H(s)]ds−Tm∑ı=1qi∫Tξi(s−ξ)α−2Γ(α−1)[F(s)−λIρT−Iσ0+H(s)]ds}+b2(t){m∑ı=1ri∫Tξi(s−ξi)α−1Γ(α)[F(s)−λIρT−Iσ0+H(s)]ds−Tm∑i=1vi∫Tξi(s−ξi)α−2Γ(α−1)[F(s)−λIρT−Iσ0+H(s)]ds}, | (2.2) |
where
b1(t)=1Δ(tS6−S7−TS9+T),b2(t)=1Δ[(1−S1)t+S2+TS4−T],Δ=(S1−1)(S7+TS9−T)−S6(S2+TS4−T),S1=m∑ı=1pi,S2=m∑ı=1piξi,S3=m∑ı=1piAi,S4=m∑ı=1qi,S5=m∑ı=1qiBi,S6=m∑ı=1ri,S7=m∑ı=1riξi,S8=m∑ı=1riAi,S9=m∑ı=1vi,S10=m∑ı=1viBi,Ai=IαT−[F(ξi)−λIρT−Iσ0+H(ξi)],Bi=−Iα−1T−[F(ξi)−λIρT−Iσ0+H(ξi)]. | (2.3) |
Proof. Applying the right fractional integral operator IαT− to the integro-differential equation in (2.1), we get
y(t)=IαT−F(t)−λIα+ρT−Iσ0+H(t)−c0−c1t, | (2.4) |
where c0 and c1 are unknown arbitrary constants. Using (2.4) in the nonlocal closed boundary conditions of (2.1), we obtain
{(S1−1)c0+(S2+TS4−T)c1=S3+TS5,S6c0+(S7+TS9−T)c1=S8+TS10, | (2.5) |
where Si,i=1,…,10, are given in (2.3).
Solving the system (2.5) for c0 and c1, we find that
c0=1Δ[(S7+TS9−T)(S3+TS5)−(S2+TS4−T)(S8+TS10)],c1=1Δ[−S6(S3+TS5)+(S1−1)(S8+TS10)], |
where Δ is given in (2.3). Substituting the above values of c0 and c1 in (2.4) together with the notation (2.3), we obtain the solution (2.2). The converse of this lemma can be obtained by direct computation. This completes the proof.
This section is devoted to our main results concerning the existence and uniqueness of solutions for the problems (1.1) and (1.2).
In order to convert the problems (1.1) and (1.2) into a fixed point problem, we define an operator V:X→X by using Lemma 2.1 as follows:
Vy(t)=∫Tt(s−t)α−1Γ(α)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds+b1(t){m∑ı=1pi∫Tξi(s−ξi)α−1Γ(α)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds−Tm∑ı=1qi∫Tξi(s−ξ)α−2Γ(α−1)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds}+b2(t){m∑ı=1ri∫Tξi(s−ξi)α−1Γ(α)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds−Tm∑i=1vi∫Tξi(s−ξi)α−2Γ(α−1)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds},t∈J, | (3.1) |
where X=C([0,T],R) denotes the Banach space of all continuous functions from [0,T]→R equipped with the norm ‖y‖=sup{|y(t)|:t∈[0,T]}. Notice that the fixed point problem Vy(t)=y(t) is equivalent to the boundary value problems (1.1) and (1.2) and the fixed points of the operator V are its solutions.
In the forthcoming analysis, we use the following estimates:
∫Tt(s−t)α+ρ−1Γ(α+ρ)Iσ0+ds=∫Tt(s−t)α+ρ−1Γ(α+ρ)∫s0(s−u)σ−1Γ(σ)duds≤Tσ(T−t)α+ρΓ(σ+1)Γ(α+ρ+1),∫Tξi(s−ξi)α+ρ−1Γ(α+ρ)Iσ0+ds=∫Tξi(s−ξi)α+ρ−1Γ(α+ρ)∫s0(s−u)σ−1Γ(σ)duds≤Tσ(T−ξi)α+ρΓ(σ+1)Γ(α+ρ+1), |
where we have used uσ≤Tσ,ρ,σ>0.
In the sequel, we set
Ω1=1Γ(α+1){Tα+¯b1[m∑ı=1|pi|(T−ξi)α+αTm∑ı=1|qi|(T−ξi)α−1]+¯b2[m∑ı=1|ri|(T−ξi)α+αTm∑ı=1|vi|(T−ξi)α−1]},Ω2=|λ|TσΓ(σ+1)Γ(α+ρ+1){Tα+ρ+¯b1[m∑ı=1|pi|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|qi|(T−ξi)α+ρ−1]+¯b2[m∑ı=1|ri|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|vi|(T−ξi)α+ρ−1]}, | (3.2) |
where
¯b1=maxt∈[0,T]|b1(t)|,¯b2=maxt∈[0,T]|b2(t)|. |
In the following, Krasnosel'ski∨i's fixed point theorem [31] is applied to prove our first existence result for the problems (1.1) and (1.2).
Theorem 3.1. Assume that:
(H1) There exists L>0 such that |f(t,x)−f(t,y)|≤L|x−y|,∀t∈[0,T],x,y∈R;
(H2) There exists K>0 such that |h(t,x)−h(t,y)|≤K|x−y|,∀t∈[0,T],x,y∈R;
(H3) |f(t,y)|≤δ(t) and |h(t,y)|≤θ(t), where δ,θ∈C([0,T],R+).
Then, the problems (1.1) and (1.2) has at least one solution on [0,T] if Lγ1+Kγ2<1, where
γ1=TαΓ(α+1),γ2=|λ|Tα+ρ+σΓ(σ+1)Γ(α+ρ+1). | (3.3) |
Proof. Introduce the ball Bη={y∈X:‖y‖≤η}, with
η≥‖δ‖Ω1+‖θ‖Ω2. | (3.4) |
Now we verify the hypotheses of Krasnosel'ski∨i's fixed point theorem in three steps by splitting the operator V:X→X defined by (3.1) on Bη as V=V1+V2, where
V1y(t)=∫Tt(s−t)α−1Γ(α)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds,t∈J,V2y(t)=b1(t){m∑ı=1pi∫Tξi(s−ξi)α−1Γ(α)[f(s,y(s))ds−λIρT−Iσ0+h(s,y(s))]ds−Tm∑ı=1qi∫Tξi(s−ξ)α−2Γ(α−1)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds}+b2(t){m∑ı=1ri∫Tξi(s−ξi)α−1Γ(α)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds−Tm∑i=1vi∫Tξi(s−ξi)α−2Γ(α−1)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds},t∈J. |
(i) For y,x∈Bη, we have
‖V1y+V2x‖≤supt∈[0,T]{∫Tt(s−t)α−1Γ(α)[|f(s,y(s))|+|λ|IρT−Iσ0+|h(s,y(s))|]ds+|b1(t)|{m∑ı=1|pi|∫Tξi(s−ξi)α−1Γ(α)[|f(s,x(s))|+|λ|IρT−Iσ0+|h(s,x(s))|]ds+Tm∑ı=1|qi|∫Tξi(s−ξi)α−2Γ(α−1)[|f(s,x(s))|+|λ|IρT−Iσ0+|h(s,x(s))|]ds}+|b2(t)|{m∑ı=1|ri|∫Tξi(s−ξi)α−1Γ(α)[|f(s,x(s))|+|λ|IρT−Iσ0+|h(s,x(s))|]ds+Tm∑ı=1|vi|∫Tξi(s−ξi)α−2Γ(α−1)[|f(s,x(s))|+|λ|IρT−Iσ0+|h(s,x(s))|]ds}}≤‖δ‖supt∈[0,T]{∫Tt(s−t)α−1Γ(α)ds+|b1(t)|[m∑ı=1|pi|∫Tξi(s−ξi)α−1Γ(α)ds+Tm∑ı=1|qi|∫Tξi(s−ξi)α−2Γ(α−1)ds]+|b2(t)|[m∑ı=1|ri|∫Tξi(s−ξi)α−1Γ(α)ds+Tm∑ı=1|vi|∫Tξi(s−ξi)α−2Γ(α−1)ds]}+‖θ‖|λ|supt∈[0,T]{∫Tt(s−t)α+ρ−1Γ(α+ρ)Iσ0+ds+|b1(t)|[m∑ı=1∫Tξi(s−ξi)α+ρ−1Γ(α+ρ)Iσ0+ds+Tm∑ı=1|qi|∫Tξi(s−ξi)α+ρ−2Γ(α+ρ−1)Iσ0+ds]+|b2(t)|[m∑ı=1|ri|∫Tξi(s−ξi)α+ρ−1Γ(α+ρ)Iσ0+ds+Tm∑ı=1|vi|∫Tξi(s−ξi)α+ρ−2Γ(α+ρ−1)Iσ0+ds]}≤‖δ‖supt∈[0,T]{(T−t)αΓ(α+1)+|b1(t)|[m∑ı=1|pi|(T−ξi)αΓ(α+1)+Tm∑ı=1|qi|(T−ξi)α−1Γ(α)]+|b2(t)|[m∑ı=1|ri|(T−ξi)αΓ(α+1)+Tm∑ı=1|vi|(T−ξi)α−1Γ(α)]}+‖θ‖|λ|TσΓ(σ+1)supt∈[0,T]{(T−t)α+ρΓ(α+ρ+1)ds+|b1(t)|[m∑ı=1(T−ξi)α+ρΓ(α+ρ+1)+Tm∑ı=1|qi|(T−ξi)α+ρ−1Γ(α+ρ)]+|b2(t)|[m∑ı=1|ri|(T−ξi)α+ρΓ(α+ρ+1)+Tm∑ı=1|vi|(T−ξi)α+ρ−1Γ(α+ρ)]}≤‖δ‖Γ(α+1){Tα+¯b1[m∑ı=1|pi|(T−ξi)α+αTm∑ı=1|qi|(T−ξi)α−1]+¯b2[m∑ı=1|ri|(T−ξi)α+αTm∑ı=1|vi|(T−ξi)α−1]}+‖θ‖|λ|TσΓ(σ+1)Γ(α+ρ+1){Tα+ρ+¯b1[m∑ı=1|pi|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|qi|(T−ξi)α+ρ−1]+¯b2[m∑ı=1|ri|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|vi|(T−ξi)α+ρ−1]}≤‖β‖Ω1+‖θ‖Ω2<η, |
where we used (3.4). Thus V1y+V2x∈Bη.
(ii) Using (H1) and (H2), it is easy to show that
‖V1y−V1x‖≤supt∈[0,T]{∫Tt(s−t)α−1Γ(α)|f(s,y(s))−f(s,x(s))|ds+|λ|∫Tt(s−t)α+ρ−1Γ(α+ρ)Iσ0+|h(s,y(s))−h(s,x(s))|ds}≤(Lγ1+Kγ2)‖y−x‖, |
which, in view of the condition Lγ1+Kγ2<1, implies that the operator V1 is a contraction.
(iii) Continuity of the functions f,h implies that the operator V2 is continuous. In addition, V2 is uniformly bounded on Bη as
‖V2y‖≤supt∈[0,T]{|b1(t)|[m∑ı=1|pi|∫Tξi(s−ξi)α−1Γ(α)|f(s,y(s))|ds+|λ|m∑ı=1|pi|∫Tξi(s−ξi)α+ρ−1Γ(α+ρ)Iσ0+|h(s,y(s))|ds+Tm∑ı=1|qi|∫Tξi(s−ξi)α−2Γ(α−1)|f(s,y(s))|ds+|λ|Tm∑ı=1|qi|∫Tξi(s−ξi)α+ρ−2Γ(α+ρ−1)Iσ0+|h(s,y(s))|ds]+|b2(t)|[m∑ı=1|ri|∫Tξi(s−ξi)α−1Γ(α)|f(s,y(s))|+|λ|m∑ı=1|ri|∫Tξi(s−ξi)α+ρ1(α+ρ)Iσ0+|h(s,y(s))|ds+Tm∑ı=1|vi|∫Tξi(s−ξi)α−2Γ(α−1)|f(s,y(s))|ds+|λ|Tm∑ı=1|vi|∫Tξi(s−ξi)α+ρ−2)Γ(α+ρ−1)|h(s,y(s))|ds]}≤‖δ‖supt∈[0,T]{|b1(t)|[m∑ı=1|pi|∫Tξi(s−ξi)α−1Γ(α)ds+Tm∑ı=1|qi|∫Tξi(s−ξi)α−2Γ(α−1)ds]+|b2(t)|[m∑ı=1|ri|∫Tξi(s−ξi)α−1Γ(α)ds+Tm∑ı=1|vi|∫Tξi(s−ξi)α−2Γ(α−1)ds]}+|λ|‖θ‖supt∈[0,T]{|b1(t)|[m∑ı=1|pi|∫Tξi(s−ξi)α+ρ−1Γ(α+ρ)Iσ0+ds+Tm∑ı=1|qi|∫Tξi(s−ξi)α+ρ−2Γ(α+ρ−1)Iσ0+ds]+|b2(t)|[m∑ı=1|ri|∫Tξi(s−ξi)α+ρ−1Γ(α+ρ)Iσ0+ds+Tm∑ı=1|vi|∫Tξi(s−ξi)α+ρ−2)Γ(α+ρ−1)Iσ0+ds]}≤‖δ‖supt∈[0,T]{|b1(t)|[m∑ı=1|pi|(T−ξi)αΓ(α+1)+Tm∑ı=1|qi|(T−ξi)α−1Γ(α)]+|b2(t)|[m∑ı=1|ri|(T−ξi)αΓ(α+1)+Tm∑ı=1|vi|(T−ξi)α−1Γ(α)]+|λ|‖θ‖TσΓ(σ+1)supt∈[0,T]{|b1(t)|[m∑ı=1|pi|(T−ξi)α+ρΓ(α+ρ+1)+Tm∑ı=1|qi|(T−ξi)α+ρ−1Γ(α+ρ)]+|b2(t)|[m∑ı=1|ri|(T−ξi)α+ρΓ(α+ρ+1)+Tm∑ı=1|vi|(T−ξi)α+ρ−1)Γ(α+ρ)]}≤‖δ‖(Ω1−γ1)+‖θ‖(Ω2−γ2), |
where Ωi, and γi, i=1,2, are defined in (3.2) and (3.3), respectively. To show the compactness of V2, we fix sup(t,y)∈[0,T]×Bη|f(t,y)|=¯f, sup(t,y)∈[0,T]×Bη|h(t,y)|=¯h. Then, for 0<t1<t2<T, we have
|(V2y)(t2)−(V2y)(t1)|≤|b1(t2)−b1(t1)|{m∑ı=1|pi|∫Tξi(s−ξi)α−1Γ(α)[|f(s,y(s))|+|λ|IρT−Iσ0+|h(s,y(s))|]ds+Tm∑ı=1|qi|∫Tξi(s−ξi)α−2Γ(α−1)[|f(s,y(s))|+|λ|IρT−Iσ0+|h(s,y(s))|]ds}+|b2(t2)−b2(t1)|{m∑ı=1|ri|∫Tξ(s−ξ)α−1Γ(α)[|f(s,y(s))|+λ|IρT−Iσ0+|h(s,y(s))|ds]+Tm∑ı=1|vi|∫Tξi(s−ξi)α−2Γ(α−1)[|f(s,y(s))|+|λ|IρT−Iσ0+|h(s,y(s))|]ds}≤|S6||t2−t1||Δ|{¯fΓ(α+1)[m∑ı=1|pi|(T−ξi)α+αTm∑i=1|qi|(T−ξi)α−1]+¯h|λ|TσΓ(σ+1)Γ(α+ρ+1)[m∑ı=1|pi|(T−ξi)α+ρ+(α+ρ)Tm∑i=1|qi|(T−ξi)α+ρ−1]}+|S1−1||t2−t1||Δ|{¯fΓ(α+1)[m∑ı=1|ri|(T−ξi)α+αTm∑i=1|vi|(T−ξi)α−1]+¯h|λ|TσΓ(σ+1)Γ(α+ρ+1)[m∑ı=1|ri|(T−ξi)α+ρ+(α+ρ)Tm∑i=1|vi|(T−ξi)α+ρ−1]}, |
which tends to zero, independent of y, as t2→t1. This shows that V2 is equicontinuous. It is clear from the foregoing arguments that the operator V2 is relatively compact on Bη. Hence, by the Arzelá-Ascoli theorem, V2 is compact on Bη.
In view of the foregoing arguments (i)–(iii), the hypotheses of the Krasnosel'ski∨i's fixed point theorem [31] are satisfied. Hence, the operator V1+V2=V has a fixed point, which implies that the problems (1.1) and (1.2) has at least one solution on [0,T]. The proof is finished.
Remark 3.1. Interchanging the roles of the operators V1 and V2 in the previous result, the condition Lγ1+Kγ2<1 changes to the following one:
L(Ω1−γ1)+K(Ω2−γ2)<1, |
where Ω1,Ω2 and γ1,γ2 are defined in (3.2) and (3.3) respectively.
The following existence result relies on Leray-Schauder nonlinear alternative [32].
Theorem 3.2. Suppose that the following conditions hold:
(H4) There exist continuous nondecreasing functions ϕ1,ϕ2:[0,∞)→(0,∞) such that ∀(t,y)∈[0,1]×R, |f(t,y)|≤ω1(t)ϕ1(‖y‖) and |h(t,y)|≤ω2(t)ϕ2(‖y‖), where ω1,ω2∈C([0,T],R+);
(H5)There exists a constant M>0 such that
M‖ω1‖ϕ1(M)Ω1+‖ω2‖ϕ2(M)Ω2>1. |
Then, the problems (1.1) and (1.2) has at least one solution on [0,T].
Proof. We firstly show that the operator V:X→X defined by (3.1) is completely continuous.
(i) V maps bounded sets into bounded sets in X.
Let y∈Br={y∈X:‖y‖≤r}, where r is a fixed number. Then, using the strategy employed in the proof of Theorem 3.1, we obtain
‖Vy‖≤‖ω1‖ϕ1(r)Γ(α+1){Tα+¯b1[m∑ı=1|pi|(T−ξi)α+αTm∑ı=1|qi|(T−ξi)α−1]+¯b2[m∑ı=1|ri|(T−ξi)α+αTm∑ı=1|vi|(T−ξi)α−1]}+|λ|Tσ‖ω2‖ϕ2(r)Γ(σ+1)Γ(α+ρ+1){Tα+ρ+¯b1[m∑ı=1|pi|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|qi|(T−ξi)α+ρ−1]+¯b2[m∑ı=1|ri|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|vi|(T−ξi)α+ρ−1]}=‖ω1‖ϕ1(r)Ω1+‖ω2‖ϕ2(r)Ω2<∞. |
(ii) V maps bounded sets into equicontinuous sets.
Let 0<t1<t2<T and y∈Br. Then, we obtain
|Vy(t2)−Vy(t1)|≤|∫Tt2(s−t2)α−1−(s−t1)α−1Γ(α)f(s,y(s))ds+∫t2t1(s−t1)α−1Γ(α)f(s,y(s))ds−λ∫Tt2(s−t2)α+ρ−1−(s−t1)α+ρ−1Γ(α+ρ)Iσ0+h(s,y(s))ds−λ∫t2t1(s−t1)α+ρ−1Γ(α+ρ)Iσ0+h(s,y(s))ds|+|b1(t2)−b1(t1)|{|m∑ı=1pi∫Tξi(s−ξ)α−1Γ(α)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds|+|Tm∑ı=1qi∫Tξi(s−ξi)α−2Γ(α−1)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds|}+|b2(t2)−b2(t1)|{|m∑ı=1ri∫Tξi(s−ξ)α−1Γ(α)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds|+|Tm∑ı=1vi∫Tξi(s−ξi)α−2Γ(α−1)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds|}≤|∫Tt2(s−t2)α−1−(s−t1)α−1Γ(α)f(s,y(s))ds+∫t2t1(s−t1)α−1Γ(α)f(s,y(s))ds|+|λ∫Tt2(s−t2)α+ρ−1−(s−t1)α+ρ−1Γ(α+ρ)Iσ0+h(s,y(s))ds+λ∫t2t1(s−t1)α+ρ−1Γ(α+ρ)Iσ0+h(s,y(s))ds|+|S6||t2−t1|Δ|{|m∑ı=1pi∫Tξi(s−ξi)α−1Γ(α)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds|+|Tm∑ı=1qi∫Tξi(s−ξi)α−2Γ(α−1)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds|}+|S1−1||t2−t1|Δ{|m∑ı=1ri∫Tξi(s−ξi)α−1Γ(α)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds|+|Tm∑ı=1vi∫Tξi(s−ξi)α−2Γ(α−1)[f(s,y(s))−λIρT−Iσ0+h(s,y(s))]ds|}≤ω1(t)Φ1(r)Γ(α+1){|(T−t2)α−(T−t1)α|+2|(t2−t1)α|+|t2−t1||Δ|[|S6|(m∑ı=1|pi|(T−ξi)α+αTm∑ı=1|qi|(T−ξi)α−1)+|S1−1|(m∑ı=1|ri|(T−ξi)α+αTm∑ı=1|vi|(T−ξi)α−1)]}+|λ|Tσω2(t)ϕ2(r)Γ(σ+1)Γ(α+ρ+1){|(T−t2)α+ρ−(T−t1)α+ρ+2|t2−t1|α+ρ+|t2−t1||Δ|[|S6|(m∑ı=1|pi|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|qi|(T−ξi)α+ρ−1)+|S1−1|(m∑ı=1|ri|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|vi|(T−ξi)α+ρ−1)]}. |
Notice that the right-hand side of the above inequality tends to 0 as t2→t1, independent of y∈Br. Thus, it follows by the Arzelá–Ascoli theorem that the operator V:X→X is completely continuous.
The conclusion of the Leray-Schauder nonlinear alternative [32] will be applicable once it is shown that there exists an open set U⊂C([0,T],R) with y≠νVy for ν∈(0,1) and y∈∂U. Let y∈C([0,T],R) be such that y=νVy for ν∈(0,1). As argued in proving that the operator V is bounded, one can obtain that
|y(t)|=|νVy(t)|≤|ω1(t)|ϕ(‖y‖)Ω1+|ω2(t)|ψ(‖y‖)Ω2, |
which can be written as
‖y‖‖ω1‖ϕ(‖y‖)Ω1+‖ω2‖ψ(‖y‖)Ω2≤1. |
On the other hand, we can find a positive number M such that ‖y‖≠M by assumption (H5). Let us set
W={y∈X:‖y‖<M}. |
Clearly, ∂W contains a solution only when ‖y‖=M. In other words, we cannot find a solution y∈∂W satisfying y=νVy for some ν∈(0,1). In consequence, the operator V has a fixed point y∈¯W, which is a solution of the problems (1.1) and (1.2). The proof is finished.
Here we apply Banach contraction mapping principle to establish the uniqueness of solutions for the problems (1.1) and (1.2).
Theorem 3.3. If the conditions (H1) and (H2) hold, then the problems (1.1) and (1.2) has a unique solution on [0,T] if
LΩ1+KΩ2<1, | (3.5) |
where Ω1 and Ω2 are defined in (3.2).
Proof. In the first step, we show that VBκ⊂Bκ, where Bκ={y∈X:‖y‖≤κ} with
κ≥f0Ω1+h0Ω21−(LΩ1+KΩ2),f0=supt∈[0,T]|f(t,0)|,h0=supt∈[0,T]|h(t,0)|. |
For y∈Bκ and using the condition (H1), we have
|f(t,y)|=|f(t,y)−f(t,0)+f(t,0)|≤|f(t,y)−f(t,0)|+|f(t,0)|≤L‖y‖+f0≤Lr+f0. | (3.6) |
Similarly, using (H2), we get
|h(t,y)|≤Kr+h0. | (3.7) |
In view of (3.6) and (3.7), we obtain
‖Vy‖≤supt∈[0,T]|Vy(t)|≤supt∈[0,T]{∫Tt(s−t)α−1Γ(α)[|f(s,y(s))|+|λ|IρT−Iσ0+|h(s,y(s))|]ds+|b1(t)|{m∑ı=1pi∫Tξi(s−ξi)α−1Γ(α)[|f(s,y(s))|+|λ|IρT−Iσ0+|h(s,y(s))|]ds+Tm∑ı=1|qi|∫Tξi(s−ξi)α−2Γ(α−1)[|f(s,y(s))|+|λ|IρT−Iσ0+|h(s,y(s))|]ds}+|b2(t)|{m∑ı=1ri∫Tξi(s−ξi)α−1Γ(α)[|f(s,y(s))|+|λ|IρT−Iσ0+|h(s,y(s))|]ds+Tm∑i=1|vi|∫Tξi(s−ξi)α−2Γ(α−1)[|f(s,y(s))|+|λ|IρT−Iσ0+|h(s,y(s))|]ds}}≤(Lr+f0)supt∈[0,T]{∫Tt(s−t)α−1Γ(α)ds+|b1(t)|[m∑i=1|pi|∫Tξi(s−ξi)α−1Γ(α)ds+Tm∑ı=1|qi|∫Tξi(s−ξi)α−2Γ(α−1)ds]+|b2(t)|[m∑i=1|ri|∫Tξi(s−ξi)α−1Γ(α)ds+Tm∑ı=1|vi|∫Tξi(s−ξi)α−2Γ(α−1)ds]}+|λ|(Kr+h0)supt∈[0,T]{∫Tt(s−t)α+ρ−1Γ(α+ρ)Iσ0+ds+|b1(t)|[m∑i=1|pi|∫Tξi(s−ξi)α+ρ−1Γ(α+ρ)Iσ0+ds+Tm∑ı=1|qi|∫Tξi(s−ξi)α+ρ−2Γ(α+ρ−1)Iσ0+ds]+|b2(t)|[m∑i=1|ri|∫Tξi(s−ξi)α+ρ−1Γ(α+ρ)Iσ0+ds+Tm∑ı=1|vi|∫Tξi(s−ξi)α+ρ−2Γ(α+ρ−1)Iσ0+ds]}≤(Lr+f0)Γ(α+1){Tα+¯b1[m∑ı=1|pi|(T−ξi)α+αTm∑ı=1|qi|(T−ξi)α−1]+¯b2[m∑ı=1|ri|(T−ξi)α+αTm∑ı=1|vi|(T−ξi)α−1]}+Tσ|λ|(Kr+h0)Γ(σ)Γ(α+ρ+1){Tα+ρ+¯b1[m∑ı=1|pi|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|qi|(T−ξi)α+ρ−1]+¯b2[m∑ı=1|ri|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|vi|(T−ξi)α+ρ−1]}=(Lr+f0)Ω1+(Kr+h0)Ω2<κ, |
which implies that Vy∈Bκ, for any y∈Bκ. Therefore, VBκ⊂Bκ.
Next, we prove that V is a contraction. For that, let x,y∈X and t∈[0,T]. Then, by the conditions (H1) and (H2), we obtain
‖Vy−Vx‖=supt∈[0,T]|(Vy)(t)−(Vx)(t)|≤supt∈[0,T]{∫Tt(s−t)α−1Γ(α)|f(s,y(s))−f(s,x(s))|ds+|λ|∫Tt(s−t)α+ρ−1Γ(α+ρ)Iσ0+|h(s,y(s))−h(s,x(s))|ds+|b1(t)|[m∑ı=1|pi|(∫Tξi(s−ξi)α−1Γ(α)|f(s,y(s))−f(s,x(s))|ds+|λ|∫Tξi(s−ξi)α+ρ−1Γ(α+ρ)Iσ0+|h(s,y(s))−h(s,x(s))|ds)+Tm∑ı=1|qi|(∫Tξi(s−ξi)α−2Γ(α−1)|f(s,y(s))−f(s,x(s))|ds+|λ|∫Tξi(s−ξi)α+ρ−2Γ(α+ρ−1)Iσ0+|h(s,y(s))−h(s,x(s))|ds)]+|b2(t)|[m∑ı=1|ri|(∫Tξi(s−ξi)α−1Γ(α)|f(s,y(s))−f(s,x(s))|ds+|λ|∫Tξi(s−ξi)α+ρ−1Γ(α+ρ)Iσ0+|h(s,y(s))−h(s,x(s))|ds)+Tm∑ı=1|vi|(∫Tξi(s−ξi)α−2Γ(α−1)|f(s,y(s))−f(s,x(s))|ds+|λ|∫Tξi(s−ξi)α+ρ−2Γ(α+ρ−1)Iσ0+|h(s,y(s))−h(s,x(s))|ds)]}≤LΓ(α+1){Tα+¯b1[m∑ı=1|pi|(T−ξi)α+αTm∑ı=1|qi|(T−ξi)α−1]+¯b2[m∑ı=1|ri|(T−ξi)α+αTm∑ı=1|vi|(T−ξi)α−1]}+Tσ|λ|KΓ(σ+1)Γ(α+ρ+1){Tα+ρ+¯b1[m∑ı=1|pi|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|qi|(T−ξi)α+ρ−1]+¯b2[m∑ı=1|ri|(T−ξi)α+ρ+(α+ρ)Tm∑ı=1|vi|(T−ξi)α+ρ−1]}‖y−x‖=(LΩ1+KΩ2)‖y−x‖, |
which shows that \mathcal{V} is a contraction in view of the condition (3.5). Therefore, we deduce by Banach contraction mapping principle that there exists a unique fixed point for the operator \mathcal{V} , which corresponds to a unique solution for the problems (1.1) and (1.2) on [0, T] . The proof is completed.
In this subsection, we construct examples for illustrating the abstract results derived in the last two subsections. Let us consider the following problem:
\begin{equation} \left\{\begin{array} {ll} &D^{9/8}_{1-}y(t)+3I^{7/3}_{1-}I^{3/4}_{0+}h(t,y(t)) = f(t,y(t)),\ t\in J: = [0,1],\\ & y(T) = \sum\limits_{ı = 1}^3 p_iy(\xi_i)+\sum\limits_{ı = 1}^3 q_iy'(\xi_i), \, y'(T) = \sum\limits_{ı = 1}^3 r_i y(\xi_i)+\sum\limits_{ı = 1}^3 v_iy'(\xi_i),\; ,\; 0 < \xi_i < 1. \end{array} \right. \end{equation} | (3.8) |
Here \alpha = 9/8, \rho = 7/3, \sigma = 3/4, \lambda = 3, \xi_1 = 3/7, \xi_2 = 2/3, \xi_3 = 4/5, p_1 = 1/2, p_2 = 1/3, p_3 = 1/4, q_1 = -2, q_2 = -3, q_3 = -4, r_1 = 1, r_2 = -1, r_3 = 3, v_1 = -2/7, v_2 = -3/7, v_3 = -4/7. Using the given data, it is found that
\overline{b}_1 = \max\limits_{t\in [0,1]}|b_1(t)| = |b_1(t)|_{t = 1} \approx 0.1112461491, \; \; \; \overline{b}_2 = \max\limits_{t\in [0,1]}|b_2(t)| = |b_2(t)|_{t = 1}\approx 0.3364235041. |
In consequence, we get \Omega_1\approx 2.517580993, \Omega_2\approx 0.3543113654 ( \Omega_1 , \Omega_2 are defined in (3.2)).
(i) For illustrating Theorem 3.1, we consider the functions
\begin{equation} f(t,y) = \frac{m_1}{2t+25}\Big(\frac{y^2}{1+y^2}+\cos 3t+1\Big), \; \; h(t,y) = \frac{m_2}{3\sqrt{t^2+64}}\bigg(2\tan^{-1}y+\sin t+e^{-t/2}\bigg), \end{equation} | (3.9) |
where m_1 and m_2 are finite positive real numbers. Observe that
|f(t,y)|\leq \delta(t) = \frac{m_1(2+\cos 3t)}{2t+25}, \; \; \; \; \; |h(t,y)|\leq \theta(t) = \frac{m_2(\pi+\sin t+e^{-t/2})}{3\sqrt{t^2+64}}, |
and f(t, y) and h(t, y) respectively satisfy the conditions (H_1) and (H_2) with \mathcal{L} = 2m_1/25 and \mathcal{K} = m_2/24. Moreover, \gamma_1 \approx0.9438765902 and \gamma_2 \approx 0.2972831604. By the condition \mathcal{L}\gamma_1+\mathcal{K}\gamma_2 < 1 , we get
\begin{equation} 0.0755101272 m_1+0.0123867984 m_2 < 1 \end{equation} | (3.10) |
For the values of m_1 and m_2 satisfying the inequality (3.10), the hypothesis of Theorem 3.1 is satisfied. Hence, it follows by the conclusion of Theorem 3.1 that the problem (3.8) with f(t, y) and h(t, y) given in (3.9) has at least one solution on [0, 1]. If the values m_1 and m_2 do not satisfy the inequality (3.10), then Theorem 3.1 does not guarantee the existence of at least one solution to the problem (3.8) with f(t, y) and h(t, y) given in (3.9) for such values of m_1 and m_2 .
(ii) In order to illustrate Theorem 3.2, we take the following functions (instead of (3.9)) in the problem (3.8):
\begin{equation} f(t,y) = \frac{e^{-3t}}{t^2+3}\big[\sin y+1/5\big], \; \; \; \; \; h(t,y) = \frac{2}{7\sqrt{t^3+1}}\Big(\frac{|y|}{1+|y|}|y|+\pi/4\Big). \end{equation} | (3.11) |
Observe that the assumption (H_4) is satisfied as |f(t, y)| \le \omega_1(t) \phi_1(\|y\|) and |h(t, y)| \le \omega_2(t) \phi_2(\|y\|), where \omega_1(t) = e^{-3t}/(t^2+3) , \phi_1(\|y\|) = (\|y\|+1/5), \omega_2(t) = 2/(7\sqrt{t^3+1}) , \phi_2(\|y\|) = (\|y\|+\pi/4). It is easy to see that \| \omega_1\| = 1/3 and \omega_2\| = 2/7. By the condition ( H_5 ), we find that M > 4.151876169 . Thus, all the conditions of Theorem 3.2 are satisfied and hence the problem (3.8) with f(t, y) and h(t, y) given by (3.11) has at least one solution on [0, 1].
(iii) The conditions (H_1) and (H_2) are respectively satisfied by f(t, y) and h(t, y) defined in (3.9) with \mathcal{L} = 2m_1/25 and \mathcal{K} = m_2/24. By the condition (3.5), we have
\begin{equation} 0.20140647944 m_1+0.0147629736 m_2 < 1. \end{equation} | (3.12) |
Clearly, all the assumptions of Theorem 3.3 hold true with the values of m_1 and m_2 satisfying the inequality (3.12). In consequence, the problem (3.8) with f(t, y) and h(t, y) given in (3.11) has a unique solution on [0, 1]. In case, we take m_1 = m_2 = m in (3.9), then the condition (3.12) implies the existence of a unique solution for the problem at hand for m < 4.62600051. One can notice that Theorem 3.1 does not guarantee the existence of a unique solution to the problem (3.8) with f(t, y) and h(t, y) given in (3.9) for the values of m_1 and m_2, which do not satisfy the inequality (3.12).
In this study, we discussed the existence and uniqueness of solutions under different assumptions for a boundary value problem involving a right Caputo fractional derivative with usual and mixed Riemann-Liouville integrals type nonlinearities, equipped with nonlocal multipoint version of the closed boundary conditions. Our results are not only new in the given configuration, but also yield some new results as special cases. Here are some examples.
● If \lambda = 0 in (1.1), then our results correspond to the fractional differential equation {^CD^{\alpha}_{T-}}y(t) = f(t, y(t)) with the boundary conditions (1.2).
● In case, we take q_i = 0, r_i = 0, \forall i = 1, \dots, m in the results of this paper, we obtain the ones for the Eq (1.1) supplemented with boundary conditions: y(T) = \sum\limits_{ı = 1}^m p_iy(\xi_i), \; \; y'(T) = \sum\limits_{ı = 1}^m v_iy'(\xi_i).
● We get the results for the Eq (1.1) complemented with boundary conditions: y(T) = T\sum\limits_{ı = 1}^m q_iy'(\xi_i), \; \; Ty'(T) = \sum\limits_{ı = 1}^m r_i y(\xi_i) by taking p_i = 0, v_i = 0, \, \forall i = 1, \dots, m in the obtained results.
The Deanship of Scientific Research (DSR) at King Abdulaziz University (KAU), Jeddah, Saudi Arabia, has funded this project under grant No. (KEP-PhD: 35-130-1443).
The authors declare no conflict of interest.
[1] |
N. Rodríguez, On the global well-posedness theory for a class of PDE models for criminal activity, Phys. D Nonlinear Phenom., 260 (2013), 191–200. https://doi.org/10.1016/j.physd.2012.08.003 doi: 10.1016/j.physd.2012.08.003
![]() |
[2] |
M. Short, M. D'Orsogna, V. Pasour, G. Tita, P. Brantingham, A. Bertozzi, et al., A statistical model of criminal behavior, Math. Mod. Meth. Appl. Sci., 18 (2008), 1249–1267. https://doi.org/10.1142/S0218202508003029 doi: 10.1142/S0218202508003029
![]() |
[3] |
M. Short, A. Bertozzi, P. Brantingham, G. Tita, Dissipation and displacement of hotspots in reaction-diffusion model of crime, Proc. Natl. Acad. Sci. USA, 107 (2010), 3961–3965. https://doi.org/0.1073/pnas.0910921107 doi: 10.1073/pnas.0910921107
![]() |
[4] |
H. Berestycki, J. Wei, M. Winter, Existence of symmetric and asymmetric spikes for a crime hotspot model, SIAM J. Math. Anal., 46 (2014), 691–719. https://doi.org/10.1137/130922744 doi: 10.1137/130922744
![]() |
[5] |
R. Cantrell, C. Cosner, R. Manásevich, Global bifurcation of solutions for crime modeling equations, SIAM J. Math. Anal., 44 (2012), 1340–1358. https://doi.org/10.1137/110843356 doi: 10.1137/110843356
![]() |
[6] |
Y. Gu, Q. Wang, G. Yi, Stationary patterns and their selection mechanism of urban crime models with heterogeneous near-repeat victimization effect, Eur. J. Appl. Math., 28 (2017), 141–178. https://doi.org/10.1017/S0956792516000206 doi: 10.1017/S0956792516000206
![]() |
[7] |
T. Kolokolnikov, M. Ward, J. Wei, The stability of steady-state hot-spot patterns for a reaction-diffusion model of urban crime, Discrete Contin. Dyn. Syst. Ser. B, 19 (2014), 1373–1410. https://doi.org/10.3934/dcdsb.2014.19.1373 doi: 10.3934/dcdsb.2014.19.1373
![]() |
[8] |
D. Lloyd, H. O'Farrell, On localised hotspots of an urban crime model, Phys. D Nonlinear Phenom., 253 (2013), 23–39. https://doi.org/10.1016/j.physd.2013.02.005 doi: 10.1016/j.physd.2013.02.005
![]() |
[9] |
D. Lloyd, N. Santitissadeekorn, M. Short, Exploring data assimilation and forecasting issues for an urban crime model, Eur. J. Appl. Math., 27 (2016), 451–478. https://doi.org/10.1017/S0956792515000625 doi: 10.1017/S0956792515000625
![]() |
[10] |
L. Mei, J. Wei, The existence and stability of spike solutions for a chemotax is system modeling crime pattern formation, Math. Models Methods Appl. Sci., 30 (2020), 1727–1764. https://doi.org/10.1142/S0218202520500359 doi: 10.1142/S0218202520500359
![]() |
[11] |
M. Short, A. Bertozzi, P. Brantingham, Nonlinear patterns in urban crime: hotspots, bifurcations, and suppression, SIAM J. Appl. Dyn. Syst., 9 (2010), 462–483. https://doi.org/10.1137/090759069 doi: 10.1137/090759069
![]() |
[12] |
W. Tse, M. Ward, Hotspot formation and dynamics for a continuum model of urban crime, Eur. J. Appl. Math., 27 (2016), 583–624. https://doi.org/10.1017/S0956792515000376 doi: 10.1017/S0956792515000376
![]() |
[13] |
N. Rodríguez, A. Bertozzi, Local existence and uniqueness of solutions to a PDE model for criminal behavior, Math. Models Methods Appl. Sci., 20 (2010), 1425–1457. https://doi.org/10.1142/S0218202510004696 doi: 10.1142/S0218202510004696
![]() |
[14] |
N. Rodríguez, M. Winkler, On the global existence and qualitative behavior of one-dimensional solutions to a model for urban crime, Eur. J. Appl. Math., 33 (2022), 919–959. https://doi.org/10.1017/S0956792521000279 doi: 10.1017/S0956792521000279
![]() |
[15] |
Q. Wang, D. Wang, Y. Feng, Global well-posedness and uniform boundedness of urban crime models: One-dimensional case, J. Differ. Equations, 269 (2020), 6216–6235. https://doi.org/10.1016/j.jde.2020.04.035 doi: 10.1016/j.jde.2020.04.035
![]() |
[16] |
M. Freitag, Global solutions to a higher-dimensional system related to crime modeling, Math. Meth. Appl. Sci., 41 (2018), 6326–6335. https://doi.org/10.1002/mma.5141 doi: 10.1002/mma.5141
![]() |
[17] |
J. Shen, B. Li, Mathematical analysis of a continuous version of statistical models for criminal behavior, Math. Meth. Appl. Sci., 43 (2020), 409–426. https://doi.org/10.1002/mma.5898 doi: 10.1002/mma.5898
![]() |
[18] |
J. Ahn, K. Kang, J. Lee, Global well-posedness of logarithmic Keller-Segel type systems, J. Differ. Equations, 287 (2021), 185–211. https://doi.org/10.1016/j.jde.2021.03.053 doi: 10.1016/j.jde.2021.03.053
![]() |
[19] |
Y. Tao, M. Winkler, Global smooth solutions in a two-dimensional cross-diffusion system modeling propagation of urban crime, Commun. Math. Sci., 19 (2021), 829–849. https://doi.org/10.4310/CMS.2021.v19.n3.a12 doi: 10.4310/CMS.2021.v19.n3.a12
![]() |
[20] |
M. Winkler, Global solvability and stabilization in a two-dimensional cross-diffusion system modeling urban crime propagation, Ann. Inst. Henri Poincaré, Anal. Non Linéaire, 36 (2019), 1747–1790. https://doi.org/10.1016/j.anihpc.2019.02.004 doi: 10.1016/j.anihpc.2019.02.004
![]() |
[21] |
Y. Jiang, L. Yang, Global solvability and stabilization in a three-dimensional cross-diffusion system modeling urban crime propagation, Acta Appl. Math., 178 (2022). https://doi.org/10.1007/s10440-022-00484-z doi: 10.1007/s10440-022-00484-z
![]() |
[22] |
N. Rodríguez, M. Winkler, Relaxation by nonlinear diffusion enhancement in a two-dimensional cross-diffusion model for urban crime propagation, Math. Models Methods Appl. Sci., 30 (2020), 2105–2137. https://doi.org/10.1142/S0218202520500396 doi: 10.1142/S0218202520500396
![]() |
[23] |
F. Heihoff, Generalized solutions for a system of partial differential equations arising from urban crime modeling with a logistic source term, Z. Für Angew. Math. Phys., 71 (2020). https://doi.org/10.1007/s00033-020-01304-w doi: 10.1007/s00033-020-01304-w
![]() |
[24] | B. Li, L. Xie, Generalized solution to a 2D parabolic-parabolic chemotaxis system for urban crime: Global existence and large time behavior, submitted for publication, 2022. |
[25] |
P. Jones, P. Brantingham, L. Chayes, Statistical models of criminal behavior: The effects of law enforcement actions, Math. Models Methods Appl. Sci., 20 (2010), 1397–1423. https://doi.org/10.1142/S0218202510004647 doi: 10.1142/S0218202510004647
![]() |
[26] |
A. Pitcher, Adding police to a mathematical model of burglary, Eur. J. Appl. Math., 21 (2010), 401–419. https://doi.org/10.1017/S0956792510000112 doi: 10.1017/S0956792510000112
![]() |
[27] |
J. Zipkin, M. Short, A. Bertozzi, Cops on the dots in a mathematical model of urban crime and police response, Discrete Contin. Dyn. Syst. Ser. B, 19 (2014), no. 5, 1479–1506. https://doi.org/10.3934/dcdsb.2014.19.1479 doi: 10.3934/dcdsb.2014.19.1479
![]() |
[28] |
W. Tse, M. Ward, Asynchronous instabilities of crime hotspots for a 1-D reaction-diffusion model of urban crime with focused police patrol, SIAM J. Appl. Dyn. Syst., 17 (2018), 2018–2075. https://doi.org/10.1137/17M1162585 doi: 10.1137/17M1162585
![]() |
[29] |
A. Buttenschoen, T. Kolokolnikov, M. Ward, J. Wei, Cops-on-the-dots: the linear stability of crime hotspots for a 1-D reaction-diffusion model of urban crime, Eur. J. Appl. Math., 31 (2020), 871–917. https://doi.org/10.1017/S0956792519000305 doi: 10.1017/S0956792519000305
![]() |
[30] |
B. Li, L. Xie, Global large-data generalized solutions to a two-dimensional chemotaxis system stemming from crime modelling, Discrete Contin. Dyn. Syst. Ser. B, 2022 (2022). https://doi.org/10.3934/dcdsb.2022167 doi: 10.3934/dcdsb.2022167
![]() |
[31] |
N. Rodríguez, Q. Wang, L. Zhang, Understanding the effects of on- and off-hotspot policing: Evidence of hotspot, oscillating, and chaotic activities, SIAM J. Appl. Dyn. Syst., 20 (2021), 1882–1916. https://doi.org/10.1137/20M1359572 doi: 10.1137/20M1359572
![]() |
[32] |
E. Keller, L. Segel, Initiation of slime mold aggregation viewed as an instability, J. Theor. Biol., 26 (1970), 399–415. https://doi.org/10.1016/0022-5193(70)90092-5 doi: 10.1016/0022-5193(70)90092-5
![]() |
[33] |
K. Fujie, Boundedness in a fully parabolic chemotaxis system with strongly singular sensitivity, J. Math. Anal. Appl., 424 (2015), 675–684. https://doi.org/10.1016/j.jmaa.2014.11.045 doi: 10.1016/j.jmaa.2014.11.045
![]() |
[34] |
M. Winkler, Global solutions in a fully parabolic chemotaxis system with singular sensitivity, Math. Methods Appl. Sci., 34 (2011), 176–190. https://doi.org/10.1002/mma.1346 doi: 10.1002/mma.1346
![]() |
[35] |
J. Lankeit, A new approach toward boundedness in a two-dimensional parabolic chemotaxis system with singular sensitivity, Math. Meth. Appl. Sci., 39 (2016), 394–404. https://doi.org/10.1002/mma.3489 doi: 10.1002/mma.3489
![]() |
[36] |
C. Stinner, M. Winkler, Global weak solutions in a chemotaxis system with large singular sensitivity, Nonliear Anal. Real Word Appl., 12 (2011), 3727–3740. https://doi.org/10.1016/j.nonrwa.2011.07.006 doi: 10.1016/j.nonrwa.2011.07.006
![]() |
[37] |
J. Lankeit, M. Winkler, A generalized solution concept for the Keller-Segel system with logarithmic sensitivity: global solvability for large nonradial data, NoDEA-Nonlinear Differ. Equ. Appl., 24 (2017). https://doi.org/10.1007/s00030-017-0472-8 doi: 10.1007/s00030-017-0472-8
![]() |
[38] |
A. Zhigun, Generalised supersolutions with mass control for the Keller-Segel system with logarithmic sensitivity, J. Math. Anal. Appl., 467 (2018), 1270–1286. https://doi.org/10.1016/j.jmaa.2018.08.001 doi: 10.1016/j.jmaa.2018.08.001
![]() |
[39] |
M. Winkler, Unlimited growth in logarithmic Keller-Segel systems, J. Differ. Equations, 309 (2022), 74–97. https://doi.org/10.1016/j.jde.2021.11.026 doi: 10.1016/j.jde.2021.11.026
![]() |
[40] |
M. Winkler, T. Yokota, Stabilization in the logarithmic Keller-Segel system, Nonlinear Anal., 170 (2018), 123–141. https://doi.org/10.1016/j.na.2018.01.002 doi: 10.1016/j.na.2018.01.002
![]() |
[41] |
J. Ahn, Global well-posedness and asymptotic stabilization for chemotaxis system with signal-dependent sensitivity, J. Differ. Equations, 266 (2019), 6866–6904. https://doi.org/10.1016/j.jde.2018.11.015 doi: 10.1016/j.jde.2018.11.015
![]() |
[42] |
Q. Hou, Z. Wang, K. Zhao, Boundary layer problem on a hyperbolic system arising from chemotaxis, J. Differ. Equations, 261 (2016), 5035–5070. https://doi.org/10.1016/j.jde.2016.07.018 doi: 10.1016/j.jde.2016.07.018
![]() |
[43] |
H. Jin, J. Li, Z. Wang, Asymptotic stability of traveling waves of a chemotaxis model with singular sensitivity, J. Differ. Equations, 255 (2013), 193–219. doilinkhttps://doi.org/10.1016/j.jde.2013.04.002 doi: 10.1016/j.jde.2013.04.002
![]() |
[44] |
H. Li, K. Zhao, Initial-boundary value problems for a system of hyperbolic balance laws arising from chemotaxis, J. Differ. Equations, 258 (2015), 302–308. https://doi.org/10.1016/j.jde.2014.09.014 doi: 10.1016/j.jde.2014.09.014
![]() |
[45] |
J. Li, T. Li, Z. Wang, Stability of traveling waves of the Keller-Segel system with logarithmic sensitivity, Math. Models Methods Appl. Sci., 24 (2014), 2819–2849. https://doi.org/10.1142/S0218202514500389 doi: 10.1142/S0218202514500389
![]() |
[46] |
Y. Tao, L. Wang, Z. Wang, Large-time behavior of a parabolic-parabolic chemotaxis model with logarithmic sensitivity in one dimension, Discrete Contin. Dyn. Syst. Ser. B, 18 (2013), 821–845. https://doi.org/10.3934/dcdsb.2013.18.821 doi: 10.3934/dcdsb.2013.18.821
![]() |
[47] |
Z. Wang, Z. Xiang, P. Yu, Asymptotic dynamics in a singular chemotaxis system modeling onset of tumor angiogenesis, J. Differ. Equations, 260 (2016), 2225–2258. https://doi.org/10.1016/j.jde.2015.09.063 doi: 10.1016/j.jde.2015.09.063
![]() |
[48] |
M. Winkler, Renormalized radial large-data solutions to the higher-dimensional Keller-Segel system with singular sensitivity and signal absorption, J. Differ. Equations, 264 (2018), 2310–2350. https://doi.org/10.1016/j.jde.2017.10.029 doi: 10.1016/j.jde.2017.10.029
![]() |
[49] |
M. Winkler, The two-dimensional Keller-Segel system with singular sensitivity and signal absorption: global large-data solutions and their relaxation properties, Math. Models Methods Appl. Sci., 26 (2016), 987–1024. https://doi.org/10.1142/S0218202516500238 doi: 10.1142/S0218202516500238
![]() |
[50] | B. Li, L. Xie, Generalized solution and its eventual smoothness to a logarithmic Keller-Segel system for criminal activities, submitted for publication, 2022. |
[51] |
M. Winkler, Large-data global generalized solutions in a chemotaxis system with tensor-valued sensitivities, SIAM J. Math. Anal., 47 (2015), 3092–3115. https://doi.org/dx.doi.org/10.1137/140979708 doi: 10.1137/140979708
![]() |
[52] |
M. Aida, K. Osaka, T. Tsujikawa, M. Mimura, Chemotaxis and growth system with sigular sensitivity function, Nonliear Anal. Real Word Appl., 6 (2005), 323–336. https://doi.org/10.1016/j.nonrwa.2004.08.011 doi: 10.1016/j.nonrwa.2004.08.011
![]() |
[53] |
X. Cao, Global bounded solutions of the higher-dimensional Keller-Segel system under smallness conditions in optimal spaces, Discrete Contin. Dyn. Syst., 35 (2015), 1891–1904. https://doi.org/10.3934/dcds.2015.35.1891 doi: 10.3934/dcds.2015.35.1891
![]() |
[54] |
M. Winkler, Aggregation vs. global diffusive behavior in the higher-dimensional Keller-Segel model, J. Differ. Equations, 248 (2010), 2889–2905. https://doi.org/10.1016/j.jde.2010.02.008 doi: 10.1016/j.jde.2010.02.008
![]() |
[55] |
J. Simon, Compact sets in the space L^p(0, T;B), Ann. Mat. Pura Appl., 146 (1996), 65–96. https://doi.org/10.1007/BF01762360 doi: 10.1007/BF01762360
![]() |
[56] | O. Ladyzhenskaya, N. Ural'tseva, Linear and Quasilinear Elliptic Equations, Translated from the Russian by Scripta Technica, Inc. Translation editor: Leon Ehrenpreis, Academic Press, New York-London, 1968. |
1. | Bashir Ahmad, Muhammed Aldhuain, Ahmed Alsaedi, Existence Results for a Right-Caputo Type Fractional Differential Equation with Mixed Nonlinearities and Nonlocal Multipoint Sub-strips Type Closed Boundary Conditions, 2024, 45, 1995-0802, 6457, 10.1134/S1995080224606969 |