
The numerical simulation of several virtual scenarios arising in cardiac mechanics poses a computational challenge that can be alleviated if traditional full-order models (FOMs) are replaced by reduced order models (ROMs). For example, in the case of problems involving a vector of input parameters related, e.g., to material coefficients, projection-based ROMs provide mathematically rigorous physics-driven surrogate ROMs. In this work we demonstrate how, once trained, ROMs yield extremely accurate predictions (according to a prescribed tolerance) – yet cheaper than the ones provided by FOMs – of the structural deformation of the left ventricular tissue over an entire heartbeat, and of related output quantities of interest, such as the pressure-volume loop, for any desired input parameter values within a prescribed parameter range. However, the construction of ROM approximations for time-dependent cardiac mechanics is not straightforward, because of the highly nonlinear and multiscale nature of the problem, and almost never addressed. Our approach relies on the reduced basis method for parameterized partial differential equations. This technique performs a Galerkin projection onto a low-dimensional space for the displacement variable; the reduced space is built from a set of solution snapshots – obtained for different input parameter values and time instances – of the high-fidelity FOM, through the proper orthogonal decomposition technique. Then, suitable hyper-reduction techniques, such as the Discrete Empirical Interpolation Method, are exploited to efficiently handle nonlinear and parameter-dependent terms. In this work we show how a fast and reliable approximation of the time-dependent cardiac mechanical model can be achieved by a projection-based ROM, taking into account both passive and active mechanics for the left ventricle providing all the building blocks of the methodology, and highlighting those challenging aspects that are still open.
Citation: Ludovica Cicci, Stefania Fresca, Stefano Pagani, Andrea Manzoni, Alfio Quarteroni. Projection-based reduced order models for parameterized nonlinear time-dependent problems arising in cardiac mechanics[J]. Mathematics in Engineering, 2023, 5(2): 1-38. doi: 10.3934/mine.2023026
[1] | Jianxin Yang, Zhipeng Qiu, Xue-Zhi Li . Global stability of an age-structured cholera model. Mathematical Biosciences and Engineering, 2014, 11(3): 641-665. doi: 10.3934/mbe.2014.11.641 |
[2] | Mingtao Li, Xin Pei, Juan Zhang, Li Li . Asymptotic analysis of endemic equilibrium to a brucellosis model. Mathematical Biosciences and Engineering, 2019, 16(5): 5836-5850. doi: 10.3934/mbe.2019291 |
[3] | Linhua Zhou, Meng Fan, Qiang Hou, Zhen Jin, Xiangdong Sun . Transmission dynamics and optimal control of brucellosis in Inner Mongolia of China. Mathematical Biosciences and Engineering, 2018, 15(2): 543-567. doi: 10.3934/mbe.2018025 |
[4] | Zongmin Yue, Yuanhua Mu, Kekui Yu . Dynamic analysis of sheep Brucellosis model with environmental infection pathways. Mathematical Biosciences and Engineering, 2023, 20(7): 11688-11712. doi: 10.3934/mbe.2023520 |
[5] | Mingtao Li, Guiquan Sun, Juan Zhang, Zhen Jin, Xiangdong Sun, Youming Wang, Baoxu Huang, Yaohui Zheng . Transmission dynamics and control for a brucellosis model in Hinggan League of Inner Mongolia, China. Mathematical Biosciences and Engineering, 2014, 11(5): 1115-1137. doi: 10.3934/mbe.2014.11.1115 |
[6] | Zhisheng Shuai, P. van den Driessche . Impact of heterogeneity on the dynamics of an SEIR epidemic model. Mathematical Biosciences and Engineering, 2012, 9(2): 393-411. doi: 10.3934/mbe.2012.9.393 |
[7] | Yaoyao Qin, Xin Pei, Mingtao Li, Yuzhen Chai . Transmission dynamics of brucellosis with patch model: Shanxi and Hebei Provinces as cases. Mathematical Biosciences and Engineering, 2022, 19(6): 6396-6414. doi: 10.3934/mbe.2022300 |
[8] | Qiuyi Su, Jianhong Wu . Impact of variability of reproductive ageing and rate on childhood infectious disease prevention and control: insights from stage-structured population models. Mathematical Biosciences and Engineering, 2020, 17(6): 7671-7691. doi: 10.3934/mbe.2020390 |
[9] | Pengyan Liu, Hong-Xu Li . Global behavior of a multi-group SEIR epidemic model with age structure and spatial diffusion. Mathematical Biosciences and Engineering, 2020, 17(6): 7248-7273. doi: 10.3934/mbe.2020372 |
[10] | Tsanou Berge, Samuel Bowong, Jean Lubuma, Martin Luther Mann Manyombe . Modeling Ebola Virus Disease transmissions with reservoir in a complex virus life ecology. Mathematical Biosciences and Engineering, 2018, 15(1): 21-56. doi: 10.3934/mbe.2018002 |
The numerical simulation of several virtual scenarios arising in cardiac mechanics poses a computational challenge that can be alleviated if traditional full-order models (FOMs) are replaced by reduced order models (ROMs). For example, in the case of problems involving a vector of input parameters related, e.g., to material coefficients, projection-based ROMs provide mathematically rigorous physics-driven surrogate ROMs. In this work we demonstrate how, once trained, ROMs yield extremely accurate predictions (according to a prescribed tolerance) – yet cheaper than the ones provided by FOMs – of the structural deformation of the left ventricular tissue over an entire heartbeat, and of related output quantities of interest, such as the pressure-volume loop, for any desired input parameter values within a prescribed parameter range. However, the construction of ROM approximations for time-dependent cardiac mechanics is not straightforward, because of the highly nonlinear and multiscale nature of the problem, and almost never addressed. Our approach relies on the reduced basis method for parameterized partial differential equations. This technique performs a Galerkin projection onto a low-dimensional space for the displacement variable; the reduced space is built from a set of solution snapshots – obtained for different input parameter values and time instances – of the high-fidelity FOM, through the proper orthogonal decomposition technique. Then, suitable hyper-reduction techniques, such as the Discrete Empirical Interpolation Method, are exploited to efficiently handle nonlinear and parameter-dependent terms. In this work we show how a fast and reliable approximation of the time-dependent cardiac mechanical model can be achieved by a projection-based ROM, taking into account both passive and active mechanics for the left ventricle providing all the building blocks of the methodology, and highlighting those challenging aspects that are still open.
Brucellosis, a zoonotic disease, is a natural epidemic disease that is not only prevalent among livestock and humans but also widely spread in wild animals [1,2]. There are human and animal brucellosis in most countries of the world, and high-risk areas of epidemics are mainly distributed in developing countries, such as Syria, Jordan, Zambia, Mongolia [3,4,5,6]. It has brought huge losses to the livestock industry worldwide and serious health problems to livestock-related practitioners, there are hundreds of thousands of new cases reported annually [6,7,8]. Infected domestic and wild animals and their excreta are the main source of infection, it is contact and pathogen infection which are the main modes of brucellosis transmission [9]. Human brucellosis is rarely transmitted to susceptible animals, and there is no infection between people reported [10,11]. Therefore, it is the eradication of animal brucellosis that is the only way to solving human health problem, and understanding the mechanisms and risk factors of the spread of brucellosis is one of the first problems that must be solved.
Brucellosis is transmitted to susceptible individuals mainly by contact with infectious individuals or by sucking pathogens from the environment [8]. In animals, the transmission of brucellosis mainly occurs between sexually mature animals, young animals may also be infected, but generally do not have any clinical manifestations and serological tests are usually non-positive, and the infection between young animals is so little as to be almost invisible [12]. In other words, sexually mature animals are very susceptible to brucellosis, latent infections may be found in young animals, but they are generally resistant [9]. It is important that many infected animals have a longer incubation period and that these animals remain serologically non-positive during this period. That is to say, latent animals may not be infectious and it is almost impossible to be found through detection, which is an important risk factor in the elimination of brucellosis [9]. Therefore, the impact of these mechanisms on brucellosis is worthy to be studied using mathematical models.
Statistical methods such as descriptive statistics, correlation analysis and time series analysis have been widely used for quantitative assessment of risk management measures for brucellosis, among which is very worthy of concern about the study of bison, elk and livestock brucellosis (see [13,14,15,16,17]). Theoretical studies on the impacts of the transmission mechanisms of brucellosis and applied studies on the assessments of risk management measures have also been studied using the kinetic model (see [18,19,20,21,22]). Especially in recent years, many kinetic models with indirect transmission have been established to analyze brucellosis transmission (see [23,24,25,26,27]). Although there have been many studies on brucellosis, there are still many transmission mechanisms and risk factors that are not considered in existing models. For example, infected individuals may have no infectivity in the early stage, during which time is different; and animals that are not sexually mature are hardly infected with direct and indirect modes [9]. Therefore, in the present work, a multi-stage dynamic model with distributed time delay is proposed involving the above risk factors and general nonlinear incidences. The existence and uniqueness of the endemic equilibrium is analyzed, and the local and global asymptotic stability of equilibria is proved.
The rest of this manuscript is constructed as follows. The dynamic model with distributed time delay and some preliminary results are given in Section 2. In Section 3, the content of research is the global dynamics of the disease-free equilibrium. The global stability of the equilibrium of persistent infection are analyzed in Section 4. In Section 5, the stability results are further explained by numerical simulation. A summary and further discussion is proposed in Section 6.
According to the transmission mechanisms of animal brucellosis, we classify animal population into three compartments: The young susceptible compartment S1(t), the adult susceptible compartment S2(t) (sexually mature), the infected compartment I(t). Similar to the definition in literature [25], B(t) can be defined as the concentration or number of pathogens in the environment. Some explanations and assumptions about the kinetic model are listed as follows. (Ⅰ) Since brucellosis can cause abortion and significantly reduce the survival rate of young animals, we suppose that infected animals have no birth rate. (Ⅱ) The supplementary rate of young susceptible animal population S1(t) is mainly derived from birth and import, it is assumed that the supplementary rate is A+bS2(t). (Ⅲ) Usually, there are two different stages for infected animals, which are no infection force in the early stage and are infectious carriers when infected animals begin to shed brucella. Therefore, we make use of a delay τ to describe the time from a susceptible individual to a infectious individual, and since the length of this time is varied, thus the delay τ is a distributed parameter in the interval [0,h], where h is the maximum value of the delay. As a result, the susceptible class S2(t) is reduced due to infection at the rate S2g(B)+S2∫h0ρ(τ)f(I(t−τ))e−(μ2+c)τdτ. (Ⅳ) Since the infected animal sheds the pathogen after the time τ, then the pathogen shedding rate of the infected animal is ∫h0ρ(τ)h(I(t−τ))e−(μ2+c)τdτ. (Ⅴ) In animal breeding, only the basic ewes breed newborns, and adult animals are widely used for slaughter or trade, then the birth rate b of adult animals is assumed to be less than the elimination rate μ2 of adult animals. Therefore, the modelling for animal brucellosis is given through the following distributed time-delay system:
{dS1dt=A+bS2−μ1S1−σ1S1,dS2dt=σ1S1−S2g(B)−S2∫h0ρ(τ)f(I(t−τ))e−(μ2+c)τdτ−μ2S2,dIdt=S2∫h0ρ(τ)f(I(t−τ))e−(μ2+c)τdτ+S2g(B)−(μ2+c)I,dBdt=∫h0ρ(τ)h(I(t−τ))e−(μ2+c)τdτ−dB. | (2.1) |
Here, b is the birth rate of adult animals. μ1 and μ2 are the elimination rates of young and adult animals, respectively. σ1 is the transfer rate from young susceptible individuals to sexually mature individuals. The culling rate is defined by c and d is the decay rate of the pathogen. f(I) and g(B) are contact and indirect infection rates, respectively. ρ(τ) is a distributed function, it is non-negative, continuous and ∫h0ρ(τ)dτ=1.
The initial conditions for system (2.1) are given as follows:
{S1(x)=ϕ1(x),S2(x)=ϕ2(x),I(x)=ϕ3(x),B(x)=ϕ4(x),x∈[−h,0],h>0,ϕ=(ϕ1,ϕ2,ϕ3,ϕ4)∈C+⊂C. | (2.2) |
Here, C denotes the Banach space C([−h,0],ℜ4) of continuous functions mapping the interval [−h,0] into ℜ4 with the sup-norm ‖ϕ‖=supx∈[−h,0]|ϕi(x)|,i=1,2,3,4 for ϕ∈C. The nonnegative cone of C is defined as C+=C([−h,0],ℜ4+).
In order to make epidemiological significance for system (2.1), f,g and h are assumed to be second-order continuous differentiable functions and satisfy the following hypotheses:
(H1) f(0)=g(0)=h(0)=0 and f(I),g(B),h(I)>0 for I,B>0;
(H2) f′(I),g′(B)>0 and f″(I),g″(B)≤0 for B,I≥0;
(H3) h′(I)>0 and h″(I)≤0 for I≥0.
The function f may be βIp or saturation incidences kln(1+λIk) and βIp1+kIp with constants β,p,λ,k>0 [28,29]. The pathogen infection rate g may be λB1+TB or η(1−e−αB) with constants λ,η,α>0 and T≥0 [30]. The function h may be kI with constants k>0 [31].
It is easy to verify that the system (2.1) always has a disease-free equilibrium E0=(S01,S02,0,0), where
S01=Aμ2μ2(μ1+σ1)−bσ1,S02=Aσ1μ2(μ1+σ1)−bσ1. |
The reproduction number of sytem (2.1) is given by the following expression:
R0=nS02fI(0)μ2+c+nS02gB(0)hI(0)d(μ2+c)=R01+R02. | (2.3) |
where n=∫h0ρ(τ)e−(μ2+c)τdτ. e−(μ2+c)τ denotes the probability of survival from a newly infected susceptible individual to an infectious individual. n is the total survival rate of infected individuals. fI(0) is the infection rate by an infectious individual. gB(0) is the infection rate of brucella. 1μ2+c is the average life span of infected individual. hI(0)d represents the total number of brucella shed by infectious individuals. Therefore, based on the explanation in the literature [26], R0 relies on direct and indirect infection and can be divided into two parts R01 and R02.
In the next section, we establish that the solution of system (2.1) is non-negative and bounded and analyze the uniqueness of the equilibrium point of persistent infection. The following qualitative results indicate the nonnegativity and boundedness of the solution of system (2.1) with initial conditions (2.2).
Theorem 2.1. (S1(t),S2(t),I(t),B(t)) is the solution of system (2.1) with the initial conditions (2.2), then S1(t),S2(t),I(t),B(t) are nonnegative and ultimately bounded.
Proof. Since (˙S1(t)+˙S2(t))=A>0 for t∈[0,+∞) when S1(t)=S2(t)=0, it implies that S1(t)+S2(t)≥0 for all t∈[0,+∞). Therefore, if S1(t)=0, then S2(t)≥0 and ˙S1(t)=A+bS2>0, so that S1(t)≥0. We denote
a1(t)=∫h0ρ(τ)f(I(t−τ))e−(μ2+c)τdτ+g(B)+μ2,a2(t)=S2(t)(∫h0ρ(τ)f(I(t−τ))e−(μ2+c)τdτ+g(B)), |
and
a3(t)=∫h0ρ(τ)h(I(t−τ))e−(μ2+c)τdτ. |
From the last three equations in (2.1), it can conclude that
S2(t)=S2(0)e−∫t0a1(η)dη+σ1∫t0e−∫tηa1(ξ)dξS1(η)dη≥0,I(t)=I(0)e−(μ2+c)t+∫t0e−(μ2+c)(t−η)a2(η)dη≥0, |
and
B(t)=B(0)e−dt+∫t0e−d(t−η)a3(η)dη≥0 |
for all t≥0. Thus, S1(t),S2(t),I(t),B(t)≥0 for t≥0.
We now analyze the boundedness of the solution of system (2.1). In fact,
(˙S1(t)+˙S2(t)+˙I(t))≤A−μ1S1−(μ2−b)S2−(μ2+c)I |
≤A−μ(S1(t)+S2(t)+I(t)), |
where μ=min{μ1,μ2−b}.
Hence, lim supt→∞(S1(t)+S2(t)+I(t))≤Aμ, it implies that lim supt→∞I(t)≤Aμ and lim supt→∞B(t)≤ndh(Aμ). Therefore, S1(t),S2(t),I(t),B(t) are ultimately bounded.
Since the solution of system (2.1) is nonnegative and ultimately bounded, the set
Ω={(S1(⋅),S2(⋅),I(⋅),B(⋅))∈C+:‖S1+S2+I‖≤Aμ,‖B‖≤ndh(Aμ)}, |
is positively invariant for system (2.1).
In order to analyze the uniqueness of the positive solution of system (2.1), the following lemma is given:
Lemma 2.1. Assume that conditions (H1)−(H3) are satisfied, the functions f(I)I, h(I)I, g(B)B and g(h(I))I are monotonic decreasing for I,B>0.
Proof. Since f″(I)≤0, it shows that f′(I) is monotonic decreasing, it follows that
f(I)I=f(I)−f(0)I−0=f′(ξ1)≥f′(I),ξ1∈(0,I), |
and
(f(I)I)′=f′(I)I−f(I)I2≤0. |
That is to say, f(I)I is a monotonically decreasing function. On the basis of the above method, it can show that h(I)I and g(B)B are also monotonic decreasing. Noting that
g(h(I))I=g(h(I))−g(h(0))h(I)−h(0)h(I)I=g′(h(ξ2))h′(ξ3)≥g′(h(I))h′(I),ξ2,ξ3∈(0,I), |
It can deduce that
(g(h(I))I)′=g′(h(I))h′(I)I−g(h(I))I2≤0. |
Therefore, g(h(I))I is also monotonic decreasing.
For system (2.1), the endemic equilibrium E∗=(S∗1,S∗2,I∗,B∗) can be derived from the following algebraic equations:
{A+bS2=μ1S1+σ1S1,σ1S1=S2(nf(I)+g(B))+μ2S2,S2(nf(I)+g(B))=(μ2+c)I,nh(I)=dB. | (2.4) |
By direct calculation, it can be written as
{S1=A+bS2μ1+σ1,m(S2)=S2(nf(I)+g(B)),S2(nf(I)+g(B))=(μ2+c)I,B=H(I). |
where m(S2)=Aσ1−(μ2(μ1+σ1)−bσ1)S2μ1+σ1 and H(I)=ndh(I).
Let us define
F1(S2,I)≜m(S2)−S2(nf(I)+g(H(I))),F2(S2,I)≜S2(nf(I)+g(H(I)))−(μ2+c)I. |
Similar to the method of analysis in literature [26], using Lemma 2.1, the following result can be summarized:
Theorem 2.2. Assume that conditions (H1)−(H3) hold. Then there is a unique positive solution E∗=(S∗1,S∗2,I∗,B∗) of system (2.1) if and only if R0>1.
In this section, we show that the global stability of the equilibrium point E0 of system (2.1) is independent of the initial value. The following conclusions are first obtained.
Lemma 3.1. Assume that conditions (H1)−(H3) hold. The disease-free equilibrium E0 is locally asymptotically stable if R0≤1 and is unstable if R0>1.
Proof. The characteristic equation at E0 is
|λ+μ1+σ1−b00−σ1λ+μ2S02fI(0)Γ(λ)S02gB(0)00λ+μ2+c−Γ(λ)S02fI(0)−S02gB(0)00−hI(0)Γ(λ)λ+d|=0, | (3.1) |
where
Γ(λ)=∫h0ρ(τ)e−(μ2+c+λ)τdτ. |
It follows from (3.1) that
(λ+μ1+σ1)(λ+μ2)((λ+d)(λ+μ2+c−Γ(λ)S02fI(0))−Γ(λ)S02gB(0)hI(0))=bσ1((λ+d)(λ+μ2+c−Γ(λ)S02fI(0))−Γ(λ)S02gB(0)hI(0)). |
That is,
(λ+μ1+σ1)(λ+μ2)(λ+d)(λ+μ2+c)−bσ1(λ+d)(λ+μ2+c) |
=((λ+μ1+σ1)(λ+μ2)−bσ1)((λ+d)Γ(λ)S02fI(0)+Γ(λ)S02gB(0)hI(0)). |
It follows that
H1(λ)((λ+d)(λ+μ2+c)−(λ+d)Γ(λ)S02fI(0)−Γ(λ)S02gB(0)hI(0))=0, |
where H1(λ)=(λ+μ1+σ1)(λ+μ2)−bσ1. Since μ1+σ1+μ2>0 and μ2(μ1+σ1)−bσ1>0, then H1(λ) consists of two roots which are negative real parts. Therefore, we only analyze the distribution of the roots of the following equation:
H2(λ)=(λ+d)(λ+μ2+c)−(λ+d)Γ(λ)S02fI(0)−Γ(λ)S02gB(0)hI(0)=0. | (3.2) |
Assume now that R0>1, then
H2(0)=d(μ2+c)−dnS02fI(0)−nS02gB(0)hI(0))<0,H2(+∞)=+∞. |
Hence H2(λ) has at least one positive root in [0,+∞), then E0 is unstable if R0>1.
From (3.2), we have
(λ+d)(λ+μ2+c)=Γ(λ)(λS02fI(0)+dS02fI(0)+S02gB(0)hI(0))=(μ2+c)F(λ)n(λR01+dR0), |
or
(λ+d)(λμ2+c+1)=R0F(λ)n(λR01R0+d). | (3.3) |
Next, considering the case R0≤1. If λ=x+yi is a solution of (3.3), one shows that x<0. Otherwise, x≥0 implies
|λ+d|>|λR01R0+d|,|λμ2+c+1|>1,|R0F(λ)n|≤1, |
and thus
|(λ+d)(λμ2+c+1)|>|R0F(λ)n(λR01R0+d)|, |
this is a contradiction to (3.3). Therefore, all roots of Eq (3.3) have no zero and positive real parts when R0≤1, this shows that E0 is locally asymptotically stable.
Theorem 3.1. Assume that conditions (H1)−(H3) are established. If R0≤1, the disease-free equilibrium E0=(S01,S02,0,0) of system (2.1) is globally asymptotically stable.
Proof. Since the functions f(I)I, g(B)B and h(I)I are decreasing, then we have
nS2f(I)(μ2+c)I≤limI→0+nS02f(I)(μ2+c)I=nS02fI(0)μ2+c≜b1,S2g(B)dB≤limB→0+S02g(B)dB=S02gB(0)d≜b2,nh(I)(μ2+c)I≤limI→0+nh(I)(μ2+c)I=nhI(0)μ2+c. |
Define
J=(10−nhI(0)μ2+c1),(a1,a2)=(b1,b2)J−1. |
We find a1=R0 and define a Lyapunov functional L as follows:
L(t)=L1(t)+L2(t)+L3(t), |
where
L1=R0(S1−S01−S01lnS1S01+S2−S02−S02lnS2S02+I)+a2B,L2=R0S02∫h0Ψ(τ)f(I(t−τ))dτ,L3=a2∫h0Ψ(τ)h(I(t−τ))dτ, |
and
Ψ(τ)=∫hτφ(s)ds,φ(s)=ρ(s)e−(μ2+c)s. |
Then the derivative of L1 along the positive solutions of system (2.1) is
dL1dt=R0(1−S01S1)dS1dt+R0(1−S02S2)dS2dt+R0dIdt+a2dBdt=R0(2A+bS02+σ1S01−μ1S1−(μ2−b)S2−S01AS1−bS01S2S1−σ1S02S1S2)+R0(S02g(B)+S02∫h0φ(τ)f(I(t−τ))dτ−(μ2+c)I)+a2(∫h0φ(τ)h(I(t−τ))dτ−dB). | (3.4) |
Calculating the derivative of L2(t) along the solutions of system (2.1), one obtains
dL2dt=R0S02∫h0Ψ(τ)df(I(t−τ))dtdτ=−R0S02∫h0Ψ(τ)df(I(t−τ))dτdτ=−R0S02(Ψ(τ)f(I(t−τ)))|h0+R0S02∫h0dΨ(τ)dτf(I(t−τ))dτ=R0S02nf(I(t))−R0S02∫h0φ(τ)f(I(t−τ))dτ. | (3.5) |
Similar to the above-used method, it can obtain that
dL3dt=a2nh(I(t))−a2∫h0φ(τ)h(I(t−τ))dτ. | (3.6) |
Combining the Eqs (3.4), (3.5) and (3.6), it follows that
dLdt=dL1dt+dL2dt+dL3dt=R0μ1S01(2−S1S01−S01S1)+R0bS02(2−S1S02S01S2−S01S2S1S02)+R0(μ2−b)S02(3−S2S02−S1S02S01S2−S01S1)+R0(nS02f(I)(μ2+c)I,S02g(B)dB)((μ2+c)I,dB)T−(R0,a2)(10−nh(I)(μ2+c)I1)((μ2+c)I,dB)T≤R0(nS02fI(0)(μ2+c),S02gB(0)d)((μ2+c)I,dB)T−(R0,a2)(10−nhI(0)μ2+c1)((μ2+c)I,dB)T=(R0−1)(b1,b2)((μ2+c)I,dB)T≤0. |
Similar to the analysis of Theorem 1 in the literature [32], it follows from Lemma 3.1 that the disease-free equilibrium E0 is globally asymptotically stable by LaSalle's Invariance Principle [33].
According to Lemma 3.1, the disease-free steady state solution E0 is unstable when R0>1. Using Theorem 4.2 in [34], the uniform permanence of system (2.1) can be proven, the process is ignored here. In this following section, by constructing a Lyapunov functional, the global stability of the equilibrium of persistent infection is proved. We first analyze its local stability.
Lemma 4.1. Assume that conditions (H1)−(H3) are ture, If R0>1, the endemic equilibrium E∗ of system (2.1) exists and is locally asymptotically stable.
Proof. For system (2.1), the characteristic equation at E∗ is
|λ+μ1+σ1−b00−σ1H3(λ)Γ(λ)S∗2fI(I∗)S∗2gB(B∗)0−nf(I∗)−g(B∗)H4(λ)−S∗2gB(B∗)00−Γ(λ)hI(I∗)λ+d|=0, |
where H3(λ)=λ+nf(I∗)+g(B∗)+μ2, H4(λ)=λ+μ2+c−Γ(λ)S∗2fI(I∗). By simple calculation, one obtains
(λ+d)(λ+μ2+c)Φ1(λ)=Φ2(λ)Γ(λ)((λ+d)S∗2fI(I∗)+hI(I∗)S∗2gB(B∗))=Φ2(λ)Γ(λ)nΛn(λS∗2fI(I∗)Λ+d), |
or
(λ+d)(λμ2+c+1)Φ1(λ)=Φ2(λ)Γ(λ)nΛnμ2+c(λS∗2fI(I∗)Λ+d), | (4.1) |
where
Φ1(λ)=(λ+μ1+σ1)(λ+μ2+nf(I∗)+g(B∗))−bσ1,Φ2(λ)=(λ+μ1+σ1)(λ+μ2)−bσ1,Λ=S∗2fI(I∗)+1dS∗2hI(I∗)gB(B∗) |
Noting that from Lemma 1
Λnμ2+c=S∗2fI(I∗)n+1dS∗2hI(I∗)gB(B∗)nμ2+c≤S∗2(f(I∗)nI∗+ndh(I∗)I∗g(B∗)B∗)μ2+c=S∗2(f(I∗)n+g(B∗))I∗(μ2+c)=1. |
Assume λ=x+yi is a solution of (4.1). If x≥0, we then have
|λ+d|>|λS∗2fI(I∗)Λ+d|,|λμ2+c+1|>1,|F(λ)n|≤1, |
and
|Φ1(λ)Φ2(λ)|=|(λ+μ1+σ1)(λ+μ2+nf(I∗)+g(B∗))−bσ1(λ+μ1+σ1)(λ+μ2)−bσ1|=|λ+μ2+nf(I∗)+g(B∗)−bσ1λ+μ1+σ1λ+μ2−bσ1λ+μ1+σ1|=|M1+nf(I∗)+g(B∗)+M2iM1+M2i|>1, |
where M1=x+μ2−bσ1(x+μ1+σ1)(x+μ1+σ1)2+y2>0, M2=y+bσ1y(x+μ1+σ1)2+y2. So it concludes that
|(λ+d)(λμ2+c+1)Φ1(λ)|>|Φ2(λ)F(λ)nΛnμ2+c(λS∗2fI(I∗)Λ+d)|, |
which is a contradiction to (4.1). Therefore, the Eq (4.1) can not have any roots with a nonnegative real part, this implies that E∗ is locally asymptotically stable if R0>1.
Theorem 4.1. Assume that conditions (H1)−(H3) hold. If R0>1, the endemic equilibrium E∗=(S∗1,S∗2,I∗,B∗) of system (2.1) is globally asymptotically stable.
Proof. Define
L1=S1−S∗1−S∗1lnS1S∗1+S2−S∗2−S∗2lnS2S∗2+I−I∗−I∗lnII∗+S∗2g(B∗)nh(I∗)(B−B∗−B∗lnBB∗). |
Finding the time derivative of L1 along the positive solutions of system (2.1) gives
dL1dt=(1−S∗1S1)dS1dt+(1−S∗2S2)dS2dt+(1−I∗I)dIdt+S∗2g(B∗)nh(I∗)(1−B∗B)dBdt=2A+bS∗2+σ1S∗1−μ1S1−(μ2−b)S2−S∗1AS1−bS∗1S2S1−σ1S∗2S1S2+S∗2g(B)+S∗2∫h0φ(τ)f(I(t−τ))dτ−(μ2+c)I−I∗IS2g(B)−I∗IS2∫h0φ(τ)f(I(t−τ))dτ+S∗2g(B∗)nh(I∗)(∫h0φ(τ)h(I(t−τ))dτ−dB+dB∗)−S∗2g(B∗)nh(I∗)B∗B∫h0φ(τ)h(I(t−τ))dτ. |
By the Eq (2.4), we get
dL1dt=μ1S∗1(2−S1S∗1−S∗1S1)+(μ2−b)S∗2(3−S2S∗2−S1S∗2S∗1S2−S∗1S1)+bS∗2(2−S1S∗2S∗1S2−S∗1S2S1S∗2)+nS∗2f(I∗)(3−S∗1S1−S1S∗2S2S∗1−II∗)+S∗2g(B∗)(3+g(B)g(B∗)−S∗1S1−S1S∗2S2S∗1−II∗−S2g(B)I∗S∗2g(B∗)I)+S∗2∫h0φ(τ)f(I(t−τ))dτ−I∗IS2∫h0φ(τ)f(I(t−τ))dτ+S∗2g(B∗)(h(I)h(I∗)−BB∗−B∗h(I)Bh(I∗)+1)+S∗2g(B∗)nh(I∗)(∫h0φ(τ)h(I(t−τ))dτ−nh(I)+B∗Bnh(I))−S∗2g(B∗)nh(I∗)B∗B∫h0φ(τ)h(I(t−τ))dτ. | (4.2) |
Since the function ν(x)=1−x+lnx is nonpositive for x>0 and ν(x)=0 if and only if x=1, so we define
L2=−S∗2f(I∗)∫h0Ψ(τ)ν(f(I(t−τ))f(I∗))dτ. |
A direct calculation shows that
dL2dt=−S∗2f(I∗)∫h0Ψ(τ)ddtν(f(I(t−τ))f(I∗))dτ=S∗2f(I∗)∫h0Ψ(τ)ddτν(f(I(t−τ))f(I∗))dτ=S∗2f(I∗)(Ψ(τ)ν(f(I(t−τ))f(I∗))∣h0+∫h0φ(τ)ν(f(I(t−τ))f(I∗))dτ)=−S∗2f(I∗)∫h0φ(τ)ν(f(I(t))f(I∗))dτ+S∗2f(I∗)∫h0φ(τ)ν(f(I(t−τ))f(I∗))dτ=S∗2f(I∗)∫h0φ(τ)(f(I)f(I∗)−f(I(t−τ))f(I∗)+lnf(I(t−τ))f(I))dτ=nS∗2f(I∗)f(I)f(I∗)+S∗2f(I∗)∫h0φ(τ)(−f(I(t−τ))f(I∗)+lnf(I(t−τ))f(I))dτ. | (4.3) |
Define
L3=−S∗2g(B∗)n∫h0Ψ(τ)ν(h(I(t−τ))h(I∗))dτ. |
Calculating the time derivative of L3(t), one obtains
dL3dt=−S∗2g(B∗)n∫h0Ψ(τ)ddtν(h(I(t−τ))h(I∗))dτ=S∗2g(B∗)n∫h0Ψ(τ)ddτν(h(I(t−τ))h(I∗))dτ=S∗2g(B∗)n∫h0φ(τ)(h(I)h(I∗)−h(I(t−τ))h(I∗)+lnh(I(t−τ))h(I))dτ=S∗2g(B∗)h(I)h(I∗)+S∗2g(B∗)n∫h0φ(τ)(−h(I(t−τ))h(I∗)+lnh(I(t−τ))h(I))dτ. | (4.4) |
For system (2.1), the following Lyapunov functional is considered:
L=L1+L2+L3. |
From (4.2), (4.3) and (4.4), we can get
dLdt=μ1S∗1(2−S1S∗1−S∗1S1)+(μ2−b)S∗2(3−S2S∗2−S1S∗2S∗1S2−S∗1S1)+bS∗2(2−S1S∗2S∗1S2−S∗1S2S1S∗2)+nS∗2f(I∗)(3+f(I)f(I∗)−S∗1S1−S1S∗2S2S∗1−II∗−S2f(I)I∗S∗2f(I∗)I)+S∗2g(B∗)(3+g(B)g(B∗)−S∗1S1−S1S∗2S2S∗1−II∗−S2g(B)I∗S∗2g(B∗)I)+S∗2g(B∗)(h(I)h(I∗)−BB∗−B∗h(I)Bh(I∗)+1)+S∗2f(I∗)∫h0φ(τ)F(τ)dτ+S∗2g(B∗)n∫h0φ(τ)H(τ)dτ, | (4.5) |
where
∫h0φ(τ)F(τ)dτ=∫h0φ(τ)(S2f(I)I∗S∗2f(I∗)I−S2I∗f(I(t−τ))S∗2f(I∗)I+lnf(I(t−τ))f(I))dτ=∫h0φ(τ)(ν(S2I∗f(I(t−τ))S∗2f(I∗)I)−ν(S2f(I)I∗S∗2f(I∗)I))dτ=∫h0φ(τ)ν(S2I∗f(I(t−τ))S∗2f(I∗)I)dτ−nν(S2f(I)I∗S∗2f(I∗)I), | (4.6) |
and
∫h0φ(τ)H(τ)dτ=∫h0φ(τ)(h(I)B∗h(I∗)B−B∗h(I(t−τ))h(I∗)B+lnh(I(t−τ))h(I))dτ=∫h0φ(τ)ν(B∗h(I(t−τ))h(I∗)B)dτ−nν(h(I)B∗h(I∗)B). | (4.7) |
it follows from (4.5), (4.6) and (4.7) that
dLdt=μ1S∗1(2−S1S∗1−S∗1S1)+(μ2−b)S∗2(3−S2S∗2−S1S∗2S∗1S2−S∗1S1)+bS∗2(2−S1S∗2S∗1S2−S∗1S2S1S∗2)+nS∗2f(I∗)(f(I)f(I∗)−1)(1−f(I∗)If(I)I∗)+nS∗2f(I∗)(ν(S∗1S1)+ν(S1S∗2S∗1S2)+ν(S2f(I)I∗S∗2f(I∗)I)+ν(f(I∗)If(I)I∗))+S∗2g(B∗)((g(B)g(B∗)−1)(1−g(B∗)Bg(B)B∗)+(h(I)h(I∗)−1)(1−h(I∗)Ih(I)I∗))+S∗2g(B∗)(ν(S∗1S1)+ν(S1S∗2S∗1S2)+ν(S2g(B)I∗S∗2g(B∗)I)+ν(g(B∗)Bg(B)B∗))+S∗2g(B∗)(ν(h(I)B∗h(I∗)B)+ν(h(I∗)Ih(I)I∗))+S∗2f(I∗)(∫h0φ(τ)ν(S2I∗f(I(t−τ))S∗2f(I∗)I)dτ−nν(S2f(I)I∗S∗2f(I∗)I))+S∗2g(B∗)n(∫h0φ(τ)ν(B∗h(I(t−τ))h(I∗)B)dτ−nν(h(I)B∗h(I∗)B))≤μ1S∗1(2−S1S∗1−S∗1S1)+(μ2−b)S∗2(3−S2S∗2−S1S∗2S∗1S2−S∗1S1)+bS∗2(2−S1S∗2S∗1S2−S∗1S2S1S∗2)+nS∗2f(I∗)(f(I)f(I∗)−1)(1−f(I∗)If(I)I∗)+S∗2g(B∗)((g(B)g(B∗)−1)(1−g(B∗)Bg(B)B∗)+(h(I)h(I∗)−1)(1−h(I∗)Ih(I)I∗)) |
By Lemma 2.1, it can conclude that
\frac{dL}{dt} = \frac{dL_1}{dt}+\frac{dL_2}{dt}+\frac{dL_3}{dt}\leq 0. |
The equality \frac{dL}{dt} = 0 suggests that \frac{S_1^*}{S_1} = 1 , \frac{S_1S_2^*}{S_1^*S_2} = 1 , \frac{f(I)}{f(I^*)} = 1 and \frac{g(B)}{g(B^*)} = 1 , it implies that E_* is the maximum invariant set of system (2.1) in the set \{\frac{dL}{dt} = 0\} . Using Lemma 4.1, the endemic steady state E_* is globally asymptotically stable.
By choosing a specific kernel function and some infection functions, system (2.1) can be evolved into different dynamic models with time delays. In this section, some such examples are used to further illustrate theoretical results.
Example: Consider the kernel function \rho(\tau) = \delta(\tau-\tau_0) , where \delta is the Dirac delta function. Then system (2.1) can be rewritten as
\begin{equation} \begin{cases} \frac{dS_1}{dt} = A+bS_2-\mu_1S_1-\sigma_1S_1,\\ \frac{dS_2}{dt} = \sigma_1S_1-\alpha S_2 B-\beta S_2f(I(t-\tau_0))e^{-(\mu_2+c)\tau_0}-\mu_2S_2, \\ \frac{dI}{dt} = \alpha S_2 B+\beta S_2f(I(t-\tau_0))e^{-(\mu_2+c)\tau_0}-(\mu_2+c)I,\\ \frac{dB}{dt} = kI(t-\tau_0)e^{-(\mu_2+c)\tau_0}-dB. \end{cases} \end{equation} | (5.1) |
Where A, b, \mu_1, \mu_2, \sigma_1, \alpha, \beta, c, k and d > 0 , g(B) = \alpha B and h(I) = kI . If f(I) = \beta I , it is easy to verify that f(I), g(B) and h(I) satisfy the assumptions (H_1)-(H_3) . Therefore, the basic reproduction number of system (5.1) can be defined as
\begin{equation*} R_0 = \frac{\beta S_2^0e^{-(\mu_2+c)\tau_0}}{\mu_2+c}+\frac{k\alpha S_2^0e^{-(\mu_2+c)\tau_0}}{d(\mu_2+c)}. \end{equation*} |
In this case, if R_0\leq1 , the disease-free equilibrium is globally asymptotically stable and the endemic equilibrium is also globally asymptotically stable if R_0 > 1 (see Figures 1 and 2). As is shown in Figure 2, initial values and time delays can not change dynamic properties of system (2.1) when R_0 is greater than one or less than or equal to one. However, from Figure 2b, it is easy to find that time delay \tau_0 has a significant impact on the positive equilibrium of system (5.1).
If f(I) = \beta I^2 [28], then f(I) does not satisfy hypothesis (H_2) , that is, \frac{f(I)}{I} is a monotonically increasing function. From Figure 3, it is easy to see that the system (5.1) appears periodic oscillation behavior under certain conditions.
In this article, a general S_1S_2IB dynamics model with distributed time delay for animal brucellosis is formulated. Under the assumptions of general biological significance, the non-negative and boundedness of the solution of system (2.1) is first proved. And then the global dynamics of the steady-state solution of system (2.1) are analyzed by constructing Lyapunov functional, it is found that the dynamic properties of equilibria depend on the basic reproduction number R_0 : If R_0\leq 1 , animal brucellosis will eventually die out regardless of the initial value; and if R_0 > 1 , the spread of animal brucellosis is persistent and it eventually reaches the endemic steady state. These results imply that distributed time delay does not change the dynamic properties of system (2.1) when R_0 is greater than one or less than or equal to one. Finally, similar to numerical methods in [35,36], the stability results and other dynamic behaviors are further illustrated through numerical simulation, it turns out that the system experiences periodic oscillations if the assumption (H_2) is not satisfied. In other words, the system (2.1) may exhibit more complex dynamic behaviors if the function \frac{f}{I}, \frac{g}{B} or \frac{h}{I} is monotonically increasing. In these cases, the impact of distributed time delay on other dynamical behaviors of system (2.1) is not completely clear, then we leave these for future work.
This research is partially supported by the National Youth Science Foundation of China (11501528), the National Sciences Foundation of China(11571324) and the Fund for Shanxi "1331KIRT".
The authors declare that they have no competing interests.
[1] |
A. Quarteroni, A. Manzoni, C. Vergara, The cardiovascular system: mathematical modelling, numerical algorithms and clinical applications, Acta Numer., 26 (2017), 365–590. http://doi.org/10.1017/S0962492917000046 doi: 10.1017/S0962492917000046
![]() |
[2] | A. Quarteroni, L. Dede', A. Manzoni, C. Vergara, Mathematical modelling of the human cardiovascular system. data, numerical approximation, clinical applications, Cambridge: Cambridge University Press, 2019. https://doi.org/10.1017/9781108616096 |
[3] |
P. Hauseux, J. Hale, S. Cotin, S. P. A. Bordas, Quantifying the uncertainty in a hyperelastic soft tissue model with stochastic parameters, Appl. Math. Model., 62 (2018), 86–102. https://doi.org/10.1016/j.apm.2018.04.021 doi: 10.1016/j.apm.2018.04.021
![]() |
[4] |
J. Campos, J. Sundnes, R. Dos Santos, B. M. Rocha, Uncertainty quantification and sensitivity analysis of left ventricular function during the full cardiac cycle, Phil. Trans. R. Soc. A, 378 (2020), 20190381. https://doi.org/10.1098/rsta.2019.0381 doi: 10.1098/rsta.2019.0381
![]() |
[5] |
L. Marx, M. A. F. Gsell, A. Rund, F. Caforio, A. J. Prassl, G. Toth-Gayor, et al., Personalization of electro-mechanical models of the pressure-overloaded left ventricle: fitting of windkessel-type afterload models, Phil. Trans. R. Soc. A, 378 (2020), 20190342. https://doi.org/10.1098/rsta.2019.0342 doi: 10.1098/rsta.2019.0342
![]() |
[6] |
R. Rodríguez-Cantano, J. Sundnes, M. Rognes, Uncertainty in cardiac myofiber orientation and stiffnesses dominate the variability of left ventricle deformation response, Int. J. Numer. Meth. Biomed. Eng., 35 (2020), e3178. https://doi.org/10.1002/cnm.3178 doi: 10.1002/cnm.3178
![]() |
[7] |
S. Pagani, A. Manzoni, Enabling forward uncertainty quantification and sensitivity analysis in cardiac electrophysiology by reduced order modeling and machine learning, Int. J. Numer. Meth. Biomed. Eng., 37 (2021), e3450. https://doi.org/10.1002/cnm.3450 doi: 10.1002/cnm.3450
![]() |
[8] | A. Quarteroni, A. Manzoni, F. Negri, Reduced basis methods for partial differential equations: an Introduction, Cham: Springer, 2016. https://doi.org/10.1007/978-3-319-15431-2 |
[9] |
P. Benner, S. Gugercin, K. Willcox, A survey of projection-based model reduction methods for parametric dynamical systems, SIAM Rev., 57 (2015), 483–531. https://doi.org/10.1137/130932715 doi: 10.1137/130932715
![]() |
[10] |
M. Morales, M. van den Boomen, C. Nguyen, J. Kalpathy-Cramer, B. R. Rosen, C. M. Stultz, et al., Deepstrain: A deep learning workflow for the automated characterization of cardiac mechanics, Front. Cardiovasc. Med., 8 (2021), 1041. https://doi.org/10.3389/fcvm.2021.730316 doi: 10.3389/fcvm.2021.730316
![]() |
[11] |
P. Di Achille, A. Harouni, S. Khamzin, O. Solovyova, J. J. Rice, V. Gurev, Gaussian process regressions for inverse problems and parameter searches in models of ventricular mechanics, Front. Physiol., 9 (2018), 1002. https://doi.org/10.3389/fphys.2018.01002 doi: 10.3389/fphys.2018.01002
![]() |
[12] |
U. Noè, A. Lazarus, H. Gao, V. Davies, B. Macdonald, K. Mangion, et al., Gaussian process emulation to accelerate parameter estimation in a mechanical model of the left ventricle: a critical step towards clinical end-user relevance, J. R. Soc. Interface, 16 (2019), 20190114. https://doi.org/10.1098/rsif.2019.0114 doi: 10.1098/rsif.2019.0114
![]() |
[13] |
S. Longobardi, A. Lewalle, S. Coveney, I. Sjaastad, E. K. S. Espe, W. E. Louch, et al., Predicting left ventricular contractile function via gaussian process emulation in aortic-banded rats, Phil. Trans. R. Soc. A, 378 (2019), 20190334. https://doi.org/10.1098/rsta.2019.0334 doi: 10.1098/rsta.2019.0334
![]() |
[14] |
S. Buoso, T. Joyce, S. Kozerke, Personalising left-ventricular biophysical models of the heart using parametric physics-informed neural networks, Med. Image Anal., 71 (2021), 102066. https://doi.org/10.1016/j.media.2021.102066 doi: 10.1016/j.media.2021.102066
![]() |
[15] |
D. Bonomi, A. Manzoni, A. Quarteroni, A matrix DEIM technique for model reduction of nonlinear parametrized problems in cardiac mechanics, Comput. Meth. Appl. Mech. Eng., 324 (2017), 300–326. https://doi.org/10.1016/j.cma.2017.06.011 doi: 10.1016/j.cma.2017.06.011
![]() |
[16] | A. Manzoni, D. Bonomi, A. Quarteroni, Reduced order modeling for cardiac electrophysiology and mechanics: New methodologies, challenges and perspectives, In: Mathematical and numerical modeling of the cardiovascular system and applications, Cham: Springer, 2018,115–166. https://doi.org/10.1007/978-3-319-96649-6_6 |
[17] |
D. Chapelle, A. Gariah, P. Moireau, J. Sainte-Marie, A Galerkin strategy with proper orthogonal decomposition for parameter-dependent problems – analysis, assessments and applications to parameter estimation, ESAIM: M2AN, 47 (2013), 1821–1843. https://doi.org/10.1051/m2an/2013090 doi: 10.1051/m2an/2013090
![]() |
[18] |
M. Pfaller, M. Cruz Varona, J. Lang, C. Bertoglio, W. A. Wall, Using parametric model order reduction for inverse analysis of large nonlinear cardiac simulations, Int. J. Numer. Meth. Biomed. Eng., 36 (2020), e3320. https://doi.org/10.1002/cnm.3320 doi: 10.1002/cnm.3320
![]() |
[19] | J. M. Guccione, A. D. McCulloch, Finite element modeling of ventricular mechanics, In: Theory of heart, New York, NY: Springer, 1991,121–144. https://doi.org/10.1007/978-1-4612-3118-9_6 |
[20] |
J. Guccione, A. McCulloch, Mechanics of active contraction in cardiac muscle: Part Ⅰ–Constitutive relations for fiber stress that describe deactivation, J. Biomech. Eng., 115 (1993), 72–81. http://doi.org/10.1115/1.2895473 doi: 10.1115/1.2895473
![]() |
[21] |
K. Costa, J. Holmes, A. McCulloch, Modelling cardiac mechanical properties in three dimensions, Phil. Trans. R. Soc. A, 359 (2001), 1233–1250. https://doi.org/10.1098/rsta.2001.0828 doi: 10.1098/rsta.2001.0828
![]() |
[22] |
G. Holzapfel, R. Ogden, Constitutive modelling of passive myocardium: a structurally based framework for material characterization, Phil. Trans. R. Soc. A, 367 (2009), 3445–3475. https://doi.org/10.1098/rsta.2009.0091 doi: 10.1098/rsta.2009.0091
![]() |
[23] |
D. Nordsletten, S. Niederer, M. Nash, P. J. Hunter, N. P. Smith, Coupling multi-physics models to cardiac mechanics, Prog. Biophys. Mol. Biol., 104 (2011), 77–88. https://doi.org/10.1016/j.pbiomolbio.2009.11.001 doi: 10.1016/j.pbiomolbio.2009.11.001
![]() |
[24] |
H. M. Wang, H. Gao, X. Y. Luo, C. Berry, B. E. Griffith, R. W. Ogden, et al., Structure-based finite strain modelling of the human left ventricle in diastole, Int. J. Numer. Meth. Biomed. Eng., 29 (2013), 83–103. https://doi.org/10.1002/cnm.2497 doi: 10.1002/cnm.2497
![]() |
[25] |
L. Barbarotta, S. Rossi, L. Dede', A. Quarteroni, A transmurally heterogeneous orthotropic activation model for ventricular contraction and its numerical validation, Int. J. Numer. Meth. Biomed. Eng., 34 (2018), e3137. https://doi.org/10.1002/cnm.3137 doi: 10.1002/cnm.3137
![]() |
[26] | L. Dede', A. Gerbi, A. Quarteroni, Segregated algorithms for the numerical simulation of cardiac electromechanics in the left human ventricle, In: The mathematics of mechanobiology, Cham: Springer, 2020, 81–116. https://doi.org/10.1007/978-3-030-45197-4_3 |
[27] |
D. Lin, F. Yin, A multiaxial constitutive law for mammalian left ventricular myocardium in steady-state barium contracture or tetanus, J. Biomech. Eng., 120 (1998), 504–517. https://doi.org/10.1115/1.2798021 doi: 10.1115/1.2798021
![]() |
[28] |
D. Ambrosi, S. Pezzuto, Active stress vs. active strain in mechanobiology: constitutive issues, J. Elast., 107 (2012), 199–212. https://doi.org/10.1007/s10659-011-9351-4 doi: 10.1007/s10659-011-9351-4
![]() |
[29] | G. A. Holzapfel, Nonlinear solid mechanics. A continuum approach for engineering, Chichester: Wiley, 2001. |
[30] |
C. G. Broyden, A class of methods for solving nonlinear simultaneous equations, Math. Comp., 19 (1965), 577–593. https://doi.org/10.1090/S0025-5718-1965-0198670-6 doi: 10.1090/S0025-5718-1965-0198670-6
![]() |
[31] |
A. Radermacher, S. Reese, POD-based model reduction with empirical interpolation applied to nonlinear elasticity, Int. J. Numer. Meth. Eng., 107 (2016), 477–495. https://doi.org/10.1002/nme.5177 doi: 10.1002/nme.5177
![]() |
[32] |
J. Edmonds, Matroids and the greedy algorithm, Mathematical Programming, 1 (1971), 127–136. https://doi.org/10.1007/BF01584082 doi: 10.1007/BF01584082
![]() |
[33] | C. Farhat, S. Grimberg, A. Manzoni, A. Quarteroni, Computational bottlenecks for proms: precomputation and hyperreduction, In: Model order reduction, Berlin: De Gruyter, 2020,181–244. https://doi.org/10.1515/9783110671490-005 |
[34] |
M. Barrault, Y. Maday, N. Nguyen, A. T. Patera, An 'empirical interpolation' method: application to efficient reduced-basis discretization of partial differential equations, C. R. Math., 339 (2004), 667–672. https://doi.org/10.1016/j.crma.2004.08.006 doi: 10.1016/j.crma.2004.08.006
![]() |
[35] |
M. Grepl, Y. Maday, N. Nguyen, A. T. Patera, Efficient reduced-basis treatment of nonaffine and nonlinear partial differential equations, ESAIM: M2AN, 41 (2007), 575–605. https://doi.org/10.1051/m2an:2007031 doi: 10.1051/m2an:2007031
![]() |
[36] |
Y. Maday, N. Nguyen, A. Patera, S. H. Pau, A general multipurpose interpolation procedure: The magic points, Commun. Pure Appl. Anal., 8 (2008), 383–404. https://doi.org/10.3934/cpaa.2009.8.383 doi: 10.3934/cpaa.2009.8.383
![]() |
[37] |
S. Chaturantabut, D. Sorensen, Nonlinear model reduction via discrete empirical interpolation, SIAM J. Sci. Comput., 32 (2010), 2737–2764. https://doi.org/10.1137/090766498 doi: 10.1137/090766498
![]() |
[38] |
D. Wirtz, D. Sorensen, B. Haasdonk. A posteriori error estimation for DEIM reduced nonlinear dynamical systems, SIAM J. Sci. Comput., 36 (2014), A311–A338. https://doi.org/10.1137/120899042 doi: 10.1137/120899042
![]() |
[39] |
K. Carlberg, R. Tuminaro, P. Boggs, Preserving lagrangian structure in nonlinear model reduction with application to structural dynamics, SIAM J. Sci. Comput., 37 (2015), B153–B184. https://doi.org/10.1137/140959602 doi: 10.1137/140959602
![]() |
[40] |
F. Negri, A. Manzoni, D. Amsallem, Efficient model reduction of parametrized systems by matrix discrete empirical interpolation, J. Comput. Phys., 303 (2015), 431–454. https://doi.org/10.1016/j.jcp.2015.09.046 doi: 10.1016/j.jcp.2015.09.046
![]() |
[41] |
P. Astrid, S. Weiland, K. Willcox, T. Backx, Missing point estimation in models described by proper orthogonal decomposition, IEEE Trans. Automat. Contr., 53 (2008), 2237–2251. https://doi.org/10.1109/TAC.2008.2006102 doi: 10.1109/TAC.2008.2006102
![]() |
[42] |
K. Carlberg, C. Bou-Mosleh, C. Farhat, Efficient non-linear model reduction via a least-squares petrov–galerkin projection and compressive tensor approximations, Int. J. Numer. Meth. Eng., 86 (2011), 155–181. https://doi.org/10.1002/nme.3050 doi: 10.1002/nme.3050
![]() |
[43] |
M. Drohmann, B. Haasdonk, M. Ohlberger, Reduced basis approximation for nonlinear parametrized evolution equations based on empirical operator interpolation, SIAM J. Sci. Comput., 34 (2012), A937–A969. https://doi.org/10.1137/10081157X doi: 10.1137/10081157X
![]() |
[44] |
D. Amsallem, M. Zahr, C. Farhat, Nonlinear model order reduction based on local reduced-order bases, Int. J. Numer. Meth. Eng., 92 (2012), 891–916. https://doi.org/10.1002/nme.4371 doi: 10.1002/nme.4371
![]() |
[45] | P. Tiso, D. J. Rixen, Discrete empirical interpolation method for finite element structural dynamics, In: Topics in nonlinear dynamics, volume 1, New York, NY: Springer, 2013,203–212. https://doi.org/10.1007/978-1-4614-6570-6_18 |
[46] |
F. Ghavamian, P. Tiso, A. Simone, POD-DEIM model order reduction for strain-softening viscoplasticity, Comput. Meth. Appl. Mech. Eng., 317 (2017), 458–479. https://doi.org/10.1016/j.cma.2016.11.025 doi: 10.1016/j.cma.2016.11.025
![]() |
[47] |
S. Land, V. Gurev, S. Arens, C. M. Augustin, L. Baron, R. Blake, et al., Verification of cardiac mechanics software: benchmark problems and solutions for testing active and passive material behaviour, Proc. R. Soc. A, 471 (2015), 20150641. https://doi.org/10.1098/rspa.2015.0641 doi: 10.1098/rspa.2015.0641
![]() |
[48] |
D. Arndt, W. Bangerth, B. Blais, T. C. Clevenger, M. Fehling, A. V. Grayver, et al., The deal.Ⅱ library, version 9.2, J. Numer. Math., 28 (2020), 131–146. https://doi.org/10.1515/jnma-2020-0043 doi: 10.1515/jnma-2020-0043
![]() |
[49] |
J. Bayer, R. Blake, G. Plank, N. A. Trayanova, A novel rule-based algorithm for assigning myocardial fiber orientation to computational heart models, Ann. Biomed. Eng., 40 (2012), 2243–2254. https://doi.org/10.1007/s10439-012-0593-5 doi: 10.1007/s10439-012-0593-5
![]() |
[50] |
J. Wong, E. Kuhl, Generating fibre orientation maps in human heart models using poisson interpolation, Comput. Meth. Biomech. Biomed. Eng., 17 (2014), 1217–1226. https://doi.org/10.1080/10255842.2012.739167 doi: 10.1080/10255842.2012.739167
![]() |
[51] |
R. Doste, D. Soto-Iglesias, G. Bernardino, A. Alcaine, R. Sebastian, S. Giffard-Roisin, et al., A rule-based method to model myocardial fiber orientation in cardiac biventricular geometries with outflow tracts, Int. J. Numer. Meth. Biomed. Eng., 35 (2019), e3185. http://doi.org/10.1002/cnm.3185 doi: 10.1002/cnm.3185
![]() |
[52] |
R. Piersanti, P. Africa, M. Fedele, C. Vergara, L. Dedè, A. F. Corno, et al., Modeling cardiac muscle fibers in ventricular and atrial electrophysiology simulations, Comput. Meth. Appl. Mech. Eng., 373 (2021), 113468. https://doi.org/10.1016/j.cma.2020.113468 doi: 10.1016/j.cma.2020.113468
![]() |
[53] |
S. Rossi, T. Lassila, R. Ruiz-Baier, A. Sequeira, A. Quarteroni, Thermodynamically consistent orthotropic activation model capturing ventricular systolic wall thickening in cardiac electromechanics, Eur. J. Mech-A/Solids, 48 (2014), 129–142. https://doi.org/10.1016/j.euromechsol.2013.10.009 doi: 10.1016/j.euromechsol.2013.10.009
![]() |
[54] | A. Saltelli, M. Ratto, T. Andres, F. Campolongo, J. Cariboni, D. Gatelli, et al., Global sensitivity analysis: the primer, John Wiley & Sons, 2008. https://doi.org/10.1002/9780470725184 |
[55] |
F. Regazzoni, L. Dede', A. Quarteroni, Machine learning of multiscale active force generation models for the efficient simulation of cardiac electromechanics, Comput. Meth. Appl. Mech. Eng., 370 (2020), 113268. https://doi.org/10.1016/j.cma.2020.113268 doi: 10.1016/j.cma.2020.113268
![]() |
[56] | Zygote solid 3d heart generation ii development report, Zygote Media Group Inc., 2014. |
[57] |
F. Regazzoni, M. Salvador, P. C. Africa, M. Fedele, L. Dedè, A. Quarteroni, A cardiac electromechanical model coupled with a lumped-parameter model for closed-loop blood circulation, J. Comput. Phys., 457 (2022), 111083. https://doi.org/10.1016/j.jcp.2022.111083 doi: 10.1016/j.jcp.2022.111083
![]() |
[58] |
F. Regazzoni, L. Dede', A. Quarteroni, Biophysically detailed mathematical models of multiscale cardiac active mechanics, PLoS Comput. Biol., 16 (2020), e1008294. https://doi.org/10.1371/journal.pcbi.1008294 doi: 10.1371/journal.pcbi.1008294
![]() |
[59] |
N. Halko, P. Martinsson, J. Tropp, Finding structure with randomness: Probabilistic algorithms for constructing approximate matrix decompositions, SIAM Rev., 53 (2011), 217–288. https://doi.org/10.1137/090771806 doi: 10.1137/090771806
![]() |