Stationary solutions and asymptotic behaviour for a chemotaxis hyperbolic model on a network

1.   Introduction

In this paper we consider a semilinear hyperbolic-parabolic system evolving on a finite planar network composed from nodes connected by $m$ oriented arcs $I_i$,

the system is complemented by initial, boundary and transmission conditions at the nodes (see Section 2).

We are interested in the study of stationary solutions and asymptotic behaviour of global solutions of the problem.

The above system has been proposed as a model for chemosensitive movements of bacteria or cells on an artificial scaffold [12]. The unknown $u$ stands for the cells concentration, $\lambda v$ is the average flux and $\phi$ is the chemo-attractant concentration. In particular, the model turns out to be useful to describe the process of dermal wound healing, when the stem cells in charge of the reparation of dermal tissue (fibroblasts) create an extracellular matrix and move along it to fill the wound, driven by chemotaxis; tissue engineers use artificial scaffolds, constituted by a network of crossed polymeric threads, inserting them within the wound to accelerate the process (see [13,20,25]). In the above mathematical model, the arcs of the graph mimic the fibers of the scaffold; each of them is characterized by a tipical velocity ${\lambda _i}$, a friction coefficient $\beta_i$, a diffusion coefficient $D_i$, and a production rate $a_i$ and a degradation one $b_i$; the functions $u_i, {\phi _i}$ are the densities of fibroblasts and chemoattractant on each arc.

Starting from the Keller-Segel paper [18] in 1970 until now, a lot of articles have been devoted to PDE models in domains of $\mathbb{R}^n$ for chemotaxis phenomena. The parabolic (or parabolic-elliptic) Patlak-Keller-Segel system is the most studied model [17,23,22]; in recent years, hyperbolic models have been introduced too, in order to avoid the unrealistic infinite speed of propagation of cells, occurring in parabolic models [8,9,14,23,24,1,15,16,11].

In [11] the Cauchy and the Neumann problems for the system in (1.1), respectively in $\mathbb{R}$ and in bounded intervals of $\mathbb{R}$, are studied, providing existence of global solutions and stability of constant states results.

Recently an interest in these mathematical models evolving on networks is arising, due to their applications in the study of biological phenomena and traffic flows, both in parabolic cases [2,6,21] and in hyperbolic ones [10,7,26,12,3].

We notice that the transmission conditions for the unknowns, at each inner node, which complement the equations on networks, are crucial characteristics of the model, since they are the coupling among the solution's components on each arc.

Most of the studies carried out until now, consider continuity conditions at each inner node for the density functions [7,6,21]; nevertheless, the eventuality of discontinuities at the nodes seems a more appropriate framework to decribe movements of individuals or traffic flows phenomena [5].

For these reason in [12], transmission conditions which link the values of the density functions at the nodes with the fluxes, without imposing any continuity, are introduced; these conditions guarantee the fluxes conservation at each inner node, and, at the same time, the m-dissipativity of the linear spatial differential operators, a crucial property in the proofs of existence of local and global solutions contained in that paper.

In this paper we focus our attention on stationary solutions to problem (1.1) complemented by null fluxes boundary conditions and by the same transmission conditions of [12] (see next section and Section 3 in [12] for details). We consider acyclic networks and we prove the existence and uniqueness of the stationary solution with fixed mass of cells $\sum\limits_{i = 1}^m {{\rm{ }}\int_{{I_i}} {{u_i}} } (x)dx$, under the assumption that the mass is suitably small. If the quantity $\frac{a_i}{b_i}$ does not vary with the index $i$, we easily show that such solution is a constant state on the whole network. We notice that, in the case of acyclic networks, although the transmission conditions do not set the continuity of the density $u$ at the inner nodes, the fluxes conservation at those nodes and the boundary null fluxes conditions imply the absence of jumps discontinuities at the inner vertexes, for the component $u$ of a stationary solution.

For general networks and the parameters $a_i$ and $b_i$ in the same range as above, it is easy to show that, for any fixed mass, a stationary solution constant on the whole network exists and the constant values of the densities are determined by the mass. In this case we also obtain a uniqueness result: in the set of stationary solutions with small density $u$ in $H^1$-norm and fixed mass of cells, the constant state on the network is the unique element.

Finally we study the large time behaviour of global solutions on general networks, when the ratio between $a_i$ and $b_i$ is constant. We consider initial data with fixed small mass, which are small perturbations of the constant state on the network with the same mass, then we prove that such state is the asymptotic profile of the solutions corresponding to the data. So, we point out that, for small global solutions to our problem, the discontinuities at the inner nodes vanish when $t$ goes to infinity, since their asymptotic profiles are continuous functions on the whole network.

The study of the asymptotic behaviour provide informations about the evolution of a small mass of individuals moving on a network driven by chemotaxis: suitable initial distributions of individuals and chemoattractant, for large time evolve towards constant distributions on the network, preserving the mass of individuals.

We recall that the stability of the constant solutions to this system, considered on bounded interval in $\mathbb{R}$, is studied in [11] and stationary solutions and asymptotic behaviour for a linear system of uncoupled conservation laws on network are studied in [19].

Finally, in [3] the authors introduce a numerical scheme to approximate the solutions to the problem (2.1); in that paper transmission conditions are set for the Riemann invariants of the hyperbolic part of the system, $w_i^{\pm} = \frac 1 2 (u_i\pm v_i)$, and are equivalent to our ones for some choices of the transmission coefficients. The tests presented there, in the case of acyclic graph and dissipative transmission coefficients, show an asymptotic behaviour of the solutions which agrees with our theoretical results.

The paper is organized as follows. In Section 2 we give the statement of the problem and, in particular, we introduce the transmission conditions, while in Section 3 we prove the results about existence and uniqueness of stationary solutions. The last section is devoted to study the asymptotic behaviour of solutions; the results obtained in this section constitute the sequel and the development of the result of existence of global solutions in [12] and the proofs are based on the same techniques and use simple modifications of the a priori estimates obtained in [12].

2.   Statement of the evolution problem

We consider a finite connected graph $\mathcal G = (\mathcal V, \mathcal A)$ composed by a set $\mathcal V$ of $n$ nodes (or vertexes) and a set $\mathcal A$ of $m$ oriented arcs, $\mathcal A = \{I_i:i\in \mathcal M = \{1, 2, ..., m\}\}$.

Each node is a point of the plane and each oriented arc $I_i$ is an oriented segment joining two nodes.

We use $e_j$, $j\in \mathcal J$, to indicate the external vertexes (or boundary vertexes) of the graph, i.e. the vertexes belonging to only one arc, and by $I_{i(j)}$ the external arc outgoing or incoming in the external vertex $e_j$.

Moreover, we use $N_\nu$, $\nu\in\mathcal P$, to denote the inner nodes; for each of them we consider the set of incoming arcs $\mathcal A_{in}^\nu = \{I_i:i\in \mathcal I^\nu\}$ and the set of the outgoing ones $\mathcal A_{out}^\nu = \{I_i:i\in \mathcal O^\nu\}$; let $\mathcal M^\nu = \mathcal I^\nu\cup\mathcal O^\nu$.

In this paper, a path in the graph is a sequence of arcs, two by two adjacent, without taking into account orientations. Moreover, we call acyclic a graph which does not contain cycles: for each couple of nodes, there exists a unique path with no repeated arcs connecting them (an example of acyclic graph is in Fig. 1).

Each arc $I_i$ is considered as a one dimensional interval $(0, L_i)$. A function $f$ defined on $\mathcal A$ is a m-tuple of functions $f_i$, $i\in \mathcal M$, each one defined on $I_i$; $ f_i(N_\nu)$ denotes $f_i(0) $ if $N_\nu$ is the initial point of the arc $I_i$ and $f_i(L_i)$ if $N_\nu$ is the end point, and similarly for $f(e_j)$.

We set $ L^p(\mathcal A): = \{f: f_i\in L^p(I_i), i\in \mathcal{M}\}$, $ H^s(\mathcal A): = \{f :f_i\in H^s(I_i), i\in\mathcal{M}\}$ and

We consider the evolution of the following one-dimensional problem on the graph $\mathcal G$

where $ a_i\geq 0, {\lambda _i}\ b_i, D_i, \beta_i >0 $.

We complement the system with the initial conditions

the boundary conditions at each outer point $e_j$ are the null flux conditions

In addition, at each inner node $N_\nu$ we impose the following transmission conditions for the unknown $\phi$

which imply the continuity of the flux at each node, for all $t>0$,

For the unknonws $v$ and $u$ we impose the transmission conditions

These conditions ensure the conservation of the flux of the density of cells at each node $N_\nu$, for $t>0$,

which corresponds to the conservation of the total mass

i.e. no death nor birth of individuals occours during the observation.

Motivations for the constraints on the coefficients in the transmission conditions can be found in [12].

Finally, we impose the following compatibility conditions

Existence and uniqueness of local solutions to problem (2.1)-(2.7),

