Next Article in Journal
Dynamics of Rössler Prototype-4 System: Analytical and Numerical Investigation
Next Article in Special Issue
A Full Description of ω-Limit Sets of Cournot Maps Having Non-Empty Interior and Some Economic Applications
Previous Article in Journal
Nonlocal Reaction–Diffusion Models of Heterogeneous Wealth Distribution
Previous Article in Special Issue
On a Retarded Nonlocal Ordinary Differential System with Discrete Diffusion Modeling Life Tables
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

About the Structure of Attractors for a Nonlocal Chafee-Infante Problem

by
Rubén Caballero
1,
Alexandre N. Carvalho
2,
Pedro Marín-Rubio
3 and
José Valero
1,*
1
Centro de Investigación Operativa, Universidad Miguel Hernández de Elche, Avenida Universidad s/n, 03202 Elche, Spain
2
Instituto de Ciências Matemáticas e de Computaçao, Universidade de São Paulo, Campus de São Carlos, Caixa Postal 668, 13560-970 São Carlos, Brazil
3
Departamento Ecuaciones Diferenciales y Análisis Numérico, Universidad de Sevilla, C/Tarfia, 41012 Sevilla, Spain
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(4), 353; https://doi.org/10.3390/math9040353
Submission received: 30 December 2020 / Revised: 4 February 2021 / Accepted: 5 February 2021 / Published: 10 February 2021
(This article belongs to the Special Issue Mathematical Methods on Economic Dynamics)

Abstract

:
In this paper, we study the structure of the global attractor for the multivalued semiflow generated by a nonlocal reaction-diffusion equation in which we cannot guarantee the uniqueness of the Cauchy problem. First, we analyse the existence and properties of stationary points, showing that the problem undergoes the same cascade of bifurcations as in the Chafee-Infante equation. Second, we study the stability of the fixed points and establish that the semiflow is a dynamic gradient. We prove that the attractor consists of the stationary points and their heteroclinic connections and analyse some of the possible connections.

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
t u α Δ u = f ( u ) c .
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]
t u α Δ u = u f ( u ) α R N g ( x y ) u ( y , t ) d y ;
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
u t a ( Ω g ( y ) u ( t , y ) d y ) Δ u = f ( t ) .
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
u t a ( Φ Ω ( u ( t ) ) Δ u = f ( t , u )
(here the functional Φ Ω may represent a general nonlocal functional acting over the whole domain Ω , for instance, u ( t ) H 0 1 2 or Ω g ( y ) u ( t , y ) d y ), 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
u t a ( u H 0 1 2 ) 2 u x 2 = λ f ( u )
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 H 0 1 . In Section 4 we prove that strong solutions generate a multivalued semiflow in H 0 1 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:
u t a ( u H 0 1 2 ) 2 u x 2 = λ f ( u ) + h ( t ) , t > 0 , x Ω , u ( t , 0 ) = u ( t , 1 ) = 0 , u ( 0 , x ) = u 0 ( x ) ,
where Ω = ( 0 , 1 ) and λ > 0 . Throughout the paper we will use the following conditions (but not all of them at the same time):
(A1)
f C ( R ) .
(A2)
f ( 0 ) = 0 .
(A3)
f ( 0 ) exists and f ( 0 ) = 1 .
(A4)
f is strictly concave if u > 0 and strictly convex if u < 0 .
(A5)
Growth and dissipation conditions: for p 2 , C i > 0 , i = 1 , . . , 4 , we have
| f ( u ) | C 1 + C 2 | u | p 1 ,
f ( u ) u C 3 C 4 | u | p , if p > 2 ,
lim sup u ± f ( u ) u 0 , if p = 2 .
(A6)
The function a C ( R + ) satisfies:
a ( s ) m > 0 .
(A7)
The function a C ( R + ) satisfies:
a ( s ) M 1 , s 0 ,
where M 1 > 0 .
(A8)
The function a C ( R + ) is non-decreasing.
(A9)
h L l o c 2 0 , + ; L 2 Ω .
(A10)
h does not depend on time and h L 2 Ω .
We define the function F ( u ) = 0 u f ( s ) d s . We observe that from (4) we have
| F ( s ) | C ˜ ( 1 + | s | p ) s R ,
whereas (5) implies
F ( s ) κ ˜ α ˜ 1 | s | p .
Additionally, from condition (6) it follows that for all ε > 0 , there exists a constant M > 0 such that f ( u ) u ε , for all | u | M . Hence, there exists m ε > 0 such that
f ( u ) u m ε + ε u 2 , u R .
In addition, it follows that
F ( u ) ε u 2 + C ε ,
where C ε > 0 . 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 h 0 .
Throughout the paper, · X will denote the norm in the Banach space X .

3. Existence of Solutions

In this section we will establish the existence of strong solutions for problem (3) with an initial condition in the phase space H 0 1 Ω . Although we will follow along the same lines as a similar result given in [18], we would like to point out that in the present case, as we are working in a one-dimensional problem, the assumptions for the function f are much weaker. In particular, we do not need to impose a growth assumption of any kind.
Definition 1.
For u 0 L 2 ( Ω ) , a weak solution to (3) is an element u L ( 0 , T ; L 2 ( Ω ) ) L 2 ( 0 , T ; H 0 1 ( Ω ) ) , for any T > 0 , such that
d d t ( u , v ) + a ( u H 0 1 2 ) ( u , v ) = λ ( f ( u ) , v ) + ( h ( t ) , v ) v H 0 1 ( Ω ) ,
where the equation is understood in the sense of distributions.
As usual, let A : D ( A ) H , D ( A ) = H 2 ( Ω ) H 0 1 Ω , be the operator A = d 2 d x 2 with Dirichlet boundary conditions. This operator is the generator of a C 0 -semigroup T ( t ) = e A t .
Definition 2.
For u 0 H 0 1 ( Ω ) , a strong solution to (3) is a weak solution with the extra regularity u L ( 0 , T ; H 0 1 ( Ω ) ) , u L 2 ( 0 , T ; D ( A ) ) and d u d t L 2 ( 0 , T ; L 2 ( Ω ) ) for any T > 0 .
Remark 1.
We observe that if u is a strong solution, then u C ( [ 0 , T ] ; H 0 1 ( Ω ) ) (see [27] p.102). This way, the initial condition makes sense.
Remark 2.
Since d u d t L 2 0 , T ; L 2 Ω for any strong solution, in this case equality (11) is equivalent to the following one:
0 T Ω d u t , x d t ξ t , x d x d t 0 T a ( u ( t ) H 0 1 2 ) Ω 2 u x 2 ξ d x d t = 0 T Ω λ f u t , x ξ t , x d x d t + 0 T Ω h t , x ξ t , x d x d t ,
for all ξ L 2 0 , T ; L 2 Ω .
Theorem 1.
Assume conditions (A1), (A6) and (A9). Assume also the existence of constants β , γ > 0 such that
f u u γ + β u 2 f o r a l l u R .
Then, for any u 0 H 0 1 ( Ω ) , problem (3) has at least one strong solution.
Remark 3.
Assumption (13) is weaker than the dissipative property (9) as the constant ε is arbitrarily small. Due to the fact that we are working in a one-dimensional domain, no growth condition of the type given in (A5) is necessary in order to prove existence of solutions. Additionally, (13) implies that
F u γ ˜ + β ˜ u 2
for some constants γ ˜ , β ˜ > 0 .
Proof. 
Consider a fixed value T > 0 . In order to use the Faedo–Galerkin method, let { w j } j 1 be the sequence of eigenfunctions of Δ in H 0 1 ( Ω ) with homogeneous Dirichlet boundary conditions, which forms a special basis of L 2 ( Ω ) . Since Ω is a bounded regular domain, it is known that { w j } H 0 1 ( Ω ) and that n N V n is dense in the spaces L 2 ( Ω ) and H 0 1 ( Ω ) , where V n = s p a n [ w 1 , , w n ] . As usual, P n will be the orthogonal projection in L 2 Ω , that is,
z n : = P n z = j = 1 n ( z , w j ) w j ,
and λ j will be the eigenvalues associated with the eigenfunctions w j . For each integer n 1 , we consider the Galerkin approximations
u n ( t ) = j = 1 n γ n j ( t ) w j ,
which are given by the following nonlinear ODE system:
d d t ( u n , w i ) + a ( u n H 0 1 2 ) ( u n , w i ) = λ ( f ( u n ) , w i ) + ( h , w i ) i = 1 , , n , u n ( 0 ) = P n u 0 .
We observe that P n u 0 u 0 in H 0 1 ( Ω ) . This Cauchy problem possesses a solution on some interval [ 0 , t n ) and by the estimates in the space L 2 ( Ω ) of the sequence { u n } given below for any T > 0 , such a solution can be extended to the whole interval [ 0 , T ] .
Firstly, multiplying the equation in (15) by γ n i ( t ) and summing from i = 1 to n, we obtain
1 2 d d t u n ( t ) L 2 2 + a ( u n H 0 1 2 ) u n ( t ) H 0 1 2 = λ ( f ( u n ( t ) , u n ( t ) ) + ( h ( t ) , u n ( t ) ) for a . e . t ( 0 , t n ) .
By using the Young and Poincaré inequalities we deduce that
( h ( t ) , u n ( t ) ) m 2 u n ( t ) H 0 1 2 + 1 2 λ 1 m h ( t ) L 2 2 ,
where m is the constant from (A6). Hence, from (A6), (13) and (16) it follows that
1 2 d d t u n ( t ) L 2 2 + m 2 u n ( t ) H 0 1 2 λ γ | Ω | + β λ u n ( t ) L 2 2 + 1 2 λ 1 m h ( t ) L 2 2 .
We infer that
u n ( t ) L 2 2 u n 0 L 2 2 e 2 β λ t + 0 t e 2 β λ t s 2 λ γ | Ω | + 1 λ 1 m h ( s ) L 2 2 d s u n 0 L 2 2 e 2 β λ T + K 1 T .
Therefore, the solution exists on any given interval [ 0 , T ] and
{ u n } is bounded in L ( 0 , T ; L 2 ( Ω ) ) .
Now, we multiply Equation (3) by d u n d t to obtain
d u n d t ( t ) L 2 2 + a ( u n H 0 1 2 ) 1 2 d d t u n H 0 1 2 = d d t Ω λ F ( u n ) d x + ( h ( t ) , d u n d t ) .
Introducing
A ( s ) = 0 s a ( r ) d r
we have
1 2 d u n d t ( t ) L 2 2 + d d t 1 2 A ( u n H 0 1 2 ) Ω λ F ( u n ) d x 1 2 h ( t ) L 2 2 .
Integrating the previous expression between 0 and t we get
1 2 A ( u n ( t ) H 0 1 2 ) + λ Ω F ( u n ( 0 ) ) d x + 1 2 0 t d d s u n ( s ) L 2 2 d s 1 2 A ( u n ( 0 ) H 0 1 2 ) + λ Ω F ( u n ( t ) ) d x + 1 2 0 t h ( s ) L 2 2 d s .
By (A6), (14) and (17) it follows that
m 2 u n ( t ) H 0 1 2 + λ Ω F ( u n ( 0 ) ) d x + 1 2 0 t d d s u n ( s ) L 2 2 d s 1 2 A ( u n ( 0 ) H 0 1 2 ) + λ β ˜ u n ( t ) L 2 2 + λ γ ˜ | Ω | + K 2 ( T ) 1 2 A ( u n ( 0 ) H 0 1 2 ) + λ β ˜ e 2 β λ T u n ( 0 ) L 2 2 + K 3 ( T ) .
Since d i m ( Ω ) = 1 , H 0 1 ( Ω ) L ( Ω ) , so u n 0 is bounded in L ( Ω ) . Thus, as f maps bounded sets of R into bounded ones, F u n 0 is bounded in L ( Ω ) as well. Therefore, we deduce that
{ u n } is bounded in L ( 0 , T ; H 0 1 ( Ω ) )
and
d u n d t is bounded in L 2 ( 0 , T ; L 2 ( Ω ) ) .
By using again the embedding H 0 1 ( Ω ) L ( Ω ) we obtain that u n is bounded in the space L ( 0 , T ; L ( Ω ) ) . Thus,
f ( u n ) is bounded in L ( 0 , T ; L ( Ω ) ) .
Additionally, we deduce that u n ( t ) H 0 1 2 is uniformly bounded in [ 0 , T ] , and then by the continuity of the function a · we get that the sequence a u n t H 0 1 2 is also uniformly bounded in [ 0 , T ] .
Finally, multiplying (15) by λ j γ n i ( t ) and summing from i = 1 to n we obtain
1 2 d d t u n H 0 1 2 + m Δ u n L 2 2 λ ( f ( u n ) , Δ u n ) + ( h ( t ) , Δ u ) .
By (23) and applying the Young inequality, we get
1 2 d d t u n H 0 1 2 + m Δ u n L 2 2 λ 2 m f ( u n ) L 2 2 + m 4 Δ u n L 2 2 + 1 m h ( t ) L 2 2 + m 4 Δ u L 2 2 .
Integrating the previous expression between 0 and t, it follows that
u n ( t ) H 0 1 2 + m 0 t Δ u n ( s ) L 2 2 d s u n ( 0 ) H 0 1 2 + 2 λ 2 m 0 t f ( u n ( s ) ) L 2 2 d s + 2 m 0 t h ( s ) L 2 2 d s .
Taking into account (23), the last inequality implies that
u n is bounded in L 2 ( 0 , T ; D ( A ) ) ,
so { Δ u n } and { a ( u n H 0 1 2 ) Δ u n } are bounded in L 2 ( 0 , T ; L 2 ( Ω ) ) .
As a consequence, there exist u L ( 0 , T ; H 0 1 ( Ω ) ) and a subsequence u n (relabelled the same) such that
u n * u in L ( 0 , T ; H 0 1 ( Ω ) ) , u n u in L 2 ( 0 , T ; D ( A ) ) , f ( u n ) * χ in L ( 0 , T ; L ( Ω ) ) , a ( u n H 0 1 2 ) * b in L ( 0 , T ) ,
where ⇀ ( * ) stands for the weak (weak star) convergence. By (22) and (24) the Aubin–Lions compactness lemma gives that u n u in L 2 ( 0 , T ; H 0 1 ( Ω ) ) , so u n ( t ) u ( t ) in H 0 1 ( Ω ) a.e. on ( 0 , T ) . Consequently, there exists a subsequence u n , relabelled the same, such that u n ( t , x ) u ( t , x ) a.e. in Ω × ( 0 , T ) .
Moreover, thanks to the inequality
u n ( t 2 ) u n ( t 1 ) L 2 2 = t 1 t 2 d d t u n ( s ) d s L 2 2 d d t u n L 2 ( 0 , T ; L 2 ( Ω ) ) 2 | t 2 t 1 | t 1 , t 2 [ 0 , T ] ,
(21), (22) and H 0 1 ( Ω ) L 2 ( Ω ) , the Ascoli-Arzelà theorem implies that { u n } converges strongly in C ( [ 0 , T ] ; L 2 ( Ω ) ) for all T > 0 . Therefore, we obtain from (21) that u n ( t ) u ( t ) in H 0 1 ( Ω ) , for any t 0 .
Additionally, by (25) we have that P n f ( u n ) ) χ in L q ( 0 , T ; L q ( Ω ) ) for any q 1 (see [28] p.224). Since f is continuous, it follows that f ( u n ( t , x ) ) f ( u ( t , x ) ) a.e. in Ω × ( 0 , T ) . Therefore, in view of (25), by ([29] Lemma 1.3) we have that χ = f ( u ) .
As a consequence, by the continuity of a we get that
a ( u n ( t ) H 0 1 2 ) a ( u ( t ) H 0 1 2 ) a . e . on ( 0 , T ) .
Since the sequence is uniformly bounded, by Lebesgue’s theorem this convergence takes place in L 2 ( 0 , T ) , so b = a ( u H 0 1 2 ) . Thus,
a ( u n H 0 1 2 ) Δ u n a ( u H 0 1 2 ) Δ u , in L 2 ( 0 , T ; L 2 ( Ω ) ) .
Therefore, we can pass to the limit to conclude that u is a strong solution.
It remains to show that u ( 0 ) = u 0 which makes sense since u C ( [ 0 , T ] ; H 0 1 ( Ω ) ) (see Remark 4). Indeed, let be ϕ C 1 ( [ 0 , T ] ; H 0 1 ( Ω ) ) with ϕ ( T ) = 0 , ϕ ( 0 ) 0 . We multiply the equation in (3) and (15) by ϕ and integrate by parts in the t variable to obtain that
0 T ( u ( t ) , ϕ ( t ) ) a ( u ( t ) H 0 1 2 ) ( Δ u ( t ) , ϕ ( t ) ) d t = 0 T ( λ f ( u ( t ) ) + h ( t ) , ϕ ( t ) ) d t + ( u ( 0 ) , ϕ ( 0 ) ) ,
0 T ( u n ( t ) , ϕ ( t ) ) a ( u n ( t ) H 0 1 2 ) ( Δ u n ( t ) , ϕ ( t ) ) d t = 0 T ( λ f ( u n ( t ) ) + h ( t ) , ϕ ( t ) ) d t + ( u n ( 0 ) , ϕ ( 0 ) ) .
In view of the previous convergences, we can pass to the limit in (27). Taking into account (26) and bearing in mind u n ( 0 ) = P n u 0 u 0 , since ϕ 0 H 0 1 ( Ω ) is arbitrary, we infer that u ( 0 ) = u 0 . □

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 d i s t X C , D the Hausdorff semidistance from C to D, that is,
d i s t X ( C , D ) = sup c C inf d D ρ c , d .
Let us consider the phase space X = H 0 1 Ω and the sets
K u 0 = { u ( · ) : u is a strong solution of ( 3 ) such that u 0 = u 0 } ,
R = u 0 X K u 0 .
Denote by P ( X ) the class of nonempty subsets of X. We define the (possibly multivalued) map G : R + × X P ( X ) by
G ( t , u 0 ) = { u ( t ) : u R and u ( 0 ) = u 0 } .
In order to study the map G let us consider the following axiomatic properties of the set R :
(K1)
For every x X there is ϕ R satisfying ϕ ( 0 ) = x .
(K2)
ϕ τ ( · ) : = ϕ ( · + τ ) R for every τ 0 and ϕ R (translation property).
(K3)
Let ϕ 1 , ϕ 2 R be such that ϕ 2 ( 0 ) = ϕ 1 ( s ) for some s > 0 . Then, the function ϕ defined by
ϕ ( t ) = ϕ 1 ( t ) 0 t s , ϕ 2 ( t s ) s t ,
belongs to R (concatenation property).
(K4)
For every sequence { ϕ n } R satisfying ϕ n ( 0 ) x 0 in X, there is a subsequence { ϕ n k } and ϕ R such that ϕ n k ( t ) ϕ ( t ) for every t 0 .
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 R fulfilling (K1) and (K2) gives rise to a multivalued semiflow G through (28) (m-semiflow for short), which means that:
  • G ( 0 , x ) = x for all x X ;
  • G ( t + s , x ) G ( t , G ( s , x ) ) for all t , s 0 and x X .
Moreover, (K3) implies that the m-semiflow is strict, that is, G ( t + s , x ) = G ( t , G ( s , x ) ) for all t , s 0 and x X .
We will show first that the m-semiflow G possesses a bounded absorbing set in the space L 2 Ω and that property (K4) is satisfied.
Lemma 1.
Assume conditions (A1), (A6), (A10) and (13). Given { u n } R , u n ( 0 ) u 0 weakly in H 0 1 ( Ω ) , there exists a subsequence of { u n } (relabelled the same) and u K ( u 0 ) such that
u n ( t ) u ( t ) i n H 0 1 ( Ω ) , t > 0 .
Additionally, if u n ( 0 ) u 0 strongly in H 0 1 ( Ω ) , then for t n 0 we get u n ( t n ) u 0 strongly in H 0 1 ( Ω ) .
Proof. 
Since d u n d t L 2 ( 0 , T ; L 2 ( Ω ) ) and u n L 2 ( 0 , T ; H 0 1 ( Ω ) ) , we have by ([27] p. 102) that
d d t u n H 0 1 2 = 2 ( Δ u n , u t n ) for a . a . t
and u n C ( [ 0 , T ] ; H 0 1 ( Ω ) ) . Additionally, as f ( u n ) L ( 0 , T ; L ( Ω ) ) , by regularization one can show that ( F ( u n ( t ) ) , 1 ) is an absolutely continuous function on [ 0 , T ] and
d d t ( F ( u n ( t ) ) , 1 ) = ( f ( u n ( t ) ) , d u n d t ) for a . a . t > 0 .
By a similar argument as in Theorem 1, there is a subsequence of u n such that
u n is bounded in L ( 0 , T ; L ( Ω ) ) , u n is bounded in L ( 0 , T ; H 0 1 ( Ω ) ) , f ( u n ) is bounded in L ( 0 , T ; L ( Ω ) ) , u n is bounded in L 2 ( 0 , T ; D ( A ) ) .
Therefore, arguing as in the proof of Theorem 1, there exist u K u 0 and a subsequence u n , relabelled the same, such that
u n * u in L ( 0 , T ; H 0 1 ( Ω ) ) u n u in L 2 ( 0 , T ; D ( A ) ) f ( u n ) * f ( u ) in L ( 0 , T ; L ( Ω ) ) d u n d t d u d t in L 2 ( 0 , T ; L 2 ( Ω ) ) a ( u n H 0 1 2 ) Δ u n a ( u H 0 1 2 ) Δ u in L 2 ( 0 , T ; L 2 ( Ω ) ) , u n u in L 2 ( 0 , T ; H 0 1 ( Ω ) ) , u n u in C ( [ 0 , T ] , L 2 ( Ω ) ) , u n ( t ) u ( t ) in H 0 1 ( Ω ) t ( 0 , T ] .
We also need to prove that u n ( t ) u ( t ) in H 0 1 ( Ω ) for all t ( 0 , T ] . To that end, we multiply (3) by u t n , and using (A10), (29) and (31) we have
1 2 d u n d t L 2 2 + d d t 1 2 A ( u n ( t ) H 0 1 2 ) C .
Thus, we obtain
A ( u n ( t ) H 0 1 2 ) A ( u n ( s ) H 0 1 2 ) + 2 C ( t s ) , t s 0 .
Since this inequality is also true for u ( · ) , the functions Q n ( t ) = A ( u n ( t ) H 0 1 2 ) 2 C t , Q ( t ) = A ( u ( t ) H 0 1 2 ) 2 C t are continuous and non-increasing in [ 0 , T ] . Moreover, from (32) we deduce that
Q n ( t ) Q ( t ) for a . e . t ( 0 , T ) .
Take 0 < t T and 0 < t j < t such that t j t and Q n ( t j ) Q ( t j ) for all j. Then
Q n ( t ) Q ( t ) Q n ( t j ) Q ( t ) | Q n ( t j ) Q ( t j ) | + | Q ( t j ) Q ( t ) | .
For any δ > 0 there exist j ( δ ) and N ( j ( δ ) ) such that Q n ( t ) Q ( t ) δ if n N . Then lim sup Q n ( t ) Q ( t ) , so lim sup u n ( t ) H 0 1 2 u ( t ) H 0 1 2 , which follows by contradiction using the continuity of the function A ( s ) . As u n ( t ) u ( t ) weakly in H 0 1 ( Ω ) implies that lim inf u n ( t ) H 0 1 2 u ( t ) H 0 1 2 , we obtain
u n ( t ) H 0 1 2 u ( t ) H 0 1 2 ,
so that u n ( t ) u ( t ) strongly in H 0 1 ( Ω ) .
Finally, if u n ( 0 ) u 0 strongly in H 0 1 ( Ω ) and we take t n 0 , then
Q n ( t n ) Q ( 0 ) Q n ( 0 ) Q ( 0 ) = A ( u n 0 H 0 1 2 ) A ( u 0 H 0 1 2 ) 0 ,
so lim sup Q n ( t n ) Q ( 0 ) . Repeating the above argument, we infer that u n ( t n ) u 0 strongly in H 0 1 ( Ω ) . □
Corollary 1.
Assume the conditions of Lemma 1. Then the set R satisfies condition ( K 4 ) .
The map t G ( t , x ) is said to be upper semicontinuous if for every x X and for an arbitrary neighbourhood O ( G ( t , x ) ) in X there is δ > 0 such that as soon as d ( y , x ) < δ , we have G ( t , y ) O .
Proposition 1.
Assume the conditions of Lemma 1. The multivalued semiflow G is upper semicontinuous for all t 0 . Additionally, it has compact values.
Proof. 
By contradiction let us assume that there exist t 0 , u 0 H 0 1 ( Ω ) , a neighbourhood O ( G ( t , u 0 ) ) and sequences { y n } , { u 0 n } such that y n G ( t , u 0 n ) , u 0 n converges strongly to u 0 in H 0 1 ( Ω ) and y n O ( G ( t , u n ) ) for all n N . Thus, there exists u n K ( u 0 n ) such that y n = u n ( t ) . From Lemma 1 there exists a subsequence of y n which converges to some y G ( t , u 0 ) . This contradicts y n O ( G ( t , u 0 ) ) for any n N . □
In order to prove the existence of an absorbing set in the space L 2 Ω 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 L 2 Ω ; that is, there exists a constant K > 0 such that for any R > 0 there is a time t 0 = t 0 ( R ) such that
y L 2 K f o r a l l t t 0 , y G ( t , u 0 ) ,
where u 0 L 2 R . Moreover, there is L > 0 such that
t t + 1 u ( s ) H 0 1 2 d s L f o r a l l t t 0 , u K u 0 .
Proof. 
By multiplying Equation (3) by u and using (A6) and (9), we get
1 2 d d t u ( t ) L 2 2 + m u ( t ) H 0 1 2 ( f ( u ) , u ) + ( h , u ) m ε | Ω | + ε u ( t ) L 2 2 + 1 2 λ 1 m h L 2 2 + λ 1 m 2 u L 2 2 .
By using the Poincaré inequality it follows that
d d t u L 2 2 2 m ε | Ω | + 2 ( ε m 2 λ 1 ) u ( t ) L 2 2 + 1 λ 1 m h L 2 2 = δ u ( t ) L 2 2 + κ ,
where δ = m λ 1 2 ε , κ = 2 m ε | Ω | + 1 λ 1 m h L 2 2 . We take a small enough ε > 0 so that δ > 0 . Then Gronwall’s lemma gives
u ( t ) L 2 2 u ( 0 ) L 2 2 e δ t + κ δ .
Hence, taking
t t 0 = 1 δ ln δ R 2 κ
we get (33) for K = 2 κ δ .
On the other hand, using again the Poincaré inequality from (35) we get
d d t u ( t ) L 2 2 + m λ 1 2 ε λ 1 u ( t ) H 0 1 2 κ
and integrating from t to t + 1 we obtain
m λ 1 2 ε λ 1 t t + 1 u ( s ) H 0 1 2 d s u ( t ) L 2 2 + κ .
Therefore, applying (33) and (34) follows. □
Further, in order to obtain an absorbing set in H 0 1 Ω we need to assume additionally that either the function a · 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 B 1 for G, which is compact in H 0 1 ( Ω ) .
Proof. 
In view of Proposition 2 we have an absorbing set B 0 in L 2 ( Ω ) . Let K > 0 be such that y K for all y B 0 .
Through multiplying (3) by u and using (9) and (36) we get
d d t u ( t ) L 2 2 + a u ( t ) H 0 1 2 u ( t ) H 0 1 2 2 m ε | Ω | + 2 ε u ( t ) L 2 2 + 1 λ 1 m h L 2 2 K 1 + K 2 u ( 0 ) L 2 2 .
Thus, integrating between t and t + r , 0 < r 1 , we deduce by using (36) again that
u ( t + r ) L 2 2 + t t + r a u ( s ) H 0 1 2 u ( s ) H 0 1 2 d s K 1 + K 2 u ( 0 ) L 2 2 + u ( t ) L 2 2 K 3 u ( 0 ) L 2 2 + K 4 .
Additionally, if p > 2 in (A5), we multiply again by (3) by u and use (5) and (A6) to obtain
1 2 d d t u ( t ) L 2 2 + m 2 u ( t ) H 0 1 2 + C 4 u t L p p C 3 + 1 2 λ 1 m h L 2 2 .
Integrating over t , t + r we have
u ( t + r ) L 2 2 + 2 C 4 t t + r u s L p p d s K 5 + u ( t ) L 2 2 K 6 + u ( 0 ) L 2 2 .
If we assume (A7), by (37) and (A6) we have that
t t + r A ( u ( s ) H 0 1 2 ) d s t t + r M 1 u ( s ) H 0 1 2 d s K 7 ( 1 + u ( 0 ) L 2 2 ) .
If we assume (A8), by (37) we obtain
t t + r A ( u ( s ) H 0 1 2 ) d s = t t + r 0 u ( s ) H 0 1 2 a r d r d s t t + r a u ( s ) H 0 1 2 u ( s ) H 0 1 2 d s K 3 u ( 0 ) L 2 2 + K 4 .
On the other hand, by (7) we get
Ω F u t d x C ˜ Ω ( 1 + | u t | p ) d x .
By using (29) and (30) we can argue as in Theorem 1 to obtain
1 2 u t L 2 2 + d d t 1 2 A ( u ( t ) H 0 1 2 Ω λ F ( u n ) d x 1 2 h L 2 2 .
Since (38)–(41) imply that
t t + r 1 2 A ( u ( s ) H 0 1 2 Ω λ F ( u s ) d x d s K 8 + K 9 u ( 0 ) L 2 2 ,
we can apply the uniform Gronwall lemma to get
1 2 A ( u ( t + r ) H 0 1 2 ) Ω λ F ( u t + r ) d x K 8 + K 9 u ( 0 ) L 2 2 r + K 10 , for all t 0 ,
so by condition ( A 6 ) , (10) and (36) it follows that
u ( t + 1 ) H 0 1 2 K 11 + K 12 u ( 0 ) L 2 2 ,
for all t 0 . In particular,
u ( 1 ) H 0 1 2 K 11 + K 12 u ( 0 ) L 2 2 ,
for any strong solution u ( · ) with initial condition u ( 0 ) .
For any u 0 H 0 1 Ω with u 0 H 0 1 R and any u R such that u 0 = u 0 , the semiflow property G ( t + 1 , u 0 ) G ( 1 , G ( t , u 0 ) ) and G ( t , u 0 ) B 0 , if t t 0 R , imply that
u t + 1 H 0 1 2 C ( 1 + K 2 ) t t 0 R .
Then there exists M > 0 such that the closed ball B M in H 0 1 Ω centred at 0 with radius M is absorbing for G.
By Lemma 1 the set B 1 = G ( 1 , B M ) ¯ is an absorbing set which is compact in H 0 1 Ω . □
Given an m-semiflow G , a set B X is said to be negatively (positively) invariant if B G ( t , B ) ( G ( t , B ) B ) for all t 0 , and strictly invariant (or, simply, invariant) if it is both negatively and positively invariant.
We recall that a set A X is called a global attractor for the m-semiflow G if it is negatively invariant and attracts all bounded subsets; i.e., d i s t X ( G ( t , B ) , A ) 0 as t + . When A 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 A .
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 H 0 1 Ω . Therefore, by ([32] Theorem 4 and Remark 8) the existence of the global invariant attractor and its compactness in H 0 1 Ω follow. □
We recall some concepts which are necessary to study the structure of the global attractor.
Definition 3.
A map ϕ : R X is a complete trajectory of R if ϕ ( · + s ) [ 0 , ) R for all s R . It is a complete trajectory of G if ϕ ( t + s ) G ( t , ϕ ( s ) ) for every s R , t 0 .
An element z X is a fixed point of R if φ ( · ) z R . We denote the set of all fixed points by R R .
An element z X is a fixed point of G if z G ( t , z ) for every t 0 .
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 R is also a fixed point (complete trajectory) of G.
Let (K1)-(K4) hold. Then the fixed points of R and G are the same. In addition, a map ϕ : R X is a complete trajectory of R 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 A . Then
A = { γ ( 0 ) : γ K } = t R { γ ( t ) : γ K } ,
where K stands for the set of all bounded complete trajectories in R .
In view of Theorem 3, as R 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 x B there exists a complete trajectory γ of R contained in B such that γ ( 0 ) = x . Characterisation (42) implies that the attractor A is weakly invariant.
The set of fixed points R R is characterised as follows.
Lemma 3.
Assume the conditions of Lemma 1. Let R be the set of z H 2 ( Ω ) H 0 1 ( Ω ) such that
a ( z H 0 1 2 ) d 2 z d x 2 = λ f ( z ) + h i n L 2 ( Ω ) .
Then R R = R .
Proof. 
If z R R , then u ( t ) z R . Thus, u ( · ) satisfies (12) and d u d t = 0 in L 2 ( 0 , T ; L 2 ( Ω ) ) , so (43) is satisfied. Let z R . Then the map u ( t ) z satisfies (43) for any t 0 and d u d t = 0 in L 2 ( 0 , T ; L 2 ( Ω ) ) , 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 that
A = M + ( R ) = M ( R ) ,
where
M + ( R ) = { z : γ ( · ) K , γ ( 0 ) = z , dist H 0 1 ( γ ( t ) , R ) 0 , t + } ,
M ( R ) = { z : γ ( · ) F , γ ( 0 ) = z , dist H 0 1 ( γ ( t ) , R ) 0 , t } ,
and F denotes the set of all complete trajectories of R (see Definition 3).
Remark 4.
In (45) it is equivalent to use K instead of F because all the solutions are bounded forward in time.
Proof. 
We consider the function E : A R
E ( y ) = 1 2 A ( y H 0 1 2 ) λ Ω F ( y ( x ) ) d x Ω h ( x ) y ( x ) d x .
Note that E ( y ) is continuous in H 0 1 ( Ω ) . Indeed, the maps y 1 2 A ( y H 0 1 2 ) and y Ω h x y x d x are obviously continuous in H 0 1 ( Ω ) . On the other hand, by the embedding H 0 1 ( Ω ) L ( Ω ) and using Lebesgue’s theorem, the continuity of y Ω F ( y ( x ) ) d x follows.
By using (29) and (30) and multiplying Equation (3) by d u d t for any u R , we can obtain the following energy equality:
s t d d r u ( r ) L 2 2 d r + E ( u ( t ) ) = E ( u ( s ) ) for all t s 0 .
Hence, E ( u ( t ) ) is non-increasing, and by ( A 6 ) , (10) and the boundedness of A , it is bounded from below. Thus E ( u ( t ) ) l , and t + , for some l R .
Let z A and u K be such that u ( 0 ) = z . By contradiction, suppose the existence of ε > 0 and u ( t n ) , where t n + , for which d i s t H 0 1 ( u ( t n ) , R ) > ε . Since A is compact in H 0 1 ( Ω ) , we can take a converging subsequence (relabelled the same) such that u ( t n ) y in H 0 1 ( Ω ) , where t n . By the continuity of the function E, it follows that E ( y ) = l . We will obtain a contradiction by proving that y R . Define v n · = u ( · + t n ) . By Lemma 1, there exist v R and a subsequence satisfying v ( 0 ) = y and v n ( t ) v ( t ) in H 0 1 ( Ω ) for t 0 . Thus, from E ( v n ( t ) ) E ( v ( t ) ) we infer that E ( v t ) = l . Additionally, v ( · ) satisfies the energy equality, so that
l + 0 t v r L 2 2 d r = E ( v t ) + 0 t v r L 2 2 d r = E ( v ( 0 ) ) = E ( y ) = l .
Therefore, d v d t ( s ) = 0 for a.a. s, and then by Lemma 3 we have y R R = R . As a consequence, A M + ( R ) . The converse inclusion follows from (42).
As before, take arbitrary z A and u K satisfying u ( 0 ) = z . Since by the embedding H 0 1 ( Ω ) C ( [ 0 , 1 ] ) the energy function is bounded from above in A , E ( u ( t ) ) l , as t , for some l R . Suppose that there are ε > 0 and u ( t n ) , where t n + , such that d i s t H 0 1 ( u ( t n ) , R ) > ε . Up to a subsequence we have that u ( t n ) y in H 0 1 ( Ω ) , E ( y ) = l . Moreover, for v n ( · ) = u ( · t n ) there are v R and a subsequence such that v ( 0 ) = y and v n ( t ) v ( t ) in H 0 1 ( Ω ) for t 0 . Therefore, E ( v n ( t ) ) E ( v ( t ) ) gives E ( v ( t ) ) = l and then by the above arguments we get a contradiction because y R . Hence, A M ( R ) and we deduce the converse inclusion from (42). □
Finally, we are able to obtain that the global attractor is compact in the space C 1 [ 0 , 1 ] . 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 w t = u α 1 t , where α ( t ) = 0 t a ( u ( s ) H 0 1 2 ) d s , which is under the conditions of Proposition 3 (see [18] for more details) a strong solution to the problem
w t 2 w x 2 = f ( w ) + h a ( w H 0 1 2 ) , in ( 0 , ) × Ω , w = 0 on ( 0 , ) × Ω , w ( 0 , x ) = u 0 ( x ) in Ω .
Let V 2 r = D ( A r ) , r 0 . We will prove first that the attractor is compact in any space V 2 r with 0 r < 1 . For this aim we will need the concept of mild solution. We consider the auxiliary problem
d v d t + A v ( t ) = g t , t > 0 , v 0 = u 0 ,
where g L l o c 2 0 , + ; L 2 Ω . The function u C ( [ 0 , + ) , L 2 Ω ) is called a mild solution to problem (48) if
v t = e A t u 0 + 0 t e A t s g ( s ) d s , t 0 .
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 g t = f ( w t ) + h / a ( w t H 0 1 2 ) .
Lemma 4.
Assume the conditions of Proposition 3. Then the global attractor A is compact in V 2 r for every 0 r < 1 .
Proof. 
Let z A be arbitrary. Since A is invariant, there exist u 0 A and u R such that z = u 1 and u t A for all t 0 . Since w t = u α 1 t is a mild solution of (48) with g t = f ( w t ) + h / a ( w t H 0 1 2 ) , the variation of constants formula (49) gives
z = w ( α 1 ) = e A α 1 u 0 + 0 α 1 e A α 1 s g ( s ) d s .
As A is bounded in H 0 1 Ω (and then in L Ω ), condition (A6) and the continuity of f imply that
u 0 L 2 C , g L 0 , α 1 ; L 2 Ω C ,
where C > 0 does not depend on z. The standard estimate e A t L ( L 2 Ω , D ( A r ) ) M r t r e a t , M r , a > 0 ([27] Theorem 37.5), implies that
A r z L 2 A r e A α 1 u 0 L 2 + 0 α 1 A r e A α 1 s g ( s ) L 2 d s M r e a α 1 α 1 r C + M r C 0 α 1 α 1 s r d s ,
so A is bounded in V 2 r for every 0 r < 1 .
From the compact embedding V α V β , for α > β , and the fact that A is closed in any V 2 r we obtain the result. □
Corollary 2.
Assume the conditions of Proposition 3. Then the global attractor A is compact in C 1 ( [ 0 , 1 ] ) .
Proof. 
We obtain by Lemma 37.8 in [27] the continuous embedding
V 2 r C 1 ( [ 0 , 1 ] ) if r > 3 4 .
Hence, the statement follows from Lemma 4. □

5. Fixed Points

In this section we are interested in studying the fixed points of problem (3) when h 0 , that is, the solutions of the boundary-value problem
a ( u H 0 1 2 ) d 2 u d x 2 = λ f ( u ) , 0 < x < 1 , u 0 = u 1 = 0 .
For this aim we will use the properties of the fixed points of the standard Chafee-Infante equation. In order to do that, for any d 0 we will study the )–(Aollowing boundary-value problem.
a ( d ) d 2 u d x 2 = λ f ( u ) , 0 < x < 1 , u 0 = u 1 = 0 ,
as it is obvious that u · is solution to problem (50) if and only if u · is a solution to problem (51) with d = u H 0 1 2 .

5.1. Dependence on the Parameters of the Fixed Points for the Chafee-Infante Equation

Denoting λ ˜ = λ a d problem (51) becomes
d 2 u d x 2 = λ ˜ f ( u ) , 0 < x < 1 , u 0 = u 1 = 0 .
Assuming conditions (A1)–(A5), it is known [35] that if n 2 π 2 < λ ˜ n + 1 2 π 2 , then this problem has exactly 2 n + 1 solutions, denoted by v 0 0 , v 1 ± , . . . , v n ± . The function v k ± has k + 1 simple zeros in [ 0 , 1 ] .
We need to study the dependence of the norm of these fixed points on the parameter λ ˜ . First, we will show that the H 1 -norm of the fixed points of problem (52) is strictly increasing with respect to the parameter λ ˜ .
Lemma 5.
Assume conditions (A1)–(A5). Let v 1 = v k , λ 1 + , v 2 = v k , λ 2 + with k 2 π 2 < λ 1 < λ 2 . Then v 1 H 0 1 < v 2 H 0 1 .
Proof. 
We consider the equivalent norm in H 0 1 Ω given by v L 2 . The fixed points are the solutions of the initial value problem
d 2 u d x 2 + λ ˜ f ( u ) = 0 , u ( 0 ) = 0 , u ( 0 ) = v 0
such that u 1 = 0 . The solutions of (53) satisfy the relation
( u ( x ) ) 2 2 + λ ˜ F ( u ( x ) ) = λ ˜ E , 0 x 1 ,
for some constant E 0 . Denote u λ ˜ = v k , λ ˜ + . By Theorem 7 in [35] we have that u λ ˜ is associated with a unique value E = E k + ( λ ˜ ) > 0 . Moreover, E k + ( λ ˜ ) is a solution of one of the following equations:
m τ + λ ˜ ( E ) + ( m 1 ) τ λ ˜ ( E ) = 1 2 , m τ λ ˜ ( E ) + ( m 1 ) τ + λ ˜ ( E ) = 1 2 , m τ + λ ˜ ( E ) + m τ λ ˜ ( E ) = 1 2 ,
where either k = 2 m 1 or k = 2 m and
τ + λ ˜ ( E ) = λ ˜ 1 / 2 0 U + ( E ) ( E F ( u ) ) 1 / 2 d u ,
τ λ ˜ ( E ) = λ ˜ 1 / 2 U ( E ) 0 ( E F ( u ) ) 1 / 2 d u ,
being U + ( E ) ( U ( E ) ) the positive (negative) inverse of F at E. It is obvious that for E fixed the functions τ + λ ˜ ( E ) , τ λ ˜ ( E ) are strictly decreasing with respect to λ ˜ . Then from (55) we deduce that the root E k + ( λ ˜ ) is strictly increasing with respect to λ ˜ . Thus, If λ 1 < λ 2 , we have
2 λ 1 ( E k + ( λ 1 ) F ( u ) ) < 2 λ 2 ( E k + ( λ 2 ) F ( u ) ) , U ( E k + ( λ 1 ) ) u U + ( E k + ( λ 1 ) ) .
We will prove now that u λ ˜ L 2 is strictly increasing in λ ˜ .
The function u λ ˜ has k + 1 simple zeros in [ 0 , 1 ] and u λ ˜ is positive in the first subinterval. Let T + ( E k + ( λ ) ) be the x-time necessary to go from the initial condition u λ ( 0 ) = 0 to the point where u λ ( T + ( E k + ( λ ) ) ) = 0 . Then the length of the first subinterval is 2 T + ( E k + ( λ ) ) [35]. By (54),
( u λ ˜ ( x ) ) 2 = 2 λ ˜ E k + ( λ ˜ ) F ( u λ ˜ ( x ) ) u λ ˜ ( x ) ,
so we have
0 T + ( E k + ( λ ˜ ) ) ( u λ ˜ ( x ) ) 2 d x = 0 T + ( E k + ( λ ˜ ) ) 2 λ ˜ E k + ( λ ˜ ) F ( u λ ˜ ( x ) ) u λ ˜ ( x ) d x .
By the change of variable v = u λ ˜ ( x ) we obtain
0 T + ( E k + ( λ ˜ ) ) ( u λ ˜ ( x ) ) 2 d x = 0 U + ( E k + ( λ ˜ ) ) 2 λ ˜ E k + ( λ ˜ ) F ( v ) d v = g ( λ ˜ ) .
Since λ ˜ U + ( E k + ( λ ˜ ) ) is strictly increasing and by using (58), we conclude that the function g ( λ ˜ ) is strictly increasing. Hence, by putting x 1 ( λ ˜ ) = 2 T + ( E k + ( λ ˜ ) ) we obtain that the norm of u λ ˜ in the first subinterval, u λ ˜ L 2 ( 0 , x 1 ( λ ˜ ) ) , is strictly increasing. By arguing in the same way as for the other subintervals, we obtain that λ ˜ u λ ˜ L 2 is strictly increasing. □
Let us prove the same result but with respect to the norm u λ ˜ L p with p 1 .
Lemma 6.
Assume conditions (A1)–(A5) and let f be odd. Let v 1 = v k , λ 1 + , v 2 = v k , λ 2 + with k 2 π 2 < λ 1 < λ 2 . Then v 1 L p < v 2 L p for any p 1 .
Proof. 
As in the previous lemma, denote u λ ˜ = v k , λ ˜ + . The function u λ ˜ has k + 1 zeros in [ 0 , 1 ] at the points 0 < x 1 < x 2 < . . . < x k 1 < 1 . When f is odd, by symmetry, the length of all subintervals has to be the same, so x j = j k regardless the value of λ ˜ .
We shall prove that in the first subinterval we have that u λ 1 x < u λ 2 x , for all x 0 , 1 k . By (54) for x [ 0 , 1 2 k ] we have
x = 0 x d s = 0 u λ ˜ x d u 2 λ ˜ E k + λ ˜ F u ,
so (58) yields
x = 0 u λ 2 x d u 2 λ 2 E k + λ 2 F u = 0 u λ 1 x d u 2 λ 1 E k + λ 1 F u > 0 u λ 1 x d u 2 λ 2 E k + λ 2 F u , if x ( 0 , 1 2 k ] .
Thus, u λ 1 x < u λ 2 x , for all x ( 0 , 1 2 k ] . By symmetry we obtain that the inequality is true in 0 , 1 k .
Repeating the same argument in the other subintervals we get that
u λ 1 x < u λ 2 x for all x 0 , 1 , x j k , j = 1 , . . . k 1 .
This implies that u λ 1 L p < u λ 2 L p for any p 1 .
Remark 5.
The statements in Lemmas 5 and 6 are also true for v k , λ ˜ , because v k , λ ˜ x = v k , λ ˜ + 1 x , so the H 0 1 and L p norms of v k , λ ˜ and v k , λ ˜ + are the same.

5.2. Nonlocal Fixed Points

Although in this paper we are mainly interested in problem (3), we will study the existence of stationary points for an elliptic problem with a more general nonlocal term than in (50). Namely, let us consider the following problem:
a l ( u ) u x x = λ f u , 0 < x < 1 , u 0 = u ( 1 ) = 0 ,
where
l u = u H 0 1 r or u L p r , p 1 , r > 0 .
Let
d k = sup { d : λ > a d ¯ π 2 k 2 d ¯ d } .
Then for any d < d k there exists the fixed point u k d of (51), where u k d is either equal to u k + or u k .
It is obvious that any solution of (59) is a solution of (51) with d = l u . Therefore, all the solutions to problem (59) have to be solutions u k d to problem (51) for a suitable d .
Theorem 5.
Assume conditions (A1)–(A6) and, additionally, that
a 0 π 2 k 2 < λ .
Then:
  • For any 1 j k there exists d j * < d k such that u j d j * is a fixed point of problem (59).
  • If λ a 0 π 2 k + 1 2 and a ( 0 ) = min s 0 { a s } , there are no fixed points for j > k .
  • If N k is the first integer such that λ inf s 0 { a s π 2 N + 1 2 } , there are no fixed points for j > N .
  • If l u = u H 0 1 r , λ a 0 π 2 k + 1 2 and a is non-decreasing, there are exactly 2 k + 1 solutions to problem (59): 0 , u 1 , d 1 * ± , . . . , u k , d k * ± .
  • If l u = u L p r , λ a 0 π 2 k + 1 2 , f is odd and a is non-decreasing, there are exactly 2 k + 1 solutions to problem (59): 0 , u 1 , d 1 * ± , . . . , u k , d k * ± .
Proof. 
For the first statement, it is enough to prove the result for j = k . By condition (60) we have that d k ( 0 , + ] .
Consider first the case where d k is finite. We need to obtain the existence of d k * < d k such that l u k d k * = d k * . When d = 0 it is clear that l u k 0 > 0 . Additionally, we know that l u k d k = 0 . Through multiplying (51) by u k d and using (9), (A6) and the Poincaré inequality we obtain
u k d L 2 2 λ a d f u k d , u k d λ m m ε + ε u k d L 2 2 K 1 + 1 2 u k d L 2 2 ,
so, by using the embedding H 0 1 Ω L Ω , l u k d is bounded in d. This implies that the function g d = l u k d has to intersect the line y d = d at some point d k * . It remains to check that d k * < d k . For this aim we prove first that u k d d d k 0 strongly in H 0 1 Ω . Indeed, as u k d is bounded in H 0 1 Ω , there exist v and a sequence { u k d j } such that u k d j v in L 2 Ω . The embedding H 0 1 Ω C [ 0 , 1 ] and the continuity of the function f u imply that { f ( u k d j ) } is bounded in C ( [ 0 , 1 ] ) , so from
u k d j L 2 λ a d j f u k d j L 2 λ m f u k d j L 2 C
we deduce that { u k d j } is bounded in H 2 Ω . Hence, u k d j v in H 0 1 Ω and C 1 ( [ 0 , 1 ] ) . Additionally, f ( u k d j ) f ( v ) in C [ 0 , 1 ] . Therefore, for any ψ H 0 1 Ω we have that
u k d j , ψ = λ a d j f u k d j , ψ v , ψ = λ a d k f v , ψ ,
which implies that v is a solution to problem (51) with d = d k . However, from u k d j v in C 1 ( [ 0 , 1 ] ) it follows that v cannot be a point with less than k + 1 simple zeros in [ 0 , 1 ] and then λ / a d k = k 2 π 2 implies that v 0 . As the limit is the same for every converging subsequence, u k d d d k 0 strongly in H 0 1 Ω . Thus, d k > 0 and lim d d k u k d L 2 = 0 imply that d k * < d k .
Second, let d k = + . Then the existence of d k * < + follows by the same argument as before.
The second and third statements are a consequence of
λ a 0 π 2 k + 1 2 a d π 2 k + 1 2 for any d 0
and
λ inf s 0 { a s } π 2 N + 1 2 a d π 2 N + 1 2 for any d 0 ,
respectively, because in such a case for problem (51) the fixed points v j ± , j > k (respectively j > N ), do not exist.
The last two statements are a consequence of the first two statements and of the fact that the points of intersection of the functions g d = l u k d and y d = d has to be unique, because if a is non-decreasing, then g ( d ) is non-increasing by Lemmas 5 and 6. □
In view of this theorem, we have exactly the same equilibria and bifurcations as in the classical Chafee-Infante equation (see [24,35]) when the function a ( d ) is non-decreasing, because in this case in view of the monotone dependence between the functions a ( d ) and g ( d ) , there is only one intersection point of the function g d with the bisector, as it is shown in Figure 1. This follows from the fact that g ( d ) d is strictly decreasing, but there may be weaker conditions on a · that would lead g ( d ) d to be strictly decreasing.
When the function a · is not assumed to be monotone, an interesting situation appears. More precisely, it is possible to have more than two equilibria with the same number of zeros. If l u = u H 0 1 2 , for the equilibria with k + 1 zeros in [ 0 , 1 ] this happens when the equation
d = 0 1 d u k d ( x ) d x 2 d x = g ( d )
has more than one solution. For instance, if a ( 0 ) = a ( d ¯ ) for some 0 < d ¯ < g ( 0 ) , then g ( 0 ) = g ( d ¯ ) . Assuming that there are 0 < d k 1 < d k 2 < d ¯ such that a ( d k 2 ) = a ( d k 1 ) = λ π 2 k 2 , there must exist 0 < d 1 * < d k 1 < d k 2 < d 2 * < d ¯ such that g ( d i * ) = d i * . Now, by the argument in Theorem 5, there must exist a d 3 * > d ¯ such that g ( d 3 * ) = d 3 * , obtaining six fixed points with k + 1 zeros in [ 0 , 1 ] . This situation is shown in Figure 2, where d 1 * , d 2 * and d 3 * are solutions of (61), that is, there are three intersection points with the bisector. We notice that when a ( d ) > λ / ( π 2 k 2 ) , the function g ( d ) is not defined since the condition for such equilibria to exist is not satisfied, but we can make this function continuous by putting g d = 0 whenever a ( d ) λ / ( π 2 k 2 ) . This procedure establishes that, having fixed a natural number k, for any j N we may construct a · in such a way that we have 2 ( 2 j + 1 ) equilibria with k + 1 zeros in [ 0 , 1 ] .
At least there is always one intersection point with the bisector, but the function g ( d ) could be even tangent to the bisector at some point or not cut it again.

5.3. Lap Number and Some Forbidden Connections

With Theorem 5 at hand we can improve the description of the global attractor given in Theorem 4.
Under conditions (A1)–(A6), (A8) and h 0 , if
a 0 π 2 n 2 < λ a 0 π 2 n + 1 2
then problem (3) possesses exactly 2 n + 1 fixed points: v 0 = 0 , u 1 , d 1 * ± , . . . , u n , d n * ± .
Let ϕ be a bounded complete trajectory. We know by Theorem 4 that
d i s t H 0 1 ( ϕ t , R ) 0 , as t ± .
As the number of fixed points is finite, we will prove that in fact the solution has to converge to one fixed point forwards and backwards. We recall the omega and alpha limit sets of ϕ , given by
ω ϕ = { y : t n + such that ϕ t n y } , α ϕ = { y : t n such that ϕ t n y } ,
are non-empty, compact and connected ([36] Lemma 3.4 and Proposition 4.1). Additionally, d i s t H 0 1 ϕ t , ω ϕ t + 0 , d i s t H 0 1 ϕ t , α ϕ t 0 . Since ω ϕ , α ϕ R and R is finite, the only possibility is that ω ϕ = z 1 R , α ϕ = z 2 R .
Thus, we have established the following result.
Theorem 6.
Let us assume conditions (A1)–(A6), (A8), (62) and h 0 . Then
A = k = 0 2 n + 1 M + v k = k = 0 2 n + 1 M v k ,
where n is given in (62) and v 0 = 0 , v 1 = u 1 , d 1 * + , v 2 = u 1 , d 1 * , . . .
In other words, the global attractor A consists of the set of stationary points R (which has 2 n + 1 elements) and the bounded complete trajectories that connect them (the heteroclinic connections).
Remark 6.
As the Lyapunov function (46) is strictly decreasing along a trajectory ϕ which is not a fixed point, then there cannot exist homoclinic connections for any fixed point. This implies in particular that if n = 0 , then A = { 0 } .
Remark 7.
If we use condition (A7) instead of (A8), then we cannot guarantee that the number of fixed points is finite. However, if we suppose that this is the case, then the result remains valid. In this situation, there could be more than two fixed points with the same number of zeros.
Lemma 7.
Let us assume conditions (A1)–(A6), h 0 and either (A7) or (A8). Let u k , d k * + , u k , d k * be a pair of fixed points corresponding to the same value d k * . Then there cannot be an heteroclinic connection between them.
Proof. 
The function v x = u k , d k * + 1 x is a fixed point corresponding to d k * as
2 v x 2 x = 2 u k , d k * + x 2 1 x = λ a d k * f u k , d k * + 1 x = λ a d k * f v x ,
so u k , d k * x = v x = u k , d k * + 1 x . The equalities
0 1 u k , d k * x x 2 d x = 0 1 u k , d k * + x 1 x 2 d x = 0 1 u k , d k * + x y 2 d y ,
0 1 0 u d k * x f s d s d x = 0 1 0 u d k * + 1 x f s d s d x = 0 1 0 u d k * + y f s d s d y
imply that E ( u k , d k * ) = E u k , d k * + , where E is the Lyapunov function (46). Since this function is strictly decreasing along a trajectory ϕ which is not a fixed point, there cannot exist a heteroclinic connection between these two points. □
Remark 8.
In the case where condition (A7) is assumed, there could be more than two equilibria with k + 1 zeros in [ 0 , 1 ] . In this case there could exist connections between fixed points with different values of the constant d.
By using the concept of lap number of the solutions we can discard some more heteroclinic connections.
We consider the function w t = u ( α 1 t ) , which is a strong solution to problem (47). For any strong solution u · conditions (A1), (A3), (A6) and u C ( [ 0 , + ) , H 0 1 ( Ω ) ) imply that the function
r t , x = λ a ( w t H 0 1 2 ) f w t , x w t , x
is continuous and w · is a solution of the linear equation
w t 2 w x 2 = r t , x w .
Thus, by Theorem A3 in the Appendix A (see also Theorem C in [37]) the number of zeros of w t in [ 0 , 1 ] is a nonincreasing function of t. Since α 1 t is an increasing function of time, the result is also true for the solution u · . Making use of this property we will prove the following result.
Lemma 8.
Let us assume conditions (A1)–(A6), h 0 and either (A7) or (A8). Then if n > k , there cannot exist a connection from the fixed point u k , d k * ± to the fixed point u n , d n * ± , that is, there cannot exist a bounded complete trajectory ϕ such that
ϕ t u n , d n * ± a s t + , ϕ t u k , d k * ± a s t .
Proof. 
By contradiction assume that such complete trajectory exists. Denote by l z the number of zeros of z in [ 0 , 1 ] . By using the compactness of the attractor in C 1 ( [ 0 , 1 ] ) (see Corollary 2) we obtain that
ϕ t u n , d n * ± in C 1 ( [ 0 , 1 ] ) as t + , ϕ t u k , d k * ± in C 1 ( [ 0 , 1 ] ) as t .
Then, as the zeros are simple, we can choose t 1 > 0 large enough such that l ϕ t 1 = l u k , d k * ± = k + 1 . Put u t = ϕ t t 1 , which is a strong solution of (3). Now we choose t 2 > 0 such that l u t 2 = l u n , d n * ± = n + 1 . Then l u 0 = k + 1 and l u ( t 2 ) = n + 1 > k + 1 . This contradicts the fact that the number of zeros of u t is non-increasing. □

6. Morse Decomposition

In this section we study in more detail the structure of the global attactor in the case where the function f is odd. More precisely, we obtain a dynamically gradient m-semiflow G, which is equivalent to saying that there is a Morse decomposition of the attractor [38], and we study the stability of the fixed points.

6.1. Approximations

We consider now the situation when conditions (A1)–(A6), h = 0 and either (A7) or (A8) are satisfied, and moreover, the function f is odd.
In this section we consider the following problems:
u t a ( u H 0 1 2 ) 2 u x 2 = λ f ε n ( u ) , t > 0 , x ( 0 , 1 ) , u ( t , 0 ) = 0 , u ( t , 1 ) = 0 , u ( 0 , x ) = u 0 ( x ) ,
where the function f ε n is defined below and ε n 0 , as n .
Let ρ ε n ( · ) be a mollifier in R . We define the function f ε n ( u ) = R ρ ε n ( s ) f ( u s ) d s . It is well known that f ε n ( · ) C ( R ) and that for any compact subset A R we have f ε n f uniformly on A. It is clear that for u > ε n the function f ε n u is strictly concave.
We need the approximation to fulfil (A2) and (A3). To that end, we consider the approximation except on the interval [ ε n , ε n ] , for any ε n > 0 . There exists a polynomial of sixth degree p ( x ) such that
p ( 0 ) = 0 , p ( ε n ) = h ( ε n ) , p ( 0 ) = 1 , p ( ε n ) = h ( ε n ) , p ( 0 ) = 0 , p ( ε n ) = h ( ε n ) , p ( 0 ) = 1 .
We choose γ > 0 such that p s < 0 for all s ( 0 , γ ] . We can assume that ε n < γ for all n.
Thus, by construction the function
f ε n ( x ) = f ε n ( x ) i f x < ε n , p ( x ) i f ε n x 0 , p ( x ) i f 0 x ε n , f ε n ( x ) i f x > ε n
approximates the function f uniformly in compact sets, that is, for any [ M , M ] and δ > 0 there exists n 0 ( M , δ ) N such that
| f ( x ) f ε n ( x ) | < δ , for all n n 0 , x [ M , M ] .
Additionally, it satisfies the following properties:
(B1)
f ε n C 2 ( R ) ;
(B2)
f ε n ( 0 ) = 0 ;
(B3)
f ε n 0 = 1 ;
(B4)
f ε n is strictly concave if u > 0 and strictly convex if u < 0 ;
(B5)
f ε n is odd.
Lemma 9.
Let f satisfy (A5). Then the functions f ε n satisfy condition (A5) and (9) with independent constants of ε n .
Proof. 
We assume without loss of generality that ε n < 1 . In order to check (4) and (5) we only need to consider u outside the interval [ 1 , 1 ] , because the sequence { f ε n } is uniformly bounded in any compact set of R . Then for u [ 1 , 1 ] the Hölder inequality and R ρ ε n s d s = 1 give
f ε n u = R f u s ρ ε n s d s R f u s ρ ε n s d s R C 1 + C 2 | u s | p 1 ρ ε n s d s C 1 + C 2 2 p 2 ε n ε n u p 1 + s p 1 ρ ε n s d s C ˜ 1 + C ˜ 2 u p 1 .
If f satisfies (5), then
f ε n u u = R f u s u s ρ ε n s d s + R f u s s ρ ε n s d s R C 3 C 4 | u s | p ρ ε n s d s + R C 1 + C 2 | u s | p 1 s ρ ε n s d s K 1 C 4 R 2 1 p u p s p ρ ε n s d s + C 2 2 p 2 R u p 1 + s p 1 s ρ ε n s d s C ˜ 3 C ˜ 4 | u | p ,
where we have used u p 2 p 1 s p + u s p and the Young inequality.
For (9) we put in the above inequality p = 2 , C 3 = m ε , C 4 = ε and obtain
f ε n u u m ˜ ε + ε u 2 ,
which obviously implies (6). □
Our next aim is to focus on the convergence of solutions of the approximations.
Theorem 7.
Let conditions (A1)–(A6), h = 0 and either (A7) or (A8) be satisfied and let, moreover, the function f be odd. If u ε n , 0 u 0 in H 0 1 ( Ω ) as ε n 0 , then for any sequence of solutions of (64) u ε n ( · ) with u ε n ( 0 ) = u ε n , 0 there exists a subsequence of ε n such that u ε n converges to some strong solution u · of (3) in the space C ( [ 0 , T ] , H 0 1 ( Ω ) ) , for any T > 0 .
Proof. 
By using (29) and (30) we can repeat the same lines of the proof of Theorem 1 and obtain the existence of a function u · and a subsequence of u ε n such that
u ε n * u in L ( 0 , T ; H 0 1 ( Ω ) ) ,
u ε n u in L 2 ( 0 , T ; D ( A ) ) ,
d u ε n d t d u d t in L 2 ( 0 , T ; L 2 ( Ω ) ) ,
u ε n u in C ( [ 0 , T ] ; L 2 ( Ω ) ) ,
u ε n u in L 2 ( 0 , T ; H 0 1 ( Ω ) ) ,
f ε n ( u n ε ) * f ( u ) in L ( 0 , T ; L ( Ω ) ) ,
a ( u ε n H 0 1 2 ) Δ u ε n a ( u H 0 1 2 ) Δ u in L 2 ( 0 , T ; L 2 ( Ω ) ) .
Additionally, in the same way we prove that u · is a strong solution to problem (3) such that u 0 = u 0 .
The uniform estimate in the space H 0 1 Ω implies also that if t n t 0 , then u ε n t n u t 0 in H 0 1 ( Ω ) . We need to prove that this convergence is in fact strong, proving then the convergence in C ( [ 0 , T ] , H 0 1 ( Ω ) ) for any T > 0 .
In the same way as in the proof of Lemma 1 we deduce that for some C > 0 the functions Q n ( t ) = A ( u ε n ( t ) H 0 1 2 ) 2 C t , Q ( t ) = A ( u ( t ) H 0 1 2 ) 2 C t are continuous and non-increasing in [ 0 , T ] . Moreover, Q n ( t ) Q ( t ) for a.e. t ( 0 , T ) . Let first t 0 > 0 ˙ . Then we take 0 < t j < t 0 such that t j t 0 and Q n ( t j ) Q ( t j ) for all j . Then
Q n ( t n ) Q ( t 0 ) Q n ( t j ) Q ( t 0 ) | Q n ( t j ) Q ( t j ) | + | Q ( t j ) Q ( t 0 ) | for t n > t j .
For any δ > 0 there exist j ( δ ) and N ( j ( δ ) ) such that Q n ( t n ) Q ( t 0 ) δ if n N , so lim sup Q n ( t n ) Q ( t 0 ) . Hence, a contradiction argument using the continuity of A s shows that lim sup u ε n ( t n ) H 0 1 2 u ( t 0 ) H 0 1 2 . This, together with lim inf u ε n ( t n ) H 0 1 2 u ( t 0 ) H 0 1 2 , implies that u ε n ( t n ) H 0 1 2 u ( t 0 ) H 0 1 2 , so that u ε n ( t n ) u ( t 0 ) strongly in H 0 1 ( Ω ) . For the case when t 0 = 0 we use the same argument as in Lemma 1. □
We denote by A ε n the global attractor for the semiflow G ε n corresponding to problem (64).
Lemma 10.
Assume the condition of Theorem 7. Then n N A ε n is bounded in H 0 1 ( Ω ) . Hence, the set n N A ε n ¯ is compact in L 2 ( Ω ) .
Proof. 
By Lemma 9 inequality (9) is satisfied for any n with constants which are independent of ε n , so inequality (36) holds true with constants independent of ε n . Thus, there a exists a common absorbing ball B 0 in L 2 Ω (with radius K > 0 ) for problems (64). Further, by repeating the same steps as in Proposition 3 we obtain a common absorbing ball in H 0 1 Ω (with radius K ˜ > 0 ), as by Lemma 9 the constants which are involved are independent of ε n . Thus, y H 0 1 K ˜ for any y n N A ε n .
Lemma 11.
Assume the condition of Theorem 7. Then n N A ε n is bounded in V 2 r for any 0 r < 1 . Hence, n N A ε n ¯ is compact in V 2 r and C 1 ( [ 0 , 1 ] ) .
Proof. 
By using Lemma 10 we obtain the boundedness of n N A ε n in V 2 r by repeating the same lines in Lemma 4. The rest of the proof follows from the compact embedding V α V β , α > β , and the continuous embedding V 2 r C 1 ( [ 0 , 1 ] ) if r > 3 4 .
Corollary 3.
Assume the condition of Theorem 7. Then any sequence ξ n A ε n with ε n 0 is relatively compact in C 1 ( [ 0 , 1 ] ) .
Lemma 12.
Assume the condition of Theorem 7. Then up to a subsequence any bounded complete trajectory u ε n of (64) converges to a bounded complete trajectory u of (3) in C ( [ T , T ] , H 0 1 ( Ω ) ) for any T > 0 . On top of that, if y n A ε n , then passing to a subsequence y n y A in H 0 1 Ω . Hence,
d i s t H 0 1 A ε n , A 0 a s n .
Proof. 
Let us fix T > 0 . By Corollary 3 u ε n ( T ) y in H 0 1 ( Ω ) up to a subsequence. Theorem 7 implies that u ε n converges in C ( [ T , T ] , H 0 1 ( Ω ) ) to some solution u of (3). If we choose successive subsequences for 2 T , 3 T and apply the standard diagonal procedure, we obtain that a subsequence u ε n converges to a complete trajectory u of (3) in C ( [ T , T ] , H 0 1 ( Ω ) ) for any T > 0 . Finally, from Lemma 10 this trajectory is bounded.
If y n A ε n , by Corollary 3 we can extract a subsequence converging to some y. If we take a sequence of bounded complete trajectories ϕ n · of (64) such that ϕ n 0 = y n , then by the previous result it converges in C ( [ T , T ] , H 0 1 ( Ω ) ) to some bounded complete trajectory ϕ · of (3), so y A .
Finally, if (67) was not true, there would exist δ > 0 and a sequence y n A ε n such that d i s t H 0 1 ( y , A ) > δ . However, passing to a subsequence y n y A , which is a contradiction. □
Lemma 13.
Assume the conditions of Theorem 7. Let τ ± d n , ε n be the functions (56)–(57) for problem (51) but replacing f by f ε n and d by d n . Let d n , E n 0 as n . Then
lim n τ ± d n , ε n ( E n ) = a 0 π 2 λ .
Proof. 
Let us consider f d n , ε n ( u ) = λ f ε n ( u ) a ( d n ) . In view of property ( B 4 ) and (66), since f ε n ( 0 ) = f ( 0 ) = 1 and f ε n ( 0 ) = f ( 0 ) = 0 , given γ ( 0 , 1 ) there exists δ > 0 (independent of ε n ) such that
( 1 γ ) u f ε n ( u ) ( 1 + γ ) u , for any u ( 0 , δ ) . 1 1 + γ u f ε n ( u ) 1 1 γ , for any u ( 0 , δ ) .
The sequence F ε n · converges uniformly to F · in compact sets. Moreover, as U + ( E ) is continuous and using ([39] p. 60), given δ > 0 , there exists η > 0 such that U + ε n ( E ) δ for any 0 < E η . Now, if we integrate the first inequality in (68) between 0 and u we obtain
1 2 ( 1 γ ) u 2 F ε n ( u ) 1 2 ( 1 + γ ) u 2 , for any 0 u δ .
By using the change of variable E n y 2 = F ε n ( u ) , we have
1 γ 2 E n 1 / 2 u y 1 + γ 2 E n 1 / 2 u , if 0 < E n η , 0 y 1 .
Dividing the previous expression by λ a d n f d n , ε n ( u ) and using (68) we obtain
a ( d n ) ( 1 γ ) 2 λ E n ( 1 + γ ) 2 1 / 2 a d n y λ f d n , ε n ( u ) a ( d n ) ( 1 + γ ) 2 λ E n ( 1 γ ) 2 1 / 2 , if 0 < E n η , 0 y 1 .
Now if we multiply by 2 E n ( 1 y 2 ) 1 2 and integrate from 0 to 1, we get
π a ( d n ) ( 1 γ ) 2 λ ( 1 + γ ) 2 1 / 2 τ + ε n ( E n ) π a ( d n ) ( 1 + γ ) 2 λ ( 1 γ ) 2 1 / 2 , if 0 < E n η .
Then the theorem follows as a d n a 0 when n . The proof for τ ε n is analogous. □
Under the conditions of Theorem 7, if (A8) is satisfied and
a 0 π 2 k 2 < λ a 0 π 2 k + 1 2 , k Z , k 0 ,
holds, then by Theorem 5 problem (64) has exactly 2 k + 1 fixed points (denoted by v 0 = 0 , v 1 , d 1 ε n ± , . . . , v k , d k ε n ± ) and v m , d m ε n ± has m + 1 zeros in [ 0 , 1 ] for each 1 m k . The same is valid for problem (3) and we denote the 2 k + 1 fixed points by v 0 = 0 , u 1 , d 1 * ± , . . . , u k , d k * ± .
Lemma 14.
Assume the conditions of Theorem 7, (A8) and (69). Let m N , 1 m k , be fixed. Then v m , d m ε n + (resp. v m , d m ε n ) do not converge to 0 in H 0 1 ( Ω ) as ε n 0 .
Proof. 
Assume that v m , d m ε n + 0 in H 0 1 0 , 1 . Then it converges to 0 in C [ 0 , 1 ] and the equality
d 2 v m , d m ε n + d x 2 x = λ f ε n v m , d m ε n + x a ( d m ε n )
implies that v m , d m ε n + 0 in C 2 [ 0 , 1 ] . In particular, d v m , d m ε n + d x 0 0 and d m ε n = v m , d m ε n + H 0 1 2 0 . The value E n corresponding to the fixed point v m , d m ε n + is equal to a d m ε n 2 λ d v m , d m ε n + d x 0 , so E n 0 . We will show that this is not possible. We know by Lemma 13 that
lim n τ ± d m ε n , ε n E n = π a ( 0 ) 2 λ .
Additionally, since v m , d m ε n + is a fixed point with d = d m ε n one of the following conditions has to be satisfied (see (55)):
j τ + d m ε n , ε n E n + j 1 τ d m ε n , ε n E n = 1 2 1 2 ,
j τ d m ε n , ε n E n + j 1 τ + d m ε n , ε n E n = 1 2 1 2 , if m = 2 j 1
j τ + d m ε n , ε n ( E n ) + j τ d m ε n , ε n ( E n ) = 1 2 1 2 , if m = 2 j .
Since E n 0 and λ > k 2 π 2 a ( 0 ) m 2 π 2 a ( 0 ) , there exists ε n 0 such that
τ ± d m ε n 0 , ε n 0 ( E n 0 ) < 1 2 m .
Hence, neither of (70)–(72) is possible. □
Lemma 15.
Assume the conditions of Theorem 7, (A8) and (69). Let m N , 1 m k , be fixed. Then v m , d m ε n + (resp. v m , d m ε n ) converges to u m , d m * + in H 0 1 ( Ω ) (resp. u m , d m * ) as ε n 0 .
Proof. 
We consider v m , d m ε n + . In view of Corollary 3, v m , d m ε n + is relatively compact in C 1 [ 0 , 1 ] , so up to a subsequence v m , d m ε n + v strongly in C 1 ( [ 0 , 1 ] ) and d m ε n d * = v H 0 1 2 . The proof will be finished if we prove that v = u m , d m * + . We observe that since in such a case every subsequence would have the same limit, the whole sequence would converge to u m , d m * + .
In view of (66) f ε n ( v m , d m ε n + ) converges to f v in C ( [ 0 , 1 ] ) . It follows that
2 v x 2 = λ f v a ( v H 0 1 2 )
and v is a solution of (50), so v is a fixed point of (3). We need to prove that v = u m , d m * + . By Lemma 14 v 0 , and then v = u j , d j * ± for some 1 j k . Since u j , d j * ± has j + 1 simple zeros, the convergence v m , d m ε n + u j , d j * ± in C 1 ( [ 0 , 1 ] ) implies that v m , d m ε n + has j + 1 zeros for n N . However, v m , d m ε n + possesses m + 1 zeros in [ 0 , 1 ] . Thus, m = j .
For the sequence v m , d m ε n the proof is analogous. □

6.2. Instability

We will prove that the fixed points 0 and u k , d k * ± , k 2 , are unstable under some additional assumptions on the functions f and a. For that aim we need to use the approximative problems (64).
Theorem 8.
Assume that the conditions (A1)–(A8), h = 0 and (69) with k 1 are satisfied; and let the function f · be odd and a · be globally Lipschitz continuous. Then the equilibria v 0 = 0 and u j , d j * ± , 2 j k (if k 2 ), are unstable.
Remark 9.
The condition that a · is globally Lipschitz continuous could be dropped, as we can replace a · in (64) by a sequence a ε n · of globally Lipschitz continuous functions.
Proof. 
Problem (64) generates a single-valued semigroup { T ε n ( t ) ; t 0 } with a finite number of fixed points: v 0 = 0 , v 1 , d 1 ε n ± , . . . , v k , d k ε n ± [26]. We know by Theorems 3.5 and 3.6 in [26] that for any v j , d j ε n + with j 2 and v 0 there exists a bounded complete trajectory u ε n such that
u ε n ( t ) v j , d j ε n + as t , for k 2 ,
so v 0 , v j , d j ε n + are unstable. The same is valid for v j , d j ε n . On the other hand, by Lemma 15 we have
v j , d j ε n ± u j , d j * ± ,
where u j , d j * ± is a fixed point of problem (3) with j + 1 zeros in [ 0 , 1 ] . We prove the result for u j , d j * + . For u j , d j * and v 0 the proof is the same.
By Lemma 12 we obtain that up to a subsequence u ε n converges to a bounded complete trajectory u of problem (3) in the space C ( [ T , T ] , H 0 1 ( Ω ) ) for every T > 0 . Thus, either u · is a fixed point v 1 or by Theorem 4 there exists a fixed point v 1 of problem (3) such that
u ( t ) v 1 as t in H 0 1 ( Ω ) .
In the second case, if v 1 = u j , d j * + , the proof would be finished, so let assume the opposite.
Assume first that either u · is not a fixed point or it is a fixed point but v 1 u j , d j * + . We consider r 0 > 0 such that the neighbourhood O 2 r 0 ( v 1 ) does not contain any other fixed point of problem (3). For any r r 0 we can choose t r and n r such that u ε n ( t r ) O r ( v 1 ) for all n n r . On the other hand, since u ε n ( t ) v j , d j ε n + , as t , and v j , d j ε n + u j , d j * + B 2 r 0 ( v 1 ) , there exists t r < t r such that
u ε n r ( t ) O r 0 ( v 1 ) for t ( t r , t r ] ,
u ε n r ( t r ) v 1 H 0 1 = r 0 .
Let first t t t r + . We define the sequence u 1 ε n r ( t ) = u ε n r ( t + t r ) , which passing to a subsequence converges to a bounded complete trajectory ϕ t such that ϕ t O r 0 ( v 1 ) for all t 0 . As there is no other fixed point in O 2 r 0 ( v 1 ) , ϕ t v 1 as t + . However, ϕ 0 v 1 = r 0 , so ϕ · is not a fixed point. Then ϕ t v 2 as t , where v 2 is a fixed point different from v 1 . Second, let t t t r C . Then put u 1 ε n r ( t ) = u ε n r ( t + t r ) . Passing to a subsequence we have that
u 1 ε n r ( 0 ) v 1 , t r t r t 0 , as r 0 .
Additionally, u 1 ε n r · converges to a bounded complete trajectory u 1 · of problem (3) such that u 1 ( 0 ) = v 1 . Let
ψ 1 ( t ) = u 1 t if t 0 , v 1 if t 0 .
We note that u 1 ( t 0 ) v 1 H 0 1 = r 0 implies that u 1 · is not a fixed point. Then ψ 1 is a bounded complete trajectory of problem (3) such that ψ 1 ( t ) v 2 v 1 as t . If v 2 = u j , d j * + , the proof is finished.
If v 2 u j , d j * + , we continue constructing by the same procedure a chain of connections in which the new fixed point is always different from the previous ones, because the existence of the Lyapunov function (46) avoids the existence of a cyclic chain of connections. Since the number of fixed points is finite, at some moment we obtain a bounded complete trajectory ϕ · such that ϕ t u j , d j * + as t , proving that u j , d j * + is unstable.
Now let u · = v 1 = u j , d j * + . Defining the neighbourhood O 2 r 0 ( v 1 ) as before, for any r r 0 we can choose n r such that u ε n ( 0 ) O r ( v 1 ) for all n n r . Additionally, since u ε n ( t ) z 0 n , as t + , where z 0 n v j , d j ε n + is a fixed point of (64), there exists t r > 0 such that
u ε n r ( t ) O r 0 ( v 1 ) for t [ 0 , t r ) ,
u ε n r ( t r ) v 1 H 0 1 = r 0 .
The sequence { t r } cannot be bounded. Indeed, if t r t 0 , then u ε n r ( t r ) u t 0 = v 1 , which is a contradiction with u ε n r ( t 0 ) v 1 H 0 1 = r 0 . Then t r + . We define the functions u 1 ε n r ( t ) = u ε n r ( t + t r ) , which satisfy that u 1 ε n r ( t ) O r 0 ( v 1 ) for all t [ t r , 0 ) . Passing to a subsequence it converges to a bounded complete trajectory ϕ · such that ϕ t O r 0 ( v 1 ) for all t 0 . This trajectory is not a fixed point as ϕ ( 0 ) v 1 H 0 1 = r 0 and ϕ t u j , d j * + as t , so u j , d j * + is unstable. □
Further, we will prove that there is also a connection from 0 to the point u k , d k * ± .
Theorem 9.
Assume the conditions of Theorem 8. Then there exists a bounded complete trajectory ϕ · such that ϕ t t 0 , ϕ t t + u k , d k * + (and the same is valid for u k , d k * ). Thus, E ( 0 ) = 0 > E ( u k , d k * ± ) .
Proof. 
We start with the case where k = 1 . We have three fixed points: 0 , u 1 , d 1 * + , u 1 , d 1 * . By Theorem 8 there exists a bounded complete trajectory ϕ · such that ϕ t t 0 , whereas Theorem 4 and Remark 6 imply that it has to converge forward to a fixed point different from 0, that is, to either u 1 , d 1 * + or u 1 , d 1 * . If, for example, ϕ t t + u 1 , d 1 * + , then as the function f is odd, ψ t = ϕ t is another bounded complete trajectory and ψ t t + u 1 , d 1 * + = u 1 , d 1 * .
Further we consider the problem
u t a ( u H 0 1 2 ) 2 u x 2 = λ f k ( u ) , t > 0 , 0 < x < 1 k , u ( t , 0 ) = u ( t , 1 k ) = 0 , u ( 0 , x ) = u 0 ( x ) ,
where f k ( u ) = k f u / k satisfies (A1)–(A5). In this problem, condition (69) implies that there are again three fixed points: 0 , u 1 , d 1 * , 1 k + , u 1 , d 1 * , 1 k . By the above argument there is a connection ϕ 1 k · from 0 to u 1 , d 1 * , 1 k + (also to u 1 , d 1 * , 1 k ). Since the function f is odd, u k , d k * + x is equal to 1 k u 1 , d 1 * , 1 k + x on [ 0 , 1 k ] , to 1 k u 1 , d 1 * , 1 k + x 1 k on [ 1 k , 2 k ] , etc. Then the function ϕ · such that ϕ t , x = ( 1 ) j k ϕ 1 k t , x j k on [ j k , j + 1 k ] , j = 0 , 1 , . . . , k 1 , is a bounded complete trajectory of problem (3) which goes from 0 to u k , d k * + . □
Remark 10.
When k = 1 the structure of the global attractor is the same as in the Chafee-Infante equation.

6.3. Gradient Structure

We will obtain that the m-semiflow G is dynamically gradient. Let us recall this concept.
A weakly invariant set M of X is isolated if there is a neighbourhood O of M such that M is the maximal weakly invariant subset on O . If M belongs to the global attractor A , then it is compact ([38] Lemma 19). In this case, it is equivalent to use a δ -neighbourhood O δ ( M ) = { y X : d i s t y , M < δ } .
Suppose that there is a finite disjoint family of isolated weakly invariant sets M = { M 1 , , M m } in A , that is, for every j { 1 , , n } there is ϵ j > 0 such that M j A is the maximal weakly invariant set on O ϵ j ( M j ) , and suppose that there exists δ > 0 such that O δ ( M i ) O δ ( M j ) = , if i j .
Definition 4.
We say the m-semiflow G : R + × X P ( X ) is dynamically gradient with respect to the disjoint family of isolated weakly invariant sets M = { M 1 , , M m } if for every complete and bounded trajectory ψ of R we have that either ψ ( R ) M j , for some j { 1 , , m } , or α ( ψ ) M i and ω ( ψ ) M j with 1 j < i m .
Let us consider the case when the conditions of Theorem 8 hold. Then (3) possesses exactly 2 k + 1 fixed points: v 0 = 0 , u 1 , d 1 * ± , . . . , u k , d k * ± . Additionally, as f is odd, u j , d j * + = u j , d j * for any j. We define the following sets:
M 1 = { u 1 , d 1 * + , u 1 , d 1 * } , . . . , M k = { u k , d k * + , u k , d k * } , M k + 1 = { 0 } .
They are weakly invariant and using Lemma 7 we deduce easily that they are isolated. Then the family M = { M 1 , , M k + 1 } is a finite disjoint family of isolated weakly invariant sets.
Proposition 4.
Assume the conditions of Theorem 8. Then G is a dynamically gradient semiflow with respect to the family (75) after (possibly) reordering them.
Proof. 
We reorder the family (75) in such a way that if the value of the Lyapunov function E given in (46) is equal to L i for the set M ˜ i , then L j L n for j < n . Then Theorem 25 in [38] implies that G is dynamically gradient with respect to this family. □
We will obtain then that the fixed points u 1 , d 1 * + , u 1 , d 1 * are asymptotically stable. The compact set M A is a local attractor for G in X if there is ε > 0 such that ω O ε ( M ) = M , where
ω B = { y : t n + , y n G ( t n , B ) such that y n y }
is the ω -limit set of B. By Lemma 14 in [38] if M is a local attractor in X, then it is stable. Thus, a local attractor is asymptotically stable.
Theorem 10.
Assume the conditions of Theorem 8. Then the stationary points u 1 , d 1 * + , u 1 , d 1 * are asymptotically stable.
Proof. 
By ([38] Theorem 23 and Lemma 15) M ˜ 1 is a local attractor in X, so it is asymptotically stable. By Theorem 8 the sets M j , j 2 , are unstable. Thus, M ˜ 1 = M 1 . As M 1 consists of the two elements u 1 , d 1 * + , u 1 , d 1 * , which are obviously disjoint, they are asymptotically stable as well. □
We will prove that there is a connection from 0 to any other fixed point u j , d j * ± .
Theorem 11.
Assume the conditions of Theorem 8. Then there exists a bounded complete trajectory ϕ · such that ϕ t t 0 , ϕ t t + u j , d j * + for all 1 j k (and the same is valid for u j , d j * ).
Proof. 
Let us consider problem (74) with k = j . The function u 1 , d j * , 1 j + ( x ) = j u j , d j * + ( x ) , x [ 0 , 1 j ] , is the unique positive fixed point of problem (74). Let X j + = { u H 0 1 0 , 1 j : u x 0 x [ 0 , 1 j ] } be the positive cone of H 0 1 0 , 1 j . If we consider the restriction of the semigroup T j ε n · of problem (64) in the interval 0 , 1 j to X j + , denoted by T j ε n , + · , then there exists a global attractor A n , j + [25]. Since 0 and v 1 , d j ε n , 1 j + = j v j , d j ε n + [ 0 , 1 j ] are the unique fixed points of T j ε n , + , A n , j + is connected, v 1 , d 1 ε n , 1 j + is stable [26] and A n , j + consists of the fixed points and their heteroclinic connections, there must exist a bounded complete trajectory ϕ j ε n · of T j ε n , + which goes from 0 to v 1 , d j ε n , 1 j + . By Lemma 12 up to a subsequence it converges to a bounded complete trajectory ϕ j · of problem (74) with k = j such that ϕ j t 0 for all t R . Since by Theorem 10 the fixed point u 1 , d j * , 1 j + is stable, the only possibility is that ϕ j t 0 , as t , ϕ j t u 1 , d j * , 1 j + , as t + . Then the function ϕ · such that ϕ t , x = 1 i j ϕ j t , x i j on [ i j , i + 1 j ] , i = 0 , 1 , . . . , j 1 , is a bounded complete trajectory of problem (3) which goes from 0 to u j , d j * + .
For u j , d j * , noting that u j , d j * = u j , d j * + , the result follows by choosing the bounded complete trajectory ϕ ˜ ( t ) = ϕ t . □
As a consequence we obtain that the order of the family M has to be the one given in (75).
Theorem 12.
The semiflow G is dynamically gradient with respect to the family M in the order given in (75), that is, M ˜ i = M i for any i.
Proof. 
As by Theorem 11 there is a connection from 0 to u j , d j * ± , 1 j k , we have proved that M ˜ k + 1 = { 0 } = M k + 1 . The fact that the order of the other sets is the one given in (75) follows from Lemma 8. □

Author Contributions

Conceptualization, R.C., A.N.C., P.M.-R. and J.V.; methodology, R.C., A.N.C., P.M.-R. and J.V.; validation, R.C., A.N.C., P.M.-R. and J.V.; formal analysis, R.C., A.N.C., P.M.-R. and J.V.; investigation, R.C., A.N.C., P.M.-R. and J.V.; writing—original draft preparation, R.C., A.N.C., P.M.-R. and J.V.; writing—review and editing, R.C., A.N.C., P.M.-R. and J.V.; supervision, R.C., A.N.C., P.M.-R. and J.V.; funding acquisition, R.C., A.N.C., P.M.-R. and J.V. All authors have read and agreed to the published version of the manuscript.

Funding

The first author is a fellow of the FPU program of the Spanish Ministry of Education, Culture and Sport, reference FPU15/03080. This work has been partially supported by the Spanish Ministry of Science, Innovation and Universities, project PGC2018-096540-B-I00; by the Spanish Ministry of Science and Innovation, project PID2019-108654GB-I00; and by the Junta de Andalucía and FEDER, project P18-FR-4509.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank also the referees for their useful comments. During the preparation of this manuscript our colleague and friend María José Garrido-Atienza passed away. She was a very kind and warm person and we will miss her a lot. We dedicate this paper to her memory. We would like also to devote this work to the memory of Encarna López, Rubén’s grandmother, with sorrow and love. She was always an example of goodness and willpower. Even in her last moments, she always had a smile in bad times. The whole family will miss her very much.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this appendix we generalise the lap number property of solutions of linear equations proved in [24] to the case when we do not have classical solutions. For this we will use a maximum principle for non-smooth functions from [40].
Let O be a region in R 2 and let t 0 , x 0 O and ρ , σ > 0 . We denote
Q ρ , σ = { ( t , x ) : t ( t 0 σ , t 0 ) , x x 0 < ρ } ,
where we assume that t 0 , x 0 , ρ , σ are such that Q ¯ ρ , σ O .
We denote by W the space of all functions from L 2 O such that
O u t , x 2 + u x t , x 2 d μ < + .
As a particular case of Theorem 6.4 in [40] we obtain the following maximum and minimum principles.
Theorem A1.
(Maximum principle) Let u W be such that
u t 2 u x 2 0
in the sense of distributions. If
sup e s s t , x Q ρ ν , σ 1 u ( t , x ) = M ,
for some ν, 0 < ν < 1 , and any σ 1 , where 0 < σ 1 < σ , then u t , x = M for a.a. t , x Q ρ , σ .
Theorem A2.
(Minimum principle) Let u W be such that
u t 2 u x 2 0
in the sense of distributions. If
inf e s s t , x Q ρ ν , σ 1 u ( t , x ) = M ,
for some ν, 0 < ν < 1 , and any σ 1 , where 0 < σ 1 < σ , then u t , x = M for a.a. t , x Q ρ , σ .
We are ready to prove the lap-number property, saying that the number of zeros is a non-increasing function of time.
Theorem A3.
Let r t , x be a continuous function and u C ( [ t 0 , t 1 ] , H 0 1 Ω ) L 2 t 0 , t 1 ; H 2 Ω be such that d u d t L 2 t 0 , t 1 ; L 2 Ω and satisfies the equation
u t 2 u x 2 = r ( t , x ) u , 0 < x < 1 , t 0 < t t 1 .
Then the number of components of
{ x : 0 < x < 1 , u t , x 0 }
is a non-increasing function of t.
Proof. 
We follow similar lines as in ([24] Theorem 6).
Denote Q t = { x 0 , 1 : u t , x 0 } . We need to show that there is an injective map from the components of Q t 1 to the components of Q t 0 if t 1 > t 0 . If we denote by C a component of Q t 1 and by S C the component of [ t 0 , t 1 ] × 0 , 1 { u t , x ) 0 } which contains C, then in order to obtain the injective map it is necessary to prove two facts:
1.
S C Q ( t 0 ) Ø ;
2.
If C 1 , C 2 are two components of Q t 1 , then S C 1 S C 2 = Ø .
Let us prove the first statement by contradiction, so assume that S C Q ( t 0 ) = Ø . We can assume without loss of generality that r t , x < 0 , because this property is satisfied for the function W t , x = u t , x e λ t with λ > 0 large enough and the components of these two functions coincide. Consider, for example, that u t , x > 0 in S C . Let M = max t , x S C u t , x . By hypothesis and the Dirichlet boundary conditions this maximum has to be attained at a point t , x such that t 0 < t t 1 , 0 < x < 1 . Additionally, there has to be an ε > 0 such that if t , x S C and t 0 < t t 0 + ε , then u t , x < M , as otherwise there would be a sequence t n , x n S C , t n > t 0 , such that t n t 0 and u t n , x n = M . By the continuity of u this would imply that u t 0 , x 0 = M for some t 0 , x 0 S C , which is a contradiction. Then we can choose t as the first time when the maximum is attained, so u t , x < M for all t , x S C , t 0 < t < t . By the continuity of u there exists a rectangle R = [ t δ , t ] × [ x γ , x + γ ] such that R belongs to S C . In order to apply Theorem A1 we put O = R and
Q γ , δ = { ( t , x ) : t ( t δ , t ) , x x < γ } .
We have that
sup ( t , x ) Q ν γ , σ 1 u ( t , x ) = M ,
for some 0 < ν < 1 and any 0 < σ 1 < δ . Since u satisfies (A1), we conclude from Theorem A1 that u t , x = M for all t , x Q ρ , σ , which is a contradiction.
For the second statement suppose the existence of two disjoints components C 1 , C 2 of Q t 1 such that S C 1 S C 2 Ø , which implies in fact that S C 1 = S C 2 . In this case we can assume that r t , x > 0 , being this justified by the function W t , x = u t , x e λ t with λ > 0 large enough. Let, for example, u t , x > 0 in S C 1 and assume that the interval C 1 has lesser values than the interval C 2 . Additionally, it is clear that between C 1 and C 2 there must exist a point t 1 , x 0 such that u t 1 , x 0 = 0 . On the other hand, the set S C 1 t 0 , t 1 × [ 0 , 1 ] is path connected. Thus, there exists a simple path ξ such that one end point is in { t 1 } × C 1 and the other one is in { t 1 } × C 2 . Let us consider the set L of all points which are above the curve ξ and such that the function u vanishes at them. This set is non-empty because t 1 , x 0 L . Since L is compact, the function g : L [ t 0 , t 1 ] given by g t , x = t attains it minimum at a certain point t , x L such that t 0 < t . Then there exists a set R = [ t δ , t ) × [ x γ , x + γ ] which belongs to S C 1 . Let O = R and
Q γ , δ = { ( t , x ) : t ( t δ , t ) , x x < γ } .
We have that
inf ( t , x ) Q ν γ , σ 1 u ( t , x ) = 0 ,
for some 0 < ν < 1 and any 0 < σ 1 < δ . Since u satisfies (A2), we conclude from Theorem A2 that u t , x = 0 for all t , x Q ρ , σ , which is a contradiction. □

References

  1. Burger, M.; Caffarelli, L.; Markowich, P.A. Partial differential equation models in the socio-economic sciences. Phil. Trans. R. Soc. A 2014, 372, 2028. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ramsey, F. A mathematical theory of saving. Econ. J. 1928, 38, 543–559. [Google Scholar] [CrossRef]
  3. Brito, P. A Bentham-Ramsey model for spatially heterogeneous growth. In Working Papers of the Department of Economics; ISEG, University of Lisboa: Lisbon, Portugal, 2001. [Google Scholar]
  4. Brito, P. The dynamics of growth and distribution in a spatially heterogeneous world. In Working Papers of the Department of Economics; ISEG, University of Lisboa: Lisbon, Portugal, 2004. [Google Scholar]
  5. Brito, P. Global endogeneous growth and distributional dynamics. In Proceedings of the Dynamics, Economic Growth, and International Trade, Milan, Italy, 13–14 September 2012. [Google Scholar]
  6. Anita, S.; Capasso, V.; Kunze, H.; La Torre, D. Dynamics and optimal control in a spatially structured economy growth model with pollution diffusion and environmental taxation. Appl. Math. Lett. 2015, 42, 36–40. [Google Scholar] [CrossRef]
  7. Deng, K.; Wu, Y. Global stability for a nonlocal reaction-diffusion population model. Nonlinear Anal. Real World Appl. 2015, 25, 127–136. [Google Scholar] [CrossRef]
  8. Achleitner, F.; Kuehn, C. On bounded positive stationary solutions for a nonlocal Fisher-KPP equation. Nonlinear Anal. 2015, 112, 15–29. [Google Scholar] [CrossRef] [Green Version]
  9. Delgado, M.; Molina-Becerra, M.; Santos Júnior, J.R.; Suárez, A. A non-local perturbation of the logistic equation in ℝN. Nonlinear Anal. 2019, 187, 147–158. [Google Scholar] [CrossRef]
  10. Chipot, M.; Lovat, B. Some remarks on nonlocal elliptic and parabolic problems. Nonlinear Anal. 1997, 30, 461–627. [Google Scholar] [CrossRef]
  11. Chipot, M.; Lovat, B. On the asymptotic behaviour of some nonlocal problems. Positivity 1999, 3, 65–81. [Google Scholar] [CrossRef]
  12. Chipot, M.; Molinet, L. Asymptotic behaviour of some nonlocal diffusion problems. Appl. Anal. 2001, 80, 273–315. [Google Scholar]
  13. Chipot, M.; Rodrigues, J.F. On a class of nonlocal nonlinear elliptic problems. Math. Model. Numer. Anal. 1992, 26, 447–467. [Google Scholar] [CrossRef]
  14. Chipot, M.; Siegwart, M. On the Asymptotic behaviour of some nonlocal mixed boundary value problems. In Nonlinear Analysis and Applications: To V. Lakshmikantham on his 80th Birthday; Agarwal, R.P., O’Regan, D., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2003; Volume 1, pp. 431–449. [Google Scholar]
  15. Chipot, M.; Valente, V.; Vergara Caffarelli, G. Remarks on a nonlocal problem involving the Dirichlet energy. Rend.Sem. Mat. Univ. Padova 2003, 110, 199–220. [Google Scholar]
  16. Zheng, S.; Chipot, M. Asymptotic behavior of solutions to nonlinear parabolic equations with nonlocal terms. Asymptot. Anal. 2005, 45, 301–312. [Google Scholar]
  17. Anh, C.T.; Tinh, L.T.; Toi, V.M. Global attractors for nonlocal parabolic equations with a new class of nonlinearities. J. Korean Math. Soc. 2018, 55, 531–551. [Google Scholar]
  18. Caballero, R.; Marín-Rubio, P.; Valero, J. Existence and characterization of attractors for a nonlocal reaction-diffusion equation with an energy functional. J. Dyn. Differ. Equ. 2021, in press. [Google Scholar] [CrossRef]
  19. Caraballo, T.; Herrera-Cobos, M.; Marín-Rubio, P. Long-time behavior of a non-autonomous parabolic equation with nonlocal diffusion and sublinear terms. Nonlinear Anal. 2015, 121, 3–18. [Google Scholar] [CrossRef] [Green Version]
  20. Caraballo, T.; Herrera-Cobos, M.; Marín-Rubio, P. Time-dependent attractors for non-autonomous non-local reaction-diffusion equations. Proc. Roy. Soc. Edinburgh Sect. A 2018, 148A, 957–981. [Google Scholar] [CrossRef]
  21. Caraballo, T.; Herrera-Cobos, M.; Marín-Rubio, P. Robustness of time-dependent attractors in H1-norm for nonlocal problems. Discret. Contin. Dyn. Syst. Ser. B 2018, 23, 1011–1036. [Google Scholar]
  22. Caraballo, T.; Herrera-Cobos, M.; Marín-Rubio, P. Global attractor for a nonlocal p-laplacian equation without uniqueness of solution. Discret. Contin. Dyn. Syst. Ser. B 2017, 17, 1801–1816. [Google Scholar]
  23. Caraballo, T.; Herrera-Cobos, M.; Marín-Rubio, P. Asymptotic behaviour of nonlocal p-Laplacian reaction-diffusion problems. J. Math. Anal. Appl. 2018, 459, 997–1015. [Google Scholar] [CrossRef]
  24. Henry, D. Some infinite-dimensional Morse-Smale systems defined by parabolic partial differential equations. J. Differ. Equations 1985, 59, 165–205. [Google Scholar] [CrossRef] [Green Version]
  25. Carvalho, A.N.; Li, Y.; Luna, T.L.M.; Moreira, E. A non-autonomous bifurcation problem for a non-local scalar one-dimensional parabolic equation. Commun. Pure Appl. Anal. 2020, 19, 5181–5196. [Google Scholar]
  26. Carvalho, A.N.; Moreira, E. Stability and hyperbolicity of equilibria for a scalar nonlocal one-dimensional quasilinear parabolic problem. arXiv 2020, arXiv:2005.07857. [Google Scholar]
  27. Sell, G.R.; You, Y. Dynamics of Evolutionary Equations; Springer: New York, NY, USA, 2002. [Google Scholar]
  28. Robinson, J.C. Infinite-Dimensional Dynamical Systems: An Introduction to Dissipative Parabolic PDEs and the Theory of Global Attractors; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  29. Lions, J.L. Quelques Méthodes de Résolution des Problèmes aux Limites non Linéaires; Gauthier-Villar: Paris, France, 1969. [Google Scholar]
  30. Caraballo, T.; Marín-Rubio, P.; Robinson, J. A comparison between two theories for multi-valued semiflows and their asymptotic behaviour. Set-Valued Anal. 2003, 11, 297–322. [Google Scholar] [CrossRef]
  31. Kapustyan, O.V.; Pankov, V.; Valero, J. On global attractors of multivalued semiflows generated by the 3D Bénard system. Set-Valued Var. Anal. 2012, 20, 445–465. [Google Scholar] [CrossRef]
  32. Melnik, V.S.; Valero, J. On attractors of multi-valued semi-flows and differential inclusions. Set-Valued Anal. 1998, 6, 83–111. [Google Scholar] [CrossRef]
  33. Kapustyan, O.V.; Kasyanov, P.O.; Valero, J. Structure and regularity of the global attractor of a reacction-diffusion equation with non-smooth nonlinear term. Discret. Contin. Dyn. Syst. 2014, 32, 4155–4182. [Google Scholar] [CrossRef]
  34. Valero, J. On Lr-regularity of global attractors generated by strong solutions of reaction-diffusion equations. Appl. Math. Nonlinear Sci. 2016, 1, 375–390. [Google Scholar] [CrossRef] [Green Version]
  35. Caballero, R.; Carvalho, A.N.; Marín-Rubio, P.; Valero, J. Robustness of dynamically gradient multivalued dynamical systems. Discret. Contin. Dyn. Syst. Ser. B 2019, 24, 1049–1077. [Google Scholar] [CrossRef] [Green Version]
  36. Ball, J.M. Continuity properties and global attractors of generalized semiflows and the Navier-Stokes equations. J. Nonlinear Sci. 1997, 7, 475–502. [Google Scholar] [CrossRef] [Green Version]
  37. Angenent, S. The zero set of a solution of a parabolic equation. J. Reine Angew. Math. 1988, 390, 79–96. [Google Scholar]
  38. Da Costa, H.B.; Valero, J. Morse decompositions and Lyapunov functions for dynamically gradient multivalued semiflows. Nonlinear Dyn. 2016, 84, 19–34. [Google Scholar] [CrossRef]
  39. De Souza, P.N.; Nuno, J. Berkeley Problems in Mathematics; Springer: New York, NY, USA, 2002. [Google Scholar]
  40. Kadlec, J. Strong maximum principle for weak solutions of nonlinear parabolic differential inequalities. Časopis Pěst. Mat. 1967, 92, 373–391. [Google Scholar] [CrossRef]
Figure 1. a ( d ) non-decreasing.
Figure 1. a ( d ) non-decreasing.
Mathematics 09 00353 g001
Figure 2. a ( d ) whatever.
Figure 2. a ( d ) whatever.
Mathematics 09 00353 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Caballero, R.; Carvalho, A.N.; Marín-Rubio, P.; Valero, J. About the Structure of Attractors for a Nonlocal Chafee-Infante Problem. Mathematics 2021, 9, 353. https://doi.org/10.3390/math9040353

AMA Style

Caballero R, Carvalho AN, Marín-Rubio P, Valero J. About the Structure of Attractors for a Nonlocal Chafee-Infante Problem. Mathematics. 2021; 9(4):353. https://doi.org/10.3390/math9040353

Chicago/Turabian Style

Caballero, Rubén, Alexandre N. Carvalho, Pedro Marín-Rubio, and José Valero. 2021. "About the Structure of Attractors for a Nonlocal Chafee-Infante Problem" Mathematics 9, no. 4: 353. https://doi.org/10.3390/math9040353

APA Style

Caballero, R., Carvalho, A. N., Marín-Rubio, P., & Valero, J. (2021). About the Structure of Attractors for a Nonlocal Chafee-Infante Problem. Mathematics, 9(4), 353. https://doi.org/10.3390/math9040353

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