Home Classical solutions to Cauchy problems for parabolic–elliptic systems of Keller-Segel type
Article Open Access

Classical solutions to Cauchy problems for parabolic–elliptic systems of Keller-Segel type

  • Michael Winkler EMAIL logo
Published/Copyright: May 11, 2023

Abstract

The Cauchy problem in R n , n 2 , for

u t = Δ u ( u S v ) , 0 = Δ v + u , ( )

is considered for general matrices S R n × n . A theory of local-in-time classical existence and extensibility is developed in a framework that differs from those considered in large parts of the literature by involving bounded classical solutions. Specifically, it is shown that for all non-negative initial data belonging to BUC ( R n ) L p ( R n ) with some p [ 1 , n ) , there exist T max ( 0 , ] and a uniquely determined u C 0 ( [ 0 , T max ) ; BUC ( R n ) ) C 0 ( [ 0 , T max ) ; L p ( R n ) ) C ( R n × ( 0 , T max ) ) such that with v Γ u , and with Γ denoting the Newtonian kernel on R n , the pair ( u , v ) forms a classical solution of ( ) in R n × ( 0 , T max ) , which has the property that

if T max < , then both limsup t T max u ( , t ) L ( R n ) = and limsup t T max v ( , t ) L ( R n ) = .

An exemplary application of this provides a result on global classical solvability in cases when S + 1 is sufficiently small, where 1 = diag ( 1 , , 1 ) .

MSC 2010: 35K55; 35B65; 35Q92; 92C17

1 Introduction

The Cauchy problem for the Keller-Segel system,

(1.1) u t = Δ u ( u S v ) , x R n , t > 0 , 0 = Δ v + u , x R n , t > 0 , u ( x , 0 ) = u 0 ( x ) , x R n ,

has been at the core of interest in numerous studies concerned with fundamental principles of self-organization and spontaneous aggregation both in biology and in astrophysics: in the case when the cross-diffusion mechanism in (1.1) represents complete attraction in the sense that the tensor S R n × n coincides with the unit matrix 1 = diag ( 1 , , 1 ) R n × n , this problem has been used since the 1940s to describe self-interaction of particles, mediated, e.g., through gravitational or electrical potentials [14]. As very similar types of interplay have been found responsible for striking experimental observations on chemotactically induced spontaneous aggregation in microbial populations, from the 1970s on system (1.1) has received considerable additional attention [2,58], with the choice S = 1 again being of predominant interest initially, but with more general classes of S appearing especially in the literature of the past few years, quite in line with experimental and refined modeling-oriented literature reflecting rotational parts in cross-diffusive movement [9].

According to these roles of (1.1) in application contexts, a predominant focus of corresponding analytical studies has ever been on phenomena related to blow-up of solutions, as widely being understood as a mathematical expression of spontaneous aggregation [6]. In fact, one of the striking technical advantages of the markedly slim parabolic–elliptic setting of the Cauchy problem (1.1) consists in the circumstance that the destabilizing action of cross-diffusion can formally be traced by methods significantly more elementary than those that appear necessary in more complex frameworks of related boundary value problems [10,11], or fully parabolic variants [1216]: nonexistence of global solutions to (1.1) with presupposed suitably favorable smoothness and integrability features, namely, can already be conjectured by means of fairly simple virial-type arguments based on moment evolution [4,1719].

For rigorous confirmations of such reasonings, having available not only a handy theory of local classical solvability, but especially also corresponding conveniently applicable criteria for extensibility, seems of considerable relevance; to the best of our knowledge, however, these issues have not been addressed comprehensively in the literature so far. Indeed, basic existence theories for (1.1) with S = 1 are fairly well developed both in frameworks of weak solutions enjoying certain additional properties related to the evolution of moments and energies [17,20,21], and within concepts of mild solutions that involve various types of requirements on spatial regularity, operating inter alia in spaces of pseudomeasures or in Morrey spaces [2,2225]. Unlike in relatives of (1.1) posed in bounded domains [26], however, settings of genuinely classical solutions seem to have not been explicitly covered in the literature.

Main results. The purpose of the present note is to briefly describe an approach toward local-in-time existence and uniqueness in (1.1) that operates in frameworks of classical solutions, and that artlessly provides some easy-to-handle basic information on behavior near maximal existence times in cases when these are finite. This will be achieved on the basis of an essentially well-established argument involving Banach’s fixed point theorem, in contrast to previous related literature [2] now arranged in a setting of solutions ( u , v ) for which u ( , t ) belongs to BUC ( R n ) L p ( R n ) for fixed times t . Here, as usual, BUC ( R n ) denotes the Banach space of bounded and uniformly continuous functions on R n , and the summability exponent p 1 will be chosen small enough such that the second equation in (1.1) can be fulfilled in a standard manner by defining v through convolution with the corresponding Newtonian kernel. By construction, at this first stage, we will obtain solvability up to some maximal T max , which, if being finite, satisfies u ( , t ) L ( R n ) + u ( , t ) L p ( R n ) as T T max (Lemma 2.2). A refined analysis relying on the local regularity properties of this solution asserted by this first step, however, will afterward enable us to make sure that actually the spatial L norm of any non-global solution must become unbounded, and to hence exclude any situation of blow-up due to a loss of sufficient spatial decay at t = T max , occurring exclusively with respect to L p norms with finite p (Lemma 3.1); apart from that, an independent additional argument will ensure that also v cannot remain bounded near any finite blow-up time (Lemma 3.2).

In summary, this will lead to the following first result of this manuscript, where as usual we let B R ( 0 ) { x R n x < R } for n 2 and R > 0 .

Proposition 1.1

Let n 2 and S R n × n , and suppose that with some p [ 1 , n ) , the function u 0 BUC ( R n ) L p ( R n ) is non-negative. Then, there exist T max ( 0 , ] and a non-negative

(1.2) u C 0 ( [ 0 , T max ) ; BUC ( R n ) ) C 0 ( [ 0 , T max ) ; L p ( R n ) ) C ( R n × ( 0 , T max ) )

such that writing v ( , t ) Γ u ( , t ) , t ( 0 , T max ) , with

(1.3) Γ ( z ) 1 2 π ln z , z R n \ { 0 } , i f n = 2 , 1 n ( n 2 ) B 1 ( 0 ) z 2 n z R n \ { 0 } , i f n 3 ,

we obtain v C ( R n × ( 0 , T max ) ) such that

(1.4) v L loc ( [ 0 , T max ) ; L ( R n ; R n ) ) ,

that ( u , v ) forms a classical solution of (1.1) in R n × ( 0 , T max ) , and that

(1.5) i f T max < , t h e n b o t h limsup t T max u ( , t ) L ( R n ) = a n d limsup t T max v ( , t ) L ( R n ) = .

This solution is uniquely determined in the sense that if T ( 0 , T max ) , and if ( u ^ , v ^ ) is a classical solution of (1.1) in R n × ( 0 , T ) fulfilling u ^ C 0 ( [ 0 , T ] ; BUC ( R n ) ) C 0 ( [ 0 , T ] ; L p ( R n ) ) C 2 , 1 ( R n × ( 0 , T ) ) and v ^ C 2 , 0 ( R n × ( 0 , T ) ) as well as v ^ L ( R n × ( 0 , T ) ; R n ) , then u ^ u in R n × ( 0 , T ) .

Moreover, in the case when p = 1 , we additionally have

(1.6) R n u ( , t ) = R n u 0 f o r a l l t ( 0 , T max ) .

Proposition 1.1 can be used to provide alternative approaches toward blow-up detections in (1.1)[2,4,17,18], under appropriate assumptions on initial concentratedness now yielding solutions known to be bounded and continuous up to some T max ( 0 , ) , and to satisfy limsup t T max u ( , t ) L ( R n ) = due to (1.5); some exemplary applications in this direction can be found detailed in [27]. Moreover, as will also be seen in [27], the second part of the blow-up characterization in (1.5) can be used to preclude the occurrence of explosions especially in some radially symmetric situations, in which suitable conditions on smallness of u 0 ensure L bounds for the functions ( r , t ) r B r ( 0 ) u ( , t ) .

In the second part of the present manuscript, we shall next apply the above general theory in order to derive a result on global solvability in (1.1), unconditional with respect to the size of the initial data, in cases when the matrix S represents cross-diffusion, which is suitably close to repulsion. In planar situations, this will rigorously confirm part of the properties presented and formally addressed in [19]; in fact, no additional requirements on finiteness of second moments, and hence on spatial decay of u 0 , will be needed here, and furthermore, the obtained solutions will be smooth and bounded. In three- or higher-dimensional frameworks, this even seems to be the first statement on global solvability in (1.1) when rotational flux components are involved. Indeed, by making essential use of the refined extensibility criterion in (1.5) we will see in Section 4 that the solutions from Proposition 1.1 actually are global and bounded when S + 1 is suitably small:

Proposition 1.2

Let n 2 . Then, there exists δ > 0 with the property that if

(1.7) S = 1 + S 1 ,

with some S 1 R n × n fulfilling

(1.8) S 1 δ ,

then for any non-negative u 0 BUC ( R n ) L 1 ( R n ) , problem (1.1) admits a global classical solution ( u , v ) with

(1.9) u C 0 ( [ 0 , ) ; BUC ( R n ) ) C 0 ( [ 0 , ) ; L 1 ( R n ) ) C ( R n × ( 0 , ) )

and v ( , t ) = Γ u ( , t ) , t > 0 , uniquely determined in the sense specified in Proposition 1.1. This solution is bounded in the sense that there exists C > 0 satisfying

(1.10) u ( , t ) L ( R n ) C f o r a l l t > 0 .

2 Local existence, uniqueness, non-negativity and mass conservation

As a starting point for our analysis, let us state a basic observation on Newtonian potentials.

Lemma 2.1

Let n 2 , and for z R n \ { 0 } , let Γ be as in (1.3). Then for each p [ 1 , n ) , there exists γ = γ ( p ) > 0 such that

(2.1) ( Γ φ ) L ( R n ) γ φ p f o r a l l φ L ( R n ) L p ( R n ) ,

where we have set

(2.2) φ p φ L ( R n ) + φ L p ( R n ) for φ L ( R n ) L p ( R n ) .

Proof

We decompose

(2.3) ( Γ φ ) ( x ) = x y 1 ( Γ ) ( x y ) φ ( y ) d y + x y > 1 ( Γ ) ( x y ) φ ( y ) d y , x R n ,

and observe that since

( Γ ) ( z ) c 1 z 1 n for all z R n \ { 0 } , with c 1 1 2 π if n = 2 , 1 n B 1 ( 0 ) if n 3 ,

we have

x y 1 ( Γ ) ( x y ) φ ( y ) d y c 1 φ L ( R n ) z 1 z 1 n d z for all x R n .

Moreover, if p > 1 , then by the Hölder inequality,

x y > 1 ( Γ ) ( x y ) φ ( y ) d y c 1 φ L p ( R n ) z > 1 z ( 1 n ) p p 1 d z p 1 p = c 1 ( p 1 ) n B 1 ( 0 ) n p p 1 p φ L p ( R n ) for all x R n ,

while in the case p = 1 , we can simply estimate

x y > 1 ( Γ ) ( x y ) φ ( y ) d y c 1 φ L 1 ( R n ) for all x R n .

Therefore, (2.1) results from (2.3).□

The inequality in (2.1) will form a crucial ingredient in an otherwise fairly straightforward derivation of local solvability in (1.1) on the basis of a contraction mapping argument.

Lemma 2.2

Let n 2 and S R n × n , let p [ 1 , n ) , and let u 0 BUC ( R n ) L p ( R n ) be non-negative. Then, there exist T max ( 0 , ] as well as a function

(2.4) u C 0 ( [ 0 , T max ) ; BUC ( R n ) ) C 0 ( [ 0 , T max ) ; L p ( R n ) ) C ( R n × ( 0 , T max ) ) ,

such that letting

(2.5) v ( , t ) Γ u ( , t ) t ( 0 , T max ) ,

defines a function v C ( R n × ( 0 , T max ) ) , which satisfies

(2.6) v L loc ( [ 0 , T max ) ; L ( R n ; R n ) ) ,

which is such that ( u , v ) solves (1.1) in the classical sense in R n × ( 0 , T max ) , and that

(2.7) i f T max < , t h e n u ( , t ) L ( R n ) + u ( , t ) L p ( R n ) a s t T max .

Proof

With T ( 0 , 1 ) to be specified below, in the Banach space

X C 0 ( [ 0 , T ] ; BUC ( R n ) ) C 0 ( [ 0 , T ] ; L p ( R n ) ) ,

we fix the closed subset

{ u X u X R + 1 } ,

where u X sup t ( 0 , T ) u ( , t ) p for u X , and where R u 0 p . Given u , with Γ as in (1.3), we then let v ( , t ) Γ u ( , t ) for t ( 0 , T ) and define

( F u ) ( , t ) e t Δ u 0 0 t e ( t s ) Δ { u ( , t ) S v ( , s ) } d s , t [ 0 , T ] ,

where ( e t Δ ) t 0 denotes the heat semigroup on R n . Then according to known smoothing properties thereof ([28, Sect. 48.2]), by using Lemma 2.1, it can readily be verified that F u belongs to X for any such u and that there exist positive constants c 1 , c 2 and c 3 such that whenever u ,

( F u ) ( , t ) L ( R n ) u 0 L ( R n ) + c 1 0 t ( t s ) 1 2 u ( , s ) S v ( , s ) L ( R n ) d s u 0 L ( R n ) + c 2 0 t ( t s ) 1 2 u ( , s ) L ( R n ) v ( , s ) L ( R n ) d s u 0 L ( R n ) + c 3 0 t ( t s ) 1 2 u ( , s ) L ( R n ) u ( , s ) p d s u 0 L ( R n ) + c 3 0 t ( t s ) 1 2 u ( , s ) p 2 d s u 0 L ( R n ) + 2 c 3 ( R + 1 ) 2 T 1 2 for all t [ 0 , T ] .

As we can similarly find c i > 0 , i { 4 , 5 , 6 } , such that for any choice of u , we have

( F u ) ( , t ) L p ( R n ) u 0 L p ( R n ) + c 4 0 t ( t s ) 1 2 u ( , s ) S v ( , s ) L p ( R n ) d s u 0 L p ( R n ) + c 5 0 t ( t s ) 1 2 u ( , s ) L p ( R n ) v ( , s ) L ( R n ) d s u 0 L p ( R n ) + c 6 0 t ( t s ) 1 2 u ( , s ) L p ( R n ) u ( , s ) p d s u 0 L p ( R n ) + c 6 0 t ( t s ) 1 2 u ( , s ) p 2 d s u 0 L p ( R n ) + 2 c 6 ( R + 1 ) 2 T 1 2 for all t [ 0 , T ] ,

we infer that if 2 ( c 3 + c 6 ) ( R + 1 ) 2 T 1 2 1 , then

( F u ) ( , t ) p u 0 L ( R n ) + u 0 L p ( R n ) + 1 = R + 1 for all t [ 0 , T ] ,

meaning that F . Since by straightforward adaptation of the above argument it can be shown that if T is further diminished, then also

( F u F u ¯ ) ( , t ) L ( R n ) + ( F u F u ¯ ) ( , t ) L p ( R n ) 1 2 u u ¯ X for all  u , u ¯  and t [ 0 , T ] ,

it follows that if we fix some T = T ( R ) ( 0 , 1 ) suitably small such that F acts as a contractive self-map on .

For its fixed point u , correspondingly existing thanks to the Banach contraction mapping theorem, it can then readily be verified by means of a standard reasoning that the pair ( u , v ) , with v as accordingly introduced above, solves (1.1) in the natural weak sense underlying the analysis in [29], and that hence, according to the well-established bootstrap arguments involving interior parabolic and elliptic Hölder and Schauder estimates [2931], both u and v actually belong to C ( R n × ( 0 , T ) ) and solve (1.1) classically. As the above choice of T = T ( R ) depends on R = u 0 p only, extensibility up to some maximal T max ( 0 , ] fulfilling (2.7) is evident.□

A uniqueness property in the style of that from Proposition 1.1 will be verified by means of an argument based on Duhamel-tape identities satisfied by differences between two solutions. In order to unambiguously guarantee validity of corresponding variation-of-constants representations for arbitrary solutions from the class of functions addressed above, our reasoning in this direction will be arranged in a suitably localized setting. In preparation of this and of several further related arguments to follow, let us fix ζ ( 0 ) C ( R ) such that 0 ζ ( 0 ) 1 on R , ζ ( 0 ) 1 on ( , 0 ) and supp ζ ( 0 ) ( , 1 ) , and for R > 1 , we let

(2.8) ζ R ( x ) ζ ( 0 ) ( x R ) , x R n ,

noting that then ζ R C 0 ( R n ) with

(2.9) 0 ζ R 1 on R n , ζ R 1 in  B R and supp ζ R B R + 1 ,

where we have abbreviated B R B R ( 0 ) R n for R > 1 . Moreover, writing K ζ ( ζ ( 0 ) ) L ( R ) + n ( ζ ( 0 ) ) L ( R ) , we see that

(2.10) ζ R + Δ ζ R K ζ on R n for all R > 1 .

Utilizing this family of cutoff functions, we can now safely derive a statement on uniqueness, which actually is slightly more comprehensive than that announced in Proposition 1.1:

Lemma 2.3

Let n 2 and p [ 1 , n ) , let S R n × n and 0 u 0 BUC ( R n ) L p ( R n ) , and let T max and ( u , v ) be as obtained in Lemma 2.2. Then, ( u , v ) is unique in the sense that whenever T ( 0 , T max ) , u ^ C 0 ( [ 0 , T ] ; BUC ( R n ) ) C 0 ( [ 0 , T ] ; L p ( R n ) ) C 2 , 1 ( R n × ( 0 , T ) ) ; and v ^ C 2 , 0 ( R n × ( 0 , T ) ) are such that v ^ L ( R n × ( 0 , T ) ; R n ) and that ( u ^ , v ^ ) forms a classical solution of (1.1) in R n × ( 0 , T ) , it follows that there exists v ˜ C 2 , 0 ( R n × ( 0 , T ) ) such that Δ v ˜ ( , t ) 0 for all t ( 0 , T ) , and that ( u ^ , v ^ ) = ( u , v + v ˜ ) in R n × ( 0 , T ) .