are achieved in [12] by means of the linear contraction semigroup theory coupled with the abstract theory of nonhomogeneous and semilinear evolution problems; in fact, the transmission conditions (2.5) and (2.6) allows us to prove that the linear differential operators in (2.1) are m-dissipative and then, to apply the Hille-Yosida-Phillips Theorem (see [4]). The existence of global solutions when the initial data are small in $(H^1(\mathcal A))^2 \times H^2(\mathcal A)$ norm is proved [12] too; this result holds under the further assumption

3.   Non negative small stationary solutions on acyclic networks

In this section we approach the question of existence and uniqueness of stationary solutions of problem (2.1)-(2.8), with fixed mass

in the case of an acyclic network (see Section 2). We look for stationary solutions $(u, v, \phi)\in (H^1(\mathcal A))^2\times H^2(\mathcal A)$.

Obviously, the flux $v$ of a stationary solution has to be constant on each arc and has to be null on the external arcs; in the case of acyclic networks, the boundary and transmission conditions (2.3), (2.6) force it to be null on each arc. In order to prove this fact we consider an internal arc $I_j$ and its initial node $N_\mu$; we consider the set

(see Fig. 1: if, for example, $j = 9$ then $\mu = 4$, $\mathcal Q = \{1, 2, 3, 5\}$ and the arcs in bold type form the path which links the nodes $N_5$ and $N_4$).

At each node the conservation of the flux of the density of cells, stated in Section 2, holds; then

Since, for all $i\in\mathcal{M}$, $v_i(x)$ is constant on $I_i$ and $v_i(x) = 0$ if $I_i$ is an external arc, the above equality reduces to

then $v_j(x) = 0$ for all $x\in I_j$.

The previous result implies that stationary solutions must have the form $(u, 0, \phi)$, where $u$ and $\phi$ have to verify the system

with the boundary condition at each outer point $e_j$, $j\in\mathcal J$,

and the transmission conditions, at each inner node $N_\nu$,

For each fixed inner node $N_\nu$, let $k\in\mathcal M^\nu$ be the index in condition (2.8) and let consider the transmission relations, for $i\in\mathcal M^\nu$, $i\neq k$,

the assumptions on $K^\nu_{kj}$ in (2.8) ensure that the matrix of the coefficients of this linear system in the unknowns $(u_j(N_\nu)-u_k(N_\nu))$, $j\neq k$, is non singular (if $k = 1$ it is immediate to check that it has strictly dominant diagonal). Then the condition (3.3) can be rewritten as

Now we fix $\mu_0\geq 0$ and we look for stationary solutions such that

notice that for the evolution problem, the quantity $ \sum\limits_{i \in \mathcal{M}} {\int_{{I_i}}{u_i}(x, t)} dx$ is preserved for all $t\geq 0$, thanks to the transmission conditions (2.6).

Integrating the first equation in (3.1) we can rewrite problem (3.1)-(3.6) as the following elliptic problem on network:

Find $C = (C_1, C_2, ..., C_m)$ and $\phi\in H^2(\mathcal A)$ such that

We consider the linear operator $A:D(A)\to L^2(\mathcal A)$,

then the equation in (3.7) and the boundary and transmission conditions for $\phi$ can be written as

where, for $i\in\mathcal{M}$, $F_i(\phi(x), C) = a_iC_i \exp\left( \frac {\phi_i(x)} {\lambda_i}\right)$.

We are going to prove the existence and uniqueness of solutions to the problem (3.7) by using the Banach Fixed Point Theorem; in order to do this we need some preliminary results about the linear equation

where $f\in H^2(\mathcal{A})$ is a given function, $C^f = (C^f_1, ..., C^f_m)$ and $C^f_i$ are non-negative given real constants.

The existence and uniqueness of the solution $\phi\in H^2(\mathcal{A})$ to the above problem (for a general $F\in L^2(\mathcal{A})$ and a general network) is showed in [12], by Lax-Milgram theorem, in the proof of Proposition 4.1; here, we need to prove some properties holding for the solution in the case of acyclic graphs, under suitable assumptions on $f$ and $C^f_i$.

The transmission conditions (2.5) imply the following equality which will be useful in the next proofs:

Let $ |\mathcal A|: = \sum\limits_{i\in\mathcal M} | I_i|$ and $\Vert g\Vert_{\infty}: = \max \{ \Vert g_i\Vert_{\infty}, i\in\mathcal{M}\}$.

Lemma 3.1. Let $\mathcal G = (\mathcal V, \mathcal{A})$ be an acyclic network, let $f\in H^2(\mathcal{A})$ and let $C_i^f$ be non-negative real numbers, for $i\in \mathcal{M}$. Then the solution $\phi$ to problem (3.8), (3.10) is non-negative. Moreover, if

then

if (3.12) holds and

then there exists a quantity $K_{\mu_0} = K_{\mu_0}(a_i, b_i, D_i, {\lambda _i}, |\mathcal A|, \mu_0)$, depending only on the parameters appearing in brackets, infinitesimal when $\mu_0$ goes to zero, such that

Proof. Let consider a function $\Gamma \in C^1(\mathbb{R})$, strictly increasing in $(0, +\infty)$, and let $\Gamma (y) = 0$ for $y\leq 0$; following standard methods for the proofs of the maximum principle for elliptic equations and setting $F_i(x) = F_i(f(x), C^f)$, we obtain

As regard to the first term, we can argue as in (3.11), taking into account the properties of $\Gamma$,

the above inequality and the non-negativity of $F_i$ imply that

so that, thanks to the properties of $\Gamma$, we can conclude that ${\phi _i}(x)\geq 0$ for all $i\in\mathcal{M}$.

By integration of the equation (3.10), taking into account (3.4) and (3.2), we obtain

which implies

In order to obtain (3.13), first we notice that, if $I_j$ is an external arc, then the following inequality holds

Then we consider an internal arc $I_j$ and its initial node $N_\mu$ and the sets

(see Fig. 1: if, for example, $j = 9$, then $\mu = 4$, $\mathcal Q = \{1, 2, 3, 5\}$, $\mathcal S = \{1, 2, 3, 4, 5, 6, 8, $ $10, 11, 12\}$); at each node the conservation of the flux, stated in Section 2 as a consequence of the transmission conditions, holds; then

Let $x$ be a point on the arc $I_j$ (see Fig. 1 for $j = 9$, $\mu = 4$) and $I_j^x$ be the part of $I_j$ which connects $N_\mu$ and $x$; then, using the above equality and the boundary conditions (2.4), we have

Then $ D_j \Vert \phi_{j_x}\Vert_{\infty}\leq 2 \mu_0 \max \{a_i\}_{i\in\mathcal{M}}$ for all $j\in\mathcal{M}$ and we obtain (3.13) which implies

and

moreover, by Sobolev embedding theorem, we obtain

where $K_1 = K_1(a_i, b_i, D_i, |\mathcal A|) $ is a suitable constant.

The estimates for the function $\phi_{xx}$ follow by using the equation (3.10); first, using (3.12) and (3.17), we obtain

then, using (3.11), we have

and the embedding of $W^{1, 1}(I_i) $ in $L^{\infty}(I_i)$ gives

where $K_2 = K_2(a_i, b_i, D_i, {{\lambda _i}})$ is a suitable constant.

Finally, the inequalities (3.18), (3.20)-(3.24) imply the inequalities (3.15) in the claim.

Now we can prove the following theorem.

Theorem 3.1. Let $\mathcal G = (\mathcal V, \mathcal A)$ be an acyclic network. There exists $\epsilon>0$ such that, if $ 0\leq\mu_0\leq \epsilon$, then problem (2.1)-(2.8) has a unique stationary solution satisfying (3.6); the solution has the form

where ${\phi _i}(x)\geq 0$ and $C_i$ are nonnegative constants such that $u_j(N_\nu) = u_i(N_\nu)$ for all $\nu\in\mathcal P$, $ i, j\in \mathcal M^\nu$.

Proof. First we notice that, if a stationary solution $(u, v, \phi)$ satisfying (3.6) exists, then $u$ is non-negative, since the costants $C_i$ in (3.7) must have the same sign, so that they have to be non-negative to satisfy the condition $\mu_0 \geq 0$; arguing as in the proof of Lemma 3.1 we prove that $\phi$ is non-negative too. If $\mu_0>0$ then $u$ and $\phi$ are positive functions.

We are going to use a fixed point technique. Given $\phi^0 \in D(A)$, we want to define a function $u^0(x)$ on the network, such that, for $i\in\mathcal{M}$,

where the constants $C_i^{\phi^0}$ satisfy the following linear system composed by the last conditions in (3.7)

The system (3.25), (3.26) has a unique solution; actually, since the network has no cycles, the system (3.25) has $\infty^1$ solutions $ C^{\phi^0, \alpha} = (\alpha, \alpha \delta_2, \alpha \delta_3, ..., \alpha \delta_m)$, $\alpha\in\mathbb{R}$, where $\delta_i$ are suitable coefficients, and the condition (3.26) determines the value of $\alpha$.

In order to give an explicit expression for the coefficients $\delta_i$ we consider an arc, $I_1$, and we define

Let $N_\mu$ one of the extreme points of $I_1$, then we define the function $u^0$ on the other arcs which are incident with $N_\mu$ in such a way to verify the equalities in (3.25) for the node $N_\mu$,

i.e. we set $ C^{\phi^0}_j = \alpha \exp\left(\frac {\phi_1^0(N_\mu)}{ \lambda_1}\right)\exp\left( -\frac {\phi_j^0(N_\mu)}{\lambda_j}\right) $, $j\in\mathcal M^\mu, j\neq 1$.

This procedure can be iterated at each node reached by one of the arcs $I_j$, $j\in \mathcal{M}^\mu$, and at the other extreme point of $I_1$, if it is an internal arc, and so on, defining in this way the function $u^0$ on each arc of the network. Notice that this construction is possible since there are no cycles in the graph. The function $u$ can be expressed, on each arc of the network, as follows (if it is the case, renumbering in suitable way the arcs and the nodes): let consider the path from the outer point $e_1$ to an inner node $N_{h-1}$, composed from the arcs $I_i$, $i = 1, ..., h-1$, (passing through the vertexes $N_i$, $i = 1, ..., h-1$), and let $I_h$ be an arc incident with the node $N_{h-1}$, not belonging to the path (see Fig. 2 where h = 5 and the highlighted arcs form the path); following the procedure described before, after setting

we define

The quantity $\alpha$ is fixed in such a way to verify the last condition in (3.7),

so that, for all $i\in \mathcal{M}$,

Let $G$ be the operator defined in $D(A)$ such that $\phi^1: = G(\phi^0)$ is the solution of problem (3.10) where $f = \phi_0$ and $ C_i^f = C_i^{\phi^0}$ for $i\in\mathcal{M}$,

let $K_{\mu_0} = K _{\mu_0}(a_i, b_i, D_i, {\lambda _i}, |\mathcal A|, {\mu_0})$ be the quantity in Lemma 3.1, and let

equipped with the distance $d$ generated by norm of $H^2(\mathcal{A})$; $(B_{\mu_0}, d)$ is a complete metric space. From the lemma we know that solutions to problem (3.8)-(3.12) have to belong to $ B_{\mu_0}$, then $G(B_{\mu_0})\subseteq B_{\mu_0}$; next we are proving that, if $\mu_0$ is small enough, then $G$ is a contraction in $B_{\mu_0}$.

We consider $\phi^0, \overline {\phi}^0 \in B_{\mu_0}$ and the corresponding $u^0, \overline {u}^0$ and $\phi^1, \overline {\phi}^1$; using the equation satisfied by $\phi^0$ and $\overline {\phi}^0$, for all $i\in\mathcal{M}$ we can write

using (3.11), from (3.28) we infer that

We set

we have

In order to treat the above quantity we have to consider that, for all $g\in B_{\mu_0}$, $E^g_i(x)\geq 1$, $ {J^g_i}\geq |I_i|$ and there exists a constant $K_6 = K_6(K_{\mu_0}, {\lambda _i})$, increasing with $\mu_0$, such that, for all $i\in\mathcal{M}$

The above inequalities can be used in (3.30) so that (3.29) implies

where $K_7$ increases with $\mu_0$; hence, for $\mu_0$ small enough, $G$ is a contraction on $B_{\mu_0}$ and we can use the Banach Fixed Point Theorem.

Let $\phi$ be the unique fixed point of $G$ in $B_{\mu_0}$ and let $C^\phi = (C_1^\phi, C_2^\phi, ..., C_m^\phi)$ where $C_i^{\phi}$, for $i\in\mathcal{M}$, are computed as in (3.27); then $(\phi, C^{\phi})$ is the unique solution to Problem $(3.7)$ and the claim is proved.

For any constant $U\geq 0$, the triple $(U, 0, \frac {a_i}{b_i} U)$ satisfies the equations in (2.1) on the arc $I_i$. Let $Q$ be a real non-negative number; if $\frac {a_i}{b_i} = Q$ for all $i\in\mathcal{M}$, then the same triple $(U, 0, Q U)$ satisfies the equations on each arc $I_i$ and it is a stationary solution to the problem (2.1)-(2.8). Then, as a consequence of the previous theorem, we have the following proposition, with $ \epsilon $ as in the theorem.

Proposition 3.1. Let $\mathcal G = (\mathcal V, \mathcal A)$ be an acyclic network. If $\frac {a_i} {b_i} = Q$ for all $i\in \mathcal{M}$ and $0\leq\mu_0\leq \epsilon$, then the unique stationary solution to problem (2.1)-(2.8), (3.6) is the constant solution $\left(\frac{\mu_0}{|\mathcal{A}|}, 0, Q \frac{\mu_0}{|\mathcal{A}|}\right)$.

Remark 3.1. For general networks, when the value of $\frac {a_i}{b_i} = Q$ on each arc, the stationary solution of Proposition 3.1 always exists. More precisely, if $\frac {a_i} {b_i} = Q$, in the class of the functions $(u, v, \phi)$ which are constant on each arc, the stationary solution $\left(\frac{\mu_0}{|\mathcal{A}|}, 0, Q \frac{\mu_0}{|\mathcal{A}|}\right)$ is the unique stationary solution with mass $\mu_0$; this fact is true without any restrictions on the value of $\mu_0$ and on the structure of the network. Actually, if we assume that $u$ is constant on each arc, then, using the equations, we infer that, on each arc, $\phi_x(x)$ is constant too, hence $\phi_{xx} = 0$ and $\phi(x)$ is constant. Then $v(x) = 0$ on each arc; hence, arguing as at the beginning of this section, we obtain that $u$ is continuous on the network.

In the next proposition we are going to prove that, in a set of small solutions, such stationary solution is the unique one with fixed mass $\mu_0$.

Proposition 3.2. Let $\frac {a_i} {b_i} = Q$ for all $i\in\mathcal{M}$ and let $(u, v, \phi)\in (H^1(\mathcal{A}))^2\times H^2(\mathcal{A}) $ be a stationary solution of problem (2.1)-(2.8), (3.6). There exists $\epsilon_0>0, $ depending on ${\lambda _i}, a_i, b_i, D_i, \beta_i, |\mathcal{A}|$, such that, if $\Vert u\Vert_{H^1}\leq \epsilon_0$, then $(u, v, \phi) = \left(\frac{\mu_0}{|\mathcal{A}|}, 0, Q \frac{\mu_0}{|\mathcal{A}|}\right)$.

Proof. We set $H: = \Vert u\Vert_{H^1}$. The transmission conditions (2.6) imply that

so, by using the first two equations in (2.1), we obtain

and

the above inequalities implies the following one

where $K_0$ is a positive constant depending on the parameters ${\lambda _i}, \beta_i$ and the Sobolev embedding costant.

The transmission conditions (2.5) imply that

moreover, the assumption (2.8) imply that, for each $\nu\in\mathcal P$, for suitable coefficients $\theta^\nu_{ij}$ and suitable $k\in \mathcal{M}^\nu$,

(see Lemma 5.9 in [12]); then, by the last equation in (2.1), arguing as in the proof of Proposition 5.8 in [12], we obtain

where $K_1$ is a positive constant depending on the parameters $D_i, a_i, b_i, \theta^\nu_{ij}$.

By inequalities (3.32) and (3.33) we deduce the following one

which, for small $H$, implies $ \Vert v\Vert_2, \Vert u_x\Vert_2 = 0$.

In the cases when $\frac {a_i} {b_i}$ depends on the arc in consideration, stationary solutions with the component $u$ constant on each arc, can be inadmissible. As we showed before, $v$ should be zero, $u$ should be constant on the whole network and $\phi$ should be constant on each arc,

Therefore the transmission conditions, for each $\nu\in\mathcal P$,

are constraints on the relations between the parameters of the problem which have to hold if the constant stationary solution exists.

For example, in the case of two arcs, if $\frac {b_2} {a_2}\neq\frac {b_1} {a_1}$ (and $0<\mu_0\leq \epsilon$), the stationary solution can not be constant on the arcs, since the trasmission condition at the node,

cannot be satisfied.

Hence, in the cases when $\frac {a_i} {b_i}$ depends on the arc in consideration, if $(u, v, \phi)$ is the stationary solution in Theorem 3.1, then $u$ is a continuous function on all the network but it is not constant on each arc.

4.   Asymptotic behaviour

In this section we are going to show that the constant stationary solutions previously introduced, provide the asymptotic profiles for a class of solutions to problem (2.1)-(2.8). We recall that existence and uniqueness of global solutions

to such problem is proved in [12], when the initial data are sufficiently small in $(H^1(\mathcal{A}))^2\times H^2(\mathcal{A})$ norm and the following condition holds

in particular it is proved that the functional $F$ defined by

is uniformly bounded for $T>0$.

Here and below we use the notations

Now we assume (4.2), we fix $\overline \mu\geq 0$ and we consider the constant stationary solution, $(\overline {u}, 0, \overline {\phi})$, to problem (2.1)-(2.8), such that $ \overline {u} |\mathcal{A}| = \overline \mu$; moreover let $(\tilde u_0, \tilde v_0, \tilde \phi_0)\in (H^1(\mathcal{A}))^2\times H^2(\mathcal{A})$ be a small perturbation of $(\overline {u}, 0, \overline {\phi})$, i.e., setting $u_0: = \tilde u_0-\overline {u}, v_0: = \tilde v_0, \phi_0: = \tilde\phi_0-\overline {\phi}\ $, the $(H^1(\mathcal{A}))^2\times H^2(\mathcal{A})$ norm of $(u_0, v_0, \phi_0)$ is bounded by a suitable small $\epsilon_0>0$.

If $(\tilde u, \tilde v, \tilde \phi)$ is the solution to problem (2.1)-(2.8) with initial data $(\tilde u_0, \tilde v_0, \tilde \phi_0)$ and $u: = \tilde u-\overline {u}, v: = \tilde v, \phi: = \tilde\phi-\overline {\phi}\ $, then $(u, v, \phi)$ is solution to the system

complemented with the conditions (2.2)-(2.8) and initial data $(u_0, v_0, \phi_0) $ defined above.

The existence and uniqueness of local solutions to this problem can be achieved by means of semigroup theory, following the method used in[12], with little modifications.

On the other hand, if we assume that $\overline {u} $ is suitably small, the method used in that paper to obtain the global existence result in the case of small initial data can be used here too, modifying the estimates in order to treat the further term in the second equation and then using the smallness of $\overline {u}$.

Below we list a priori estimates holding for the solutions to problem (4.4), (2.2)-(2.7); we don't give the proofs since they are equal to those in [12], in Section 5, except for easy added calculations to treat the term $\overline {u} {\phi_{i}}_x$.

Proposition 4.1. Let $(u, v, \phi)$ be a local solution to problem (4.4), (2.2)-(2.7),

then

a)

b)

c)

d)

e)

f)

g)

h) if (2.8) and (4.2) hold, then

for suitable costants $C$.

The estimates in the previous proposition allow to prove the following theorem about the existence of global solutions to problem (4.4), (2.2)-(2.8).

Let $F_T(u, v, \phi)$ be the functional defined in (4.3).

Theorem 4.1. Let (4.2) hold. There exists $\epsilon_0, \epsilon_1>0$ such that, if

then there exists a unique global solution $(u, v, \phi)$ to problem (4.4), (2.2)-(2.8),

Moreover, $F_T(u, v, \phi)$ is bounded, uniformly in $T$.

Proof. It is sufficient to prove that the functional $F_T(u, v, \phi)$ is bounded, uniformly in $T$.

We notice that each term in $F_T^2(u, v, \phi)$ is in the left hand side of one of the estimates in Proposition 4.1, therefore, arranging all the estimates, we can prove the following inequality

taking into account also that, on the right hand side of the estimates, the quadratic terms (not involving initial data) which have not the coefficient $\overline {u}$, can be bounded by means of cubic ones.

If $\overline {u}$ is sufficiently small, the previous inequality implies

for suitable positive constants $c_i$.

It is easy to verify that, if $y_0$ is a sufficiently small positive real number and $ h(y) = c_5 y^3-y^2 +c_4 y_0$ then there exists $0<\overline y<\frac 2 {3 c_5}$ such that $h(y)>0$ in $[0, \overline y)$ and $h(y)<0$ in $(\overline y, \frac 2 {3 c_5}]$.

Then we can conclude that, if $F_0$ is suitably small, then $F_T$ remains bounded for all $T>0$.

The above result, in particular the uniform, in time, boundedness of the functional $F_T$, allow us to prove the theorem below.

Let (4.2) hold and let $(\overline {u}, 0, \overline {\phi})$ be the constant stationary solution to problem (2.1)-(2.8) such that $ \overline {u} |\mathcal{A}| = \overline \mu$; moreover, let $\mathcal C(\mathcal{A})$ be the set of the funcions $f$ such that $f_i\in\mathcal C(\overline I_i)$ for $i\in \mathcal{M}$.

Theorem 4.2. Let (4.2) hold. There exist $\epsilon_0, \epsilon_2>0$ such that, if

then problem (2.1)-(2.8) has a unique global solution $(u, v, \phi)$,

and, for all $i\in\mathcal{M}$,

Proof. Let $(u, v, \phi)$ be the local solution to problem (2.1)-(2.8) and let

