1. Introduction
Ordinary and partial differential equations play a key role in modelling for all sciences: engineering, physics, chemistry, biology, medicine, economics and many others. The right understanding of the behaviour of solutions (in particular, well-posedness versus blow-up) means not only to predict the future of trajectories but also to establish strategies for control (i.e., optimisation). Concerning PDE and economics, it is interesting to cite the nice survey [
1] and the references therein on many different problems dealing with effects such as aggregation and repulsion, optimal control, mean-field games and so on as applications.
Parabolic PDE models reflect the diffusion phenomena due to local touching of molecules and dissipation of energy, and when different internal and external factors come into play, they link naturally to some reaction-diffusion models, such as the growth versus capacity of the environment in biology or the endogenous growth versus the neoclassical theories in economics. In particular, capital accumulation distribution in space and time following spatial extensions of the continuous Ramsey model [
2] by Brito [
3,
4,
5] and others later use the semilinear parabolic PDE
This spatiality introduces important issues about the steady states distribution and the dynamic evolution, convergence, local interaction among local agents and so on.
Not for the sake of generality but for real modelling purposes, in the last two decades the increment of nonlocal PDE models that attempt to capture in a more accurate way the real spreading of the problem (density of population, capital accumulation, consumption or prices and innovation indexes and so on) has been very important. Firstly we might comment about extensions by using some nonlocal operators acting in the right-hand side of the PDE and/or the boundary conditions as integral operators, leading to integro-differential equations. Among others, we can cite [
6] for a system coupling capital and the pollution stock model; a population dynamic model in [
7]
the elliptic (stationary) counterpart in population/physics models of the Fischer-KPP [
8]; and a logistic model [
9]. Secondly, we wish to point out that the nonlocal extensions have also been performed on the diffusion operators as well. The literature about fractional Laplacian is vast nowadays. However, let us concentrate in an intermediate step. Coming originally from modelling of bacterial populations in biology, the introduction of a nonlocal viscosity in front of the Laplacian has become an interesting problem for different applications and for its mathematical study, as for example occurs in the equation
In this way, the spreading (or aggregating/concentrating) effects are given by the increasing (resp. non-increasing) function
a as a viscosity nonlocal coefficient. One should cite Prof. Chipot and his collaborators [
10,
11,
12,
13,
14,
15,
16] among others for a detailed analysis, including existence, uniqueness, steady states and convergence of evolutionary solutions to equilibria.
When the reaction term
f depends on the unknown
u(here the functional
may represent a general nonlocal functional acting over the whole domain
for instance,
or
), equilibria are difficult to analyse. Oppositely to ordinary differential equations, the analysis of the existence of stationary states for the above problem is much more involved. Additionally, comparing the reaction-diffusion equations with local diffusion, another difficulty is that in general a Lyapunov functional is not known to exist in most cases.
The dynamical analysis of problem (
1) and in particular the existence of global attractors has been established till now in several papers (cf. [
17,
18,
19,
20,
21]). Other differential operators such as the
p-Laplacian coupled with nonlocal viscosity has also been considered (cf. [
21,
22,
23]). However, in general little is known about the internal structure of the attractor, which is very important as it gives us a deep insight into the long-term dynamics of the problem. When we manage to obtain a Lyapunov functional some insights can be obtained.
If we consider the non-local equation
with Dirichlet boundary conditions, then it is possible to define a suitable Lypaunov functional. In [
18] it is shown that regular and strong solutions generate (possibly) multivalued semiflows having a global attractor which is described by the unstable set of the stationary points. Although this is already a good piece of information, our goal is to describe the structure of the attractor as accurately as possible. For this aim we need to study the particular situation where the domain is one-dimensional and the function
f is of the type of the standard Chafee-Infante problem, for which the dynamics inside the attractor has been completely understood [
24].
The first step when studying the structure of the attractor consists of analysing the stationary points. In the case where the function
f is odd and Equation (
2) generates a continuous semigroup, the existence of fixed points of the type given in the Chafee-Infante problem was established in [
25]. Moreover, if
a is non-decreasing, then they coincide with the ones in the Chafee-Infante problem, and moreover, in [
26] the stability and hyperbolicity of the fixed points was studied. In this paper we extend these results for a more general function
f (not necessarily odd and for which we do not known whether the Cauchy problem has a unique solution or not), showing that Equation (
2) undergoes the same cascade of bifurcations as the Chafee-Infante equation. Moreover, when we allow the function
a to decrease, though the problem possesses at least the same fixed point as in the Chafee-Infante problem, we show that more equilibria can appear. For a non-decreasing function
a and an odd function
f we prove also that even when uniqueness fails, the stability of the fixed points is the same as for the corresponding ones in the Chafee-Infante problem. Finally, we are able to prove that in this last case we have a dynamically gradient semiflow with respect to the disjoint family of isolated weakly invariant sets generated by the equilibria, which is ordered by the number of zeros of the fixed points. More precisely, the attractor consists of the set of equilibria and their heteroclinic connections and a connection from a fixed point to another is allowed only if the number of zeros of the first one is greater.
In
Section 3 we study the existence of strong solutions of the Cauchy problem in the space
. In
Section 4 we prove that strong solutions generate a multivalued semiflow in
having a global attractor which is equal to the unstable set of the stationary points. In
Section 5 we study the existence and properties of equilibria. In
Section 6 we analyse the stability of the fixed points and establish that the semiflow is dynamically gradient.
2. Setting of the Problem
Let us consider the following problem:
where
and
. Throughout the paper we will use the following conditions (but not all of them at the same time):
- (A1)
.
- (A2)
- (A3)
exists and .
- (A4)
f is strictly concave if and strictly convex if .
- (A5)
Growth and dissipation conditions: for
, we have
- (A6)
The function
satisfies:
- (A7)
The function
satisfies:
where
- (A8)
The function is non-decreasing.
- (A9)
- (A10)
h does not depend on time and
We define the function
. We observe that from (
4) we have
whereas (
5) implies
Additionally, from condition (
6) it follows that for all
, there exists a constant
such that
for all
. Hence, there exists
such that
In addition, it follows that
where
. These two inequaities are also true under condition (
5).
The main aim of this paper consists of describing in as much detail as possible the internal structure of the global attractor in a similar way as for the classical Chafee-Infante equation.
Some of these conditions will be used all the time, whereas other ones will be used only in certain results. In particular, the function
h will be considered as a time-dependent function satisfying (A9) only for establishing the existence of solution for problem (
3). However, since we will study the asymptotic behaviour of solutions in the autonomous situation, for the second part concerning the existence and properties of global attractors, the function
h will be time-independent, so assumption (A10) will be used instead. Finally, in order to study the structure of the global attractors in terms of the stationary points and their possible heteroclinic connections, we will assume that
Throughout the paper, will denote the norm in the Banach space
4. The Existence and Structure of Attractors
In this section, we will prove the existence of a global attractor for the semiflow generated by strong solutions in the autonomous case. Thus, the function
h will be an independent of time function satisfying (A10) instead of (A9). Additionally, we will establish that the attractor is equal to the unstable set of the stationary points (see the definition in (
45)).
Throughout this section, for a metric space
X with metric
d we will denote by
the Hausdorff semidistance from
C to
D, that is,
Let us consider the phase space
and the sets
Denote by
the class of nonempty subsets of
X. We define the (possibly multivalued) map
by
In order to study the map G let us consider the following axiomatic properties of the set :
- (K1)
For every there is satisfying .
- (K2)
for every and (translation property).
- (K3)
Let
be such that
for some
. Then, the function
defined by
belongs to
(concatenation property).
- (K4)
For every sequence satisfying in X, there is a subsequence and such that for every .
Assuming conditions (A1), (A6), (A10) and (
13), property (K1) follows from Theorem 1, whereas (K2) and (K3) can be proved easily using equality (
12). By ([
30] Proposition 2) or ([
31] Lemma 9) we know that
fulfilling (K1) and (K2) gives rise to a multivalued semiflow
G through (
28) (m-semiflow for short), which means that:
for all
for all and .
Moreover, (K3) implies that the m-semiflow is strict, that is, for all and .
We will show first that the m-semiflow G possesses a bounded absorbing set in the space and that property (K4) is satisfied.
Lemma 1. Assume conditions (A1), (A6), (A10) and (13). Given weakly in , there exists a subsequence of (relabelled the same) and such thatAdditionally, if strongly in , then for we get strongly in . Proof. Since
and
, we have by ([
27] p. 102) that
and
. Additionally, as
, by regularization one can show that
is an absolutely continuous function on
and
By a similar argument as in Theorem 1, there is a subsequence of
such that
Therefore, arguing as in the proof of Theorem 1, there exist
and a subsequence
, relabelled the same, such that
We also need to prove that
in
for all
. To that end, we multiply (
3) by
, and using (A10), (
29) and (
31) we have
Thus, we obtain
Since this inequality is also true for
, the functions
are continuous and non-increasing in
. Moreover, from (
32) we deduce that
Take
and
such that
and
for all
j. Then
For any
there exist
and
such that
if
. Then
so
, which follows by contradiction using the continuity of the function
. As
weakly in
implies that
, we obtain
so that
strongly in
.
Finally, if
strongly in
and we take
, then
so
. Repeating the above argument, we infer that
strongly in
. □
Corollary 1. Assume the conditions of Lemma 1. Then the set satisfies condition .
The map is said to be upper semicontinuous if for every and for an arbitrary neighbourhood in X there is such that as soon as , we have .
Proposition 1. Assume the conditions of Lemma 1. The multivalued semiflow G is upper semicontinuous for all . Additionally, it has compact values.
Proof. By contradiction let us assume that there exist , a neighbourhood and sequences such that , converges strongly to in and for all . Thus, there exists such that . From Lemma 1 there exists a subsequence of which converges to some . This contradicts for any . □
In order to prove the existence of an absorbing set in the space
we need to use the stronger condition (A5) instead of (
13).
Proposition 2. Assume that conditions (A1), (A5), (A6) and (A10) hold. Then the m-semiflow G has a bounded absorbing set in ; that is, there exists a constant such that for any there is a time such thatwhere Moreover, there is such that Proof. By multiplying Equation (
3) by
u and using (A6) and (
9), we get
By using the Poincaré inequality it follows that
where
,
. We take a small enough
so that
. Then Gronwall’s lemma gives
Hence, taking
we get (
33) for
On the other hand, using again the Poincaré inequality from (
35) we get
and integrating from
t to
we obtain
Therefore, applying (
33) and (
34) follows. □
Further, in order to obtain an absorbing set in we need to assume additionally that either the function is bounded above or that it is non-decreasing.
Proposition 3. Assume the conditions in Proposition 2 and that either (A7) or (A8) holds true. Then there exists an absorbing set for G, which is compact in .
Proof. In view of Proposition 2 we have an absorbing set in . Let be such that for all .
Through multiplying (
3) by
u and using (
9) and (
36) we get
Thus, integrating between
t and
,
, we deduce by using (
36) again that
Additionally, if
in (A5), we multiply again by (
3) by
u and use (
5) and (A6) to obtain
Integrating over
we have
If we assume (A7), by (
37) and (A6) we have that
If we assume (A8), by (
37) we obtain
On the other hand, by (
7) we get
By using (
29) and (
30) we can argue as in Theorem 1 to obtain
Since (
38)–(
41) imply that
we can apply the uniform Gronwall lemma to get
so by condition
, (
10) and (
36) it follows that
for all
. In particular,
for any strong solution
with initial condition
.
For any
with
and any
such that
, the semiflow property
and
, if
imply that
Then there exists
such that the closed ball
in
centred at 0 with radius
M is absorbing for
G.
By Lemma 1 the set is an absorbing set which is compact in . □
Given an m-semiflow a set is said to be negatively (positively) invariant if () for all , and strictly invariant (or, simply, invariant) if it is both negatively and positively invariant.
We recall that a set
is called a global attractor for the m-semiflow
G if it is negatively invariant and attracts all bounded subsets; i.e.,
as
. When
is compact, it is the minimal closed attracting set ([
32] Remark 5).
Theorem 2. Assume the conditions of Proposition 3. Then the multivalued semiflow G possesses a global compact invariant attractor .
Proof. From Propositions 1 and 3 we deduce that the multivalued semiflow
G is upper semicontinuous with closed values and the existence of an absorbing which is compact in
. Therefore, by ([
32] Theorem 4 and Remark 8) the existence of the global invariant attractor and its compactness in
follow. □
We recall some concepts which are necessary to study the structure of the global attractor.
Definition 3. A map is a complete trajectory of if for all . It is a complete trajectory of G if for every ,
An element is a fixed point of if . We denote the set of all fixed points by .
An element is a fixed point of G if for every .
Several properties concerning fixed points, complete trajectories and global attractors are summarised in the following results [
33].
Lemma 2. Let (K1)-(K2) hold. Then each fixed point (complete trajectory) of is also a fixed point (complete trajectory) of G.
Let (K1)-(K4) hold. Then the fixed points of and G are the same. In addition, a map is a complete trajectory of if and only if it is continuous and a complete trajectory of G.
The standard well-known result in the single-valued case for describing the attractor as the union of bounded complete trajectories reads in the multivalued case as follows.
Theorem 3. Suppose that (K1) and (K2) are satisfied and that either (K3) or (K4) holds true. The semiflow G is assumed to have a compact global attractor . Then where stands for the set of all bounded complete trajectories in .
In view of Theorem 3, as
satisfies (K3) and (K4) (by Corollary 1), the global attractor is characterised in terms of bounded complete trajectories, so (
42) follows.
The set
B is said to be weakly invariant if for any
there exists a complete trajectory
of
contained in
B such that
. Characterisation (
42) implies that the attractor
is weakly invariant.
The set of fixed points is characterised as follows.
Lemma 3. Assume the conditions of Lemma 1. Let be the set of such thatThen . Proof. If
then
Thus,
satisfies (
12) and
in
, so (
43) is satisfied. Let
. Then the map
satisfies (
43) for any
and
in
, so (
12) holds true. □
Finally, we shall obtain the characterisation of the global attractor in terms of the unstable and stable sets of the stationary points.
Theorem 4. Assume the conditions of Proposition 3. Then it holds thatwhereand denotes the set of all complete trajectories of (see Definition 3). Remark 4. In (45) it is equivalent to use instead of because all the solutions are bounded forward in time. Proof. We consider the function
Note that
is continuous in
. Indeed, the maps
and
are obviously continuous in
. On the other hand, by the embedding
and using Lebesgue’s theorem, the continuity of
follows.
By using (
29) and (
30) and multiplying Equation (
3) by
for any
, we can obtain the following energy equality:
Hence,
is non-increasing, and by
, (
10) and the boundedness of
, it is bounded from below. Thus
, and
, for some
Let
and
be such that
. By contradiction, suppose the existence of
and
, where
, for which
Since
is compact in
, we can take a converging subsequence (relabelled the same) such that
in
, where
. By the continuity of the function
E, it follows that
. We will obtain a contradiction by proving that
. Define
. By Lemma 1, there exist
and a subsequence satisfying
and
in
for
. Thus, from
we infer that
. Additionally,
satisfies the energy equality, so that
Therefore,
for a.a.
s, and then by Lemma 3 we have
. As a consequence,
. The converse inclusion follows from (
42).
As before, take arbitrary
and
satisfying
. Since by the embedding
the energy function is bounded from above in
,
, as
, for some
. Suppose that there are
and
, where
, such that
. Up to a subsequence we have that
in
,
. Moreover, for
there are
and a subsequence such that
and
in
for
. Therefore,
gives
and then by the above arguments we get a contradiction because
. Hence,
and we deduce the converse inclusion from (
42). □
Finally, we are able to obtain that the global attractor is compact in the space . This property will be important in order to study a more precise structure of the global attractor in terms of the stationary points and their heteroclinic connections.
We define the function
, where
, which is under the conditions of Proposition 3 (see [
18] for more details) a strong solution to the problem
Let
,
. We will prove first that the attractor is compact in any space
with
For this aim we will need the concept of mild solution. We consider the auxiliary problem
where
. The function
is called a mild solution to problem (
48) if
In the same way as in Lemma 2 in [
34] we obtain that a strong solution to problem (
47) is a mild solution to problem (
48) with
.
Lemma 4. Assume the conditions of Proposition 3. Then the global attractor is compact in for every
Proof. Let
be arbitrary. Since
is invariant, there exist
and
such that
and
for all
. Since
is a mild solution of (
48) with
, the variation of constants formula (
49) gives
As
is bounded in
(and then in
), condition (A6) and the continuity of
f imply that
where
does not depend on
z. The standard estimate
([
27] Theorem 37.5), implies that
so
is bounded in
for every
From the compact embedding for , and the fact that is closed in any we obtain the result. □
Corollary 2. Assume the conditions of Proposition 3. Then the global attractor is compact in
Proof. We obtain by Lemma 37.8 in [
27] the continuous embedding
Hence, the statement follows from Lemma 4. □