Mathematical modeling of cardiac function to evaluate clinical cases in adults and children
Authors:
Selim Bozkurt aff001
Authors place of work:
Institute of Cardiovascular Science, University College London, London, United Kingdom
aff001
Published in the journal:
PLoS ONE 14(10)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0224663
Summary
Time-varying elastance models can simulate only the pressure and volume signals in the heart chambers while the diagnosis of clinical cases and evaluation of different treatment techniques require more information. In this study, an extended model utilizing the geometric dimensions of the heart chambers was developed to describe the cardiac function. The new cardiac model was evaluated by simulating a healthy and dilated cardiomyopathy (DCM) condition for adults and children. The left ventricular ejection fraction, end-diastolic volume, end-diastolic diameter and diastolic sphericity index were 53.60%, 125 mL, 5.08 cm and 1.82 in the healthy adult cardiovascular system model and 23.70%, 173 mL, 6.60 cm and 1.40 in the DCM adult cardiovascular system model. In the healthy child cardiovascular system model, the left ventricular ejection fraction, end-diastolic volume, end-diastolic diameter and diastolic sphericity index were 59.70%, 92 mL, 4.10 cm and 2.26 respectively and 30.70%, 125 mL, 4.94 cm and 1.87 in the DCM child cardiovascular system model. The developed cardiovascular system model simulates the hemodynamic variables and clinical diagnostic indicators within the physiological range for healthy and DCM conditions proving the feasibility of this new model to evaluate clinical cases in adults and children.
Keywords:
Simulation and modeling – Child health – blood pressure – Ejection fraction – Cardiovascular physiology – Cardiac ventricles – Cardiac atria
Introduction
Numerical modeling of cardiac physiology plays an important role to evaluate clinical scenarios and outcomes of different treatment techniques before experimental or clinical applications. Nonetheless, the information obtained from a numerical simulation depends on the modeled relations in the utilized model. For instance, time-varying elastance models describe the pressure-volume relations in the heart chambers using an elastance function which changes with respect to time over a cardiac cycle [1]. The time-varying elastance models have been used for different clinical scenarios such as evaluating heart failure [2], left ventricular assist device support [3] or simulating the interaction between the cardiovascular and respiratory systems [4,5]. Although time-varying elastance models are widely used to model the cardiac function, they simulate only the blood pressure and volume in heart chambers.
More detailed numerical models such as single fiber contraction models describe the cardiac function by simulating the contraction of a sarcomere over a cardiac cycle [6]. The single fiber contraction models allow simulating the fiber tension in a heart chamber, in addition to the blood pressure and volume [7]. The single fiber contraction models have also been used to simulate physiological and clinical scenarios such as the influence of intra-myocardial pressure on the coronary arterial blood flow rate [8], simulation of fetal heart rate variability [9] or evaluation of rotary blood pump support [10] etc. Although these models are driven by the fiber contraction and provide more information about the cardiac physiology, they also remain insufficient to simulate different mechanisms occurring at each level of cardiac contraction and do not provide information about clinical indicators such as heart chamber dimensions for the physiological cases.
Multi-scale models of the cardiac dynamics simulating the physiological processes at the cellular, protein and organ levels also were developed and used to understand the pathophysiology of heart failure [11]. However, increasing the complexity of a numerical model may not improve the outcome of a simulation. As shown in [12] relatively simple models describing the cardiac function, simulate the physiological and clinical cases such as heart failure and mechanical circulatory support more accurately with respect to the multiscale models because the small discrepancies between the real physiology and the numerical model in each scale cause large deviations at the organ level. Moreover, lumped parameter modeling mostly was used to simulate the adult cardiovascular system. Numerical simulation of cardiac function in children remains to work further to evaluate different physiological scenarios as only a small number of studies focus on pediatric cases [13]. A detailed review of lumped parameter models is given in [14].
The similarities and differences between children and adult physiologies have been subject to developmental theories [15]. Hemodynamic parameters and size of the heart in children reaches similar values of the adult cardiovascular system in time. Therefore, a numerical model capable of simulating the dimensions of the heart chambers and related parameters used to diagnose clinical cases and describe the cardiac function not only for adults but also for children can be utilized to evaluate clinical cases and effects of different treatment in pediatric patients at different ages.
The aim of this study is to develop a lumped parameter model describing the cardiac function and circulatory system and simulating the diagnostic criteria for the clinical cases, therefore, can be used to evaluate physiological scenarios and treatment techniques for cardiovascular system diseases in adults and children at different ages.
Methods
The mathematical model proposed in this study describes the relationship between the pressure and radius in the heart chambers. Left ventricular geometry resembles a truncated ellipsoid [16], therefore, an ellipsoidal geometry was used for the modeling of the left ventricle. The left ventricular volume (Vlv) was expressed using the left ventricular radius (rlv), long axis length (llv) and an additional coefficient (Klv) which allows to include effects of the contraction in the long axis and scales the proportion between the left ventricular radius and volume over a cardiac cycle.
Change of the left ventricular volume with respect to time (dVlv/dt) can also be described as the difference between the mitral and aortic valve flow rates (Qmv, Qav).
Change of the left ventricular radius with respect to time (drlv/dt) becomes as given below.
The left ventricular pressure during the active contraction and the relaxation phases (plv,a) is described using an activation function (fact,lv), end-systolic elastance (Ees,lv), the left ventricular radius (rlv) and left ventricular radius at zero pressure volume (rlv,0).
The relationship between the left ventricular pressure and the radius for the passive left ventricular contraction is adapted from [13] using the Eq 1 as the left ventricular volume.
Right ventricle has a more complex shape with respect to left ventricle and resembles a crescentic cross-section and a triangle from the side view [18]. The crescentic cross-section of the right ventricle is more extensive with respect to the circular cross-section of the left ventricle, and the triangular side view of right ventricle seems like an ellipsoid which has a larger volume occupied by left ventricle [19]. Therefore, an ellipsoidal volume trimmed at the long axis and the basal axis of a right ventricle provides a quite good approximation to model the right ventricular volume (Vrv).
In the equations above, Qtv and Qpv are the flow rates through the tricuspid and pulmonary valves. The right ventricular active and passive pressure components (prv,a, prv,p) are modeled as given below. The right ventricular pressure is obtained using active and passive pressure components.
The left and right atrial geometries were modeled using a similar model to the left ventricular geometry. Therefore, the same relations for the left and right atrial volumes and radiuses are used as in the left ventricle model with different parameter values. The atrial pressure-radius relationship is modeled adopting a time-varying elastance pressure-volume and expressing the atrial volume as a truncated ellipsoidal shape.
The right atrial pressure-radius relations are described in the same way as the left atrial pressure-radius relations. The parameter values used in the adult atria and ventricle models are given in Table 1.
The longitudinal lengths of the left and right ventricles and atria were selected considering the anatomical ranges given in [20–23] for the adult cardiovascular system model. K values were adjusted to simulate the end-systolic and end-diastolic heart chamber radiuses (r) within the reference range for the physiological heart chamber pressures and volumes. The zero pressure volume values of the heart chambers (V0) reported in the literature [17,24] used in the healthy adult cardiovascular system model. The left and right ventricular end-systolic elastance values (Ees) were selected as 2.5 mmHg/mL and 1 mmHg/mL similar to the reported values in [25,26]. The coefficients in the passive component of the ventricular pressure signals (A, B) were 1 and 0.02 [13]. The parameters describing the onset of the ventricular and atrial contraction and relaxation times (T1, T2, Ta, T, D) were selected considering the duration of the cardiac cycle phases [27].
The circulation loop was described using a lumped parameter model which includes electrical analogues for the blood vessel resistance, compliance and inertance. The heart valves were modeled as ideal diodes allowing one-way blood flow. The time-varying geometric models of the heart chambers and electric-analogue of the circulatory loop are given in Fig 1. The resistance, compliance and inertance values in the circulatory system were adjusted manually within the reported physiological range for the healthy adult cardiovascular system [28–30]. The parameter values used in the adult circulatory loop are given in Table 2.
The healthy child cardiovascular system model was simulated by modifying the parameters used in the healthy adult cardiovascular system model. The model parameters were adjusted to simulate children with 1 m2 body surface area and 8–12 age interval [31]. The systolic elastance values in the child cardiovascular system model were selected, considering the data given in [32]. The ventricular long axis length is taken from [33] for the simulated age and body surface area. The zero pressure volumes of the heart chambers (V0) were decreased in the child cardiovascular system model [34]. The compliance values of the blood vessels (C) were reduced in the children as the vessel compliance increases with the increasing body surface area and age [35]. The resistances of the blood vessels (R) are kept the same in the child cardiovascular system model as in the adult cardiovascular system model because the difference between the adult vascular resistances and children vascular resistances for the considered ages is not very high and the vascular resistances vary in a very high variation range in children [36]. Duration of a cardiac cycle (T) was decreased to 0.75 bpm in the child cardiovascular system model as the heart rate increases with the decreasing age [36]. The total blood volume (Vtot) is significantly less in children with relatively small body surface area with respect to the adult cardiovascular system [37]. The total blood volume in the children having 1 m2 body surface area is around 2000 mL [37]. The parameter values used in the heart chambers and circulatory loop of the child cardiovascular system model are given in Tables 3 and 4.
DCM is a condition which left ventricle is dilated and weakened due to contractile dysfunction [38]. Diastolic dysfunction may also coexist in patients with ejection fraction less than 45 per cent [39]. The systemic arterial pressure is increased to maintain the perfusion pressure in DCM patients. The contractile dysfunction was simulated in the DCM left ventricle reducing end-systolic elastance of the left ventricle (Ees,lv) in both adult and child cardiovascular system models. The diastolic dysfunction in the left ventricle was simulated by decreasing the coefficient (A) used in the Eq 8. The systemic arteriolar resistance was increased to 1.6 mmHg/mL/s from 0.95 mmHg/mL/s to simulate the increased resistance in the systemic circulation in both adult and child cardiovascular system models. The scaling coefficient (Klv) was also increased to simulate the remodeling in the left ventricular wall and the enlargement of the left ventricular cavity. The modified parameters in the left ventricle and circulatory system for both adult and child cardiovascular system models are given in the Tables 1–4 in the brackets.
Fractional shortening (fs) and sphericity index (si) were used to evaluate healthy and DCM conditions in the adult and child cardiovascular system models along with the ejection fraction, volumes and dimensions of the heart chambers.
Here, EDD and ESD are the diameters at the end of diastole and systole, l is the length of a heart chamber and Dmid is the diameter of the mid-cavity of the heart chambers. The right ventricular basal diameter was assumed to be the same as the right ventricular radius because of the anatomical geometry of right ventricle. The simulations were performed using Matlab Simulink R2017a. The set of equations was solved using the ode15s solver. The maximum step size was 1e-3 s, relative tolerance was set to 1e-3.
Results
First, the simulations were performed for the healthy adult and child cardiovascular system models. The ventricular, atrial and aortic pressures, the ventricular and atrial volumes and diameters in the healthy adult and child cardiovascular system models are given in Fig 2.
The left ventricular and aortic systolic pressures in the healthy adult cardiovascular system were around 120 mmHg. In both adult and child cardiovascular system models, the right ventricular volumes remain higher with respect to the left ventricular volumes over a cardiac cycle. The right atrial volumes are higher with respect to the left atrial volumes over a cardiac cycle in both models. Similarly, the right atrial diameters are higher as well with respect to the left atrial diameters over a cardiac cycle. The right ventricular basal diameter remains lower with respect to the left ventricular basal diameter similar to the normal physiological conditions in the adults and children. The ventricular and atrial pressures, aortic pressure, ventricular and atrial volumes and the diameters in the DCM adult and child cardiovascular system models are given in Fig 3.
The left ventricular and aortic systolic pressures decreased in the DCM cases for the adult and child cardiovascular system models. The left ventricular volume increased significantly simulating the enlargement because of the impaired contractility in the DCM left ventricles. Left atrial volumes in DCM models increased with respect to the left atrial volumes in the healthy cardiovascular system models. The left ventricular and atrial diameters increased as well with the increased volumes in the left atria and ventricles of the adult and child DCM cardiovascular system models. The pressure-volume loops for the healthy and DCM models are given in Fig 4.
The left ventricular pressure-volume loop shifted to the right and became smaller in the cardiovascular system models simulating DCM for both adults and children. The pressure-volume areas of the left and right ventricles are smaller in the cardiovascular system models simulating the DCM condition with respect to the cardiovascular system models simulating the healthy condition. Additionally, the pressure-volume loop area is smaller in the child cardiovascular system model with respect to the adult cardiovascular system model. The end-systolic and end-diastolic volumes, diameters and pressures in the heart chambers for the adult and child cardiovascular systems with the healthy and DCM conditions are given in Fig 5, the maximal and minimal dimensions of the heart chambers and the hemodynamic parameters used to evaluate the healthy and DCM conditions are given in Table 5.
The hemodynamic parameters shown in Fig 5 and Table 5 were within the normal physiological ranges for the healthy and DCM condition in the adult and child cardiovascular system models. The left ventricular end-diastolic and end-systolic diameters increased in the DCM cardiovascular system models with respect to the healthy cardiovascular system models for an adult and a child. Left atrial maximal and minimal diameters increased in the DCM cardiovascular system models as well. Right ventricular end-diastolic diameters reduced while end-systolic diameters in the right ventricle do not change significantly for the DCM adult and child cardiovascular system models. Right atrial diameters reduced at the end-systole and end-diastole in the DCM adult and child cardiovascular system models. The cardiac output reduced significantly in the DCM condition for both adult child cardiovascular system models as well as the left ventricular stroke volume and ejection fraction. Fractional shortening, diastolic and systolic sphericity indices and aortic mean pressure decreased in the DCM conditions as well. The mean pulmonary arterial pressure remained at similar levels in the DCM and healthy conditions for both adult and child cardiovascular system models.
Discussion
Time-varying elastance models reveal only the relationship between the pressure and volume signals in the ventricles by lumping all the contraction dynamics in the elastance term over a cardiac cycle. These models have been widely used to evaluate clinical devices or clinical scenarios numerically despite their limited capability [30,40]. The cardiovascular system model proposed in this study describes the relationship between pressure and dimensions of the heart chambers. Also, systolic elastances were used to describe pressures in the ventricle models. Therefore, the proposed numerical model extends the capability of the existing relatively simple time-varying elastance models without losing information. Single fiber contraction models are driven by contraction of muscle fibers, and they simulate length and tension on the fibers, along with the pressure and volume relation in heart chambers. Although single fiber contraction models have not been used as widely as the time-varying elastance models, they have been utilized to evaluate medical devices numerically [7,41]. These models also have limited capabilities in simulating indicators for clinical scenarios as ejection fraction is the only clinical indicator can be obtained. Therefore, their capability also remains limited with respect to the developed cardiovascular system model in this study. Multiscale models simulating the action potentials and Calcium kinetics provide detailed information about the cardiac dynamics at cellular and protein levels for a healthy and DCM heart [11]. However, at the organ level ejection fraction is the only parameter that such models can provide. Moreover, calcium dynamics in children at early ages are quite different [42], therefore, substantial modifications in the equations may be required to simulate heart failure when these models are used to simulate children cardiovascular system. The developed cardiovascular system model simulates cardiac function at organ level and allows simulation of additional clinical indicators such as sphericity index and fractional shortening, which are used in the diagnosis of heart failure. Moreover, the presented cardiovascular system model shows that children have a relatively long longitudinal axis with respect to the minor axis in their left ventricle.
Scaling coefficient (K) which include the effects of the longitudinal contraction and irregularities in the heart chamber geometries and allow simulation of the physiological heart chamber diameter signals was used in the proposed cardiovascular system models. The scaling coefficients (K) also relate the diameter and volume signals in the heart chambers. Modifying these parameters for a heart chamber does not influence the volume and pressure signals while the diameter changes. A higher value of the scaling coefficients (K) shows a stronger longitudinal contraction and a more complex geometry for the heart chambers. Right ventricle has the most complex geometry among all the heart chambers, and right ventricular longitudinal contraction contributes to the ejection of the blood through the pulmonary valve significantly [43]. Therefore, higher values for the right ventricular scaling constant (Krv) with respect to the scaling constants in other heart chambers were used in the adult and child cardiovascular system models. Also, higher scaling constant values for the heart chambers in the child cardiovascular system model with respect to the adult cardiovascular system show stronger longitudinal contractions in children. In the adult and child DCM models, the scaling coefficients (Klv) were modified as well as the systolic elastances (Esys,lv) and zero pressure-volume values (Vlv,0) for the left ventricle because of the remodeling in the myocardial tissue causing in an increase in the left ventricular volume and diameter. The modifications in the systolic elastance (Esys,lv) and zero pressure-volume (Vlv,0) were done to simulate the impaired contractile behavior and enlarged cavity volume in the DCM left ventricles. The scaling constant (Klv) was modified to simulate the changes in the contractile strength of the longitudinal contractions.
Ventricular long axis length in the model was assumed to be constant because the ventricular long axis signal over a cardiac cycle is similar and seems to be proportional to the ventricular volume signal [44]. A similar change in the ventricular diameter is also seen. Taking into account both signals increases the complexity of the model. Scaling coefficients (K) for each heart chambers were used to express the effects of longitudinal contraction and irregularities in the heart chamber geometries and reduce the complexity of the model.
Accuracy of the right ventricular volume estimation from 2D echocardiography requires assumptions on the geometry of right ventricle [45]. Trimmed ellipsoidal geometries similar to the presented right ventricular model in this study to estimate right ventricular volume have already been used in 2D echocardiography [46,47]. Empirical coefficients expressing the irregularities in right ventricular shape are also used in 2D echocardiographic area-length methods estimating right ventricular volume [45]. Using a trimmed ellipsoidal geometry and an additional coefficient (Krv) as in 2D echocardiographic cardiac imaging describes the right ventricular volume signals for the healthy and DCM models accurately.
The pressure and volume signals simulated in the adult cardiovascular system model for the healthy condition remain within the physiological range [27,48–51]. The simulated maximal and minimal heart chamber dimensions in the healthy adult cardiovascular system model are also within the healthy physiological ranges and corresponds to the reported average values for the adults [52–58]. As the dilatation occurs in a DCM left ventricle, the left ventricular volume increases significantly, internal ventricular diameter at the diastole and systole increase around 6.5 cm and 5.5 cm and ejection fraction may reduce below 30 per cent in adults [39,59,60]. The simulation results were in the reported physiological range for the adult DCM cardiovascular system model. Left ventricular geometry assumes a more spherical shape with the remodeling due to the dilation [61]. Reduced sphericity indexes at the end-systole and end-diastole in DCM (Table 5) indicate a more globular shape for the left ventricle [39]. Ejection fraction and fractional shortening are parameters used for diagnosis and risk stratification in DCM patients [39]. Moreover, fractional shortening is used for risk assessment in DCM [62]. Reduced ejection fraction and fractional shortening in the left ventricle (Table 5) indicate the existence of a dilated cardiomyopathy in the adult cardiovascular system model. [39]. The changes occur in the left ventricle also affect the other heart chambers and may increase the left atrial volume over a cardiac cycle or worsening the right ventricular function as in the simulations [63–65].
In pediatric patients, fractional shortening and ejection fraction are the most commonly used parameters to evaluate the left ventricular function [66]. Measurement of fractional shortening is based on the assumption of the cylindrical shape of the left ventricle. However, the assumption of cylindrical shape can be a limitation for using ejection fraction [66]. Relatively left ventricular high sphericity index values in the child cardiovascular system model (Table 5) indicate that the ratio between the long axis and mid-cavity diameter is relatively high. This indicates a close shape to a cylinder near the base of the left ventricle in the child cardiovascular system model. Children’s cardiac output reaches the similar levels to those in the adult cardiovascular system around 15 years while mean arterial pressure remains considerably lower in children at the same age with respect to mean arterial pressure level in adults’ cardiovascular system [36]. In this study, body surface area was selected as 1 m2 corresponding to 8–12 years children’s body surface area [31] and the hemodynamic parameters obtained from the simulations for the child cardiovascular system models compared to this age interval. The average cardiac output at these ages is around 4.5 L/min and changes within a large variation range [36]. The systolic pressure in children’s left ventricle generally remains below 110 mmHg while mean arterial pressure is around 80 mmHg in children circulation [67]. Reported right ventricular end-diastolic and systolic pressures are around 5 mmHg, 28 mmHg respectively and mean pulmonary arterial pressure can reach up to 22 mmHg being slightly lower with respect to those of adult cardiovascular system in the children around 9 years old [68]. The simulated pressures in the healthy child cardiovascular system were similar to those values reported in the literature. Mean of the left and right ventricular end-diastolic volume indexes in children are around 80 mL/m2 with slightly higher right ventricular volume and may reach to values around 100 mL/m2 while the end-systolic volume in the ventricles is around 30 mL/m2 and with 45 mL/m2 upper reference limit [51]. The volume indexes obtained from the simulations for 1 m2 body surface area corresponds well with the reported volume indexes in literature. The left atrial volume signal over a cardiac cycle in the healthy child cardiovascular system corresponds to a reported average left atrial volume signal over a cardiac cycle quite well [34]. Moreover, the reported range of the atrial maximal and minimal volume index values [51] matches with the simulated left and right atrial maximal and minimal volume indexes in the healthy child cardiovascular system for 1 m2 body surface area. Typically, the left ventricular end-diastolic and end-systolic diameters are around 4.1 cm and 2.6 cm, the right ventricular basal diameter at the end of the diastole is around 2 cm and the left atrial maximal diameter 2.9 cm and changing again within a physiologic variation range for the studied age interval [69–71]. The simulated dimensions in the cardiac chambers remain within the physiologic range and match with the reported values in the literature.
The DCM cardiovascular system model simulating the child circulation responded to the changes in parameters in a similar way as the adult DCM cardiovascular system model. Cardiac output and ejection fraction reduced from healthy values to the values characterized as impaired [72,73]. The left ventricular volume signal over a cardiac cycle increased as a result of dilatation and cause dilation in the left atrial volume, as reported in [74]. Left ventricular end-diastolic and end-systolic diameters and right ventricular diameter were similar to patient-specific values reported in [73] and remained within the reported range in the children for the studied age interval [75]. Moreover, the left ventricular fractional shortening in the child DCM cardiovascular system model matched closely to the reported average value of 19 per cent and remained within the physiologic variation range for the DCM condition in children heart [75].
Mitral regurgitation may occur in DCM because increased left ventricular size cause improper leaflet coaptation [38]. Simulation of the regurgitant flow rate through the mitral valve requires a numerical model describing the motion of leaflets [10,17] or zener diode analogue models allowing reverse blood flow. Flow rate signals through the heart valves were simulated using ideal diodes in the developed cardiovascular system model. Therefore, the change of the volume signals in the heart chambers depends only forward flow rate signals. Utilizing an elaborated mitral valve model will allow simulating effects regurgitant flow rate on the left ventricular function as the change of the left ventricular volume signal depends on mitral and aortic valve flow rates regardless of the direction of the blood flow (Eq 4). Nonetheless, ideal diode models have also been used widely to simulate clinical scenarios such as DCM and LVAD therapies in cardiovascular system [76,77].
The developed cardiovascular system model is driven by an activation function which utilizes time signal to simulate the hemodynamic variables. Clinically, wall stress and strain are also commonly used physiologic parameters. Simulation of these parameters requires performing finite element analyses since wall stress and strain have orientations in space and also depend on the geometry of myocardium. Additionally, the interaction between the ventricles through the intraventricular septum can be simulated using finite element analyses. The ventricles interact in the developed cardiovascular system model due to change in the parameters such as contractility as afterload and preload of the heart chambers are influenced because of variation of this parameter. However, it should be noted that the developed cardiovascular system model simulated the hemodynamic parameters for normal and DCM conditions in adults and children within the physiological range. Moreover, lumped parameter models have the advantage of performing simulations in shorter computational times, and they are powerful clinical decision-making tools [14].
The DCM in adult and child cardiovascular system models was simulated by adjusting the parameters in the cardiovascular system models manually. Progression of heart failure depends on complex and interacting mechanisms which are not entirely understood [78]. Modeling the progression of heart failure requires identification of these mechanisms. Therefore, progression of the DCM was not included in the simulations. Similarly, the interacting compensatory mechanisms in heart failure can be implemented after the identification of the exact mechanisms [79].
The cardiovascular system model presented in this study can be used to simulate different clinical cases for different ages and conditions as it allows to simulate clinical indicators extracted from the geometry of heat chambers. In the future, it is aimed to use this cardiovascular system model to simulate heart failure for pediatric patients in patient-specific studies to evaluate mechanical circulatory support. An important portion of the pediatric patients implanted with a heart pump exhibit complications associated with mechanical circulatory support [80]. Continuous and pulsatile flow LVADs are being tested to use in children. Moreover, LVAD success criteria in adults include evaluation of parameters such as ejection fraction, left ventricular end-diastolic volume, left ventricular end-diastolic diameter, left ventricular end-systolic diameter, left atrial volume, right atrial pressure [81]. The developed cardiovascular system model can simulate the dimensions of heart chambers along with the other parameters. The parameter values in the developed model can be tuned using parameter estimation and system identification techniques to simulate patient-specific cardiovascular system dynamics as presented in [4]. Therefore, it could be used to simulate patient-specific scenarios and evaluate devices such as LVADs before experimental tests.
Conclusions
In this study, a mathematical model simulating hemodynamic parameters and clinical indicators in the cardiovascular system was developed to evaluate clinical cases in adults and children. The developed numerical model demonstrated its feasibility to simulate clinical cases and can be used to evaluate clinical cases and treatment techniques in adults and children at different ages.
Supporting information
S1 File [zip]
Cardiovascular system Model.slx, Matlab Simulink Code for model simulation.
Zdroje
1. Lankhaar J-W, Rövekamp FA, Steendijk P, Faes TJC, Westerhof BE, Kind T, et al. Modeling the instantaneous pressure-volume relation of the left ventricle: a comparison of six models. Ann Biomed Eng. 2009;37: 1710–1726. doi: 10.1007/s10439-009-9742-x 19554450
2. Faes TJ, Kerkhof PL. The Volume Regulation Graph versus the Ejection Fraction as Metrics of Left Ventricular Performance in Heart Failure with and without a Preserved Ejection Fraction: A Mathematical Model Study. Clin Med Insights Cardiol. 2015;9: 73–91. doi: 10.4137/CMC.S18748 26052232
3. Song Z, Gu K, Gao B, Wan F, Chang Y, Zeng Y. Hemodynamic effects of various support modes of continuous flow LVADs on the cardiovascular system: A numerical study. Med Sci Monit Int Med J Exp Clin Res. 2014;20: 733–741. doi: 10.12659/MSM.890824 24793178
4. Ellwein LM, Pope SR, Xie A, Batzel JJ, Kelley CT, Olufsen MS. Patient-specific modeling of cardiovascular and respiratory dynamics during hypercapnia. Math Biosci. 2013;241: 56–74. doi: 10.1016/j.mbs.2012.09.003 23046704
5. Gaudenzi F, Avolio AP. Lumped parameter model of cardiovascular-respiratory interaction. Conf Proc Annu Int Conf IEEE Eng Med Biol Soc IEEE Eng Med Biol Soc Annu Conf. 2013;2013: 473–476. doi: 10.1109/EMBC.2013.6609539 24109726
6. van de Vosse FN, Stergiopulos N. Pulse Wave Propagation in the Arterial Tree. Annu Rev Fluid Mech. 2011;43: 467–499. doi: 10.1146/annurev-fluid-122109-160730
7. Cox LGE, Loerakker S, Rutten MCM, de Mol BAJM, van de Vosse FN. A mathematical model to evaluate control strategies for mechanical circulatory support. Artif Organs. 2009;33: 593–603. doi: 10.1111/j.1525-1594.2009.00755.x 19558561
8. Bovendeerd PHM, Borsje P, Arts T, van De Vosse FN. Dependence of Intramyocardial Pressure and Coronary Flow on Ventricular Loading and Contractility: A Model Study. Ann Biomed Eng. 2006;34: 1833–1845. doi: 10.1007/s10439-006-9189-2 17048105
9. Jongen GJLM, van der Hout-van der Jagt MB, Oei SG, van de Vosse FN, Bovendeerd PHM. Simulation of fetal heart rate variability with a mathematical model. Med Eng Phys. 2017;42: 55–64.28196652
10. Bozkurt S, Bozkurt S. In-silico evaluation of left ventricular unloading under varying speed continuous flow left ventricular assist device support. Biocybern Biomed Eng. 2017;37: 373–387. doi: 10.1016/j.bbe.2017.03.003
11. Bhattacharya-Ghosh B, Bozkurt S, Rutten MCM, van de Vosse FN, Díaz-Zuccarini V. An in silico case study of idiopathic dilated cardiomyopathy via a multi-scale model of the cardiovascular system. Comput Biol Med. 2014;53: 141–153. doi: 10.1016/j.compbiomed.2014.06.013 25147131
12. Bozkurt S. In-silico modeling of left ventricle to simulate dilated cardiomyopathy and cf-lvad support. J Mech Med Biol. 2017;17: 1750034. doi: 10.1142/S0219519417500348
13. Sheffer L, Santamore WP, Barnea O. Cardiovascular simulation toolbox. Cardiovasc Eng Dordr Neth. 2007;7: 81–88. doi: 10.1007/s10558-007-9030-z 17570062
14. Shimizu S, Une D, Kawada T, Hayama Y, Kamiya A, Shishido T, et al. Lumped parameter model for hemodynamic simulation of congenital heart diseases. J Physiol Sci JPS. 2018;68: 103–111. doi: 10.1007/s12576-017-0585-1 29270856
15. MacGregor J. By Janet MacGregor—Introduction to the Anatomy and Physiology of Children: A Guide for Students of Nursing, Child Care and Health. 2 edition. ROUTLEDGE; 2008.
16. Udelson JE. Left Ventricular Shape: The Forgotten Stepchild of Remodeling Parameters. JACC Heart Fail. 2017;5: 179–181. doi: 10.1016/j.jchf.2017.01.005 28254123
17. Korakianitis T, Shi Y. A concentrated parameter model for the human cardiovascular system including heart valve dynamics and atrioventricular interaction. Med Eng Phys. 2006;28: 613–628.16293439
18. Walker LA, Buttrick PM. The right ventricle: biologic insights and response to disease: updated. Curr Cardiol Rev. 2013;9: 73–81. doi: 10.2174/157340313805076296 23092273
19. Denslow S, Wiles HB. Right ventricular volumes revisited: a simple model and simple formula for echocardiographic determination. J Am Soc Echocardiogr Off Publ Am Soc Echocardiogr. 1998;11: 864–873.
20. Holland DJ, Sharman JE, Leano RL, Marwick TH. Gender differences in systolic tissue velocity: role of left ventricular length. Eur J Echocardiogr J Work Group Echocardiogr Eur Soc Cardiol. 2009;10: 941–946. doi: 10.1093/ejechocard/jep103 19720628
21. Kossaify A. Echocardiographic Assessment of the Right Ventricle, from the Conventional Approach to Speckle Tracking and Three-Dimensional Imaging, and Insights into the “Right Way” to Explore the Forgotten Chamber. Clin Med Insights Cardiol. 2015;9: 65–75. doi: 10.4137/CMC.S27462
22. Peverill RE, Chou B, Donelan L. Left ventricular long axis tissue Doppler systolic velocity is independently related to heart rate and body size. PloS One. 2017;12: e0173383. doi: 10.1371/journal.pone.0173383 28288162
23. Rudski LG, Lai WW, Afilalo J, Hua L, Handschumacher MD, Chandrasekaran K, et al. Guidelines for the echocardiographic assessment of the right heart in adults: a report from the American Society of Echocardiography endorsed by the European Association of Echocardiography, a registered branch of the European Society of Cardiology, and the Canadian Society of Echocardiography. J Am Soc Echocardiogr Off Publ Am Soc Echocardiogr. 2010;23: 685–713; quiz 786–788. doi: 10.1016/j.echo.2010.05.010 20620859
24. Ursino M. Interaction between carotid baroregulation and the pulsating heart: a mathematical model. Am J Physiol. 1998;275: H1733–1747. doi: 10.1152/ajpheart.1998.275.5.H1733 9815081
25. Heldt T, Shim EB, Kamm RD, Mark RG. Computational modeling of cardiovascular response to orthostatic stress. J Appl Physiol Bethesda Md 1985. 2002;92: 1239–1254. doi: 10.1152/japplphysiol.00241.2001 11842064
26. Pinsky MR. The right ventricle: interaction with the pulmonary circulation. Crit Care. 2016;20: 266. doi: 10.1186/s13054-016-1440-0 27613549
27. Hall JE. Guyton and Hall Textbook of Medical Physiology, 13e. 13 edition. Philadelphia, PA: Saunders; 2015.
28. Blanco PJ, Feijóo RA. A 3D-1D-0D Computational Model for the Entire Cardiovascular System. Mecánica Comput. 2010;29: 5887–5911.
29. Albanese A, Cheng L, Ursino M, Chbat NW. An integrated mathematical model of the human cardiopulmonary system: model development. Am J Physiol Heart Circ Physiol. 2016;310: H899–921. doi: 10.1152/ajpheart.00230.2014 26683899
30. Shi Y, Lawford P, Hose R. Review of Zero-D and 1-D Models of Blood Flow in the Cardiovascular System. Biomed Eng OnLine. 2011;10: 33. doi: 10.1186/1475-925X-10-33 21521508
31. El Edelbi R, Lindemalm S, Eksborg S. Estimation of body surface area in various childhood ages—validation of the Mosteller formula. Acta Paediatr Oslo Nor 1992. 2012;101: 540–544. doi: 10.1111/j.1651-2227.2011.02580.x
32. Kiani A, Shakibi JG. Normal value of left ventricular end-systolic elastance in infants and children. Iran J Med Sci. 2003;28: 169–172.
33. Regen DM, Graham TP, Wyse RK, Deanfield J, Franklin RC. Left-ventricular cavity dimensions in children with normal and dilated hearts. Pediatr Cardiol. 1988;9: 17–24. doi: 10.1007/BF02279878 2964589
34. Poutanen T, Jokinen E, Sairanen H, Tikanoja T. Left atrial and left ventricular function in healthy children and young adults assessed by three dimensional echocardiography. Heart Br Card Soc. 2003;89: 544–549.
35. Gardner AW, Parker DE. Association between arterial compliance and age in participants 9 to 77 years old. Angiology. 2010;61: 37–41. doi: 10.1177/0003319709339588 19638351
36. Cattermole GN, Leung PYM, Ho GYL, Lau PWS, Chan CPY, Chan SSW, et al. The normal ranges of cardiovascular parameters measured using the ultrasonic cardiac output monitor. Physiol Rep. 2017;5. doi: 10.14814/phy2.13195 28320891
37. Raes A, Van Aken S, Craen M, Donckerwolcke R, Vande Walle J. A reference frame for blood volume in children and adolescents. BMC Pediatr. 2006;6: 3. doi: 10.1186/1471-2431-6-3 16503982
38. Japp AG, Gulati A, Cook SA, Cowie MR, Prasad SK. The Diagnosis and Evaluation of Dilated Cardiomyopathy. J Am Coll Cardiol. 2016;67: 2996–3010. doi: 10.1016/j.jacc.2016.03.590 27339497
39. Mathew T, Williams L, Navaratnam G, Rana B, Wheeler R, Collins K, et al. Diagnosis and assessment of dilated cardiomyopathy: a guideline protocol from the British Society of Echocardiography. Echo Res Pract. 2017;4: G1–G13. doi: 10.1530/ERP-16-0037 28592613
40. Jelenc M, Jelenc B, Vrtovec B, Kneževič I. Mitral regurgitation and axial flow left ventricular assist device: a computer simulation study. ASAIO J Am Soc Artif Intern Organs 1992. 2013;59: 405–409. doi: 10.1097/MAT.0b013e31829ff57d 23820280
41. Bozkurt S. Effect of Cerebral Flow Autoregulation Function on Cerebral Flow Rate Under Continuous Flow Left Ventricular Assist Device Support. Artif Organs. 2018;42: 800–813. doi: 10.1111/aor.13148 29726017
42. Porter GA, Makuck RF, Rivkees SA. Intracellular calcium plays an essential role in cardiac development. Dev Dyn. 2003;227: 280–290. doi: 10.1002/dvdy.10307 12761855
43. Brown SB, Raina A, Katz D, Szerlip M, Wiegers SE, Forfia PR. Longitudinal shortening accounts for the majority of right ventricular contraction and improves after pulmonary vasodilator therapy in normal subjects and patients with pulmonary arterial hypertension. Chest. 2011;140: 27–33. doi: 10.1378/chest.10-1136 21106653
44. Veronesi F, Corsi C, Caiani EG, Sarti A, Lamberti C. Tracking of left ventricular long axis from real-time three-dimensional echocardiography using optical flow techniques. IEEE Trans Inf Technol Biomed Publ IEEE Eng Med Biol Soc. 2006;10: 174–181.
45. Ostenfeld E, Flachskampf FA. Assessment of right ventricular volumes and ejection fraction by echocardiography: from geometric approximations to realistic shapes. Echo Res Pract. 2015;2: R1–R11. doi: 10.1530/ERP-14-0077 26693327
46. Jorstig S, Waldenborg M, Lidén M, Wodecki M, Thunberg P. Determination of Right Ventricular Volume by Combining Echocardiographic Distance Measurements. Echocardiogr Mt Kisco N. 2016;33: 844–853. doi: 10.1111/echo.13173 26841195
47. Jorstig S, Waldenborg M, Lidén M, Thunberg P. Right ventricular ejection fraction measurements using two-dimensional transthoracic echocardiography by applying an ellipsoid model. Cardiovasc Ultrasound. 2017;15: 4. doi: 10.1186/s12947-017-0096-5 28270161
48. Fukuta H, Little WC. The cardiac cycle and the physiologic basis of left ventricular contraction, ejection, relaxation, and filling. Heart Fail Clin. 2008;4: 1–11. doi: 10.1016/j.hfc.2007.10.004 18313620
49. Kovalova S, Necas J, Vespalec J. What is a “normal” right ventricle? Eur J Echocardiogr. 2006;7: 293–297. doi: 10.1016/j.euje.2005.06.010
50. Hudsmith LE, Petersen SE, Francis JM, Robson MD, Neubauer S. Normal human left and right ventricular and left atrial dimensions using steady state free precession magnetic resonance imaging. J Cardiovasc Magn Reson Off J Soc Cardiovasc Magn Reson. 2005;7: 775–782.
51. Kawel-Boehm N, Maceira A, Valsangiacomo-Buechel ER, Vogel-Claussen J, Turkbey EB, Williams R, et al. Normal values for cardiovascular magnetic resonance in adults and children. J Cardiovasc Magn Reson Off J Soc Cardiovasc Magn Reson. 2015;17: 29. doi: 10.1186/s12968-015-0111-7 25928314
52. El Missiri AM, El Meniawy KAL, Sakr SAS, Mohamed ASE. Normal reference values of echocardiographic measurements in young Egyptian adults. Egypt Heart J. 2016;68: 209–215. doi: 10.1016/j.ehj.2016.01.002
53. Gibson PH, Becher H, Choy JB. Classification of left ventricular size: diameter or volume with contrast echocardiography? Open Heart. 2014;1: e000147. doi: 10.1136/openhrt-2014-000147 25525505
54. Ho SY, Nihoyannopoulos P. Anatomy, echocardiography, and normal right ventricular dimensions. Heart. 2006;92: i2–i13. doi: 10.1136/hrt.2005.077875 16543598
55. Kossaify A. Echocardiographic Assessment of the Right Ventricle, from the Conventional Approach to Speckle Tracking and Three-Dimensional Imaging, and Insights into the “Right Way” to Explore the Forgotten Chamber. Clin Med Insights Cardiol. 2015;9: 65–75. doi: 10.4137/CMC.S27462
56. Kou S, Caballero L, Dulgheru R, Voilliot D, De Sousa C, Kacharava G, et al. Echocardiographic reference ranges for normal cardiac chamber size: results from the NORRE study. Eur Heart J Cardiovasc Imaging. 2014;15: 680–690. doi: 10.1093/ehjci/jet284 24451180
57. Lang RM, Badano LP, Mor-Avi V, Afilalo J, Armstrong A, Ernande L, et al. Recommendations for cardiac chamber quantification by echocardiography in adults: an update from the American Society of Echocardiography and the European Association of Cardiovascular Imaging. J Am Soc Echocardiogr Off Publ Am Soc Echocardiogr. 2015;28: 1–39.e14. doi: 10.1016/j.echo.2014.10.003 25559473
58. Lang RM, Bierig M, Devereux RB, Flachskampf FA, Foster E, Pellikka PA, et al. Recommendations for chamber quantification. Eur J Echocardiogr J Work Group Echocardiogr Eur Soc Cardiol. 2006;7: 79–108. doi: 10.1016/j.euje.2005.12.014 16458610
59. Gavazzi A, De Maria R, Renosto G, Moro A, Borgia M, Caroli A, et al. The spectrum of left ventricular size in dilated cardiomyopathy: clinical correlates and prognostic implications. SPIC (Italian Multicenter Cardiomyopathy Study) Group. Am Heart J. 1993;125: 410–422. doi: 10.1016/0002-8703(93)90020-a 8427135
60. Matsumura Y, Elliott PM, Mahon NG, Virdee MS, Doi Y, McKenna WJ. Familial dilated cardiomyopathy: assessment of left ventricular systolic and diastolic function using Doppler tissue imaging in asymptomatic relatives with left ventricular enlargement. Heart. 2006;92: 405–406. doi: 10.1136/hrt.2005.065474 16501207
61. Weintraub RG, Semsarian C, Macdonald P. Dilated cardiomyopathy. The Lancet. 2017;390: 400–414. doi: 10.1016/S0140-6736(16)31713-5
62. Niemeijer MN, Leening MJG, van den Berg ME, Hofman A, Franco OH, Deckers JW, et al. Subclinical Abnormalities in Echocardiographic Parameters and Risk of Sudden Cardiac Death in a General Population: The Rotterdam Study. J Card Fail. 2016;22: 17–23. doi: 10.1016/j.cardfail.2015.06.007 26093333
63. Ferreira F, Galrinho A, Soares R, Branco L, Abreu J, Feliciano J, et al. Prognostic value of left atrial volume in patients with dilated cardiomyopathy. Rev Port Cardiol Orgao Of Soc Port Cardiol Port J Cardiol Off J Port Soc Cardiol. 2013;32: 865–872. doi: 10.1016/j.repc.2012.12.017 24119867
64. Gopal D, Wang J, Han Y. Determinants of Normal Left Atrial Volume in Heart Failure with Moderate-to-Severely Reduced Ejection Fraction. In: Cardiology Research and Practice [Internet]. 2018 [cited 4 Nov 2018]. doi: 10.1155/2018/7512758 29850229
65. Parcharidou DG, Giannakoulas G, Efthimiadis GK, Karvounis H, Papadopoulou KN, Dalamanga E, et al. Right ventricular function in ischemic or idiopathic dilated cardiomyopathy. Circ J Off J Jpn Circ Soc. 2008;72: 238–244.
66. Tissot C, Singh Y, Sekarski N. Echocardiographic Evaluation of Ventricular Function—For the Neonatologist and Pediatric Intensivist. Front Pediatr. 2018;6. doi: 10.3389/fped.2018.00006
67. de Simone G, Devereux RB, Kimball TR, Roman MJ, Palmieri V, Celentano A, et al. Relation of heart rate to left ventricular dimensions in normotensive, normal-weight children, adolescents and adults. Ital Heart J Off J Ital Fed Cardiol. 2001;2: 599–604.
68. Krovetz LJ, McLoughlin TG, Mitchell MB, Schiebler GL. Hemodynamic findings in normal children. Pediatr Res. 1967;1: 122–130. doi: 10.1203/00006450-196703000-00006 6029810
69. Koestenberger M, Nagel B, Ravekes W, Avian A, Burmas A, Grangl G, et al. Reference values and calculation of z-scores of echocardiographic measurements of the normal pediatric right ventricle. Am J Cardiol. 2014;114: 1590–1598. doi: 10.1016/j.amjcard.2014.08.028 25248810
70. Majonga ED, Rehman AM, McHugh G, Mujuru HA, Nathoo K, Patel MS, et al. Echocardiographic reference ranges in older children and adolescents in sub-Saharan Africa. Int J Cardiol. 2017;248: 409–413. doi: 10.1016/j.ijcard.2017.06.109 28711335
71. Oran B, Bodur AS, Arslan D, Çımen D, Güvenç O. Normal M mode values in healthy Turkish children. Turk J Med Sci. 2014;44: 756–763. 25539541
72. Cain PA, Ahl R, Hedstrom E, Ugander M, Allansdotter-Johnsson A, Friberg P, et al. Age and gender specific normal values of left ventricular mass, volume and function for gradient echo magnetic resonance imaging: a cross sectional study. BMC Med Imaging. 2009;9: 2. doi: 10.1186/1471-2342-9-2 19159437
73. Friedberg MK. Chapter 9—Echocardiographic Quantitation of Ventricular Function. In: Jefferies JL, Chang AC, Rossano JW, Shaddy RE, Towbin JA, editors. Heart Failure in the Child and Young Adult. Boston: Academic Press; 2018. pp. 105–124. doi: 10.1016/B978-0-12-802393-8.00009-0
74. Muscogiuri G, Ciliberti P, Mastrodicasa D, Chinali M, Rinelli G, Santangelo TP, et al. Results of Late Gadolinium Enhancement in Children Affected by Dilated Cardiomyopathy. Front Pediatr. 2017;5: 13. doi: 10.3389/fped.2017.00013 28220144
75. Lee CK, Margossian R, Sleeper LA, Canter CE, Chen S, Tani LY, et al. Variability of M-mode versus two-dimensional echocardiography measurements in children with dilated cardiomyopathy. Pediatr Cardiol. 2014;35: 658–667. doi: 10.1007/s00246-013-0835-9 24265000
76. Bozkurt S, van de Vosse FN, Rutten MCM. Improving arterial pulsatility by feedback control of a continuous flow left ventricular assist device via in silico modeling. Int J Artif Organs. 2014;37: 773–785. doi: 10.5301/ijao.5000328 24970558
77. Granegger M, Masetti M, Laohasurayodhin R, Schloeglhofer T, Zimpfer D, Schima H, et al. Continuous Monitoring of Aortic Valve Opening in Rotary Blood Pump Patients. IEEE Trans Biomed Eng. 2016;63: 1201–1207. doi: 10.1109/TBME.2015.2489188 26461795
78. Mann Douglas L, Bristow Michael R. Mechanisms and Models in Heart Failure. Circulation. 2005;111: 2837–2849. doi: 10.1161/CIRCULATIONAHA.104.500546 15927992
79. Jackson G, Gibbs CR, Davies MK, Lip GYH. Pathophysiology. BMJ. 2000;320: 167–170. doi: 10.1136/bmj.320.7228.167 10634740
80. Su JA, Menteer J. Outcomes of Berlin Heart EXCOR® pediatric ventricular assist device support in patients with restrictive and hypertrophic cardiomyopathy. Pediatr Transplant. 2017;21. doi: 10.1111/petr.13048
81. Anwer LA, Poddi S, Tchantchaleishvili V, Daly RC, Toeg H, Strand JJ, et al. Left Ventricular Assist Devices: How Do We Define Success? ASAIO J Am Soc Artif Intern Organs 1992. 2019;65: 430–435. doi: 10.1097/MAT.0000000000000894 30312213
Článek vyšel v časopise
PLOS One
2019 Číslo 10
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
Nejčtenější v tomto čísle
- Correction: Low dose naltrexone: Effects on medication in rheumatoid and seropositive arthritis. A nationwide register-based controlled quasi-experimental before-after study
- Combining CDK4/6 inhibitors ribociclib and palbociclib with cytotoxic agents does not enhance cytotoxicity
- Risk factors associated with IgA vasculitis with nephritis (Henoch–Schönlein purpura nephritis) progressing to unfavorable outcomes: A meta-analysis
- Prevalence of pectus excavatum (PE), pectus carinatum (PC), tracheal hypoplasia, thoracic spine deformities and lateral heart displacement in thoracic radiographs of screw-tailed brachycephalic dogs