we already noticed that $(\hat u, \hat v, \hat \phi)$ is the local solution to system (4.4) complemented by the initial condition $(u_0-\overline {u}, v_0, \phi_0-\overline {\phi})$ and the same boundary and transmission condition given by (2.3)-(2.8) for system (2.1).

For suitable $\epsilon_0, \epsilon_2$ the assumptions of Theorem 4.1 are satisfied, then we obtain the uniform boundedness of the functional $F_T(\hat u, \hat v, \hat \phi)$, for $T\in [0, +\infty)$. Hence the set $\{\hat u(t), \hat v(t), \hat \phi(t)\}_{t\in[0, +\infty)}$ is uniformly bounded in $(H^1(\mathcal{A}))^2\times H^2(\mathcal{A})$; thus, if we call $E_s$ the set of accumulation points of $\{\hat u(t), \hat v(t), \hat \phi(t)\}_{t\geq s}$ in $(C(\mathcal{A}))^2\times C^1(\mathcal{A})$, then $E_s$ is not empty and $ E: = \cap_{s\geq 0} E_s \neq \emptyset$.

Let $(U(x), V(x), \Phi(x))$ be such that, for a sequence $t_n\to +\infty$,

In order to identify the limit functions we notice that $ \sum\limits_{i \in \mathcal{M}} \int_{{I_i}} U_i(x) dx = 0$, since $ \sum\limits_{i \in \mathcal{M}} \int_{{I_i}} \hat u(x, t_n), dx = 0$ for all $t_n$.

Moreover, since $ \hat v_i\in H^1(I_i\times (0, +\infty))$ for all $i\in\mathcal{M}$, if we set

then $\omega_i\in H^1((0, +\infty))$ and, as a consequence, $ \lim\limits_{t\to+\infty} \omega_i(t) = 0$.

As $ \lim\limits_{n\to +\infty}\Vert \hat v_i(\cdot, t_n)\Vert_2 = \Vert V _i(\cdot)\Vert_2$, we obtain $\Vert V\Vert_2 = 0$.

The same argument can be applied to the functions $\hat \phi_{i_x}$, for all $i\in\mathcal{M}$, since they belongs to $H^1(I_i\times (0, +\infty))$. Finally, it can be applied to the functions $a_i\hat u_i-b_i\hat \phi_i$ since $\hat \phi_{i_t}, \hat \phi_{i_{xx}}, \hat u_{i_x}, \hat \phi_{i_x}\in L^2(I_i\times (0, +\infty))$, thanks to the uniform boundedness of $F_T(\hat u, \hat v, \hat \phi) $ and to estimate f) in Proposition 4.1.

As a consequence we have that

where $\overline \Phi_i$ are real numbers, so that the limit function is given by $(\frac {b_i}{a_i} \Phi_i, 0, \Phi_i)$ in each interval $I_i$, for all $i\in\mathcal{M}$. It is easily seen that such function is a stationary solution to problem (2.1)-(2.8), which is constant in each arc $I_i$; in particular it verifies the transmission conditions since $(\hat u, \hat v, \hat \phi)$ verifies them and the convergence result (4.5) holds.

The condition $ \sum\limits_{i \in \mathcal{M}} \int_{{I_i}} U_i(x) dx = 0$ and Remark 3.1 imply that $\Phi_i = 0$ for all $i\in\mathcal{M}$, so that we can conclude that the unique possible limit function is $(U(x), V(x), \Phi(x)) = (0, 0, 0)$; this fact proves the claimed convergence results.

5.   Conclusions and perspectives for the future

The main features of the present work are: $a)$ the proof of the existence and uniqueness of stationary solutions with fixed small mass to problem (2.1)-(2.6) considered on acyclic networks; $b)$ the proof of the stability of particular stationary solutions, the constant states on the whole network, when the their masses are small and the quantity $\frac {a_i} {b_i} $ does not vary with the index $i$, for general networks. We can conclude that, in this range of parameters, although the transmission conditions do not impose continuity of the densities at the internal nodes, for suitable initial data the asymptotic profiles of the solutions are continuous functions, constant on the network.

These results are useful in describing the large time behaviour of small masses of individuals moving on networks driven by chemotaxis.

For the future, our aim is approaching the same questions when the system (2.1) is complemented by non-null fluxes conditions at the boundaries, which provide models for different situations at the outer nodes, in order to describe the features of the behaviour of cells moving along the arcs searching food. We notice that the condition $v_i(x) = 0$ for $x\in{I_i}$, for all $i\in \mathcal{M}$, which prevents the presence of jumps for the density $u$ at the inner nodes, in this case is no longer necessary for stationary solutions on acyclic networks.