Proof

According to the regularity features asserted by Lemma 2.2, as well as the current hypothesis, we can find c 1 > 0 and c 2 > 0 such that

(2.11) v ( , t ) L ( R n ) c 1 for all t ( 0 , T ) ,

and that

(2.12) u ^ ( , t ) p c 2 for all t ( 0 , T ) ,

and in order to make sure that these inequalities imply that

(2.13) u ( , t ) u ^ ( , t ) p = 0 for all t ( 0 , T ) ,

we take ( ζ R ) R > 1 from (2.8) and note that then for each R > 1 , thanks to (1.1),

t { ζ R ( u u ^ ) } = ζ R Δ ( u u ^ ) ζ R ( u S v ) + ζ R ( u ^ S v ^ ) = Δ { ζ R ( u u ^ ) } 2 { ( u u ^ ) ζ R } + ( u u ^ ) Δ ζ R { ζ R ( u u ^ ) S v } + ( u u ^ ) ( S v ) ζ R { ζ R u ^ S ( v v ^ ) } + u ^ ( S ( v v ^ ) ) ζ R

holds in the classical pointwise sense in R n × ( 0 , T ) . Since supp ζ R is bounded, we may thus rely on a corresponding variation-of-constants representation to see that again due to known regularization features of the heat semigroup ( e t Δ ) t 0 , there exists c 3 > 0 such that for both choices of q from the set { p , } ,

(2.14) ζ R ( u ( , t ) u ^ ( , t ) ) L q ( R n ) c 3 0 t ( t s ) 1 2 ( u ( , s ) u ^ ( , s ) ) ζ R L q ( R n ) d s + c 3 0 t ( u ( , s ) u ^ ( , s ) ) Δ ζ R L q ( R n ) d s + c 3 0 t ( t s ) 1 2 ζ R ( u ( , s ) u ^ ( , s ) ) S v ( , s ) L q ( R n ) d s + c 3 0 t ( u ( , s ) u ^ ( , s ) ) ( S v ( , s ) ) ζ R L q ( R n ) d s + c 3 0 t ( t s ) 1 2 ζ R u ^ ( , s ) S ( v ( , s ) v ^ ( , s ) ) L q ( R n ) d s + c 3 0 t u ^ ( , s ) ( S ( v ( , s ) v ^ ( , s ) ) ) ζ R L q ( R n ) d s

for all t ( 0 , T ) , because u ( , 0 ) u ^ ( , 0 ) 0 in R n by assumption. Here, thanks to (2.9), (2.10), and the rough estimate 1 T 1 2 ( t s ) 1 2 for 0 < s < t < T ,

(2.15) c 3 0 t ( t s ) 1 2 ( u ( , s ) u ^ ( , s ) ) ζ R L q ( R n ) d s + c 3 0 t ( u ( , s ) u ^ ( , s ) ) Δ ζ R L q ( R n ) d s c 3 K ζ 0 t ( t s ) 1 2 u ( , s ) u ^ ( , s ) L q ( R n ) d s + c 3 K ζ 0 t u ( , s ) u ^ ( , s ) L q ( R n ) d s c 3 K ζ ( 1 + T 1 2 ) 0 t ( t s ) 1 2 u ( , s ) u ^ ( , s ) p d s for all t ( 0 , T ) ,

while a combination of (2.9) and (2.10) with (2.11) shows that, similarly,

(2.16) c 3 0 t ( t s ) 1 2 ζ R ( u ( , s ) u ^ ( , s ) ) S v ( , s ) L q ( R n ) d s + c 3 0 t ( u ( , s ) u ^ ( , s ) ) ( S v ( , s ) ) ζ R L q ( R n ) d s c 3 S 0 t ( t s ) 1 2 u ( , s ) u ^ ( , s ) L q ( R n ) v ( , s ) L ( R n ) d s + c 3 K ζ S 0 t u ( , s ) u ^ ( , s ) L q ( R n ) v ( , s ) L ( R n ) d s c 1 c 3 S ( 1 + K ζ T 1 2 ) 0 t ( t s ) 1 2 u ( , s ) u ^ ( , s ) p d s for all t ( 0 , T ) .

Furthermore, using (2.9) and (2.10) along with (2.12) and Lemma 2.1 we obtain that with γ > 0 taken from the latter we have

c 3 0 t ( t s ) 1 2 ζ R u ^ ( , s ) S ( v ( , s ) v ^ ( , s ) ) L q ( R n ) d s + c 3 0 t u ^ ( , s ) ( S ( v ( , s ) v ^ ( , s ) ) ) ζ R L q ( R n ) d s c 3 S 0 t ( t s ) 1 2 u ^ ( , s ) L q ( R n ) ( v ( , s ) v ^ ( , s ) ) L ( R n ) d s + c 3 K ζ S 0 t u ^ ( , s ) L q ( R n ) ( v ( , s ) v ^ ( , s ) ) L ( R n ) d s c 2 c 3 S ( 1 + K ζ T 1 2 ) 0 t ( t s ) 1 2 ( v ( , s ) v ^ ( , s ) ) L ( R n ) d s c 2 c 3 γ S ( 1 + K ζ T 1 2 ) 0 t ( t s ) 1 2 u ( , s ) u ^ ( , s ) p d s for all t ( 0 , T ) ,

so that from (2.14)–(2.16) we infer that with some c 4 > 0 we have

(2.17) ζ R ( u ( , t ) u ^ ( , t ) ) p c 4 0 t ( t s ) 1 2 u ( , s ) u ^ ( , s ) p d s for all t ( 0 , T ) .

Since ζ R 1 a.e. in R n as R by (2.9), in view of Fatou’s lemma this entails that

u ( , t ) u ^ ( , t ) p c 4 0 t ( t s ) 1 2 u ( , s ) u ^ ( , s ) p d s for all t ( 0 , T ) ,

and thereby, due to the continuity of [ 0 , T ] t u ( , t ) u ^ ( , t ) p guaranteed by the regularity features of u and u ^ , indeed implies (2.13) upon an application of a Gronwall-type inequality [32, Exercise 4 , p. 190].

As thus u ^ u and, in particular, Δ v ^ Δ v , the existence of v ˜ with the claimed properties is evident.□

In quite a straightforward manner, the functions from (2.8) can also be used in the course of a localized testing procedure asserting sign preservation in the first solution components:

Lemma 2.4

Let n 2 and S R n × n , and for p [ 1 , n ) and 0 u 0 BUC ( R n ) L p ( R n ) , let T max ( 0 , ] and ( u , v ) be as obtained in Lemma 2.2. Then,

(2.18) u ( x , t ) 0 f o r a l l x R n a n d t ( 0 , T max ) .

Proof

We fix T ( 0 , T max ) and can then find c 1 ( T ) > 0 and c 2 ( T ) > 0 such that

(2.19) u ( x , t ) c 1 ( T ) and v ( x , t ) c 2 ( T ) for all  x R n  and t ( 0 , T ) .

We moreover take any q > 2 such that q p , and note that then R z z q max { z , 0 } belongs to C 2 ( R ) , whence for each R > 1 we may use (1.1) to see that for all t ( 0 , T ) ,

(2.20) 1 q d d t R n ζ R 2 u q = R n ζ R 2 u q 1 { u u S v }

(2.21) = ( q 1 ) R n ζ R 2 u q 2 u 2 + 2 R n ζ R u q 1 u ζ R ( q 1 ) R n ζ R 2 u q 1 u ( S v ) + 2 R n ζ R u q ( S v ) ζ R .

Here, by Young’s inequality, (2.9) and (2.10), with K ζ as introduced there we have

(2.22) 2 R n ζ R u q 1 u ζ R q 1 2 R n ζ R 2 u q 2 u 2 + 2 q 1 R n u q ζ R 2 q 1 2 R n ζ R 2 u q 2 u 2 + 2 K ζ 2 q 1 B R + 1 \ B R u q for all t ( 0 , T ) ,

whereas combining Young’s inequality with (2.19) shows that

(2.23) ( q 1 ) R n ζ R 2 u q 1 u ( S v ) q 1 2 R n ζ R 2 u q 2 u 2 + q 1 2 R n ζ R 2 u q S v 2 q 1 2 R n ζ R 2 u q 2 u 2 + q 1 2 c 2 2 ( T ) S 2 R n ζ R 2 u q

for all t ( 0 , T ) . As furthermore, again by (2.9), (2.10), and (2.19),

2 R n ζ R u q ( S v ) ζ R 2 c 2 ( T ) S K ζ B R + 1 \ B R u q for all t ( 0 , T ) ,

from (2.20)–(2.23) we obtain that

(2.24) d d t R n ζ R 2 u q c 3 ( T ) R n ζ R 2 u q + c 4 ( T ) B R + 1 \ B R u q for all t ( 0 , T ) ,

