Next Article in Journal
Joint Resource Allocation for Frequency-Domain Artificial Noise Assisted Multiuser Wiretap OFDM Channels with Finite-Alphabet Inputs
Previous Article in Journal
A Numerical Solution of Fredholm Integral Equations of the Second Kind Based on Tight Framelets Generated by the Oblique Extension Principle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Priori Estimates for a Nonlinear System with Some Essential Symmetrical Structures

1
Department of Mathematics, Jincheng College of Sichuan University, Chengdu 611731, China
2
College of Mathematics, Sichuan University, Chengdu 610064, China
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(7), 852; https://doi.org/10.3390/sym11070852
Submission received: 18 May 2019 / Revised: 17 June 2019 / Accepted: 20 June 2019 / Published: 2 July 2019

Abstract

:
In this paper, we are concerned with a nonlinear system containing some essential symmetrical structures (e.g., cross-diffusion) in the two-dimensional setting, which is proposed to model the biological transport networks. We first provide an a priori blow-up criterion of strong solution of the corresponding Cauchy problem. Based on this, we also establish a priori upper bounds to strong solution for all positive times.

1. Introduction and Main Results

Driven by the need to understand the biological transportation networks (for instance, leaf venation in plants, angiogenesis of blood vessels and neural networks which transport electric charge), biologists and physicists have expressed great interest in investigating the qualitative properties of network structures in the last few decades (see for instance [1,2,3,4,5] and references therein). Recently, Hu and Cai [6] introduced a purely local dynamic adaptation model based on mechanical laws on a graph, which was extended to a continuum one in [7,8,9] that was subsequently studied in the series of papers [10,11,12,13,14,15]. This continuum model, posed in spatial domain Ω, can be read as
t m D 2 Δ m c 2 ( m · p ) p + | m | 2 ( γ 1 ) m = 0 , x Ω , t > 0 , · I + m m p = S , x Ω , t > 0 .
The unknown function p = p ( t , x ) denotes a scalar pressure of the fluid transported within the network which satisfies Darcy’s type equation due to Darcy’s law for slow flow in the network being valid, and thus p represents a driving force for the evolution of the vector-valued conductance m = m 1 ( t , x ) , , m n ( t , x ) that describes the dynamics of networks by using a reaction–diffusion equation consisting of three different mechanisms—pressure effect, diffusion (representing microscopic Brownian process) and an algebraic conductance-relaxation; the given function S = S ( x ) models the sources and sinks. Values of the parameters D 2 > 0 (diffusion coefficient), c 2 > 0 (activation parameter) and γ R (relaxation exponent) are determined by the particular physical applications. For instance, we get from the known experimental studies (see [6,16] and ([8], Section 2) for details) that γ = 1 2 can be used to describe blood vessel systems in the human body and that γ = 1 corresponds to leaf venation.
From a mathematical perspective, system (1) exhibits two rather peculiar nonlinear structures: c 2 ( m · p ) p in system (1)1 and · m m p in (1)2, which may result in yielding several difficulties in the mathematical analysis. For instance, the absence of a priori L -bound for m may cause the elliptic coefficients in (1)2 to be singular, and thus the solution of (1)2 is too weak to control the nonlinear term c 2 ( m · p ) p in (1)1. These challenges make the cross-diffusion system (1) of interest. We would like to mention that, when system (1) is posed in a bounded domain Ω , it should be supplemented with the initial condition
m ( 0 , x ) = m 0 ( x ) in Ω
and the homogeneous Dirichlet boundary conditions
p ( t , x ) = 0 and m ( t , x ) = 0 on Ω .
The initial-boundary value problems (1)–(3) has attracted a lot of interest already—see, for instance, ref. [10] for the existence of global weak solution and of local mild solution when γ 1 that was extended to the case of γ 1 2 in [11], ref. [13] for the partial regularity of weak solution, ref. [14] for the regularity of stationary weak solution in two space dimensions and ref. [15] for the existence of local classical large-data solution and of global classical small-data solution.
To our best knowledge, the Cauchy problems (1) and (2) was only studied in [17] when Ω = R 3 , in which local existence as well as blow-up criterion for large initial data and global existence for small initial data were established based on H ˙ ( R 3 ) being a Hilbert space; however, there is no result on the Cauchy problems (1) and (2) on R 2 . Our goal is to fill this gap. It is worthwhile pointing out that the important problem of the mathematical theory concerning the Cauchy problems (1) and (2) is whether or not the global in time smooth solution exists for any prescribed smooth initial data. Hence, in the absence of a global well-posedness theory, establishing a priori estimates is of major importance for both theoretical and practical purposes, which motivates us to investigate an a priori blow-up criterion and a priori upper bounds of strong solutions to system (1) and (2) with Ω = R 2 . We would like remark that, in contrast to system (1) and (2) with Ω = R 3 , H ˙ ( R 2 ) not being a Hilbert space poses obstacles with adopting the same strategies in [17] to study the Cauchy problems (1) and (2) with Ω = R 2 . To overcome this challenge, we will derive the a priori L -bound for m to rule out the degeneracy and develop some new dissipation mechanisms hidden in the system (1) and (2) by fully utilizing some essential symmetrical structures of system when Ω = R 2 .
Before stating our main results, we need to layout some notations. H k ( R 2 ) and L q ( R 2 ) denote the usual Hilbert space and Lebesgue space with the norms · H k ( R 2 ) and · L q ( R 2 ) (or · H k and · L q for short), respectively. The functions in these spaces are usually understood to be real valued. If no confusion is likely, we shall use the same notation for similar spaces of vector-valued functions and of matrix-valued functions, for instance,
m L q = i = 1 2 m i L q and m L q = i , j = 1 2 j m i L q .
Throughout this paper, we will use A B to denote A C B for some uniform constant C > 0 . Unless specified, the values of the constants may vary line by line according to the context.
With the aforementioned notations, we now state our first result as follows.
Theorem 1.
Let γ 1 , m 0 H 2 ( R 2 ) and w i S L 2 ( R 2 ) , i = 0 , 1 , 2 , where w ( x ) : = ( 1 + | x | ) 1 + ln ( 1 + | x | ) . Assume that there exist a time T > 0 and a strong solution ( m , p ) to systems (1) and (2) on R 2 × [ 0 , T ] satisfying
m C [ 0 , T ] ; H 2 ( R 2 ) , m L 2 [ 0 , T ] ; H 2 ( R 2 ) and p C [ 0 , T ] ; H 2 ( R 2 ) .
Then, there exists a maximal time T max ( T , ] such that ( m , p ) is a strong solution of the Cauchy problem (1) and (2) on R 2 × [ 0 , T max ) . Moreover, if T max < , it holds that
lim sup t T max m ( t ) H 2 ( R 2 ) =
if and only if
0 T max m ( t ) L q ( R 2 ) β d t = ,
for any q ( 2 , + ] , where
β = max 4 q q 2 , 4 q ( γ 2 ) 2 q 3 .
Remark 1.
The blow-up criterion (4) in Theorem 1 is similar to the blow-up criterion of the strong solution of the Cauchy problem (1) and (2) on R 3 [17] and to the Serrin-type criterion of Leray–Hopf weak solution of the Navier–Stokes equations [18].
Invoking the blow-up criterion (4) obtained in Theorem 1, we can further present that the local strong solution of the Cauchy problem (1) and (2) can be extended to a global one.
Theorem 2.
Assume that all conditions in Theorem 1 hold. Then, there exists a positive constant h ( 4 ) (see (32) below) such that, if D h ( 4 ) , then the local strong solution of the Cauchy problem (1) and (2) with Ω = R 2 is global indeed in the sense that
m ( t ) H 2 ( R 2 ) + p H 2 ( R 2 ) + 0 t m ( s ) H 2 ( R 2 ) 2 d s C e k t
for all t > 0 and for some positive constants k and C being independent of t.
Remark 2.
As far as we are concerned, even though there are substantial results regarding the system (1) on R 3 or a bounded domain Ω R n , there are not yet any on R 2 . Theorem 1 and Theorem 2 seem to be the first rigorous theoretical analysis on the initial-value problem (1) and (2) on R 2 and are a first step toward filling this gap.
Remark 3.
Compared to the initial-value problem (1) and (2) on R 3 in [17], establishing an a priori blow-up criterion (4) and a priori upper bounds (5) are nontrivial in the sense that the benefit emanating from H ˙ ( R 3 ) being a Hilbert space will be not granted and that our strategy depends on the weight Hardy inequality over R 2 (see (12)) as well as seeking some new estimates.
The rest of this paper is organized as follows. In Section 2, an a priori blow-up criterion is established. In Section 3, invoking this a priori blow-up criterion, we present that the local strong solution of the Cauchy problem (1) and (2) is global indeed in the two-dimensional setting.

2. Blow-Up Criterion. Proof of Theorem 1

Assume that there exist a time T 0 > 0 and a unique strong solution ( m , p ) to systems (1) and (2) on [ 0 , T 0 ] such that
m C [ 0 , T 0 ] ; H 2 ( R 2 ) , m L 2 [ 0 , T 0 ] ; H 2 ( R 2 ) and p C [ 0 , T 0 ] ; H 2 ( R 2 ) .
Clearly, by employing a standard bootstrap argument, we can extend this strong solution ( m , p ) to the maximal interval of existence [ 0 , T max ) , where either T max = or T max < . If T max < , then we have
lim t T max m ( · , t ) H 2 ( R 2 ) =
and vice versa. The goal of this section is to further establish a more precise blow-up criterion for such strong solution, which allows us to extend the local strong solution to a global one in the next section. To this end, for any T ( 0 , T max ) , we abbreviate
X ( 0 , T ) : = { ( m , p ) | m C [ 0 , T ] ; H 2 ( R 2 ) , m L 2 [ 0 , T ] ; H 2 ( R 2 ) , t m L 2 [ 0 , T ] ; H 1 ( R 2 ) , p L [ 0 , T ] ; H 2 ( R 2 ) }
for simplicity. We now begin with establishing the following a priori estimates to the strong solution of systems (1) and (2) on R 2 .
Lemma 1.
Let ( m , p ) be a solution to system (1) and (2) in X ( 0 , T ) . We have
m ( t ) L 2 2 + m ( t ) L 2 γ 2 γ + p ( t ) L 2 2 + m · p ( t ) L 2 2 + 0 t t m ( s ) L 2 2 d s 1 + m 0 H 2 2 γ + w S L 2 2 ,
m ( t ) L 2 2 + 0 t m ( s ) L 2 2 + m ( s ) L 2 γ 2 γ + p ( s ) L 2 2 + m · p ( s ) L 2 2 d s C T w S L 2 2 + m 0 L 2 2
for any t ( 0 , T ) , where w ( x ) : = ( 1 + | x | ) ( 1 + ln ( 1 + | x | ) ) and the positive constant C depend only on D, c and γ.
Proof. 
Without loss of generality, we may assume that ( m , p ) is sufficiently smooth. The general case can be dealt with by taking an approximation procedure.
We proceed along the lines of the proof of ([17], Lemma 3.1) and have
d d t D 2 m L 2 2 + 1 γ m L 2 γ 2 γ + c 2 p L 2 2 + c 2 m · p L 2 2 + 2 t m L 2 2 = 0 .
For any t ( 0 , T ) , by integrating (9) from 0 to t, we infer that
D 2 m ( t ) L 2 2 + 1 γ m ( t ) L 2 γ 2 γ + c 2 p ( t ) L 2 2 + c 2 m · p ( t ) L 2 2 + 2 0 t t m ( s ) L 2 2 d s = D 2 m 0 L 2 2 + 1 γ m 0 L 2 γ 2 γ + c 2 p 0 L 2 2 + c 2 m 0 · p 0 L 2 2 ,
where the function p 0 = p ( 0 , x ) solves the following Poisson equation
· I + m 0 m 0 p 0 = S , x R 2 .
To estimate the terms related to p 0 on the right-hand side of (10), by using p 0 as a test function in (11), we obtain from the integration by parts that
R 2 | p 0 | 2 d x + R 2 ( m 0 · p 0 ) 2 d x = R 2 p 0 S d x .
Recalling the Hardy inequality over R 2 ([19])
R 2 f w ( x ) 2 d x C R 2 | f | 2 d x for w ( x ) = ( 1 + | x | ) ( 1 + ln ( 1 + | x | ) )
and employing Hölder’s inequality, we arrive at
R 2 p 0 S d x p 0 w L 2 w S L 2 p 0 L 2 w S L 2 1 2 p 0 L 2 2 + C w S L 2 2 .
From this, we conclude that
p 0 L 2 2 + m 0 · p 0 L 2 2 C w S L 2 2 .
This together with (10) yields that
m ( t ) L 2 2 + m ( t ) L 2 γ 2 γ + p ( t ) L 2 2 + m · p ( t ) L 2 2 + 0 t t m ( s ) L 2 2 d s m 0 L 2 2 + m 0 L 2 γ 2 γ + w S L 2 2 1 + m 0 H 2 2 γ + w S L 2 2 ,
where we used Sobolev’s embedding m 0 L 2 γ m 0 H 2 , Young’s inequality and the fact 2 γ 2 in the last inequality.
Next, taking the L 2 inner product of ( ) 1 with m and using the integration by parts, we can obtain
1 2 d d t R 2 | m | 2 d x + D 2 R 2 | m | 2 d x + R 2 | m | 2 γ d x = c 2 R 2 ( m · p ) 2 d x .
Taking the L 2 inner product of ( ) 2 with p, we can obtain from the integration by parts that
R 2 | p | 2 d x + R 2 ( m · p ) 2 d x = R 2 p S d x .
Similarly, we conclude from Hölder’s inequality, Sobolev’s inequality and Hardy’s inequality (12) that
R 2 | p | 2 d x + 2 R 2 ( m · p ) 2 d x C w S L 2 2 ,
which together with (13) yields that
1 2 d d t R 2 | m | 2 d x + R 2 D 2 | m | 2 + | m | 2 γ + c 2 | p | 2 + c 2 ( m · p ) 2 d x C w S L 2 2 .
For any t ( 0 , T ) , we integrate the above inequality from 0 to t and thus obtain
m ( t ) L 2 2 + 0 t D 2 m ( s ) L 2 2 + m ( s ) L 2 γ 2 γ + c 2 p ( s ) L 2 2 + c 2 m · p ( s ) L 2 2 d s C T w S L 2 2 + m 0 L 2 2 .
This completes the proof of Lemma 1.  □
To obtain the higher-order estimates of the solution component m , we should establish the higher-order estimates of the solution component p.
Lemma 2.
Let ( m , p ) be a solution to systems (1) and (2) in X ( 0 , T ) . Then, there exists a positive constant C depending only on D 2 , c 2 and γ such that
2 p L 2 2 + m · p L 2 2 C m L q 2 q q 2 + w S L 2 2 ,
3 p L 2 2 + 2 m · p L 2 2 1 2 3 m L 2 2 + C 2 m L 2 2 m L q 4 q q 2 + w S L 2 4 + C m L q 4 q q 2 + w S L 2 4 + w 2 S L 2 2
for any q ( 2 , + ] and all t ( 0 , T ) .
Proof. 
We begin with establishing the estimate (16). Differentiating ( ) 2 with respect to x k for k = 1 , 2 , we have
Δ k p · k ( m · p ) m + ( m · p ) k m = k S .
By taking the L 2 inner product of (18) with k p , we obtain from the integration by parts that
k p L 2 2 + R 2 k ( m · p ) m · k p + R 2 ( m · p ) k m · k p = R 2 k S k p .
Hence, it reduces to
k p L 2 2 + k m · p L 2 2 = R 2 k ( m · p ) k m · p R 2 ( m · p ) k m · k p + R 2 k S k p .
Invoking Hölder’s inequality and Hardy’s inequality (12), the third term on the right-hand side of (19) can be controlled as follows:
R 2 k S k p d x k p w L 2 w k S L 2 C k p L 2 w k S L 2 .
Based on this, we get from Hölder’s inequality and Young’s inequality that
k p L 2 2 + k m · p L 2 2 k ( m · p ) L 2 k m · p L 2 + ( m · p ) k m L 2 k p L 2 + C k p L 2 w S L 2 1 2 k m · p L 2 2 + k p L 2 2 + k m · p L 2 2 + ( m · p ) k m L 2 2 + C w S L 2 2 .
From this, one arrives at
k p L 2 2 + k m · p L 2 2 2 k m · p L 2 2 + 2 ( m · p ) k m L 2 2 + C w S L 2 2 ,
which implies that
2 p L 2 ( R 2 ) 2 + m · p L 2 ( R 2 ) 2 2 p · m L 2 ( R 2 ) 2 + ( m · p ) m L 2 ( R 2 ) 2 + C w S L 2 2 .
Note that, for any q ( 2 , + ] , Hölder’s inequality, Sobolev’s embedding, Young’s inequality and Lemma 1 yield that
p · m L 2 ( R 2 ) 2 m L q 2 p L 2 q q 2 2 m L q 2 p L 2 2 1 2 q 2 p L 2 4 q m L q 2 2 p L 2 4 q 1 4 2 p L 2 2 + C m L q 2 q q 2 ,
and similarly that
( m · p ) m L 2 ( R 2 ) 2 1 4 ( m · p ) L 2 2 + C m L q 2 q q 2 .
Hence, we can conclude from (20) that, for any q ( 2 , + ] ,
2 p L 2 ( R 2 ) 2 + m · p L 2 ( R 2 ) 2 C m L q 2 q q 2 + w S L 2 2 .
This indicates that (16) holds.
We now turn to the higher-order estimate (17). To achieve this, differentiating (18) with respect to x j for j = 1 , 2 , we obtain
Δ k j p · ( m · p ) k j m + j ( m · p ) k m + k ( m · p ) j m + k j ( m · p ) m = k j S .
Then, multiplying (22) by i j 2 p and using the integration by parts, we have
k j p L 2 2 + ( m · p ) k j m · k j p + j ( m · p ) k m + k ( m · p ) j m · k j p + k j m · p m · k j p = k j S k j p .
Since
m · k j p = k j ( m · p ) k j m · p j m · k p k m · j p ,
we have
k j p L 2 2 + k j m · p L 2 2 = ( m · p ) k j m · k j p j ( m · p ) k m + k ( m · p ) j m · k j p + k j m · p + j m · k p + k m · j p k j ( m · p ) + k j S k j p .
We now use Hölder’s inequality and Hardy’s inequality (12) to estimate the rightmost term of the above equation as follows:
R 2 k j S k j p k j p w L 2 w k j S L 2 C k j p L 2 w k j S L 2 .
One has
k j p L 2 2 + k j m · p L 2 2 k j p L 2 ( m · p ) k j m L 2 + j ( m · p ) k m L 2 + k ( m · p ) j m L 2 + k j m · p L 2 k j m · p L 2 + j m · k p L 2 + k m · j p L 2 + C k j p L 2 w 2 S L 2 .
By employing Young’s inequality, we can deduce that
k j p L 2 2 + k j m · p L 2 2 ( m · p ) k j m L 2 2 + j ( m · p ) k m L 2 2 + k ( m · p ) j m L 2 2 + k j m · p L 2 2 + j m · k p L 2 2 + k m · j p L 2 2 + w 2 S L 2 2 ,
and thus that
3 p L 2 2 + 2 m · p L 2 2 ( m · p ) 2 m L 2 2 + ( m · p ) · m L 2 2 + j = 1 2 p · j m L 2 2 + m 2 p L 2 2 + w 2 S L 2 2 .
We estimate the terms on the right-hand side of (23) one by one. Firstly, we conclude from Hölder’s inequality, Sobolev’s embedding, the interpolation, Young’s inequality and Lemma 1 that
( m · p ) 2 m L 2 2 ( m · p ) L 6 2 2 m L 3 2 m · p L 2 2 3 ( m · p ) L 2 4 3 2 m L 2 4 3 3 m L 2 2 3 1 8 3 m L 2 2 + C 2 m L 2 2 ( m · p ) L 2 2 .
Secondly, we get from a similar procedure that, for any q ( n , + ] ,
( m · p ) · m L 2 2 ( m · p ) L 2 q q 2 2 m L q 2 ( m · p ) L 2 2 1 2 q 2 ( m · p ) L 2 4 q m L q 2 1 2 2 ( m · p ) L 2 2 + C ( m · p ) L 2 2 m L q 2 q q 2 .
Similarly, we have
p · j m L 2 2 1 8 3 m L 2 2 + C 2 m L 2 2 2 p L 2 2 , m 2 p L 2 2 1 2 3 p L 2 2 + C 2 p L 2 2 m L q 2 q q 2
for any q ( 2 , + ] . It follows from inserting the above estimates into (23) that
3 p L 2 2 + 2 m · p L 2 2 1 2 3 m L 2 2 + C 2 m L 2 2 2 p L 2 2 + ( m · p ) L 2 2 + C 2 p L 2 2 + ( m · p ) L 2 2 m L q 2 q q 2 + C w 2 S L 2 2
for any q ( 2 , + ] . This implies (17) and thus completes the proof of Lemma 2.  □
With the improved a priori estimates on p at hand, we can directly derive the higher-order estimates of the solution component m under the same assumptions as above.
Lemma 3.
Let ( m , p ) be a solution to systems (1) and (2) in X ( 0 , T ) . Then, there exists a positive constant C depending only on D, c and γ such that for all 0 < t < T ,
d d t 2 m L 2 2 + D 2 2 3 m L 2 2 2 m L 2 2 m L q 4 q q 2 + w S L 2 4 + m L q 4 q q 2 + m L 2 8 ( γ 1 ) 9 2 γ + w S L 2 4 + w 2 S L 2 2 + 1
for 1 γ < 3 , and
d d t 2 m L 2 2 + D 2 2 3 m L 2 2 2 m L 2 2 m L q 4 q q 2 + w S L 2 4 + m L q 4 q q 2 + w S L 2 4 + w 2 S L 2 2 + 1 + m L 2 1 + 1 2 q 3 m L q 2 q ( γ 2 ) 2 q 3 2 m L 2
for γ 3 .
Proof. 
Similar to ([17], Lemma 3.3), we get that there exists a positive constant C 0 independent of t such that
d d t 2 m L 2 2 + D 2 3 m L 2 2 C 0 ( m · p ) p L 2 2 + ( | m | 2 ( γ 1 ) m ) L 2 2 .
We estimate the terms on the right-hand side of (24) one by one. For the first term on the right-hand side of (24), we obtain from Hölder’s inequality, the interpolation inequality, Young’s inequality and Lemma 1 that
( m · p ) p L 2 2 ( m · p ) p L 2 2 + ( m · p ) 2 p L 2 2 ( m · p ) L 3 2 p L 6 2 + m · p L 6 2 2 p L 3 2 ( m · p ) L 2 4 3 2 ( m · p ) L 2 2 3 p L 2 2 3 2 p L 2 4 3 + m · p L 2 2 3 ( m · p ) L 2 4 3 2 p L 2 4 3 3 p L 2 2 3 D 2 2 C 0 2 ( m · p ) L 2 2 + 3 p L 2 2 + C ( m · p ) L 2 2 2 p L 2 2 + ( m · p ) L 2 2 2 p L 2 2 ,
which together with Lemma 2 yields that
( m · p ) p L 2 2 D 2 4 C 0 3 m L 2 2 + C 2 m L 2 2 m L q 4 q q 2 + w S L 2 4 + C m L q 4 q q 2 + w S L 2 4 + w 2 S L 2 2 + 1 .
To estimate the second term on the right-hand side of (24), by using the Leibniz’s product rule and Hölder’s inequality, we obtain
( | m | 2 ( γ 1 ) m ) L 2 2 | m | 2 ( γ 1 ) m L 2 2 m L 4 ( γ 1 ) m L 2 2 .
For 1 γ < 3 , a straightforward application of the interpolation inequality yields that
m L 4 ( γ 1 ) m L 6 24 ( γ 1 ) 7 3 m L 2 4 ( γ 1 ) 7 m L 2 8 ( γ 1 ) 7 m L 2 16 ( γ 1 ) 7 3 m L 2 4 ( γ 1 ) 7
and thus that
( | m | 2 ( γ 1 ) m ) L 2 2 m L 2 8 ( γ 1 ) 7 m L 2 2 + 16 ( γ 1 ) 7 3 m L 2 4 ( γ 1 ) 7 D 2 4 C 0 3 m L 2 2 + C m L 2 8 ( γ 1 ) 9 2 γ m L 2 14 + 16 ( γ 1 ) 9 2 γ .
On the other hand, for γ 3 , we have from Hölder’s inequality and the interpolation inequality that
m L 4 ( γ 1 ) m L 4 m L 4 ( γ 2 ) m L 6 3 2 m L 2 m L 6 3 ( q 2 ) 2 ( 2 q 3 ) m L q q 2 ( 2 q 3 ) 4 ( γ 2 ) m L 2 m L 2 2 2 m L 2 m L 2 q 2 2 ( 2 q 3 ) m L 2 q 2 2 q 3 m L q q 2 ( 2 q 3 ) 4 ( γ 2 ) m L 2 1 + 1 2 q 3 m L 2 4 ( q 2 ) ( γ 2 ) 2 q 3 + 2 m L q 2 q ( γ 2 ) 2 q 3 2 m L 2
for any q ( 2 , + ] and thus
( | m | 2 ( γ 1 ) m ) L 2 2 m L 2 1 + 1 2 q 3 m L 2 4 ( q 2 ) ( γ 2 ) 2 q 3 + 4 m L q 2 q ( γ 2 ) 2 q 3 2 m L 2 .
Substituting (25), (26) and (27) into (24), we deduce from Lemma 1 that, for 1 γ < 3 ,
d d t 2 m L 2 2 + D 2 2 3 m L 2 2 2 m L 2 2 m L q 4 q q 2 + w S L 2 4 + m L q 4 q 2 q 3 + m L 2 8 ( γ 1 ) 9 2 γ + w S L 2 4 + w 2 S L 2 2 + 1 ,
and that, for γ 3 ,
d d t 2 m L 2 2 + D 2 2 3 m L 2 2 2 m L 2 2 m L q 4 q q 2 + w S L 2 4 + m L q 4 q q 2 + w S L 2 4 + w 2 S L 2 2 + 1 + m L 2 1 + 1 2 q 3 m L q 2 q ( γ 2 ) 2 q 3 2 m L 2 .
This completes the proof of Lemma 3.  □
Invoking Lemmas 1–3, we now establish a more precise blow-up criterion for the strong solution as follows.
Lemma 4.
Let ( m , p ) be a solution to systems (1) and (2) in X ( 0 , T ) and T max be the maximal existence time. If T max < , it holds that
lim sup t T max m ( t ) H 2 ( R 2 ) =
if and only if
0 T max m ( t ) L q ( R 2 ) β d t =
for any q ( 2 , + ] , where
β = max 4 q q 2 , 4 q ( γ 2 ) 2 q 3 .
Proof. 
We can first deduce from Lemma 3 that, for 1 γ < 3 and all 0 < t < T max ,
d d t 2 m L 2 2 + D 2 2 3 m L 2 2 2 m L 2 2 m L q 4 q q 2 + 1 + m L q 4 q q 2 + m L 2 8 ( γ 1 ) 9 2 γ + 1 ,
and that for γ 3 and all 0 < t < T max ,
d d t 2 m L 2 2 + D 2 2 3 m L 2 2 2 m L 2 2 m L q 4 q q 2 + 1 + m L q 4 q q 2 + m L q 4 q ( γ 2 ) 2 q 3 m L 2 4 ( q 1 ) 2 q 3 + 1 .
Hence, Gronwall’s inequality implies that, for 1 γ < 3 ,
2 m ( t ) L 2 2 e C 0 t m ( s ) L q 4 q q 2 + 1 d s 2 m 0 L 2 2 + C 0 t m ( s ) L q 4 q q 2 + m ( s ) L 2 8 ( γ 1 ) 9 2 γ + 1 d s ,
and that, for γ 3 ,
2 m ( t ) L 2 2 e C 0 t m ( s ) L q 4 q q 2 + 1 d s 2 m 0 L 2 2 + C 0 t m L q 4 q q 2 + m ( s ) L q 4 q ( γ 2 ) 2 q 3 m ( s ) L 2 4 ( q 1 ) 2 q 3 + 1 d s .
By employing Lemma 1 again, we conclude that
lim sup t T max m ( t ) H 2 =
if and only if
0 T max m ( s ) L q 4 q q 2 d s =
for 1 γ < 3 , and
0 T max m ( s ) L q 4 q q 2 d s + 0 T max m ( s ) L q 4 q ( γ 2 ) 2 q 3 d s =
for γ 3 , where q ( 2 , + ] . This completes the proof of Lemma 4.  □
Based on Lemma 4, we now prove Theorem 1.
Proof of Theorem 1.
The blow-up criterion is a by-product of Lemma 4. Thus, we have completed the proof of Theorem 1.  □

3. Global Existence. Proof of Theorem 2

In this section, invoking blow-up criterion (4), we show that the local strong solution of the Cauchy problem (1) and (2) is global indeed in the two-dimensional setting.
Lemma 5.
Assume that all conditions in Theorem 2 hold and the systems (1) and (2) possess a solution ( m , p ) in X ( 0 , T ) . For fixed r ( 2 , ) , then there exists a positive constant h ( r ) (see (32) below) such that, for all D h ( r ) ,
m ( t ) L r r C 1 + t
for any 0 t T and certain positive constant C independent of t.
Proof. 
We begin with r 4 . We now apply the operator i ( i = 1 , 2 ) to ( ) 1 :
t i m D 2 Δ i m + i | m | 2 ( γ 1 ) m = c 2 i ( m · p ) p .
Then, taking the L 2 inner product of (29) with i m | i m | r 2 and using the integration by parts, we obtain
1 r d d t R 2 | i m | r d x + D 2 R 2 i m · i m | i m | r 2 d x + R 2 i | m | 2 ( γ 1 ) m · i m | i m | r 2 d x = R 2 i m · p p · i m | i m | r 2 d x + R 2 m · p i p · i m | i m | r 2 d x .
The second and third terms on the left-hand side of (30) can be calculated as follows:
R 2 i m · i m | i m | r 2 d x = R 2 | i m | 2 | i m | r 2 d x + ( r 2 ) R 2 i m · i m 2 | i m | r 4 d x 4 ( r 1 ) r 2 R 2 | i m | r 2 2 d x ,
and
R 2 i | m | 2 ( γ 1 ) m · i m | i m | r 2 d x = R 2 | m | 2 ( γ 1 ) | i m | r d x + 2 ( γ 1 ) R 2 m · i m 2 | m | 2 ( γ 2 ) | i m | r 2 d x .
Due to r 2 > 0 and γ 1 0 , we then get from Hölder’s inequality and (30) that
1 r d d t i m L r r + 4 ( r 1 ) D 2 r 2 | i m | r 2 L 2 2 i m · p L 2 p · i m | i m | r 2 L 2 + i p L 2 m · p i m | i m | r 2 L 2 i m · p L 2 p L 3 i m L 6 ( r 1 ) r 1 + i p L 2 m · p L 3 i m L 6 ( r 1 ) r 1 .
We use the interpolation inequality and Lemma 1 to show
1 r d d t i m L r r + 4 ( r 1 ) D 2 r 2 | i m | r 2 L 2 2 i m · p L 2 p L 2 2 3 2 p L 2 1 3 i m L 6 ( r 1 ) r 1 + i p L 2 m · p L 2 2 3 m · p L 2 1 3 i m L 6 ( r 1 ) r 1 i m · p L 2 2 p L 2 1 3 i m L 6 ( r 1 ) r 1 + i p L 2 m · p L 2 1 3 i m L 6 ( r 1 ) r 1 .
Hence, by employing (16) with q = 4 , it arrives at
1 r d d t i m L r r + 4 D 2 ( r 1 ) r 2 | i m | r 2 L 2 2 i m L 6 ( r 1 ) r 1 m L 4 4 + 1 2 3 i m L 6 ( r 1 ) r 1 m L 4 8 3 + i m L 6 ( r 1 ) r 1 .
On the other hand, a direct application of the interpolation inequality yields that for r 4
m L r r 2 = | m | r 2 L 2 | m | r 2 L 1 1 2 | m | r 2 L 2 1 2 m L r 2 r 4 | m | r 2 L 2 1 2 m L 2 r 2 ( r 2 ) m L r r ( r 4 ) 4 ( r 2 ) | m | r 2 L 2 1 2 ,
which leads to
m L r r 2 4 ( r 2 ) m L 2 r 2 ( r 2 ) | m | r 2 L 2 1 2 .
Here, we used the fact that r 2 2 , due to r 4 . From this, one then has
m L 4 m L 2 r 4 2 ( r 2 ) m L r r 2 ( r 2 ) m L 2 r 4 2 ( r 2 ) + 1 r 2 | m | r 2 L 2 1 r .
Similarly, we also have
m L 6 ( r 1 ) r 2 = | m | r 2 L 12 ( r 1 ) r | m | r 2 L 2 r 6 ( r 1 ) | m | r 2 L 2 5 r 6 6 ( r 1 ) m L r r 2 12 ( r 1 ) | m | r 2 L 2 5 r 6 6 ( r 1 ) m L 2 r 6 ( r 1 ) | m | r 2 L 2 r 2 6 ( r 1 ) + 5 r 6 6 ( r 1 ) .
Here, we used the fact that 12 ( r 1 ) r > 2 due to r 4 . Whereupon, invoking the above inequalities, we infer from (31) and Lemma 1 that
1 r d d t m L r r + 4 D 2 ( r 1 ) r 2 | m | r 2 L 2 2 m L 2 1 3 | m | r 2 L 2 2 ( 3 r 4 ) 3 r m L 2 4 3 | m | r 2 L 2 8 3 r + m L 2 1 3 | m | r 2 L 2 2 ( 3 r 4 ) 3 r | m | r 2 L 2 2 ( 3 r 4 ) 3 r + 8 3 r + | m | r 2 L 2 2 ( 3 r 4 ) 3 r .
Note that 0 < 2 ( 3 r 4 ) 3 r < 2 and 2 ( 3 r 4 ) 3 r + 8 3 r = 2 , due to r 4 . Thus, Young’s inequality implies from the above inequality that
1 r d d t m L r r + 4 D 2 ( r 1 ) r 2 | m | r 2 L 2 2 C 1 | m | r 2 L 2 2 + C 2
for two positive constants C 1 , C 2 independent of t. Furthermore, taking
h ( r ) = C 1 r 2 4 ( r 1 ) 1 2 ,
under our assumptions, then one arrives at
1 r d d t m L r r C 2 .
Finally, a time integration of the above inequality shows that, for all r 4 ,
m ( t ) L r r m 0 L r r + C 2 t .
Next, we turn to the case of r ( 2 , 4 ) . Recalling the interpolation inequality
m L r m L 2 4 r r m L 4 2 ( r 2 ) r for all r ( 2 , 4 ) ,
we can get from (7), (33) and Young’s inequality that
m ( t ) L r r m 0 L 4 4 + C 2 t 2 ( r 2 ) 4 1 + t .
Combining with (33) and (34) yields (28). This finishes the proof of Lemma 5.  □
By means of Lemma 5 and blow-up criterion (4), we can proceed to prove Theorem 2.
Proof of Theorem 2.
By taking r = 4 in (32), we can get from Lemma 5 that
m ( t , · ) L 4 4 C ( 1 + t ) for all t > 0 .
Based on this, the local strong solution ( m , p ) can be extended to a global one by employing the blow-up criterion (4) obtained in Theorem 2. Moreover, the boundedness (5) directly follows from Lemma 1, Lemma 2, Lemma 3 and Lemma 5. This completes the proof of Theorem 2.  □

4. Discussion

In this section, let us begin with presenting a short summary of our main results for the nonlinear PDE system (1) and (2) in the two-dimensional setting:
  • By developing some new dissipation mechanisms hidden in the system (1), we provide an a priori blow-up criterion of a strong solution of the Cauchy problems (1) and (2);
  • Invoking an a priori blow-up criterion and seeking some suitable estimates, we establish a priori upper bounds to a strong solution for all positive times.
We conclude this section by providing two unsolved problems which will be the subject of our future research:
  • How to establish the local existence of a strong solution of the Cauchy problems (1) and (2);
  • We establish an a priori blow-up criterion and a priori upper bounds to a strong solution of systems (1) and (2) when γ 1 . Can these be generalized to γ > 1 2 ?

Author Contributions

Conceptualization, J.S. and B.L.; Methodology, J.S. and B.L.; Validation, J.S.; Writing—original draft, B.L.; and Writing—review and editing, J.S.

Funding

This research received no external funding.

Acknowledgments

The authors are very grateful to the anonymous referees for their valuable comments and suggestions, which greatly improved the manuscript, and to editor for his/her kind help.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nancy, D.; Julie, K. Vascular patterning and leaf shape. Curr. Opin. Plant Biol. 2001, 4, 50–56. [Google Scholar]
  2. Malinowski, R. Understanding of leaf development-the science of complexity. Plants 2013, 2, 396–415. [Google Scholar] [CrossRef] [PubMed]
  3. Sedmera, D. Function and form in the developing cardiovascular system. Cardiovasc. Res. 2011, 91, 252–259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Anne, E.; Ferdinand, L.N.; Monica, A.; Peter, C. Guidance of vascular and neural network formation. Curr. Opin. Neurobiol. 2005, 15, 108–115. [Google Scholar]
  5. Michel, O.; Biondi, J. Morphogenesis of neural networks. Neural Process. Lett. 1995, 2, 9–12. [Google Scholar] [CrossRef]
  6. Hu, D.; Cai, D. Adaptation and optimization of biological transport networks. Phys. Rev. Lett. 2013, 111, 138701. [Google Scholar] [CrossRef] [PubMed]
  7. Albi, G.; Artina, M.; Foransier, M.; Markowich, P. Biological transportation networks: Modeling and simulation. Anal. Appl. 2016, 14, 185–206. [Google Scholar] [CrossRef]
  8. Albi, G.; Burger, M.; Haskovec, J.; Markowich, P.; Schlottbom, M. Continuum Modeling of Biological Network Formation. In Active Particles, Volume 1. Modeling and Simulation in Science and Technology; Bellomo, N., Degond, P., Tamdor, T., Eds.; Birkhäuser: Boston, MA, USA, 2017; pp. 1–48. [Google Scholar]
  9. Hu, D. Optimization, Adaptation, and Initialization of Biological Transport Networks. Notes Lect. 2014. [Google Scholar]
  10. Haskovec, J.; Markowich, P.; Perthame, B. Mathematical analysis of a PDE system for biological network formation. Commun. Part. Differ. Equ. 2015, 40, 918–956. [Google Scholar] [CrossRef]
  11. Haskovec, J.; Markowich, P.; Perthame, B.; Schlottbom, M. Notes on a PDE system for biological network formation. Nonlinear Anal. 2016, 138, 127–155. [Google Scholar] [CrossRef] [Green Version]
  12. Li, B. Long time behavior of the solution to a parabolic-elliptic system. Comput. Math. Appl. 2019. [Google Scholar] [CrossRef]
  13. Liu, J.; Xu, X. Partial regularity of weak solutions to a PDE system with cubic nonlinearity. J. Differ. Equ. 2018, 264, 5489–5526. [Google Scholar] [CrossRef] [Green Version]
  14. Xu, X. Regularity theorems for a biological network formulation model in two space dimensions. Kinet. Relat. Models 2018, 11, 397–408. [Google Scholar] [CrossRef] [Green Version]
  15. Xu, X. Life-span of smooth solutions to a PDE system with cubic nonlinearity. arXiv 2017, arXiv:1706.06057. [Google Scholar]
  16. Katifori, E.; Szöllosi, G.J.; Magnasco, M.O. Damage and fluctuations induce loops in optimal transport networks. Phys. Rev. Lett. 2010, 104, 048704. [Google Scholar] [CrossRef] [PubMed]
  17. Li, B. On the blow-up criterion and global existence of a nonlinear PDE system in biological transport networks. Kinet. Relat. Models 2019, in press. [Google Scholar]
  18. Serrin, J. On the interior regularity of weak solutions of Navier-Stokes equations. Arch. Ration. Mech. Anal. 1962, 9, 187–195. [Google Scholar] [CrossRef]
  19. Duan, R.; Lorz, A.; Markowich, P. Global solutions to the coupled chemotaxis-fluid equations. Commun. Part Differ. Equ. 2010, 35, 1635–1673. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Shen, J.; Li, B. A Priori Estimates for a Nonlinear System with Some Essential Symmetrical Structures. Symmetry 2019, 11, 852. https://doi.org/10.3390/sym11070852

AMA Style

Shen J, Li B. A Priori Estimates for a Nonlinear System with Some Essential Symmetrical Structures. Symmetry. 2019; 11(7):852. https://doi.org/10.3390/sym11070852

Chicago/Turabian Style

Shen, Jieqiong, and Bin Li. 2019. "A Priori Estimates for a Nonlinear System with Some Essential Symmetrical Structures" Symmetry 11, no. 7: 852. https://doi.org/10.3390/sym11070852

APA Style

Shen, J., & Li, B. (2019). A Priori Estimates for a Nonlinear System with Some Essential Symmetrical Structures. Symmetry, 11(7), 852. https://doi.org/10.3390/sym11070852

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop