- Open Access
Closed-loop real-time simulation model of hemodynamics and oxygen transport in the cardiovascular system
BioMedical Engineering OnLine volume 12, Article number: 69 (2013)
Computer technology enables realistic simulation of cardiovascular physiology. The increasing number of clinical surgical and medical treatment options imposes a need for better understanding of patient-specific pathology and outcome prediction.
A distributed lumped parameter real-time closed-loop model with 26 vascular segments, cardiac modelling with time-varying elastance functions and gradually opening and closing valves, the pericardium, intrathoracic pressure, the atrial and ventricular septum, various pathological states and including oxygen transport has been developed.
Model output is pressure, volume, flow and oxygen saturation from every cardiac and vascular compartment. The model produces relevant clinical output and validation of quantitative data in normal physiology and qualitative directions in simulation of pathological states show good agreement with published data.
The results show that it is possible to build a clinically relevant real-time computer simulation model of the normal adult cardiovascular system. It is suggested that understanding qualitative interaction between physiological parameters in health and disease may be improved by using the model, although further model development and validation is needed for quantitative patient-specific outcome prediction.
The increasing number of surgical and medical treatment options in cardiovascular disease imposes a need for better understanding of patient-specific pathology and outcome prediction. Due to the complexity of the cardiovascular system, it is not always possible to understand the role of every cardiovascular parameter; additionally, prediction of the hemodynamic features can be even more challenging when several of those parameters change simultaneously. In many cases the use of simulation models can improve the understanding of both normal physiology and pathophysiology. Moreover, clinical training in simulation centres, where simulation models have been connected to mannequins responding to therapeutic measures, has proved to further increase diagnostic and therapeutic accuracy .
Other authors  have published complex hemodynamic models based on a reductionistic approach, where contractions are based on cellular properties including electrophysiology, ion channels, myofilament structure and fiber orientation. These models are relevant for exploring normal and disturbed contractile mechanisms as well as integrated organ function. The computational load in 3D models is however presently too large to allow clinical and educational use in studies of central hemodynamics. Simpler models based on cellular properties  tend to be less relevant on higher explanatory levels when studying integrated properties of the entire cardiovascular system, although being computationally less demanding.
Models focusing on particular aspects of hemodynamics such as mitral valve dynamics , remodeling of heart and vessels , individual drug effects  and effects of therapeutic measures such as positive end-expiratory pressure (PEEP)  offer good examples of how cardiovascular simulation models can be validated in specific situations and used clinically. It is however difficult to know to what extent these models can be used in clinical situations other than those described in the articles.
The aim of this paper is to present a closed-loop, real-time, lumped parameter cardiovascular simulation model, which includes the behavior of the four cardiac chambers, interatrial and interventricular septum, pericardium, intrathoracic pressure changes, heart valves, intracardiac shunts and vascular system. The current model is based on previous publications [8, 9]. The circulatory system interacts with pericardial as well as intrathoracic pressures as in one of the previous models  and several improvements have been made. The vascular system has been expanded and the circulatory system closed in order to reach more realistic steady state output when parameters in the model are changed. The systolic cardiac model is replaced in order to more realistically simulate cardiac work and valvular timing. Passive diastolic properties are now non-linear  resulting in more realistic changes of filling pressures with ventricular dilatation. Atrial septal interactions are included in addition to a modified ventricular septum . The valves are gradually opening and closing  and adjustment of valvular properties therefore allow more realistic simulation of both stenotic and regurgitant valves. Potential atrial septal defects, ventricular septal defects and persistent ductus arteriosus have been included. An oxygen transport model allowing studies of oxygen delivery and uptake including myocardial oxygen balance have been developed. The most unique feature of this model is however not any particular advanced modeling solutions in each sub-part, but rather the holistic approach giving a comprehensive overview of cardiovascular hemodynamics in health and disease. The possibility to change the model parameters one by one facilitate understanding in educational sessions, but real-time simulation and fast feed-back after parameter changes are also important features in a clinical decision support scenario in cardiac surgery and intensive care.
The model is presented from a normal cardiovascular physiological point of view, but a limited number of pathological scenarios will also be shown. A secondary aim is to show that the model generates valid output, which is consistent with earlier published clinical and experimental data, focusing on the behavior of the entire system rather than on single variables.
General overview of the model
The model (Figure 1) is based on the following compartments: the four cardiac chambers and corresponding valves, the pericardium, the systemic circulation consisting of the aortic root, ascending aorta, proximal aortic arch, right carotid/subclavian artery, distal aortic arch, left carotid/subclavian artery, descending aorta and the representation of a peripheral arterial system (See also electrical analogue in Additional file 1: Figure S1). The peripheral arterial system contains resistance arteries (arteriolae), capillaries and small veins emptying into the inferior caval vein, while the two carotid/subclavian arteries in a similar way end up in the superior caval vein. The pulmonary circulation is represented by the main pulmonary trunk, pulmonary resistance arteries (arteriolae), capillaries, small veins and pulmonary veins emptying into the left atrium. Two coronary vessels are also included. Input to the vascular model are physical dimensions [12, 13], tissue properties and number of parallel vessels in each compartment modified to fit the anatomical structure of the present vascular model; while pressure, flow, volume and oxygen saturation are simulated outputs and can be displayed. The cardiac, valvular, pericardial and vascular properties can be modified in detail by changing the parameters presented in the supplement (Additional file 2: Table S1-5). Model complexity is chosen to enable future simulations of central hemodynamics in major congenital and acquired circulatory pathologies and extra-corporeal circulatory support.
The heart chambers
The four cardiac chambers passive and active properties are modeled with separate time-varying elastance functions based on the “Double-Hill equation” as previously described [11, 14] including a slightly modified Starling mechanism (Figure 2) . The model includes internal chamber flow resistance and viscous wall properties of the heart, in accordance with previous models [9, 11]. Details are shown in the cardiac model section of the supplement.
By using a continuous function through the cardiac cycle manual setting of valvular timing  is avoided. Each cardiac chamber has its own elastance function of realistic shape (Additional file 3: Figure S3). Timing and duration of contractions in the cardiac chambers can be set independently, in order to create variable atrioventricular (AV) and interventricular (VV) intervals as well as realistic timing of valve opening and closure. Thereby the model offers additional opportunities to simulate arrhythmias as well as interventricular dyssynchronicity. As an example atrial fibrillation can be simulated by turning off atrial contractions and randomizing the time between ventricular contractions.
The electrocardiogram (ECG) (Additional file 4: Figure S4; Additional file 5: Figure S5) is simulated by summing the four cardiac elastance functions with a linear weight factor, based on the approximate distance from each cardiac chamber to the chosen surface ECG lead assuming an electromechanical delay of 60 ms. The ECG does not by itself influence the simulated physiology and is included as an output of the model only to add clinical realism during educational sessions, therefore the detailed algorithms are not described.
The heart valves
Valve pressure gradients are composed of a Bernoulli resistance and an inertial term (Additional file 2: Equation S4 and Figure S6 of Additional file 6) as earlier described [8, 9]. Minimum and maximum valvular areas can be set to simulate valvular stenosis and regurgitation. The valve areas change gradually (Additional file 4: Figure S4) as a function of pressure gradients and valve inertial constants  as is shown in the supplement (Additional file 2: Table S4). Blood flow inertance L is determined by density ρ, the variable inflow length l and the valve area A (Equation 1). Inflow length l is set to a value identical to the instantaneous diameter of the valve in contrast to the constant value in Mynard et al..
The atrial and ventricular septum
The interventricular septal displacement is simulated using pressure transmission through the muscular septum as previously described [9, 15] (Equation 2). As a further development the ventricular septal stiffness E sv is modeled in order to increase in proportion to the left ventricular systolic elastance e lv (t) (Equation 3), thereby stiffening during contraction, in accordance with known physiological characteristics. Additional file 5: Figure S5 shows the ventricular septal shift at different stiffness.
Atrial septal interaction is modeled in a similar way. The septal interaction does not give any significant contribution during normal physiology in apnea, but is important when simulating interventricular dyssynchrony, variations in intrathoracic pressures including artificial ventilation and changes in ventricular loading conditions.
In order to resemble pericardial function in health and disease a user-defined exponential function relating intra-pericardial pressure to total heart volume is adopted from the literature  (Additional file 7: Figure S7). As a further development of the previous model the minimal pericardial pressure is allowed to be a negative value in agreement with what is found experimentally when measured during intrathoracic pressure changes or hypovolemia . Further details can be found in the supplement.
Intracardiac and arteriovenous shunts
Atrial septal defects, ventricular septal defects and persistent ductus arteriosus flows are calculated according to Equation 4, with instantaneous flow Q proportional to the area A of the defect, the square-root of the pressure difference ∆P and the Gorlin empirical constant C assumed to be 1 . Acceleration g due to gravity is set to 980 cm/s2.
The systemic circulation
Transmural pressure p in each vascular compartment is related to the variable segmental vascular volume v according to an exponential relation (Additional file 2: Equation S7, Additional file 8: Figure S12)  creating an increase in stiffness with progressive distension of the vascular wall.
Resistance R 0 , inertia I 0 , and elastance E 0 at normal mean pressure P 0 are calculated from published and estimated vascular properties [12, 13], and relations  described in Equations 5–7. These vascular properties are all approximately inversely proportional to n, which is the number of parallel vessels, lumped into each compartment . The parameter η is blood viscosity, l, r 0 and h is the length, radius and thickness of a single vessel respectively and ρ is blood density.
The same value of Young’s modulus Y inc is used in every vascular segment, implying that vascular stiffness properties are solely explained by this constant and vascular dimensions. The actual volume-dependent segmental elastances, resistances and inertias are updated in each calculation step based on the volume and radius assuming constant vessel length and thickness. Details can be found in the blood vessels section of the supplement (Additional file 2: Table S5-7).
The resistance Ω, in series with the capacitor situated in parallel to the blood flow in every vascular segment, corresponds to the viscous property of the vascular wall, dampening the flow pulse. The parameter value assigned to this component is calculated as in Equation 8 with addition of a scaling factor λ as compared to the characteristic impedance in a classical 3-component Windkessel model [19, 20].
The same formula was used to calculate Ω in every vascular segment as previously suggested , although the name characteristic impedance is less well suited in a distributed model. A scaling factor λ between 0.0 and 1.0 is used to tune the vascular damping properties. The effects of changing the scaling factor λ can be seen in Figure 3.
The systemic capillaries and veins are described by the same principal vascular model as the arteries although inertial effects are of less importance due to lower pulsatility  and a larger cross-sectional area. Every single vessel has high resistance and stiffness; however the large amounts of vessels in parallel lumped into each compartment result in both a low resistance (Equation 5) and a low elastance (Equation 7), corresponding to high compliance and being related to their physiological volume reservoir function. The distribution of the blood volume between arteries, capillaries and veins in the systemic and pulmonary circulation in the model is shown in Additional file 2: Table S7.
The pulmonary circulation
The pulmonary artery, resistance vessels, capillaries and veins are modeled as in the systemic circulation. An adjustable pulmonary shunt mimics the circulation without gas exchange due to anatomical anomalies or ventilation-perfusion disturbances as in atelectatic lung tissue.
The coronary circulation
Coronary circulation is simulated with a left and right coronary artery emptying in the right atrium as described in the supplement (See also Additional file 9: Figure S8). Myocardial oxygen consumption is calculated according to the pressure-volume area concept of Suga et al. on a beat-by-beat basis and estimations of myocardial oxygen balance can therefore be calculated based on both supply and demand.
Blood oxygen content in every compartment is calculated as a constant multiplied by the hemoglobin (Hb ), the oxygen saturation (Sat ) and the volume of the compartment (Volume ) as in clinical routine and most catheterization laboratories  (Equation 9).
Oxygen saturation within each compartment is considered homogenous and oxygen transport between compartments is proportional to blood flow, hemoglobin level and saturation. Calculations concerning blood mixing are simplified as effects of blood jets, incomplete mixing and physically dissolved oxygen not are taken into account.
A baroreceptor reflex can be activated in the model that affects heart rate, cardiac contractility (maximum elastance) and arterial vascular resistance as in Sun et al. and described in the baroreceptor section of the supplement.
The simulator was compiled as a stand-alone software developed in Visual Basic.NET 2010 (Microsoft Corporation, Redmond, WA, USA) and .NET Framework 4.0 (Microsoft Corporation, Redmond, WA, USA). A time step of 0.25 ms was used to assure stability. The implemented 62 differential equations were solved with an implicit Euler numerical method . Mean values were calculated as running arithmetical means to simplify calculations and conserve memory resources. The program version used was Aplysia CorVascSim 184.108.40.206 (Aplysia Medical AB, Stockholm, Sweden). The model can be run in real-time on a standard personal computer (example: Intel processor Core i3 3.06 GHz, RAM 3 Gb, Graphics ATI Radeon HD4550, Windows 7 32/64-bit).
Parameter settings in simulations of normal physiology, systolic heart failure, diastolic heart failure, aortic stenosis, aortic regurgitation, a Valsalva maneuver, exercise and progressive arteriosclerosis
All hemodynamic data were presented with zero intrathoracic pressure in order to resemble end-expiratory values. With the basic set of parameters presented in the supplement an adult man weighing 70 kg is described with a blood volume of 5600 ml (80 ml/kg), blood viscosity of 0.00024 mmHg · s and blood density of 1.060 g/cm3. Systolic left ventricular heart failure was simulated by decreasing the maximum elastance from 2.8 to 1.0 mmHg/ml. Diastolic left ventricular heart failure was simulated by increasing the basal passive diastolic elastance E min from 0.05 to 0.12 mmHg/ml. Aortic stenosis was simulated by decreasing the open aortic valve area from 5.0 to 0.7 cm2. Aortic regurgitation was simulated by increasing the closed aortic valve area from 0.0 to 0.2 cm2. Simulation of a Valsalva maneuver was accomplished by increasing intrathoracic pressure from 0 to 10 mmHg for approximately 20 seconds with and without baroreceptors activated. All other parameters and settings were identical.
The stepwise simulation of exercise was accomplished by increasing heart rate from 72 to 144 min-1, increased left and right ventricular contractility by 50% and increased systemic and pulmonary resistance arterial diameters by 50%. The cardiac elastance parameters α 1 and α 2 were scaled with the square-root of the factor changing the heart rate to preserve a realistic relation between systolic and diastolic times. The vascular damping scaling factor λ was set to 0.5 in all vascular segments and presented simulations. The oxygen saturation in the pulmonary capillaries was set to 99.4%, corresponding to a normal alveolar pO2 of 13.3 kPa and the pulmonary shunt was set to 10% of cardiac output in all presented simulations. Hemoglobin level was set to 140 g/l and systemic oxygen consumption to 250 ml/min in addition to myocardial oxygen extraction. A Young’s modulus Y inc of 3000 mmHg (= 0.40 MPa) similar to the value in Wang et al. was used to calculate the vascular elastance E 0 at the vessel-specific normal mean pressure P 0 (see Additional file 2: Table S5) in all simulations above. A progressive arteriosclerotic process was simulated by changing vascular stiffness in all vascular segments through variations in Young’s modulus Y inc between 2000 and 6000 mmHg. Note that comparisons with real clinical cases will not be performed due to unknown patient-specific secondary adaptations.
The relation between changes in input parameters x and model output variables y was studied with a sensitivity analysis during simulation of normal physiology (Additional file 10). Input parameters were increased 10% one by one. Seven model output variables [left ventricular end-systolic pressure (LVESP), right ventricular end-systolic pressure (RVESP), left atrial pressure (LAP), right atrial pressure (RAP), left ventricular stroke work (LVSW), right ventricular stroke work (RVSW) and cardiac output (CO)] focusing on cardiac function were evaluated after 60 seconds assuring a hemodynamic steady-state. Sensitivity S was calculated according to Equation 10.
Overall, the direction and magnitude of the changes in output variables and curve shapes were within a plausible clinical range in the simulated cases. The specific data for each case is presented below.
Running the model with basal settings representative of a 30 years old healthy man resulted in cardiac output 5.09 l/min, systemic arterial blood pressure 112/61(79) mmHg, pulmonary arterial pressure 24/8(11) mmHg, mean right atrial pressure 4 mmHg and mean left atrial pressure 5 mmHg. The model output of pressures (Figure 4 and Additional file 11: Figure S9), volume-flow rates (Figure 5 and Additional file 12: Figure S10) and timing within the cardiac chambers and the systemic and pulmonary circulatory parts created by the model resembled in pattern and magnitude the ones found in clinical practise (Table 1).
The shapes and areas (stroke work) of the ventricular pressure-volume loops (Additional file 13: Figure S11a) as well as right and left ventricular ejection fraction were within normal limits  (Table 1). The atrial loops (Additional file 13: Figure S11b) were similar to the scarce published data .
Systolic heart failure
The simulation of systolic heart failure resulted in decreased left ventricular ejection fraction, stroke work, stroke volume and systolic pressure, while left ventricular end-diastolic volume and filling pressures increased (Figure 6 and Table 2).
Diastolic heart failure
The simulation of diastolic heart failure showed preserved left ventricular ejection fraction, but decreased stroke work, end-diastolic volume, stroke volume and systolic pressure, while left ventricular filling pressures increased (Figure 7 and Table 2).
The simulation of aortic stenosis resulted in decreased left ventricular ejection fraction, stroke volume and systemic blood pressure, although left intraventricular pressure and stroke work increased (pressure-load). Left ventricular end-diastolic volume and filling pressures increased (Figure 8 and Table 2).
Simulated data of aortic regurgitation showed apparent left ventricular ejection fraction and stroke volume increase, systemic blood pressure decrease although stroke work increased (volume-load). Left ventricular end-diastolic volume and filling pressures increased (Figure 8 and Table 2).
The simulation of a Valsalva maneuver without and with activated baroreceptor reflex is shown in Figure 9. The sudden decrease in afterload with preserved preload when increasing intrathoracic pressure  is the explanation for the initial rise in arterial pressure as indicated in the figure. The opposite effect was seen immediately after decreasing intrathoracic pressure. Arterial blood pressure was as expected better preserved when the baroreceptor reflex was activated due to increase in heart rate, ventricular contractility and systemic arterial resistance.
The stepwise simulation of exercise is shown in Figure 10. Starting from a resting state, heart rate and cardiac contractility is increased followed by reductions in systemic and pulmonary vascular resistance resulting in a total cardiac output above 10 l/min.
The presence of arteriosclerosis, simulated by increased vascular stiffness, resulted in increased systolic aortic pressure and decreased diastolic aortic pressure. A decrease in systolic ejection time during this process is indicated by the premature arrival of the reflected pressure wave changing the shape of the pressure curves in a characteristic way (Figure 11).
Blood volume was the single most important parameter in the sensitivity analysis (Additional file 10). The pericardial volume constant v_pc0 was second most important. Pulmonary arteriolae radius (radius_pur) was more important than systemic arteriolae radius (radius_r). Right ventricular passive diastolic stiffness (e0_rv) and systolic contractility (emax_rv) was more important than their left sided equivalents (e0_lv and emax_lv, respectively). Parameters with less than 10% sensitivity are not shown.
This work presents a closed-loop real-time lumped parameter simulation model of the cardiovascular system where normal physiology and oxygen transport as well as systolic heart failure, diastolic heart failure, valve failure, intracardiac shunts and pulmonary failure can be simulated separately or in combination.
The validation and appreciation of clinical relevance of the model is based on empirically relevant changes in simulated clinical parameters when being compared with measurements made by clinically available methods such as intravascular pressure catheters and ultrasound. It is not possible to validate the large number of parameters unambiguously with available clinical and/or experimental data, but when taking the complexity of the system into account and the variety of disease states that can be simulated, the clinical relevance becomes a powerful validation and hypothesis creating tool. Another approach to challenge the adequacy of a cardiovascular simulation model is to study if model responses are adequate when evaluating model output in new scenarios with pathological states not previously tested during construction and tuning of the model. We have chosen to present model output illuminating a range of normal and pathophysiological states, illustrating the wide scope of the model. The normal resting state, the Valsalva maneuver and physiological changes during exercise are parts of normal physiology that should be handled correctly by a multi-purpose cardiovascular simulation model. The pathophysiological states chosen represent important clinical problems with various backgrounds. Notable is that very few parameter changes are needed to simulate the basic features of these states. Further parameter changes are needed to simulate secondary compensatory adaptations such as e.g. ventricular hypertrophy in aortic stenosis. This can be automatized in a simulation model , but manual stepwise adaptation often enables a better understanding of the relation between the primary problem and its consequences. The lack of secondary changes such as blood volume increase and systemic vasoconstriction usually seen in clinical cases with low cardiac output can also explain the small increase in filling pressures seen in simulation of heart failure and aortic valve failure in the present work.
The combined effect of elastic recoil in large elastic arteries and reflection of pressure waves in the peripheral vasculature is of importance in coronary physiology as well as in determining left ventricular afterload. A distributed lumped parameter model as presented in this article is usually not considered adequate for describing wave transmission and reflection in the vascular tree, although the use of an inductor/inertia in each vascular compartment and a high resistance in peripheral components create pressure and flow curves with realistic amplitudes and time delays resembling effects of reflected waves (Figures 4, 5 and 11). The advantage with our approach is less demanding calculations enabling real-time simulation and also the possibility to connect artificial circulatory support to any node in the model as will be explored in future works.
Aortic characteristic impedance is originally defined as a single purely resistive component describing the properties of the entire vascular tree in the high frequency domain, when seen from the aortic valve . Placing this resistor in line with the blood flow, as is usually done in the 3-compartment Windkessel model, gives adequate results at high frequencies, but obscures the meaning of the more basic concept "total peripheral resistance" in the low frequency domain . When placing the resistor in line with the capacitor in parallel to the blood flow in each vascular compartment as is done in the present model (Additional file 14: Figure S2), this component still fulfills the role of the characteristic impedance in the high-frequency domain without affecting the total peripheral resistance. This is the reason why the characteristic impedance is calculated according to Equation 8. Furthermore it has a physical meaning representing the damping effect imposed by the viscous properties of the vascular wall (non-linear Maxwell viscoelastance model), although the detailed relation to tissue properties remains to be elucidated.
The algorithms for scaling of the cardiac elastance functions with changing heart rate seem efficient up to a frequency of 160 min-1. Above this heart rate cardiac output decreases suggesting either suboptimal diastolic heart function (in the model) or suboptimal scaling of the elastance functions. It is also possible that atrioventricular (AV)-plane motion, not included in the present model, is necessary to improve ventricular filling during tachycardia in real life physiology.
Sensitivity analysis in a non-linear model should ideally be performed for every possible set of input parameters. With many input parameters it is neither feasible to perform the calculations nor to interpret them in every possible model state. The choice of relevant parameters therefore must rely on physical and physiological model relations; however the relevance of the model is supported by the performed sensitivity analysis in normal physiology, where parameters determining cardiac filling, afterload and ventricular systolic performance are most important for cardiac flow and pressure generation (Additional file 10).
Limitations and future development of the current model
The scientific, clinical and educational value of the model is the possibility to study combinations of physiological and pathophysiological cardiac and vascular properties as well as the effects of potential modifications or treatments. This holistic approach enables hypothesis generation regarding the direction of changes and qualitative relations between parameters in a wide variety of normal and pathological states suggesting possible usefulness as a clinical decision support tool. Detailed quantitative conclusions are however not justified without further validating studies. It must also be kept in mind that differences in the magnitude of changes between model output and patients will persist due to model simplifications and non-modelled parameters. Patient related predictions or hypotheses created by any simulation model have to be weighed against other sources of knowledge such as clinical studies and experience.
The mathematical functions creating the time-varying elastance functions are not based directly on cellular or mechanical properties of the myocardium, but should rather be seen as a convenient and well-established way to mimic a realistic physiological ventricular filling and contraction pattern during variable loading conditions in a computer model. The chosen functions are scalable with changes in cardiac mechanical properties and heart rate.
Future improvements of the cardiac model aim at adjustment to known physiological relations, including the AV-plane movement contributions to the cardiac pumping mechanism  and scaling of the model to neonatal and pediatric size. A specific aim is also to develop and validate simulated hemodynamic responses in patients treated with extra-corporeal circulation during intensive care, such as left ventricular unloading during support with mechanical left ventricular assist devices and systemic oxygen delivery and pulmonary hypertension in congenital cardiopulmonary anomalies. Further development of the model will also aim at implementing a semi-automatic adaptation to available patient-specific clinical data and from this enable prediction of the hemodynamic results from possible surgical and pharmacological interventions.
During extreme conditions such as cardiac arrest, extreme hypovolemia or volume overload, when pressures are outside a normal range an improved model including simulation of vascular collapse and overdistension needs to be developed. Moreover, due to the lack of true 3D properties some pathology, e.g. myocardial ischemia and infarctions with regional cardiac dysfunction, cannot be represented in detail. The same problem is inherent in many cardiac arrhythmias involving aberrant activation of the myocardium.
The results of the present publication show that it is possible to build a real-time computer simulation model illustrating the cardiovascular system, including oxygen transport, as evidenced by the base-line values almost entirely within reference limits of normal physiology and by relevant direction of changes in clinical parameters in a variety of pathophysiological states. The understanding of complex hemodynamic states may be further improved when known model parameter changes are studied and simultaneously compared in detail by the user with clinical data.
Blood flow inertance
Valve inflow length
Ventricular septal stiffness
Ventricular septal stiffness constant
Left ventricular systolic elastance
Left ventricular pressure
Left ventricular volume
Right ventricular pressure
Gorlin empirical constant
Resistance at normal mean pressure
Inertia at normal mean pressure
Elastance at normal mean pressure
Normal mean pressure
Number of parallel vessels in vascular compartment
Vessel radius at normal mean pressure
Vascular wall resistance
Vascular wall resistance scaling factor
Model compartment volume
Basal cardiac chamber passive stiffness
Left ventricular end-systolic pressure
Right ventricular end-systolic pressure
Left atrial pressure
Right atrial pressure
Left ventricular stroke work
Right ventricular stroke work
Pericardial volume constant
Pulmonary arteriolae radius
Systemic arteriolae radius
Right ventricular passive diastolic stiffness
Right ventricular systolic contractility
Duration of right ventricular relaxation
Systemic arterial radius
Left ventricular passive diastolic stiffness
Normal mean pressure in small systemic veins
Pulmonary capillary radius
Left ventricular systolic contractility
Right carotid arterioli radius
Left carotid arterioli radius
Pulmonary small vein radius
Pulmonary arterioli number
Pulmonary arterioli length
Systemic arterioli length
Right ventricular stiffness constant
Systemic arterioli number
Maximum pressure derivative
Minimum pressure derivative
Ventricular relaxation constant Tau
Left ventricular end-diastolic volume
- LA volume:
Left atrial volume
Tei index = Myocardial performance index
- RA volume min --- max:
Right atrial volume
Systolic arterial pressure
Left ventricular ejection fraction
Left ventricular stroke volume
Left ventricular end-diastolic volume
Pulmonary arterial pressure
Okuda Y, Bryson EO, DeMaria S, Jacobson L, Quinones J, Shen B, Levine AI: The utility of simulation in medical education: what is the evidence? Mt Sinai J Med 2009, 76(4):330–43. 10.1002/msj.20127
Hunter PJ, Pullan AJ, Smaill BH: Modeling total heart function. Annu Rev Biomed Eng 2003, 5: 147–77. 10.1146/annurev.bioeng.5.040202.121537
Shim EB, Jun HM, Leem CH, Matusuoka S, Noma A: A new integrated method for analyzing heart mechanics using a cell-hemodynamics-autonomic nerve control coupled model of the cardiovascular system. Prog Biophys Mol Biol 2008, 96(1–3):44–59.
Paeme S, Moorhead KT, Chase JG, Lambermont B, Kolh P, D'orio V, Pierard L, Moonen M, Lancellotti P, Dauby PC, Desaive T: Mathematical multi-scale model of the cardiovascular system including mitral valve dynamics. Application to ischemic mitral insufficiency. Biomed Eng Online 2011, 10: 86. 10.1186/1475-925X-10-86
Arts T, Delhaas T, Bovendeerd P, Verbeek X, Prinzen FW: Adaptation to mechanical load determines shape and properties of heart and circulation: the CircAdapt model. Am J Physiol Heart Circ Physiol 2005, 288(4):H1943–54.
Chase JG, Starfinger C, Hann CE, Revie JA, Stevenson D, Shaw GM, Desaive T: Model-based prediction of the patient-specific response to adrenaline. Open Med Inform J 2010, 4: 149–63. 10.2174/1874431101004010149
Starfinger C, Chase JG, Hann CE, Shaw GM, Lambert P, Smith BW, Sloth E, Larsson A, Andreassen S, Rees S: Model-based identification of PEEP titrations during different volemic levels. Comput Methods Programs Biomed 2008, 91(2):135–44. 10.1016/j.cmpb.2008.03.005
Sun Y, Sjöberg BJ, Ask P, Loyd D, Wranne B: Mathematical model that characterizes transmitral and pulmonary venous flow velocity patterns. Am J Physiol 1995, 268(1 PT 2):H476–89.
Sun Y, Beshara M, Lucariello RJ, Chiaramida SA: A comprehensive model for right-left heart interaction under the influence of pericardium and baroreflex. Am J Physiol 1997, 272(3 PT 2):H1499–515.
Chung DC, Niranjan SC, Clark JW, Bidani A, Johnston WE, Zwischenberger JB, Traber DL: A dynamic model of ventricular interaction and pericardial influence. Am J Physiol 1997, 272(6 PT 2):H2942–62.
Mynard JP, Davidson MR, Penny DJ, Smolich JJ: A simple versatile model of valve dynamics for use in lumped parameter and one-dimensional cardiovascular models. International Journal for Numerical Methods in Biomedical Engineering 2011, 28(6–7):626–41.
Wang JJ, Parker KH: Wave propagation in a model of the arterial circulation. J Biomech 2004, 37(4):457–70. 10.1016/j.jbiomech.2003.09.007
Milnor WR: Hemodynamics. 2nd edition. Baltimore: William and Wilkins; 1989.
Stergiopulos N, Meister JJ, Westerhof N: Determinants of stroke volume and systolic and diastolic aortic pressure. Am J Physiol 1996, 270(6 PT 2):H2050–9.
Maughan WL, Sunagawa K, Sagawa K: Ventricular systolic interdependence: volume elastance model in isolated canine hearts. Am J Physiol 1987, 253(6 PT 2):H1381–90.
Applegate RJ, Santamore WP, Klopfenstein HS, Little WC: External pressure of undisturbed left ventricle. Am J Physiol 1990, 258(4 PT 2):H1079–86.
Gorlin R, Gorlin SG: Hydraulic formula for calculation of the area of the stenotic mitral valve, other cardiac valves, and central circulatory shunts. Am Heart J 1951, 41(1):1–29. 10.1016/0002-8703(51)90002-6
Westerhof N, Stergiopulos N, Noble M: Snapshots of Hemodynamics. An aid for clinical research and graduate education. Boston: Springer; 2005.
Arts T, Reesink K, Kroon W, Delhaas T: Simulation of adaptation of blood vessel geometry to flow and pressure: Implications for arterio-venous impedance. Mech Res Commun 2012, 42: 15–21.
Murgo JP, Westerhof N, Giolma JP, Altobelli SA: Aortic input impedance in normal man: relationship to pressure wave forms. Circulation 1980, 62(1):105–16. 10.1161/01.CIR.62.1.105
Lee JJ, Tyml K, Menkis AH, Novick RJ, Mckenzie FN: Evaluation of pulsatile and nonpulsatile flow in capillaries of goat skeletal muscle using intravital microscopy. Microvasc Res 1994, 48(3):316–27. 10.1006/mvre.1994.1058
Suga H: Total mechanical energy of a ventricle model and cardiac oxygen consumption. Am J Physiol 1979, 236(3):H498–505.
Hufner G: Neue versuche zur bestimmung der sauerstoffcapacitat der blutyfarbstoffs. Arch Physiol 1902, 17(22):130–76.
Butcher JC Numerical methods for ordinary differential equations. In Chapter 23: Runge–Kutta Methods. Chichester, UK: John Wiley & Sons; 2003:93–104.
Lentner C Geigy Scientific Table. In Volume 5: Heart and Circulation.8th revised and enlarged edition. 5th edition. West Caldwell, New Jersey: CIBA-GEIGY; 1990.
Dernellis JM, Stefanadis CI, Zacharoulis AA, Toutouzas PK: Left atrial mechanical adaptation to long-standing hemodynamic loads based on pressure-volume relations. Am J Cardiol 1998, 81(9):1138–43. 10.1016/S0002-9149(98)00134-9
Harrigan PWJ, Pinsky MR: Heart-lung interactions. Part II: effects of intrathoracic pressure. Intensive Care 2001, Part II: 99–106. Summer
Lundbäck S: Cardiac pumping and function of the ventricular septum. Acta Physiol Scand Suppl 1986, 550: 1–101.
MBs research during the years 2013–2015 is founded by The Swedish Research Council Grant 2012–2800. BJS have a strategic professor position in Medical Technology at Karolinska Institutet.
Michael Broomé is the founder and owner of the company Aplysia Medical AB developing the simulation software Aplysia CorVascSim. There are no other conflicts of interest.
MB constructed the model, performed programming, sensitivity analysis and drafted the manuscript. EM and AB participated in model construction and adaptation of the model to engineering standards as well as in manuscript drafting. BJS participated in model development in relation to previously published models and drafting of the manuscript. MB, BF and BJS all participated in evaluation of the clinical relevance of the model as being clinically active medical doctors. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1.: Electrical analogue sketch showing both the heart and the vascular system. The dark yellow area is the pericardium containing the cardiac chambers and coronary vessels. The light yellow area is the intrathoracic space containing the pericardium, the pulmonary circulation and the thoracic aorta. (TIFF 254 KB)
Additional file 3: Figure S3: Left (black thick) and right (gray thick) ventricular time-varying elastance functions during two heart cycles. Atrial elastance functions are shown with thin lines. The amplitude is closely related to contractile function while the volume-dependent basal level describes passive stiffness. (TIFF 143 KB)
Additional file 4: Figure S4: The ECG and valvular areas simulated during two normal heart beats. ECG (light blue) and area-changes during the heart cycle for aortic (red), pulmonary (yellow), mitral (brown), tricuspid (blue) valves are shown. (TIFF 251 KB)
Additional file 5: Figure S5: Modeled volume changes (black) due to pressure dependent ventricular septal volume shift from left to the right ventricle during two heart cycles. Septal stiffness values of 40, 30, 20 and 10 mmHg/ml are shown. Lower septal stiffness increases septal shift. ECG (blue) is shown as reference. (TIFF 151 KB)
Additional file 6: Figure S6: Simulated pressure gradients over the mitral valve during diastole. The total mitral pressure gradient (thick black) is decomposed into a dominating Bernoulli gradient (grey) and an inertial component (dotted black). (TIFF 204 KB)
Additional file 7: Figure S7: Pericardial pressure-volume relations. The black curve illustrates the relation used in the present model with the range used by a single heart beat in orange. Pericardial pressure is slightly negative. (TIFF 96 KB)
Additional file 8: Figure S12: Pressure-volume relations in arterial vascular segments illustrating non-linear stiffness in the actual working pressure range of each vascular segment in accordance with Equation 10. (TIFF 68 KB)
Additional file 9: Figure S8: Blood flow in left (black) and right (gray) coronary artery during two heart cycles. Left coronary artery blood flow decreases during systole due to vascular compression caused by high left ventricular wall stress. (TIFF 134 KB)
Additional file 11: Figure S9: Pressure changes simulated during two heart cycles in the aortic root (black), ascending aorta, proximal and distal aortic arch, descending aorta and a peripheral artery (light gray). The segments between the aortic root and the peripheral artery are shown in different shades of gray. (TIFF 107 KB)
Additional file 12: Figure S10: Flow pulsatility in pulmonary capillaries (pink) is larger than in the systemic capillaries (green). Pulmonary valve flow (orange) and tricuspid valve flow (blue) resembling left-sided flows are also shown. (TIFF 221 KB)
Additional file 13: Figure S11: a. Intracardiac pressure-volume loops from the left ventricle (black) and right ventricle (gray). b. Intracardiac pressure-volume loops from the left atrium (black) and right atrium (gray). (TIFF 81 KB)
Additional file 14: Figure S2: Electrical analogue sketches showing details of the heart chambers, valves and vascular compartments. The arrows show directions of blood flow. e(t); time-varying elastance representing systolic and diastolic chamber properties, Rwall; a resistance representing viscous chamber wall properties, Routflow; a linear chamber outflow resistance, B; a non-linear Bernoulli valve resistance, Ivalve; an inductance representing flow inertia,, R; a linear resistance, I; an inductance representing flow inertia, C; a non-linear capacitor representing vascular elasticity, Ω; a non-linear resistance term representing viscous vascular wall properties. (TIFF 70 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Broomé, M., Maksuti, E., Bjällmark, A. et al. Closed-loop real-time simulation model of hemodynamics and oxygen transport in the cardiovascular system. BioMed Eng OnLine 12, 69 (2013). https://doi.org/10.1186/1475-925X-12-69
- Cardiovascular simulation
- Time-varying elastance functions
- Lumped parameter model
- Valve model
- Oxygen transport model
- Computer simulation