Beat-to-beat estimation of the continuous left and right cardiac elastance from metrics commonly available in clinical settings
© Stevenson et al.; licensee BioMed Central Ltd. 2012
Received: 23 December 2011
Accepted: 30 July 2012
Published: 21 September 2012
Functional time-varying cardiac elastances (FTVE) contain a rich amount of information about the specific cardiac state of a patient. However, a FTVE waveform is very invasive to directly measure, and is thus currently not used in clinical practice. This paper presents a method for the estimation of a patient specific FTVE, using only metrics that are currently available in a clinical setting.
Correlations are defined between invasively measured FTVE waveforms and the aortic and pulmonary artery pressures from 2 cohorts of porcine subjects, 1 induced with pulmonary embolism, the other with septic shock. These correlations are then used to estimate the FTVE waveform based on the individual aortic and pulmonary artery pressure waveforms, using the “other” dysfunction’s correlations as a cross validation.
The cross validation resulted in 1.26% and 2.51% median errors for the left and right FTVE respectively on pulmonary embolism, while the septic shock cohort had 2.54% and 2.90% median errors.
The presented method accurately and reliably estimated a patient specific FTVE, with no added risk to the patient. The cross validation shows that the method is not dependent on dysfunction and thus has the potential for generalisation beyond pulmonary embolism and septic shock.
Cardiac disturbances are difficult to diagnose in critical care, which can lead to poor management [1, 2]. Inadequate diagnosis can be common, and plays a significant role in increased length of stay, cost and mortality [3–5]. Currently, measurements are obtained via catheters placed around the heart. However, this limited set of data can severely restrict clinical diagnostic capability, and, as a result, these catheters are not necessarily associated with improved outcomes [6–8].
Overall, much of the data currently available to clinicians in an intensive care unit (ICU) that could have significant clinical value is under utilised. For example, acute cardiovascular dysfunction, like pulmonary embolism and septic shock, severely alter cardiovascular system hemodynamics around the heart. These changes can be seen directly by catheter measurements as changes in pressure and flow that reflect dynamic changes in the balance of pre-load and afterload, resulting in an altered cardiac energetic state [9, 10]. Detailed cardiac energetics are too invasive to measure in an ICU setting. However, if the relevant energetics could be captured using data from a nearby catheter, the clinical potential of such measurements could be realised. To date, no such method achieves this aim.
All models are approximations to capture observed physics, and thus all models offer advantages in understanding these observations as well as limitations of any necessary assumptions made. This paper presents a model for estimating the functional time-varying elastance from observed ventricular pressure-volume behaviour, to be used in the modelling of broader cardiovascular system [11–15]. Thus, this is the elastance required to deliver the pressure-volume measurements observed, which implicitly contain the pre-load and afterload information. This is in contrast to the “true” elastance which is a physical property of the heart muscle alone and can be tested only ex-vivo. To make this distinction clear, the elastance presented in this paper is referred to as the “functional elastance”.
where, V d is assumed to be equal to V 0 for simplicity, V 0 is the intercept of the end-systolic pressure-volume relation with the volume axis , P v (t) is the ventricle pressure and V v (t) is the ventricle volume. The waveform e(t) is typically normalised to a value of 1.0 . This model of load-dependent functional time-varying elastance should not be confused with the measurement of end-systolic elastance, which is regarded as load insensitive within physiological range. However, the load-independence of end-systolic elastance is an approximation made to that model, as the true end-systolic elastance is dependent on loading condition, if only weakly within physiological range . Load conditions are diagnostic and thus the model for FTVE presented monitors these changes as the time varying ability of the heart to pump blood. It thus provides a measure of heart function and energetics [18, 20, 21].
There have been several attempts to estimate FTVE [20, 22–25]. Most studies present a method using the FTVE to estimate a specific parameter, most commonly end-systolic elastance [20, 24, 25] and ejection fraction . However, their validation is based on these metrics, not on the resulting FTVE waveform, which contains unique clinical information. This study is focused not on the absolute values of the elastance (such as end-systolic elastance). Rather, the focus is on the shape, and change of shape of the time-varying waveform within a specific patient as dysfunction occurs and the patient’s state changes. This intra-patient variability is reflected in the (relative) shape.
It is unclear how much specific clinical information can be obtained from the FTVE waveform, other than the highly sought-after end-systolic elastance , and a measure of cardiac work (the area under the FTVE waveform is analogous to work done by the ventricle). The absolute value of the work done is lost with the normalisation of the waveform, but the relative changes to work over time remain. This holds true for a constant inotropic state, as the end-systolic elastance will remain the same. For a varying end-systolic elastance, the relative change in cardiac power will still be visible to the FTVE through the altered shape from an increased heart rate.
Beyond cardiac work a number of FTVE features have been highly correlated to clinically relevant parameters in a previously developed circulatory model . Equally, it contains similar information to that of the pressure-volume (P-V) loops, which are known to contain information on cardiac function  including cardiac work [29, 30], contractility [18, 31], O2 consumption [30, 32], and all the states of filling, contraction ejection and relaxation . Overall, cardiac elastances reflect cardiac state, cardiac output or blood volume, and net pre-load and afterload, all of which change with different cardiac dysfunction. Hence, the ability to directly measure FTVE, which does not currently exist without significantly invasive added testing, should yield clinically useful insight and diagnostics.
This research is unique in that the end goal is to produce the FTVE function in its own right, validating the FTVE waveform on its own accuracy compared to the invasively and directly measured waveform. This is done through the use of only commonly available metrics in an ICU setting, and thus, by non-invasive we mean, no more invasive than traditional care.
Cardiac elastance cannot be directly measured without invasive procedures. However, the information contained within the cardiac elastance waveform can be seen in other parts of the closely intertwined cardiovascular system. Therefore, with the right knowledge and transformations, that information can be collected and translated into an estimation of a metric (cardiac elastance) that would not normally be possible to obtain.
The best accessible source of indirect information about cardiac elastance is directly down stream of the ventricles at the aorta and the pulmonary arteries. They are the first to see any results of a change in elastance during systole. Thus, the aortic and pulmonary artery pressure waveforms (P ao and P pa , respectively) are key sources of information, which contain surrounding load conditions, to reconstruct cardiac elastances. It can be assumed, and seen from the data presented in this paper, that the shape of the functional time-varying cardiac elastance is load dependent.
Figure 1 shows the overall method of this paper. Known points on the aortic pressure waveform are correlated with points on the left ventricle FTVE (described in Section “Cardiac elastance correlations”), giving equations that enable estimations of these left ventricle points from the aortic pressure alone. A function is then used to draw a smooth line through these estimated points (described in Section “Waveform construction”) giving a continuous approximation of the left ventricle FTVE. The right ventricle is estimated in similar way. Known points from the estimated left ventricle FTVE, time values taken from the pulmonary artery pressure, and the global end-diastolic volume (GEDV) are correlated with points from the right ventricle FTVE, again creating equations used to estimated a continuous right ventricle FTVE.
It should be noted that the two catheter locations (aorta and pulmonary artery) can only “see” the ventricles during systole, inferring that information about the cardiac elastance can only be known during systole. However, this limitation is not clinically significant as cardiac elastance is primarily defined by what happens during systole.
The method presented was developed on a cohort of 9 porcine subjects divided into two groups. The first group was induced with pulmonary embolism, (5 subjects) [33, 34]. The second group (4 subjects) were induced with septic shock . Both trials were under the control of the ethics committee of the medical faculty of Liege, Belgium.
Measurements were taken every 30 minutes. In total the first group had 51 sets of measurements across the 5 pigs, while the septic shock cohort had 34 over 4 pigs. For each pig, the following relevant measurements were taken: P ao , P pa , left and right ventricle volume (V lv , V rv ), left and right ventricle pressure (P lv , P rv ). Aortic and pulmonary artery pressure were measured using catheters, while right and left ventricle pressures and volumes were directly measured using 7F, 12 electrodes conductance micro manometer tipped catheters.
Left and right cardiac elastances were calculated directly with the respective pressure and volume waveforms via (1), while setting the dead space volume to zero. GEDV was calculated as the sum of the individual maximum volumes of the left and right ventricles.
The data contains varying degrees of cardiac dysfunction from healthy to the fully developed disease state. It thus provides a good test of the presented methods. The use of two dysfunction states is designed to illustrate the method’s potential robustness.
The shape of the cardiac elastance waveform is formed from 3 known locations on the waveform; time, slope and height of the maximum rising gradient () and maximum falling gradient (), (analogous to MPG and MNG respectively in Figure 1), and the height of the waveform ′ s left shoulder, x 2, (analogous to LS, in Figure 1). These are 7 of the 8 values needed to reproduce the waveforms. The eighth is the height of the right shoulder (MX, in Figure 1) which is always unity. The 7 variable values are found with correlations from points on the P ao and P pa waveforms, as described in Section “Cardiac elastance correlations”.
where a β , b β and c β are similarly defined by replacing subscript 1 with 3 and setting x 2 = 1.
Cardiac elastance correlations
The right and left ventricle elastances are constructed using Equations (3) to (5). However, the 8 values that identify the separate elastances are found in different ways.
Reconstruction formulae for the three different types of correlations used to estimate the cardiac elastance
y = mx + c
y = m 1 x 1 + m 2 x 2 + c
y = m 1 x 1 + m 2 x 2+ m 3 x 3 + c
The left cardiac elastance values are relatively straightforward to find, and all use the standard correlations, α, from Table 1. The right cardiac elastance is slightly harder to estimate, partly because there is no direct down stream access to the right ventricle similar to that of the left ventricle. This is because only timing information is taken from P pa , which can be replaced with electrocardiography (ECG) and central venous pressure (CVP), as the P pa waveform is not commonly measured in an ICU. This self imposed limitation makes the method more clinically applicable.
This limitation necessitates other means to estimate the required values on the right side, which is mostly done through the use of multi-variable correlations involving GEDV, and the estimated left cardiac elastance. GEDV is closely connected to the average volume in the right ventricle, and thus serves to add additional scaling. The specific correlations found of the left and right elastances are shown in the Results.
The cardiac elastances were estimated for every measurement taken at 30 minute intervals during the onset of the porcine cohort’s cardiac dysfunction, in two different ways. First, using the correlations created using the entire cohort of 9 pigs (both dysfunctions). Second, using the correlations determined only from the five (or four) pigs from the “other" dysfunction. These estimations were then evaluated against the true cardiac elastance to determine the error, as in (6).
Left ventricle cardiac elastance correlations, all points use correlation α from Table 1
MP G x (x 1)
m = 2.521×10−5
c = 3.732×10−1
MP G t (t 1)
P ao (DMPG) t
m = 1.007
c = 1.615×10−3
m = 4.126×10−3
c = 8.665
MN G x (x 3)
P ao (MNG) t
m = 3.370×10−1
c = 4.886×10−1
MN G t (t 3)
P ao (DN) t
m = 1.028
c = −9.148×10−3
m = −5.510
c = 1.871×10+ 1
L S x (x 2)
m = 9.627×10−2
c = 5.867×10−1
Right cardiac elastance correlations
MPG s x (x 1)
GEDV, e lv (MPG) x
m 1 = −1.137×10−3
m 2 = 3.521×10−1
c = 4.524×10−1
MP G t (t 1)
P pa (MN) t
m = 1.023
c = −4.102×10−3
m 1 = −1.286×10−2
m 2 = 6.446×10−1
c = 6.910
MN G x (x 3)
GEDV, e lv (MNG) x
m 1 = 1.233×10−3
m 2 = 1.824×10−1
c = 2.348×10−1
MN G t (t 3)
P pa (DN) t
m = 9.751×10−1
c = −2.081×10−2
GEDV, e lv (MPG) x
m 1 = −6.096×10−2
m 2 = 1.446
c = −4.331
L S x (x 2)
GEDV, e lv (LS) x ,e lv (LS) t
m 1 = −2.054×10−3
m 2 = 3.209×10−1
m 3 = −3.538×10−1
c = 9.045×10−1
Reconstruction errors for pulmonary embolism and septic shock
Cross validation reconstruction errors for pulmonary embolism and septic shock
(based on septic shock)
(based on PE)
Table 5 shows the errors using the correlations derived from one dysfunction to approximate the driver function of the ‘other’ dysfunction. This last assessment cancels out any inherent dependence of the method, outlined in this paper, on the dysfunction itself and is a rigorous robustness test. The results show the method is independent of dysfunction, in this case, which is an extremely important result. For a method to reliably give insight into a specific patients cardiac dysfunction, it must not rely on the patient having a specific cardiac dysfunction. It is also another validation of the accuracy and potential of this method, which gives further confidence that the method will generalise over a wider set of dysfunctions.
A method for identifying time varying elastance of the left and right ventricles, was developed that doesn’t require an occlusion manoeuvre to find V d , and avoids invasive measurements of left and right ventricle volumes and pressure. These time varying elastance functions can be used as the input driver function into a lumped parameter cardiovascular model. The method tracks the main, clinically important time features on the driver functions, such as the start and end of ejection. This method was validated on two sets of porcine data, which included multiple measurements of five pigs that were induced with pulmonary embolism.
The approach taken to chose the mapping from the known metrics to the cardiac elastance, was based both on physiology and accuracy of the resulting cardiac elastance waveform. Thus some mappings may not be the most obvious choice based on purely physiological reasoning. Some of these mapping may appear to be counter-intuitive. For example, although MPG and MNG appear to have symmetry of meaning in the aortic pressure waveform, this does not transfer into symmetrical meaning for cardiac energetics, and hence do not have symmetry in the chosen mappings.
The correlations that were used in the reconstructions presented in this paper are generally strong. The time based correlations are especially strong (see the 2 left most plots in Figure 5), and lead to very good estimations of the waveform timing. There are also some rather poor correlations, the worst of which predicts the value of the right ventricle FTVE at the point MNG. This correlation has R = 0.31, and a correlation plot that barely shows any correlation at all. However, the final estimations of the waveform are not very sensitive to this value, and hence can produce reasonably accurate waveform estimations from values of MNG x that error prone. Note that the main results for this paper are the errors of the reconstructed waveforms, rather than the presence of correlations.
There is also a marked difference in the correlation’s R values between the left and right sides. This is mostly due to not have direct access to the right side of the heart (since all but the timing of the P pa waveform have been neglected), and instead are inferring information about the right side from the left side. Although this results in significantly lower R valued correlations (compared to the left), the waveform reconstructions are still good, with the errors averaging only 5% higher than the left side, for the cross dysfunction validation.
As would be expected from the correlations in Table 2, the timings of the reconstructed waveforms were very good, capturing the ascending and descending exponential sections very closely across a large change in cardiac condition, including an approximate doubling in heart rate. The biggest difference in the real and reconstructed waveforms usually appeared in the flat section (B of Figure 3). This error could be reduced by using a higher order polynomial or model. However, simulation has shown that this error does not have a significant effect on outputs of the existing cardiovascular model . In addition the philosophy of this research is to only add complexity as required to capture and predict clinical data, and to keep the model as simple as possible. The straight line approach would also have the benefit of simplifying any semi-analytical solution approach to speeding up the simulations which is planned in future work .
The specific waveform for FTVE was chosen as a balance between accuracy and complexity. A model (waveform shape) was needed that provided enough information to capture the necessary cardiac dynamics, while remaining as simple as possible. The simplest waveform shape is that of three straight lines. To fit this model, four points would be required (or equivalently, two points with associated gradients and two heights (points with free ‘x’ coordinate) - the currently used data, see Figure 1). However, for the same amount of information, and no additional computational complexity, a model of two exponential curves and one straight line (the chosen model) can be fitted. Furthermore, this model is continuously differentiable, which can be computationally advantageous in numerical simulation. This more complex shape, provides more accurate results, without adding complexity to the process of fitting. Beyond this model structure, more complex models would require additional information to be found from the pressure waveforms.
For simplicity of the model, V d was set to zero, which is a limitation of the current approach. However, even with this discarded parameter, the model is still able to capture the necessary dynamics, and provide useful additional information to an ICU clinician base on the trends observed in FTVE shape relative to this assumption. Having an accurate V d value in addition would add useful information, but is not identifiable from the data available without added invasive sensors.
Another limitation of the paper is the small number (N = 9) of porcine subjects (although there were a large number (N = 85) of data sets used) and only two dysfunctions. However, this paper serves as a proof of concept as its main goal. In particular, the results are intended to justify the methods and thus a larger, complete prospective validation.
Although this study used the pulmonary artery pressure, the points found in the pulmonary artery pressure waveform could be estimated using the central venous pressure and ECG. Thus, there is the potential to construct the driver functions without any knowledge of the right side of the heart. This option would avoid using pulmonary artery pressure, which is less commonly available than aortic pressure in an ICU. This study did not have access to CVP and ECG data for the porcine cohort so used pulmonary artery pressure instead. Future work will investigate these clinically relevant extensions to minimise invasiveness, including the application of the methods onto ICU patients.
This paper describes a method that has huge clinical potential. Currently, disease states like pulmonary embolism and septic shock are difficult to diagnose accurately and reliably due to the lack of useful information accessible to the clinician. Knowing the continuous cardiac elastance can gather this otherwise under utilised information into a more clinically applicable data set. This is achieved through a method that uses only currently available metrics in the ICU, and thus is no more invasive (no additional risk) than standard clinical practice. Preliminary results for pulmonary embolism indications from the cardiac elastance show promise.
This work was supported by the Foundation for Research Science and Technology (FRST) of New Zealand, the Fonds de la Recherche Scientifique (FNRS, Belgium) and the French Community of Belgium (Actions de Recherches Concertées Académie Wallonie-Europe).
- Guyton A, Hall J: Textbook of Medical Physiology. Philadelphia: W B Saunders Company; 2000.Google Scholar
- Grenvik A, Ayres SM, Holbrook PR: Textbook of Critical Care. Philadelphia: W B Saunders Company; 1989.Google Scholar
- Angus DC, Linde-Zwirble WT, Lidicker J, Clermont G, Carcillo J, Pinsky MR: Epidemiology of severe sepsis in the United States: analysis of incidence, outcome, and associated costs of care. Crit Care Med 2001, 29(7):1303–10. 10.1097/00003246-200107000-00002View ArticleGoogle Scholar
- Kearon C: Diagnosis of pulmonary embolism. CMAJ 2003, 168(2):183–94.Google Scholar
- Pineda LA, HV S, Grant BJB: Clinical suspicion of fatal pulmonary embolism. Chest 2001, 120(3):791–5. 10.1378/chest.120.3.791View ArticleGoogle Scholar
- Frazier S, Skinner GJ: Pulmonary artery catheters: state of the controversy. J Cardiovasc Nursing 2008, 23: 113–21.View ArticleGoogle Scholar
- Chatterjee K: The Swan-Ganz Catheters: Past, Present, and Future. Circulation 2009, 119: 147–152. 10.1161/CIRCULATIONAHA.108.811141View ArticleGoogle Scholar
- Cooper A, Doig WJ, Sibbald GS: Pulmonary artery catheters in the critically ill. Crit Care Med 1996, 12: 777–94.Google Scholar
- Weber KT, Janicki JS: The heart as a muscle-pump system and the concept of heart-failure. Am Heart J 1979, 98(3):371–84. 10.1016/0002-8703(79)90051-6View ArticleGoogle Scholar
- Ross JJ: Afterload mismatch and preload reserve: a conceptual framework for the analysis of ventricular function. Prog Cardiovasc Dis 1976, 18: 255–64. 10.1016/0033-0620(76)90021-9View ArticleGoogle Scholar
- Smith BW, Chase JG, Nokes RI, Shaw GM, Wake G: Minimal haemodynamic system model including ventricular interaction and valve dynamics. Med Eng Phys 2004, 26(2):131–9. 10.1016/j.medengphy.2003.10.001View ArticleGoogle Scholar
- Starfinger C, Chase JG, Hann CE, Shaw GM, Lambermont B, Ghuysen A, Kolh P, Dauby PC, Desaive T: Model-based identification and diagnosis of a porcine model of induced endotoxic shock with hemofiltration. Math Biosci 2008, 216(2):132–9. 10.1016/j.mbs.2008.08.014MathSciNetView ArticleGoogle Scholar
- Chung DC, Niranjan SC, Clark J, 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.Google Scholar
- Revie J, Stevenson D, Chase G, Hann CE, Lambermont B, Ghuysen A, Kolh P, Morimont P, Shaw GM, Desaive T: Clinical detection and monitoring of acute pulmonary embolism: proof of concept of a computer-based method. Ann of Intensive Care 2011., 1(33):Google Scholar
- Revie J, Stevenson D, Chase JG, Hann CE, Lambermont B, Ghuysen A, Kolh P, Shaw GM, Heldmann S, Desaive T: Validation of subject-specific cardiovascular system models from porcine measurements. Comput Methods Programs Biomed 2011. in-pressGoogle Scholar
- Suga H, Sagawa K: Instantaneous pressure-volume relationships and their ratio in the excised, supported canine left ventricle. Circ Res 1974, 35: 117–26. 10.1161/01.RES.35.1.117View ArticleGoogle Scholar
- Sagawa K: The end-systolic pressure-volume relation of the ventricle: definition, modifications and clinical use. Circulation 1981, 63(6):1223–7. 10.1161/01.CIR.63.6.1223View ArticleGoogle Scholar
- Suga H, Sagawa K, Shoukas AA: Load independence of the instantaneous pressure-volume ratio of the canine left ventricle and effects of epinephrine and heart rate on the ratio. Circ Res 1973, 32(3):314–22. 10.1161/01.RES.32.3.314View ArticleGoogle Scholar
- Burkhoff D, Sugiura S, Yue DT, Sagawa K: Contractility-dependent curvilinearity of end-systolic pressure-volume relations. Am J Physiol 1987, 252(6 Pt 2):H1218–27.Google Scholar
- Senzaki H, Chen CH, Kass DA: Single-beat estimation of end-systolic pressure-volume relation in humans. A new method with the potential for noninvasive application. Circulation 1996, 94(10):2497–506. 10.1161/01.CIR.94.10.2497View ArticleGoogle Scholar
- Sunagawa K, Sagawa K, Maughan WL: Ventricular interaction with the loading system. Ann Biomed Eng 1984, 12(2):163–89. 10.1007/BF02584229View ArticleGoogle Scholar
- Guarini M, Urzúa J, Cipriano A, González W: Estimation of caridac function from computer analysis of the arterial pressure waveform. IEEE Trans on Biomed Eng 1998, 45(12):1420–8. 10.1109/10.730436View ArticleGoogle Scholar
- Swamy G, Kuiper J, Gudur MSR, Bari ON, Mukkamala R: Continuous left ventricular ejection fraction monitoring by aortic pressure waveform analysis. Ann Biomed Eng 2009, 37(6):1055–68. 10.1007/s10439-009-9675-4View ArticleGoogle Scholar
- Shishido T, Hayashi K, Shigemi K, Sato T, Sugimachi M, Sunagawa K: Single-beat estimation of end-systolic elastance using bilinearly approximated time-varying elastance curve. Circulation 2000, 102(16):1983–9. 10.1161/01.CIR.102.16.1983View ArticleGoogle Scholar
- ten Brinke EA, Klautz RJ, Verwey HF, van der Wall EE, Dion RA, Steendijk P: Single-beat estimation of the left ventricular end-systolic pressure-volume relationship in patients with heart failure. Acta Physiol (Oxf) 2010, 198: 37–46. 10.1111/j.1748-1716.2009.02040.xView ArticleGoogle Scholar
- Burkhoff D: Chasing the elusive pressure-volume relationship. J Am Coll Cardiol Img 2009, 2: 1282–4.View ArticleGoogle Scholar
- Starfinger C, Hann CE, Chase JG, Desaive T, Ghuysen A, Shaw GM: Model-based cardiac diagnosis of pulmonary embolism. Comput Methods Programs Biomed 2007, 87: 46–60. 10.1016/j.cmpb.2007.03.010View ArticleGoogle Scholar
- Sagawa K: The ventricular pressure-volume diagram revisited. Circ Res 1978, 43(5):677–87. 10.1161/01.RES.43.5.677View ArticleGoogle Scholar
- Suga H: Ventricular energetics. Physiol Rev 1990, 70(2):247–77.Google Scholar
- Burkhoff D, Sagawa K: Ventricular efficiency predicted by an analytical model. Regu Physiol 1986, 250(6):R1021–7.Google Scholar
- Broscheit JA, Weidemann F, Strotmann J, Steendijk P, Karle H, Roewer N, Greim CA: Time-varying elastance concept applied to the relation of carotid arterial flow velocity and ventricular area. J Cardiothorac Vasc Anesth 2006, 20(3):340–6. 10.1053/j.jvca.2006.01.014View ArticleGoogle Scholar
- Suga H: Total mechanical energy of a ventricle model and cardiac oxygen consumption. Heart 1979, 236(3):H498–5.Google Scholar
- Desaive T, Dutron S, Lambermont B, Kolh P, Hann CE, Chase JG, Dauby PC, Ghuysen A: Close-loop model of the cardiovascular system including ventricular interaction and valve dynamics: application to pulmonary embolism. In IFMBE Proceedings, 12th Intl Conference on Biomedical Engineering (ICBME): Dec 7–10;. Singapore; 2005:4. ISSN 1727–1983Google Scholar
- Ghuysen A, Lambermont B, Kolh P, Tchana-Sato V, Magis D, Gerard P, Mommens V, Janssen N, Desaive T, D’Orio V: Alteration of right ventricular-pulmonary vascular coupling in a porcine model of progressive pressure overloading. Shock 2008, 29(2):197–204.Google Scholar
- Lambermont B, Ghuysen A, Kolh P, Tchana-Sato V, Segers P, Gerard P, Morimont P, Magis D, Dogne JM, Masereel B, D’Orio V: Effects of endotoxic shock on right ventricular systolic function and mechanical efficiency. Cardiovasc Res 2003, 59(2):412–8. 10.1016/S0008-6363(03)00368-7View ArticleGoogle Scholar
- Hann CE, Chase JG, Shaw GM: Efficient implementation of non-linear valve law and ventricular interaction dynamics in the minimal cardiac model. Comput Methods Programs Biomed 2005, 80: 65–74. 10.1016/j.cmpb.2005.06.003View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.