Loading [MathJax]/jax/output/SVG/jax.js

A class of infinite horizon mean field games on networks

  • We consider stochastic mean field games for which the state space is a network. In the ergodic case, they are described by a system coupling a Hamilton-Jacobi-Bellman equation and a Fokker-Planck equation, whose unknowns are the invariant measure $ m $, a value function $ u $, and the ergodic constant $ \rho $. The function $ u $ is continuous and satisfies general Kirchhoff conditions at the vertices. The invariant measure $ m $ satisfies dual transmission conditions: in particular, $ m $ is discontinuous across the vertices in general, and the values of $ m $ on each side of the vertices satisfy special compatibility conditions. Existence and uniqueness are proven under suitable assumptions.

    Citation: Yves Achdou, Manh-Khang Dao, Olivier Ley, Nicoletta Tchou. A class of infinite horizon mean field games on networks[J]. Networks and Heterogeneous Media, 2019, 14(3): 537-566. doi: 10.3934/nhm.2019021

    Related Papers:

    [1] T. A. Shaposhnikova, M. N. Zubova . Homogenization problem for a parabolic variational inequality with constraints on subsets situated on the boundary of the domain. Networks and Heterogeneous Media, 2008, 3(3): 675-689. doi: 10.3934/nhm.2008.3.675
    [2] Ken-Ichi Nakamura, Toshiko Ogiwara . Periodically growing solutions in a class of strongly monotone semiflows. Networks and Heterogeneous Media, 2012, 7(4): 881-891. doi: 10.3934/nhm.2012.7.881
    [3] Junlong Chen, Yanbin Tang . Homogenization of nonlinear nonlocal diffusion equation with periodic and stationary structure. Networks and Heterogeneous Media, 2023, 18(3): 1118-1177. doi: 10.3934/nhm.2023049
    [4] Benjamin Contri . Fisher-KPP equations and applications to a model in medical sciences. Networks and Heterogeneous Media, 2018, 13(1): 119-153. doi: 10.3934/nhm.2018006
    [5] Thomas Geert de Jong, Georg Prokert, Alef Edou Sterk . Reaction–diffusion transport into core-shell geometry: Well-posedness and stability of stationary solutions. Networks and Heterogeneous Media, 2025, 20(1): 1-14. doi: 10.3934/nhm.2025001
    [6] Iryna Pankratova, Andrey Piatnitski . Homogenization of convection-diffusion equation in infinite cylinder. Networks and Heterogeneous Media, 2011, 6(1): 111-126. doi: 10.3934/nhm.2011.6.111
    [7] Feiyang Peng, Yanbin Tang . Inverse problem of determining diffusion matrix between different structures for time fractional diffusion equation. Networks and Heterogeneous Media, 2024, 19(1): 291-304. doi: 10.3934/nhm.2024013
    [8] Xavier Blanc, Claude Le Bris . Improving on computation of homogenized coefficients in the periodic and quasi-periodic settings. Networks and Heterogeneous Media, 2010, 5(1): 1-29. doi: 10.3934/nhm.2010.5.1
    [9] Avner Friedman . PDE problems arising in mathematical biology. Networks and Heterogeneous Media, 2012, 7(4): 691-703. doi: 10.3934/nhm.2012.7.691
    [10] Bendong Lou . Self-similar solutions in a sector for a quasilinear parabolic equation. Networks and Heterogeneous Media, 2012, 7(4): 857-879. doi: 10.3934/nhm.2012.7.857
  • We consider stochastic mean field games for which the state space is a network. In the ergodic case, they are described by a system coupling a Hamilton-Jacobi-Bellman equation and a Fokker-Planck equation, whose unknowns are the invariant measure $ m $, a value function $ u $, and the ergodic constant $ \rho $. The function $ u $ is continuous and satisfies general Kirchhoff conditions at the vertices. The invariant measure $ m $ satisfies dual transmission conditions: in particular, $ m $ is discontinuous across the vertices in general, and the values of $ m $ on each side of the vertices satisfy special compatibility conditions. Existence and uniqueness are proven under suitable assumptions.



    In this paper, we study the positive solutions of the periodic-parabolic problem

    $ {ut=μk(x,t)Δu+m(x,t)uc(x,t)up, in Ω×R,uν=0, on Ω×R,u(x,0)=u(x,T), in Ω, $ (1.1)

    where $ \Omega $ is a bounded domain of $ \mathbb{R}^N(N\geq1) $ with smooth boundary $ \partial\Omega $, $ \nu $ is the outward normal vector of $ \partial\Omega $, $ \mu > 0 $ and $ p > 1 $ is constant, $ m(x, t)\in C^{\alpha, \frac{\alpha}{2}}(\bar\Omega\times\mathbb{R})\; (0 < \alpha < 1) $ is $ T $-periodic in $ t $, $ k(x, t) $, $ c(x, t)\in C^{\alpha, 1}(\bar\Omega\times\mathbb{R}) $ are positive and $ T $-periodic in $ t $. It is known that the periodic reaction-diffusion equation (1.1) can be accurately used to describe different diffusion phenomena in infectious diseases, microbial growth, and population ecology, see [1,2,3,4]. From a biological point of view, $ \Omega $ represents the habitat of species $ u $ and $ \mu k(x, t) $ stands for the diffusion rate, which is time and space dependent. The function $ m(x, t) $ represents the growth rate of species. In this situation, in the subset $ \{(x, t)\in\Omega\times\mathbb{R}:m(x, t) > 0\} $, the species will increase, while in $ \{(x, t)\in\Omega\times\mathbb{R}:m(x, t) < 0\} $, species will decrease. The coefficient $ c(x, t) $ means that environment $ \Omega $ can accommodate species $ u $. There are many interesting conclusions about the study of the reaction-diffusion equation, see [5,6,7,8] for the elliptic problems and [9,10,11,12,13,14] for the periodic problems.

    In particular, if $ k(x, t)\equiv k(t) $ for $ x\in\bar\Omega $, problem (1.1) has been well investigated by Hess [2], Cantrell and Cosner [1]. Let $ \lambda(\mu) $ be the unique principal eigenvalue of the eigenvalue problem

    $ {utμk(t)Δum(x,t)u=λ(μ)u, in Ω×R,uν=0, on Ω×R,u(x,0)=u(x,T), in Ω. $

    It follows from [2,15] that Eq (1.1) has a positive periodic solution $ \theta_\mu(x, t) $ if and only if $ \lambda(\mu) < 0 $. In addition, Dancer and Hess [16] and Daners and López-Gómez [17] studied the effect of $ \mu $ on the positive periodic solution of Eq (1.1) with various boundary conditions. The most interesting conclusion of [11,16,17] is that

    $ limμ0+θμ(x,t)=θ(x,t) locally uniformly  in Ω×[0,T], $

    here $ \theta(x, t) $ is the maximum nonnegative periodic solution of

    $ {ut=m(x,t)uc(x,t)up,tR,u(x,0)=u(x,T). $

    However, there is little result on the associated large diffusion and the effect of large diffusion on positive solutions.

    Our goal is to study the existence and uniqueness of positive periodic solutions of Eq (1.1) and the asymptotic behavior of positive periodic solutions when the diffusion rate $ \mu $ is large. To this end, let $ \lambda(\mu; m) $ be the principal eigenvalue of

    $ {utμk(x,t)Δum(x,t)u=λ(μ;m)u, in Ω×R,uν=0, on Ω×R,u(x,0)=u(x,T), in Ω. $ (1.2)

    {It is well known that $ \lambda(\mu; m) $ plays a major role in the study of the positive periodic solution of Eq (1.1). The properties of $ \lambda(\mu; m) $ will be established in Section 2. In addition, let $ W^{2, p}_{\nu}(\Omega) = \{u\in W^{2, p}(\Omega):\frac{\partial u}{\partial\nu} = 0\}(N < p < \infty) $. If $ u_0\in W^{2, p}_{\nu}(\Omega) $, it follows from [2] that the semilinear initial value problem}

    $ {ut=μk(x,t)Δu+m(x,t)uc(x,t)up, in Ω×R,uν=0, on Ω×R,u(x,0)=u0(x), in Ω, $

    has a unique solution $ U(x, t) = U(x, t;u_0) $ satisfying

    $ U(x,t)C1+α,1+α2(ˉΩ×[0,T])C2+α,1+α2(ˉΩ×(0,T]). $

    Our first result is the existence and uniqueness of positive periodic solutions of Eq (1.1). For simplicity, in the rest of this paper, we use the following notations:

    $ k(t)=Ω1k(x,t)dx,m(t)=Ωm(x,t)k(x,t)dx,c(t)=Ωc(x,t)k(x,t)dx. $

    Theorem 1.1. Suppose that $ \int_0^T m_*(t)\, dt > 0 $. Then Eq (1.1) admits a unique positive periodic solution $ \theta_\mu(x, t) $ for all $ \mu > 0 $.

    Remark 1.1. By the result of Section 3, we know that there exists a unique positive solution to Eq (1.1) if and only if $ \lambda(\mu; m) < 0 $. In the case $ \int_0^Tm_*(t)\, dt > 0 $, we obtain that $ \lambda(\mu; m) < 0 $.

    Next, we study the asymptotic behavior of positive periodic solutions when the diffusion rate is large.

    Theorem 1.2. Suppose that

    $ m(t)>0 for t[0,T]. $ (1.3)

    Let $ \theta_\mu(x, t) $ be the unique positive periodic solution of Eq (1.1) for $ \mu > 0 $. Then we have

    $ limμθμ(x,t)=ω(t) in C1,12(ˉΩ×[0,T]), $ (1.4)

    where $ \omega(t) $ is the unique positive periodic solution of

    $ {ut=m(t)k(t)uc(t)k(t)up,tR,u(0)=u(T). $ (1.5)

    Remark 1.2. With the approach of local upper-lower solutions developed by Daners and López-Gómez [17] in the study of classical periodic-parabolic logistic equations, we can prove that

    $ limμ0+θμ(x,t)=θ(x,t) locally uniformly  in Ω×[0,T], $

    provided $ \max_{\Omega}\int_0^Tm(x, t)\, dt > 0 $. It also shows that when $ m_*(t) < 0 < \int_0^T\max_{\Omega}m(x, t)\, dt $, populations with small dispersal rates survive, while populations with large dispersal rates perish. This means that a small diffusion rate is a better strategy than a large diffusion rate under appropriate circumstances.

    The rest of this paper is arranged as follows: In Section 2, we study the properties of principal eigenvalues for the periodic eigenvalue problems. In Section 3, we mainly study the existence, uniqueness and stability of the positive solution to Eq (1.1). Moreover, we investigate the asymptotic profiles of the positive periodic solution to Eq (1.1) as $ \mu\to\infty $ in Section 4.

    In this section, we consider the principal eigenvalue of Eq (1.2). To this end, we first study the linear initial value problem

    $ {utμk(x,t)Δua(x,t)u=0, in Ω×(τ,T],uν=0, on Ω×(τ,T],u(x,τ)=u0(x), in Ω, $ (2.1)

    where $ 0\leq\tau < T $, $ u_0\in W^{2, p}_{\nu}(\Omega)(N < p < \infty) $ and $ a(x, t)\in C^{\alpha, \frac{\alpha}{2}}(\bar\Omega\times[\tau, \infty)) $. It is well known that there is a one-to-one correspondence between Eq (2.1) and the evolution operator $ U_\mu(t, \tau) $. Then we can define that $ u(x, t) = U_\mu(t, \tau)u_0 $ is the solution of Eq (2.1). For simplicity, let $ X = L^p(\Omega)(N < p < \infty) $, $ X_1 = W^{2, p}_{\nu}(\Omega) $ and

    $ F={uCα,α2(ˉΩ×R):u(,t+T)=u(,t) in R}. $

    Inspired by the classical works of Hess [2], we first give some important results of Eq (2.1), which will be used in the rest of this paper.

    Lemma 2.1. If $ u_0\in X $ is positive, then $ U_\mu(t, \tau)u_0 > 0 $ in $ C^1_{\nu}(\bar\Omega) $ for $ 0\le\tau < t\leq T $.

    Proof. Note that $ X_1 $ is compactly embedded in $ X $. The operator $ U_\mu(t, \tau)/_{X_1}:X_1\to X_1 $ can be continuously extended to the positive operator $ U_\mu(t, \tau)\in\mathcal{L}(X, X_1) $. Thus $ U_\mu(t, \tau)u_0\geq0 $ in $ X_1 $. Since $ s\mapsto U_\mu(s, \tau)u_0 $ is continuous from $ [\tau, T] $ to $ X_1 $ and $ U_\mu(\tau, \tau)u_0 = u_0\neq0 $, we can get that $ U_\mu(s, \tau)u_0 > 0 $ in $ X_1 $ as $ s > \tau $ goes to $ \tau $. In addition, we have

    $ Uμ(t,τ)u0=Uμ(t,s)Uμ(s,τ)u0, $

    for $ \tau < s < t $. Thus it can be obtained that $ U_\mu(t, \tau)u_0 > 0 $ for $ 0\le\tau < t\leq T $.

    We now study the periodic-parabolic eigenvalue problem

    $ {utμk(x,t)Δua(x,t)u=λ(μ;a)u, in Ω×(0,T],uν=0, on Ω×(0,T],u(x,0)=u(x,T), in Ω. $ (2.2)

    If there is a nontrivial solution $ u(x, t) $ of Eq (2.2), then $ \lambda(\mu; a) $ is called the eigenvalue. In particular, if $ u(x, t) $ is positive, then $ \lambda(\mu; a) $ is the principal eigenvalue.

    Theorem 2.1. Let $ K_\mu: = U_\mu(T, 0) $ and $ r $ be the spectral radius of $ K_\mu $. Then $ r $ is the principal eigenvalue of $ K_\mu $ with positive eigenfunction $ u_0 $ if and only if $ \lambda(\mu; a) = -\frac{1}{T}\ln r $ is the principal eigenvalue of Eq (2.2) with positive eigenfunction $ u(x, t) = e^{\lambda(\mu; a) t}U_\mu(t, 0)u_0 $.

    Proof. It can be proved by the similar arguments as in [2, Proposition 14.4]. For the completeness, we provide a proof in the following. Suppose that $ r $ is the principal eigenvalue of $ K_\mu $ with positive eigenfunction $ u_0\in X_1 $. Let $ u(x, t) = e^{\lambda(\mu; a) t}U_\mu(t, 0)u_0 $. Then $ u(x, t) $ satisfies

    $ {utμk(x,t)Δua(x,t)u=λ(μ;a)u, in Ω×(0,T],uν=0, on Ω×(0,T],u(x,0)=u0=1rKμu0=eλ(μ;a)TKμu0=u(x,T), in Ω. $

    According to the regularity results, we have $ u(x, t)\in C^{2+\alpha, 1+\frac{\alpha}{2}}(\bar\Omega\times\mathbb{R}) $. This means that $ \mu = -\frac{1}{T}\ln r $ is the principal eigenvalue of Eq (2.2), while $ u(x, t) = e^{\lambda(\mu; a) t}U_\mu(t, 0)u_0 $ is the corresponding positive eigenfunction.

    On the contrary, suppose that $ \lambda(\mu; a) = -\frac{1}{T}\ln r $ is the eigenvalue of Eq (2.2) with positive eigenfunction $ u(x, t) $. Set $ v(x, t) = e^{-\lambda(\mu; a) t}u(x, t) $. Then $ v(x, t) $ is the solution of

    $ {vtμk(x,t)Δva(x,t)v=0, in Ω×(0,T],vν=0, on Ω×(0,T],v(x,0)=u(x)=:u0, in Ω. $

    Thus, we obtain $ v(x, t) = U_\mu(t, 0)u_0 $ for $ 0\leq t\leq T $ and $ u_0\in X_1 $ is positive. Hence,

    $ v(T)=eλ(μ;a)Tu0=Kμu0. $

    It follows from Krein-Rutman theorem that $ e^{-\lambda(\mu; a)T} = r $.

    Remark 2.1. For $ \tau < t $, it follows that $ U_\mu(t, \tau) $ is a compact and strongly positive operator on $ X_1 $. Moreover, by Krein-Rutman theorem, we obtain $ r > 0 $, and $ r $ is the unique principal eigenvalue of $ K_\mu $. This implies that Eq (2.2) has the unique principal eigenvalue $ \lambda(\mu; a) $ for any $ \mu > 0 $.

    Lemma 2.2. Let $ a_1(x, t) $, $ a(x, t)\in F $ satisfy

    $ a1(x,t)<a2(x,t) in ˉΩ×[0,T]. $

    Then $ \lambda(\mu; a_2) < \lambda(\mu; a_1) $ for any $ \mu > 0 $.

    Proof. Assume that there exists $ \mu_1 > 0 $ such that $ \lambda(\mu_1;a_2)\geq\lambda(\mu_1;a_1) $. Let $ u_1(x, t) $ and $ u_2(x, t) $ be corresponding positive eigenfunctions, chosen in such a way that

    $ 0<u1(x,t)<u2(x,t) in ˉΩ×[0,T]. $

    Then $ \omega(x, t) = u_2(x, t)-u_1(x, t) $ satisfies

    $ {ωtμ1k(x,t)Δωa1(x,t)ω>λ(μ1;a1)ω, in Ω×(0,T],ων=0, on Ω×(0,T],ω(x,0)=ω(x,T), in Ω. $

    Set $ \phi(x, t) = e^{-\lambda(\mu_1;a_1)t}\omega(x, t) $, then we have

    $ {ϕtμ1k(x,t)Δϕa1(x,t)ϕ>0, in Ω×(0,T],ϕν=0, on Ω×(0,T],ϕ(x,0)=ω(x,0)=ω(x,T), in Ω. $

    Thus, for any $ x\in\Omega $, we can obtain

    $ ϕ(x,T)>Kμ1ω(x,0) and ϕ(x,T)=eλ(μ1;a1)Tω(x,0). $

    Hence, we obtain

    $ (eλ(μ1;a1)TKμ1)ω(x,0)>0 in X1. $

    Note that $ \omega(x, 0) > 0 $, it follows from [2, Theorem 7.3] that

    $ eλ(μ1;a1)T=rμ1<eλ(μ1;a1)T, $

    where $ r_{\mu_1} $ is the principal eigenvalue of $ K_{\mu_1} $. This is a contradiction.

    Lemma 2.3. Suppose that for any $ n\in\mathbb{N} $, $ a_n(x, t)\in F $ satisfies

    $ limnan(x,t)=a(x,t) in C1(ˉΩ×[0,T]). $

    Then for fixed $ \mu > 0 $, we have

    $ limnλ(μ;an)=λ(μ;a). $

    Proof. For any given $ \varepsilon > 0 $, there exists $ n_\varepsilon\in\mathbb{N} $ such that for any $ n > n_\varepsilon $, there holds

    $ a(x,t)ε<an(x,t)<a(x,t)+ε in ˉΩ×[0,T]. $

    Notice that $ \lambda(\mu; a\pm\varepsilon) = \lambda(\mu; a)\mp\varepsilon $. From Lemma 2.2, we have

    $ λ(μ;a)ε<λ(μ;an)<λ(μ;a)+ε, $

    for any $ n > n_\varepsilon $.

    Lemma 2.4. Let $ \lambda(\mu; a) $ be the principal eigenvalue of Eq (2.2) for $ \mu > 0 $. Then we have

    $ λ(μ;a)T0a(t)dtT0k(t)dt, $ (2.3)

    here $ a_*(t) = \int_{\Omega}\frac{a(x, t)}{k(x, t)}\, dx $.

    Proof. First, we consider the case

    $ T0Ωkt(x,t)k2(x,t)dxdt0. $

    Let $ \varphi(x, t) $ be the positive eigenfunction corresponding to the principal eigenvalue $ \lambda(\mu; a) $. Taking $ \alpha > 0 $ satisfies

    $ lnα=T0Ωkt(x,t)lnφ(x,t)k2(x,t)dxdtT0Ωkt(x,t)k2(x,t)dxdt. $

    Then $ \varphi_\alpha: = \alpha\varphi(x, t) $ is also the principal eigenfunction of Eq (2.2). It is easy to obtain

    $ λ(μ;a)T0Ω1k(x,t)dxdt=T0Ωa(x,t)k(x,t)dxdtμT0ΩΔφαφαdxdt=T0Ωa(x,t)k(x,t)dxdtμT0Ω|Dφα|2φ2αdxdt. $ (2.4)

    This implies that Eq (3.2) holds.

    Next, we consider the case of

    $ T0Ωkt(x,t)k2(x,t)dxdt=0. $

    We can find smooth $ T $-periodic functions $ \{k_n(x, t)\} $ such that

    $ limnkn(x,t)=k(x,t) in C(ˉΩ×[0,T]), $

    and

    $ T0Ω(kn(x,t))tk2n(x,t)dxdt0. $

    It follows from Lemma 2.3 that

    $ limnλn(μ;a)=λ(μ;a), $

    where $ \lambda_n(\mu; a) $ is the principal eigenvalue of Eq (2.2) with $ k(x, t) $ replaced by $ k_n(x, t) $. It is clear from Eq (2.4) that

    $ λn(μ;a)=T0Ωa(x,t)kn(x,t)dxdtμT0Ω|Dφα|2φ2αdxdt. $

    Letting $ n\to\infty $, we have Eq (3.2).

    Remark 2.2. In Eq (3.2), we obtain upper estimates for the principal eigenvalue of the Neumann problem Eq (2.2). Indeed, let $ \lambda_D $ be the principal eigenvalue of the eigenvalue problem

    $ {utμk(x,t)Δua(x,t)u=λDu, in Ω×(0,T],u=0, on Ω×(0,T],u(x,0)=u(x,T), in Ω. $

    By a similar way as in [2], we can show

    $ λD1TT0[μk(x,s)+a(s)]ds, $

    for any $ \mu > 0 $.

    In this section, we study the existence and uniqueness of positive solutions of Eq (1.1). First, we show that if Eq (1.2) has negative principal eigenvalues, then Eq (1.1) has a unique positive solution. To this end, we recall the upper-lower solutions of Eq (1.1). For the sake of convenience, let

    $ QT=Ω×(0,T],Q1=Ω×(0,T]. $

    Definition 3.1. The continuous function $ \bar u(x, t) $ is called the upper-solution of Eq (1.1), if

    $ {utμk(x,t)Δu+m(x,t)uc(x,t)up,inQT,uν0,onQ1,u(x,0)u(x,T),inΩ, $

    is satisfied.

    The definition of the lower-solution is similar. We then can prove the following result, see [2,4,5,15].

    Theorem 3.1. Suppose that $ \bar u(x, t) $, $ \underline u(x, t) $ are a pair of ordered bounded upper-lower solutions of Eq (1.1). Then Eq (1.1) has a unique positive periodic solution $ \theta_\mu(x, t)\in C^{1+\alpha, (1+\alpha)/2}(\bar Q_T) $ that satisfies

    $ u_(x,t)θμ(x,t)ˉu(x,t) in ˉQT. $

    Proof. Let

    $ f(x,t,u)=m(x,t)uc(x,t)up. $

    Then there exists a constant $ K > 0 $ such that

    $ |f(x,t,u2)f(x,t,u1)|K|u2u1|, $

    for any $ (x, t, u_i)\in\bar Q_T\times[\underline u(x, t), \bar u(x, t)] $, $ i = 1, 2 $. It follows from $ L^p $ theory that for any $ u\in C^{\alpha, \frac{\alpha}{2}}(\bar Q_T) $ satisfying $ [\underline u(x, t), \bar u(x, t)] $, the linear initial value problem

    $ {vtμk(x,t)Δv+Kv=Ku+f(x,t,u), in QT,vν=0, on Q1,v(x,0)=u(x,T), in Ω, $

    admits a unique solution $ v $. Thus, a nonlinear operator $ v = \mathcal{F}u $ is defined. We will prove that there is a fixed point for $ \mathcal{F} $ in four steps.

    $ \bf{Step\; 1.} $ In this step, we prove that if $ \underline u\leq u_1\leq u_2\leq\bar u $, then $ \underline u\leq v_1 = \mathcal{F} u_1\leq v_2 = \mathcal{F} u_2\leq\bar u $.

    Take $ \omega_1 = v_2-v_1 $, then $ \omega_1 $ satisfies

    $ {[ω1]tμk(x,t)Δω1+Kω1=K(u2u1)+f(x,t,u2)f(x,t,u1)0, in QT,ω1ν=0, on Q1,ω1(x,0)=u2(x,T)u1(x,T)0, in Ω. $

    By the comparison principle, we obtain $ \omega_1\geq0 $. This implies $ \mathcal{F}u_2\geq\mathcal{F}u_1 $. Similarly, let $ \omega_2 = v_1-\underline u $, then $ \omega_2 $ satisfies

    $ {[ω2]tμk(x,t)Δω2+Kω2=K(u1u_)+f(x,t,u1)f(x,t,u_)0, in QT,ω2ν=0, on Q1,ω1(x,0)=u1(x,T)u_(x,T)0, in Ω. $

    Thus, $ \underline u\leq v_1 $. Similarly, $ v_2\leq\bar u $.

    $ \bf{Step\; 2.} $ In this step, we construct a convergent monotone sequence.

    The iterative sequences $ \{u_n\} $ and $ \{v_n\} $ are constructed as follows:

    $ u1=Fˉu,u2=Fu1,,un=Fun1, $
    $ v1=Fu_,v2=Fv1,,vn=Fvn1. $

    Since $ \underline u\leq\bar u $ and $ \mathcal{F} $ is monotonically non-decreasing, then

    $ u_v1u1ˉu. $

    Similarly, we obtain

    $ u_vnunˉu. $

    And since $ \underline u\leq u_1\leq\bar u $,

    $ u_u2u1ˉu. $

    By induction, we have $ u_{n+1}\leq u_n $. In the same way, we obtain $ v_n\leq v_{n+1} $. Thus, we have

    $ u_v1v2u2u1ˉu. $

    This also implies that $ \{u_n\} $ and $ \{v_n\} $ are monotonically bounded sequences, so there are $ u_0(x, t) $ and $ v_0(x, t) $ such that

    $ limnun(x,t)=u0(x,t),limnvn(x,t)=v0(x,t). $

    Thus

    $ u_(x,t)v0(x,t)u0(x,t)ˉu(x,t) in ˉQT. $

    $ \bf{Step\; 3.} $ In this step, we prove that $ u_0(x, t) $ and $ v_0(x, t) $ are solutions of Eq (1.1).

    Take $ E = W^{2, 1}_p(Q_T)(p > n+2) $. First, we prove that $ \mathcal{F}:D\to C(\bar Q_T) $ is a compact operator, where

    $ D={u(x,t)E:u_(x,t)u(x,t)ˉu(x,t) in ˉQT}. $

    For $ u_1 $, $ u_2\in E $, let $ v_1 = \mathcal{F}u_1 $ and $ v_2 = \mathcal{F}u_2 $, then $ \omega_3 = v_2-v_1 $ satisfies

    $ {[ω3]tμk(x,t)Δω3+Kω3=K(u2u1)+f(x,t,u2)f(x,t,u1), in QT,ω3ν=0, on Q1,ω1(x,0)=u2(x,T)u1(x,T), in Ω. $

    By the $ L^p $ estimate and embedding theorem, it follows that

    $ ω3C1+α,1+α2(ˉQT)Cω3W2,1p(QT)C1(u2u1Lp(QT)+u2(x,T)u1(x,T)Lp(Ω)), $

    here $ C $ and $ C_1 $ are positive constants. Thus $ \mathcal{F}:D\to C(\bar Q_T) $ is continuous. It is known from the embedding theorem that if $ u $ is bounded in $ W^{2, 1}_p(Q_T) $, then $ \mathcal{F}u $ is bounded in $ C^{1+\alpha, (1+\alpha)/2}(\bar Q_T) $. This means that $ \mathcal{F}:D\to C(\bar Q_T) $ is a compact operator.

    Since $ u_n $ is bounded, $ u_n = \mathcal{F}u_{n-1} $ has a convergent subsequence in $ C(\bar Q_T) $. By the monotonicity of $ u_n $,

    $ limnun(x,t)=u0(x,t) in C(ˉQT). $

    Thus $ u_0(x, t) $ is the solution of Eq (1.1) in $ W^{2, 1}_p(Q_T) $. The embedding theorem is used again to get $ u_0(x, t)\in C^{1+\alpha, (1+\alpha)/2}(\bar Q_T) $. In the same way, we get that $ v_0(x, t) $ is the classical solution of Eq (1.1).

    $ \bf{Step\; 4.} $ In this step, we prove the uniqueness and periodicity of the solution of Eq (1.1).

    Since $ k(x, t) $, $ m(x, t) $ and $ c(x, t) $ are periodic about $ t $, then $ \tau(x, t) = u_0(x, t+T)-u_0(x, t) $ satisfies

    $ {τt(x,t)μk(x,t)Δτ(x,t)=m(x,t)τ(x,t)pc(x,t)˜up1(x,t)τ(x,t), in QT,τν=0, on Q1,τ(x,0)=0, in Ω, $ (3.1)

    here $ \tilde{u}(x, t) $ is between $ u_0(x, t+T) $ and $ u_0(x, t) $. It is well known that the solution of Eq (3.1) is unique, thus $ u_0(x, t+T)\equiv u_0(x, t) $ in $ \bar Q_T $.

    The uniqueness of the solution is based on the results of [2,4] and can also be found in recent research results [11,15]. Assume that $ v_1 $ and $ v_2 $ are two positive periodic solutions of Eq (1.1). We first prove that there exists a large constant $ M > 1 $ such that

    $ M1v1<v2<Mv1 in QT. $

    Indeed, it is clear that there exists $ M_1 > 1 $ such that

    $ v2(x,0)=v2(x,T)<M1v1(x,T)=M1v1(x,0) in Ω. $

    This implies $ v_2(x, 0)\not\equiv M_1v_1(x, 0) $ on $ \bar\Omega $. Let $ \eta(x, t): = M_1v_1(x, t)-v_2(x, t) $, then $ \eta(x, t) $ satisfies

    $ ηtμk(x,t)Δη=m(x,t)ηc(x,t)[M1vp1vp2]>m(x,t)ηc(x,t)[(M1v1)pvp2]=m(x,t)ηpc(x,t)ςp1(x,t)η, $

    where $ \varsigma(x, t) $ is lying between $ M_1v_1(x, t) $ and $ v_2(x, t) $. Notice that $ \frac{\partial\eta}{\partial\nu} = 0 $ on $ Q_1 $. By the maximum principle, we have $ \eta > 0 $ in $ \bar Q_T $. Similarly, we can obtain that there exists $ M_2 > 0 $ such that $ v_1 < M_2v_2 $ in $ \bar Q_T $. Take $ M = \max\{M_1, M_2\} $, then we have

    $ M1v1<v2<Mv1 in QT. $

    We know that $ Mv_1(x, t) $ and $ M^{-1}v_1(x, t) $ are a pair of upper-lower solutions of Eq (1.1). According to the second step, Eq (1.1) has a minimum solution $ u_* $ and a maximum solution $ u^* $, which satisfies $ u_*\leq v\leq u^* $ in $ \bar Q_T $ for all solution $ v $ satisfying $ M^{-1}v_1\leq v\leq Mv_1 $. Thus, we obtain $ u_*\leq v_i\leq u^* $ for $ i = 1, 2 $. Hence, $ u_*\equiv u^* $ implies the uniqueness of the solution to Eq (1.1). Set

    $ α=inf{α>0|u(x,t)αu(x,t) in ˉQT}. $

    It is clear that $ \alpha_{*}\geq1 $. Note that if $ \alpha_{*} = 1 $, then $ u_*(x, t)\equiv u^*(x, t) $ in $ \bar Q_T $. Assume that $ \alpha_{*} > 1 $. Let $ \pi(x, t) = \alpha_{*} u_*-u^* $. It is known from the maximum principle that $ \pi(x, t) > 0 $ in $ \bar Q_T $. On the other hand, we know that

    $ π(x,0)=π(x,T)α1u(x,T)=α1u(x,0) on ˉΩ, $

    for some small $ \alpha_1 > 0 $. We can use the previous method to prove the existence of $ M $ to show that

    $ π(x,t)α1u(x,t) on ˉQT. $

    Then we have $ u^*(x, t)\leq (\alpha_*-\alpha) u_*(x, t) $. This is in contradiction with the definition of $ \alpha_* $. Thus, we obtain $ \alpha_{*} = 1 $. The uniqueness is proved.

    Lemma 3.1. If $ \lambda(\mu; m) < 0 $, then Eq (1.1) admits a unique positive periodic solution $ \theta_\mu(x, t)\in C^{1+\alpha, (1+\alpha)/2}(\bar Q_T) $. Moreover, $ \theta_\mu(x, t) $ is globally asymptotically stable.

    Proof. Let $ \theta(x, t) $ be a principal eigenfunction of Eq (1.2) normalized by $ \|\theta(x, t)\|_{C(\bar Q_T)} = 1 $. Then $ \underline{u} = \varepsilon\theta(x, t) $ is a lower-solution of Eq (1.1) for any

    $ 0<ε[λ(μ;m)maxˉQTc(x,t)]1p. $

    Take

    $ M>max{1,[λ(μ;m)minˉQTc(x,t)]1p}. $

    Then we have $ \bar u = M\theta(x, t) $ is an upper-solution of Eq (1.1). From Theorem 3.1, we get that Eq (1.1) has a unique positive solution $ \theta_\mu(x, t) $.

    Since $ \theta_\mu(x, t) $ is the solution of Eq (1.1), then $ \lambda(\mu; m-c\theta_\mu^{p-1}) = 0 $. Let $ \lambda_1 $ be the eigenvalue of the linear problem

    $ {utμk(x,t)Δu[m(x,t)pc(x,t)θp1μ]u=λ1u, in Ω×(0,T],uν=0, on Ω×(0,T],u(x,0)=u(x,T), in Ω. $

    Due to

    $ m(x,t)c(x,t)ˆup1>m(x,t)pc(x,t)ˆup1, $

    for $ u > 0 $. Then $ \lambda_1 > 0 $. It follows from Theorem 2.1 that $ r < 1 $. Thus, $ \theta_\mu(x, t) $ is locally asymptotically stable. In addition, we can choose $ \varepsilon $ small enough and $ M $ large enough such that $ \varepsilon\theta(x, t) $ and $ M\theta(x, t) $ are a pair of ordered bounded upper-lower solutions of Eq (1.1). Then we know that $ \theta_\mu(x, t) $ is globally asymptotically stable by the standard iteration argument as in [2].

    Lemma 3.2. If (1.1) has a positive periodic solution, then $ \lambda(\mu; m) < 0 $.

    Proof. Let $ \theta_\mu(x, t) $ be a positive periodic solution of Eq (1.1). Thanks to [2], we can get that Eq (1.1) is equivalent to

    $ (IKμ)θμ(x,0)=T0Uμ(T,τ)c(x,τ)θpμ(x,τ)dτ in X1. $

    Notice that $ \theta_\mu(x, t) > 0 $. We now apply [2, Theorem 7.3] to obtain

    $ eλ(μ;m)T>1. $

    Thus, $ \lambda(\mu; m) < 0 $.

    Proposition 3.1. If $ \int_0^Tm_*(t)\, dt > 0 $, then Eq (1.1) admits a unique positive periodic solution $ \theta_\mu(x, t) $ for all $ \mu > 0 $.

    Proof. According to Lemma 2.4, we know that

    $ λ(μ;m)T0m(t)dtT0k(t)dt. $ (3.2)

    Due to $ \int_0^Tm_*(t)\, dt > 0 $, $ \lambda(\mu; m) < 0 $. This together with Lemma 3.1 implies that Eq (1.1) admits a unique positive periodic solution for all $ \mu > 0 $.

    In this section, we study the asymptotic behavior of the positive periodic solution of Eq (1.1) when the diffusion rate is large. Here we use regularity estimates together with the perturbation technique to prove our main result. To do this, we first consider the perturbation equation

    $ {ut=μk(x,t)Δu+m(x,t)(u+ε)c(x,t)up, in QT,uν=0, on Q1,u(x,0)=u(x,T), in Ω, $ (4.1)

    where the parameter $ \varepsilon > 0 $.

    Lemma 4.1. Assume that Eq (1.3) holds. Then Eq (4.1) has a positive periodic solution $ \theta^{\varepsilon}_{\mu}(x, t) $ for $ \mu > 0 $, provided $ \varepsilon > 0 $ is small. Moreover, we can find $ \mu_1 > 0 $ such that

    $ limε0+θεμ(x,t)=θμ(x,t) in C1,12(ˉQT), $ (4.2)

    for $ \mu\geq\mu_1 $.

    Proof. Through a similar argument as in Theorem 3.1, we know the existence of the positive solution $ \theta^{\varepsilon}_{\mu}(x, t) $ to Eq (4.1). We only need to prove Eq (4.2). Let $ \sigma(x, t) = \theta^{\varepsilon}_{\mu}(x, \frac{t}{\mu}) $. Then $ \sigma(x, t) $ satisfies

    $ {σt=k(x,tμ)Δσ+1μ[m(x,tμ)(σ+ε)c(x,tμ)σp], in QT,σν=0, on Q1,σ(x,0)=θεμ(x,0), in Ω. $

    It is known from the $ L^p $ estimate that there exists $ \mu_1 > 0 $ such that $ \sigma(x, t) $ is bounded in $ W^{2, 1}_p(\Omega\times(0, \mu T]) $ for any $ \mu > \mu_1 $. It then follows that $ \theta^{\varepsilon}_{\mu}(x, t) $ is bounded in $ W^{2, 1}_p(Q_T) $ for any $ \mu > \mu_1 $. Then, taking $ p $ large enough, we know from the embedding theorem that $ \theta^{\varepsilon}_{\mu}(x, t) $ is compact in $ C^{1, \frac{1}{2}}(\bar Q_T) $. Thus there is a subsequence, still denoted by $ \theta_\mu^{\varepsilon}(x, t) $, such that

    $ limε0+θεμ(x,t)=ν(x,t) in C1,12(ˉQT), $ (4.3)

    for some nonnegative periodic function $ \nu(x, t)\in C(\bar Q_T) $. It follows from the argument of Lemma 3.1 that $ \varepsilon\theta(x, t) $ is a lower-solution of Eq (4.1). Thus we have $ \nu(x, t) > 0 $ for $ (x, t)\in\bar\Omega\times[0, T] $. Since $ \theta^{\varepsilon}_{\mu}(x, t) $ is bounded and Eq (4.3), $ \nu $ satisfies

    $ ν(x,t)=ν(x,0)+μt0[k(x,s)Δν(x,s)+m(x,s)νc(x,s)νp]ds. $

    It is easy to obtain

    $ {νt=μk(x,t)Δν+m(x,t)νc(x,t)νp, in QT,νν=0, on Q1,ν(x,0)=ν(x,T), in Ω. $

    By standard parabolic regularity, we know that $ \nu(x, t)\in C^{1+\alpha, (1+\alpha)/2}(\bar Q_T) $. The uniqueness of the solution means that Eq (4.2) holds.

    At the end of this section, we prove Theorem 1.2.

    Proof of Theorem 1.2. We divide the proof into the following three steps.

    $ \bf{Step\; 1.} $ In this step, we prove that $ \theta_\mu(x, t) $ has a convergent subsequence as $ \mu\to\infty $.

    It follows from a similar argument to Lemma 4.1 that there exists $ \mu_1 > 0 $ such that $ \theta_\mu(x, t) $ is compact in $ C^{1, \frac{1}{2}}(\bar Q_T) $ for any $ \mu > \mu_1 $. Thus, by passing to a subsequence, there is a nonnegative periodic function $ \theta\in C(\bar Q_T) $ such that

    $ limμθμ(x,t)=θ(x,t) in C1,12(ˉQT). $

    $ \bf{Step\; 2.} $ In this step, we show that $ \theta(x, t) $ is independent of $ t $.

    Let $ f(t) $ be a smooth $ T $-periodic function, then we have

    $ T0θμ(x,t)ft(t)dt=μT0k(x,t)f(t)Δθμ(x,t)dt+T0m(x,t)θμ(x,t)f(t)dtT0c(x,t)θpμ(x,t)f(t)dt. $

    By dividing $ \mu $ and making $ \mu\to\infty $, we have

    $ T0k(x,t)f(t)Δθ(x,t)dt=0. $

    Since $ f(t) $ is arbitrary, we obtain

    $ Δθ(x,t)=0. $

    Then we derive

    $ Ω|θ(x,t)|2dx=0. $

    Thus we have $ \theta(x, t)\equiv\theta(t) $ for $ x\in\bar\Omega $.

    $ \bf{Step\; 3.} $ In this step, we show that $ \theta(t) = \omega(t) $ in $ [0, T] $.

    First, we assert that $ \theta(t)\in C^{1}((0, \infty)) $. Indeed, it is easy to obtain from Eq (1.1) that

    $ t+εtΩus(x,s)k(x,s)dxds=t+εtΩm(x,s)k(x,s)u(x,s)dxdst+εtΩc(x,s)k(x,s)up(x,s)dxds. $

    Then we have

    $ Ωu(x,t+ε)k(x,t+ε)dxΩu(x,t)k(x,t)dx+t+εtΩkt(x,s)k2(x,s)u(x,s)dxds=t+εtΩm(x,s)k(x,s)u(x,s)dxdst+εtΩc(x,s)k(x,s)up(x,s)dxds. $

    Taking $ \mu\to\infty $, we obtain

    $ Ωθ(t+ε)k(x,t+ε)dxΩθ(t)k(x,t)dx+t+εtΩkt(x,s)k2(x,s)θ(s)dxds=t+εtΩm(x,s)k(x,s)θ(s)dxdst+εtΩc(x,s)k(x,s)θp(s)dxds. $

    Thus, we derive

    $ [θ(t)Ω1k(x,t)dx]t=Ωm(x,t)k(x,t)dxθ(t)Ωc(x,t)k(x,t)dxθp(t)Ωkt(x,t)k2(x,t)dxθ(t). $

    Hence,

    $ θt(t)=1k(t)Ωm(x,t)k(x,t)dxθ(t)1k(t)Ωc(x,t)k(x,t)dxθp(t)1k(t)Ωkt(x,t)k2(x,t)dxθ(t)[k(t)]tk(t)θ(t), $

    for $ t > 0 $. Thus $ \theta(t)\in C^{1}((0, \infty)) $ holds.

    We then prove that $ \theta(t) $ satisfies Eq (1.5). It is obvious from Eq (1.1) that

    $ Ωut(x,t)k(x,t)dx=Ωm(x,t)k(x,t)u(x,t)dxΩc(x,t)k(x,t)up(x,t)dx. $ (4.4)

    Similarly, suppose that $ f(t) $ is a smooth $ T $-periodic function. Multiplying $ f(t) $ on both sides of Eq (4.4) and integrating over $ [0, T] $ gives

    $ T0Ωu(x,t)[f(t)k(x,t)]tdxdt=T0Ωm(x,t)k(x,t)u(x,t)f(t)dxdtT0Ωc(x,t)k(x,t)up(x,t)f(t)dxdt. $

    Letting $ \mu\to\infty $, we know

    $ T0Ωθ(t)[f(t)k(x,t)]tdxdt=T0Ωm(x,t)k(x,t)θ(t)f(t)dxdtT0Ωc(x,t)k(x,t)θp(t)f(t)dxdt. $

    This implies

    $ T0Ωf(t)k(x,t)θt(t)dxdt=T0Ωm(x,t)k(x,t)θ(t)f(t)dxdtT0Ωc(x,t)k(x,t)θp(t)f(t)dxdt. $

    By the arbitrary of $ f(t) $, it follows that

    $ Ω1k(x,t)dxθt(t)=Ωm(x,t)k(x,t)dxθ(t)Ωc(x,t)k(x,t)dxθp(t). $

    Thus, we have

    $ {θt=˜M(t)˜k(t)θ˜C(t)˜k(t)θp,tR,θ(0)=θ(T). $

    Finally, we prove that $ \theta(t) > 0 $ in $ t\in[0, T] $. We define $ \theta^{\varepsilon}_{\mu}(x, t) $ as the unique positive periodic solution of Eq (4.1) for small $ \varepsilon > 0 $ and large $ \mu $. Similarly to the previous argument, it can be shown that

    $ limμθεμ(x,t)=ωε(t) in C1,12(ˉQT), $ (4.5)

    where $ \omega^{\varepsilon}(t) $ satisfies

    $ {ωεt=M(t)k(t)(ωε+ε)c(t)k(t)(ωε)p,tR,ωε(0)=ωε(T). $ (4.6)

    Since

    $ Ωm(x,t)k(x,t)dx>0 in [0,T], $

    we can obtain that Eq (1.5) admits a unique periodic positive solution $ \omega(t) $. Note that $ \omega(t) $ is a lower solution of Eq (4.6). Then there exists a unique positive periodic solution $ \omega^{\varepsilon}(t) $ to Eq (4.6). In addition, $ \omega^{\varepsilon}(t) $ is monotonically increasing about $ \varepsilon $, and $ \omega_{\varepsilon}(t)\geq\omega(t) > 0 $. We obtain that there exists a positive continuous function $ \omega_{0}(t) $ such that

    $ limε0+ωε(t)=ω0(t) in [0,T]. $

    The uniqueness of the positive solution of Eq (4.6) implies that

    $ limε0+ωε(t)=ω(t) in [0,T]. $

    It follows from Lemma 4.1 that

    $ limε0+θεμ(x,t)=θμ(x,t) in C1,12(ˉQT). $

    This means that $ \theta(t) $ is positive, together with (4.4)–(4.6). Thus, we must have

    $ θ(t)ω(t) in [0,T]. $

    This also implies that

    $ limμθμ(x,t)=ω(t) in C1,12(ˉQT), $

    holds for the entire sequence.

    We consider the positive solutions of the periodic-parabolic logistic equation with indefinite weight function and nonhomogeneous diffusion coefficient. When the dispersal rate is small, we can obtain a similar result as in the homogeneous diffusion coefficient. Here we are interested in the case of large diffusion coefficient with nonhomogeneous diffusion coefficient.

    In Theorem 1.1, we obtain the condition of $ m(x, t) $ to guarantee a positive periodic solution for all $ \mu > 0 $. Then we investigate the effect of large $ \mu $ on the positive solution and establish that the limiting profile is determined by the positive solution of Eq (1.5). More precisely, we prove that the positive periodic solution tends to the unique positive solution of the corresponding non-autonomous logistic equation when the diffusion rate is large.

    M. Fan was responsible for writing the original draft. J. Sun handled the review and supervision.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The authors would like to thank the anonymous reviewers for the careful reading and several valuable comments to revise the paper. Fan was supported by Gansu postgraduate scientifc research (20230XZX-055) and Sun was supported by NSF of China (12371170) and NSF of Gansu Province of China (21JR7RA535, 21JR7RA537).

    The authors declare there is no conflict of interest.



    [1]

    Y. Achdou, Finite difference methods for mean field games, in Hamilton-Jacobi Equations: Approximations, Numerical Analysis And Applications, Lecture Notes in Mathematics, 2074, Springer, Heidelberg, 2013, 1–47.

    [2] Hamilton–Jacobi equations constrained on networks. NoDEA Nonlinear Differential Equations Appl. (2013) 20: 413-445.
    [3] Mean field games: numerical methods. SIAM J. Numer. Anal. (2010) 48: 1136-1162.
    [4] Hamilton-Jacobi equations for optimal control on junctions and networks. ESAIM Control Optim. Calc. Var. (2015) 21: 876-899.
    [5] Convergence of a finite difference scheme to weak solutions of the system of partial differential equation arising in mean field games. SIAM J. Numer. Anal. (2016) 54: 161-186.
    [6]

    L. Boccardo, F. Murat and J.-P. Puel, Existence de solutions faibles pour des équations elliptiques quasi-linéaires à croissance quadratique, in Nonlinear Partial Differential Equations and Their Applications, Res. Notes in Math, 84, Pitman, Boston, Mass.-London, 1983, 19–73.

    [7] Stationary mean field games systems defined on networks. SIAM J. Control Optim. (2016) 54: 1085-1103.
    [8] The vanishing viscosity limit for Hamilton-Jacobi equations on networks. J. Differential Equations (2013) 254: 4122-4143.
    [9]

    P. Cardaliaguet, Notes on mean field games, Preprint, 2011.

    [10]

    R. Carmona and F. Delarue, Probabilistic Theory of Mean Field Games with Applications. I-II, Springer, 2017.

    [11]

    M.-K. Dao, Ph.D. thesis, 2018.

    [12] Vertex control of flows in networks. Netw. Heterog. Media (2008) 3: 709-722.
    [13]

    A. Ern and J.-L. Guermond, Theory and Practice of Finite Elements, Applied Mathematical Sciences, 159, Springer-Verlag, New York, 2004.

    [14] Homogenization of second order discrete model with local perturbation and application to traffic flow. Discrete Contin. Dyn. Syst. (2017) 37: 1437-1487.
    [15] Diffusion processes on graphs: Stochastic differential equations, large deviation principle. Probab. Theory Related Fields (2000) 116: 181-220.
    [16] Diffusion processes on graphs and the averaging principle. Ann. Probab. (1993) 21: 2215-2245.
    [17]

    M. Garavello and B. Piccoli, Traffic Flow on Networks, AIMS Series on Applied Mathematics, 1, American Institute of Mathematical Sciences, Springfield, MO, 2006.

    [18]

    D. Gilbarg and N. S. Trudinger, Elliptic Partial Differential Equations of Second Order, Reprint of the 1998 edition, Classics in Mathematics, Springer-Verlag, Berlin, 2001.

    [19] On the existence of classical solutions for stationary extended mean field games. Nonlinear Anal. (2014) 99: 49-79.
    [20] Local regularity for mean-field games in the whole space. Minimax Theory Appl. (2016) 1: 65-82.
    [21] Time-dependent mean-field games in the subquadratic case. Comm. Partial Differential Equations (2015) 40: 40-76.
    [22]

    O. Guéant, J.-M. Lasry and P.-L. Lions, Mean field games and applications, in Paris-Princeton Lectures on Mathematical Finance 2010, Lecture Notes in Mathematics, 2003, Springer, Berlin, 2011,205–266.

    [23] An invariance principle in large population stochastic dynamic games. J. Syst. Sci. Complex. (2007) 20: 162-172.
    [24] Large-population cost-coupled LQG problems with nonuniform agents: Individual-mass behavior and decentralized ϵ-Nash equilibria. IEEE Trans. Automat. Control (2007) 52: 1560-1571.
    [25] Large population stochastic dynamic games: Closed-loop McKean-Vlasov systems and the Nash certainty equivalence principle. Commun. Inf. Syst. (2006) 6: 221-251.
    [26] A Hamilton-Jacobi approach to junction problems and application to traffic flows. ESAIM Control Optim. Calc. Var. (2013) 19: 129-166.
    [27]

    C. Imbert and R. Monneau, Flux-limited solutions for quasi-convex Hamilton-Jacobi equations on networks, Ann. Sci. Éc. Norm. Supér., 50 (2017), 357–448.

    [28] Jeux à champ moyen. Ⅰ. Le cas stationnaire. C. R. Math. Acad. Sci. Paris (2006) 343: 619-625.
    [29] Jeux à champ moyen. Ⅱ. Horizon fini et contrôle optimal. C. R. Math. Acad. Sci. Paris (2006) 343: 679-684.
    [30] Mean field games. Jpn. J. Math. (2007) 2: 229-260.
    [31] Viscosity solutions for junctions: Well posedness and stability. Atti Accad. Naz. Lincei Rend. Lincei Mat. Appl. (2016) 27: 535-545.
    [32] Well-posedness for multi-dimensional junction problems with Kirchoff-type conditions. Atti Accad. Naz. Lincei Rend. Lincei Mat. Appl. (2017) 28: 807-816.
    [33] Weak solutions to Fokker-Planck equations and mean field games. Arch. Rational Mech. Anal. (2015) 216: 1-62.
    [34] On the weak theory for mean field games systems. Boll. U.M.I. (2017) 10: 411-439.
    [35] Classical solvability of linear parabolic equations on networks. J. Differential Equations (1988) 72: 316-337.
  • This article has been cited by:

    1. Peng E, Tingting Xu, Linhua Deng, Yulin Shan, Miao Wan, Weihong Zhou, Solutions of a class of higher order variable coefficient homogeneous differential equations, 2025, 20, 1556-1801, 213, 10.3934/nhm.2025011
  • Reader Comments
  • © 2019 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(3674) PDF downloads(296) Cited by(6)

Figures and Tables

Figures(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog