Evacuation dynamics influenced by spreading hazardous material

  • In this article, an evacuation model describing the egress in case of danger is considered. The underlying evacuation model is based on continuous network flows, while the spread of some gaseous hazardous material relies on an advection-diffusion equation. The contribution of this work is twofold. First, we introduce a continuous model coupled to the propagation of hazardous material where special cost functions allow for incorporating the predicted spread into an optimal planning of the egress. Optimality can thereby be understood with respect to two different measures: fastest egress and safest evacuation. Since this modeling approach leads to a pde/ode-restricted optimization problem, the continuous model is transferred into a discrete network flow model under some linearity assumptions. Second, it is demonstrated that this reformulation results in an efficient algorithm always leading to the global optimum. A computational case study shows benefits and drawbacks of the models for different evacuation scenarios.

    Citation: Simone Göttlich, Sebastian Kühn, Jan Peter Ohst, Stefan Ruzika, Markus Thiemann. Evacuation dynamics influenced by spreading hazardous material[J]. Networks and Heterogeneous Media, 2011, 6(3): 443-464. doi: 10.3934/nhm.2011.6.443

    Related Papers:

    [1] Tao Xie, Xing Xiong . Finite-time synchronization of fractional-order heterogeneous dynamical networks with impulsive interference via aperiodical intermittent control. AIMS Mathematics, 2025, 10(3): 6291-6317. doi: 10.3934/math.2025287
    [2] Chengqiang Wang, Xiangqing Zhao, Yang Wang . Finite-time stochastic synchronization of fuzzy bi-directional associative memory neural networks with Markovian switching and mixed time delays via intermittent quantized control. AIMS Mathematics, 2023, 8(2): 4098-4125. doi: 10.3934/math.2023204
    [3] Qinghua Zhou, Li Wan, Hongbo Fu, Qunjiao Zhang . Exponential stability of stochastic Hopfield neural network with mixed multiple delays. AIMS Mathematics, 2021, 6(4): 4142-4155. doi: 10.3934/math.2021245
    [4] Zhifeng Lu, Fei Wang, Yujuan Tian, Yaping Li . Lag synchronization of complex-valued interval neural networks via distributed delayed impulsive control. AIMS Mathematics, 2023, 8(3): 5502-5521. doi: 10.3934/math.2023277
    [5] Dong Pan, Huizhen Qu . Finite-time boundary synchronization of space-time discretized stochastic fuzzy genetic regulatory networks with time delays. AIMS Mathematics, 2025, 10(2): 2163-2190. doi: 10.3934/math.2025101
    [6] Rakkiet Srisuntorn, Wajaree Weera, Thongchai Botmart . Modified function projective synchronization of master-slave neural networks with mixed interval time-varying delays via intermittent feedback control. AIMS Mathematics, 2022, 7(10): 18632-18661. doi: 10.3934/math.20221025
    [7] Zhengqi Zhang, Huaiqin Wu . Cluster synchronization in finite/fixed time for semi-Markovian switching T-S fuzzy complex dynamical networks with discontinuous dynamic nodes. AIMS Mathematics, 2022, 7(7): 11942-11971. doi: 10.3934/math.2022666
    [8] Tian Xu, Jin-E Zhang . Intermittent control for stabilization of uncertain nonlinear systems via event-triggered mechanism. AIMS Mathematics, 2024, 9(10): 28487-28507. doi: 10.3934/math.20241382
    [9] Xingxing Song, Pengfei Zhi, Wanlu Zhu, Hui Wang, Haiyang Qiu . Exponential synchronization control of delayed memristive neural network based on canonical Bessel-Legendre inequality. AIMS Mathematics, 2022, 7(3): 4711-4734. doi: 10.3934/math.2022262
    [10] Saravanan Shanmugam, Mohamed Rhaima, Hamza Ghoudi . Exponential synchronization analysis for complex dynamical networks with hybrid delays and uncertainties under given control parameters. AIMS Mathematics, 2023, 8(12): 28976-29007. doi: 10.3934/math.20231484
  • In this article, an evacuation model describing the egress in case of danger is considered. The underlying evacuation model is based on continuous network flows, while the spread of some gaseous hazardous material relies on an advection-diffusion equation. The contribution of this work is twofold. First, we introduce a continuous model coupled to the propagation of hazardous material where special cost functions allow for incorporating the predicted spread into an optimal planning of the egress. Optimality can thereby be understood with respect to two different measures: fastest egress and safest evacuation. Since this modeling approach leads to a pde/ode-restricted optimization problem, the continuous model is transferred into a discrete network flow model under some linearity assumptions. Second, it is demonstrated that this reformulation results in an efficient algorithm always leading to the global optimum. A computational case study shows benefits and drawbacks of the models for different evacuation scenarios.


    In recent years, an enormous advance of complex networks (CNS) has been witnessed on account of their extensive applications in many areas such as unmanned aerial vehicles, communication systems, power systems, transportation networks, ecological networks, and so forth [1,2,3,4]. A large number of interesting achievements about dynamical properties of CNS have sprung up [5,6,7,8,9]. It can be observed that CNS in the majority of the above findings are supposed to be single-link. In contrast to single-link complex, multilayer complex networks (MCNS) can more accurately simulate the real systems, which are more general and practical due to possessing sophisticated structures. For instance, a transportation network is composed of highways, ships, railways, and airplanes, and it can be seen as one type of multilayer complex network. Apparently, the transmission speed and cost are different among these ways, which signifies that each interaction mode has different weights. On the other hand, great importance can also be attached to social networks. There are various kinds of contact ways for people to communicate with others including phone, mail, and internet, which leads to multilayer networks. Particularly, large group decision-making plays a significant role in social networks. In [10], large group decision-making with a rough integrated asymmetric cloud model under a multi-granularity linguistic environment was considered, and large group emergency decision-making with bi-directional trust in social networks was studied based on a probabilistic hesitant fuzzy integrated cloud approach in [11], where strong uncertainty and randomness were tackled effectively. Since the MCNS may better depict the connection and interaction among different layer networks, the investigation of MCNS has important theoretical and realistic significance.

    Synchronization is a typical dynamic characteristic of MCNS, which means all the states of different nodes of networks can evolve in one common mode with various different initial data, and plenty of relevant results about MCNS have been reported. Particularly, in [12], global exponential synchronization of multi-link memristive neural networks with time-varying delays was analyzed by designing state feedback controllers. In [13], the practical fixed-time synchronization issue of MCNS was addressed based on an intermittent event-triggered control strategy by utilizing the proposed new practical fixed-time stability lemma. In [14], bipartite synchronization of multilayer signed networks was investigated under aperiodic intermittent-based adaptive dynamic event-triggered control, where the control gains and the triggering parameters were adjusted with the system states. Furthermore, the discussions of various synchronizations were extended to fractional MCNS with the help of impulsive control [15,16]. On the other hand, in the real world, due to the existence of random disturbances and abrupt variations in systems' structures and parameters, stochastic Markovian switching systems were introduced to simulate such systems. Recently, stability and synchronization of stochastic MCNS with Markovian switching have aroused scholars' interests [17,18,19,20,21,22,23]. In [17], $ p $th moment exponential stability for a class of stochastic complex multi-link networks with Markovian switching and multi-delayed impulsive effects was examined by virtue of the Razumikhin approach. In [18], combining with Kirchhoff's matrix tree theorem, finite-time synchronization of stochastic MCNS with Markovian switching was studied via a novel quantized aperiodically intermittent control. Moreover, in [19], the $ p $th moment exponential synchronization issue of stochastic delayed MCNS with semi-Markov jump under aperiodically intermittent control was solved based on the Lyapunov method and graph theory, and several sufficient criteria were derived. Nevertheless, the preceding existing results mainly focus on $ p $th moment synchronization of stochastic MCNS with Markovian switching rather than almost sure exponential synchronization, which sparks the appearance of this work.

    It is worth noting that stochastic noise can be used as the control input [24,25]. For instance, in social networks, by virtue of public opinion control, which can be viewed as noise control, a certain behavior mode of the entire society can keep the consensus in a special direction [26]. Particularly, white noise was incorporated initially to stabilize the neural networks by Liao and Mao [27]. Soon afterward, stochastic stabilization of nonlinear systems was considered by utilizing an intermittent Brownian noise or aperiodical intermittent Brownian noise in [28,29]. Meanwhile, when the discrete-time feedback control was introduced to the drift term of systems, mean-square exponential stabilization of stochastic nonlinear systems was discussed in [30]. While the discrete-time feedback control was introduced to the diffusion term of systems, the almost sure exponential stabilization issue was solved in [31]. Subsequently, the studies about discrete feedback control and noise stabilization were further developed in [32,33,34]. Based on the previous findings, a new periodically intermittent discrete observation control scheme was proposed, by which the issue of mean-square exponential synchronization of stochastic neural networks was coped with in [35]. Furthermore, under the similar control strategy in [35], almost sure exponential stabilization of deterministic neural networks and hybrid neural networks was analyzed in [36,37]. More recently, aperiodically intermittent discrete observation noise control (AIDONC) as one kind of discontinuous control scheme was designed by adjusting the distributions of work/rest intervals. Since feedback control based on intermittent discrete observations on current states are imposed, the control time and control frequency can be reduced. Consequently, energy consumption and control cost can be greatly saved. Besides, efficiency and feasibility of implementation of the control strategy will be improved significantly. By using aperiodically intermittent discrete observation noise control (AIDONC), synchronization of complex networks was investigated in [38,39]. For MCNS, it can be observed that in [21,22], the intermittent control was imposed on the drift term and several criteria on $ p $th moment exponential synchronization were presented. Additionally, in [23], some sufficient conditions on almost sure exponential synchronization of MCNS were given, but the control strategy AIDONC was not adopted.

    Inspired by the above discussions, this paper will investigate the almost sure exponential synchronization of MCNS with Markovian switching via AIDONC. The main contributions are summarized below:

    1) In this paper, Markovian switching and multilayer interaction factors are taken into account simultaneously, which make our system more general compared with the networks in literature [35,36,37,38,39]. Meanwhile, the network architecture may be undirected or directed, and consequently, the adjacency matrix could be symmetrical and asymmetrical.

    2) Different from the work about $ p $th moment exponential synchronization of network systems in [21,22], by utilizing M-matrix theory and various stochastic analysis techniques including the Itô formula, the Gronwall inequality, and the Borel-Cantelli lemma, almost sure exponential synchronization of multilayer coupled networks with Markovian switching is analyzed and the upper bound of the duration time is also estimated.

    3) Compared with the work in [23], one kind of control strategy based on AIDONC is adopted, where the average control rate is integrated to depict the distributions of work/rest intervals of the control strategy from an overall perspective.

    The remainder of this article is organized as follows. In Section 2, some necessary preliminaries and model descriptions are presented. The main findings about almost sure exponential synchronization for multilayer coupled networks with Markovian switching via aperiodically intermittent noise are derived in Section 3. In Section 4, a multilayer Chua's circuits network with Markovian switching is considered and some numerical simulations are carried out to validate the effectiveness of our theoretical results, and conclusions are drawn in the last section.

    For convenience and simplicity, let us introduce some standard notations utilized in the following context. Let $ \mathbb{R} ^n $ be the $ n $-dimensional Euclidean space. If $ x\in \mathbb{R} ^n $, then $ \left | x \right | $ represents the Euclidean norm of $ x $. Let $ I_{N} $ be an $ N\times N $ identity matrix. Besides, the superscript $ T $ is defined as the transpose of a vector or a matrix. For any real symmetric matrix, $ \lambda _{\max}(\cdot) $ and $ \lambda _{\min}(\cdot) $ denote the maximum and minimum eigenvalues of the given matrix, respectively. Let $ (\Omega, \mathcal{F}, \mathbb{P}) $ be a complete probability space with a filtration $ \{\mathcal{F}_{t} \}_{t\ge 0} $ satisfying the usual condition, and $ B(t) $ represents one-dimensional Brownian motion defined on the above probability space.

    In this section, we consider the multilayer coupled networks with Markovian switching and $ M $ nodes as follows:

    $ dδk(t)=[(δk(t),t,r(t))+Hh=1Mm=1ρh(r(t))bhkm(r(t))Υδm(t)]dt+uk(t,r(t)),k=1,2,,M,
    $
    (2.1)

    where $ \delta_{k}(t) = ((\delta_{k1}(t), \delta_{k2}(t), \dots, \delta_{kn}(t))^T\in \mathbb{R} ^n $ denotes the vector of the $ \mathit{k} $th node at time $ t $. $ u_{i}(t, r(t)) $ is the controller to be devised. $ G(\delta_{k}(t), t, r(t)) \in \mathbb{R}^n: \mathbb{R} ^n\times \mathbb{R} ^+\times \mathbb{S}\rightarrow \mathbb{R}^n $ represents the nonlinear activation function with $ G(\delta_{k}(0), r(0), 0) = 0 $. $ \rho _h(r(t)) > 0 $ is the coupling strength of the $ h $th layer. $ \Upsilon $ stands for an inner coupling matrix. The adjacency matrix $ (b_{km}^{h})_{M\times M} $ is the outer coupling configuration matrix, where $ b_{km}^h > 0 $ holds if node $ m $ directs links to $ k $, and otherwise, $ b_{km}^h = 0\; (k\ne m) $. Meanwhile, the diagonal elements $ b_{kk}^h = - { \sum\limits_{m = 1, m\ne k}^{M}}b_{km} ^h $, which leads to $ { \sum\limits_{m = 1}^{M}}b_{km} ^h = 0 $. Let $ r(t), t\ge 0 $ signify a right-continuous Markov chain on the probability space $ (\Omega, \mathcal{F}, \mathbb{P}) $ with initial date $ r(0) = r_{0}\in \mathbb{S} = \{1, 2, \cdots\cdots, N \} $, and generator $ \Psi = (\psi _{ij})_{N\times N} $ given by

    $ P(r(t+Δ)=j|r(t)=i)={ψijΔ+o(Δ),ifij,1+ψijΔ+o(Δ),ifi=j,
    $
    (2.2)

    where $ \Delta > 0 $, $ \psi_{ii} = -{ \sum_{j\ne i}^{}}\psi _{ij} $, and $ \psi _{ij} \ge 0 $ denotes the transition rate from $ i $ to $ j $.

    Subsequently, the isolated node of the network is expressed as

    $ ds(t)=G(s(t),t,r(t))dt,
    $
    (2.3)

    where $ s(t) = (s_{1}(t), s_{2}(t), \dots, s_{n}(t))^T\in\mathbb{R} ^n $. Apparently, the isolated node is independent of other nodes. $ s(t) $ can be considered as the equilibrium state or stable objective trajectory, and in multi-agent dynamics, it stands for the leader. Let $ \eta_{k}(t) = \delta_{k}(t)-s(t) $ be the synchronization error. Accordingly, the error system can be depicted as

    $ dηk(t)=[ˆG(ηk(t),t,r(t))+Hh=1Mm=1ρh(r(t))bhkm(r(t))Υηi(t)]dt+uk(t,r(t)),k=1,2,,M,
    $
    (2.4)

    where $ \hat{G} (\eta _{k}(t), t, r(t)) = G(\eta _{k}(t), t, r(t))- G (s(t), t, r(t)) $. In order to realize synchronization, the aperiodically intermittent feedback control based on discrete-time state observations is designed below:

    $ uk(t)={βk(ηk(μ(t)),t,r(t))dw(t),t[ti,si),0,t[si,ti+1),
    $
    (2.5)

    where $ \mu(t) = \left [ t/\tau \right ]\tau, \tau > 0 $, $ \tau $ denotes the duration between two consecutive observations, and $ \beta_{k}(\eta_{k}(\mu(t), t, r(t)) = C_{k}(r(t))\eta_{k}(\mu(t)) $ represents the diffusion coefficient column vector. Time interval $ [t_{i}, s_{i}) $ is the working time of the $ k $th period while $ [s_{i}, t_{i+1}) $ stands for the rest time. Furthermore, under aperiodically intermittent noise control, error system (2.4) can be rewritten by

    $ {dη(t)=[˜G(η(t),t,r(t))+Hh=1ρh(r(t))(Bh(r(t))Υ)η(t)]dt+˜β(η(μ(t)),t,r(t))dw(t),t[ti,si),i=0,1,dη(t)=˜G(η(t),t,r(t))+Hh=1ρh(r(t))(Bh(r(t))Υ)η(t),t[si,ti+1),i=0,1,
    $
    (2.6)

    where $ \eta (t) = \Big(\eta_{1}^{T} (t), \eta_{2}^{T} (t), \cdots\cdots, \eta_{M}^{T} (t)\Big)^T $, $ \hat{G} (\eta (t)) = \Big(\tilde{G}^{T}(\eta_{1} (t)), \tilde{G}^{T}(\eta_{2} (t)), \cdots\cdots, \tilde{G}^{T}(\eta_{M} (t))\Big)^{T}\in \mathbb{R}^{nM} $, and $ \tilde{\beta}(\eta(\mu(t)), t, r(t)) = \Big(\eta_{1}^{T}(\mu(t))C_{1}^{T}(r(t)), \eta_{2}^{T}(\mu(t))C_{2}^{T}(r(t)), \cdots\cdots, \eta_{M}^{T}(\mu(t))C_{M}^{T}(r(t))\Big)^{T} $. Since multilayer networks, Markovian switching, and aperiodically intermittent controller based on discrete-time observation noise exist simultaneously, this makes them more complex to analyze the synchronization feature. In order to overcome the difficulties, the auxiliary systems with continuous-time observations are presented below,

    $ {dz(t)=[˜G(z(t),t,r(t))+Hh=1ρh(r(t))(Bh(r(t))Υ)z(t)]dt+˜β(z(t),t,r(t))dw(t),t[ti,si),i=0,1,dz(t)=˜G(z(t),t,r(t))+Hh=1ρh(r(t))(Bh(r(t))Υ)z(t),t[si,ti+1),i=0,1,
    $
    (2.7)

    where $ z(t) = \delta(t)-s(t) $.

    In what follows, several necessary definitions and assumptions are provided, which play a fundamental role in acquiring the theoretical results.

    Definition 1. The multilayer coupled networks (2.1) can be almost surely exponentially synchronized if for any initial data $ \eta(t_{0}) = \eta_0 $, there exists a positive constant satisfying the following inequality

    $ limtsup1tlog|η(t)|<κ,a.s.
    $
    (2.8)

    Assumption 1. Suppose that there exist two positive constants $ \alpha\in(0, 1) $ and $ N_{0} $ such that

    $ N(t,s)α(ts)N0,t>st0,
    $
    (2.9)

    where $ N(t, s) $ denotes the total control time length on $ [s, t) $, and $ N_0 $ is the elasticity number.

    Assumption 2. Suppose that there are positive constants $ l_{i} $ and $ \sigma_{ki} $ such that following inequalities hold

    $ |G(u1,t,i)G(u2,t,i)|li|u1u2|,
    $
    (2.10)

    for $ \forall u_{1}, u_{2}\in\mathbb{R}^n, $ and $ i\in \mathbb{S} $.

    In this section, some novel lemmas are first established by utilizing stochastic analysis techniques. Based on the derived lemmas, almost sure exponential synchronization of MCNS with Markovian switching are further examined via AIDONC, and some special cases are also discussed.

    Lemma 3.1. Let $ Q = \rm{diag}\{\chi_{1}, \chi_{2}, \cdots, \chi_{N}\}-\Psi $ denote an M-matrix, where $ q\in(0, 1) $, $ \chi_{i} = 0.5q\Big[(2-q)d_{i}-\sigma_{i}\Big]-q l_{i}-q\Bigg(\sum\limits_{h = 1}^{H}\rho_{h}(i)\lambda_{\mathop{\rm{max}}}\Big(\Big[B^{h}(i)\otimes\Upsilon\Big]_{s}\Big)\Bigg) $, and $ \Psi = (\psi_{ij})_{M\times M} $. Under Assumptions 1 and 2, if $ \epsilon = \tilde{\lambda}\alpha-\tilde{\zeta}(1-\alpha) > 0 $, then the solution $ z(t) $ of Eq (2.7) satisfies that

    $ E|z(t)|qK0E|z0|qeϵt,
    $
    (3.1)

    where $ \sigma_{i} = \max\Big\{\lambda_{\max}\Big(C_{1}^{T}(i)C_{1}(i)\Big), \lambda_{\max}\Big(C_{2}^{T}(i)C_{2}(i)\Big), \cdots\cdots, \lambda_{\max}\Big(C_{M}^{T}(i)C_{M}(i)\Big)\Big\}, $ $ d_{i} = \min\Big\{\lambda_{\min}^2\Big(C_{1}(i)\Big), \lambda_{\min}^2\Big(C_{2}(i)\Big), \cdots\cdots, \lambda_{\min}^2\Big(C_{M}(i)\Big)\Big\} $, $ Q\theta = Q(\theta_{1}, \theta_{2}, \cdots\cdots, \theta_{N})^{T} = (\lambda_{1}, \lambda_{2}, \cdots\cdots, \lambda_{N})^{T} > 0 $, $ \tilde{\lambda} = \min\limits_{1\leq i\leq N}\Big\{\frac{\lambda_{i}}{\theta_{i}}\Big\} > 0 $, $ \zeta_{i} = 0.5q\Big[(2-q)d_{i}-\sigma_{i}\Big]\theta_{i}-\lambda_{i} $, $ \tilde{\zeta} = \max\limits_{1\leq i\leq N}\Big\{\frac{\zeta_{i}}{\theta_{i}}\Big\} > 0 $, $ \theta _{m} = \underset{i\in S}{ \min}\; \{\theta _{i}\} $, $ \theta _{M} = \underset{i\in S}{ \max}\; \{\theta _{i}\} $, and $ K_{0} = \frac{\theta_{M}e^{N_{0}(\tilde{\lambda}+\tilde{\zeta})}}{\theta_{m}} $. In other words, the trivial solution of Eq (2.7) is $ q $th moment exponentially stable.

    Proof. Noting that $ Q $ is an M-matrix, there exists a vector $ \theta = (\theta_{1}, \theta_{2}, \cdots\cdots, \theta_{N})^{T} > 0 $ such that $ Q\theta = (\lambda_{1}, \lambda_{2}, \cdots\cdots, \lambda_{N})^{T} > 0 $. Choose the Lyapunov function $ V(z(t), t, i) = \theta _{i} |z (t)|^{q}, i\in \mathbb{S} $. Obviously, $ \theta _{m} |z (t)|^{q}\le V(z (t), t, i)\le \theta _{M}|z (t)|^{q} $. By utilizing the Itô formula, when $ t\in [t_{i}, s_{i}), i = 0, 1, 2, \cdots\cdots $, and we have that

    $ LV(z(t),t,i)=qθi|z(t)|q2zT(t)[˜G(z(t),t,i)+Hh=1ρh(i)(Bh(i)Υ)z(t)]+12qθi|z(t)|q2|˜β(z(t),t,i)|212q(2q)θi|z(t)|q4|zT(t)˜β(z(t),t,i)|2+Nj=1ψijθj|z(t)|q.
    $
    (3.2)

    According to Assumption 2, the following inequalities are calculated:

    $ qθi|z(t)|q2zT(t)ˆG(z(t),t,i)qθi|z(t)|q1|ˆG(z(t),i,t)|qθili|z(t)|q,
    $
    (3.3)
    $ qθi|z(t)|q2zT(t)(Hh=1ρh(i)(Bh(i)Υ)z(t))qθi(Hh=1ρh(i)λmax([Bh(i)Υ]s))|z(t)|q,
    $
    (3.4)
    $ |z(t)|q2|˜β(z(t),t,i)|2max(λmax(CT1(i)C1(i)),λmax(CT2(i)C2(i)),,λmax(CTM(i)CM(i)))|z(t)|q=σi|z(t)|q,
    $
    (3.5)
    $ |z(t)|q4|zT(t)˜β(z(t),t,i)|2min{λ2min(C1(i)),λ2min(C2(i)),,λ2min(CM(i))}|z(t)|q=di|z(t)|q.
    $
    (3.6)

    Substituting inequalities (3.3)–(3.6) into equality (3.2) yields that

    $ LV(z(t),t,i)[qθili+qθi(Hh=1ρh(i)λmax([Bh(i)Υ]s))+12qθiσi12q(2q)θidi+Nj=1ψijθj]|z(t)|q.
    $
    (3.7)

    Since $ Q\theta = (\lambda_{1}, \lambda_{2}, \cdots, \lambda_{N})^{T} > 0 $, when $ t\in [t_{i}, s_{i}), i = 0, 1, 2, \cdots\cdots $, we have that

    $ LV(z(t),t,i)λi|z(t)|qλiθiV(z(t),i,t)min1iN{λiθi}V(z(t),i,t)=˜λV(z(t),i,t).
    $
    (3.8)

    Moreover, when $ t\in [s_{i}, t_{i+1}), i = 0, 1, 2, \cdots\cdots $, it can be calculated that

    $ LV(z(t),t,i)[0.5q((2q)d(i)σ(i))θiλi]|z(t)|qζiθiV(z(t),t,i)max1iN{ζiθi}V(z(t),i,t)=˜ζV(z(t),i,t).
    $
    (3.9)

    Particularly, according to inequality (3.8), if $ t\in[t_{0}, s_{0}), t_{0} = 0 $, we can infer that

    $ EV(z(t),t,r(t))θME|z0|qe˜λt.
    $
    (3.10)

    Moreover, combining inequalities (3.9) and (3.10), it can be deduced that

    $ EV(z(t),t,r(t))θME|z(s0)|qe˜ζ(ts0)θME|z0|qe˜λ(s0t0)+˜ζ(ts0).
    $
    (3.11)

    Repeating the iteration leads to the following result,

    $ EV(z(t),t,r(t))θME|z0|qe˜λi1j=0(sjtj)+˜ζij=1(tjsj1)˜λ(tti),t[ti,si),
    $
    (3.12)

    and

    $ EV(z(t),t,r(t))θME|z0|qe˜λij=0(sjtj)+˜ζij=1(tjsj1)+˜ζ(tsi),t[si,ti+1),
    $
    (3.13)

    which means that

    $ EV(z(t),t,r(t))θME|z0|qe˜λNc(0,t)+˜ζNr(0,t).
    $
    (3.14)

    In light of Assumption 1, one has that $ N_{c}(t, 0)\ge \alpha t-N_0 $ and $ N_{r}(t, 0)\ge (1-\alpha) t+ N_0 $. Consequently,

    $ EV(z(t),t,r(t))θME|z0|qe˜λ[αtN0]+˜ζ[(1α)t+N0]eN0(˜λ+˜ζ)θME|z0|qe[˜λα˜ζ(1α)]t.
    $
    (3.15)

    Therefore, one has that

    $ E|z(t)|qK0E|z0|qeϵt,
    $
    (3.16)

    where $ K_{0} = \frac{\theta_{M}e^{N_{0}(\tilde{\lambda}+\tilde{\zeta})}}{\theta_{m}} $, and $ \epsilon = \tilde{\lambda}\alpha-\tilde{\zeta}(1-\alpha) > 0 $.

    Lemma 3.2. Let $ q\in(0, 1) $ and Assumption 2 hold. The following estimations are presented:

    $ E|η(t)|2|η0|2e(2ˆl+ˆσ)t,
    $
    (3.17)
    $ E|η(t)η(μ(t))|2K1(τ)e(2ˆl+ˆσ)t,
    $
    (3.18)
    $ E|η(t)z(t)|q|η0|q[4(ˆl+ˆσ)K1(τ)2ˆl+ˆσ]q2|η0|qe(ˆl+ˆσ)qt,
    $
    (3.19)

    where $ \hat{l} = \max\limits_{1\le i\le N}\left\{l_{i}+\sum\limits_{h = 1}^{H}\rho_{h}(i)\lambda_{\max}\left(\left [ B^{h}(i)\otimes\Upsilon \right ] _{s} \right) \right\} $, $ \hat{\sigma} = \max\limits_{0\le i\le N}\left\{\sigma _{i}\right\} $,

    $ \sigma_{i} = \max\left\{\lambda_{\max}\left(C_{1}^{T}(i)C_{1}(i)\right), \lambda_{\max}\left(C_{2}^{T}(i)C_{2}(i)\right), \cdots, \lambda_{\max}\left(C_{M}^{T}(i)C_{M}(i)\right)\right\} $,

    $ K_{1}(\tau) = \left\{4\tau \left[\left(\max\limits_{1\le i\le N}\left\{l_{i}\right\}\right)^{2}+H\max\limits_{1\le i\le N}\left \{ \sum\limits_{h = 1}^{H}\rho _{h}^{2}(i)\lambda_{\max}\left(\left (B^{h}(i)\otimes\Upsilon\right)^{T}\left(B^{h}(i)\otimes\Upsilon\right)\right)\right\}\right]+2\hat{\sigma}\right\}\tau $.

    Proof. Error system (2.6) in matrix form can be formulated as

    $ dη(t)=[˜G(η(t),t,r(t))+Hh=1ρh(r(t))(Bh(r(t))Υ)η(t)]dt+I(t)ˆβ(η(μ(t)),t,r(t))dB(t),
    $
    (3.20)

    where $ I(t) = \left\{1,t[ti,si),0,t[si,ti+1),i=0,1,2,.

    \right. $ By employing the Itô formula, one has that

    $ E|η(t)|2=|η0|2+Et0{2ηT(s)[˜G(η(s),s,r(s))+Hh=1ρh(r(s))(Bh(r(s))Υ)η(s)]+|ˆβ(η(μ(s)),s,r(s))|2I(s)}ds|η0|2+t02max1iN{li+Hh=1ρh(i)λmax([Bh(i)Υ]s)}E|η(s)|2ds+max1iN{σi}t0E|η(μ(s))|2ds.
    $
    (3.21)

    Moreover, it can be inferred that

    $ sup0stE|η(s)|2|η0|2+(2ˆl+ˆσ)t0(sup0stE|η(s)|2)ds.
    $
    (3.22)

    Applying the Gronwall inequality yields that

    $ sup0stE|η(s)|2|η0|2e(2ˆl+ˆσ)t.
    $
    (3.23)

    Accordingly, it can be concluded that

    $ sup0stE|η(s)|2|η0|2e(2ˆl+ˆσ)t,
    $
    (3.24)

    which indicates that

    $ E|η(t)|2|η0|2e(2ˆl+ˆσ)t.
    $
    (3.25)

    Obviously, it follows from Eq (3.20) that

    $ η(t)η(μ(t))=tμ(t)[˜G(η(s),s,r(s))+Hh=1ρh(r(s))(Bh(r(s))Υ)]η(s)ds
    $
    (3.26)
    $ +tμ(t)I(s)ˆβ(η(μ(s)),s,r(s))dw(s).
    $
    (3.27)

    Together with the Hölder inequality, it can be calculated as

    $ E|η(t)η(μ(t))|22E|tμ(t)[˜G(η(s),s,r(s))+Hh=1ρh(r(s))(Bh(r(s))Υ)]η(s)ds|2+2Etμ(t)|ˆβ(η(μ(s)),s,r(s))|2ds4τ[(max1iN{li})2+Hmax1iN{Hh=1ρ2h(i)λmax((Bh(i)Υ)T(Bh(i)Υ))}]tμ(t)E|η(s)|2ds+2ˆσtμ(t)E|η(μ(s))|2dsK1(τ)τtμ(t)sup0stE|η(s)|2duK1(τ)|η0|2e(2ˆl+ˆσ)t.
    $
    (3.28)

    On the other hand, by utilizing the Itô formula again, one gets that

    $ E|η(t)z(t)|2=2Et0{[η(s)z(s)]T[ˆG(η(s),s,r(s))ˆG(z(s),s,r(s))+Hh=1ρh(r(s))(Bh(r(s))Υ)(η(s)z(s))]+|ˆβ(η(μ(s)),s,r(s))ˆβ(η(s),s,r(s))|2I(s)}ds2[max1iN{li}+max1iN{Hh=1ρh(i)λmax([Bh(i)Υ]s)}]t0E|η(s)z(s)|2ds+max1iN{λmax(CT1(i)C1(i)),λmax(CT2(i)C2(i)),,λmax(CTM(i)CM(i))}×t0E|η(μ(s))z(s)|2ds2ˆlt0E|η(s)z(s)|2ds+ˆσt0E|η(μ(s)z(s)|2ds(2ˆl+2ˆσ)t0E|η(s)z(s)|2ds+2ˆσt0E|η(μ(s)η(s)|2ds.
    $
    (3.29)

    Recalling assertion (3.18), it can be deduced that

    $ t0E|η(μ(s)η(s)|2dsK1(τ)|η0|2t0e(2ˆl+ˆσ)sdsK1(τ)|η0|22ˆl+ˆσ[e(2ˆl+ˆσ)t1].
    $
    (3.30)

    Substituting inequality (3.30) into inequality (3.29) leads to

    $ E|η(t)z(t)|2(2ˆl+2ˆσ)t0E|η(s)z(s)|2ds+2ˆσK1(τ)|η0|22ˆl+ˆσ[e(2ˆl+ˆσ)t1].
    $
    (3.31)

    Subsequently, by using the Gronwall inequality [38], it can be readily obtained that

    $ E|η(t)z(t)|22ˆσK1(τ)|η0|22ˆl+ˆσ[e(2ˆl+ˆσ)t1]+4ˆσ(ˆl+ˆσ)K1(τ)|η0|22ˆl+ˆσt0e2(ˆl+ˆσ)(ts)[e(2ˆl+ˆσ)s1]ds[4(ˆl+ˆσ)K1(τ)2ˆl+ˆσ]|η0|2e2(ˆl+ˆσ)t.
    $
    (3.32)

    Hence, in light of the Hölder inequality, one can further acquire the following assertion:

    $ E|η(t)z(t)|q(E|η(t)z(t)|2)q2[4(ˆl+ˆσ)K1(τ)2ˆl+ˆσ]q2|η0|qeq(ˆl+ˆσ)t.
    $
    (3.33)

    Theorem 3.1. Suppose that all of the conditions in Lemma 3.1 are satisfied. Let $ \xi\in (0, 1) $. $ \tau ^{\ast} > 0 $ is the unique positive root of the following equation:

    $ ξ+K2(τ)e(ˆl+ˆσ)q(τ+1ϵlog(K0ξ))=1,
    $
    (3.34)

    where $ \theta_{m} = \min\limits_{0\le i\le N}\left\{\theta_{i}\right\} $, $ \theta_{M} = \max\limits_{0\le i\le N}\left\{\theta_{i}\right\} $, $ K_{2}(\tau) = \left[\frac{4(\hat{l}+\hat{\sigma}) K_{1}(\tau)}{2\hat{l}+\hat{\sigma}}\right]^{\frac{q}{2}} $, and $ \epsilon, \; \theta_{i}, \; \hat{\rho}, \; \hat{\sigma}, \; K_{0} $, and $ K_{1} $ are the same as in Lemmas 3.1 and 3.2. If $ 0 < \tau < \tau^{\ast} $, then the multilayer complex network system (2.1) can achieve almost sure exponential synchronization with isolated node (2.3) via aperiodically intermittent noise control, i.e.,

    $ limtsup1tlog|η(t)|2ˆδ5q<0,a.s.
    $
    (3.35)

    Proof. For given positive constant $ \xi\in (0, 1) $, we can choose sufficiently large positive integer $ m_{0} $ such that

    $ log(K0ξ)ϵτm0<1+log(K0ξ)ϵτ.
    $
    (3.36)

    Combining the above inequality with Lemma 3.1, one gets that

    $ K0em0ϵτξ,E|z(m0τ)|qK0em0ϵτ|η0|qξ|η0|q.
    $
    (3.37)

    Noting that

    $ ξ+K2(τ)e(ˆl+ˆσ)qm0τ<ξ+K2(τ)e(ˆl+ˆσ)q(τ+1ϵlog(K0ξ))=1,
    $
    (3.38)

    there exists a constant $ \hat{\delta} > 0 $ such that

    $ ξ+K2(τ)e(ˆl+ˆσ)qm0τ<eˆδm0τ.
    $
    (3.39)

    For simplicity, let $ z_{m_{0}} = z(m_{0}\tau, \eta_{0}, r_{0}) $ and $ \eta_{m_{0}} = \eta(m_{0}\tau, \eta_{0}, r_{0}) $. Combing Lemma 3.2 and the elementary inequality $ (u+v)^{q} \le u^{q} +v^{q} $, $ u\ge 0, v\ge 0, q\in (0, 1) $, one can compute that

    $ E|ηm0|qE|zm0|q+E|ηm0zm0|qK0em0ϵτ|η0|q+|η0|qK2(τ)e(ˆl+ˆσ)qm0τ|η0|q[ξ+K2(τ)e(ˆl+ˆσ)qm0τ]|η0|qeˆδm0τ.
    $
    (3.40)

    When $ t\ge m\tau $, it can be deduced that

    $ E(|η2m0|q|Fm0τ)|ηm0|qeˆδm0τ|η0|qe2ˆδm0τ.
    $
    (3.41)

    Similarly, one can derive that

    $ E|ηim0|qE|η(i1)m0|qeˆδm0τE|ηm0|qeˆδ(i1)m0τ|η0|qeˆδim0τi=1,2,3,.
    $
    (3.42)

    Let $ \tau\in(0, \tau^{\ast}) $. It follows from Eq (3.14) that

    $ η(t)=η0+t0[˜G(η(s),s,r(s))+Hh=1ρh(r(s))(Bh(r(s))Υ)η(s)]ds+I(s)ˆβ(η(μ(s)),s,r(s))dB(s).
    $

    By exploiting the Hölder inequality and the Burkholder-Davis-Gundy inequality, it can be estimated that

    $ E(sup0tm0τ|η(t)|2)3|η0|2+3E{sup0tm0τ|t0˜G(η(s),s,r(s))+Hh=1ρh(r(s))(Bh(r(s))Υ))η(s)ds|2}+3E{sup0tm0τ|t0I(s)ˆβ(η(μ(s)),s,r(s))dB(s)|2}3|η0|2+6mτ[(max1iN{li})2+Hmax1iN{Hi=1ρ2h(i)λmax(Bh(i)Υ)T(Bh(i)Υ))}].E(sup0tm0τt0|η(s)|2ds)+12ˆσE(sup0tm0τt0|η(μ(s))|2ds)3|η0|2+6m0τ[(max1iN{li})2+Hmax1iN{Hi=1ρ2h(i)λmax((Bh(i)Υ)T(Bh(i)Υ))}+12ˆσ]m0τ0E(sup0ts|η(t)|2)ds.
    $
    (3.43)

    Moreover, by the Gronwall inequality, it can be attained that

    $ E(sup0tm0τ|η(t)|2)3|η0|2eK3(τ)m0τ,
    $
    (3.44)

    where $ K_{3}(\tau) = 6m_{0}\tau\bigg[(\max\limits_{1\leq i\leq N}\{l_i\})^2 +H\max\limits_{1\leq i\leq N}\bigg\{\sum_{i = 1}^H\rho_h^2(i)\lambda_{\max}\bigg((B^h(i)\otimes \Upsilon) ^{T} (B^h(i)\otimes\Upsilon)\bigg)\bigg\} +12\hat{\sigma}\bigg] $. By virtue of the Hölder inequality, one gets that

    $ E(sup0τtm0τ|η(t)|q)3q2|η0|qeqK3(τ)mτ2=K4(τ)|η0|q,
    $
    (3.45)

    where $ K_{4}(\tau) = 3^{\frac{q}{2} }e^{\frac{qK_{3}(\tau)m\tau}{2}} $. Accordingly, based on the time stationarity of stochastic differential equations, it can be derived that

    $ E(supim0τt(i+1)m0τ|η(t)|q|Fim0τ)K4(τ)|ηim0τ|q,i=0,1,2,.
    $
    (3.46)

    Combining (3.42) and (3.46) gives

    $ E(supim0τt(i+1)m0τ|η(t)|q)K4(τ)E|ηim0τ|qK4(τ)|η0|qeˆδim0τ.
    $
    (3.47)

    According to the Markov inequality, one gets that

    $ P{supim0τt(i+1)m0τ|η(t)|qe0.4m0ˆδiτ}K4(τ)|η0|qe0.6m0ˆδiτ,i=1,2,.
    $
    (3.48)

    As a result, by the Borel-Cantelli lemma, for almost all $ \omega\in\Omega $, there exists a stochastic integer $ i_{0} = i_{0}(\omega) $ satisfying

    $ supim0τt(i+1)m0τ|η(t)|q<e0.4m0ˆδiτ,
    $
    (3.49)

    for $ i^{*} > i_{0}(\omega) $. Hence, it can be concluded that

    $ limtsup1tlog|η(t)|2ˆδ5q<0,a.s.
    $
    (3.50)

    which signifies that the multilayer complex network system (2.1) realizes almost sure exponential synchronization with isolated node (2.3) via aperiodically intermittent noise control.

    Remark 1. Recently, stability and synchronization of multilayer complex networks have been discussed by developing various control approaches such as state feedback control, event-triggered control, impulsive control, and intermittent control in [12,14,16,21]. Compared with the previous findings, this paper has dealt with the almost sure exponential synchronization issue of MCNS by adopting aperiodically intermittent discrete observation noise control, where the stochastic noise can be viewed as control input in the diffusion term rather than the drift term [22]. In the considered control scheme, the average control rate has been integrated to depict the distributions of work/rest intervals of the control strategy from an overall perspective. In general, as the average control rate become larger, the control performance will be better.

    Remark 2. Different from the work about $ p $th moment exponential synchronization of network systems in [21,22], by utilizing M-matrix theory and various stochastic analysis techniques including the Itô formula, the Gronwall inequality, and the Borel-Cantelli lemma, almost sure exponential synchronization of multilayer coupled networks with Markovian switching has been analyzed via AIDONC. Although some sufficient conditions on almost sure exponential synchronization of MCNS have been given in [23], the control strategy AIDONC has not been adopted. Additionally, in [35,36], almost sure exponential synchronization and stabilization of neural networks have been investigated under AIDONC and the theoretical work has been extended to hybrid neural networks and complex networks [37,38,39]. In this paper, Markovian switching and multilayer interaction factors have been taken into account simultaneously, which make our network system more general.

    Remark 3. In this paper, we utilize the stochastic operator $ \mathcal{L}V $ and matrix multiplication to investigate the almost sure exponential synchronization of stochastic multilayer complex networks. Actually, fast algorithms about matrix multiplication are an interesting and hot topic which can be applied in contemporary Intel Xeon microprocessors. Particularly, in [40], a broad series of algorithms taking advantage of the efficiency of fast matrix multiplication algorithms in other mathematical and computer science operations have been reported. Extending the applications of these algorithms to complex networks remains a current challenge.

    In particular, when $ \tau = 0 $ and $ \; H = 1 $, the multilayer networks are reduced to the following single-layered networks.

    $ dδk(t)=[G(δk(t),t,r(t))+ρ(r(t))Mm=1bkm(r(t))Υδm(t)]dt+uk(t,r(t)),k=1,2,,M.
    $
    (3.51)

    Corollary 1. Let $ Q = \rm{diag}\{\chi_{1}, \chi_{2}, \cdots\cdots, \chi_{N}\}-\Psi $ denote an M-matrix, where $ q\in(0, 1) $, $ \chi_{i} = 0.5q\Big((2-q)d(i)-\sigma (i)\Big)-q l_{i}-q\Bigg(\rho(i)\lambda_{\mathop{\rm{max}}}\Big(\Big[B(i)\otimes\Upsilon\Big]_{s}\Big)\Bigg) $, and $ \Psi = (\psi_{ij})_{N\times N} $. Under Assumptions 1 and 2, if $ \epsilon = \tilde{\lambda}\alpha-\tilde{\zeta}(1-\alpha) > 0 $, complex network (3.51) can achieve almost sure exponential synchronization with the isolated node (2.3), where $ \sigma_{i} = \max\Big\{\lambda_{\max}\Big(C_{1}^{T}(i)C_{1}(i)\Big), \lambda_{\max}\Big(C_{2}^{T}(i)C_{2}(i)\Big), \cdots\cdots, \lambda_{\max}\Big(C_{M}^{T}(i)C_{M}(i)\Big)\Big\} $, $ d_{i} = \min\Big\{\lambda_{\min}^2\Big(C_{1}(i)\Big), \lambda_{\min}^2\Big(C_{2}(i)\Big), \cdots\cdots, \lambda_{\min}^2\Big(C_{M}(i)\Big)\Big\} $, $ Q\theta = Q(\theta_{1}, \theta_{2}, \cdots\cdots, \theta_{N})^{T} = (\lambda_{1}, \lambda_{2}, \cdots\cdots, \lambda_{N})^{T} > 0 $, $ \tilde{\lambda} = \min\limits_{1\leq i\leq N}\Big\{\frac{\lambda_{i}}{\theta_{i}}\Big\} > 0 $, $ \zeta_{i} = 0.5q\Big[(2-q)d_{i}-\sigma_{i}\Big]\theta_{i}-\lambda_{i} $, $ \tilde{\zeta} = \max\limits_{1\leq i\leq N}\Big\{\frac{\zeta_{i}}{\theta_{i}}\Big\} > 0 $, $ \theta _{m} = \underset{i\in S}{ \min}\; \{\theta _{i}\} $, $ \theta _{M} = \underset{i\in S}{ \max}\; \{\theta _{i}\} $, and $ K_{0} = \frac{\theta_{M}e^{N_{0}(\tilde{\lambda}+\tilde{\zeta})}}{\theta_{m}} $.

    In particular, when $ \theta = (\theta_{1}, \theta_{2}, \cdots\cdots, \theta_{n})^{T} = (1, 1, \cdots\cdots, 1)^{T} and\; H = 1 $, immediately, we acquire the following assertion.

    Corollary 2. Let all the conditions in Theorem 1 hold except $ (\lambda_{1}, \lambda_{2}, \cdots\cdots, \lambda_{n})^{T} = (1, 1, \cdots\cdots, 1)^{T} $, $ \chi_{i} = 0.5q\Big[(2-q)d_{i}-\sigma_{i}\Big]-q l_{i}-q\Bigg(\rho(i)\lambda_{\mathop{\rm{max}}}\Big(\Big[B(i)\otimes\Upsilon\Big]_{s}\Big)\Bigg) $, and $ \hat{l} = \max\limits_{1\le i\le N}\left\{l_{i}+\rho(i)\lambda_{\max}\left(\left [ B(i)\otimes\Upsilon \right ] _{s} \right) \right\} $. If $ 0 < \tau < \tau^{\ast} $, then the multilayer coupled network system (3.51) can achieve almost sure exponential synchronization with isolated node (2.3) via aperiodically intermittent noise control, i.e.,

    $ limtsup1tlog|η(t)|2ˆδ5q<0,a.s.
    $
    (3.52)

    Remark 4. Fractional differences can be seen as generalization or as an extension of classical calculus while fractional order derivatives are non-locally distributed, which can describe the memory and hereditary effects of complex processes and materials accurately. In recent years, many mathematical models in engineering and biological sciences have been proposed by using discrete or continuous fractional equations [41,42,43,44,45]. Particularly, in [41,42], some dynamic behaviors such as the existence and stability of solutions of worms in a wireless sensor model in the sense of the fractal fractional and fractional nabla difference COVID-19 model have been analyzed through function analysis and the Ulam-Hyers stability technique. Furthermore, some properties such as existence, uniqueness, controllability, stability of solutions of some equations including Hilfer fractional evolution equations, ABC-fuzzy-Volterra integro-differential equations, and coupled pantograph discrete fractional order difference equations have been discussed in [44,45]. Actually, fractional complex networks have been a current attractive topic. In the future, our work can be further extended to stochastic fractional complex networks.

    In order to exhibit the feasibility of our theoretical findings, a famous Chua's circuits network is considered here, which has been extensively applied to various areas as an essential nonlinear electronic oscillator model. Initially, the single Chua's circuit is described as follows:

    $ {˙V1(t)=1R1C1V1(t)+1R1C1V2(t)1C1Γ(V1(t));˙V2(t)=1R1C2V1(t)1R1C2V2(t)+1C2I3(t);˙I3(t)=1L0(V2(t)+R2I3(t));
    $
    (4.1)

    where $ V_{1}(t) $ and $ V_{2}(t) $ denote the voltages across the capacitors $ C_{1}^{*} $ and $ C_{2}^{*} $. Meanwhile, $ I_{3} $ signifies the current through the inductance $ L_{0} $, and $ R_{1} $ and $ R_{2} $ stands for the linear resistors. $ \Gamma(V_1(t)) $ denotes the current through the nonlinear resistor $ N_{R_{1}} $, which can be recasted as $ \Gamma(V_1(t)) = \varpi_{1}V_{1}+0.5(\varpi_{2}-\varpi_{1})(|V_{1}+1|-|V_{1}-1|), \varpi_{1} = -0.05, \varpi_{2} = -0.07 $. Soon afterward, we consider a multilayer coupled network of Chua's circuits with Markovian switching as follows:

    $ (dδk1(t)dδk2(t)dδk3(t))=(˜F11(r(t))˜F12(r(t))0˜F21(r(t))˜F22(r(t))˜F23(r(t))0˜F32(r(t))˜F33(r(t)))(δk1(t)δk2(t)δk3(t))dt+(˜Γk(r(t))00)dt+Hh=1Mm=1ρh(r(t))bhkm(r(t))Υδm(t)dt+uk(t,r(t)),
    $
    (4.2)

    where $ \delta_{k}(t) = (\delta_{k1}(t), \delta_{k2}(t), \delta_{k3}(t))^{T} = (V_{k1}, V_{k2}, I_{k3})^{T} $, $ (k = 1, 2, 3, 4, 5, 6) $, $ \tilde{F}_{11}(r(t)) = -\frac{1}{R_{1}(r(t))C_{1}^{*}(r(t))}, \tilde{F}_{12}(r(t)) = -\tilde{F}_{11}(r(t)) $, $ \tilde{F}_{21}(r(t)) = \frac{1}{R_{1}(r(t))C_{2}^{*}(r(t))}, \tilde{F}_{22}(r(t)) = -\frac{1}{R_{1}(r(t))C_{2}^{*}(r(t))}, \tilde{F}_{23}(r(t)) = \frac{1}{C_{2}^{*}(r(t))} $, $ \tilde{F}_{32}(r(t)) = -\frac{1}{L_{0}(r(t))} $, $ \tilde{F}_{33}(r(t)) = -\frac{R_{2}(r(t))}{L_{0}(r(t))} $, and $ \tilde{\Gamma}_{k}(r(t)) = -\frac{\Gamma(\delta_{k1}(t))}{C_{1}^{*}(r(t))} $. Accordingly, the isolated node of the network is expressed as

    $ ds(t)=G(s(t),t,r(t))dt.
    $
    (4.3)

    Meanwhile, controller $ u_{k}(t, r(t)) $ is designed as

    $ uk(t)={Ck(r(t))ηk(μ(t))dw(t),t[ti,si),0,t[si,ti+1),
    $
    (4.4)

    where $ \eta_{k}(t) = \delta_{k}(t)-s(t) $ and $ \mu(t) = \left [ t/\tau \right ]\tau, \tau > 0 $, and $ \tau $ denotes the duration between two consecutive observations. It is supposed that Markov jump $ r(t)\in\mathbb{S} = \{1, 2\} $ with generator

    $ Ψ=(4411).
    $

    Parameters $ M $ and $ H $ can be selected as $ M = 6 $ and $ H = 3 $. When $ r(t) = 1 $, let $ C_{1}^{*}(1) = 2 $, $ C_{2}^{*}(1) = 4 $, $ R_{1}^{*}(1) = 2 $, $ R_{2}^{*}(1) = 0.5 $, $ L_{0}^{*}(1) = 4 $, $ \rho_{1}(1) = 0.6 $, $ \rho_{2}(1) = 0.5 $, and $ \rho_{3}(1) = 0.4 $. When $ r(t) = 2 $, let $ C_{1}^{*}(2) = 4 $, $ C_{2}^{*}(2) = 2 $, $ R_{1}^{*}(2) = 2.5 $, $ R_{2}^{*}(2) = 1 $, $ L_{0}^{*}(2) = 5 $, $ \rho_{1}(2) = 0.4 $, $ \rho_{2}(2) = 0.6 $, and $ \rho_{3}(2) = 0.5 $. The values of six adjacency matrices $ B^{1}(1), B^{2}(1), B^{3}(1) $, $ B^{1}(2) $, $ B^{2}(2) $, and $ B^{3}(2) $ are given as

    $ B1(1)=[0.30.10.1000.100.20.100.10000.10.1000.10.100.30.100.10000.20.1000.10.100.2],B2(1)=[0.20.100.10000.30.100.10.10.100.20.10000.100.200.10.10000.100.100.10.100.3],B3(1)=[0.40.10.100.10.10.10.300.10.10000.20.100.10000.10.100.100.100.2000.100.100.2],B1(2)=[0.30.100.100.100.20.100.100.100.30.100.100.100.20.10000.100.100.100.10.100.3],B2(2)=[0.30.100.100.100.10.1000000.20.10.1000.10.10.30.100.10000.20.100.100.100.2],B3(2)=[0.200.100.100.10.3000.10.10.100.20.10000.100.20.100.100.100.30.100.100.100.2].
    $

    Furthermore, choose $ q = 0.5, \Upsilon = 0.5I, \alpha = 0.8, C_{k}(r(t)) = 2I, \xi = 0.9 $, and $ N_{0} = 0.01 $. By calculation, we can obtain $ l_{1} = 0.3215, l_{2} = 0.3561 $, $ \hat{l} = 0.3633, \hat{\sigma} = 4 $, $ \chi_{1} = 0.3349, \chi_{2} = 0.3184 $, $ \lambda_{1} = \lambda_{2} = 1, \theta_{1} = 3.1009, \theta_{2} = 3.1105, \hat{\lambda} = 0.3215, \hat{\zeta} = 0.1785 $, and $ \epsilon = \tilde{\lambda}\alpha-\tilde{\zeta}(1-\alpha) = 0.2215 > 0 $. According to $ \xi+K_{2}(\tau)e^{(\hat{l}+\hat{\sigma})q(\tau +\frac{1}{\epsilon}\log_{}{(\frac{K_{0}}{\xi }}))} = 1 $, the upper bound of duration time $ \tau $ is estimated as $ \tau < \tau* = 3.6213\times 10^{-5} $. Therefore, all of the conditions of Theorem 1 are satisfied, and the almost sure exponential synchronization between network Eq (4.2) and isolated node Eq (4.3) is realized. Figure 1 shows a right continuous Markov chain with initial data $ r(0) = 2 $. Meanwhile, Figures 24 illustrate the synchronization sample trajectories $ \eta_{i1}, \eta_{i2}, \eta_{i3}, (i = 1, 2, 3, 4, 5, 6) $ between network Eq (4.2) and isolated node Eq (4.3), respectively. It can be observed from Figures 24 that the error system is stable under aperiodically intermittent discrete observation noise control, which means that almost sure exponential synchronization is achieved. Therefore, the numerical simulations validate the effectiveness of the theoretical findings. By increasing the average control rate and the control intensity, the synchronization performance including the exponential convergent rate will be further improved. In comparison, this method is one discontinuous control, which can save cost and promote efficiency.

    Figure 1.  The Markov chain.
    Figure 2.  Synchronization errors $ \eta_{i1}(t) $ between network Eq (4.2) and isolated node Eq (4.3).
    Figure 3.  Synchronization errors $ \eta_{i2}(t) $ between network Eq (4.2) and isolated node Eq (4.3).
    Figure 4.  Synchronization errors $ \eta_{i3}(t) $ between network Eq (4.2) and isolated node Eq (4.3).

    In this paper, almost sure exponential synchronization of multilayer complex networks with Markovian switching via AIDONC is investigated. Different from the work about $ p $th moment exponential synchronization of network systems in [21,22], by utilizing M-matrix theory and various stochastic analysis techniques including the Itô formula, the Gronwall inequality, and the Borel-Cantelli lemma, some criteria on almost sure exponential synchronization of multilayer coupled networks with Markovian switching are constructed and the upper bound of the duration time is also estimated. It is noted that the control strategy is based on aperiodically intermittent discrete observation noise, where the average control rate is integrated to depict the distributions of work/rest intervals of the control strategy from an overall perspective. Finally, some numerical simulations are exhibited to illustrate the effectiveness and feasibility of our analytical findings. Although the aperiodically intermittent discrete observation noise control strategy can reduce energy consumption and save control cost, the derived algebraic criteria on almost sure exponential synchronization generally are comparatively difficult to calculate. Meanwhile, the duration time is comparatively small, and its length needs to be extended by designing several optimization algorithms in order to demonstrate the feasibility in practise. In the future, based on the work in [46,47], the issue of stabilization and synchronization of coupled networks with semi-Markovian switching based on AIDONC will be further explored.

    Li Liu: conceptualization, investigation, writing-original draft; Yinfang Song: funding acquisition, methodology, supervision; Hong Yu: software; Gang Zhang: writing-review and editing. All authors have read and approved the final version of the manuscript for publication.

    This work was supported by the National Natural Science Foundation of China (62076039, 62273059).

    The authors declare no conflict of interest.

    [1] R. K. Ahuja, T. L. Mananti and J. B. Orlin, "Network Flows: Theory, Algorithms, and Applications," Prentice Hall, Englewood Cliffs, NJ, 1993.
    [2] D. Armbruster, P. Degond and C. Ringhofer, A model for the dynamics of large queuing networks and supply chains, SIAM J. on Applied Mathematics, 66 (2006), 896-920. doi: 10.1137/040604625
    [3] R. E. Burkard, K. Dlaska and B. Klinz, The quickest flow problem, Z. Oper. Research, 37 (1993), 31-58.
    [4] L. G. Chalmet, R. L. Francis and P. B. Saunders, Network models for building evacuation, Fire Technology, 18 (1982), 90-113. doi: 10.1007/BF02993491
    [5] W. Choi, H. W. Hamacher and S. Tufekci, Modeling of building evacuation problems by network flows with side constraints, European Journal of Operational Research, 35 (1988), 98-110. doi: 10.1016/0377-2217(88)90382-7
    [6] G. M. Coclite, M. Garavello and B. Piccoli, Traffic flow on a road network, SIAM Journal on Mathematical Analysis, 36 (2005), 1862-1886. doi: 10.1137/S0036141004402683
    [7] C. D'Apice, S. Göttlich, M. Herty and B. Piccoli, "Modeling, Simulation and Optimization of Supply Chains: A Continuous Approach," SIAM, Philadelphia, PA, 2010.
    [8] L. Fleischer and É. Tardos, Efficient continuous-time dynamic network flow algorithms, Operations Research Letters, 23 (1998), 71-80. doi: 10.1016/S0167-6377(98)00037-6
    [9] L. R. Ford and D. R. Fulkerson, Constructing maximal dynamic flows from static flows, Operations Research, 6 (1958), 419-433. doi: 10.1287/opre.6.3.419
    [10] A. Fügenschuh, S. Göttlich, M. Herty, A. Klar and A. Martin, A discrete optimization approach to large scale supply networks based on partial differential equations, SIAM Journal on Scientific Computing, 30 (2008), 1490-1507.
    [11] A. Fügenschuh, M. Herty, A. Klar and A. Martin, Combinatorial and continuous models for the optimization of traffic flows on networks, SIAM Journal on Optimization, 16 (2006), 1155-1176. doi: 10.1137/040605503
    [12] S. Göttlich, M. Herty and A. Klar, Network models for supply chains, Communications in Mathematical Sciences, 3 (2005), 545-559.
    [13] E. Graat, C. Midden and P. Bockholts, Complex evacuation; effects of motivation level and slope of stairs on emergency egress time in a sports stadium, Safety Science, 31 (1999), 127-141. doi: 10.1016/S0925-7535(98)00061-7
    [14] S. Gwynne, E. R. Galea, M. Owen, P. J. Lawrence and L. Filippidis, A review of the methodologies used in evacuation modelling, Fire and Materials, 23 (1999), 383-388. doi: 10.1002/(SICI)1099-1018(199911/12)23:6<383::AID-FAM715>3.0.CO;2-2
    [15] H. W. Hamacher, S. Heller, G. Köster and W. Klein, A Sandwich Approach for Evacuation Time Bounds, in "Pedestrian and Evacuation Dynamics" (eds. R.D. Peacock, E.D. Kuligowski, and J.D. Averill), Springer US, (2011), 503-514. doi: 10.1007/978-1-4419-9725-8_45
    [16] H. W. Hamacher, K. Leiner and S. Ruzika, Quickest Cluster Flow Problems, in "Pedestrian and Evacuation Dynamics" (eds. R.D. Peacock, E.D. Kuligowski, and J.D. Averill), Springer US, (2011), 327-336. doi: 10.1007/978-1-4419-9725-8_30
    [17] H. W. Hamacher and S. A. Tjandra, Eariest arrival flows with time dependent capacity for solving evacuation problems, in "Pedestrian and Evacuation Dynamics" (eds. M. Schreckenberger and S.D. Sharma), Springer, Berlin, (2002), 267-276.
    [18] H. W. Hamacher and S. A. Tjandra, Mathematical modelling of evacuation problems-a state of the art, in "Pedestrian and Evacuation Dynamics" (eds. M. Schreckenberger and S.D. Sharma), Springer, Berlin, (2002), 227-266.
    [19] D. Helbing, A mathematical model for the behavior of pedestrians, Behavioral Science, 36 (1991), 298-310. doi: 10.1002/bs.3830360405
    [20] M. Herty and A. Klar, Modeling, simulation, and optimization of traffic flow networks, SIAM Journal on Scientific Computing, 25 (2003), 1066-1087.
    [21] H. Holden and N. H. Risebro, A mathematical model of traffic flow on a network of unidirectional roads, SIAM Journal on Mathematical Analysis, 26 (1995), 999-1017. doi: 10.1137/S0036141093243289
    [22] S. P. Hoogendoorn and P. H. L. Bovy, Gas-kinetic modeling and simulaton of pedestrian flows, Transportation Research Record, (2000), 28-36. doi: 10.3141/1710-04
    [23] R. Hughes, A continuum theory for the flow of pedestrians, Transportation Research Part B, 36 (2002), 507-535. doi: 10.1016/S0191-2615(01)00015-7
    [24] C. Kirchner, M. Herty, S. Göttlich and A. Klar, Optimal control for continuous supply network models, Networks Heterogenous Media, 1 (2006), 675-688. doi: 10.3934/nhm.2006.1.675
    [25] H. Klüpfel, M. Schreckenberg and T. Meyer-König, "Models for Crowd Movement and Egress Simulation," Traffic and Granular Flow '03, (2005), 357-372.
    [26] A. Kneidl, M. Thiemann, A. Borrmann, S. Ruzika, H. W. Hamacher, G. Köster and E. Rank, Bidirectional Coupling of Macroscopic and Microscopic Approaches for Pedestrian Behavior Prediction, in "Pedestrian and Evacuation Dynamics" (eds. R.D. Peacock, E.D. Kuligowski, and J.D. Averill), Springer US, (2011), 459-470. doi: 10.1007/978-1-4419-9725-8_41
    [27] E. Köhler, K. Langkau and M. Skutella, "Time-Expanded Graphs for Flow-Dependent Transit Times," Lecture Notes in Computer Science, 2461 Springer, Berlin, (2002), 599-611.
    [28] E. Köhler and M. Skutella, Flows over time with load-dependent transit times, SIAM Journal on Optimization, 15 (2005), 1185-1202. doi: 10.1137/S1052623403432645
    [29] C. D. Laird, L. T. Biegler and B. G. van Bloemen Waanders, Real-time, large-scale optimization of water network systems using a subdomain approach, in "Real-Time PDE-Constrained Optimization" (eds. L.T. Biegler, O. Ghattas, M. Heinkenschloss, D. Keyes and B. G. van Bloemen Waanders), SIAM Series in Computational Science and Engineering, Philadelphia, PA, (2007), 289-306. doi: 10.1137/1.9780898718935.ch15
    [30] C. D. Laird, L. T. Biegler, B. G. van Bloemen Waanders and R. A. Bartlett, Contaminant source determination for water networks, Journal of Water Resources Planning and Management, 131 (2005), 125-134. doi: 10.1061/(ASCE)0733-9496(2005)131:2(125)
    [31] O. Østerby, Five ways of reducing the Crank-Nicolson oscillations, BIT Numerical Mathematics, 43 (2003), 811-822. doi: 10.1023/B:BITN.0000009942.00540.94
    [32] C. E. Pearson, Impulsive end condition for diffusion equation, Mathematics of Computation, 19 (1965), 570-576. doi: 10.1090/S0025-5718-1965-0193765-5
    [33] B. Rajewsky, "Strahlendosis und Strahlenwirkung," Thieme, Stuttgart, 1954.
    [34] G. Santos and B. Aguirre, "A Critical Review of Emergency Evacuation Simulation," Proceedings of Building Occupant Movement during Fire Emergencies, June 10-11, NIST/BFRL Publications Online, Gaithersburg, (2004), 27-52.
    [35] A. Schadschneider, W. Klingsch, H. Klüpfel, T. Kretz, C. Rogsch and A. Seyfried, Evacuation dynamics: Empirical results, modeling and applications, in "Encyclopedia of Complexity and System Science" (ed. B. Meyers), Springer, New York, (2009), 3142-3176.
    [36] J. G. Siek, L.-Q. Lee and A. Lumsdaine, "The Boost Graph Library: User Guide and Reference Manual (C++ In-Depth Series)," Addison-Wesley, Boston, 2001.
    [37] M. Skutella, "An Introduction to Network Flows Over Time," Research Trends in Combinatorial Optimization, Springer, Berlin, (2009), 451-482.
    [38] F. Southworth, "Regional Evacuation Modeling: A State-of-the-Art Review," ORNL/TAM-11740, Oak Ridge National Laboratory, Energy Division, Oak Ridge, 1991. doi: 10.2172/814579
    [39] IBM ILOG CPLEX Optimization Studio,. Cplex version 12 (2010) .
  • This article has been cited by:

    1. Aslı Akyol İnada, Samaneh Arman, Babak Safaei, A novel review on the efficiency of nanomaterials for solar energy storage systems, 2022, 55, 2352152X, 105661, 10.1016/j.est.2022.105661
    2. Muneesh Sethi, R.K. Tripathi, Birajashis Pattnaik, Sushil Kumar, Rohit Khargotra, Satish Chand, Abhishek Thakur, Recent developments in design of evacuated tube solar collectors integrated with thermal energy storage: A review, 2022, 52, 22147853, 1689, 10.1016/j.matpr.2021.11.324
    3. Berrin YILMAZ, Evren OSMA, FARKLI DUYULUR ISI DEPOLAMA MADDELERİNİN DEŞARJ PERİYODUNA ETKİLERİNİN DENEYSEL ARAŞTIRILMASI, 2023, 24, 2147-0308, 39, 10.59314/tujes.1378892
    4. Sally F. Naji, Hosham S. Anead, Khalid F. Sultan, 2024, 3092, 0094-243X, 050008, 10.1063/5.0199869
    5. Lukman Ahmed Omeiza, Muhammad Abid, Yathavan Subramanian, Anitha Dhanasekaran, Saifullah Abu Bakar, Abul Kalam Azad, Challenges, limitations, and applications of nanofluids in solar thermal collectors—a comprehensive review, 2023, 1614-7499, 10.1007/s11356-023-30656-9
  • Reader Comments
  • © 2011 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(5133) PDF downloads(259) Cited by(13)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog