|
|
||||||||
a Department of Cardiovascular Surgery, Kitasato University School of Medicine, Sagamihara, Japan
b Department of Pediatrics, Kitasato University School of Medicine, Sagamihara, Japan
c Department of Cardiac Surgery, The University of Tokyo, Graduate School of Medicine, Tokyo, Japan
Accepted for publication April 28, 2009.
* Address correspondence to Dr Miyaji, Kitasato University School of Medicine, 1-15-1 Kitasato, Sagamihara, Kanagawa, 225-8555, Japan (Email: kagami111{at}aol.com).
Presented at the Forty-fifth Annual Meeting of The Society of Thoracic Surgeons, San Francisco, CA, Jan 26–28, 2009.
| Abstract |
|---|
|
|
|---|
Methods: For the Fontan operation, a conduit from 14 to 22 mm was created based on angiographic data from 17 Fontan patients (mean age, 36.0 months; mean body surface area, 0.53 m2). Respiratory-driven flow of the superior and inferior vena cava was determined at rest and during exercise on two levels (0.5 and 1.0 W/kg) by magnetic resonance imaging flow studies. Flow stagnation was defined as the volume of the region where flow velocity was less than 0.01 m/second at both the expiratory and inspiratory phases.
Results: In larger conduits, backward flow at the expiratory phase was prominent. Energy loss was small even during exercise, but the change was slightly larger between 14 and 16 mm than other conduit sizes (14 mm, 5.759 mW; 16 mm, 4.881 mW; and 22 mm, 4.199 mW during 1.0 W/kg exercise). Stagnation volume at the expiratory phase increased with an increase of conduit size (14 mm, 9.20% vs 22 mm, 33.9% conduit volume at rest).
Conclusions: Fontan circulation is a low-energy system even during exercise. Larger conduits were proven to have redundant spaces, thus 16 and 18 mm conduits were optimal.
| Introduction |
|---|
|
|
|---|
There are a number of reports regarding Fontan hemodynamics based on computational flow dynamic (CFD) models [8–20]. As the fundamental Fontan physiology is a dissociation of the venous return from a ventricular power source, all of the previous numeric models have assumed that pulmonary blood flow is steady and continuous [3, 8–19] except the numeric study by Marsden and colleagues [20] that simulates the effects of exercise and respiration in two cases. However, clinical observations have proved that the pulmonary blood flow in the Fontan physiology is mostly dependent on the respiratory motion [2–4] and often referred to as a "respiratory pump" [2, 4], indicating that pulmonary blood flow is sucked into the lung during inspiration. In Fontan circulation, the decreased exercise capacity is one of the significant morbidities and detailed studies have been reported in regard to flow dynamics of the superior vena cava (SVC) and the inferior vena cava (IVC) of Fontan patients at rest, and during exercise in relation to respiratory flow fluctuation [3, 4, 21–24]. These studies demonstrated that augmentation of pulmonary blood flow with inspiration has important physiologic effects [2]. Therefore, we created respiratory-driven transient flow models using these flow study data.
Previous reports of CFD Fontan models have discussed, mainly, energy expenditure at the anastomosis site and pulmonary blood flow distribution on various types of the Fontan procedure [12, 14–20]. Because the Fontan circulation is supposed to be a low-energy system without ventricular power source, we assumed that the flow stagnation in addition to flow energy expenditure plays an essential role on evaluating flow dynamic efficiency. We schemed to evaluate flow stagnation and defined stagnation volume as the volume of the region where flow velocity was extremely low on the respiratory cycle in the conduit.
| Patients and Methods |
|---|
|
|
|---|
|
|
Computational Simulation
Three-dimensional geometries were divided into computational meshes (nodal points and elements). Nodal points ranged from 10,903 to 16,736, and the number of elements ranged from 61,007 to 95,170. Tetrahedral meshes were created using grid generator code GAMBIT 2.3.16 (Fluent Asia Pacific, Tokyo, Japan), which automatically checked the mesh quality. We took note to adopt identical mesh-interval lengths between the geometries in order to equalize the analytic accuracy. The finite element solver package FIDAP8.7.4 (Fluent Asia Pacific, Tokyo, Japan) was used to solve the Navier-Stokes equation of incompressible transient Newtonian fluid. The parameters were the following: density
= 1,060 kg/m3 and viscosity µ = 4.0 x 10–3 kg/m/s.
Energy Loss and Flow Stagnation
The formula of energy loss was justified by previous papers [8–19].
|
| (1) |
|
| (2) |
|
| (3) |
| Results |
|---|
|
|
|---|
|
|
Energy Loss
The results of the mean energy loss of each graft are listed in Table 2. Energy loss tended to increase as the conduit size decreased. The change was slightly larger between 14 mm (1.549 mW at rest, 3.529 mW during 0.5 W/kg exercise, 5.759 mW during 1.0 W/kg exercise) and 16 mm (1.423 mW at rest, 3.077 mW during 0.5 W/kg exercise, 4.881 mW during 1.0 W/kg exercise) than other conduit size intervals. Although energy loss became larger with the increase of the exercise level, we found these absolute values were significantly small even during strenuous exercise. The energy loss (5.759 mW) calculated during 1 W/kg exercise with the 14 mm graft was only equivalent to 0.118 kcal/day, which is trivial compared with the total energy intake or cardiac energy expenditure of young children 2 to 3 years old.
|
Stagnation Volume
Calculated results of stagnation volume of any conduit size in the expiratory and inspiratory phase of all exercise levels are listed in Table 3, and the stagnated regions in the conduits are dark blue in Figures 2, 4, and 5.
The stagnation volume tended to be larger in the expiratory phase than in the inspiratory phase according to the decrease of total flow conditions. With the increase of the exercise level, stagnation decreased (Fig 6). Differences of the stagnation volume between the conduit sizes were prominent. At rest, the 22 mm conduit had 33.9% of the stagnation volume in the expiratory phase (Table 3; Fig 6), whereas the 14 mm conduit had only 9.20% of the stagnation volume under the same conditions. Flow stagnation occurred at the lateral portion of the large size conduit (Figs 2; 4), where backflow at the expiratory phase and dilatory flow at the inspiratory phase occurred on the flow dynamics.
|
|
|
|
| Comment |
|---|
|
|
|---|
The flow collision between the SVC and IVC has been considered to be problematic, such as early circulatory failure after the Fontan operation [19, 25]. In the present study, nonsteady transient simulation with the respiration revealed the new findings: a penetrating flow from the SVC to the conduit, and the backflow in the lateral portion of the curved conduit in the expiratory phase. The results revealed that the larger conduit had larger flow penetration and larger backflow in the lateral portion of the conduit. This may be attributed not only to the conduit size itself but also to the geometric features of the extracardiac Fontan using large conduits. Although we set the offset between the center of the conduit and the center of the Glenn anastomosis site to be the length of one SVC-radius, the larger conduit had a larger facing area with the Glenn anastomosis site and smaller conduit curvatures. This may be the cause of the flow penetration. Because the previous Fontan simulation studies performed under the steady and continuous flow conditions did not reveal the flow penetration or the stagnated backflow at the lateral portion of the conduits, we thought this result revealed the mechanism of the "respiratory pump" (ie, the conduit acted as a reservoir in the expiratory phase and disgorged the reserved blood in the inspiratory phase). Also from this result, and from the fact that the backflow in the lateral portion of the conduit coincided well with the stagnation regions, we conclude that the penetrating flow and the backflow were the causes of the stagnation.
Efficiency of energy of the "respiratory pump" in the Fontan circulation during exercise has been discussed on clinical studies [3] and on numerical models [19, 20]. Our concern is the amount of energy dissipation due to flow collision of the SVC and IVC. Our results indicated that although energy dissipation at the anastomosis site increased during exercise, its hemodynamic attribution was sufficiently low compared with the value of the energy expenditure calculated by the previous CFD models [8–20]. Considering the form of equation (3) and the fact that Fontan circulation is a low-pressure and slow-flow system, energy loss could depend on pressure more than flow velocity. Pressure relaxation occurs in a large conduit space resulting in calculating low-energy loss. This idea may be supported by the previous numeric studies by de Leval and colleagues [9] and Bove and colleagues [13] that proved patch augmentation in the Fontan anastomosis site reduced the energy loss. On the contrary, the stenotic site where a pressure drop occurs might dissipate a large amount of energy, as evidenced in this calculation method that was first reported by Pekkan and colleagues [17]. The stenotic vessel may have a relatively large amount of wall friction in this low Reynolds number system. The difference between 14 mm and 16 mm was slightly larger than other conduit size intervals. The 14 mm conduit might not always seem to be optimal because it might soon become stenotic with the growth in body size and result in much more energy loss.
In order to evaluate stagnation, we defined an index "stagnation volume": the volume of the region where the absolute value of the velocity vector was small enough, less than the threshold value; 0.01 m/second. The stagnated region could be thought to mean high-risk area of thrombosis formation. This method of evaluation of thrombosis formation risk using a slow velocity region has been adopted in other CFD models such as thrombosis formation of coil-embolized cerebral arterial aneurysm [26, 27].
We adopted the threshold of 0.01 m/second for convenience. Because the hemodynamic mechanism of thrombosis formation has not yet fully been revealed, there is little evidence to determine the stagnation threshold.
Our results showed that a larger conduit had a larger stagnation volume. The stagnation region and flow pattern analysis discussed above revealed that a larger conduit had larger redundant volume on the lateral side. According to our results, grafts of less than 18 mm are thought to be optimal from the viewpoint of flow stagnation. However, the 14 mm conduit may be predicted to increase energy loss with body growth as it already has slightly higher energy loss than other conduits at this body size. The 16 mm and 18 mm conduits proved to be optimal for children of these ages.
Recent CFD models reconstruct patient-specific three-dimensional geometry from MRI or computed tomography to increase simulation accuracy [15, 18–20]. The present study is not a patient-specific one but is the idealized model. An idealized model has advantages for obtaining the essential or practical knowledge such as the indication or the strategies of the operation [17]. Our idealized geometry is based on the averaged data and has better agreement with the clinical observation compared with the previously reported idealized geometries [8–12, 17].
For giving the boundary conditions we adopted the previously reported MRI flow study data [3]. The patients studied in that report were older than those in the present study and older than the mean age of patients undergoing the Fontan operation in many institutions. Flow data during exercise of young patients (2- or 3- year-old patients) after having undergone the Fontan operation are not available. However, in the present study the given flow boundary data were corrected with the body surface area, and the flow ratio of the SVC and IVC (Table 1) were in good agreement with the clinical studies in the patients of these ages [28] and with the inlet boundary flow ratio of the previous CFD models [8–20].
In numeric models, several assumptions may oversimplify the multifactor system and become the study limitations. We did not consider the effect of fenestration, which would naturally change the flow fields. In the boundary conditions only respiratory fluctuations were considered, but actual flow patterns of the Fontan circulation contain pulsatile fluctuations which would make the flow character more complex. We also assumed the vessel wall to be rigid. We did not know the effect of compliances of the pulmonary arteries or the compressive effects of the conduit during the respiratory cycle. However, to confirm that these simplifications are acceptable, further studies are warranted.
We did not simulate adult Fontan cases. At the time of the operation, 20 mm or larger conduits could be considered too large; however, these sizes might become optimal after the patient grows larger. Simulation of adult size geometry, based on adult Fontan patient body-size data, should be performed using the same method.
According to our results with the Fontan operation, relative to body size, larger conduits proved to have redundant space. The 16 mm and 18 mm conduits proved to be optimal for children 2 to 3 years old. However, the optimal size considering the body growth warrants further study.
| Discussion |
|---|
|
|
|---|
For all the simulation you performed, you purposely built in flow competition between the SVC [superior vena cava] and the IVC [inferior vena cava] pathways. Do you think that this might have resulted majority of the energy loss within the TCPC [total cavopulmonary connection], and therefore hide the energy losses differences between various conduit sizes?
And my second question is, what are the Fontan pressures that you obtain in your simulations with exercise? This is a particularly important point because exercise modeling is not just increasing the flow at the boundaries from MRI [magnetic resonance imaging], as you have done. It is important to know what the corresponding pressure conditions are. With exercise, there are hormonal and neural influences that will change pulmonary and venous vascular elasticity and vascular compliance, which are not accounted for in your modeling. So I would be very interested to see what your pressure conditions are in your modeling as well.
DR ITATANI: I'll answer at first the second question. And actually pressure conditions during exercise were not available, so we defined the pressure condition of 11 mm Hg also during exercise. And we assumed that RP, resistance of pulmonary arteries, were not changed at rest or during exercise. So we assumed 11 mm Hg is the same pressure condition for convenience. And please repeat the first question.
DR HSIA: The first question is that your model's built in competitive flow between the SVC and IVC pathway. And we know that this generally results in the greatest energy loss. Therefore, if you have the majority of your energy loss occurring there, could that have artificially minimized the conduit size effects on energy loss?
DR ITATANI: Actually, in this study, SVC/IVC flow conflict is difficult to evaluate. But I think Fontan system is a low pressure and slow velocity system, and by the definition of the energy loss, which is also called Bernoulli's energy, dependent on pressure plus square of the velocity, energy loss depends much more on pressure gradient than flow conflict. So, in stenosis energy loss becomes drastically increased, so in our model only SVC/IVC flow conflict caused smaller energy loss.
DR JOHN HAWKINS (Salt Lake City, UT): I'd like to ask one quick question. Have these findings impacted your clinical care? Do you now place slightly smaller conduits than you would have four years ago or before you knew this data.
DR ITATANI: We had a case of small children operated Fontan operation, to whom we inserted a relatively large conduit. And after the Fontan operation, CT [computed tomographic] exam was taken and in this case thrombosis was found only at the lateral portion in the conduit.
DR CHRISTOPHER CALDARONE (Toronto, Ontario, Canada): Do you know, is there a relationship between 0.1 m/s blood flow speed and thrombosis.
DR ITATANI: I think this is a very good question. We define the threshold of 0.01 m/s for convenience because the fluid dynamic or hemodynamic mechanism of thrombosis formation was not revealed yet. So we wish to study further about the threshold. In this study we define the threshold for convenience.
DR CALDARONE: Just because we do have another minute, could we just have an informal survey of the audience. If you have, let's say, a 14 kg, 1-1/2 to 2-year old who is going to undergo an extracardiac Fontan, how many would use a 16-mm conduit? 18? 20? 22? (Show of hands.) It sounds like mostly 18s. It looked like a bell curve. All right. Thanks.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. Itatani, K. Miyaji, Y. Nakahata, K. Ohara, S. Takamoto, and M. Ishii The lower limit of the pulmonary artery index for the extracardiac Fontan circulation J. Thorac. Cardiovasc. Surg., July 1, 2011; 142(1): 127 - 135. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. M. Di Donato Invited Commentary Ann. Thorac. Surg., April 1, 2011; 91(4): 1246 - 1247. [Full Text] [PDF] |
||||
![]() |
S. Sathanandam, A. C. Polimenakos, C. Blair, C. El Zein, and M. N. Ilbawi Hypoplastic Left Heart Syndrome: Feasibility Study for Patients Undergoing Completion Fontan at or Prior to Two Years of Age Ann. Thorac. Surg., September 1, 2010; 90(3): 821 - 829. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ANN THORAC SURG | ASIAN CARDIOVASC THORAC ANN | EUR J CARDIOTHORAC SURG |
| J THORAC CARDIOVASC SURG | ICVTS | ALL CTSNet JOURNALS |