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

  • Received: 01 November 2016 Revised: 01 October 2017
  • Primary: 35R02; Secondary: 35M33, 35L50, 35B40, 35Q92

  • This paper approaches the question of existence and uniqueness of stationary solutions to a semilinear hyperbolic-parabolic system and the study of the asymptotic behaviour of global solutions. The system is a model for some biological phenomena evolving on a network composed by a finite number of nodes and oriented arcs. The transmission conditions for the unknowns, set at each inner node, are crucial features of the model.

    Citation: Francesca R. Guarguaglini. 2018: Stationary solutions and asymptotic behaviour for a chemotaxis hyperbolic model on a network, Networks and Heterogeneous Media, 13(1): 47-67. doi: 10.3934/nhm.2018003

    Related Papers:

    [1] Francesca R. Guarguaglini . Stationary solutions and asymptotic behaviour for a chemotaxis hyperbolic model on a network. Networks and Heterogeneous Media, 2018, 13(1): 47-67. doi: 10.3934/nhm.2018003
    [2] Raul Borsche, Axel Klar, T. N. Ha Pham . Nonlinear flux-limited models for chemotaxis on networks. Networks and Heterogeneous Media, 2017, 12(3): 381-401. doi: 10.3934/nhm.2017017
    [3] Sara Monsurrò, Carmen Perugia . Homogenization and exact controllability for problems with imperfect interface. Networks and Heterogeneous Media, 2019, 14(2): 411-444. doi: 10.3934/nhm.2019017
    [4] Piotr Biler, Grzegorz Karch, Jacek Zienkiewicz . Morrey spaces norms and criteria for blowup in chemotaxis models. Networks and Heterogeneous Media, 2016, 11(2): 239-250. doi: 10.3934/nhm.2016.11.239
    [5] Avner Friedman . PDE problems arising in mathematical biology. Networks and Heterogeneous Media, 2012, 7(4): 691-703. doi: 10.3934/nhm.2012.7.691
    [6] James Nolen . A central limit theorem for pulled fronts in a random medium. Networks and Heterogeneous Media, 2011, 6(2): 167-194. doi: 10.3934/nhm.2011.6.167
    [7] Magali Tournus, Aurélie Edwards, Nicolas Seguin, Benoît Perthame . Analysis of a simplified model of the urine concentration mechanism. Networks and Heterogeneous Media, 2012, 7(4): 989-1018. doi: 10.3934/nhm.2012.7.989
    [8] Marek Fila, John R. King . Grow up and slow decay in the critical Sobolev case. Networks and Heterogeneous Media, 2012, 7(4): 661-671. doi: 10.3934/nhm.2012.7.661
    [9] Jacek Banasiak, Proscovia Namayanja . Asymptotic behaviour of flows on reducible networks. Networks and Heterogeneous Media, 2014, 9(2): 197-216. doi: 10.3934/nhm.2014.9.197
    [10] Jérôme Coville, Nicolas Dirr, Stephan Luckhaus . Non-existence of positive stationary solutions for a class of semi-linear PDEs with random coefficients. Networks and Heterogeneous Media, 2010, 5(4): 745-763. doi: 10.3934/nhm.2010.5.745
  • This paper approaches the question of existence and uniqueness of stationary solutions to a semilinear hyperbolic-parabolic system and the study of the asymptotic behaviour of global solutions. The system is a model for some biological phenomena evolving on a network composed by a finite number of nodes and oriented arcs. The transmission conditions for the unknowns, set at each inner node, are crucial features of the model.



    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$,

    $ \left\{ tui+λixvi=0,tvi+λixui=uixϕiβivi,t0,xIi,i=1,...,m;tϕi=Dixxϕi+aiuibiϕi,
    \right.$
    (1.1)

    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].

    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).

    Figure 1. 

    Example of acyclic network; the highlighted arcs form the path linking the nodes $N_4$ and $N_5$.

    .

    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

    $\Vert f\Vert_2: = \sum\limits_{i \in \mathcal M}\Vert f_i\Vert_2, \ \Vert f\Vert_{H^s}: = \sum\limits_{i \in \mathcal M} \Vert f_i\Vert_{H^s}.$

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

    $ \left\{ tui+λixvi=0,tvi+λixui=uixϕiβivi,t0, xIi, iM,tϕi=Dixxϕi+aiuibiϕi,
    \right.$
    (2.1)

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

    We complement the system with the initial conditions

    $ u_{i0}, v_{i0}\in H^1(I_i), \ \phi_{i0}\in H^2(I_i) \ \textrm{ for } i\in\mathcal{M}\ ;$ (2.2)

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

    $ v_{i(j)}(e_j, t) = 0, \;\;\;\; t > 0, \ j\in\mathcal J, $ (2.3)
    $ {\phi_{i(j)}}_x(e_j, t) = 0 \;\;\;\; t > 0 , \ j\in \mathcal J .$ (2.4)

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

    $\left\{Diϕix(Nν,t)=jMνανij(ϕj(Nν,t)ϕi(Nν,t)),    iIν, t>0,Diϕix(Nν,t)=jMνανij(ϕj(Nν,t)ϕi(Nν,t)),  iOν, t>0,ανij0, ανij=ανji  for all i,jMν,
    \right.$
    (2.5)

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

    $\sum\limits_{i \in {\mathcal{I}^\nu }} {{D_i}} {\phi _{{i_x}}}({N_\nu },t) = \sum\limits_{i \in {\mathcal{O}^\nu }} {{D_i}} {\phi _{{i_x}}}({N_\nu },t)\;.$

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

    $\left\{λivi(Nν,t)=jMνKνij(uj(Nν,t)ui(Nν,t)), iIν, t>0,λivi(Nν,t)=jMνKνij(uj(Nν,t)ui(Nν,t)),   iOν, t>0,Kνij0, Kνij=Kνji  for all i,jMν .
    \right.\ $
    (2.6)

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

    $ \sum\limits_{i \in {\mathcal{I}^\nu }} {\lambda _i}{v_i}(N_\nu, t) = \sum\limits_{i \in {\mathcal{O}^\nu }}{\lambda _i}{v_i}(N_\nu, t), $

    which corresponds to the conservation of the total mass

    $ \sum\limits_{i \in \mathcal{M}} {\int_{{I_i}} {{u_i}(x,t)\;dx} } = \sum\limits_{i \in \mathcal{M}} {\int_{{I_i}} {{u_{0i}}(x)\;dx} } , $

    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

    $ u_{i0}, v_{i0}, \phi_{i0}\textrm{ satisfy conditions (2.3)-(2.6) for all } i\in\mathcal{M}\ . $ (2.7)

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

    $u, v \in C([0, T];H^1(\mathcal A))\cap C^1([0, T];L^2(\mathcal A)), \phi \in C([0, T];H^2(\mathcal A))\cap C^1([0, T];L^2(\mathcal A))\ $

    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

    $ \textrm{ for all }\nu\in\mathcal P, \textrm{ for some } k\in\mathcal{M}^\nu, K_{ik}^\nu\neq 0 \textrm{ for all } i\in\mathcal{M}^\nu, i\neq k\ . $ (2.8)

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

    $ \mu : = \sum\limits_{i \in \mathcal{M}} {\int_{{I_i}} {{u_i}(x)\;dx} } \ge 0, $

    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

    $\mathcal Q = \{ \nu\in \mathcal P: N_\nu \textrm{ is linked to } N_\mu \textrm{ by a path not covering } I_j\} \ $

    (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

    $ \sum\limits_{\nu\in\mathcal Q\cup \{\mu\}} \left( \sum\limits_{i\in I^\nu} {\lambda _i}{v_i}(N_\nu) -\sum\limits_{i\in O^\nu} {\lambda _i}{v_i}(N_\nu)\right) = 0\ .$

    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

    $v_j(N_\mu) = 0\ ;$

    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

    $\left\{ λiuix=uiϕix,Diϕixx+biϕi=aiui,
    \right.x\in {{I}_{i}},\ i\in \mathcal{M},\ t > 0,\text{ } $
    (3.1)

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

    $ {\phi_{i(j)}}_x(e_j, t) = 0 \;\;\;\; t > 0 , $ (3.2)

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

    $ \sum\limits_{j \in {\mathcal{M}^\nu }} {K_{ij}^\nu ({u_j}({N_\nu }) - {u_i}({N_\nu }))} = 0,\;\;\;\;i \in {\mathcal{M}^\nu }, $ (3.3)
    $ \label{stcf}Diϕix(Nν)=jMνανij(ϕj(Nν)ϕi(Nν)),iIν,Diϕix(Nν)=jMνανij(ϕj(Nν)ϕi(Nν)),iOν.
    $
    (3.4)

    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$,

    $ 0=jMν,jiKνij(uj(Nν)ui(Nν))=jMν,ji,kKνij(uj(Nν)uk(Nν))(jMν,jiKνij)(ui(Nν)uk(Nν)) ;
    $
    (3.5)

    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

    $ u_j(N_\nu) = u_k(N_\nu)\;\;\;\; \textrm{ for all } j\in\mathcal{M}^\nu.$

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

    $ \sum\limits_{i \in \mathcal{M}} {\int_{{I_i}} {{u_i}(x)dx} } = {\mu _0}\;; $ (3.6)

    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

    $\left\{Diϕixx+biϕi=aiuixIi,iM,ui(x)=Ciexp(ϕi(x)λi)xIi,iM,ϕi(j)x(ej)=0,jJ,Diϕix(Nν)=jMνανij(ϕj(Nν)ϕi(Nν)),iIν,νP,Diϕix(Nν)=jMνανij(ϕj(Nν)ϕi(Nν)),iOν,νP,Cjexp(ϕj(Nν)λj)=Ciexp(ϕi(Nν)λi),i,jMν,νP,iMCiIiexp(ϕi(x)λi)dx=μ0.
    \right. $
    (3.7)

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

    $ \label{A2}D(A)={ϕH2(A):(3.2),(3.4) hold },A(ϕ)={Diϕixx+biϕi}iM ;
    $
    (3.8)

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

    $ A \phi = F (\phi, C), $ (3.9)

    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

    $A \phi = F (f, C^f), $ (3.10)

    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:

    $ iMIiDi(ϕi(x)ϕix(x))x dx=νP(iIνDiϕi(Nν)ϕix(Nν)iOνDiϕi(Nν))ϕix(Nν))=νP ijMνανijϕi(Nν)(ϕj(Nν)ϕi(Nν)) =12νP ijMνανij(ϕj(Nν)ϕi(Nν))2 .
    $
    (3.11)

    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

    $ \sum\limits_{i \in \mathcal{M}} C_i^f \int_{{I_i}}\exp{ \left( \frac {f_i(x)} {\lambda_i}\right)} dx = \mu_0 , $ (3.12)

    then

    $\Vert \phi_{x}\Vert_{\infty} \leq \frac {2 \max \{a_i\}_{i\in\mathcal{M}}}{\min \{D_i\}_{i\in\mathcal{M}}} \mu_0 ; $ (3.13)

    if (3.12) holds and

    $ \Vert f_{x}\Vert_{\infty} \leq \frac {2 \max\{a_i\}_{i\in\mathcal{M}} }{\min\{D_i\}_{i\in\mathcal{M}} }\mu_0 , $ (3.14)

    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

    $ \label{12} \Vert \phi \Vert_{W^{2, 1}\;\;(\mathcal{A})}, \Vert \phi \Vert_{H^2(\mathcal{A})} \leq K_{\mu_0}\ . $ (3.15)

    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

    $iMIi(Di(ϕix(x)Γ(ϕi(x)))xDiΓ(ϕi(x))ϕ2ix(x)+biϕi(x)Γ(ϕi(x))Fi(x)Γ(ϕi(x))) dx=0 .
    $

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

    $ iMIiDi(Γ(ϕi)ϕix)x=12νP ijMνανij(ϕj(Nν)ϕi(Nν))(Γ(ϕj(Nν))Γ(ϕi(Nν)))0 ;
    $
    (3.16)

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

    $\sum\limits_{i \in \mathcal{M}} b_i\int_{{I_i}} {\phi _i}(x) \Gamma(-{\phi _i}(x)) dx \geq 0, $

    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

    $\sum\limits_{i \in \mathcal{M}} b_i\int_{{I_i}} {\phi _i}(x) dx = \sum\limits_{i \in \mathcal{M}} \int_{{I_i}} F_i(f(x), C^f) dx $ (3.17)

    which implies

    $ \Vert \phi\Vert_1 \leq \frac {\max\{a_i\}}{\min\{b_i\}} \mu_0\ .$ (3.18)

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

    $|{D_j}{\phi _j}_x(x)| \le \int_{{I_j}} {{D_j}} \left| {{\phi _j}_{yy}(y)} \right|dy \le \int_{{I_j}} {\left( {{b_j}{\phi _j}(y) + C_j^f{a_j}\exp \left( {\frac{{{f_j}(y)}}{{{\lambda _j}}}} \right)} \right)} dy\;.$

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

    $\mathcal Q = \{ \nu\in \mathcal P: N_\nu \textrm{ is linked to } N_\mu \textrm{ by a path not covering } I_j\} , $
    $\mathcal S = \{ i\in \mathcal M: I_i \textrm{ is incident with } N_l \textrm{ for some }l \in\mathcal Q\}\ $

    (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

    $ \sum\limits_{\nu\in\mathcal Q\cup \{\mu\}} \left( \sum\limits_{i\in I^\nu} D_i{\phi _{{i_x}}}(N_\nu) -\sum\limits_{i\in O^\nu} D_i{\phi _{{i_x}}}(N_\nu)\right) = 0\ .$

    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

    $ |Djϕjx(x)|=|Djϕjx(x)+νQ{μ}(iIνDiϕix(Nν)iOνDiϕix(Nν))||iSIiDiϕiyy(y)dy+IxjDjϕjyy(y)dy||iSIi(biϕi(y)Cfiaiexp(fi(y)λi))dy+Ixj(bjϕj(y)Cfjajexp(fj(y)λj))dy|.
    $
    (3.19)

    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

    $ \Vert \phi_x\Vert_1 \leq \frac {2 \max \{a_i\}_{i\in\mathcal{M}}}{\min \{D_i\}_{i\in\mathcal{M}}} |\mathcal A| \mu_0 , $ (3.20)

    and

    $ \Vert \phi_{x}\Vert_2 \leq \frac {2 \max \{a_i\}_{i\in\mathcal{M}}}{\min \{D_i\}_{i\in\mathcal{M}}} |\mathcal{A}|^{\frac 1 2} \mu_0\ ; $ (3.21)

    moreover, by Sobolev embedding theorem, we obtain

    $\Vert \phi\Vert_{\infty} \leq K_1 \mu_0, $ (3.22)

    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

    $ \Vert \phi_{xx}\Vert_1 \leq \frac{2 \max\{a_i\}_{i\in\mathcal{M}}}{\min\{D_i\}_{i\in\mathcal{M}}}\mu_0\ ; $ (3.23)

    then, using (3.11), we have

    $\sum\limits_{i \in \mathcal{M}} \frac {D_i^2}{b_i} \int_{{I_i}}{\phi _{{i_{xx}}}}^2(x) \ dx \leq \sum\limits_{i \in \mathcal{M}} \frac {\Vert F_i\Vert_{\infty}}{b_i} \int_{{I_i}} F_i(x)\ dx \leq \frac { \max \{a_i \Vert F_i\Vert_{\infty}\}}{\min\{b_i\}}\mu_0\ $

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

    $ \sum\limits_{i \in \mathcal{M}} \int_{{I_i}}{\phi _{{i_{xx}}}}^2(x) \ dx \leq K_2 (1+\Vert {f}_x\Vert_\infty) \mu_0^2, $ (3.24)

    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

    $\left(C_i {\exp \left( {\frac{{{\phi _i}(x)}}{{{\lambda _i}}}} \right)}, 0, {\phi _i}(x)\right)\ \ i\in \mathcal{M}, $

    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}$,

    $u_i^0(x) = C^{\phi^0}_i \exp \left(\frac {\phi_i^0(x)} {\lambda_i}\right), $

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

    $ C^{\phi^0}_j \exp \left( {\frac{{\phi _j^0\left( {{N_\nu }} \right)}}{{{\lambda _j}}}} \right) = C^{\phi^0}_i\exp \left( {\frac{{\phi _i^0\left( {{N_\nu }} \right)}}{{{\lambda _i}}}} \right) , \;\;\;\; i, j \in\mathcal M^\nu, \;\;\;\; \nu\in\mathcal P, $ (3.25)
    $\sum\limits_{i \in \mathcal{M}} C^{\phi^0}_i\int_{{I_i}} \exp \left( {\frac{{\phi _i^0\left( x \right)}}{{{\lambda _i}}}} \right) dx = \mu_0\ . $ (3.26)

    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

    $u^0_1(x): = \alpha \exp\left(\frac {\phi_1^0(x)} {\lambda_1}\right)\ .$

    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$,

    $u^0_j(x): = \alpha \exp\left(\frac {\phi_1^0(N_\mu)}{ \lambda_1}\right)\exp\left( -\frac {\phi_j^0(N_\mu)}{\lambda_j}\right)\exp\left(\frac {\phi_j^0(x)} {\lambda_j}\right) \ \textrm{ for all } j\in\mathcal M^\mu, j\neq 1;$

    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

    $ \mathcal E_h(\phi^0) : = \frac {\Pi_{i = 1, ..., h-1} \exp \left( \frac {\phi^0_i(N_i)} {{{\lambda _i}}} \right) } {\Pi_{i = 1, ..., h-1} \exp \left( \frac {\phi^0_{i+1}(N_i)} {\lambda_{i+1}} \right) }, $
    Figure 2. 

    Example: the highlighted arcs form the path from the outer point $e_1$ to the inner node $N_4$ and $I_5$ is an arc incident with $N_4$, not belonghing to the path.

    .

    we define

    $u^0_h(x) : = \alpha\mathcal E_h(\phi^0) \exp\left(\frac { \phi^0_h(x)} {\lambda_h} \right) \ . $

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

    $\alpha\sum\limits_{i \in \mathcal{M}} \mathcal E_i(\phi^0)\int_{{I_i}} \exp \left( {\frac{{\phi _i^0\left( x \right)}}{{{\lambda _i}}}} \right) dx = \mu_0, $

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

    $ u^0_i(x) = C_i^{\phi^0} \exp\left(\frac {\phi_i^0(x)} {\lambda_i}\right), \ \ C_i^{\phi^0}: = \mu_0 \frac {\ \mathcal E_i(\phi^0) } { \ \sum\limits_{j \in \mathcal{M}} \mathcal E_j(\phi^0) \int_{I_j} \exp \left( {\frac{{\phi _j^0\left( x \right)}}{{{\lambda _j}}}} \right) dx \ }. $ (3.27)

    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}$,

    $A\phi^1 = F(C^{\phi^0}, \phi^0)\ ;$

    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

    $B_{\mu_0}: = \{ \phi\in D(A) : \phi\geq 0, \Vert \phi_{x}\Vert_{\infty} \leq \frac {2 \max \{a_i\}_{i\in\mathcal{M}}}{\min \{D_i\}_{i\in\mathcal{M}}} \mu_0, \Vert \phi \Vert_{H^2} \leq K_{ \mu_0}\ \}$

    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

    $ biIi(ϕ1i(x)¯ϕ1i(x))2,dx+DiIi(ϕ1ix(x)¯ϕ1ix(x))2,dxDiIi((ϕ1ix(x)¯ϕ1ix(x))(ϕ1i(x)¯ϕ1i(x)))xdx=aiIi(u0i(x)¯u0i(x))(ϕ1i(x)¯ϕ1i(x))dx ;
    $
    (3.28)

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

    $\sum\limits_{i \in \mathcal{M}} \Vert {{\phi _i}}^1 -\overline {{\phi _i}}^1\Vert_{H^2} \leq K(a_i, b_i, D_i) \sum\limits_{i \in \mathcal{M}} \Vert u^0_i -\overline {u}_i^0\Vert_{2} , $ (3.29)

    We set

    $ J_i^{\phi^0}: = \int_{I_i}\exp \left( {\frac{{\phi _i^0\left( x \right)}}{{{\lambda _i}}}} \right) dx, \;\;\;\;E_i^{\phi^0}(x): = \exp \left( {\frac{{\phi _i^0\left( x \right)}}{{{\lambda _i}}}} \right)\ \ ;$

    we have

    $ \left| {u_i^0(x) - \bar u_i^0(x)} \right| = {\mu _0}\left| {\frac{{\;{\varepsilon _i}({\phi ^0})E_i^{{\phi ^0}}(x)\;}}{{\;\sum\limits_{j \in \mathcal{M}} {{\varepsilon _j}({\phi ^0})J_j^{{\phi ^0}}} \;}} - \frac{{\;{\varepsilon _i}({{\bar \phi }^0})E_i^{{{\bar \phi }^0}}(x)\;}}{{\;\sum\limits_{j \in \mathcal{M}} {{\varepsilon _j}({{\bar \phi }^0})J_j^{{{\bar \phi }^0}}} \;}}} \right|. $ (3.30)

    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}$

    $ \mathop {\max }\limits_{{I_i}} E_i^g(x) \le {K_6},\;\;\;\;J_i^g \le {K_6}\left| {{I_i}} \right|, $
    $ \left| {E_i^{{\phi ^0}}(x) - E_i^{{{\bar \phi }^0}}(x)} \right| \le {K_6}\left| {\phi _i^0(x) - \bar \phi _i^0(x)} \right|, $
    $ \left| {J_i^{{\phi ^0}} - J_i^{{{\bar \phi }^0}}} \right| \le {K_6}\int_{{I_i}} {\left| {\phi _i^0(x) - \bar \phi _i^0(x)} \right|{dx}} . $

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

    $ \sum\limits_{i \in \mathcal{M}} \Vert {{\phi _i}}^1 -\overline {{\phi _i}}^1\Vert_{H^2} \leq \mu_0 K_7(a_i, b_i, D_i, K_{\mu_0}, |\mathcal{A}|) \sum\limits_{i \in \mathcal{M}} \Vert \phi_i^0 -\overline {\phi}_i^0\Vert_{H^1} , $ (3.31)

    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

    $ \sum\limits_{\nu\in\mathcal P}\left( \sum\limits_{i\in\mathcal I^\nu} {\lambda _i} {u_i}(N^\nu){v_i}(N^\nu) - \sum\limits_{i\in\mathcal O^\nu} {\lambda _i} {u_i}(N^\nu){v_i}(N^\nu)\right)\geq 0, $

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

    $2 \sum\limits_{i \in \mathcal{M}} \beta_i \int_{{I_i}} v_i^2(x)dx \leq \sum\limits_{i \in \mathcal{M}} \Vert u_i\Vert_{\infty} \int_{{I_i}} (v_i^2(x) +\phi^2_{i_x}(x) ) \ dx $

    and

    $ \sum\limits_{i \in \mathcal{M}} {\lambda _i}\int_{{I_i}} u^2_{i_x}(x) \ dx \leq \sum\limits_{i \in \mathcal{M}} \Vert u_i\Vert_{\infty} \int_{{I_i}} (u_{i_x}^2(x) +\phi^2_{i_x}(x) ) \ dx \ + \sum\limits_{i \in \mathcal{M}} \frac {\beta_i^2}{{\lambda _i}} \int_{{I_i}} v^2_i(x)\ dx \ ;$

    the above inequalities implies the following one

    $ \Vert v\Vert_2^2 + \Vert u_x\Vert_2^2\leq K_0 H \left( \Vert \phi_x\Vert_2^2 + \Vert v\Vert_2^2 +\Vert u_x\Vert_2^2 \right), $ (3.32)

    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

    $ -\sum\limits_{\nu\in\mathcal P}\left( \sum\limits_{i\in\mathcal I^\nu} D_i{\phi _i}(N^\nu){\phi _{{i_x}}}(N^\nu) - \sum\limits_{i\in\mathcal O^\nu} D_i {\phi _i}(N^\nu){\phi _{{i_x}}}(N^\nu)\right)\geq 0\ ;$

    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$,

    $ u_j(N_\nu) = u_k(N_\nu)+\sum\limits_{i\in \mathcal{M}^\nu, i\neq k} \theta^\nu_{ij} {v_i}(N_\nu) \;\;\;\;\textrm{ for all } j\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

    $ \Vert \phi_x\Vert_2^2 + \Vert \phi_{xx}\Vert_2^2\leq K_1 \left( \Vert v\Vert_2^2 +\Vert u_x\Vert_2^2 \right), $ (3.33)

    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

    $ \Vert v\Vert_2^2 + \Vert u_x\Vert_2^2\leq K_0(1+K_1) H \left( \Vert v\Vert_2^2 +\Vert u_x\Vert_2^2 \right), $

    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,

    $ {u_i}(x) = \frac {\mu_0}{|\mathcal{A}|} , \ \ \ \phi_i(x) = \frac{a_i}{b_i} \frac{\mu_0}{|\mathcal{A}|}, \;\;\;i\in\mathcal{M}\ .$

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

    $ \sum\limits_{j\in\mathcal{M}^\nu}\alpha^\nu_{ij} \frac{\mu_0}{|\mathcal{A}|} \left( \frac {a_j} {b_j}-\frac {a_i} {b_i}\right) = 0 , \;\;\;i\in\mathcal{M}^\nu , $

    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,

    $ \alpha_{11} \frac{\mu_0}{|\mathcal{A}|} \left( \frac {b_2} {a_2}-\frac {b_1} {a_1}\right) = 0 , $

    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.

    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

    $ u,vC([0,+);H1(A))C1([0,+);L2(A)),ϕC([0,+);H2(A))C1([0,+);L2(A)),ϕxH1(A×(0,+)),
    $
    (4.1)

    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

    $\label {Q} \frac {a_i}{b_i} = Q\ \ \textrm{for all }i\in \mathcal{M}\ ;$ (4.2)

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

    $F2T(u,v,ϕ):=iM(supt[0,T]ui(t)2H1+supt[0,T]vi(t)2H1+supt[0,T]ϕix(t)2H1)+T0(ux(t)22+v(t)2H1+vt(t)22+ϕx(t)2H1+ϕxt(t)22) dt
    $
    (4.3)

    is uniformly bounded for $T>0$.

    Here and below we use the notations

    $\Vert f_i(t)\Vert_2: = \Vert f_i(\cdot, t)\Vert_{L^2(I_i)}, \ \ \Vert f_i(t)\Vert_{H^s}: = \Vert f_i(\cdot, t)\Vert_{H^s(I_i)} \ . $

    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

    $ \left\{ tui+λixvi=0tvi+λixui=(ui+¯u)xϕiβivixIi,t0,iM,tϕi=Dixxϕi+aiuibiϕi,
    \right.$
    (4.4)

    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),

    $ u,vC([0,T];H1(A))C1([0,T];L2(A)),ϕC([0,T];H2(A))C1([0,T];L2(A)), ϕxH1(A.×(0,T)) ;
    $

    then

    a)

    $ iM(sup[0,T]ui(t)22+sup[0,T]vi(t)22+βiT0vi(t)22dt)CiM(u0i22+v0i22)+CiM(sup[0,T]ui(t)H1+ˉu)To(ϕix(t)22+vi(t)22)dt;
    $

    b)

    $iM(sup[0,T]vix(t)22+sup[0,T]vit(t)22+T0vit(t)22,dt)C(v02H1+u02H1ϕ02H2)+CiM(sup[0,T]ui(t)H1+ˉu)T0(ϕixt(t)22+vit(t)22) dt+CiMsup[0,T]ϕx(t)H1T0(vit(t)22+vi(t)2H1)dt ;
    $

    c)

    $iMsup[0,T]uix(t)22CiM(sup[0,T]vit(t)22+sup[0,T]vi(t)22)+CiM(sup[0,T]ui(t)H1+ˉu)(sup[0,T]uix(t)22+sup[0,T]ϕix(t)22) ;
    $

    d)

    $iMT0uix(t)22,dtCiMT0(vit(t)22+vi(t)22)dt+CiM(sup[0,T]ui(t)H1+ˉu)T0(uix(t)22+ϕix(t)22)dt ;
    $

    e)

    $iMT0vix(t)22,dtCiM(v0i22+u0i2H1(1+ϕ0i2H1))+CiM(T0vit(t)22,dt+sup[0,T]vit(t)22)+CiM(sup[0,T]ui(t)H1+sup[0,T]ϕix(t)H1+ˉu)×T0(vi(t)2H1+ϕixt(t)22)dt ;
    $

    f)

    $iM(sup[0,T]ϕit(t)22+T0(ϕit(t)22+ϕitx(t)22) dt)CiM(ϕ0i2H2+u0i22+T0uit(t)22) ;
    $

    g)

    $iM(sup[0,T]ϕixx(t)22+sup[0,T]ϕix(t)22)CiM(sup[0,T]ϕit(t)22+sup[0,T]ui(t)22) ;
    \ $

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

    $iMT0(ϕix(t)22+ϕixx(t)22) dtCiMT0(uix(t)22+vi(t)2H1+ϕit(t)22) dt,
    $

    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

    $\overline {u}\leq \epsilon_1, \ \ \Vert u_0\Vert_{H^1}, \Vert v_0\Vert_{H^1}, \Vert \phi_{0}\Vert_{H^2}\leq \epsilon_0 , $

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

    $u, v \in C([0, +\infty);H^1(\mathcal{A}))\cap C^1([0, +\infty);L^2(\mathcal{A})) , $
    $\phi\in C([0, +\infty);H^2(\mathcal{A}))\cap C^1([0, +\infty);L^2(\mathcal{A})), \ \phi_{x}\in H^1\left( \mathcal A\times (0, +\infty)\right) \ .$

    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

    $F_T^2(u, v, \phi)\leq c_1 F_0^2(u, v, \phi) +c_2 \overline {u} F^2_T(u, v, \phi) +c_3 F_T^3(u, v, \phi), $

    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

    $F_T^2(u, v, \phi)\leq c_4 F_0^2(u, v, \phi) +c_5 F_T^3(u, v, \phi)$

    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

    $\overline {u} \leq \epsilon_2, \;\;\;\sum\limits_{i \in \mathcal{M}}\int_{{I_i}} u_0(x) = \overline \mu, \;\;\;\Vert( u_0-\overline {u}, v_0, \phi_0-\overline {\phi})\Vert_{(H^1)^2\times H^2} \leq \epsilon_0, $

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

    $u, v \in C([0, +\infty);H^1(\mathcal{A}))\cap C^1([0, +\infty);L^2(\mathcal{A})) , $
    $\phi\in C([0, +\infty);H^2(\mathcal{A}))\cap C^1([0, +\infty);L^2(\mathcal{A})), $

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

    $ \lim\limits_{t\to+\infty}\Vert {u_i}(\cdot, t)-\overline {u}\Vert_{C(\overline I_i)}, \lim\limits_{t\to+\infty}\Vert {v_i}(\cdot, t)\Vert_{C(\overline I_i)}, \lim\limits_{t\to+\infty}\Vert {\phi _i}(\cdot, t)-\overline {\phi}\Vert_{C^1(\overline I_i)} = 0\ . $

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

    $ \hat u: = u - \bar u,\;\hat v: = v,\;\hat \phi : = \phi - \bar \phi ; $

    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$,

    $limn+iMˆui(,tn)Ui()C(¯Ii)=0,limn+iMˆvi(,tn)Vi()C(¯Ii)=0,limn+iMˆϕi(,tn)Φi()C1(¯Ii)=0 .
    $
    (4.5)

    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

    $\omega_i(t): = \Vert \hat v_i(t, \cdot)\Vert_{L^2(I_i)}$

    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

    $ V_i(x) = 0, \;\;\;\; a_i U_i(x)- b_i \Phi_i(x) = 0, \;\;\;\; \Phi_i(x) = \overline \Phi_i, \;\;\;\; \ x\in I_i, $

    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.

    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.

    [1] Stability results for a diffusion equation with functional drift approximating a chemotaxis model. Trans. Amer. Math. Soc. (1987) 300: 235-258.
    [2] The scalar Keller-Segel model on networks. Math. Models Methods Appl. Sci. (2014) 24: 221-247.
    [3] A hyperbolic model of chemotaxis on a network: A numerical study. ESAIM: Mathematical Modelling and Numerical Analysis (2014) 48: 231-258.
    [4]

    T. Cazenave and A. Haraux, An Introduction to Semilinear Evolution Equations, Clarendon Press-Oxford, 1998.

    [5] Traffic flow on a road network. SIAM J. Math. Anal. (2005) 36: 1862-1886.
    [6] Parabolic models for chemotaxis on weighted nerworks. J. Math. Pures Appl. (2017) 108: 459-480.
    [7]

    R. Dager and E. Zuazua, Wave Propagation, Observation and Control in 1-d Flexible Multi-Structures, Vol. 50 of Mathematiques & Applications [Mathematics & Applications] Springer-Verlag, Berlin, 2006.

    [8] Cattaneo models for chemosensitive movement. Numerical solution and pattern formation. J. Math. Biol. (2003) 46: 153-170.
    [9] Derivation of hyperbolic model for chemosensitive movement. J. Math. Biol. (2005) 50: 189-207.
    [10]

    M. Garavello and B. Piccoli, Traffic Flow on Networks -Conservation Laws Models, AIMS Series on Applied Mathematics, Vol. 1, American Institute of Mathematical Sciences (AIMS), Springfield, MO, 2006.

    [11] Stability of constant states and qualitative behavior of solutions to a one dimensional hyperbolic model of chemotaxis. Discrete Contin. Dyn. Syst. Ser. B (2009) 12: 39-76.
    [12] Global smooth solutions for a hyperbolic chemotaxis model on a network. SIAM J. Math. Anal. (2015) 47: 4652-4671.
    [13] Microarchitecture of three-dimensional scaffold influences cell migration behavior via junction interaction. Biophysical Journal (2008) 29: 4013-4024.
    [14] Hyperbolic models for chemosensitive mevement. Math. Models Methods Appl. Sci. (2002) 12: 1007-1034.
    [15] Existence of weak solutions for a hyperbolic model of chemosensitive movement. J. Math. Anal. Appl. (2001) 260: 173-199.
    [16] Hyperbolic model for chemotaxis in 1-D. Nonlinear Anal. Real World Appl. (2000) 1: 409-433.
    [17] From 1970 until present: The Keller-Segel model in chemotaxis and its consequences. I.Jahresber.Deutsch Math-Verein (2003) 105: 103-165.
    [18] Initiation of slime mold aggregation viewed as an instability. J. Theor. Biol. (1970) 26: 399-415.
    [19] Spectral properties and asymptotic periodicity of flows in networks. Mathematische Zeitschrift (2005) 249: 139-162.
    [20] Cell proliferation and migration in silk broin 3D scaffolds. Biomaterials (2009) 30: 2956-2965.
    [21]

    D. Mugnolo, Simigroup Methods for Evolutions Equations on Networks, Springer, Berlin, 2014.

    [22]

    J. D. Murray, Mathematical Biology. I An Introduction, Third edition. Interdisciplinary Applied Mathematics, 17 Springer Verlag, New York, 2002; Mathematical Biology. Ⅱ Spatial models and biomedical applications, Third edition. Interdisciplinary Applied Mathematics, 18 Springer Verlag, New York, 2003.

    [23]

    B. Perthame, Transport Equations in Biology, Frontiers in Mathematics, Birkhauser, 2007.

    [24] A theoretical study of receptor mechanisms in bacterial chemotaxis. SIAM J.Appl.Math. (1977) 32: 653-665.
    [25] A G-CSF functionalized PLLA scaffold for wound repair: An in vitro preliminary study. Conf. Proc. IEEE Eng.Med.Biol.Soc. (2010) 843-846.
    [26] Stabilization of the wave equation on 1-D networks. SIAM J. Control Optim. (2009) 48: 2771-2797.
  • This article has been cited by:

    1. Francesca Romana Guarguaglini, Marco Papi, Flavia Smarrazzo, Local and global solutions for a hyperbolic–elliptic model of chemotaxis on a network, 2019, 29, 0218-2025, 1465, 10.1142/S021820251950026X
    2. Yuta Ishii, Multi-peak solutions for the Schnakenberg model with heterogeneity on star shaped graphs, 2023, 446, 01672789, 133679, 10.1016/j.physd.2023.133679
    3. Francesca R. Guarguaglini, Roberto Natalini, Vanishing viscosity approximation for linear transport equations on finite star-shaped networks, 2021, 21, 1424-3199, 2413, 10.1007/s00028-021-00688-0
    4. Yafeng Li, Chunlai Mu, Qiao Xin, Global existence and asymptotic behavior of solutions for a hyperbolic‐parabolic model of chemotaxis on network, 2022, 45, 0170-4214, 6739, 10.1002/mma.8204
  • Reader Comments
  • © 2018 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5622) PDF downloads(202) Cited by(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog