Numerical analysis of biothermal-fluids and cardiac thermal pulse of abdominal aortic aneurysm

EYK Ng; Leonard Jun Cong Looi; EYK Ng; Leonard Jun Cong Looi

doi:10.3934/mbe.2022479

Mathematical Biosciences and Engineering

2022, Volume 19, Issue 10: 10213-10251. doi: 10.3934/mbe.2022479

Previous Article Next Article

Research article Special Issues

Numerical analysis of biothermal-fluids and cardiac thermal pulse of abdominal aortic aneurysm

EYK Ng ^,,
Leonard Jun Cong Looi

School of Mechanical and Aerospace Engineering, College of Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798

Academic Editor: Wenrui Hao

Received: 19 May 2022 Revised: 12 July 2022 Accepted: 12 July 2022 Published: 20 July 2022

Abdominal aortic aneurysms are serious and difficult to detect, conditions can be deadly if they rupture. In this study, the heat transfer and flow physics of Abdominal Aortic Aneurysm (AAA) were discussed and associated with cardiac cycle to illustrate the cardiac thermal pulse (CTP) of AAA. A CTP and infrared thermography (IRT) evaluation-based on AAA and abdomen skin surface detection method was proposed, respectively. Infrared thermography (IRT) is a promising imaging technique that may detect AAA quicker and cheaper than other imaging techniques (as biomarker). From CFD rigid-wall and FSI Analysis, the transient bioheat transfer effect resulted in a distinct thermal signature (circular thermal elevation) on the temperature profile of midriff skin surface, at both regular body temperature and supine position, under normal clinical temperature. However, it is important to note that thermography is not a perfect technology, and it does have some limitations, such as lack of clinical trials. There is still work to be done to improve this imaging technique and make it a more viable and accurate method for detecting abdominal aortic aneurysms. However, thermography is currently one of the most convenient technologies in this field, and it has the potential to detect abdominal aortic aneurysms earlier than other techniques. CTP, on the other hand, was used to examine the thermal physics of AAA. In CFD rigid-wall Analysis, AAA had a CTP that only responded to systolic phase at regular body temperature. In contrast, a healthy abdominal aorta displayed a CTP that responded to the full cardiac cycle, including diastolic phase at all simulated cases. Besides, the findings from FSI Analysis suggest the influence of numerical simulation techniques on the prediction of thermal physics behaviours of AAA and abdominal skin surface. Lastly, this study correlated the relationship between natural convective heat transfer coefficient with AAA and provided reference for potential clinical diagnostic using IRT in clinical implications.

Keywords:

Citation: EYK Ng, Leonard Jun Cong Looi. Numerical analysis of biothermal-fluids and cardiac thermal pulse of abdominal aortic aneurysm[J]. Mathematical Biosciences and Engineering, 2022, 19(10): 10213-10251. doi: 10.3934/mbe.2022479

Related Papers:

[1]	Muhammad Asif Zahoor Raja, Adeeba Haider, Kottakkaran Sooppy Nisar, Muhammad Shoaib . Intelligent computing knacks for infected media and time delay impacts on dynamical behaviors and control measures of rumor-spreading model. AIMS Biophysics, 2024, 11(1): 1-17. doi: 10.3934/biophy.2024001
[2]	Julia C. Arciero, Paola Causin, Francesca Malgaroli . Mathematical methods for modeling the microcirculation. AIMS Biophysics, 2017, 4(3): 362-399. doi: 10.3934/biophy.2017.3.362
[3]	Mario Lefebvre . A Wiener process with jumps to model the logarithm of new epidemic cases. AIMS Biophysics, 2022, 9(3): 271-281. doi: 10.3934/biophy.2022023
[4]	Ramesh Rijal, Hari Prasad Lamichhane, Kiran Pudasainee . Molecular structure, homo-lumo analysis and vibrational spectroscopy of the cancer healing pro-drug temozolomide based on dft calculations. AIMS Biophysics, 2022, 9(3): 208-220. doi: 10.3934/biophy.2022018
[5]	Dafina Xhako, Niko Hyka, Elda Spahiu, Suela Hoxhaj . Medical image analysis using deep learning algorithms (DLA). AIMS Biophysics, 2025, 12(2): 121-143. doi: 10.3934/biophy.2025008
[6]	T. M. Medvedeva, A. K. Lüttjohann, M. V. Sysoeva, G. van Luijtelaar, I. V. Sysoev . Estimating complexity of spike-wave discharges with largest Lyapunov exponent in computational models and experimental data. AIMS Biophysics, 2020, 7(2): 65-75. doi: 10.3934/biophy.2020006
[7]	Kieran Greer . New ideas for brain modelling 6. AIMS Biophysics, 2020, 7(4): 308-322. doi: 10.3934/biophy.2020022
[8]	Dina Falah Noori Al-Sabak, Leila Sadeghi, Gholamreza Dehghan . Predicting the correlation between neurological abnormalities and thyroid dysfunction using artificial neural networks. AIMS Biophysics, 2024, 11(4): 403-426. doi: 10.3934/biophy.2024022
[9]	Soumyajit Podder, Abhishek Mallick, Sudipta Das, Kartik Sau, Arijit Roy . Accurate diagnosis of liver diseases through the application of deep convolutional neural network on biopsy images. AIMS Biophysics, 2023, 10(4): 453-481. doi: 10.3934/biophy.2023026
[10]	Arjun Acharya, Madan Khanal, Rajesh Maharjan, Kalpana Gyawali, Bhoj Raj Luitel, Rameshwar Adhikari, Deependra Das Mulmi, Tika Ram Lamichhane, Hari Prasad Lamichhane . Quantum chemical calculations on calcium oxalate and dolichin A and their binding efficacy to lactoferrin: An in silico study using DFT, molecular docking, and molecular dynamics simulations. AIMS Biophysics, 2024, 11(2): 142-165. doi: 10.3934/biophy.2024010

Abstract

1. Introduction

Numerous studies on recurrent neural networks (RNNs), including those on bidirectional associative memory neural networks, cellular neural networks, and Hopfield neural networks, among others, have been proposed. In the meantime, neural networks (NNs) have been employed as a tool to address issues that have arisen in associative memory, signal processing, image processing, static image treatment, pattern recognition, and optimization. Additionally, it seems that time delays are crucial to networked control systems, chemical reactions, communication systems, etc. The authors of [1], [2] investigated time-delayed artificial neural network electronic implementations. Several studies have looked into different sorts of delays inside NNs in this context; see the papers [38]–[40] and references therein. The passivity theory is widely applied in many engineering problems. Indeed, it is intimately related to the circuit analysis which is a useful and significant tool to analyze the stability of nonlinear systems, signal processing and chaos control. Thus, it has been employed in various fields of science and engineering [7], [8], [44]–[48]. In [9], the authors proposed neural adaptive output feedback control based on passivity with adaptive set-point regulation of nonlinear uncertain non-negative and compartmental systems. In the last decade, great attention has been paid to the passivity analysis of RNNs with delay-independent [10], [49] and delay-dependent [11]–[14], [36]. In [13], the author studied the passivity analysis of Markovian jump neural network with leakage time varying delay, discrete and distributed time varying delays. We can see the discussion of the extended dissipativity of discrete-time neural networks NNs with time-varying delay in [55]. However, comparatively less interest has been shown towards the passivity analysis of mode dependent delay on neutral type neural networks (NNNs) with Markovian jumping (Mj). On the other hand, a Markov jump system is a special classes of a hybrid system. Indeed, they have great ability to model the dynamical systems and their application can be found in manufacturing systems, economic systems, network control systems, modeling production system, communication systems and so on. In the recent years, several results are reported on the stability analysis for neural networks with Markovian jump parameters, see the references [16]–[19], [41]–[43]. In [20], the authors studied global exponential estimates of delayed stochastic NNs with Markovian switching by constructing with positive definite matrices in stochastic Lyapunov functional which are dependent on the system mode and a triple-integral term. The ℋ_∞ synchronization issue for singularly perturbed coupled neural networks (SPCNNs) affected by both nonlinear constraints and gain uncertainties was explored in [57] using a novel double-layer switching regulation containing Markov chain and persistent dwell-time switching regulation (PDTSR). Convolutional Neural Networks (CNNs) are efficient tools for pattern recognition applications. More on this topic can be seen in [58], [60]. An exponential synchronization problem for the multi-weighted complex dynamical network (MCDN) with hybrid delays on a time scale is investigated in [61]. We can see the establishment of fixed-point and coincidence-point consequences in generalized metric spaces in [59]. The nonfragile ℋ_∞ synchronization issue for a class of discrete-time Takagi–Sugeno (T–S) fuzzy Markov jump systems was investigated in [56]. In [21], the authors have studied stochastic NNNs with mixed time-delays under adaptive synchronization. Additionally, the problem of state estimation of RNNs with Mj parameters and mixed delays based on mode-dependent approach was investigated in [22].

However, some of the researchers discussed the robust passive filtering for NNNs with delays in [25]. In [26], the authors investigated the global asymptotic stability of NNNs with delays by utilizing the Lyapunov-Krasovskii functional (LKF) and the linear matrix inequality approach. While employing the method of Lyapunov–Krasovskii functional, we necessarily need these three steps for the derivation of a global asymptotic stability criterion: constructing a Lyapunov–Krasovskii functional, estimating the derivative of the Lyapunov–Krasovskii functional, and formulating a global asymptotic stability criterion. You will get an overview of recent developments in each of the above steps if you refer to [54].The author in [52] studied passivity and exponential passivity for NNNs with various delays. In [51], the authors investigated the robust passivity analysis of mixed delayed NNs with distributed time-varying delays. The exponential passivity of discrete-time switched NNs with transmission delay was studied in [53]. Recently, the authors in [5] studied passivity analysis for NNNs with Mj parameters and time delay in the leakage term. New delay-dependent passivity conditions are derived in terms of LMIs with a proper construction of LKF, and it can be checked easily via standard numerical packages. However, triple and quadruple integrals have not been taken into account to derive the passivity conditions and, moreover, the mode-dependent time delays have not been included in [5]. Recently, a novel Lyapunov functional with some terms involving triple or quadruple integrals are taken into account to study the state estimation problem with mode-dependent approach in [22]. Motivated by the above discussion, the main purpose of this paper is to study the global passivity of Mj for NNNs with leakage and mode-dependent delay terms. By construction of a new LKF involving mode-dependent Lyapunov matrices, some sufficient conditions are derived in terms of LMIs. For the sake of illustration, a numerical example is given to demonstrate the usefulness and effectiveness of the presented results. Unlike previous results, we will introduce an improved Lyapunov–Krasovskii functional with triple and quadruple integrals for deriving the reported stability results in this paper. Based on this discussion, our technique not only provides different approach but also gives less conservative conditions than those studied in [5], [22]. The rest of this paper is organized as follows. The problem and some preliminaries are introduced in Section 2. In Section 3, the main results are stated and proved. Some sufficient conditions for global passivity results are developed here. In Section 4, an illustrative example is provided to demonstrate the effectiveness of the proposed criteria. We conclude the results of this paper in Section 5.

Notations: Throughout this paper the following notations are used: $ℝ^{n}$ - n-dimensional Euclidean space; $ℝ^{n \times n}$ - the set of all $n \times n$ real matrices; $d i a g (\dots)$ - a block diagonal matrix; I -the identity matrix with compatible dimensions; $C^{T}$ - the transpose of C; X and Y are symmetric matrices, where $X \geq Y$ (similarly $X > Y$ )- $X - Y$ is a positive semi-definite (similarly positive definite); $(Ω, ℱ, 𝒫)$ - a complete probability space with a natural filtration ${ℱ_{t}}_{t \geq 0}$ ; $E [\cdot]$ - expectation operator with respect to the given probability measure $𝒫$ ; $𝒞 ([- d,0]; ℝ^{n})$ - the family of continuously differentiable function; $‖ φ ‖ = m a x {m a x_{- τ \leq θ \leq 0} | φ (θ) |, m a x_{- d \leq θ \leq 0} | φ^{'} (θ) |}$ ; $𝒞_{ℱ_{0}}^{2} ([- d,0]; ℝ^{n})$ - the family of bounded $ℱ_{0}$ -measurable; $𝒞 ([- d,0]; ℝ^{n})$ -valued stochastic variables $ξ = {ξ (θ) : - d \leq θ \leq 0}$ such that $\int_{- d}^{0} E | ξ (θ {) |}^{2} d s < \infty$ ; *-the symmetric block in one symmetric matrix.

2. Problem formulation and preliminaries

Let ${\nabla (t), t \in ℤ^{+}}$ be right-continuous on $(Ω, ℱ, 𝒫)$ . Here $(Ω, ℱ, 𝒫)$ is a complete probability space with Markov chain, $\nabla (t)$ takes the values from a finite state space $𝒮 = {1,2, \dots, N}$ with generator $Γ = (π_{i j})_{N \times N}$ given by

$P {\nabla (t + Δ t) = j | \nabla (t) = i} = {\begin{array}{l} π_{i j} Δ t + o (Δ t), i \neq j, \\ 1 + π_{i i} Δ t + o (Δ t), i = j, \end{array}$

where $Δ t > 0$ and ${l i m}_{Δ t \to 0} \frac{o (Δ t)}{Δ t} = 0$ , $π_{i j} \geq 0$ $(i \neq j)$ is the transition rate from i to j, and $π_{i i} = - \sum_{j = 1, j \neq i}^{N} π_{i j}$ .

We consider Mode-dependent Markov jump $N N N s$ with mixed time-delays:

${\begin{array}{l} \dot{x} (t) & = & - 𝒜 (\nabla (t)) x (t - σ (\nabla (t))) + ℬ (\nabla (t)) ℋ (x (t)) + 𝒞 (\nabla (t)) ℋ (x (t - τ (t, \nabla (t)))) \\ + 𝒟 (\nabla (t)) \dot{x} (t - h (t, \nabla (t))) + ℰ (\nabla (t)) \int_{t - d (t, \nabla (t))}^{t} ℋ (x (s)) d s + u (t) \\ y (t) & = & ℋ (x (t)) \end{array}$ (2.1)

where $x (t) = [x_{1} (t), x_{2} (t), \dots, x_{n} (t {)]}^{T} \in ℝ^{n}$ is the state vector linked with n neurons. The diagonal matrix $𝒜 (\nabla (t)) = d i a g (A_{1} (\nabla (t)), A_{2} (\nabla (t)), \dots, A_{n} (\nabla (t)))$ has positive entries $A_{i} (\nabla (t)) > 0$ $(i = 1,2, \dots, n)$ . $ℬ (\nabla (t))$ , $𝒞 (\nabla (t))$ , $𝒟 (\nabla (t))$ , $ℰ (\nabla (t))$ are known appropriate dimensional constant matrices. Here the neuron activation function is $ℋ (x (t)) = [ℋ_{1} (x_{1} (t)), ℋ_{2} (x_{2} (t)), \dots, ℋ_{n} (x_{n} (t {))]}^{T}$ . $u (t)$ denotes a constant input. $τ (t, \nabla (t))$ , $h (t, \nabla (t))$ , $d (t, \nabla (t))$ are mode dependent discrete, neutral and distributed delays, respectively and $σ (\nabla (t))$ is the mode dependent leakage delay.

Throughout this paper, we assume the following.

Assumption 1. For any $j = 1,2, \dots, n$ , $ℋ_{j} (0) = 0$ and there exist constants ${\hat{l}}_{j}^{-}$ and ${\hat{l}}_{j}^{+}$ such that

${\hat{l}}_{j}^{-} \leq \frac{ℋ_{j} (γ_{1}) - ℋ_{j} (γ_{2})}{γ_{1} - γ_{2}} \leq {\hat{l}}_{j}^{+},$ (2.2)

where $γ_{1}, γ_{2} \in ℝ$ , and $γ_{1} \neq γ_{2}$ .

For the sake of convenience, we denote $𝒜 (\nabla (t) = i) = 𝒜_{i}, ℬ (\nabla (t) = i) = ℬ_{i}, 𝒞 (\nabla (t) = i) = 𝒞_{i}$ , $𝒟 (\nabla (t) = i) = 𝒟_{i}$ , $ℰ (\nabla (t) = i) = ℰ_{i}$ , respectively.

System (2.1) can be rewritten as

${\begin{array}{l} \dot{x} (t) & = & - 𝒜_{i} x (t - σ_{i}) + ℬ_{i} ℋ (x (t)) + 𝒞_{i} ℋ (x (t - τ_{i} (t))) + 𝒟_{i} \dot{x} (t - h_{i} (t)) \\ + ℰ_{i} \int_{t - d_{i} (t)}^{t} ℋ (x (s)) d s + u (t) \\ y (t) & = & ℋ (x (t)) \end{array}$ (2.3)

and the parameters associated with time delays are assumed to satisfy following:

$0 \leq τ_{i} (t) \leq τ_{i}, {\dot{τ}}_{i} (t) \leq τ_{µ_{i}}, 0 \leq h_{i} (t) \leq h_{i}, {\dot{h}}_{i} (t) \leq h_{µ_{i}}, 0 \leq d_{i} (t) \leq d_{i}, {\dot{d}}_{i} (t) \leq d_{µ_{i}}, σ_{i} > 0$ (2.4)

where $τ_{i}, h_{i}, d_{i}, τ_{µ_{i}}, h_{µ_{i}}$ and $d_{µ_{i}}$ are some real constants and $τ = \underset{i \in S}{m a x} {τ_{i}}$ , $h = \underset{i \in S}{m a x} {h_{i}}$ , $d = \underset{i \in S}{m a x} {d_{i}}$ , $σ = \underset{i \in S}{m a x} {σ_{i}}$ .

Now we can see a few necessary lemmas and a definition.

$(b - a) [\int_{a}^{b} y {(s)}^{T} M y (s) d s] \geq {[\int_{a}^{b} y (s) d s]}^{T} M [\int_{a}^{b} y (s) d s] .$

holds.

Proof. By Schur complement,

$[\begin{matrix} y {(s)}^{T} M y (s) & y {(s)}^{T} \\ y (s) & M^{- 1} \end{matrix}] \geq 0, s \in [a, b] .$

On integration from a to b yields, $[\begin{matrix} \int_{a}^{b} y {(s)}^{T} M y (s) d s & \int_{a}^{b} y {(s)}^{T} d s \\ \int_{a}^{b} y (s) d s & (b - a) M^{- 1} \end{matrix}]$ $\geq 0, s \in [a, b] .$ Now using the Schur complement on this inequality, we obtain our desired result. □

$\pm 2 x^{T} y \leq x^{T} ℳ x + y^{T} ℳ^{- 1} y .$

Lemma 3. [31] (Schur complement) Given $Ω_{1}, Ω_{2}$ and $Ω_{3}$ are constant matrices with appropriate dimensions, where $Ω_{1}, Ω_{2} > 0$ are symmetric matrices, then

$Ω_{1} + Ω_{3}^{T} Ω_{2}^{- 1} Ω_{3} < 0 \Leftrightarrow [\begin{matrix} Ω_{1} & Ω_{3}^{T} \\ * & - Ω_{2} \end{matrix}] < 0, o r [\begin{matrix} - Ω_{2} & Ω_{3} \\ * & Ω_{1} \end{matrix}] < 0.$

Definition 1. [32] If there exists a scalar $ν \geq 0$ such that $\forall t_{p} \geq 0$ and for all solutions of (2.1), the following inequality holds under zero initial conditions,

$2 \int_{0}^{t_{p}} E {y {(s)}^{T} u (s)} d s \geq - γ \int_{0}^{t_{p}} E {u {(s)}^{T} u (s)} d s,$ (2.5)

then the system (2.3) is said to be passive.

3. Main results

Now, we denote

${\hat{L}}_{1} = d i a g {{\hat{l}}_{1}^{-} {\hat{l}}_{1}^{+}, {\hat{l}}_{2}^{-} {\hat{l}}_{2}^{+}, \dots, {\hat{l}}_{m}^{-} {\hat{l}}_{m}^{+}}, a n d {\hat{L}}_{2} = d i a g {\frac{{\hat{l}}_{1}^{-} + {\hat{l}}_{1}^{+}}{2}, \frac{{\hat{l}}_{2}^{-} + {\hat{l}}_{2}^{+}}{2}, \dots, \frac{{\hat{l}}_{m}^{-} + {\hat{l}}_{m}^{+}}{2}} .$

$[\begin{matrix} X_{i} & N_{i} \\ X_{i} \end{matrix}] \geq 0,$ (3.1)

$\sum_{j = 1}^{N} π_{i j} V_{j} - V \leq 0,$ (3.2)

$\sum_{j = 1, j \neq i}^{N} π_{i j} 𝒢_{j} - 𝒢 \leq 0,$ (3.3)

$Φ = [\begin{matrix} Ω & Γ^{T} \\ - \frac{1}{π_{i j}} K_{j} \end{matrix}] < 0,$ (3.4)

where $𝒢_{j}$ in (3.3) respectively represents $Q_{j}$ , $W_{j}$ , $R_{j}$ , $S_{j}$ , $U_{j}$ , $X_{j}$ , $Y_{j}$ , $T_{j}$ , $Z_{j}$ , $L_{j}$ and correspondingly 𝒢 represents $Q$ , $W$ , $R$ , $S$ , $U$ , $X$ , $Y$ , $T$ , $Z$ , (e.q.,when 𝒢_j is $Q_{j}$ , $𝒢$ is $Q$ ) and

$Ω = (ϑ_{i, j})_{15 \times 15},$

$ϑ_{1,1} = - P_{i} 𝒜_{i} - 𝒜_{i}^{T} P_{i} + π_{i i} P_{i} + \sum_{j \neq i} π_{i j} P_{j} + \sum_{j \neq i} π_{i j} K_{j} + Q_{1 i} + τ Q_{1} + W_{i} + τ W + R_{i} + σ R$

$+ τ_{i} U_{i} + \frac{τ^{2}}{2} U - \frac{1}{τ_{i}} X_{i} - \frac{1}{h_{i}} Y_{i} - 2 T_{i} - 2 Z_{i} - 2 L_{i} - {\hat{L}}_{1} H_{1 i} - {\hat{L}}_{1} H_{3 i} - {\hat{L}}_{1} H_{4 i}$

$- {\hat{L}}_{1} H_{5 i} - {\hat{L}}_{1} H_{6 i} - 2 σ^{2} M, ϑ_{1,2} = - \frac{1}{τ_{i}} N_{i}^{T} + \frac{1}{τ_{i}} X_{i}, ϑ_{1,3} = \frac{1}{τ_{i}} N_{i}^{T} + {\hat{L}}_{2} H_{3 i},$

$ϑ_{1,4} = \frac{2}{τ_{i}} T_{i}, ϑ_{1,5} = P_{i} ℬ_{i} + Q_{2 i} + Q_{2} + {\hat{L}}_{2} H_{1 i} + J_{1} ℬ_{i}, ϑ_{1,6} = P_{i} 𝒞_{i} + J_{1} 𝒞_{i},$

$ϑ_{1,7} = - J_{1} 𝒜_{i} - 𝒜_{i}^{T} J_{1}^{T} + {\hat{L}}_{2} H_{5 i}, ϑ_{1,8} = 𝒜_{i}^{T} P_{i} 𝒜_{i} - π_{i i} P_{i} 𝒜_{i} - \frac{2}{σ_{i}} L_{i}, ϑ_{1,9} = 2 σ M,$

$ϑ_{1,10} = P_{i}, ϑ_{1,11} = P_{i} 𝒟_{i} + J_{1} 𝒟_{i}, ϑ_{1,12} = \frac{1}{h_{i}} Y_{i} + {\hat{L}}_{2} H_{4 i}, ϑ_{1,13} = \frac{2}{h_{i}} Z_{i},$

$ϑ_{1,14} = J_{1} ℰ_{i}, ϑ_{1,15} = - J_{1} + {\hat{L}}_{2} H_{6 i}, ϑ_{2,2} = - (1 - τ_{µ_{i}}) Q_{1 i} - {\hat{L}}_{1} H_{2 i} - \frac{2}{τ_{i}} X_{i} + \frac{1}{τ_{i}} N_{i}^{T},$

$ϑ_{2,3} = \frac{1}{τ_{i}} X_{i} - \frac{1}{τ_{i}} N_{i}^{T}, ϑ_{2,6} = - (1 - τ_{µ_{i}}) Q_{1 i} + {\hat{L}}_{2} H_{2 i}, ϑ_{3,3} = - W_{i} - \frac{1}{τ_{i}} X_{i} - H_{3 i},$

$ϑ_{4,4} = - \frac{1}{τ_{i}} U_{i} - \frac{2}{τ_{i}^{2}} T_{i}, ϑ_{5,5} = Q_{3 i} + Q_{3} + d_{i} V_{i} + \frac{d^{2}}{2} V - H_{1 i}, ϑ_{5,8} = - ℬ_{i}^{T} P_{i} 𝒜_{i},$

$ϑ_{5,15} = ℬ_{i}^{T} J_{2}^{T}, ϑ_{6,6} = - (1 - τ_{µ_{i}}) Q_{3 i} - H_{2 i}, ϑ_{6,8} = - 𝒞_{i}^{T} P_{i} 𝒜_{i}, ϑ_{6,15} = 𝒞_{i}^{T} J_{2}^{T},$

$ϑ_{7,7} = - R_{i} - H_{5 i}, ϑ_{7,15} = - 𝒜_{i}^{T} J_{2}^{T}, ϑ_{8,8} = π_{i i} 𝒜_{i}^{T} P_{i} 𝒜_{i} - \frac{2}{σ_{i}^{2}} L_{i}, ϑ_{8,10} = - 𝒜_{i}^{T} P_{i},$

$ϑ_{8,11} = - 𝒜_{i}^{T} P_{i} 𝒟_{i}, ϑ_{8,14} = - 𝒜_{i}^{T} P_{i} ℰ_{i}, ϑ_{9,9} = \sum_{j \neq i} π_{i j} 𝒜_{j}^{T} P_{j} 𝒜_{j} - 2 M, ϑ_{10,10} = - γ I,$

$ϑ_{10,15} = J_{2}^{T}, ϑ_{11,11} = - S_{i} (1 - h_{µ_{i}}), ϑ_{11,15} = 𝒟_{i}^{T} J_{2}^{T}, ϑ_{12,12} = - \frac{1}{h_{i}} Y_{i} - H_{4 i},$

$ϑ_{13,13} = - \frac{2}{h_{i}^{2}} Z_{i}, ϑ_{14,14} = - \frac{(1 - d_{µ_{i}})}{d_{i}} V_{i}, ϑ_{14,15} = ℰ_{i}^{T} J_{2}^{T}, ϑ_{15,15} = S_{i} + h S + τ_{i} X_{i} + \frac{τ^{2}}{2} X$

$+ h_{i} Y_{i} + \frac{h^{2}}{2} Y + \frac{τ_{i}^{2}}{2} T_{i} + \frac{τ^{3}}{6} T + \frac{h_{i}^{2}}{2} Z_{i} + \frac{h^{3}}{6} Z + \frac{σ_{i}^{2}}{2} L_{i} + \frac{σ^{3}}{6} L + \frac{σ^{4}}{2} Z - (J_{2} + J_{2}^{T}) - H_{6 i},$

$Γ^{T} = [{\underset{}{\underset{︸}{0.....0}}}_{8 e l e m e n t s} \sum_{j \neq i} π_{i j} {(𝒜_{j}^{T} P_{j})}^{T} {\underset{}{\underset{︸}{0.....0}}}_{15 e l e m e n t s}]^{T},$

and the other coefficients are zero.

Proof. Here, we consider LKF candidate:

$V (x_{t}, i, t) = \sum_{κ = 1}^{13} V_{κ} (x_{t}, i, t),$ (3.5)

where

$V_{1} (x_{t}, i, t) = {[x (t) - 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s]}^{T} P_{i} [x (t) - 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s],$

$V_{2} (x_{t}, i, t) = \int_{t - τ_{i} (t)}^{t} ζ^{T} (s) Q_{i} ζ (s) d s + \int_{- τ}^{0} \int_{t + θ}^{t} ζ^{T} (s) Q ζ (s) d s d θ,$

$V_{3} (x_{t}, i, t) = \int_{t - τ_{i}}^{t} x^{T} (s) W_{i} x (s) d s + \int_{- τ}^{0} \int_{t + θ}^{t} x^{T} (s) W x (s) d s d θ,$

$V_{4} (x_{t}, i, t) = \int_{t - σ_{i}}^{t} x^{T} (s) R_{i} x (s) d s + \int_{- σ}^{0} \int_{t + θ}^{t} x^{T} (s) R x (s) d s d θ,$

$V_{5} (x_{t}, i, t) = \int_{t - h_{i}}^{t} x^{T} (s) S_{i} x (s) d s + \int_{- h}^{0} \int_{t + θ}^{t} x^{T} (s) S x (s) d s d θ,$

$V_{6} (x_{t}, i, t) = \int_{- d_{i} (t)}^{0} \int_{t + θ}^{t} ℋ^{T} (x (s)) V_{i} ℋ (x (s)) d s d θ + \int_{- d}^{0} \int_{θ}^{0} \int_{t + β}^{t} ℋ^{T} (x (s)) V ℋ (x (s)) d s d β d θ,$

$V_{7} (x_{t}, i, t) = \int_{- τ_{i}}^{0} \int_{t + θ}^{t} x^{T} (s) U_{i} x (s) d s d θ + \int_{- τ}^{0} \int_{θ}^{0} \int_{t + β}^{t} x^{T} (s) U x (s) d s d β d θ,$

$V_{8} (x_{t}, i, t) = \int_{- τ_{i}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) X_{i} \dot{x} (s) d s d θ + \int_{- τ}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) X \dot{x} (s) d s d β d θ,$

$V_{9} (x_{t}, i, t) = \int_{- h_{i}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) Y_{i} \dot{x} (s) d s d θ + \int_{- h}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) Y \dot{x} (s) d s d β d θ,$

$V_{10} (x_{t}, i, t) = \int_{- τ_{i}}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) T_{i} \dot{x} (s) d s d β d θ + \int_{- τ}^{0} \int_{θ}^{0} \int_{β}^{0} \int_{t + α}^{t} {\dot{x}}^{T} (s) T \dot{x} (s) d s d α d β d θ,$

$V_{11} (x_{t}, i, t) = \int_{- h_{i}}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) Z_{i} \dot{x} (s) d s d β d θ + \int_{- h}^{0} \int_{θ}^{0} \int_{β}^{0} \int_{t + α}^{t} {\dot{x}}^{T} (s) Z \dot{x} (s) d s d α d β d θ,$

$V_{12} (x_{t}, i, t) = \int_{- σ_{i}}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) L_{i} \dot{x} (s) d s d β d θ + \int_{- σ}^{0} \int_{θ}^{0} \int_{β}^{0} \int_{t + α}^{t} {\dot{x}}^{T} (s) L \dot{x} (s) d s d α d β d θ,$

$V_{13} (x_{t}, i, t) = σ^{2} \int_{- σ}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) M \dot{x} (s) d s d β d θ,$

where $ζ^{T} (t) = [x^{T} (t), ℋ^{T} (x (t {))]}^{T} .$

From (2.1), we get

$ℒ V (x_{t}, i, t) = \sum_{κ = 1}^{13} ℒ V_{κ} (x_{t}, i, t),$ (3.6)

where

$ℒ V_{1} (x_{t}, i, t) = 2 {[x (t) - 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s]}^{T} P_{i} \frac{d}{d t} [x (t) - 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s]$

$+ \sum_{j = 1}^{N} π_{i j} {[x (t) - 𝒜_{j} \int_{t - σ_{j}}^{t} x (s) d s]}^{T} P_{j} [x (t) - 𝒜_{j} \int_{t - σ_{j}}^{t} x (s) d s],$

$\leq 2 {[x (t) - 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s]}^{T} P_{i} [- 𝒜_{i} x (t) + ℬ_{i} ℋ (x (t)) + 𝒞_{i} ℋ (x (t - τ_{i} (t)))$

$+ 𝒟_{i} \dot{x} (t - h_{i} (t)) + ℰ_{i} \int_{t - d_{i} (t)}^{t} ℋ (x (s)) d s + u (t)]$

$+ π_{i i} {[x (t) - 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s]}^{T} P_{i} [x (t) - 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s]$

$+ \sum_{j \neq i} π_{i j} [x^{T} (t) P_{j} x (t) + \int_{t - σ_{j}}^{t} x^{T} (s) d s 𝒜_{j} P_{j} K_{j}^{- 1} P_{j} 𝒜_{j} \int_{t - σ_{j}}^{t} x (s) d s$

$+ \int_{t - σ}^{t} x^{T} (s) d s 𝒜_{j}^{T} P_{j} 𝒜_{j} \int_{t - σ}^{t} x (s) d s],$

$ℒ V_{2} (x_{t}, i, t) \leq ζ^{T} (t) Q_{i} ζ (t) - ζ^{T} (t - τ_{i} (t)) Q_{i} ζ (t - τ_{i} (t)) (1 - τ_{µ_{i}}) + \sum_{j = 1}^{N} π_{i j} \int_{t - τ_{j} (t)}^{t} ζ^{T} (s) Q_{j} ζ (s) d s$

$+ τ ζ^{T} (t) Q ζ (t) - \int_{t - τ}^{t} ζ^{T} (s) Q ζ (s) d s,$

$ℒ V_{3} (x_{t}, i, t) = x^{T} (t) W_{i} x (t) - x^{T} (t - τ_{i}) W_{i} x (t - τ_{i}) + \sum_{j = 1}^{N} π_{i j} \int_{t - τ_{j}}^{t} x^{T} (s) W_{j} x (s) d s$

$+ τ x^{T} (t) W x (t) - \int_{t - τ}^{t} x^{T} (s) W x (s) d s,$

$ℒ V_{4} (x_{t}, i, t) = x^{T} (t) R_{i} x (t) - x^{T} (t - σ_{i}) R_{i} x (t - σ_{i}) + \sum_{j = 1}^{N} π_{i j} \int_{t - σ_{j}}^{t} x^{T} (s) R_{j} x (s) d s$

$+ σ x^{T} (t) R x (t) - \int_{t - σ}^{t} x^{T} (s) R x (s) d s,$

$ℒ V_{5} (x_{t}, i, t) \leq {\dot{x}}^{T} (t) S_{i} \dot{x} (t) - {\dot{x}}^{T} (t - h_{i} (t)) S_{i} \dot{x} (t - h_{i} (t)) (1 - h_{µ_{i}}) + \sum_{j = 1}^{N} π_{i j} \int_{t - h_{j} (t)}^{t} {\dot{x}}^{T} (s) S_{j} \dot{x} (s) d s$

$+ h {\dot{x}}^{T} (t) S \dot{x} (t) - \int_{t - h}^{t} {\dot{x}}^{T} (s) S \dot{x} (s) d s,$

$ℒ V_{6} (x_{t}, i, t) = d_{i} (t) ℋ^{T} (x (t)) V_{i} ℋ (x (t)) - (1 - d_{µ_{i}}) \int_{t - d_{i} (t)}^{t} ℋ^{T} (x (s)) V_{i} ℋ (x (s)) d s$

$+ \sum_{j = 1}^{N} π_{i j} \int_{- d_{j} (t)}^{0} \int_{t + θ}^{t} ℋ^{T} (x (s)) V_{j} ℋ (x (s)) d s d θ$

$+ \frac{d^{2}}{2} ℋ^{T} (x (t)) V ℋ (x (t)) - \int_{- d}^{0} \int_{t + θ}^{t} ℋ^{T} (x (s)) V ℋ (x (s)) d s d θ,$

$ℒ V_{7} (x_{t}, i, t) = τ_{i} x^{T} (t) U_{i} x (t) - \int_{t - τ_{i}}^{t} x^{T} (s) U_{i} x (s) d s + \sum_{j = 1}^{N} π_{i j} \int_{- τ_{j}}^{0} \int_{t + θ}^{t} x^{T} (s) U_{j} x (s) d s d θ$

$+ \frac{τ^{2}}{2} x^{T} (t) U x (t) - \int_{- τ}^{0} \int_{t + θ}^{t} x^{T} (s) U x (s) d s d θ,$

$ℒ V_{8} (x_{t}, i, t) = τ_{i} {\dot{x}}^{T} (t) X_{i} \dot{x} (t) - \int_{t - τ_{i}}^{t} {\dot{x}}^{T} (s) X_{i} \dot{x} (s) d s + \sum_{j = 1}^{N} π_{i j} \int_{- τ_{j}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) X_{j} \dot{x} (s) d s d θ$

$+ \frac{τ^{2}}{2} {\dot{x}}^{T} (t) X \dot{x} (t) - \int_{- τ}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) X \dot{x} (s) d s d θ,$

$ℒ V_{9} (x_{t}, i, t) = h_{i} {\dot{x}}^{T} (t) Y_{i} \dot{x} (t) - \int_{t - h_{i}}^{t} {\dot{x}}^{T} (s) Y_{i} \dot{x} (s) d s + \sum_{j = 1}^{N} π_{i j} \int_{- h_{j}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) Y_{j} \dot{x} (s) d s d θ$

$+ \frac{h^{2}}{2} {\dot{x}}^{T} (t) Y \dot{x} (t) - \int_{- h}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) Y \dot{x} (s) d s d θ,$

$ℒ V_{10} (x_{t}, i, t) = \frac{τ_{i}^{2}}{2} {\dot{x}}^{T} (t) T_{i} \dot{x} (t) - \int_{- τ_{i}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) T_{i} \dot{x} (s) d s d θ + \frac{τ^{3}}{6} {\dot{x}}^{T} (t) T \dot{x} (t)$

$+ \sum_{j = 1}^{N} π_{i j} \int_{- τ_{j}}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) T_{j} \dot{x} (s) d s d β d θ - \int_{- τ}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) T \dot{x} (s) d s d β d θ,$

$ℒ V_{11} (x_{t}, i, t) = \frac{h_{i}^{2}}{2} {\dot{x}}^{T} (t) Z_{i} \dot{x} (t) - \int_{- h_{i}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) Z_{i} \dot{x} (s) d s d θ + \frac{h^{3}}{6} {\dot{x}}^{T} (t) Z \dot{x} (t)$

$+ \sum_{j = 1}^{N} π_{i j} \int_{- h_{j}}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) Z_{j} \dot{x} (s) d s d β d θ - \int_{- h}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) Z \dot{x} (s) d s d β d θ,$

$ℒ V_{12} (x_{t}, i, t) = \frac{σ_{i}^{2}}{2} {\dot{x}}^{T} (t) L_{i} \dot{x} (t) - \int_{- σ_{i}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) L_{i} \dot{x} (s) d s d θ + \frac{σ^{3}}{6} {\dot{x}}^{T} (t) L \dot{x} (t)$

$+ \sum_{j = 1}^{N} π_{i j} \int_{- σ_{j}}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) L_{j} \dot{x} (s) d s d β d θ - \int_{- σ}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) L \dot{x} (s) d s d β d θ,$

$ℒ V_{13} (x_{t}, i, t) = \frac{σ^{4}}{2} {\dot{x}}^{T} (t) M \dot{x} (t) - σ^{2} \int_{- σ}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) M \dot{x} (s) d s d θ .$

Using the upper bounds of discrete, neutral, distributed time-varying delays, leakage delays, Lemma 2 and with $π_{i i} < 0,$ the following relationship is obtained

$\begin{array}{l} \sum_{j \neq i} (- 2 π_{i j} x^{T} (t) P_{j} 𝒜_{j} \int_{t - σ_{j}}^{t} x (s) d s) \\ \leq \sum_{j \neq i} π_{i j} (x^{T} (t) K_{j} x (t) + \int_{t - σ_{j}}^{t} x^{T} (s) d s 𝒜_{j}^{T} P_{j} K_{j}^{- 1} P_{j} 𝒜_{j} \int_{t - σ_{j}}^{t} x (s) d s) \\ \leq \sum_{j \neq i} π_{i j} (x^{T} (t) K_{j} x (t) + \int_{t - σ}^{t} x^{T} (s) d s 𝒜_{j}^{T} P_{j} K_{j}^{- 1} P_{j} 𝒜_{j} \int_{t - σ}^{t} x (s) d s), \end{array}$ (3.7)

$\begin{array}{l} \sum_{j \neq i} π_{i j} (\int_{t - σ_{j}}^{t} x^{T} (s) d s) 𝒜_{j}^{T} P_{j} 𝒜_{j} (\int_{t - σ_{j}}^{t} x (s) d s) \\ \leq \sum_{j \neq i} π_{i j} (\int_{t - σ}^{t} x^{T} (s) d s) 𝒜_{j}^{T} P_{j} 𝒜_{j} (\int_{t - σ}^{t} x (s) d s) . \end{array}$ (3.8)

Similarly,

$\begin{array}{l} \sum_{j = 1}^{N} π_{i j} \int_{t - τ_{j} (t)}^{t} ζ^{T} (s) Q_{j} ζ (s) d s \leq \sum_{j = 1, j \neq i}^{N} π_{i j} \int_{t - τ_{j} (t)}^{t} ζ^{T} (s) Q_{j} ζ (s) d s \\ \leq \sum_{j = 1, j \neq i}^{N} π_{i j} \int_{t - τ}^{t} ζ^{T} (s) Q_{j} ζ (s) d s \\ \leq \int_{t - τ}^{t} ζ^{T} (s) Q ζ (s) d s, \end{array}$ (3.9)

$\sum_{j = 1}^{N} π_{i j} \int_{t - τ_{j}}^{t} x^{T} (s) W_{j} x (s) d s \leq \int_{t - τ}^{t} x^{T} (s) W x (s) d s,$ (3.10)

$\sum_{j = 1}^{N} π_{i j} \int_{t - σ_{j}}^{t} x^{T} (s) R_{j} x (s) d s \leq \int_{t - σ}^{t} x^{T} (s) R x (s) d s,$ (3.11)

$\sum_{j = 1}^{N} π_{i j} \int_{t - h_{j} (t)}^{t} {\dot{x}}^{T} (s) S_{j} \dot{x} (s) d s \leq \int_{t - h}^{t} {\dot{x}}^{T} (s) S \dot{x} (s) d s,$ (3.12)

$\sum_{j = 1}^{N} π_{i j} \int_{- d_{j} (t)}^{0} \int_{t + θ}^{t} ℋ^{T} (x (s)) V_{j} ℋ (x (s)) d s d θ \leq - \int_{- d}^{0} \int_{t + θ}^{t} ℋ^{T} (x (s)) V ℋ (x (s)) d s d θ,$ (3.13)

$\sum_{j = 1}^{N} π_{i j} \int_{- τ_{j}}^{0} \int_{t + θ}^{t} x^{T} (s) U_{j} x (s) d s d θ \leq \int_{- τ}^{0} \int_{t + θ}^{t} x^{T} (s) U x (s) d s d θ,$ (3.14)

$\sum_{j = 1}^{N} π_{i j} \int_{- τ_{j}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) X_{j} \dot{x} (s) d s d θ \leq \int_{- τ}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) X \dot{x} (s) d s d θ,$ (3.15)

$\sum_{j = 1}^{N} π_{i j} \int_{- h_{j}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) Y_{j} \dot{x} (s) d s d θ \leq \int_{- h}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) Y \dot{x} (s) d s d θ,$ (3.16)

$\sum_{j = 1}^{N} π_{i j} \int_{- τ_{j}}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) T_{j} \dot{x} (s) d s d β d θ \leq \int_{- τ}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) T \dot{x} (s) d s d β d θ,$ (3.17)

$\sum_{j = 1}^{N} π_{i j} \int_{- h_{j}}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) Z_{j} \dot{x} (s) d s d β d θ \leq \int_{- h}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) Z \dot{x} (s) d s d β d θ,$ (3.18)

$\sum_{j = 1}^{N} π_{i j} \int_{- σ_{j}}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) L_{j} \dot{x} (s) d s d β d θ \leq \int_{- σ}^{0} \int_{θ}^{0} \int_{t + β}^{t} {\dot{x}}^{T} (s) L \dot{x} (s) d s d β d θ .$ (3.19)

By using Lemma 1, we get

$- (1 - d_{µ_{i}}) \int_{t - d_{i} (t)}^{t} ℋ^{T} (x (s)) V_{i} ℋ (x (s)) d s \leq - \frac{(1 - d_{µ_{i}})}{d_{i}} \int_{t - d_{i} (t)}^{t} ℋ^{T} (x (s)) d s V_{i} \int_{t - d_{i} (t)}^{t} ℋ (x (s)) d s,$

$- \int_{t - τ_{i}}^{t} x^{T} (s) U_{i} x (s) d s \leq - \frac{1}{τ_{i}} \int_{t - τ_{i}}^{t} x^{T} (s) d s U_{i} \int_{t - τ_{i}}^{t} x (s) d s,$

$- \int_{t - h_{i}}^{t} {\dot{x}}^{T} (s) Y_{i} \dot{x} (s) d s \leq - \frac{1}{h_{i}} \int_{t - h_{i}}^{t} {\dot{x}}^{T} (s) d s Y_{i} \int_{t - h_{i}}^{t} \dot{x} (s) d s,$

$- \int_{- τ_{i}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) T_{i} \dot{x} (s) d s d θ \leq - \frac{2}{τ_{i}^{2}} \int_{- τ_{i}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) d s d θ T_{i} \int_{- τ_{i}}^{0} \int_{t + θ}^{t} \dot{x} (s) d s d θ,$

$- \int_{- h_{i}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) Z_{i} \dot{x} (s) d s d θ \leq - \frac{2}{h_{i}^{2}} \int_{- h_{i}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) d s d θ Z_{i} \int_{- h_{i}}^{0} \int_{t + θ}^{t} \dot{x} (s) d s d θ,$

$- \int_{- σ_{i}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) Z_{i} \dot{x} (s) d s d θ \leq - \frac{2}{σ_{i}^{2}} \int_{- σ_{i}}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) d s d θ Z_{i} \int_{- σ_{i}}^{0} \int_{t + θ}^{t} \dot{x} (s) d s d θ,$

$\int_{- σ}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) M \dot{x} (s) d s d θ \leq - \frac{2}{σ^{2}} \int_{- σ}^{0} \int_{t + θ}^{t} {\dot{x}}^{T} (s) d s d θ Z_{i} \int_{- σ}^{0} \int_{t + θ}^{t} \dot{x} (s) d s d θ .$

Note that from (3.1) and using the reciprocally convex technique in [34], we obtain

$- \int_{t - τ_{i}}^{t} {\dot{x}}^{T} (s) X_{i} \dot{x} (s) d s \leq - \int_{t - τ_{i}}^{t - τ_{i} (t)} {\dot{x}}^{T} (s) X_{i} \dot{x} (s) d s - \int_{t - τ_{i} (t)}^{t} {\dot{x}}^{T} (s) X_{i} \dot{x} (s) d s$

$\leq - \frac{1}{τ_{i}} ϖ^{T} (t) [\begin{matrix} R_{i} & X_{i} \\ R_{i} \end{matrix}] ϖ (t),$

where $ϖ (t) = [x^{T} (t - τ_{i} (t)) - x^{T} (t - τ_{i}), x^{T} (t) - x^{T} (t - τ_{i} (t))]$ . For positive diagonal matrices $H_{1 i}, H_{2 i}, H_{3 i}, H_{4 i}, H_{5 i}, H_{6 i}$ , by Assumption 1, we get

${[\begin{matrix} x (t) \\ ℋ (x (t)) \end{matrix}]}^{T} [\begin{matrix} {\hat{L}}_{1} H_{1 i} & - {\hat{L}}_{2} H_{1 i} \\ * & H_{1 i} \end{matrix}] [\begin{matrix} x (t) \\ ℋ (x (t)) \end{matrix}] \leq 0,$ (3.20)

${[\begin{matrix} x (t - τ_{i} (t)) \\ ℋ (x (t - τ_{i} (t))) \end{matrix}]}^{T} [\begin{matrix} {\hat{L}}_{1} H_{2 i} & - {\hat{L}}_{2} H_{2 i} \\ * & H_{2 i} \end{matrix}] [\begin{matrix} x (t - τ_{i} (t)) \\ ℋ (x (t - τ_{i} (t))) \end{matrix}] \leq 0,$ (3.21)

${[\begin{matrix} x (t) \\ x (t - τ_{i}) \end{matrix}]}^{T} [\begin{matrix} {\hat{L}}_{1} H_{3 i} & - {\hat{L}}_{2} H_{3 i} \\ * & H_{3 i} \end{matrix}] [\begin{matrix} x (t) \\ x (t - τ_{i}) \end{matrix}] \leq 0,$ (3.22)

${[\begin{matrix} x (t) \\ x (t - h_{i}) \end{matrix}]}^{T} [\begin{matrix} {\hat{L}}_{1} H_{4 i} & - {\hat{L}}_{2} H_{4 i} \\ * & H_{4 i} \end{matrix}] [\begin{matrix} x (t) \\ x (t - h_{i}) \end{matrix}] \leq 0,$ (3.23)

${[\begin{matrix} x (t) \\ x (t - σ_{i}) \end{matrix}]}^{T} [\begin{matrix} {\hat{L}}_{1} H_{5 i} & - {\hat{L}}_{2} H_{5 i} \\ * & H_{5 i} \end{matrix}] [\begin{matrix} x (t) \\ x (t - σ_{i}) \end{matrix}] \leq 0,$ (3.24)

${[\begin{matrix} x (t) \\ \dot{x} (t) \end{matrix}]}^{T} [\begin{matrix} {\hat{L}}_{1} H_{6 i} & - {\hat{L}}_{2} H_{6 i} \\ * & H_{6 i} \end{matrix}] [\begin{matrix} x (t) \\ \dot{x} (t) \end{matrix}] \leq 0.$ (3.25)

Hence, for any matrices J₁, J₂ of appropriate dimensions, we get

$\begin{array}{l} 0 = 2 [x^{T} (t) J_{1} + {\dot{x}}^{T} (t) J_{2}] [- 𝒜_{i} x (t - σ_{i}) + ℬ_{i} ℋ (x (t)) + 𝒞_{i} ℋ (x (t - τ_{i} (t))) \\ + 𝒟_{i} \dot{x} (t - h_{i} (t)) + ℰ_{i} \int_{t - d_{i} (t)}^{t} ℋ (x (s)) d s + u (t) - \dot{x} (t)] . \end{array}$ (3.26)

Using (3.6) and adding (3.20)-(3.26), we have

$ℒ V (x_{t}, i, t) - 2 y^{T} (t) u (t) - γ u^{T} (t) u (t) \leq ℒ {η^{T} (t) Φ η (t)},$ (3.27)

where

$η^{T} (t) = [x^{T} (t) x^{T} (t - τ_{i} (t)) x^{T} (t - τ_{i}) \int_{t - τ_{i}}^{t} x^{T} (s) d s ℋ^{T} (x (t)) ℋ^{T} (x (t - τ_{i} (t)))$

$x^{T} (t - σ_{i}) \int_{t - σ_{i}}^{t} x^{T} (s) d s \int_{t - σ}^{t} x^{T} (s) d s u^{T} (t) {\dot{x}}^{T} (t - h_{i} (t)) x^{T} (t - h_{i})$

$\int_{t - h_{i}}^{t} x^{T} (s) d s \int_{t - d_{i} (t)}^{t} ℋ^{T} (x (s)) d s {\dot{x}}^{T} (t)] .$

Hence from equation (3.4) we have,

$ℒ V (x_{t}, i, t) - 2 y {(t)}^{T} u (t) - γ u {(t)}^{T} u (t) \leq 0.$

Now, to show the passivity of the delayed NNs in (2.3), we take

$J (t_{p}) = 𝔼 {\int_{0}^{t_{p}} [- γ u {(t)}^{T} u (t) - 2 y {(t)}^{T} u (t)] d t}$ (3.28)

where $t_{p} \geq 0$ .

From Dynkin's formula, we get

$𝔼 [\int_{0}^{t_{p}} ℒ V (x_{t}, i, t) d t] = 𝔼 [V (x_{t_{p}}, i, t_{p})] - 𝔼 [V (x_{0}, \nabla (0),0)] .$

Therefore,

$\begin{array}{l} J (t_{p}) = 𝔼 {\int_{0}^{t_{p}} [- γ u {(t)}^{T} u (t) - 2 y {(t)}^{T} u (t) + ℒ V (x_{t}, i, t)] d t} - 𝔼 [\int_{0}^{t_{p}} ℒ V (x_{t}, i, t) d t] \\ = 𝔼 {\int_{0}^{t_{p}} [- γ u {(t)}^{T} u (t) - 2 y {(t)}^{T} u (t) + ℒ V (x_{t}, i, t)] d t} \\ - 𝔼 [V (x_{t}, i, t)] + 𝔼 [V (x_{0}, \nabla (0),0)] . \end{array}$ (3.29)

By applying lemma 3 to (3.4), we have

$Φ < 0.$ (3.30)

Thus, if (3.30) holds, then $𝔼 [V (x_{t_{p}}, i, t_{p})] \geq 0$ and $V (x_{0}, \nabla (0),0) = 0$ holds with zero initial conditions. From (3.30), it follows that $J (t_{p}) \leq 0$ for any $t_{p} \geq 0$ , which implies (2.5) is satisfied and hence the delayed NNs (2.3) is locally passive.

Now, we prove the global passivity of the system.

By taking expectation of (3.27) and then integration from 0 to t we get,

$\int_{0}^{t} 𝔼 [ℒ V (x_{s}, r (s), s)] d s - 2 \int_{0}^{t} 𝔼 [y^{T} (s) u (s)] d s - γ \int_{0}^{t} 𝔼 [u^{T} (s) u (s)] d s \leq \int_{0}^{t} 𝔼 [η^{T} (s) Φ η (s)] d s .$

Then by Dynkin's formula,

$𝔼 [ℒ V (x_{t}, i, t)] - 𝔼 [ℒ V (x_{0}, \nabla (0),0)] - 2 \int_{0}^{t} 𝔼 [y^{T} (s) u (s)] d s - γ \int_{0}^{t} 𝔼 [u^{T} (s) u (s)] d s$

$\leq \int_{0}^{t} 𝔼 [η^{T} (s) Φ η (s)] d s .$

Hence,

$\begin{matrix} 𝔼 [ℒ V (x_{t}, i, t)] - \int_{0}^{t} 𝔼 [η^{T} (s) Φ η (s)] d s \leq 𝔼 [ℒ V (x_{0}, \nabla (0),0)] + 2 \int_{0}^{t} 𝔼 [y^{T} (s) u (s)] d s \\ + γ \int_{0}^{t} 𝔼 [u^{T} (s) u (s)] d s \\ < \infty, t \geq 0. \end{matrix}$ (3.31)

By Jenson's inequality and (3.6), we get