where c 3 ( T ) q ( q 1 ) 2 c 2 2 ( T ) S 2 and c 4 ( T ) 2 q K ζ 2 q 1 + 2 q c 2 ( T ) S K ζ . Since u 0 0 , an integration of (2.24) leads to the inequality

(2.25) R n ζ R 2 u q ( , t ) 0 t e c 3 ( T ) ( t s ) c 4 ( T ) B R + 1 \ B R u ( , s ) q d s c 4 ( T ) e c 3 ( T ) T 0 T B R + 1 \ B R u q for all t ( 0 , T ) ,

where using that q p , by means of the dominated convergence theorem we find that

0 T B R + 1 \ B R u q c 1 q p ( T ) 0 T B R + 1 \ B R u p 0 as R

because of the inclusion u L p ( R n × ( 0 , T ) ) implied by Lemma 2.2. Once more relying on (2.9), according to Fatou’s lemma we thus conclude from (2.25) on letting R that

R n u q ( , t ) = 0 for all t ( 0 , T ) ,

and that hence (2.18) results after taking T T max .□

A close relative of the above procedure leads, in its first step, to a localized counterpart of an identity describing the evolution of spatial L q norms. Besides asserting mass conservation in Lemma 2.6, this will form a crucial ingredient in our refinement of the extensibility criterion (3.1) in Lemma 3.1, and will moreover be applied later on in our analysis of the approximately repulsive version of (1.1) addressed in Proposition 1.2 (cf. Lemma 4.1).

Lemma 2.5

Let n 2 , S R n × n , and p [ 1 , n ) , let 0 u 0 BUC ( R n ) L p ( R n ) , and let T max and ( u , v ) as well as ( ζ R ) R > 1 be as in Lemma 2.2and (2.8). Then whenever q 1 and R > 1 ,

(2.26) d d t R n ζ R 2 u q + q ( q 1 ) R n ζ R 2 u q 2 u 2 = 2 q R n ζ R u q 1 u ζ R ( q 1 ) R n ζ R 2 u q ( S v ) + 2 R n ζ R u q ( S v ) ζ R f o r a l l t ( 0 , T max ) .

Proof

This can be verified upon integrating by parts in (1.1) in a straightforward manner.□

Indeed, in the case when in Lemma 2.2 we have p = 1 , constancy of spatial L 1 norms can be obtained as a straightforward consequence of this.

Lemma 2.6

Let n 2 and S R n × n , and let 0 u 0 BUC ( R n ) L 1 ( R n ) . Then, the solution ( u , v ) of (1.1) from Lemma 2.2has the property that R n u ( , t ) = R n u 0 for all t ( 0 , T max ) .

Proof

For T ( 0 , T max ) , we again pick c 1 ( T ) > 0 such that v c 1 ( T ) in R n × ( 0 , T ) , and then obtain that in the identity

(2.27) d d t R n ζ R 2 u = h R ( t ) 2 R n ζ R u ζ R + 2 R n ζ R u ( S v ) ζ R , t ( 0 , T ) ,

valid according to Lemma 2.5, we can again rely on (2.9) and (2.10) to verify that

(2.28) 2 R n ζ R u ( S v ) ζ R c 2 ( T ) B R + 1 \ B R u for all t ( 0 , T ) ,

where c 2 ( T ) 2 c 1 ( T ) K ζ S . In the first summand on the right of (2.27), we first integrate by parts once more to see that thereafter we may employ (2.9) and (2.10) in estimating

(2.29) 2 R n ζ R u ζ R = 2 R n ζ R u Δ ζ R + 2 R n u ζ R 2 c 3 B R + 1 \ B R u for all t ( 0 , T )

with c 3 2 K ζ + 2 K ζ 2 . Now since u L 1 ( R n × ( 0 , T ) ) by Lemma 2.2, in view of the dominated convergence theorem the inequalities in (2.28) and (2.29) ensure that

0 T h R ( t ) d t ( c 2 ( T ) + c 3 ) 0 T B R + 1 \ B R u 0 as R ,

whence (2.27) implies that for each t ( 0 , T ) ,

R n ζ R 2 u ( , t ) R n ζ R 2 u 0 = 0 t h R ( s ) d s 0 as R .

Since (2.9) entails that as R , we have R n ζ R 2 u 0 R n u 0 and R n ζ R 2 u ( , t ) R n u ( , t ) for all t ( 0 , T ) , this establishes (1.6).□

3 Refined criteria for extensibility

Again on the basis of Lemma 2.5, we can now exclude the possibility that the solution from Lemma 2.2 ceases to exist in finite time, but that this is purely due to a lack of sufficient decay at spatial infinity near the blow-up time:

Lemma 3.1

Let n 2 , S R n × n , and p [ 1 , n ) , and given a non-negative u 0 BUC ( R n ) L p ( R n ) , let T max ( 0 , ] and ( u , v ) be as in Lemma 2.2. Then, in fact,

(3.1) e i t h e r T max = , o r limsup t T max u ( , t ) L ( R n ) = .

Proof

If the claim was false, then for some non-negative u 0 BUC ( R n ) L p ( R n ) , we would have T max < but

(3.2) u ( , t ) L ( R n ) M for all t ( 0 , T max )

with some M > 0 . To see that this is impossible, we fix T ( 0 , T max ) and recall (2.6) to pick c 1 ( T ) > 0 such that

(3.3) v ( , t ) L ( R n ) c 1 ( T ) for all t ( 0 , T ) .

In the context of an application of (2.26) to q p , for each R > 1 hence leading to the identity

(3.4) d d t R n ζ R 2 u p + p ( p 1 ) R n ζ R 2 u p 2 u 2 = 2 p R n ζ R u p 1 u ζ R ( p 1 ) R n ζ R 2 u p ( S v ) + 2 R n ζ R u p ( S v ) ζ R , t ( 0 , T max ) ,

this enables us to estimate the rightmost summand on the basis of (2.9) and (2.10) according to

(3.5) 2 R n ζ R u p ( S v ) ζ R 2 c 1 ( T ) K ζ S B R + 1 \ B R u p for all t ( 0 , T ) .

Apart from that, twice employing Young’s inequality shows that for all t ( 0 , T max ) , again by (2.9) and (2.10),

(3.6) 2 p R n ζ R u p 1 u ζ R p ( p 1 ) R n ζ R 2 u p 2 u 2 + p p 1 R n u p ζ R 2 p ( p 1 ) R n ζ R 2 u p 2 u 2 + p K ζ 2 p 1 B R + 1 \ B R u p

and

( p 1 ) R n ζ R 2 u p ( S v ) ( p 1 ) R n ζ R 2 u p + 1 + ( p 1 ) R n ζ R 2 ( S v ) p + 1 ( p 1 ) R n ζ R 2 u p + 1 + ( p 1 ) S p + 1 R n D 2 v p + 1 ,

where in line with a Calderón-Zygmund estimate [33,34], with some c 2 > 0 , we have

( p 1 ) S p + 1 R n D 2 v p + 1 c 2 R n u p + 1 for all t ( 0 , T max ) .

Since thus (3.3) guarantees that

( p 1 ) R n ζ R 2 u p ( S v ) ( p 1 + c 2 ) M R n u p for all t ( 0 , T max ) ,

from (3.4)-(3.6) we obtain that

d d t R n ζ R 2 u p c 3 R n u p + c 4 ( T ) B R + 1 \ B R u p for all t ( 0 , T ) ,

where c 3 ( p 1 + c 2 ) M and c 4 ( T ) 2 c 1 ( T ) K ζ S + p K ζ 2 p 1 . Upon integration, this implies that

(3.7) Ω ζ R 2 u p ( , t ) R n ζ R 2 u 0 p + c 3 0 t R n u p + c 4 ( T ) 0 t B R + 1 \ B R u p for all t ( 0 , T ) ,

where thanks to (2.9) and Beppo Levi’s theorem,

R n ζ R 2 u p ( , t ) R n u p ( , t ) for all t ( 0 , T max ) and R n ζ R 2 u 0 p R n u 0 as R ,

and where due to the dominated convergence theorem,

c 4 ( T ) 0 t B R + 1 \ B R u p 0 for all t ( 0 , T ) as R ,

because u L p ( R n × ( 0 , t ) ) for each t ( 0 , T max ) by Lemma 2.2. In the limit R , from (3.7) we thus infer that

R n u p ( , t ) R n u 0 p + c 3 0 t R n u p for all t ( 0 , T ) ,

whence relying on the continuity of [ 0 , T max ) t R n u p ( , t ) , as entailed by Lemma 2.2, we may invoke a Gronwall lemma to see that

R n u p ( , t ) R n u 0 p e c 3 t for all t ( 0 , T ) .

Noting that c 3 is independent of T , we may take T T max here to obtain a contradiction to (2.7).□

Independently of the latter, by once more employing semigroup estimates, we can finally make sure that signal gradients cannot be bounded near a finite blow-up time:

Lemma 3.2

Let n 2 and S R n × n , let p [ 1 , n ) , and u 0 BUC ( R n ) L p ( R n ) be non-negative, and let T max ( 0 , ] and ( u , v ) be as in Lemma 2.2. Then, it follows that

(3.8) i f T max < , t h e n limsup t T max v ( , t ) L ( R n ) = .

Proof

Let us assume that T max < , but that

(3.9) v ( , t ) L ( R n ) M for all t ( 0 , T max )

with some M > 0 . Then, with ( ζ R ) R > 1 as in (2.8), from (1.1) we obtain that

t ( ζ R u ) = Δ ( ζ R u ) 2 ( u ζ R ) + u Δ ζ R ( ζ R u S v ) + u ( S v ) ζ R in R n × ( 0 , T max ) .

Therefore, drawing on (2.9) and (2.10) and on standard smoothing estimates for the heat semigroup, we find c 1 > 0 such that

ζ R u ( , t ) L ( R n ) ζ R u 0 L ( R n ) + c 1 0 t ( t s ) 1 2 u ( , s ) ζ R L ( R n ) d s + c 1 0 t u ( , s ) Δ ζ R L ( R n ) d s + c 1 0 t ( t s ) 1 2 ζ R u ( , s ) S v ( , s ) L ( R n ) d s + c 1 0 t u ( , s ) ( S v ( , s ) ) ζ R L ( R n ) d s u 0 L ( R n ) + c 1 K ζ 0 t ( t s ) 1 2 u ( , s ) L ( R n ) d s + c 1 K ζ 0 t u ( , s ) L ( R n ) d s + c 1 S M 0 t ( t s ) 1 2 u ( , s ) L ( R n ) d s + c 1 K ζ S M 0 t u ( , s ) L ( R n ) d s u 0 L ( R n ) + c 2 0 t ( t s ) 1 2 u ( , s ) L ( R n ) d s for all t ( 0 , T max ) ,

where c 2 c 1 K ζ + c 1 K ζ T max 1 2 + c 1 S M + c 1 K ζ S M T max 1 2 . Taking R here shows that

u ( , t ) L ( R n ) u 0 L ( R n ) + c 2 0 t ( t s ) 1 2 u ( , s ) L ( R n ) d s for all t ( 0 , T max ) ,

whence again utilizing the Gronwall inequality from [32, Exercise 4 , p. 190] we obtain that

sup t ( 0 , T max ) u ( , t ) L ( R n ) < ,

contrary to what has been asserted by Lemma 3.1.□

Our local theory of classical solvability in (1.1) has thereby been completed:

Proof of Proposition 1.1

We only need to combine Lemma 2.2 with Lemmas 2.4, 2.3, 3.1, 3.2, and 2.6.□

4 Global existence in essentially repulsive versions of (1.1)

As a favorable consequence of Proposition 1.1 and especially of (1.5), verifying Proposition 1.2 essentially reduces to a derivation of L bounds for u under an appropriate assumption on smallness of S + 1 . This will be accomplished in a two-step argument, the first part of which again relies on Lemma 2.5 to obtain L q bounds for arbitrarily large but finite q in such close-to-repulsive frameworks:

Lemma 4.1

Given n 2 and q > 1 , one can find δ 0 ( q ) > 0 with the property that if S 1 R n × n is such that S 1 δ 0 ( q ) , and if u 0 BUC ( R n ) L 1 ( R n ) is non-negative, then there exists C = C ( q , u 0 ) > 0 such that for the solution ( u , v ) of (1.1) with S = 1 + S 1 , as obtained in Proposition 1.1, we have

(4.1) R n u q ( , t ) C f o r a l l t ( 0 , T max ) .

Proof

Again using a Calderón-Zygmund estimate for Newtonian potentials, we fix c 1 = c 1 ( q ) > 0 such that with Γ as in Lemma 2.1,

(4.2) Ω D 2 ( Γ φ ) q + 1 c 1 Ω φ q + 1 for all φ L q + 1 ( R n ) ,

and we let

(4.3) δ 0 ( q ) ( 2 c 1 ) 1 q + 1 .

Then, assuming that S 1 , S , u 0 , and ( u , v ) are as indicated above, we again rely on Lemma 2.5 to see that since

( q 1 ) R n ζ R 2 u q { ( 1 + S 1 ) v } = ( q 1 ) R n ζ R 2 u q Δ v ( q 1 ) R n ζ R 2 u q ( S 1 v ) = ( q 1 ) R n ζ R 2 u q + 1 ( q 1 ) R n ζ R 2 u q ( S 1 v )

for all t ( 0 , T max ) by (1.1), we have

(4.4) d d t R n ζ R 2 u q + q ( q 1 ) R n ζ R 2 u q 2 u 2 + ( q 1 ) R n ζ R 2 u q + 1 = 2 q R n ζ R u q 1 u ζ R ( q 1 ) R n ζ R 2 u q ( S 1 v ) + 2 R n ζ R u q ( S v ) ζ R for all t ( 0 , T max ) .

Here once more by Young’s inequality, (2.9) and (2.10),

(4.5) 2 q R n ζ R u q 1 u ζ R q ( q 1 ) R n ζ R 2 u q 2 u 2 + q K ζ 2 q 1 B R + 1 \ B R u q for all t ( 0 , T max ) ,

and combining (2.9) with (2.10) we find that for each fixed T ( 0 , T max ) ,

(4.6) 2 R n ζ R u q ( S v ) ζ R 2 c 2 ( T ) K ζ S B R + 1 \ B R u q for all t ( 0 , T ) ,

where c 2 ( T ) sup t ( 0 , T ) v ( , t ) L ( R n ) is finite by Proposition 1.1. To estimate the second summand on the right-hand side in (4.4), we now use (4.2) to see that, again due to Young’s inequality and (2.9),

(4.7) ( q 1 ) R n ζ R 2 u q ( S 1 v ) q ( q 1 ) q + 1 R n ζ R 2 u q + 1 + q 1 q + 1 R n ζ R 2 ( S 1 v ) q + 1 q ( q 1 ) q + 1 R n u q + 1 + ( q 1 ) S 1 q + 1 q + 1 R n D 2 v q + 1 q ( q 1 ) q + 1 R n u q + 1 + ( q 1 ) S 1 q + 1 q + 1 c 1 R n u q + 1 q ( q 1 ) q + 1 R n u q + 1 + q 1 2 ( q + 1 ) R n u q + 1 = ( q 1 ) ( 2 q + 1 ) 2 ( q + 1 ) R n u q + 1 for all t ( 0 , T max ) ,

because S 1 q + 1 c 1 δ 0 q + 1 ( q ) c 1 1 2 according to (4.3).

In summary, (4.4)–(4.7) imply that if for T ( 0 , T max ) we let c 3 ( T ) q K ζ 2 q 1 + 2 c 2 ( T ) K ζ S , then

(4.8) d d t R n ζ R 2 u q + ( q 1 ) R n ζ R 2 u q + 1 ( q 1 ) ( 2 q + 1 ) 2 ( q + 1 ) R n u q + 1 + c 3 ( T ) B R + 1 \ B R u q for all t ( 0 , T ) .

Here we note that c 4 ( q 1 ) ( 2 q + 1 ) 2 ( q + 1 ) has the property that c 5 q 1 c 4 q 1 2 ( q + 1 ) is positive, and that thanks to the Hölder inequality, (2.9), (1.6) and Young’s inequality, writing c 6 R n u 0 we have

R n ζ R 2 u q = ζ R 2 q u L q ( R n ) q ζ R 2 q u L q + 1 ( R n ) ( q + 1 ) ( q 1 ) q ζ R 2 q u L 1 ( R n ) 1 q ζ R 2 q + 1 u L q + 1 ( R n ) ( q + 1 ) ( q 1 ) q u L 1 ( R n ) 1 q = c 6 1 q R n ζ R 2 u q + 1 q 1 q c 6 1 q R n ζ R 2 u q + 1 + c 6 1 q for all t ( 0 , T max ) ,

and hence

c 5 R n ζ R 2 u q + 1 c 7 R n ζ R 2 u q c 5 for all t ( 0 , T max ) ,

with c 7 c 5 c 6 1 q . From (4.8) it therefore follows that if T ( 0 , T max ) and R > 1 , and if for t [ 0 , T ] we let

y R ( t ) R n ζ R 2 u q ( , t ) and h R ( t ) c 5 + c 4 R n ( 1 ζ R 2 ) u q + 1 ( , t ) + c 3 ( T ) B R + 1 \ B R u q ( , t ) ,

then

(4.9) y R ( t ) + c 7 y R ( t ) h R ( t ) for all t ( 0 , T ) ,

so that

y R ( t ) y R ( 0 ) e c 7 t + 0 t e c 7 ( t s ) h R ( s ) d s for all t ( 0 , T ) .

Since Beppo Levi’s theorem along with (2.9) ensures that

y R ( t ) R n u q ( , t ) for all t ( 0 , T ) and y R ( 0 ) R n u 0 q as R ,

and that

0 t e c 7 ( t s ) c 4 R n ( 1 ζ R 2 ) u q + 1 ( , s ) d s 0 for all t ( 0 , T ) as R ,

and since by dominated convergence we have

0 t e c 7 ( t s ) c 3 ( T ) B R + 1 \ B R u q ( , s ) d s 0 for all t ( 0 , T ) as R ,

taking R in (4.9) leads to the inequality

R n u q ( , t ) R n u 0 q e c 7 t + 0 t e c 7 ( t s ) c 5 d s R n u 0 q + c 5 c 7 for all t ( 0 , T ) .

As T ( 0 , T max ) was arbitrary, this establishes the claim.□

Once again thanks to heat semigroup regularization, made accessible by regularization using the functions from (2.8), an application of Lemma 4.1 to suitably large q enables us to establish the second step of our reasoning:

Lemma 4.2

Let n 2 . Then, there exists δ > 0 with the following property: Whenever S 1 R n × n satisfies S 1 δ , for any non-negative u 0 BUC ( R n ) L 1 ( R n ) , the solution ( u , v ) of (1.1) with S = 1 + S 1 is global and bounded; that is, in Proposition 1.1we have T max = , and there exists C > 0 such that (1.10) holds.

Proof

We fix any q > n , and with δ 0 ( ) taken from Lemma 4.1, we let δ δ 0 ( q ) . Then assuming that S = 1 + S 1 with some S 1 R n × n fulfilling S 1 δ , and that u 0 BUC ( R n ) L 1 ( R n ) is non-negative, from Lemma 4.1 we obtain c 1 > 0 such that

(4.10) u ( , t ) L q ( R n ) c 1 for all t ( 0 , T max ) .

Recalling from Lemma 2.1 that with γ as provided there we have v ( , t ) L ( R n ) γ u ( , t ) p for all t ( 0 , T max ) , and using Young’s inequality to estimate u ( , t ) L p ( R n ) u ( , t ) L q ( R n ) + u ( , t ) L 1 ( R n ) for all t ( 0 , T max ) , by means of (1.6) and (4.10) we obtain c 2 > 0 such that

(4.11) v ( , t ) L ( R n ) c 2 u ( , t ) L ( R n ) + c 2 for all t ( 0 , T max ) .

Moreover, choosing c 3 > 0 and c 4 > 0 such that in line with known smoothing properties of the heat equation we have

(4.12) e t Δ φ L ( R n ) c 3 t n 2 q φ L q ( R n ) for all t > 0 and each φ C 0 ( R n ) ,

and that

(4.13) e t Δ φ L ( R n ) c 4 t 1 2 n 2 q φ L q ( R n ) for all  t > 0  and each φ C 0 ( R n ; R n ) ,

we use that q > n to fix τ ( 0 , 1 ] suitably small such that

(4.14) c 1 c 2 c 4 S 1 2 n 2 q τ 1 2 n 2 q 1 4 and c 1 c 2 c 3 K ζ S 1 n 2 q τ 1 n 2 q 1 4 .

For arbitrary T ( 0 , T max ) , we then estimate the finite number

(4.15) M ( T ) sup t ( 0 , T ) u ( , t ) L ( R n )

by relying on a Duhamel representation associated with the identity

t ( ζ R u ) = Δ ( ζ R u ) 2 ( u ζ R ) + u Δ ζ R ( ζ R u S v ) + u ( S v ) ζ R ,

valid for each R > 1 according to (1.1) (cf. also the proof of Lemma 3.2). With τ as above, namely, this implies that thanks to (4.12) and (4.13),

(4.16) ζ R u ( , t ) L ( R n ) = e [ t ( t τ ) + ] Δ { ζ R u ( , ( t τ ) + ) } 2 ( t τ ) + t e ( t s ) Δ { u ( , s ) ζ R } d s + ( t τ ) + t e ( t s ) Δ { u ( , s ) Δ ζ R } d s ( t τ ) + t e ( t s ) Δ { ζ R u ( , s ) S v ( , s ) } d s + ( t τ ) + t e ( t s ) Δ { u ( , s ) ( S v ( , s ) ) ζ R } d s L ( R n ) e [ t ( t τ ) + ] Δ { ζ R u ( , ( t τ ) + ) } L ( R n ) + 2 c 4 ( t τ ) + t ( t s ) 1 2 n 2 q u ( , s ) ζ R L q ( R n ) d s + c 3 ( t τ ) + t ( t s ) n 2 q u ( , s ) Δ ζ R L q ( R n ) d s + c 4 ( t τ ) + t ( t s ) 1 2 n 2 q ζ R u ( , s ) S v ( , s ) L q ( R n ) d s + c 3 ( t τ ) + t ( t s ) n 2 q u ( , s ) ( S v ( , s ) ) ζ R L q ( R n ) d s for all t ( 0 , T ) .

Here if t ( τ , T ) , then combining (4.12) with (4.10) shows that

(4.17) e [ t ( t τ ) + ] Δ { ζ R u ( , ( t τ ) + ) } L ( R n ) = e τ Δ { ζ R u ( , t τ ) } L ( R n ) c 3 τ n 2 q ζ R u ( , t τ ) L q ( R n ) c 1 c 3 τ n 2 q ,

because ζ R 1 ; if t ( 0 , τ ] , however, then by a simple application of the maximum principle we obtain that

(4.18) e [ t ( t τ ) + ] Δ { ζ R u ( , ( t τ ) + ) } L ( R n ) = e t Δ { ζ R u 0 } L ( R n ) ζ R u 0 L ( R n ) u 0 L ( R n ) .

Apart from that, using (2.9) and (2.10) together with (4.10) and the inequalities q > n and τ 1 we can estimate

(4.19) 2 c 4 ( t τ ) + t ( t s ) 1 2 n 2 q u ( , s ) ζ R L q ( R n ) d s 2 c 1 c 4 K ζ ( t τ ) + t ( t s ) 1 2 n 2 q d s 2 c 1 c 4 K ζ 0 1 σ 1 2 n 2 q d σ = 2 c 1 c 4 K ζ 1 2 n 2 q for all t ( 0 , T )

and, similarly,

(4.20) c 3 ( t τ ) + t ( t s ) n 2 q u ( , s ) Δ ζ R L q ( R n ) d s c 1 c 3 K ζ ( t τ ) + t ( t s ) n 2 q d s c 1 c 3 K ζ 1 n 2 q for all t ( 0 , T ) .

In view of the definition of M ( T ) and our restriction in (4.14), finally, the two last summands in (4.16) can be controlled, again on the basis of (2.9), (2.10), and (4.10), and of (4.11) and (4.15), according to

(4.21) c 4 ( t τ ) + t ( t s ) 1 2 n 2 q ζ R u ( , s ) S v ( , s ) L q ( R n ) d s c 4 S ( t τ ) + t ( t s ) 1 2 n 2 q u ( , s ) L q ( R n ) v ( , s ) L ( R n ) d s c 1 c 2 c 4 S ( M ( T ) + 1 ) ( t τ ) + t ( t s ) 1 2 n 2 q d s c 1 c 2 c 4 S 1 2 n 2 q τ 1 2 n 2 q M ( T ) + c 1 c 2 c 4 S 1 2 n 2 q 1 4 M ( T ) + c 1 c 2 c 4 S 1 2 n 2 q for all t ( 0 , T )

and, in much the same fashion,

(4.22) c 3 ( t τ ) + t ( t s ) n 2 q u ( , s ) ( S v ( , s ) ) ζ R L q ( R n ) d s c 1 c 2 c 3 K ζ S ( M ( T ) + 1 ) ( t τ ) + t ( t s ) n 2 q d s c 1 c 2 c 3 K ζ S 1 n 2 q τ 1 n 2 q M ( T ) + c 1 c 2 c 3 K ζ S 1 n 2 q 1 4 M ( T ) + c 1 c 2 c 3 K ζ S 1 n 2 q for all t ( 0 , T ) ,

again because τ 1 . A combination of (4.16) with (4.17)–(4.22) now shows that if we let

c 5 c 1 c 3 τ n 2 q + u 0 L ( R n ) + 2 c 1 c 4 K ζ 1 2 n 2 q + c 1 c 3 K ζ 1 n 2 q + c 1 c 2 c 4 S 1 2 n 2 q + c 1 c 2 c 3 K ζ S 1 n 2 q ,

then whenever T ( 0 , T max ) and R > 1 ,

ζ R u ( , t ) L ( R n ) 1 2 M ( T ) + c 5 for all t ( 0 , T ) .

In the limit R , this yields the inequality

u ( , t ) L ( R n ) 1 2 M ( T ) + c 5 for all t ( 0 , T ) ,

from which it follows that M ( T ) 1 2 M ( T ) + c 5 and hence

u ( , t ) L ( R n ) 2 c 5 for all t ( 0 , T max ) ,

as T ( 0 , T max ) was arbitrary. In light of Proposition 1.1 and especially the extensibility criterion (1.5) therein, this shows that indeed ( u , v ) exists globally and possesses the claimed boundedness property.□

Our final result on global existence and boundedness in (1.1) for small perturbations S of 1 has thereby been achieved already:

Proof of Proposition 1.2

The statement is an immediate consequence of Proposition 1.1, Lemmas 4.1 and 4.2.□

Acknowledgement

The author acknowledges support of the Deutsche Forschungsgemeinschaft (Project No. 462888149).

  1. Conflict of interest: The author states no conflict of interest.

  2. Data availability statement: Data sharing is not applicable to this article as no datasets were generated or analyzed during this study.

References

[1] S. Chandrasekhar, Principles of Stellar Dynamics, University of Chicago Press, Chicago, 1942. Search in Google Scholar

[2] P. Biler, Singularities of Solutions to Chemotaxis Systems, Series in Mathematics and Life Sciences, De Gruyter, Berlin, 2020. 10.1515/9783110599534Search in Google Scholar

[3] T. Suzuki, Free Energy and Self-Interacting Particles, Progress in Nonlinear Differential Equations and Their Applications Vol. 62, Birkhäuser, Boston, 2005. 10.1007/0-8176-4436-9Search in Google Scholar

[4] P. Biler and T. Nadzieja, Existence and nonexistence of solutions for a model of gravitational interaction of particles I, Colloq. Math. 66 (1993), 319–334. 10.4064/cm-66-2-319-334Search in Google Scholar

[5] E. F. Keller and L. A. Segel, Initiation of slime mold aggregation viewed as an instability, J. Theoret. Biol. 26 (1970), 399–415. 10.1016/0022-5193(70)90092-5Search in Google Scholar PubMed

[6] D. Horstmann, From 1970 until present: the Keller-Segel model in chemotaxis and its consequences. I, Jahresber. Deutsch. Math.-Verein. 105 (2003), 103–165. Search in Google Scholar

[7] T. Senba, Blowup behavior of radial solutions to Jäger-Luckhaus system in high dimensional domains, Funkcial. Ekvac. 48 (2005), 247–271. 10.1619/fesi.48.247Search in Google Scholar

[8] T. Senba, A fast blow-up solution to an elliptic-parabolic system related to chemotaxis, Adv. Differential Equations 1 (2006), 981–1030. 10.57262/ade/1355867610Search in Google Scholar

[9] C. Xue and H. G. Othmer, Multiscale models of taxis-driven patterning in bacterial populations, SIAM J. Appl. Math. 70 (2009), 133–167. 10.1137/070711505Search in Google Scholar PubMed PubMed Central

[10] T. Nagai, Blow-up of radially symmetric solutions to a chemotaxis system, Adv. Math. Sci. Appl. 5 (1995), 581–601. Search in Google Scholar

[11] T. Nagai, Blowup of nonradial solutions to parabolic–elliptic systems modeling chemotaxis in two-dimensional domains, J. Inequal. Appl. 6 (2001), 37–55. 10.1155/S1025583401000042Search in Google Scholar

[12] M. A. Herrero and J. J. L. Velázquez, A blow-up mechanism for a chemotaxis model, Ann. Sc. Norm. Super. Pisa Cl. Sci. 24 (1997), 633–683. Search in Google Scholar

[13] M. Winkler, Finite-time blow-up in the higher-dimensional parabolic-parabolic Keller-Segel system, J. Math. Pures Appl. 100 (2013), 748–767, arXiv:1112.4156v1. 10.1016/j.matpur.2013.01.020Search in Google Scholar

[14] T. Senba and T. Suzuki, Parabolic system of chemotaxis: blowup in a finite and the infinite time, Methods Appl. Anal. 8 (2001), 349–368. 10.4310/MAA.2001.v8.n2.a9Search in Google Scholar

[15] T. Nagai, T. Senba, and T. Suzuki, Chemotactic collapse in a parabolic system of mathematical biology, Hiroshima Math. J. 30 (2000), 463–497. 10.32917/hmj/1206124609Search in Google Scholar

[16] M. Winkler, Single-point blow-up in the Cauchy problem for the higher-dimensional Keller-Segel system, Nonlinearity 33 (2020), 5007. 10.1088/1361-6544/ab9247Search in Google Scholar

[17] A. Blanchet, J. Dolbeault, and B. Perthame, Two-dimensional Keller-Segel model: optimal critical mass and qualitative properties of the solutions, Electron. J. Differential Equations 2006 (2006), no. 44, 1–33. Search in Google Scholar

[18] P. Biler and T. Nadzieja, Growth and accretion of mass in an astrophysical model II, Appl. Math. (Warsaw) 23 (1995), 351–361. 10.4064/am-23-3-351-361Search in Google Scholar

[19] E. Espejo and H. Wu, Optimal critical mass for the two-dimensional Keller-Segel model with rotational flux terms, Commun. Math. Sci. 18 (2020), 379–394. 10.4310/CMS.2020.v18.n2.a5Search in Google Scholar

[20] B. Perthame, Transport Equations in Biology, Birkhäuser, Basel, 2007. 10.1007/978-3-7643-7842-4Search in Google Scholar

[21] A. Blanchet, J. A. Carrillo, and N. Masmoudi, Infinite time aggregation for the critical Patlak-Keller-Segel model in R2, Commun. Pure Appl. Math. 61 (2008), 1449–1481. 10.1002/cpa.20225Search in Google Scholar

[22] D. Wei, Global well-posedness and blow-up for the 2-D Patlak-Keller-Segel equation, J. Funct. Anal. 274 (2018), 388–401. 10.1016/j.jfa.2017.10.019Search in Google Scholar

[23] P. Biler, The Cauchy problem and self-similar solutions for a nonlinear parabolic equation, Studia Math. 114 (1995), 181–205. 10.4064/sm-114-2-181-205Search in Google Scholar

[24] A. Raczyński, Stability property of the two-dimensional Keller-Segel model, Asymptot. Anal. 61 (2009), 35–59. 10.3233/ASY-2008-0907Search in Google Scholar

[25] P. Biler, M. Cannone, I. A. Guerra, and G. Karch, Global regular and singular solutions for a model of gravitating particles, Math. Ann. 330 (2004), 693–708. 10.1007/s00208-004-0565-7Search in Google Scholar

[26] N. Bellomo, A. Bellouquid, Y. Tao, and 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. 10.1142/S021820251550044XSearch in Google Scholar

[27] M. Winkler, Radial data with slow spatial decay in the Cauchy problem for the parabolic–elliptic Keller-Segel system, preprint. Search in Google Scholar

[28] P. Quittner and P. Souplet, Superlinear Parabolic Problems. Blow-up, Global Existence and Steady States, Birkhäuser, Basel, 2007. Search in Google Scholar

[29] O. A. Ladyzenskaja, V. A. Solonnikov, and N. N. Ural’ceva, Linear and Quasi-linear Equations of Parabolic Type, American Mathematical Society, Translations of Mathematical Monographs, Vol. 23, Providence, 1968. 10.1090/mmono/023Search in Google Scholar

[30] M. M. Porzio and V. Vespri, Holder estimates for local solutions of some doubly nonlinear degenerate parabolic equations, J. Differential Equations 103 (1993), 146–178. 10.1006/jdeq.1993.1045Search in Google Scholar

[31] D. Gilbarg and N. S. Trudinger, Elliptic Partial Differential Equations of Second Order, Springer-Verlag, Berlin, 2001. 10.1007/978-3-642-61798-0Search in Google Scholar

[32] D. Henry, Geometric Theory of Semilinear Parabolic Equations, Lecture Notes in Mathematics, Vol. 840, Springer-Verlag, Berlin, 1981. 10.1007/BFb0089647Search in Google Scholar

[33] E. M. Stein, Singular Integrals and Differentiability Properties of Functions, Vol. 30, Princeton University Press, Princeton, 1970. 10.1515/9781400883882Search in Google Scholar

[34] J. Bedrossian, N. Rodriguez, and A. L. Bertozzi, Local and global well-posedness for aggregation equations and Patlak-Keller-Segel models with degenerate diffusion, Nonlinearity 24 (2011), 1683–1714. 10.1088/0951-7715/24/6/001Search in Google Scholar

Received: 2023-01-18
Revised: 2023-03-28
Accepted: 2023-03-28
Published Online: 2023-05-11

© 2023 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Special Issue on Future Directions of Further Developments in Mathematics
  2. What will the mathematics of tomorrow look like?
  3. On H 2-solutions for a Camassa-Holm type equation
  4. Classical solutions to Cauchy problems for parabolic–elliptic systems of Keller-Segel type
  5. Control of multi-agent systems: Results, open problems, and applications
  6. Logical perspectives on the foundations of probability
  7. Subharmonic solutions for a class of predator-prey models with degenerate weights in periodic environments
  8. A non-smooth Brezis-Oswald uniqueness result
  9. Luenberger compensator theory for heat-Kelvin-Voigt-damped-structure interaction models with interface/boundary feedback controls
  10. Special Issue on Fractional Problems with Variable-Order or Variable Exponents (Part II)
  11. Positive solution for a nonlocal problem with strong singular nonlinearity
  12. Analysis of solutions for the fractional differential equation with Hadamard-type
  13. Hilfer proportional nonlocal fractional integro-multipoint boundary value problems
  14. A comprehensive review on fractional-order optimal control problem and its solution
  15. The θ-derivative as unifying framework of a class of derivatives
  16. Review Articles
  17. On the use of L-functionals in regression models
  18. Minimal-time problems for linear control systems on homogeneous spaces of low-dimensional solvable nonnilpotent Lie groups
  19. Regular Articles
  20. Existence and multiplicity of solutions for a new p(x)-Kirchhoff problem with variable exponents
  21. An extension of the Hermite-Hadamard inequality for a power of a convex function
  22. Existence and multiplicity of solutions for a fourth-order differential system with instantaneous and non-instantaneous impulses
  23. Relay fusion frames in Banach spaces
  24. Refined ratio monotonicity of the coordinator polynomials of the root lattice of type Bn
  25. On the uniqueness of limit cycles for generalized Liénard systems
  26. A derivative-Hilbert operator acting on Dirichlet spaces
  27. Scheduling equal-length jobs with arbitrary sizes on uniform parallel batch machines
  28. Solutions to a modified gauged Schrödinger equation with Choquard type nonlinearity
  29. A symbolic approach to multiple Hurwitz zeta values at non-positive integers
  30. Some results on the value distribution of differential polynomials
  31. Lucas non-Wieferich primes in arithmetic progressions and the abc conjecture
  32. Scattering properties of Sturm-Liouville equations with sign-alternating weight and transmission condition at turning point
  33. Some results for a p(x)-Kirchhoff type variation-inequality problems in non-divergence form
  34. Homotopy cartesian squares in extriangulated categories
  35. A unified perspective on some autocorrelation measures in different fields: A note
  36. Total Roman domination on the digraphs
  37. Well-posedness for bilevel vector equilibrium problems with variable domination structures
  38. Binet's second formula, Hermite's generalization, and two related identities
  39. Non-solid cone b-metric spaces over Banach algebras and fixed point results of contractions with vector-valued coefficients
  40. Multidimensional sampling-Kantorovich operators in BV-spaces
  41. A self-adaptive inertial extragradient method for a class of split pseudomonotone variational inequality problems
  42. Convergence properties for coordinatewise asymptotically negatively associated random vectors in Hilbert space
  43. Relating the super domination and 2-domination numbers in cactus graphs
  44. Compatibility of the method of brackets with classical integration rules
  45. On the inverse Collatz-Sinogowitz irregularity problem
  46. Positive solutions for boundary value problems of a class of second-order differential equation system
  47. Global analysis and control for a vector-borne epidemic model with multi-edge infection on complex networks
  48. Nonexistence of global solutions to Klein-Gordon equations with variable coefficients power-type nonlinearities
  49. On 2r-ideals in commutative rings with zero-divisors
  50. A comparison of some confidence intervals for a binomial proportion based on a shrinkage estimator
  51. The construction of nuclei for normal constituents of Bπ-characters
  52. Weak solution of non-Newtonian polytropic variational inequality in fresh agricultural product supply chain problem
  53. Mean square exponential stability of stochastic function differential equations in the G-framework
  54. Commutators of Hardy-Littlewood operators on p-adic function spaces with variable exponents
  55. Solitons for the coupled matrix nonlinear Schrödinger-type equations and the related Schrödinger flow
  56. The dual index and dual core generalized inverse
  57. Study on Birkhoff orthogonality and symmetry of matrix operators
  58. Uniqueness theorems of the Hahn difference operator of entire function with a Picard exceptional value
  59. Estimates for certain class of rough generalized Marcinkiewicz functions along submanifolds
  60. On semigroups of transformations that preserve a double direction equivalence
  61. Positive solutions for discrete Minkowski curvature systems of the Lane-Emden type
  62. A multigrid discretization scheme based on the shifted inverse iteration for the Steklov eigenvalue problem in inverse scattering
  63. Existence and nonexistence of solutions for elliptic problems with multiple critical exponents
  64. Interpolation inequalities in generalized Orlicz-Sobolev spaces and applications
  65. General Randić indices of a graph and its line graph
  66. On functional reproducing kernels
  67. On the Waring-Goldbach problem for two squares and four cubes
  68. Singular moduli of rth Roots of modular functions
  69. Classification of self-adjoint domains of odd-order differential operators with matrix theory
  70. On the convergence, stability and data dependence results of the JK iteration process in Banach spaces
  71. Hardy spaces associated with some anisotropic mixed-norm Herz spaces and their applications
  72. Remarks on hyponormal Toeplitz operators with nonharmonic symbols
  73. Complete decomposition of the generalized quaternion groups
  74. Injective and coherent endomorphism rings relative to some matrices
  75. Finite spectrum of fourth-order boundary value problems with boundary and transmission conditions dependent on the spectral parameter
  76. Continued fractions related to a group of linear fractional transformations
  77. Multiplicity of solutions for a class of critical Schrödinger-Poisson systems on the Heisenberg group
  78. Approximate controllability for a stochastic elastic system with structural damping and infinite delay
  79. On extremal cacti with respect to the first degree-based entropy
  80. Compression with wildcards: All exact or all minimal hitting sets
  81. Existence and multiplicity of solutions for a class of p-Kirchhoff-type equation RN
  82. Geometric classifications of k-almost Ricci solitons admitting paracontact metrices
  83. Positive periodic solutions for discrete time-delay hematopoiesis model with impulses
  84. On Hermite-Hadamard-type inequalities for systems of partial differential inequalities in the plane
  85. Existence of solutions for semilinear retarded equations with non-instantaneous impulses, non-local conditions, and infinite delay
  86. On the quadratic residues and their distribution properties
  87. On average theta functions of certain quadratic forms as sums of Eisenstein series
  88. Connected component of positive solutions for one-dimensional p-Laplacian problem with a singular weight
  89. Some identities of degenerate harmonic and degenerate hyperharmonic numbers arising from umbral calculus
  90. Mean ergodic theorems for a sequence of nonexpansive mappings in complete CAT(0) spaces and its applications
  91. On some spaces via topological ideals
  92. Linear maps preserving equivalence or asymptotic equivalence on Banach space
  93. Well-posedness and stability analysis for Timoshenko beam system with Coleman-Gurtin's and Gurtin-Pipkin's thermal laws
  94. On a class of stochastic differential equations driven by the generalized stochastic mixed variational inequalities
  95. Entire solutions of two certain Fermat-type ordinary differential equations
  96. Generalized Lie n-derivations on arbitrary triangular algebras
  97. Markov decision processes approximation with coupled dynamics via Markov deterministic control systems
  98. Notes on pseudodifferential operators commutators and Lipschitz functions
  99. On Graham partitions twisted by the Legendre symbol
  100. Strong limit of processes constructed from a renewal process
  101. Construction of analytical solutions to systems of two stochastic differential equations
  102. Two-distance vertex-distinguishing index of sparse graphs
  103. Regularity and abundance on semigroups of partial transformations with invariant set
  104. Liouville theorems for Kirchhoff-type parabolic equations and system on the Heisenberg group
  105. Spin(8,C)-Higgs pairs over a compact Riemann surface
  106. Properties of locally semi-compact Ir-topological groups
  107. Transcendental entire solutions of several complex product-type nonlinear partial differential equations in ℂ2
  108. Ordering stability of Nash equilibria for a class of differential games
  109. A new reverse half-discrete Hilbert-type inequality with one partial sum involving one derivative function of higher order
  110. About a dubious proof of a correct result about closed Newton Cotes error formulas
  111. Ricci ϕ-invariance on almost cosymplectic three-manifolds
  112. Schur-power convexity of integral mean for convex functions on the coordinates
  113. A characterization of a ∼ admissible congruence on a weakly type B semigroup
  114. On Bohr's inequality for special subclasses of stable starlike harmonic mappings
  115. Properties of meromorphic solutions of first-order differential-difference equations
  116. A double-phase eigenvalue problem with large exponents
  117. On the number of perfect matchings in random polygonal chains
  118. Evolutoids and pedaloids of frontals on timelike surfaces
  119. A series expansion of a logarithmic expression and a decreasing property of the ratio of two logarithmic expressions containing cosine
  120. The 𝔪-WG° inverse in the Minkowski space
  121. Stability result for Lord Shulman swelling porous thermo-elastic soils with distributed delay term
  122. Approximate solvability method for nonlocal impulsive evolution equation
  123. Construction of a functional by a given second-order Ito stochastic equation
  124. Global well-posedness of initial-boundary value problem of fifth-order KdV equation posed on finite interval
  125. On pomonoid of partial transformations of a poset
  126. New fractional integral inequalities via Euler's beta function
  127. An efficient Legendre-Galerkin approximation for the fourth-order equation with singular potential and SSP boundary condition
  128. Eigenfunctions in Finsler Gaussian solitons
  129. On a blow-up criterion for solution of 3D fractional Navier-Stokes-Coriolis equations in Lei-Lin-Gevrey spaces
  130. Some estimates for commutators of sharp maximal function on the p-adic Lebesgue spaces
  131. A preconditioned iterative method for coupled fractional partial differential equation in European option pricing
  132. A digital Jordan surface theorem with respect to a graph connectedness
  133. A quasi-boundary value regularization method for the spherically symmetric backward heat conduction problem
  134. The structure fault tolerance of burnt pancake networks
  135. Average value of the divisor class numbers of real cubic function fields
  136. Uniqueness of exponential polynomials
  137. An application of Hayashi's inequality in numerical integration
Downloaded on 5.10.2025 from https://www.degruyterbrill.com/document/doi/10.1515/math-2022-0578/html
Scroll to top button