Non-invasive evaluation of energy loss in the pulmonary arteries using 4D phase contrast MR measurement: a proof of concept
© Lee et al.; licensee BioMed Central Ltd. 2013
Received: 8 August 2013
Accepted: 12 September 2013
Published: 23 September 2013
The repair surgery of congenital heart disease (CHD) associated with the right ventricular (RV)-pulmonary artery (PA) pathophysiology often left patients with critical post-operative lesions, leading to regurgitation and obstruction in the PAs. These lesions need longitudinal (with time) assessment for monitoring the RV function, in order for patients to have appropriate treatment before irreversible RV dysfunction occurs. In this research, we computed energy loss in the branch PAs using blood flow and pressure drop data obtained from 4D phase contrast (PC) MRI, to non-invasively quantify the RV-PA pathophysiology.
4D PC MRI was acquired for a CHD patient with abnormal RV-PA physiology, including pulmonary regurgitation and PA stenosis, and a subject with normal RV-PA physiology. The blood velocity, flow rate, and pressure drop data, obtained from 4D PC MRI, were used to compute and compare the energy loss values between the patient and normal subjects.
The pressure drop in the branch PAs for the patient was −1.3 mmHg/s and −0.2 mmHg/s for the RPA and LPA, respectively, and was larger (one order of magnitude) than that for the control. Similarly, the total energy loss in the branch PAs for the patient, -96.9 mJ/s and −16.4 mJ/s, for the RPA and LPA, respectively, was larger than that for the control.
The amount of energy loss in the pulmonary blood flow for the patient was considerably larger than the normal subject due to PA regurgitation and PA stenosis. Thus, we believe that the status of RV-PA pathophysiology for CHD patients can be evaluated non-invasively using energy loss endpoint.
Keywords4D phase contrast MRI Pulmonary insufficiency Energy-based endpoint
With the excellent survival rate of palliated congenital heart disease (CHD) patients, monitoring residual lesions has become increasingly important [1, 2]. Particularly, patients with right ventricular (RV) or pulmonary valve lesions, such as tetralogy of Fallot (TOF), aortic valve disease requiring the Ross procedure (aortic autograft with RV-PA homograft), or complex transposition of the great arteries (TGA), are often left with pulmonary insufficiency (PI) leading to progressive RV dilatation and occasionally resulting in pressure overload due to residual pulmonary stenosis. Those sequelae can result in progressive RV myocardial dysfunction, increasing the risk of sudden death [3–6].
Surgical or catheter-based pulmonary valve replacement is often required to rectify severe RV myocardial dysfunction . Since the appropriate time for intervention is critical, the patho-physiology of the RV and the pulmonary arteries (PAs) has to be carefully monitored throughout the patient’s lifetime. However, due to the complexity of symptoms, it is sometimes difficult to accurately assess the disease progression with existing cardiac indices alone. Commonly used metrics include body-surface-area (BSA) indexed RV end-diastolic and end-systolic volumes (EDVI and ESVI, respectively), ejection fraction, and RV end-systolic pressure (ESP).
Recently, energy-based endpoints have been investigated in our research group [8–10] to help determine timing for surgical interventions. In particular, BSA indexed RV stroke work (SWI) and energy transfer ratio (e MPA ) between the RV and main PA (MPA) were proposed in our previous study to evaluate the hemodynamic status of the RV and the PA for repaired CHD patients. These energy-based endpoints had the advantage of incorporating RV volume, pressure data, and flow conditions in the PAs into one single index, and differentiated the hemodynamics of RV and PA of patients from those of normal subjects with statistical significance (p <0.05). In addition, energy-based endpoints correlated well with current indices, such as RV EDVI and RV ESP . Further, other energy-based endpoints, such as energy dissipation and power loss, also have been evaluated in single-ventricle (Fontan) physiology from other group [11–13].
However, energy-based endpoints require invasive pressure measurement, i.e., cardiac catheterization, which limits their applicability to only those patients undergoing catheterization. With recent development of 4D phase contrast magnetic resonance imaging (PC MRI), three dimensional and three directional velocity data over the cardiac cycle can be obtained for the conduits and chambers of the entire heart [14–16]. Therefore, the pressure data can be estimated non-invasively for any of the heart’s conduits and chambers by using the time varying 3D velocity vector field from 4D PC MRI data [16, 17].
In this research we computed the pressure drop along the PA non-invasively, from the MPA to the branch PAs (RPA and LPA, right and left PA, respectively), using 4D PC MRI data. This enabled us to calculate energy loss in the branch PAs over the cardiac cycle leading to quantification of localized PI caused by obstruction in the RVOT and PA. The novelty of this research is that the use of 4D PC MRI will allow non-invasive assessment of energy loss in the RV-PA pathophysiology; thus, avoiding the need for catheterization. We believe that non-invasive assessment of pressure drop along the PAs and subsequent energy loss calculation in the PAs will help in longitudinal monitoring of RV-PA physiology for CHD patients. This allows accurate assessment of progression of the disease, which may lead to improving the timing of intervention, and patient outcomes.
Demographics of the subject in the study
Control (normal RV-PA physiology)
Patient (abnormal RV-PA physiology)
Detail of 4D PC MRI acquisition parameters
Spatial resolution (mm)
128 × 128
128 × 128
Pixel size (mm)
2.34 × 2.34
Slice spacing (mm)
Acquisition volume (cm)
32.0 × 32.0 × 6.0
30.0 × 30.0 × 10.0
Repetition time (ms)
Echo time (ms)
Flip angle (°)
Diameter, average flow rate and velocity in the PAs for the control and patient
The streamlines of each pulmonary flow, MPA-RPA and MPA-LPA flows, and the velocity vectors in the respective subplanes during a systole phase are also shown in Figure 5. The streamlines that passed through the predominant portion of each pulmonary flow were chosen for the further pressure drop computation. The streamlines were started from uniformly distributed 3 × 3 seed points with a distance of approximately 2 mm on the plane. The time varying pressure drop calculation between the branch PAs and the MPA over the cardiac cycle was performed along those streamlines (described in the following section). Multiple streamlines were used for the pressure drop calculation since the calculated pressure drop may be dependent upon the selected streamline. Thus, the average pressure drop values for each pulmonary flow were reported in this study.
Pressure drop calculation
where, u is blood velocity data, ρ the blood density (= 1,050 kg/m3), μ the blood viscosity (= 0.00345 Pa/s).
Thus, there are two dPBranch PAs; one is dPRPA and the other is dPLPA.
The calculation procedures (Figure 2) to compute pressure drop were verified with the pressure drop computed for a simplified 2D stenosis model using a finite difference CFD solver (Fluent, ANSYS, Inc., Canonsburg, PA, USA). The detail of the verification procedures and results are provided in the Appendix section.
Energy loss calculation
The is the difference in the rate of the pressure-flow and kinetic energies transferred between the branch PAs and the MPA as shown in Eq. 5a:
where, dP BranchPA is the time varying pressure drop in the branch PA computed previously. The first two terms of the right-hand side of Eq. 5e account for the rate of pressure-flow and kinetic energy losses in the RPA, respectively, and the next two terms represent those energy losses in the LPA.
In this study, the patient had abnormal RV-PA physiology compared to the control subject. To highlight this, the PA diameters measured from 4D MRI along with the spatially averaged blood velocity and flow rate at each PAs are presented in the following results section. Subsequently, the results for non-invasively computed pressure drop and energy loss in the branch PAs for the control and the patient are provided.
Average and peak pressure drop, and average and net total energy loss in the PAs for the control and patient
Ave. pressure drop (mmHg/s)
Peak pressure drop (mmHg)
Ave. total energy loss (mJ/s)
Net total energy loss (mJ)
Time averaged blood velocity and flow rate at the PAs
The time averaged blood velocities at the PAs for both the subjects are presented in Table 4. The time average blood velocities for the control were 18.4 cm/s, 25.9 cm/s, and 24.6 cm/s for the MPA, RPA, and LPA, respectively. The time average blood velocity at the MPA was lower than that at the branch PAs, because the area of MPA is larger (diameter of the MPA, DMPA = 2.5 cm; Table 4) compared to the sum of the areas of the two branch PAs (DRPA = 1.5 cm and DLPA = 1.7 cm). The time average blood velocities for the patient were 18.2 cm/s, 15.6 cm/s, and 9.7 cm/s for the MPA, RPA, and LPA, respectively. The time average blood velocities at the PAs for the patient were lower than the control due to the negative flow in the patient’s PAs. Further, the average blood velocity was lower at the LPA than the RPA due to considerably larger regurgitation that occurred at the LPA during the diastole phase of the pulse for the patient.
In Figure 6, the blood flow rate versus time curves at the PAs for both the subjects are shown. The control had the average blood flow volumes of 88.6 ml/s, 47.7 ml/s, and 56.3 ml/s for the MPA, RPA, and LPA, respectively (Figure 6A and Table 4). The summation of RPA and LPA average flow was 103.9 ml/s and was somewhat more than the MPA average blood flow (88.6 ml/s) by 17.3%. The small amount of reverse blood flow in the MPA was observed during the early and late diastole phases. No reverse flow was observed in the branch PAs, RPA and LPA.
In contrast, the time average blood flow volumes for the patient were 96.0 ml/s, 48.5 ml/s, and 44.9 ml/s for the MPA, RPA, and LPA, respectively (Figure 6B and Table 4). The summation of RPA and LPA flow was 93.4 ml/s and was within 2.7% of the MPA flow. As shown in Figure 6B, the large amount of reverse blood flow in the MPA and the LPA was observed during diastole phase. During the early systole over 70% of blood from the MPA directly flows to the LPA; however, the blood flow in the RPA increased with time until it reached the late systole phase. Reverse flow in the RPA was insignificant during the diastole phase.
Pressure drop along the branch PAs
Figure 7 shows the pressure drop between the branch PA and the MPA (dP Branch PA = PBranch PA − PMPA) versus time. The average pressure drop values at any instant of time calculated from nine streamlines for each pulmonary flow was used to construct the time varying pressure drop curves with standard error. As presented in Table 3, the average pressure drop was minimal in the control (−0.09 mmHg/s and −0.04 mmHg/s, for the RPA and LPA, respectively) compared to that of the patient (−1.3 mmHg/s and −0.2 mmHg/s, for the RPA and LPA, respectively). Similarly, the peak pressure drop was smaller for the control (−0.4 mmHg and −0.2 mmHg, for the RPA and LPA) than that for the patient (−4.7 mmHg and −0.8 mmHg, for the RPA and LPA).
Energy loss in the branch PAs
The rate of total energy loss along the branch PA () versus time curves are shown in Figure 8. The average total energy loss was larger for the patient (−96.9 mJ/s and −16.4 mJ/s, for the RPA and LPA, respectively) compared to that for the control (−2.1 mJ/s and −1.1 mJ/s, for the RPA and LPA). Also, the net total energy loss in the PAs over the cardiac cycle was much larger for the patient (−70.3 mJ and −11.9 mJ, for the RPA and LPA) than that for the control (−1.6 mJ and −0.8 mJ, for the RPA and LPA).
Pulmonary regurgitation, residual RVOT and PA obstruction are the main post-operative lesions seen in many cases of repaired CHD [2, 6, 19]. The abnormal RV-PA hemodynamics, resulting from pulmonary regurgitation and residual obstruction, directly affect the efficiency of pulmonary blood flow. Further, chronic pulmonary regurgitation and PA obstruction can cause progressive RV volume and pressure overload resulting in RV dilatation and hypertrophy. As these symptoms worsen, RV dysfunction can occur, leading to often sudden death.
The patient had significantly larger energy loss compared to the control subject as shown in Figure 8. The average energy loss in the patient was −96.9 mJ/s and −16.4 mJ/s, for the RPA and LPA, respectively. These values for the patient are one order of magnitude higher than those for the control subject (Table 3). The energy loss in the RPA of the patient was particularly higher due to the stenosis as mentioned in the previous paragraph. The higher pressure drop in the patient leads to the increased energy loss in the branch PAs, as shown in Figure 8.
The additional energy loss due to blood flow separation in the MPA bifurcation, in Eq. 7, which is the minor energy loss from the perspective of a fluid mechanics, was calculated by the multiplication of a flow resistance coefficient (ζ) and kinetic energy loss and was added to the total energy loss. A flow resistance coefficient (ζ) was determined based on the ratio of flow rates in the PAs, Q Branch PA /Q MPA , for each time point during the cardiac cycle . The contribution by minor energy loss to the total energy loss was less significant for the control subject (2.1% and 4.1% for the RPA and LPA, respectively) than the patient (19.0% and 7.1% for the RPA and LPA, respectively). Particularly, the minor energy loss for the patient’s RPA was relatively large (19.0%) due to flow obstruction in the RPA caused by the stenosis.
Inherent noise exists on phase contrast velocity mapping due to the motion of PAs during the contraction and relaxation of the RV . It is observed that noise resulting from PA movement was manifested near the boundary of artery. Therefore, the streamlines which passed adjacent to the PA boundary were not considered in computing the pressure drop along the PA.
The elastic energy storage in the MPA was estimated to be about 1% of the total energy transferred to the MPA . The majority of the energy transfer occurs through the form of the pressure-flow energy and the kinetic energy, as discussed in our previous study . Thus, elastic energy loss in the branch PAs would be insignificant. However, other possible energy losses in the branch PAs, such as friction loss between blood and tissues as blood flows in the PAs and local effects due to narrowing or stenosis in the daughter branches of the branch PAs, are difficult to measure under a clinical setting. Therefore, the energy loss calculation in this study may somewhat underestimate the actual energy loss in the branch PAs while such losses need to be assessed in future. We believe such losses are order of magnitude lower than pressure-flow and kinetic energy losses, which are accounted for in this study.
The accurate measurement of cardiac pressure data is necessary for the previously proposed energy-based endpoints, such as RV stroke work index (SWI) and energy transfer ratio (e MPA ). However, diagnostic cardiac catheterization is not a part of standard care for typical CHD patients. Therefore, continuous quantification and monitoring of the change in the pressure using catheterization for a patient is not trivial. In this regard, 4D PC MRI has an advantage as it can provide 4D (u, v, w, t) velocity information needed to compute pressure drop, flow, and energy-based endpoints non-invasively . Further, the methodology used for this research does not include a complicated numerical computation, such as computational fluid dynamics (CFD), which is time consuming and requires significant lead time to pre-process, compute, post process, and analyze the numerical data. Therefore, the proposed technique would extend the clinical applicability of the energy-based endpoints for longitudinal assessment of patient on a regular basis. A study that extends the proposed methodology and involves statistically relevant control subjects and CHD patients is needed. We believe such a study would reveal the significant difference in the pressure drop and energy loss in the branch PAs between control subjects and CHD patients.
We have used non-invasive 4D PC MRI data to calculate the energy loss in the branch PAs and compared between a CHD patient with abnormal PA physiology and a normal subject with normal PA physiology. Based on our results, the patient had considerably larger energy loss in the branch PAs due to pulmonary regurgitation and PA obstruction compared to the normal subject. Therefore, we believe that non-invasively obtained energy loss in the branch PAs can be a useful clinical measure for evaluating the longitudinal changes of RV-PA pathophysiology of repaired CHD patients.
CHD = congenital heart disease, RV = right ventricle, PA = pulmonary artery, 4D PC MRI = 4 dimensional phase contrast magnetic resonance imaging, TOF = tetralogy of Fallot, TGA = transposition of the great arteries, PI; pulmonary insufficiency, BSA = body surface area, EDV = end-diastolic volume [ml], ESV = end-systolic volume [ml], SW = stroke work [J], MPA = main pulmonary artery, RPA = right pulmonary artery, LPA = left pulmonary artery, e MPA = energy transfer ratio, CFD = computational fluid dynamics, CCHMC = Cincinnati Children’s Hospital Medical Center, VENC = velocity encoding, AP = anterior to posterior, RL = right to left, FH = foot to head, ∇P = pressure gradient field, MRA = magnetic resonance angiogram, u = blood velocity [cm/s], Q = blood flow rate [ml/s], ρ = blood density [kg/m3], μ = blood viscosity [Pa/s], dP = pressure drop [mmHg], , = the rate of energy loss in the branch PA [mJ/s], = the rate of major energy loss [mJ/s], = the rate of minor energy loss [mJ/s], = the rate of energy in the PA [mJ/s].
The verification of the pressure drop computation
Physiologic pressure data in the subject's PA was not available because catheterization was not performed in this study. Thus, the pressure drop calculation in this study was verified with CFD result using a simplified geometry as following.
Pressure drop computation
After CFD analysis, the 2D velocity information was extracted from the converged solution. The pressure gradient field () was computed using the extracted velocity field (Eq. 1). Then, in Ensight six streamlines originating from the inlet were generated with a distance of approximately 2.5 mm and the pressure gradient along each streamlines () was obtained (Eq. 2). The pressure drop at the outlet was computed by integrating the pressure gradient along each streamline (Eq. 3). The average pressure drop value was obtained from six streamlines (0.22 mmHg) using the proposed method (Figure 2) and was compared with the pressure drop computed from CFD solution (0.24 mmHg). The difference in the pressure drop between two methodologies was 7.3% [= (0.24-0.22)/0.24 × 100].
The comparison of pressure drop and energy loss in the branch PAs
Comparison of average pressure drop and total energy loss in the branch PA in this study with physiological data
Total energy loss
4D PC MRI
4D PC MRI
This work was partially supported by “TOF study” funded by an internal project in CCHMC. We appreciate Dr. Michael Markl for helping us in developing an initial version of Matlab function to convert 4D PC MRI data into the Ensight format.
- Ooi A, Moorjani N, Baliulis G, Keeton BR, Salmon AP, Monro JL, Haw MP: Medium term outcome for infant repair in tetralogy of fallot: indicators for timing of surgery. Eur J Cardiothorac Surg 2006, 30(6):917–922. 10.1016/j.ejcts.2006.08.022View ArticleGoogle Scholar
- Murphy JG, Gersh BJ, Mair DD, Fuster V, McGoon MD, Ilstrup DM, McGoon DC, Kirklin JW, Danielson GK: Long-term outcome in patients undergoing surgical repair of tetralogy of fallot. N Engl J Med 1993, 329: 593–599. 10.1056/NEJM199308263290901View ArticleGoogle Scholar
- Taussig HB: Diagnosis of the tetralogy of fallot and medical aspects of the surgical treatment. Bull N Y Acad Med 1947, 23: 705–718.Google Scholar
- d’Udekem Y, Ovaert C, Grandjean F, Gerin V, Cailteux M, Shango-Lody P, Vliers A, Sluysmans T, Robert A, Rubay J: Tetralogy of Fallot: transannular and right ventricular patching equally affect late functional status. Circulation 2000, 102: III116–122.Google Scholar
- Davlouros PA, Niwa K, Webb G, Gatzoulis MA: The right ventricle in congenital heart disease. Heart 2006, 92(Suppl 1):i27–38.View ArticleGoogle Scholar
- Oosterhof T, Tulevski II, Vliegen HW, Spijkerboer AM, Mulder BJ: Effects of volume and/or pressure overload secondary to congenital heart disease (tetralogy of fallot or pulmonary stenosis) on right ventricular function using cardiovascular magnetic resonance and B-type natriuretic peptide levels. Am J Cardiol 2006, 97: 1051–1055. 10.1016/j.amjcard.2005.10.047View ArticleGoogle Scholar
- Harrild DM, Berul CI, Cecchin F, Geva T, Gauvreau K, Pigula F, Walsh EP: Pulmonary valve replacement in tetralogy of fallot: impact on survival and ventricular tachycardia. Circulation 2009, 119: 445–451. 10.1161/CIRCULATIONAHA.108.775221View ArticleGoogle Scholar
- Das A, Banerjee RK, Gottliebson WM: Right ventricular inefficiency in repaired tetralogy of fallot: proof of concept for energy calculations from cardiac MRI data. Ann Biomed Eng 2010, 38: 3674–3687. 10.1007/s10439-010-0107-2View ArticleGoogle Scholar
- Lee N, Das A, Gottliebson W, Banerjee RK: Comparison of stroke work between repaired tetralogy of fallot and normal right ventricular physiologies. Heart Vessels 2011. DOI: 10.1007/s00380–011–0212–7Google Scholar
- Lee N, Das A, Taylor M, Hor K, Banerjee RK: Energy transfer ratio as a metric of right ventricular efficiency in repaired congenital heart disease. Congenit Heart Dis 2012. DOI: 10.1111/chd.12034Google Scholar
- Dasi LP, Pekkan K, Katajima HD, Yoganathan AP: Functional analysis of fontan energy dissipation. J Biomech 2008, 41: 2246–2252. 10.1016/j.jbiomech.2008.04.011View ArticleGoogle Scholar
- Dasi LP, Pekkan K, de Zelicourt D, Sundareswaran KS, Krishnankutty R, Delnido PJ, Yoganathan AP: Hemodynamic energy dissipation in the cardiovascular system: generalized theoretical analysis on disease states. Ann Biomed Eng 2009, 37: 661–673. 10.1007/s10439-009-9650-0View ArticleGoogle Scholar
- Whitehead KK, Pekkan K, Kitajima HD, Paridon SM, Yoganathan AP, Fogel MA: Nonlinear power loss during exercise in single-ventricle patients after the fontan: insights from computational fluid dynamics. Circulation 2007, 116: I165–171.View ArticleGoogle Scholar
- Eriksson J, Dyverfeldt P, Engvall J, Bolger AF, Ebbers T, Carlhall CJ: Quantification of presystolic blood flow organization and energetics in the human left ventricle. Am J Physiol Heart Circ Physiol 2011, 300: H2135–2141. 10.1152/ajpheart.00993.2010View ArticleGoogle Scholar
- Fredriksson AG, Zajac J, Eriksson J, Dyverfeldt P, Bolger AF, Ebbers T, Carlhall CJ: 4-D blood flow in the human right ventricle. Am J Physiol Heart Circ Physiol 2011, 301: H2344–2350. 10.1152/ajpheart.00622.2011View ArticleGoogle Scholar
- Markl M, Kilner PJ, Ebbers T: Comprehensive 4D velocity mapping of the heart and great vessels by cardiovascular magnetic resonance. J Cardiovasc Magn Reson 2011, 13: 7. 10.1186/1532-429X-13-7View ArticleGoogle Scholar
- Ebbers T, Wigstrom L, Bolger AF, Engvall J, Karlsson M: Estimation of relative cardiovascular pressures using time-resolved three-dimensional phase contrast MRI. Magn Reson Med 2001, 45: 872–879. 10.1002/mrm.1116View ArticleGoogle Scholar
- Idelchik IE: Handbook of Hydraulic Resistance. Connecticut: Begell House; 2001.Google Scholar
- Pigula FA, Khalil PN, Mayer JE, del Nido PJ, Jonas RA: Repair of tetralogy of fallot in neonates and young infants. Circulation 1999, 100: II157–161.View ArticleGoogle Scholar
- Andersen AH, Kirsch JE: Analysis of noise in phase contrast MR imaging. Med Phys 1996, 23: 857–869. 10.1118/1.597729View 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.