$\begin{array}{l} E {∥ 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s ∥}^{2} = E {[𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s]}^{T} [𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s] \\ \leq λ_{m a x} (𝒜_{i}^{2}) E {[\int_{t - σ_{i}}^{t} x (s) d s]}^{T} [\int_{t - σ_{i}}^{t} x (s) d s] \\ \leq \frac{λ_{m a x} (𝒜_{i}^{2})}{λ_{m i n} (R_{i})} {[\int_{t - σ_{i}}^{t} E x (s) d s]}^{T} R_{i} [\int_{t - σ_{i}}^{t} E x (s) d s] \\ \leq σ_{i} \frac{λ_{m a x} (𝒜_{i}^{2})}{λ_{m i n} (R_{i})} {\int_{t - σ_{i}}^{t} E x^{T} (s) R_{i} x (s) d s} \\ \leq σ_{i} \frac{λ_{m a x} (𝒜_{i}^{2})}{λ_{m i n} (R_{i})} E V_{4} (x_{t}, i, t) \\ \leq σ \frac{λ_{m a x} (𝒜_{i}^{2})}{λ_{m i n} (R_{i})} E V (x_{t}, i, t) \\ \leq σ \frac{λ_{m a x} (𝒜_{i}^{2})}{λ_{m i n} (R_{i})} E V (x_{0}, 𝒜 (0),0), t \geq 0. \end{array}$ (3.32)

Similarly, it follows from the definition of $V_{1} (x_{t}, i, t)$ that

$E {∥ x (t) - 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s ∥}^{2} = E {[𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s]}^{T} [𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s]$

$\leq \frac{E V_{1} (x_{t}, i, t)}{λ_{m i n} (P_{i})}$

$\leq \frac{E V (x_{t}, i, t)}{λ_{m i n} (P_{i})}$

$\leq \frac{E V (x_{0}, \nabla (0),0)}{λ_{m i n} (P_{i})}, t \geq 0.$

Hence, it can be obtained that

$\begin{array}{l} E {∥ x (t) ∥}^{2} = E {∥ x (t) - 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s + 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s ∥}^{2} \\ \leq 2 E {∥ A_{i} \int_{t - σ_{i}}^{t} x (s) d s ∥}^{2} + 2 E {∥ x (t) - 𝒜_{i} \int_{t - σ_{i}}^{t} x (s) d s ∥}^{2} \\ \leq 2 σ \frac{λ_{m a x} (𝒜_{i}^{2})}{λ_{m i n} (R_{i})} E V (x_{0}, \nabla (0),0) + 2 \frac{E V (x_{0}, \nabla (0),0)}{λ_{m i n} (P_{i})} < \infty, t \geq 0, \end{array}$ (3.33)

$\begin{array}{l} E V (x_{0}, \nabla (0),0) \\ \leq {2 λ_{\underset{i \in S}{m a x}} (P_{i}) (1 + σ_{i}^{2} \underset{i \in S}{m a x} 𝒜_{i}) + τ \underset{i \in S}{m a x} {λ_{m a x} (Q_{i})} + τ^{2} λ_{m a x} (Q) + τ \underset{i \in S}{m a x} {λ_{m a x} (W_{i})} \\ + τ^{2} λ_{m a x} (W) + σ \underset{i \in S}{m a x} {λ_{m a x} (R_{i})} + σ^{2} λ_{m a x} (R) + h \underset{i \in S}{m a x} {λ_{m a x} (S_{i})} + h^{2} λ_{m a x} (S) \\ + d^{2} \underset{i \in S}{m a x} {λ_{m a x} (V_{i})} + d^{3} λ_{m a x} (V) + τ^{2} \underset{i \in S}{m a x} {λ_{m a x} (U_{i})} + τ^{3} λ_{m a x} (U) + τ^{2} \underset{i \in S}{m a x} {λ_{m a x} (X_{i})} \\ + τ^{3} λ_{m a x} (X) + h^{2} \underset{i \in S}{m a x} {λ_{m a x} (Y_{i})} + h^{3} λ_{m a x} (Y) + τ^{3} \underset{i \in S}{m a x} {λ_{m a x} (T_{i})} + τ^{4} λ_{m a x} (T) \\ + h^{3} \underset{i \in S}{m a x} {λ_{m a x} (Z_{i})} + h^{4} λ_{m a x} (Z) + σ^{3} \underset{i \in S}{m a x} {λ_{m a x} (L_{i})} + σ^{4} λ_{m a x} (L) + σ^{5} λ_{m a x} (M)} < \infty . \end{array}$ (3.34)

From (3.33) and (3.34), we get that the solution of the system (2.3) is locally passive. Then the solutions $x (t) = x (t,0, φ)$ of system (2.3) is bounded on $[0, \infty)$ . The solution $x (t)$ on $[0, \infty)$ is uniformly continuous because $\dot{x} (t)$ is bounded on $[0, \infty)$ . Further, from the equation (3.31), the following holds:

$λ_{m i n} (Φ) \int_{0}^{t} E [x^{T} (s) x (s)] d s \leq 𝔼 [ℒ V (x_{t}, i, t)] - \int_{0}^{t} 𝔼 [η^{T} (s) Φ η (s)] d s$

$\leq 𝔼 [ℒ V (x_{0}, \nabla (0),0)] + 2 \int_{0}^{t} 𝔼 [y^{T} (s) u (s)] d s$

$+ γ \int_{0}^{t} 𝔼 [u^{T} (s) u (s)] d s$

$< \infty, t \geq 0.$

From Barbalats lemma [30], $𝔼 [{‖ x (t) ‖}^{2}] \to 0$ as $t \to \infty$ holds. Hence the proof is complete.

Remark 1. Recently, studies on passivity analysis for neural networks of neutral type with Markovian jumping parameters and time delay in the leakage term were conducted in [5]. By constructing proper Lyapunov–Krasovskii functional, new delay-dependent passivity conditions are derived in terms of LMIs and it can easily be checked using standard numerical packages. Moreover, it is well known that the passivity behaviour of neural networks is very sensitive to the time delay in the leakage term. Triple and quadruple integrals have not been taken into account to derive the passivity conditions in [5]. Mode-dependent time delays were not included in [5]. Very recently, a mode-dependent approach is proposed by constructing a novel Lyapunov functional, where some terms involving triple or quadruple integrals are taken into account to study the state estimation problem in [22]. Motivated by this reason, we have introduced improved Lyapunov–Krasovskii functional with triple and quadruple integrals for deriving the reported stability results in this paper. Based on this discussion, our results will give less conservative results than those studied in [5], [22].

4. Numerical example

In this section, a numerical example is provided to demonstrate the validity of the proposed theorems.

Example 1. Consider a 2-D Mode-dependent Markov jump NNNs with mixed time-delays (2.3) with the following parameters

$𝒜_{1} = [\begin{matrix} 8.4 & 0 \\ 0 & 9 \end{matrix}], 𝒜_{2} = [\begin{matrix} 7.8 & 0 \\ 0 & 8.5 \end{matrix}], ℬ_{1} = [\begin{matrix} - 0.21 & - 0.19 \\ - 0.24 & 0.1 \end{matrix}], ℬ_{2} = [\begin{matrix} 0.9 & - 0.9 \\ 0.5 & - 0.8 \end{matrix}], 𝒞_{1} = [\begin{matrix} - 0.09 & - 0.2 \\ 0.2 & 0.1 \end{matrix}], 𝒞_{2} = [\begin{matrix} 0.1 & 0.1 \\ 0.2 & 0.3 \end{matrix}], 𝒟_{1} = [\begin{matrix} - 0.2 & 0 \\ 0.2 & - 0.09 \end{matrix}], 𝒟_{2} = [\begin{matrix} 0.1 & 0 \\ 0.5 & - 0.1 \end{matrix}], ℰ_{1} = [\begin{matrix} - 0.5 & 0 \\ 0 & - 0.5 \end{matrix}], ℰ_{2} = [\begin{matrix} 0.1 & - 0.02 \\ - 0.2 & 0.07 \end{matrix}], {\hat{L}}_{1} = [\begin{matrix} 0 & 0 \\ 0 & 0 \end{matrix}], {\hat{L}}_{2} = [\begin{matrix} 0.25 & 0 \\ 0 & 0.25 \end{matrix}] .$

Take $ℋ_{1} (s) = ℋ_{2} (s) = t a n h (s)$ , $τ_{1} (t) = τ_{2} (t) = h_{1} (t) = h_{2} (t) = d_{1} (t) = d_{2} (t) = 0.1 c o s t + 0.4$ , $σ_{1} = σ_{2} = 0.1$ , $τ_{µ_{1}} = τ_{µ_{2}} = h_{µ_{1}} = h_{µ_{2}} = d_{µ_{1}} = d_{µ_{2}} = 0.1$ . $Γ = [\begin{matrix} - 7 & 7 \\ 6 & - 6 \end{matrix}] .$

By using the MATLAB LMI toolbox, we can obtain the following feasible solution for the LMIs (3.1)–(3.4):

$P_{1} = [\begin{matrix} 0.0157 & 0.0008 \\ 0.0008 & 0.0143 \end{matrix}], P_{2} = [\begin{matrix} 0.0068 & 0.0004 \\ 0.0004 & 0.0060 \end{matrix}], Q_{11} = [\begin{matrix} 0.2401 & 0.0007 \\ 0.0007 & 0.2592 \end{matrix}], Q_{12} = [\begin{matrix} 0.3770 & - 0.0192 \\ - 0.0192 & 0.3992 \end{matrix}], Q_{21} = [\begin{matrix} - 0.4101 & 0.0043 \\ 0.0201 & - 0.4457 \end{matrix}], Q_{22} = [\begin{matrix} - 0.6428 & 0.0773 \\ 0.0231 & - 0.6959 \end{matrix}], Q_{31} = [\begin{matrix} 1.3786 & 0.0133 \\ 0.0133 & 1.4063 \end{matrix}], Q_{32} = [\begin{matrix} 1.9101 & - 0.0284 \\ - 0.0284 & 1.9903 \end{matrix}], W_{1} = [\begin{matrix} 0.1039 & 0.0088 \\ 0.0088 & 0.1018 \end{matrix}], W_{2} = [\begin{matrix} 0.1517 & 0.0121 \\ 0.0121 & 0.1487 \end{matrix}], R_{1} = [\begin{matrix} 0.1758 & 0.0160 \\ 0.0160 & 0.1721 \end{matrix}], R_{2} = [\begin{matrix} 0.5675 & 0.0333 \\ 0.0333 & 0.5563 \end{matrix}], S_{1} = [\begin{matrix} 1.0817 & - 0.0075 \\ - 0.0075 & 0.7204 \end{matrix}], S_{2} = [\begin{matrix} 1.2231 & 0.0093 \\ 0.0093 & 0.9256 \end{matrix}], V_{1} = [\begin{matrix} 10.7113 & 0.5692 \\ 0.5692 & 9.9973 \end{matrix}], V_{2} = [\begin{matrix} 10.8483 & 0.5953 \\ 0.5953 & 10.2619 \end{matrix}], U_{1} = [\begin{matrix} 0.5120 & 0.0419 \\ 0.0419 & 0.5056 \end{matrix}], U_{2} = [\begin{matrix} 0.6238 & 0.0492 \\ 0.0492 & 0.6133 \end{matrix}], X_{1} = [\begin{matrix} 21.0470 & - 0.1153 \\ - 0.1153 & 21.2070 \end{matrix}], X_{2} = [\begin{matrix} 18.5642 & - 0.3139 \\ - 0.3139 & 18.6620 \end{matrix}], Y_{1} = [\begin{matrix} 27.5946 & - 0.0019 \\ - 0.0019 & 26.6557 \end{matrix}], Y_{2} = [\begin{matrix} 28.4115 & 0.1096 \\ 0.1096 & 27.4131 \end{matrix}], T_{1} = [\begin{matrix} 3.0928 & 0.0943 \\ 0.0943 & 3.0505 \end{matrix}], T_{2} = [\begin{matrix} 3.8182 & 0.0431 \\ 0.0431 & 3.7941 \end{matrix}], Z_{1} = [\begin{matrix} 3.0867 & 0.1080 \\ 0.1080 & 3.0591 \end{matrix}], Z_{2} = [\begin{matrix} 3.8647 & 0.0504 \\ 0.0504 & 3.8446 \end{matrix}], L_{1} = [\begin{matrix} 0.0068 & 0.0001 \\ 0.0001 & 0.0068 \end{matrix}], L_{2} = [\begin{matrix} 0.0059 & 0.0004 \\ 0.0004 & 0.0062 \end{matrix}], K_{1} = [\begin{matrix} 1.3461 & - 0.1410 \\ - 0.1410 & 1.3280 \end{matrix}], K_{2} = [\begin{matrix} 0.0919 & 0.0135 \\ 0.0135 & 0.0907 \end{matrix}], Q_{1} = [\begin{matrix} 4.0913 & - 0.1850 \\ - 0.1850 & 4.4495 \end{matrix}], Q_{2} = [\begin{matrix} - 6.6332 & 0.7686 \\ 0.3444 & - 7.3948 \end{matrix}], γ = 0.2272.$

Continuing in this way, the remaining feasible matrices are obtained. This shows that the given system (2.3) is globally passive in the mean square.

5. Conclusion

In this paper, passivity analysis of Markovian jumping NNNs with time delays in the leakage term is considered. Delay-mode-dependent passivity conditions are derived by taking the inherent characteristic of such kinds of NNs into account. An improved LKF, with the triple integral terms and quadruple integrals, is constructed and the results are derived in terms of linear matrix inequalities. The information of the mode-dependent of all delays have been taken into account in the constructed LKF and derived novel stability criterion. Theoretical results are validated through a numerical example.

References

[1]	F. A. Lederle, J. A. Freischlag, T. C. Kyriakides, F. T. Padberg, J. S. Matsumura, T. R. Kohler, et al., Outcomes following endovascular vs open repair of abdominal aortic aneurysm: A randomized trial, JAMA, 302 (2009), 1535-1542. https://doi.org/10.1001/jama.2009.1426 doi: 10.1001/jama.2009.1426
[2]	A. Chervu, G. P. Clagett, R. J. Valentine, S. I. Myers, P. J. Rossi, Role of physical examination in detection of abdominal aortic aneurysms, Surgery, 117 (1995), 454-457. https://doi.org/10.1016/S0039-6060(05)80067-4 doi: 10.1016/S0039-6060(05)80067-4
[3]	B. Keisler, C. Carter, Abdominal aortic aneurysm, Am. Fam. Phys., 91 (2015), 538-543.
[4]	P. M. Shaw, J. Loree, R. C. Gibbons, Abdominal Aortic Aneurysm, StatPearls Publishing, Treasure Island (FL), 2022.
[5]	T. Canchi, A. Saxena, E. Y. K. Ng, E. C. H. Pwee, S. Narayanan, Application of fluid-structure interaction methods to estimate the mechanics of rupture in asian abdominal aortic aneurysms, Bionanoscience, 8 (2018), 1035-1044. https://doi.org/10.1007/s12668-018-0554-z doi: 10.1007/s12668-018-0554-z
[6]	L. Smith-Burgess, Early identification and detection of abdominal aortic aneurysms, Nurs. Times, 113 (2017), 36-39.
[7]	D. C. Brewster, J. L. Cronenwett, J. W. Hallett, K. W. Johnston, W. C. Krupski, J. S. Matsumura, Guidelines for the treatment of abdominal aortic aneurysms: Report of a subcommittee of the Joint Council of the American Association for Vascular Surgery and Society for Vascular Surgery, J. Vasc. Surg., 37 (2003), 1106-1117. https://doi.org/10.1067/mva.2003.363 doi: 10.1067/mva.2003.363
[8]	K. C. Kent, R. M. Zwolak, N. N. Egorova, T. S. Riles, A. Manganaro, A. J. Moskowitz, et al., Analysis of risk factors for abdominal aortic aneurysm in a cohort of more than 3 million individuals, J. Vasc. Surg., 52 (2010), 539-548. https://doi.org/10.1016/j.jvs.2010.05.090 doi: 10.1016/j.jvs.2010.05.090
[9]	H. Bengtsson, D. Bergqvist, N. H. Sternby, Increasing prevalence of abdominal aortic aneurysms: A necropsy study, Eur. J. Surg., 158 (1992), 19-23.
[10]	F. L. Moll, J. T. Powell, G. Fraedrich, F. Verzini, S. Haulon, M. Waltham, et al., Management of abdominal aortic aneurysms clinical practice guidelines of the European society for vascular surgery, Eur. J. Vasc. Endovasc. Surg., 41 (2011), S1-S58. https://doi.org/10.1016/j.ejvs.2010.09.011 doi: 10.1016/j.ejvs.2010.09.011
[11]	J. Lieberg, L. L. Pruks, M. Kals, K. Paapstel, A. Aavik, J. Kals, Mortality after elective and ruptured abdominal aortic aneurysm surgical repair: 12-year single-center experience of Estonia, Scand. J. Surg., 107 (2018), 152-157. https://doi.org/10.1177/1457496917738923 doi: 10.1177/1457496917738923
[12]	R. E. Brightwell, A. M. Choong, A. G. Barnett, P. J. Walker, Changes in temperature affect the risk of abdominal aortic aneurysm rupture, ANZ J. Surg., 84 (2014), 871-876. https://doi.org/10.1111/ans.12446 doi: 10.1111/ans.12446
[13]	A. Saxena, E. Y. K. Ng, C. Manchanda, T. Canchi, Cardiac thermal pulse at the neck-skin surface as a measure of stenosis in the carotid artery, Therm. Sci. Eng. Prog., 19 (2020), 100603. https://doi.org/10.1016/j.tsep.2020.100603 doi: 10.1016/j.tsep.2020.100603
[14]	E. Y. K. Ng, E. Y. L. Pang, Thermal elevation on midriff skin surface as a potential diagnostic feature for abdominal aortic aneurysm using infrared thermography (IRT), Int. J. Therm. Sci., 172 (2022), 107305. https://doi.org/10.1016/j.ijthermalsci.2021.107305 doi: 10.1016/j.ijthermalsci.2021.107305
[15]	C. M. Scotti, A. D. Shkolnik, S. C. Muluk, E. A. Finol, Fluid-structure interaction in abdominal aortic aneurysms: Effects of asymmetry and wall thickness, Biomed. Eng. Online, 4 (2005), 64. https://doi.org/10.1186/1475-925x-4-64 doi: 10.1186/1475-925x-4-64
[16]	W. Hao, S. Gong, S. Wu, J. Xu, M. Go, A. Friedman, et al., A mathematical model of aortic aneurysm formation, PLoS One, 12 (2017), e0170807. https://doi.org/10.1371/journal.pone.0170807 doi: 10.1371/journal.pone.0170807
[17]	S. Lin, X. Han, Y. Bi, S. Ju, L. Gu, Fluid-structure interaction in abdominal aortic aneurysm: Effect of modeling techniques, Biomed Res. Int., 2017 (2017), 7023078. https://doi.org/10.1155/2017/7023078 doi: 10.1155/2017/7023078
[18]	G. A. Holzapfel, T. C. Gasser, R. W. Ogden, A new constitutive framework for arterial wall mechanics and a comparative study of material models, J. Elast., 61 (2000), 1-48. https://doi.org/10.1023/A:1010835316564 doi: 10.1023/A:1010835316564
[19]	J. Humphrey, Introduction, World Dev., 23 (1995), 1-7. https://doi.org/10.1016/0305-750X(95)90011-O doi: 10.1016/0305-750X(95)90011-O
[20]	G. A. Holzapfel, R. W. Ogden, Constitutive modelling of arteries, Proc. Math. Phys. Eng. Sci., 466 (2010), 1551-1597. https://doi.org/10.1098/rspa.2010.0058 doi: 10.1098/rspa.2010.0058
[21]	C. Reeps, M. Gee, A. Maier, M. Gurdan, H. H. Eckstein, W. A. Wall, The impact of model assumptions on results of computational mechanics in abdominal aortic aneurysm, J. Vasc. Surg., 51 (2010), 679-688. https://doi.org/10.1016/j.jvs.2009.10.048 doi: 10.1016/j.jvs.2009.10.048
[22]	C. M. Scotti, J. Jimenez, S. C. Muluk, E. A. Finol, Wall stress and flow dynamics in abdominal aortic aneurysms: finite element analysis vs. fluid-structure interaction, Comput. Methods Biomech. Biomed. Eng., 11 (2008), 301-322. https://doi.org/10.1080/10255840701827412 doi: 10.1080/10255840701827412
[23]	J. H. Leung, A. R. Wright, N. Cheshire, J. Crane, S. A. Thom, A. D. Hughes, et al., Fluid structure interaction of patient specific abdominal aortic aneurysms: A comparison with solid stress models, Biomed. Eng. Online, 5 (2006), 33. https://doi.org/10.1186/1475-925X-5-33 doi: 10.1186/1475-925X-5-33
[24]	Y. Mesri, H. Niazmand, A. Deyranlou, Numerical study on fluid-structure interaction in a patient-specific abdominal aortic aneurysm for evaluating wall heterogeneity and material model effects on its rupture, J. Appl. Fluid Mech., 10 (2017), 1699-1709. https://doi.org/10.18869/acadpub.jafm.73.243.27678 doi: 10.18869/acadpub.jafm.73.243.27678
[25]	A. Grytsan, P. N. Watton, G. A. Holzapfel, A thick-walled fluid-solid-growth model of abdominal aortic aneurysm evolution: application to a patient-specific geometry, J. Biomech. Eng., 137 (2015). https://doi.org/10.1115/1.4029279 doi: 10.1115/1.4029279
[26]	H. Schmid, A. Grytsan, E. Poshtan, P. N. Watton, M. Itskov, Influence of differing material properties in media and adventitia on arterial adaptation—Application to aneurysm formation and rupture, Comput. Methods Biomech. Biomed. Eng., 16 (2013), 33-53. https://doi.org/10.1080/10255842.2011.603309 doi: 10.1080/10255842.2011.603309
[27]	K. Y. Volokh, Modeling failure of soft anisotropic materials with application to arteries, J. Mech. Behav. Biomed. Mater., 4 (2011), 1582-1594. https://doi.org/10.1016/j.jmbbm.2011.01.002 doi: 10.1016/j.jmbbm.2011.01.002
[28]	P. N. Watton, N. A. Hill, M. Heil, A mathematical model for the growth of the abdominal aortic aneurysm, Biomech. Model. Mechanobiol., 3 (2004), 98-113. https://doi.org/10.1007/s10237-004-0052-9 doi: 10.1007/s10237-004-0052-9
[29]	M. L. Raghavan, D. A. Vorp, Toward a biomechanical tool to evaluate rupture potential of abdominal aortic aneurysm: identification of a finite strain constitutive model and evaluation of its applicability, J. Biomech., 33 (2000), 475-482. https://doi.org/10.1016/S0021-9290(99)00201-8 doi: 10.1016/S0021-9290(99)00201-8
[30]	M. J. Riedl, Optical Design Fundamentals for Infrared Systems, SPIE press, Bellingham, 2001.
[31]	A. Saxena, E. Y. K. Ng, M. Mathur, C. Manchanda, N. A. Jajal, Effect of carotid artery stenosis on neck skin tissue heat transfer, Int. J. Therm. Sci., 145 (2019), 106010. https://doi.org/10.1016/j.ijthermalsci.2019.106010 doi: 10.1016/j.ijthermalsci.2019.106010
[32]	G. Varjú, C. F. Pieper, J. B. Renner, V. B. Kraus, Assessment of hand osteoarthritis: Correlation between thermographic and radiographic methods, Rheumatology, 43 (2004), 915-919. https://doi.org/10.1093/rheumatology/keh204 doi: 10.1093/rheumatology/keh204
[33]	K. Woźniak, L. Szyszka-Sommerfeld, G. Trybek, D. Piątkowska, Assessment of the sensitivity, specificity, and accuracy of thermography in identifying patients with TMD, Med. Sci. Monit., 21 (2015), 1485-1493. https://doi.org/10.12659/MSM.893863 doi: 10.12659/MSM.893863
[34]	S. V. Patankar, Numerical Heat Transfer and Fluid Flow, CRC Press, Boca Raton, (1980), 214.
[35]	C. Childs, H. Soltani, Abdominal cutaneous thermography and perfusion mapping after caesarean section: A scoping review, Int. J. Environ. Res. Public Health, 17 (2020), 8693. https://doi.org/10.3390/ijerph17228693 doi: 10.3390/ijerph17228693
[36]	R. B. Barnes, Thermography of the human body, Science, 140 (1963), 870-877. https://doi.org/10.1126/science.140.3569.870 doi: 10.1126/science.140.3569.870
[37]	O. Ley, T. Kim, Determination of atherosclerotic plaque temperature in large arteries, Int. J. Therm. Sci., 47 (2008), 147-156. https://doi.org/10.1016/j.ijthermalsci.2007.01.034 doi: 10.1016/j.ijthermalsci.2007.01.034
[38]	T. Canchi, S. D. Kumar, E. Y. K. Ng, S. Narayanan, A review of computational methods to predict the risk of rupture of abdominal aortic aneurysms, Biomed. Res. Int., 2015 (2015), 861627. https://doi.org/10.1155/2015/861627 doi: 10.1155/2015/861627
[39]	J. Xu, A. Psikuta, J. Li, S. Annaheim, R. M. Rossi, Influence of human body geometry, posture and the surrounding environment on body heat loss based on a validated numerical model, Build. Environ., 166 (2019), 106340. https://doi.org/10.1016/j.buildenv.2019.106340 doi: 10.1016/j.buildenv.2019.106340
[40]	Y. Kurazumi, T. Tsuchikawa, N. Matsubara, T. Horikoshi, Convective heat transfer area of the human body, Eur. J. Appl. Physiol., 93 (2004), 273-285. https://doi.org/10.1007/s00421-004-1207-1 doi: 10.1007/s00421-004-1207-1
[41]	Y. Kurazumi, T. Tsuchikawa, J. Ishii, K. Fukagawa, Y. Yamato, N. Matsubara, Radiative and convective heat transfer coefficients of the human body in natural convection, Build. Environ., 43 (2008), 2142-2153. https://doi.org/10.1016/j.buildenv.2007.12.012 doi: 10.1016/j.buildenv.2007.12.012
[42]	K. Ouriel, R. M. Green, C. Donayre, C. K. Shortell, J. Elliott, J. A. DeWeese, An evaluation of new methods of expressing aortic aneurysm size: Relationship to rupture, J. Vasc. Surg., 15 (1992), 12-20. https://doi.org/10.1016/0741-5214(92)70008-9 doi: 10.1016/0741-5214(92)70008-9
[43]	C. M. Scotti, E. A. Finol, Compliant biomechanics of abdominal aortic aneurysms: A fluid-structure interaction study, Comput. Struct., 85 (2007), 1097-1113. https://doi.org/10.1016/j.compstruc.2006.08.041 doi: 10.1016/j.compstruc.2006.08.041
[44]	S. Bernad, E. Bernad, T. Barbat, C. Brisan, V. Albulescu, An analysis of blood flow dynamics in AAA, in Etiology, Pathogenesis and Pathophysiology of Aortic Aneurysms and Aneurysm Rupture, IntechOpen, 2011.
[45]	W. S. Cobb, J. M. Burns, K. W. Kercher, B. D. Matthews, H. J. Norton, B. T. Heniford, Normal intraabdominal pressure in healthy adults, J. Surg. Res., 129 (2005), 231-235. https://doi.org/10.1016/j.jss.2005.06.015 doi: 10.1016/j.jss.2005.06.015
[46]	B. Hernández-Gascón, E. Peña, H. Melero, G. Pascual, M. Doblaré, M. P. Ginebra, et al., Mechanical behaviour of synthetic surgical meshes: Finite element simulation of the herniated abdominal wall, Acta Biomater., 7 (2011), 3905-3913. https://doi.org/10.1016/j.actbio.2011.06.033 doi: 10.1016/j.actbio.2011.06.033
[47]	E. S. Di Martino, G. Guadagni, A. Fumero, G. Ballerini, R. Spirito, P. Biglioli, et al., Fluid-structure interaction within realistic three-dimensional models of the aneurysmatic aorta as a guidance to assess the risk of rupture of the aneurysm, Med. Eng. Phys., 23 (2001), 647-655. https://doi.org/10.1016/S1350-4533(01)00093-5 doi: 10.1016/S1350-4533(01)00093-5
[48]	M. H. Cardoso, Experimental study of the human anterolateral abdominal wall: Biomechanical properties of fascia and muscles, 2012. Available from: https://repositorio-aberto.up.pt/bitstream/10216/65576/1/000154315.pdf.
[49]	ANSYS Mechanical Material Reference, Release 2021 R1 ed., ANSYS, Inc., Canonsburg, (2021), 93-106.
[50]	J. Xiang, V. M. Tutino, K. V. Snyder, H. Meng, CFD: Computational fluid dynamics or confounding factor dissemination? The role of hemodynamics in intracranial aneurysm rupture risk assessment, Am. J. Neuroradiol., 35 (2014), 1849-1857. https://doi.org/10.3174/ajnr.A3710 doi: 10.3174/ajnr.A3710
[51]	F. Gijsen, F. van de Vosse, J. D. Janssen, The influence of the non-Newtonian properties of blood on the flow in large arteries: Steady flow in a carotid bifurcation model, J. Biomech., 32 (1999), 601-608. https://doi.org/10.1016/S0021-9290(99)00015-9 doi: 10.1016/S0021-9290(99)00015-9
[52]	K. Perktold, M. Resch, H. Florian, Pulsatile non-Newtonian flow characteristics in a three-dimensional human carotid bifurcation model, J. Biomech. Eng., 113 (1991), 464-475. https://doi.org/10.1115/1.2895428 doi: 10.1115/1.2895428
[53]	J. P. V. Geest, M. S. Sacks, D. A. Vorp, The effects of aneurysm on the biaxial mechanical behavior of human abdominal aorta, J. Biomech., 39 (2006), 1324-1334. https://doi.org/10.1016/j.jbiomech.2005.03.003 doi: 10.1016/j.jbiomech.2005.03.003
[54]	D. A. Vorp, M. L. Raghavan, M. W. Webster, Mechanical wall stress in abdominal aortic aneurysm: Influence of diameter and asymmetry, J. Vasc. Surg., 27 (1998), 632-639. https://doi.org/10.1016/S0741-5214(98)70227-7 doi: 10.1016/S0741-5214(98)70227-7
[55]	S. Pasta, A. Rinaudo, A. Luca, M. Pilato, C. Scardulla, T. G. Gleason, et al., Difference in hemodynamic and wall stress of ascending thoracic aortic aneurysms with bicuspid and tricuspid aortic valve, J. Biomech., 46 (2013), 1729-1738. https://doi.org/10.1016/j.jbiomech.2013.03.029 doi: 10.1016/j.jbiomech.2013.03.029
[56]	J. Pearce, S. Thomsen, Blood vessel architectural features and their effects on thermal phenomena in Matching the Energy Source to the Clinical Need: A Critical Review, SPIE, (2000), 122-168.
[57]	F. A. Duck, Chapter 2—Thermal properties of tissue, in Physical Properties of Tissues (ed. F. A. Duck), Academic Press, London, (1990), 9-42.
[58]	G. Giannakoulas, G. Giannoglou, J. Soulis, T. Farmakis, S. Papadopoulou, G. Parcharidis, et al., A computational model to predict aortic wall stresses in patients with systolic arterial hypertension, Med. Hypotheses, 65 (2005), 1191-1195. https://doi.org/10.1016/j.mehy.2005.06.017 doi: 10.1016/j.mehy.2005.06.017
[59]	G. J. Müller, A. Roggan, Laser-induced Interstitial Thermotherapy, SPIE Optical Engineering Press, 1995.
[60]	M. Halabian, B. Beigzadeh, A. Karimi, H. A. Shirazi, M. H. Shaali, A combination of experimental and finite element analyses of needle-tissue interaction to compute the stresses and deformations during injection at different angles, J. Clin. Monit. Comput., 30 (2016), 965-975. https://doi.org/10.1007/s10877-015-9801-9 doi: 10.1007/s10877-015-9801-9
[61]	S. Mukhopadhyay, M. S. Mandal, S. Mukhopadhyay, Heat transfer in pulsatile blood flow obeying Cross viscosity model through an artery with aneurysm, J. Eng. Math., 131 (2021), 6. https://doi.org/10.1007/s10665-021-10172-w doi: 10.1007/s10665-021-10172-w
[62]	Y. G. Stergiou, A. G. Kanaris, A. A. Mouza, S. V. Paras, Fluid-structure interaction in abdominal aortic aneurysms: Effect of haematocrit, Fluids, 4 (2019), 11. https://doi.org/10.3390/fluids4010011 doi: 10.3390/fluids4010011
[63]	J. D. Humphrey, G. A. Holzapfel, Mechanics, mechanobiology, and modeling of human abdominal aorta and aneurysms, J. Biomech., 45 (2012), 805-814. https://doi.org/10.1016/j.jbiomech.2011.11.021 doi: 10.1016/j.jbiomech.2011.11.021
[64]	Y. Chu, D. S. Bilal, M. Hajizadeh, Hybrid ferrofluid along with MWCNT for augmentation of thermal behavior of fluid during natural convection in a cavity, Math. Methods Appl. Sci., 2020. https://doi.org/10.1002/mma.6937 doi: 10.1002/mma.6937

Reader Comments

Your name:*

Email:*
© 2022 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)