This paper is concerned with a chemotaxis system in a two-dimensional setting as follows:
{ut=Δu−χ∇⋅(u∇lnv)−κuv+ru−μu2+h1,vt=Δv−v+uv+h2,
with the parameters χ,κ,μ>0 and r∈R, and with the given functions h1,h2≥0. This model was originally introduced by Short et al for urban crime with the particular values χ=2,r=0 and μ=0, and the logistic source term ru−μu2 was incorporated into (⋆) by Heihoff to describe the fierce competition among criminals. Heihoff also proved that the initial-boundary value problem of (⋆) possesses a global generalized solution in the two-dimensional setting. The main purpose of this paper is to show that such a generalized solution becomes bounded and smooth at least eventually. In addition, the long-time asymptotic behavior of such a solution is discussed.
Citation: Zixuan Qiu, Bin Li. Eventual smoothness of generalized solutions to a singular chemotaxis system for urban crime in space dimension 2[J]. Electronic Research Archive, 2023, 31(6): 3218-3244. doi: 10.3934/era.2023163
[1] | Huawei Huang, Xin Jiang, Changwen Peng, Geyang Pan . A new semiring and its cryptographic applications. AIMS Mathematics, 2024, 9(8): 20677-20691. doi: 10.3934/math.20241005 |
[2] | Hailun Wang, Fei Wu, Dongge Lei . A novel numerical approach for solving fractional order differential equations using hybrid functions. AIMS Mathematics, 2021, 6(6): 5596-5611. doi: 10.3934/math.2021331 |
[3] | Shousheng Zhu . Double iterative algorithm for solving different constrained solutions of multivariate quadratic matrix equations. AIMS Mathematics, 2022, 7(2): 1845-1855. doi: 10.3934/math.2022106 |
[4] | Miloud Sadkane, Roger Sidje . Computing a canonical form of a matrix pencil. AIMS Mathematics, 2024, 9(5): 10882-10892. doi: 10.3934/math.2024531 |
[5] | Zahra Pirouzeh, Mohammad Hadi Noori Skandari, Kamele Nassiri Pirbazari, Stanford Shateyi . A pseudo-spectral approach for optimal control problems of variable-order fractional integro-differential equations. AIMS Mathematics, 2024, 9(9): 23692-23710. doi: 10.3934/math.20241151 |
[6] | Yinlan Chen, Min Zeng, Ranran Fan, Yongxin Yuan . The solutions of two classes of dual matrix equations. AIMS Mathematics, 2023, 8(10): 23016-23031. doi: 10.3934/math.20231171 |
[7] | Waleed Mohamed Abd-Elhameed, Abdullah F. Abu Sunayh, Mohammed H. Alharbi, Ahmed Gamal Atta . Spectral tau technique via Lucas polynomials for the time-fractional diffusion equation. AIMS Mathematics, 2024, 9(12): 34567-34587. doi: 10.3934/math.20241646 |
[8] | Emrah Polatlı . On some properties of a generalized min matrix. AIMS Mathematics, 2023, 8(11): 26199-26212. doi: 10.3934/math.20231336 |
[9] | Jin Li . Barycentric rational collocation method for semi-infinite domain problems. AIMS Mathematics, 2023, 8(4): 8756-8771. doi: 10.3934/math.2023439 |
[10] | Samia Bushnaq, Kamal Shah, Sana Tahir, Khursheed J. Ansari, Muhammad Sarwar, Thabet Abdeljawad . Computation of numerical solutions to variable order fractional differential equations by using non-orthogonal basis. AIMS Mathematics, 2022, 7(6): 10917-10938. doi: 10.3934/math.2022610 |
This paper is concerned with a chemotaxis system in a two-dimensional setting as follows:
{ut=Δu−χ∇⋅(u∇lnv)−κuv+ru−μu2+h1,vt=Δv−v+uv+h2,
with the parameters χ,κ,μ>0 and r∈R, and with the given functions h1,h2≥0. This model was originally introduced by Short et al for urban crime with the particular values χ=2,r=0 and μ=0, and the logistic source term ru−μu2 was incorporated into (⋆) by Heihoff to describe the fierce competition among criminals. Heihoff also proved that the initial-boundary value problem of (⋆) possesses a global generalized solution in the two-dimensional setting. The main purpose of this paper is to show that such a generalized solution becomes bounded and smooth at least eventually. In addition, the long-time asymptotic behavior of such a solution is discussed.
The importance of existence and uniqueness theorems for initial value and boundary value problems (IVP and BVP) involving the classical derivative operator is indisputable because, without them, one cannot understand modeled systems correctly and make predictions how they behave. Recently, with the popularity of fractional derivative operators such as Riemann-Liouville (R-L), Caputo (C), etc., the equations involving these operators have begun to be studied in detail (See, [1,2,3,4,5,6]). However, such a generalization leads to some difficulties and differences especially in R-L case. For instance, unlike the initial value problems involving the classical derivative, the existence of continuous solution to some IVPs in the sense of R-L derivative strictly depends on the initial values and smoothness conditions on nonlinear functions in right-hand side of equations in IVPs. For example, as shown in [7], the following initial value problem for R-L derivative of order σ∈(0,1)
{Dσω(x)=f(x,ω(x)),x>0ω(0)=b≠0 | (1.1) |
has no continuous solution when f is continuous on [0,T]×R. Therefore, this equation was investigated when subjected to different initial values (See, Theorem 3.3 in [1]). By considering the aforementioned remark, we investigate the existence and uniqueness of solutions to the following problem
{Dσω(x)=f(x,ω(x),Dσ−1ω(x)),x>0ω(0)=0,Dσ−1ω(x)|x=0=b, | (1.2) |
where σ∈(1,2), b∈R, Dσ represents the Riemann-Liouville fractional derivative of order σ, which is given by
Dσω(x)=1Γ(2−σ)d2dx2∫x0ω(t)(x−t)σ−1dt, |
and f fulfills the following condition:
(C1) Let f(x,t1,t2)∈C((0,T]×R×R) and xσ−1f(x,t1,t2)∈C([0,T]×R×R),
where C(X) represents the class of continuous functions defined on the subspace X of R3.
The equation in (1.2) was first considered by Yoruk et. al. [8], when the second initial value is also homogeneous (b=0) and the right-hand side function is continuous on [0,T]×R×R. They gave Krasnoselskii-Krein, Roger and Kooi-type uniqueness results. As seen from condition (C1), the equation in (1.2) has a singularity at x=0. Such singular equations involving R-L and Caputo derivatives were recently investigated in [7,9,10,11,12] and they proved local or global existence theorem for initial value problems involving singular equations. In the investigation of existence of solutions to these problems, converting the problem into the Volterra-type integral equation is one of the referenced tools. These integral equations have weakly and double singularities since the corresponding fractional differential equations are singular. For this reason, some new techniques or lemmas had to be developed to reveal the existence and uniqueness of solutions to integral equations (See for example [10,11]). In this paper, we also encounter a Volterra-type integral equations having a single singularity when we transform the problem (1.2) into a integral equation. For the uniqueness of solutions to the integral equation we need to generalize the definitions and methods previously given in [7,8,13] and use the tools of Lebesgue spaces such as Hölder inequality. These uniqueness theorems are the type of Nagumo-type, Krasnoselskii-Krein and Osgood which are well-known in the literature (See, [13,14,15,16]) Moreover, we give a Peano-type existence theorem for this problem as well.
We first give definition of the space of functions where we investigate the existence of solutions to the problem (1.2). The space is defined below as in [17] :
Theorem 2.1. The space of continuous functions defined on [0,T], whose R-L fractional derivative of order σ−1, 1<σ<2 are continuous on [0,T] is a Banach space when endowed with the following norm:
||ω||σ−1=||ω||∞+||Dσ−1ω||∞, |
where ||.||∞ is the supremum norm defined on the class of continuous functions. This space will be denoted by Cσ−1([0,T]).
According to this space, we define R-L integral and R-L derivative of higher order [11]. The lower terminal points of integrals in their formulas will be taken as zero.
Definition 2.1. The R-L integral of order σ>0 of a function ω(x)∈C0[0,T] is defined for all x∈[0,T] by
Iσω(x):=1Γ(σ) ∫x0ω(t)(x−t)1−σdt. | (2.1) |
Definition 2.2. For 1<σ<2 and ω(x)∈Cσ−1[0,T] with Dσ−1ω(x)∈C1(0,T]∩L1[0,T] the R-L fractional derivative Dβω is defined for all x∈(0,T] by
Dσω(x)=1Γ(2−σ)d2dx2∫x0ω(t)(x−t)σ−1dt. | (2.2) |
The local existence of solutions to problem (1.2) will be proved with the aid of Schauder fixed point theorem [18]:
Theorem 2.2. Let C be a closed, bounded, convex subset of a Banach space X:={u:I→R continuous:I⊂R closed and bounded interval}. If operator S:C→C is continuous and, if S(C) is an equicontinuous set on I, then S has at least one fixed point in C.
The one of mathematical tools used for showing the existence and uniqueness of the desired type of solution to a given initial or boundary value problem is first to convert them into an integral equation. One investigates the existence and uniqueness of the solution to the integral equation instead of the associated problem. Here, we follow this way by taking the aid of the lemma given below:
Lemma 3.1. Under condition (C1), if ω∈Cσ−1[0,T] is a solution of problem (1.2), then ω∈Cσ−1[0,T] is a solution of the following integral equation
ω(x)=bΓ(σ)xσ−1+1Γ(σ)∫x0f(t,ω(t),Dσ−1ω(t))(x−t)1−σdt | (3.1) |
and, vice versa.
Proof. We assume that ω∈Cσ−1[0,T] is a solution of problem (1.2). By condition (C1), we have f(x,ω(x),Dσ−1ω(x)) is continuous on (0,T] and xσ−1f(x,ω(x),Dσ−1ω(x)) is continuous on [0,T]. It means that f(x,ω(x),Dσ−1ω(x)) is integrable, i.e f(x,ω(x),Dσ−1ω(x))∈C(0,T]∩L1[0,T]. Then, by integrating the both sides of the equation in (1.2) and using the relation IDσ=IDDσ−1
Dσ−1ω(x)=Dσ−1ω(0)+If(x,ω(x),Dσ−1ω(x))∈C[0,T] |
is obtained. From here, by integration of the both sides of last equation and by use of IDσ−1=IDI2−σ we have
I2−σω(x)=bx+I2f(x,ω(x),Dσ−1ω(x))∈C1[0,T] | (3.2) |
where we used Dσ−1ω(0)=b and I2−σω(0)=0 since ω(0)=0.
If the operator Iσ−1 is applied to the both sides of (3.2), then by semigroup and commutative property of R-L derivative we get
Iω(x)=bxσΓ(1+σ)+I[Iσf(x,ω(x),Dσ−1ω(x))] | (3.3) |
for all x∈[0,T]. Differentiating the both sides of (3.3), we have
ω(x)=bxσ−1Γ(σ)+Iσf(x,ω(x),Dσ−1ω(x)) |
which is the equivalent to the integral Eq (3.1).
Now we suppose that ω∈Cσ−1[0,T] is a solution of integral Eq (3.1), and let us show that ω is a solution of the problem (1.2). If Dσ is applied to the both sides of (3.1), and then, if
DσIσω(x)=ω(x) for all ω∈Cσ−1[0,T] |
is used, then one can observe that ω∈Cσ−1[0,T] satisfies the equation in (1.2). Moreover, let us prove that ω∈Cσ−1[0,T] also fulfills initial value conditions. By change of variables and condition (C1) we have
ω(0)=limx→0+ω(x)=bΓ(σ)xσ−1+1Γ(σ)limx→0+∫x0f(t,ω(t),Dσ−1ω(t))(x−t)1−σdt=1Γ(σ)limx→0+∫x0tσ−1f(t,ω(t),Dσ−1ω(t))tσ−1(x−t)1−σdt=1Γ(σ)limx→0+x∫10(xτ)σ−1f(xτ,ω(xτ),Dσ−1ω(xτ))τσ−1(1−τ)1−σdτ=0, | (3.4) |
since the integral is finite. Thus ω satisfies the first initial condition in (1.2).
Now let us show that ω provides the second initial condition in (1.2). If Dσ−1 is applied to both sides of (3.1), and if the relation Dσ−1Iσh(x)=Ih(x) is used, then we can first get
Dσ−1ω(x)=b+∫x0f(t,ω(t),Dσ−1ω(t))dt. | (3.5) |
From here, by passing the limit as x→0+, we then obtain
Dσ−1ω(0)=b+limx→0+∫x0f(t,ω(t),Dσ−1ω(t))dt=b+limx→0+∫x01tσ−1tσ−1f(t,ω(t),Dσ−1ω(t))dt=b+limx→0+x2−σ∫101τσ−1(xτ)σ−1f(xτ,ω(xτ),Dσ−1ω(xτ))dτ=b, | (3.6) |
since 2−σ>0 and the integral is finite due to the continuity of tσ−1f(t,ω(t),Dσ−1ω(t)) on [0,T]. Consequently, it has been shown that any solution of (3.1) provides the problem (1.2) if condition (C1) is assumed to be satisfied.
Theorem 3.1 (Existence). Let condition (C1) be satisfied, and assume that there exist positive real numbers r1, r2 and M such that |xσ−1f(x,ω,v)|≤Mfor all(x,ω,v)∈I=[0,T]×[−r1,r1]×[b−r2,b+r2]. Then problem (1.2) admits at least one solution in Cσ−1[0,T0 ], where
T0 ={TifT<rC(b,σ,M)rC(b,σ,M)ifT≥rC(b,σ,M)≥1[rC(b,σ,M)]σ−1,ifT≥rC(b,σ,M),1≥rC(b,σ,M)and1<σ≤1.5[rC(b,σ,M)]2−σ,ifT≥rC(b,σ,M),1≥rC(b,σ,M)and1.5<σ<2, | (3.7) |
and
r=r1+r2andC(b,σ,M)=[|b|Γ(σ)+M(1+Γ(3−σ)2−σ)]. | (3.8) |
Proof. As it is known from Lemma 3.1, solutions of problem (1.2) are solutions of integral equation (3.1) as well. Moreover, the fixed points of the opeator S:Cσ−1[0,T0 ]→Cσ−1[0,T0 ] defined by
Sω(x)=bΓ(σ)xσ−1+∫x0f(t,ω(t),Dσ−1ω(t))(x−t)1−σdt | (3.9) |
interfere with solutions of the integral equation. For this reason, it is sufficient to prove that operator S admits at least one fixed point. For this, it will be verified that operator S satisfies the hypotheses of Schauder fixed-point theorem. Let us start with showing the following inclusion to be valid:
S(Br)⊂Br |
where
Br={ω∈Cσ−1[0,T0 ]:||ω||∞+||Dσ−1ω−b||∞≤r} |
is a closed compact subset of Cσ−1[0,T0 ]. Accordingly to the norm on Cσ−1[0,T0 ], upper bounds of ‖Sω(x)‖∞ and ‖Dσ−1Sω(x)−b‖∞ can be determined as follows:
|Sω(x)|≤|b|Γ(σ)xσ−1+1Γ(σ)∫x0|tσ−1f(t,ω(t),Dσ−1ω(t)|tσ−1(x−t)1−σdt≤|b|Γ(σ)xσ−1+MΓ(σ)∫10xτσ−1(1−τ)1−σdτ≤|b|Γ(σ)xσ−1+Γ(2−σ)Mx | (3.10) |
and
|Dσ−1Sω(x)−b|≤∫x0|tσ−1f(t,ω(t),Dσ−1ω(t)|tσ−1dt≤Mx2−σ∫10τ1−σdτ=Mx2−σ2−σ. | (3.11) |
From (3.10)) and (3.11),
|Sω(x)|+|Dσ−1Sω(x)−b|≤|b|Γ(σ)xσ−1+Γ(2−σ)Mx+Mx2−σ2−σ | (3.12) |
is obtained. Taking supremum over [0,T0] for a T0>0 for the right hand-side of the above equation,
|Sω(x)|+|Dσ−1Sω(x)−b|≤C(b,σ,M)Tα0 | (3.13) |
can be written, where α∈Ω={σ−1,1,2−σ}. α depends on values of b,M,σ,r. To determine T0 and α, let
C(b,σ,M)Tα0=r. |
If Tα0=rC(b,σ,M)<1, then it is observed that T0<1 for any α∈Ω. If Tα0=rC(b,σ,M)≥1, it must be T0≥1 for any α∈Ω. Thus,
supx∈[0,T0][|Sω(x)|+|Dσ−1Sω(x)−b|]≤C(b,σ,M)Tα0=r, | (3.14) |
where
T0:=[rC(b,σ,M)]1/α |
and
α={1ifrC(b,σ,M)≥1σ−1ifrC(b,σ,M)<1and1<σ≤1.52−σifrC(b,σ,M)<1and1.5≤σ<2. | (3.15) |
As a result, for all cases we obtain
||Sω||∞+||Dσ−1Sω−b||∞≤r, |
which is the desired.
Now, let us prove the equicontinuity of S(Br)⊂Cσ−1[0,T0 ]. Since the composition of uniformly continuous functions is so as well, the function xσ−1f(x,ω(x),Dσ−1ω(x)) is uniformly continuous on [0,T0 ]. Because for any ω∈Br, both ω(x) and Dσ−1ω(x) and xσ−1f(x,ω,v) are uniformly continuous on I, respectively. Therefore, for given any ϵ>0, one can find a δ=δ(ϵ)>0 so that for all x1,x2∈[0,T0 ] with |x1−x2|<δ it is
|xσ−11f(x1,ω(x1),Dσ−1ω(x))−xσ−12f(x2,ω(x2),Dσ−1ω(x2))|<Kϵ, |
where K=max(1T0 Γ(2−σ),2−σT0 2−σ). It follows that
|Sω(x1)−Sω(x2)|+|Dσ−1Sω(x1)−Dσ−1Sω(x2)|≤∫10|h(ηx1)−h(ηx2)|Γ(σ)η1−σ(1−η)σ−1xdη+∫10|h(ηx1)−h(ηx2)|η1−σx2−σdη<T0 Γ(2−σ)Kϵ+T0 2−σ2−σKϵ=ϵ, |
where h(x)=xσ−1f(x,ω(x),Dσ−1ω(x)). This implies that S(Br) is an equicontinuous set of Cσ−1[0,T0 ].
Finally, the continuity of S on Br will be proven. Assume that {ωk}∞k=1⊂Br is a sequence with ωkCσ[0,T0 ]→ω as k→∞. Then, one can easily conclude that ωk and Dσ−1ωk(t) converges uniformly to ω and Dσ−1ω(t), respectively. With these and the uniform continuity of xσ−1f(x,ω,v) on I=[0,T]×[−r1,r1]×[b−r2,b+r2], it leads to
‖Sωk−Sω‖σ−1=supx∈[0,T0 ]|1Γ(σ)∫x0[f(t,ωk(t),Dσ−1ωk(t))−f(t,ω(t),Dσ−1ω(t))](x−t)1−σdt|+supx∈[0,T0 ]|∫x0[f(t,ωk(t),Dσ−1ωk(t))−f(t,ω(t),Dσ−1ω(t))]dt|≤supηx∈[0,T0 ]∫10(ηx)σ−1|f(ηx,ωk(ηx),Dσ−1ωk(ηx))−f(ηx,ω(ηx),Dσ−1ω(ηx))|Γ(σ)ησ−1(1−η)1−σxdη+supηx∈[0,T0 ]∫10(ηx)σ−1|f(ηx,ωk(ηx),Dσ−1ωk(ηx))−f(ηx,ω(ηx),Dσ−1ω(ηx))|ησ−1x2−σdη→0ask→∞. |
In conclusion, since hypotheses of Theorem 2.2 are fulfilled, it implies that operator S admits at least one fixed point in Cσ−1[0,T0 ], which is a solution of problem (1.2) as well.
The mean value theorem for R-L derivative of order σ∈(0,1) was correctly given by [7]. Now, its counterparts for order σ∈(1,2) is given as follows:
Lemma 3.2. Let σ∈(1,2) and ω∈Cσ−1([0,T]). Then, there is a function μ:[0,T]→[0,T] with 0<μ(x)<x so that
ω(x)=Dσ−1ω(0)xσ−1Γ(σ)+Γ(2−σ)x(μ(x))σ−1Dσ−1ω(μ(x)), |
is satisfied.
The lemma can be proved by following the way used in [7] and so we omit it here. With the aid of this lemma we can obtain the Nagumo-type uniqueness:
Theorem 3.2. (Nagumo type uniqueness) Let 1<σ<2, 0<T<∞ and let condition (C1) be satisfied. Moreover, assume that there exists a positive real number L≤2−σmax(T,T2−σ)(1+Γ(3−σ)) such that the inequality
xσ−1|f(x,ω1,v1)−f(x,ω2,v2)|≤L(|ω1−ω2|+|v1−v2|) | (3.16) |
is fulfilled for all x∈[0,T] and for all ωi,vi∈R with i=1,2. Then, (1.2) has at most one solution in the space of Cσ−1([0,T0]).
Proof. We have just showed the existence of the solution to problem (1.2) in the previous theorem. For the uniqueness, we first assume that (1.2) admits two different solutions such as ω1 and ω2 in the space of Cσ−1([0,T0]). Let us define a function Φ(x) to be in the form
Φ(x):={|ω1(x)−ω2(x)|+|Dσ−1ω1(x)−Dσ−1ω2(x)|,x>0 0 ,x=0. |
Since ω1,ω2∈Cσ−1([0,T]), the continuity of Φ(x) on x∈(0,T0] can obviously be seen. For its continuity at x=0,
0≤limx→0+Φ(x)=limx→0+1Γ(σ)|∫x0f(t,ω1(t),Dσ−1ω1(t))−f(t,ω2(t),Dσ−1ω2(t))(x−t)1−σdt|+limx→0+|∫x0f(t,ω1(t),Dσ−1ω1(t))−f(t,ω2(t),Dσ−1ω2(t))dt|≤∫10limx→0+x|H(xη,ω1(xη))−H(xη,ω2(xη))|ησ−1(1−η)1−σdη+∫10limx→0+x2−σ|H(xη,ω1(xη))−H(xη,ω2(xη))|ησ−1dη=0, |
where H(x,ω(x))=xσ−1f(x,ω(x),Dσ−1ω(x)) and we made the change of variable t=xη and used condition (C1), respectively. Consequently, limx→0+Φ(x)=0=Φ(0).
The fact that Φ(x)≥0 on [0,T] allows us to choose a point x0∈(0,T] so that
0<Φ(x0)=|ω1(x0)−ω2(x0)|+|Dσ−1ω1(x0)−Dσ−1ω2(x0)|. |
By using the mean value theorem given by Lemma 3.2
|ω1(x0)−ω2(x0)|=Γ(2−σ)x0|xσ−11Dσ−1(ω1−ω2)(x1)|=Γ(2−σ)x0xσ−11|f(x1,ω1(x1),Dσ−1ω1(x1))−f(x1,ω2(x1),Dσ−1ω2(x1))| | (3.17) |
is obtained for x1∈(0,x0).
Secondly, for the estimation of |Dσ−1ω1(x0)−Dσ−1ω2(x0)|, we have from the well-known integral mean theorem for the classical calculus
|Dσ−1ω1(x0)−Dσ−1ω2(x0)|=∫x00tσ−1|f(t,ω1(t),Dσ−1ω1(t)−f(t,ω2(t),Dσ−1ω2(t)|tσ−1dt=x2−σ02−σxσ−12|f(x2,ω1(x2),Dσ−1ω1(x2))−f(x2,ω2(x2),Dσ−1ω2(x2))|, | (3.18) |
where x2∈(0,x0).
We assign x3 as one of the points x1 and x2 so that |H(x3,ω1(x3))−H(x3,ω2(x3))| :=max(|H(x1,ω1(x1))−H(x1,ω2(x1))|,|H(x2,ω1(x2))−H(x2,ω2(x2))|).
Thus, from (3.17) and (3.18), we have
0<Φ(x0)≤(Γ(2−σ)x0+x2−σ02−σ)|H(x3,ω1(x3))−H(x3,ω2(x3))|≤max(T,T2−σ)(1+Γ(3−σ)2−σ)xσ−13|f(x3,ω1(x3))−f(x3,ω2(x3))|≤xσ−13(|ω1(x3)−ω2(x3)|+|Dσ−1ω1(x3)−Dσ−1ω2(x3)|)=Φ(x3) |
since L≤2−σmax(T,T2−σ) (1+Γ(3−σ)). Repeating the same procedure for the point x3, it enables us to find some points x6∈(0,x3) so that 0<Φ(x0)≤Φ(x3)≤Φ(x6). Continuing in the same way, the sequence {x3n}∞n=1⊂[0,x0) can be constructed so that x3n→0 and
0<Φ(x0)≤Φ(x3)≤Φ(x6)≤...≤Φ(x3n)≤... | (3.19) |
However, the fact that Φ(x) is continuous at x=0 and x3n→0 leads to Φ(x3n)→Φ(0)=0, and this contradicts with (3.19). Consequently, IVP (1.2) possesses a unique solution.
Theorem 3.3. (Krasnoselskii-Krein type uniqueness) Let 1<σ<2 and T∗0=min{T0 ,1}, where T0 is defined by (3.7). Let condition (C1) be fulfilled. Furthermore, suppose that there exist a L>0 and an α∈(0,1) so that
x1−α(σ−1)|f(x,ω1,v1)−f(x,ω2,v2)|≤Γ(σ)α(σ−1)+L2(|ω1−ω2|+|v1−v2|) | (3.20) |
holds for all x∈[0,T] and for all ωi,vi∈R with i=1,2, and that there exist C>0 satisfying (1−σ)(1−α)−L(1−α)+1>0 such that
xσ−1|f(x,ω1,v1)−f(x,ω2,v2)|≤C(|ω1−ω2|α+xα(σ−1)|v1−v2|α) | (3.21) |
holds for all x∈[0,T] and for all ωi,vi∈R with i=1,2. Then, problem (1.2) has a unique solution in the space of Cσ−1([0,T∗0]).
Proof. As claimed in the previous theorem, we first assume that problem (1.2) has two different solutions such as ω1(x) and ω2(x) in Cσ−1([0,T∗0]). However, by contradiction, we will show that they are indeed equal. For this, let us first define Φ1(x)=|ω1(x)−ω2(x)| and Φ2(x)=|Dσ−1ω1(x)−Dσ−1ω2(x)| and try to find estimates for each functions by using condition (C1) and inequality (3.21). Hence, we first have
Φ1(x)≤1Γ(σ)∫x0|f(t,ω1(t),Dσ−1ω1(t))−f(t,ω2(t),Dσ−1ω2(t))|(x−t)1−σdt≤1Γ(σ)∫x0tσ−1|f(t,ω1(t),Dσ−1ω1(t))−f(t,ω2(t),Dσ−1ω2(t))|tσ−1(x−t)1−σdt≤CΓ(σ)∫x0[Φα1(x)+tα(σ−1)Φα2(x)]tσ−1(x−t)1−σdt≤CΓ(σ)(∫x0(1tσ−1(x−t)1−σ)qdt)1/q(∫x0[Φα1(t)+tα(σ−1)Φα2(t)]pdt)1/p≤CΓ(1+(1−σ)q)Γ(1+(σ−1)q))Γ(σ)x1/qΩ1/p(x) |
where we used Hölder inequality with q>1 satisfying (1−σ)q+1>0 and p=q/(q−1), and Ω(x) is defined by
Ω(x)=∫x0[Φα1(t)+tα(σ−1)Φα2(t)]pdt. |
From here, we have the following estimation
Φp1(x)≤Cxp/qΩ(x), | (3.22) |
where C is not specified here and throughout the proof. In addition to this, the upper bound for Φ2(x) can be found as follows:
Φ2(x)≤∫x0|f(t,ω1(t),Dσ−1ω1(t))−f(t,ω2(t),Dσ−1ω2(t))|dt≤∫x0tσ−1|f(t,ω1(t),Dσ−1ω1(t))−f(t,ω2(t),Dσ−1ω2(t))|tσ−1dt≤CΓ(σ)∫x0[Φα1(x)+tα(σ−1)Φα2(x)]tσ−1(x−t)1−σdt≤C(∫x0(1tσ−1)qdt)1/q(∫x0[Φα1(t)+tα(σ−1)Φα2(t)]pdt)1/p≤CΓ(1+(1−σ)q)Γ(2+(σ−1)q))x(1+q(1−σ))/qΩ1/p(x). |
From here,
Φp2(x)≤Cx(1+q(1−σ))p/qΩ(x) | (3.23) |
is then obtained. By using estimations in (3.22) and (3.23) in the derivative of Ω(x) we have
Ω′(x)=[Φα1(x)+tα(σ−1)Φα2(x)]p≤2p−1[(Φα1(x))p+tpα(σ−1)(Φα2(x))p]=2p−1[(Φp1(x))α+xpα(σ−1)(Φp2(x))α]≤2p−1[Cαxαp/qΩα(x)+Cαxpα(σ−1)x(1+q(1−σ))αp/qΩα(x)]≤Cxαp/qΩα(x). | (3.24) |
If we multiply the both sides of the above inequality with (1−α)Ω−α(x),
(1−α)Ω−α(x)Ω′(x)=ddx[Ω1−α(x)]≤Cxαp/q. |
is then obtained. Integrating the both sides of the inequality over [0,x], we get
Ω1−α(x)≤Cx(αp+q)/q, |
since Ω(0)=0. Consequently, this leads to the following estimation on Ω(x)
Ω(x)≤Cx(αp+q)/(1−α)q. | (3.25) |
By considering (3.22) and (3.23) together with (3.25), one can conclude that
Φp1(x)≤Cxp/qΩ(x)≤Cxp/qx(αp+q)/(1−α)pq=Cxp+q/(1−α)q. |
or
Φ1(x)≤Cx(p+q)/(1−α)pq=x1/(1−α), | (3.26) |
and
Φp2(x)≤Cxp(1+q(1−σ))/qΩ(x)≤Cxp(1+q(1−σ))/qx(αp+q)/(1−α)pq=Cx(1−α)(1−σ)pq+p+q(1−α)q. |
or
Φ2(x)≤Cx(1−σ)+p+q(1−α)pq=x(1−σ)+1(1−α), | (3.27) |
since p+qpq=1.
Let us now define
Ψ(x)=x−Lmax{Φ1(x),Φ2(x)}, |
where L(1−α)<1+(1−σ)(1−α). If Φ1(x)=max{Φ1(x),Φ2(x)}, then from (3.26) we get
0≤Ψ(x)≤x11−α−L, |
or in the case of Φ2(x)=max{Φ1(x),Φ2(x)}, by the inequality (3.27) we have the following
0≤Ψ(x)≤x(1−σ)+11−α−L=x(1−σ)(1−α)−L(1−α)+11−α. |
In both cases, Ψ(x) is continuous on [0,T∗0] and Ψ(0)=0. Let us now show that Ψ(x)≡0 on [0,T∗0]. For this, suppose otherwise and let Ψ(x)≢ This means and from its continuity one can say that there exists a point so that takes its maximum value at that point. Thus, let
By assuming
(3.28) |
is obtained for However, we get to a contradiction.
On the other hand, when we get
(3.29) |
which is contraction as well.
Consequently, vanishes identically on This gives us the uniqueness of solutions to the considered problem.
Theorem 3.4 (Osgood-type uniqueness). Let and be defined by (3.7), and let condition (C1) be satisfied. Furthermore, suppose that the equality
(3.30) |
is fulfilled for all and for all with where is conjugate of satisfying and
Moreover, assume that is a continuous, non-negative and non-decreasing function in so that and it satisfies
(3.31) |
for any Then, has a unique solution in the space
Proof. As made in previously given uniqueness theorems, we assume that there exist two different solutions such as and to problem (1.2) in Moreover, let and At first, we get the estimation on as follows:
(3.32) |
where we used the inequality (3.30), Hölder inequality and the assumption on respectively. From here, it follows that
(3.33) |
Similarly to above, we have
(3.34) |
This leads to
(3.35) |
Now, set
and assume that for We will show that it can not be possible under assumptions.
From the definition of one easily conclude that for each it is and there exists a so that Then, from estimations (3.33)–(3.35) and from the fact that is non-decreasing function
(3.36) |
is then obtained. It can be seen that Moreover, we have
for all From this fact, for sufficiently small we have
Furthermore, by changing variables in the above integral and by using the continuity of and , we have
for sufficiently small with and for However, this contradicts with the assumption on given in (3.31). Consequently, for i.e.
Remark 3.1. It must be pointed out that, as noted in Theorem 1.4.3 in [13], the condition that function is non-decreasing can be dropped.
Example 3.1. Let us consider the following problem
(3.37) |
with initial conditions and Let in Theorem 3.1. Then, and Hence, problem has a solution in Now, we investigate the uniqueness of the solution to the problem in Let the function in the previous theorem defined by
(3.38) |
It is obvious that is positive for . Since for it is non-decreasing. Also,
(3.39) |
for any which can be seen by considering the inequality in the neighborhood of and the divergence of the integral and by applying comparison test. Moreover, by using the concavity of the nonlinear function with respect to second and third variables, we have
(3.40) |
where for and
Hence, assumptions of Theorem 3.4 are satisfied. So, the problem has a unique solution in
In this research, we gave some sufficient conditions for the existence and uniqueness of a problem involving a nonlinear differential equations in the sense of R-L derivative when the right-hand side function has a discontinuity at zero. We presented an example associated with two theorems. Considering the literature, these results can be generalized and improved. Besides, one can obtain another uniqueness results for this problem as well.
There is no conflict of interest.
[1] |
M. Short, M. Drsogna, V. Pasour, G. Tita, P. Brantingham, A. Bertozzi, et al., A statistical model of criminal behavior, Math. Models Methods Appl. Sci., 18 (2008), 1249–1267. https://doi.org/10.1142/S0218202508003029 doi: 10.1142/S0218202508003029
![]() |
[2] |
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
![]() |
[3] |
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), 80. https://doi.org/10.1007/s00033-020-01304-w doi: 10.1007/s00033-020-01304-w
![]() |
[4] |
N. Bellomo, F. Colasuonno, D. Knopoff, J. Soler, From a systems theory of sociology to modeling the onset and evolution of criminality, Networks Heterog. Media, 10 (2015), 421–441. https://doi.org/10.3934/nhm.2015.10.421 doi: 10.3934/nhm.2015.10.421
![]() |
[5] |
H. Berestycki, J. Nadal, Self-organised critical hot spots of criminal activity, Eur. J. Appl. Math., 21 (2010), 371–399. https://doi.org/10.1017/S0956792510000185 doi: 10.1017/S0956792510000185
![]() |
[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] |
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
![]() |
[8] |
M. Short, G. Mohler, P. Brantingham, G. Tita, Gang rivalry dynamics via coupled point process networks, Discrete Contin. Dyn. Syst. Ser. B, 19 (2014), 1459–1477. https://doi.org/10.3934/dcdsb.2014.19.1459 doi: 10.3934/dcdsb.2014.19.1459
![]() |
[9] |
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
![]() |
[10] |
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), 1479–1506. https://doi.org/10.3934/dcdsb.2014.19.1479 doi: 10.3934/dcdsb.2014.19.1479
![]() |
[11] |
N. Bellomo, N. Outada, J. Soler, Y. Tao, M. Winkler, Chemotaxis and cross-diffusion models in complex environments: Models and analytic problems toward a multiscale vision, Math. Models Methods Appl. Sci., 32 (2022), 713–792. https://doi.org/10.1142/S0218202522500166 doi: 10.1142/S0218202522500166
![]() |
[12] |
M. D'Orsogna, M. Perc, Statistical physics of crime: a review, Phys. Life Rev., 12 (2015), 1–21. https://doi.org/10.1016/j.plrev.2014.11.001 doi: 10.1016/j.plrev.2014.11.001
![]() |
[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), 11. https://doi.org/10.1007/s10440-022-00484-z doi: 10.1007/s10440-022-00484-z
![]() |
[22] | 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. |
[23] |
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
![]() |
[24] |
L. Yang, X. Yang, Global existence in a two-dimensional nonlinear diffusion model for urban crime propagation, Nonliear Anal., 224 (2022), 113086. https://doi.org/10.1016/j.na.2022.113086 doi: 10.1016/j.na.2022.113086
![]() |
[25] | M. Fuest, F. Heihoff, Unboundedness phenomenon in a reduced model of urban crime, preprint, arXiv: 2109.01016. |
[26] | 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. https://doi.org/10.3934/dcdsb.2022167 |
[27] |
B. Li, Z. Wang, L. Xie, Regularization effect of the mixed-type damping in a higher-dimensional logarithmic Keller-Segel system related to crime modeling, Math. Biosci. Eng., 24 (2023), 4532–4559. https://doi.org/10.3934/mbe.2023210 doi: 10.3934/mbe.2023210
![]() |
[28] |
R. Manásevich, Q. Phan, P. Souplet, Global existence of solutions for a chemotaxis-type system arising in crime modelling, Eur. J. Appl. Math., 24 (2013), 273–296. https://doi.org/10.1017/S095679251200040X doi: 10.1017/S095679251200040X
![]() |
[29] |
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
![]() |
[30] |
D. Wang, Y. Feng, Global well-posedness and uniform boundedness of a higher dimensional crime model with a logistic source term, Math. Meth. Appl. Sci., 45 (2022), 4727–4740. https://doi.org/10.1002/mma.8066 doi: 10.1002/mma.8066
![]() |
[31] |
T. Black, Eventual smoothness of generalized solutions to a singular chemotaxis-Stokes system in 2D, J. Differ. Equations, 265 (2018), 2296–2339. https://doi.org/10.1016/j.jde.2018.04.035 doi: 10.1016/j.jde.2018.04.035
![]() |
[32] |
T. Black, C. Wu, Prescribed signal concentration on the boundary: eventual smoothness in a chemotaxis-Navier-Stokes system with logistic proliferation, Calc. Var., 61 (2022), 96. https://doi.org/10.1007/s00526-022-02201-y doi: 10.1007/s00526-022-02201-y
![]() |
[33] |
M. Ding, J. Lankeit, Generalized solutions to a chemotaxis-Navier-Stokes system with arbitrary superlinear degradation, SIAM J. Math. Anal., 54 (2022), 1022–1052. https://doi.org/10.1137/21M140907X doi: 10.1137/21M140907X
![]() |
[34] |
Y. Tao, M. Winkler, Eventual smoothness and stabilization of large-data solutions in a three-dimensional chemotaxis system with consumption of chemoattractant, J. Differ. Equations, 252 (2012), 2520–2543. https://doi.org/10.1016/j.jde.2011.07.010 doi: 10.1016/j.jde.2011.07.010
![]() |
[35] | B. Li, L. Xie, Generalized solution and eventual smoothness in a logarithmic Keller-Segel system for criminal activities, Math. Models Methods Appl. Sci., 2023. https://doi.org/10.1142/S0218202523500306 |
[36] |
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
![]() |
[37] |
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
![]() |
[38] |
T. Xiang, Finite time blow-up in the higher dimensional parabolic-elliptic-ODE minimal chemotaxis-haptotaxis system, J. Differ. Equations, 336 (2022), 44–72. https://doi.org/10.1016/j.jde.2022.07.015 doi: 10.1016/j.jde.2022.07.015
![]() |
[39] |
T. Hillen, K. Painter, M. Winkler, Convergence of a cancer invasion model to a logistic chemotaxis model, Math. Models Methods Appl. Sci., 23 (2013), 165–198. https://doi.org/10.1142/S0218202512500480 doi: 10.1142/S0218202512500480
![]() |
[40] | 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. |
[41] |
N. Bellomo, A. Bellouquid, Y. Tao, M. Winkler, Toward a mathematical theory of Keller-Segel models of pattern formation in biological tissues, Math. Models Methods Appl. Sci., 25 (2015), 1663–1763. https://doi.org/10.1142/S021820251550044X doi: 10.1142/S021820251550044X
![]() |
[42] |
T. Xiang, Boundedness and global existence in the higher-dimensional parabolic-parabolic chemotaxis system with/with-out growth source, J. Differ. Equations, 258 (2015), 4275–4323. https://doi.org/10.1016/j.jde.2015.01.032 doi: 10.1016/j.jde.2015.01.032
![]() |
[43] |
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
![]() |
[44] |
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
![]() |
[45] | O. Ladyzhenskaya, V. Solonnikov, N. Ural'tseva, Linear and Quasilinear Equations of Parabolic Type, Translated from the Russian by S. Smith. Translations of Mathematical Monographs, American Mathematical Society, Providence, R.I., 1968. |
[46] |
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
![]() |
[47] |
K. Fujie, A. Ito, M. Winkler, T. Yokota, Stabilization in a chemotaxis model for tumor invasion, Discrete Contin. Dyn. Syst., 36 (2016), 151–169. https://doi.org/10.3934/dcds.2016.36.151 doi: 10.3934/dcds.2016.36.151
![]() |
1. | Peng E, Tingting Xu, Linhua Deng, Yulin Shan, Miao Wan, Weihong Zhou, Solutions of a class of higher order variable coefficient homogeneous differential equations, 2025, 20, 1556-1801, 213, 10.3934/nhm.2025011 |